“…The models are mostly based on energy and mass conservation to quantify the coupling between the reactions of heat generation and the transport of heat and mass in porous solid materials. The model equations in general forms ,, are the energy balance equation for heat transfer, and mass balance equation for transport of the species, for example, oxygen, where ρ, C p , λ, and ε are the bulk density (kg/m 3 ), specific heat capacity (J/(kg·K)), effective thermal conductivity (W/(m·K)), and porosity (unitless) of the bulk material, respectively, ρ g , C p g , and u are the density (kg/m 3 ), specific heat capacity (J/(kg·K)), and flow velocity (m/s) of the gases in the pores, respectively, T is the temperature (K), t the time (s), C j is the concentration (mol/m 3 ) of species j , D j is the effective diffusivity (m 2 /s) of species j inside the bulk material, q i is the heat source i accounting for the rate of heat generation (W/m 3 ) from the physical, biological, or chemical process, and R j is the mass source or sink, representing the rate of production or consumption of species j (mol/(m 3 ·s)). In eqs and , the left-hand side is the accumulation term of heat or mass, the first term on the right-hand side represents the thermal conduction or species diffusion, and the second term on the right-hand side denotes the convective heat or gas transport.…”