Numerical simulation of the mass loss process in pyrolizing char materials

We present here a decoupling technique to tackle the entanglement of the nonlinear boundary condition and the movement of the char/virgin front for a thermal pyrolysis model for charring materials. Standard numerical techniques to solve moving front problems — often referred to as Stefan problems — encounter difficulties when dealing with nonlinear boundaries. While special integral methods have been developed to solve this problem, they suffer from several limitations which the technique described here overcomes. The newly developed technique is compared with the exact analytical solutions for some simple ideal situations which demonstrate that the numerical method is capable of producing accurate numerical solutions. The pyrolysis model is also used to simulate the mass loss process from a white pine sample exposed to a constant radiative flux in a nitrogen atmosphere. Comparison with experimental results demonstrates that the predictions of mass loss rates and temperature profile within the solid material are in good agreement with the experiment.