Identifying a modeling procedure of biomass thermal decomposition that is not only simple enough to implement and use, and computationally efficient, but also sufficiently accurate for engineering design activities, and with a spectrum of applications as broad as possible is a very difficult task. The authors propose a procedure which consists of two main stages: (a) the static modeling phase with the purpose of generating the algorithm (macro functions) that supplies a Computational Fluid Dynamics (CFD) model with specific input data (source/sink terms and local material properties) and (b) the dynamic modeling phase, where the CFD model is bi-directionally coupled to the external biomass decomposition model in the form of a User-Defined Function (UDF). The modeling approach was successfully validated against data obtained from single particle decomposition experiments, demonstrating its applicability even to large biomass particles, under high heating rates and combusting conditions.