8.1. Fundamentals#

As stated in section 1.1, pyrolysis describes the emission of gases from the solid material. It is the thermal decomposition of materials under inert atmosphere. Whereby thermal decomposition is defined as a process of extesive chemical species change caused by heat. The gases released from the solid may afterwards react with oxygen, this reaction is not part of the pyrolysis.

Studying pyrolysis with the TGA#

A thermogravimetric analyser (TGA) is often used to study pyrolysis. In a tga the mass of the sample is so small (usually between \(\mf 5~mg\) to \(\mf 10~mg\)) that it is assumed that heat and mass transfer can be neglected. A picture of a tga device can be seen in Fig. 8.1.

../../../_images/TGA.jpg

Fig. 8.1 TGA#

In a TGA a small amount of sample mass is placed in a crucible. This crucible is placed in the TGA. During the measurement the device will heat with a specified heating rate (e.g. \(\mf 5~K/min\)), while the mass loss of the sample is measured. This can be done under different atmospheres, typically experiments are done either under inert atmosphere (e.g. nitrogen or argon) or under synthetic air. The measurement will give the mass loss as function of time. This can be used to determine the massloss rate as function of temperature, which can be used to determine the kinetic parameters of the material.

../../../_images/TGA1.PNG

Fig. 8.2 Mass loss curve from TGA data.#

../../../_images/TGA2.PNG

Fig. 8.3 Mass loss rate curve from TGA data#

Arrhenius equation#

To describe a pyrolysis process Arrhenius equations are usually used. Suppose an amount of one specific material (X) is completly converted to gaseous fuel with one reaction. This pyrolysis reaction can than be described by a rate equation and the Arrhenius equation:

\[\mf \frac{dX}{dt} = - kX^n\]
\[\mf k = A e^{\frac{-E_a}{RT}}\]

where \(\mf k\) is the rate constant and \(\mf X\) is the amount of material that takes part in the reaction. The kinetic parameters are

  • \(\mf E_a\): activation energy

  • \(\mf A\): pre-exponential factor

  • \(\mf n\): reaction order

Modeling with FDS#

There are several ways to model pyrolysis in fds. The simplest way is to specify the heat release rate. In this case the pyrolising sample will be modeled as ejection of gaseous fuel from the surface.

However fds, also allows for more complex pyrolysis modelling. Every material described in the fds input file can undergo multiple reactions. The amount of reactions should be specified by N_REACTIONS for every material. For every single reaction the reaction parameters, the created solid material, the released gas species, and others should be specified in the input file.

For each reaction the kinetic parameters need to be specified. The general evolution equation is given by:

\[\mf \frac{d Y_{s,i}}{d t} = - \sum^{N_{r,i}}_{j=1} r_{ij} + \sum^{N_{m}}_{i'=1} \sum^{N_{r,i'}}_{j=1} \nu_{s,i'j} r_{i'j} \;\;\;\;\; (i' \neq i) \]
\[\mf r_{ij} = A_{ij}Y_{s,i}^{n_{s,ij}} exp \left(-\frac{E_{ij}}{ R T_s} \right) X_{O_2}^{n_{O_2,ij}}\]
\[\mf Y_{s,i} = \left( \frac{\rho_{s,i}}{\rho_s(0)}\right)\]

where:

  • \(\mf r_{ij}\) : rate of reaction at the tempeature T_s of the ith material undergoing the jth reaction

  • \(\mf \nu_{s, i'j}\) : yield of production of material i by other reactions

  • \(\mf \rho_{s,i}\) : density of the ith matierial of the layer

  • \(\mf \rho_s(0)\) : original density of the layer

  • \(\mf X_{O_2}\) : local oxygen volume fraction