Numerical Study of Cylindrical Tropical Woods Pyrolysis Using Python Tool

: In this paper, the thermal behavior of large pieces of wood pyrolysis has been modeled. Two mathematical models coupling heat transfer equations to chemical kinetics were used to predict the pyrolytic degradation of a 25 mm radius wood sample, assumed to be dry in the ﬁrst model and wet in the second, when heated to 973.15 K. The reactions involved in the pyrolysis process are assumed to be endothermic. The diffusion of bounded water during the process is taken into account in the second model, where the heat transfer equation has been coupled to that of the diffusion of moisture. This model, although simple, provides more information on the drying and pyrolysis processes during the heating of wood, which is its originality. It can therefore be advantageously used to calculate the temperature distribution in a pyrolysis bed. The equations of the two models, discretized by an explicit ﬁnite difference method, were solved numerically by a program written in Python. The validation of both models against experimental work in the literature is satisfactory. The two models allow examination of the temperature proﬁle in the radial direction of wood samples and highlighting of the effect of temperature on some thermal, physical and physicochemical characteristics.


Introduction
Renewable energies are of growing importance from an environmental point of view as a substitute for fossil fuels. These renewable energies use inexhaustible sources of energy of natural origin such as the sun, the wind, without forgetting the biomass which in the energy field of this document, will concern only the products of vegetable origin, or more precisely the wood. Wood from tropical forests is the object of a strong attraction both in tropical and temperate countries [1]. This wood is used not only to satisfy the construction needs of buildings, works of art and furniture, but also for the production of energy [1].
Africa is currently experiencing a period of sustained economic growth and transformation, with a rapidly increasing population and growing and diversifying economies. To be sustainable, such growth requires massive investment in the energy sector. In Africa, fossil fuels are scarce, costly to exploit as a source of energy or for the production of raw Many authors have simply considered a global reaction, converting wood into volatile products and charcoal [7][8][9][10], while others have used kinetic models involving two or more parallel or successive, irreversible, first-order reactions [11][12][13]. Most of these works were carried out on small wood samples and do not include, with a few exceptions, the effect of the water content of the wood. They assume initially dry wood samples. However, initially dry wood generally has at best an equilibrium moisture content of between 0.08 and 0.1 kg/kg in North Africa and between 0.2 and 0.12 kg/kg in Sub-Saharan Africa. It would therefore be somewhat incorrect to study pyrolysis from an anhydrous state of the wood.
In addition, when pyrolysis is conducted on small particles of low mass, physical limitations are overcome, chemical kinetics being the only limiting step. On the other hand, for massive samples, relevant to applications that are artisanal as well as industrial, coupling between the chemical kinetics and the physical limitations takes place. It is therefore necessary to take into account, during process modelling, the complex nature of wood (anisotropy, hygroscopy, etc.), the coupling between chemical kinetics, the heat transfer and the mass transfer which take place during the pyrolysis process.
Furthermore, experimental and theoretical studies have been carried out for large wet biomass particles by Shen et al. [12] and Chan et al. [14]. They translated the evaporation of bounded water during the heating process into a chemical drying reaction. Saastamoinen and Richard [15] proposed a mathematical model to evaluate the drying and pyrolysis of wood samples. They found that for large particles, drying and pyrolysis overlap during thermal decomposition, and that moisture in the particles significantly delays pyrolysis and reduces its rate. Furthermore, the volatile materials released during the pyrolysis process can lead to convective heat transfer near the surface, thus delaying the passage of the conductive wave into the solid [15,16].
Taking into account previous studies in literature reviews, and in order to optimize the heat treatment, we proposed, in this work, to model the thermal behavior of a large piece of tropical wood subjected to pyrolysis.

a Assumptions
To model the pyrolysis of solid wood samples, a number of assumptions must be made. For the purpose of this study, we will assume that:

•
The studied sample is of cylindrical geometry, of infinite length, heated through its surface by convection and radiation; • Samples are directly exposed to the heat fixed numerically, thus the heating rates are neglected. • Pyrolysis does not generate any change in the volume of the sample and the latter is considered isotropic; • For simplicity, heat is assumed to be transmitted inside the solid particle by conduction only; • The heat and mass transfer of the volatile products and the vapor inside the solid are ignored; • The volatile products and the vapor leave the solid as soon as they are produced. They are in thermal equilibrium with the solid matrix and therefore the rate of solid mass loss is taken as the mass flux of the volatiles; • Secondary reactions of the volatile products are not taken into account; • The chemical reactions of pyrolysis are described by the Arrhenius law of first order.
b Setting in equation of the pyrolysis of wood We propose to establish the equations that govern the pyrolysis of a cylindrical wood sample. The study of the process will be undertaken by means of two kinetic approaches: a uni-reactionary approach and a multi-reactionary one. The first is a global reaction, developing uniformly in the particle and converting the wood into volatile products and charcoal. The second, proposed by Shen et al. [12], is a one-step multi-reaction model formed by three competitive reactions (one reaction for the drying step and two others composing the pyrolysis step).
Based on the above assumptions, let us establish the material and energy balances.
Uni-reactional approach [7]: In this approach, wood is assumed to decompose into charcoal and volatile products in a single, endothermic, irreversible reaction of order 1. In addition to the assumptions made above, the physical properties of the sample, including C p and λ are assumed to be constant. The energy balance translating the heat accumulation, coming from the contribution of the diffusive heat flow and the reactions over time is coupled to the material balance which is the speed of the reaction. Bamford et al. [7] are written as follows (model 1): Multi-reactional approach [12]: In this approach, we consider all the phenomena that take place during the pyrolysis process of a wet wood sample, including the formation of water vapor, gases and charcoal [12]. The reactions involved in the process are characterized by three rate constants (associated with the formation of charcoal, gases and steam) obeying the Arrhenius law. For this second model, in addition to the coupling between heat transfer and chemical kinetics, we have taken into account the diffusion of bounded water during the heating process. This has allowed us to follow the spatio-temporal evolution of the humidity and to evaluate its influence. The heat balance and mass conservation equations for solids, gases, steam and liquid are written as follows (model 2): where: λ, in (W/m. K), is the thermal conductivity of the solid. It is calculated, taking into account the effect of moisture (H) and carbonaceous residue (λ c ), as follows: with: The humidity (H) is obtained by solving the following equation: D h = exp −9.9 + 9.8χ eq − 4300 T a (10) χ eq , is the moisture content at equilibrium. It is calculated according to the relation of Simpson and TenWolde taken from Simo-Tagne [17]: where RH is the relative humidity (%/100) and the parameters A, B, C and D are calculated as follows: with T a , the temperature of the surrounding air expressed in ( • C).
• Mass balances of wood, charcoal, gas, liquid and steam: • Initial and boundary conditions The initial and boundary conditions of the problem are: T is equal to the ambient air temperature. The density of the gases and steam is initially 0 kg/m 3 .
-Due to the cylindrical symmetry, the condition at the center gives: -At the surface of the sample, in r = R, we have: with: Under these conditions, the equations including the moisture content of the wood are associated with the multi-reaction approach. To solve the equations established in this study, we propose to use the explicit finite difference method. The differentiation scheme Sustainability 2021, 13, 13892 6 of 23 used is centered. The fictitious node approach has been used to discretize the boundary conditions [18]. The obtained equations were solved numerically by a program written in Python. c Numerical values of the model constants The numerical values of the parameters used in this work are summarized in the Table 1. In the absence of Comorian wood samples, the parameters (thermo-physical and kinetic) used in this study were all taken from the literature. To discretize all equations obtained above, a global explicit scheme has been used using the finite differences method. Thus, the uni-reactional approach gives: with: Multi-reaction approach Pyrolysis reaction Drying reaction −420 −420 −2440 [12] C p (J.kg −1 .K −1 ) κ a (m 2 .s −1 ) 2.77 × 10 −5 In this work In the case of multi reactional approach, we have: with: Equations given as boundary conditions have also been discretized. At the axis level (r = 0), we have: Case of uni-reactional approach: i = −1, is a fictive position with a fictive temperature T n −1 located at the left of the axis. Using Equation (21) with i = 0, we have: where: For the exterior surface (r = R for i = m), we have: i = m + 1, is a fictitious node at the fictitious temperature T n m+1 located at the right of the radial axis. Replacing the expression for T n m+1 in the equation (Equation

Results and Discussion
Based on the mathematical formulation of the pyrolysis process of solid wood samples, established in the previous sections, Python language has been used to generate numerical results. This model allowed us to examine the impact of the oven temperature, the heat conduction, the moisture content and the particle size on the pyrolysis process of a cylindrical wood sample.

Validation of the Model
The pyrolysis model was validated against the experiments of Larfeldt et al. [21]. The latter presented measurements of temperature distribution and mass loss on dry, solid wood cylindrical samples during pyrolysis in an inert atmosphere heated to 973.15 K. They considered cylindrical birch wood samples with a radius of 25 mm, a length of 300 mm and an average density of 410 kg/m 3 . Figure 1 presents the wood temperature evolutions at the surface (r = 25 mm), at 5 mm, and 15 mm to the axis of the cylindrical wood. In general, temperature increases from the surface to the center because parts near the surface are more heated than parts far from the surface. Model 1 (one-reactional model) gives acceptable results that are more satisfactory than model 2 (multi-reactional model).  Figure 2 shows relative density evolution versus time of process. Model 1 gives more satisfactory results than Model 2. Humidity in wood explains the poor results given by Model 2. In effect, when humidity is taken into account, water content absorbs an amount of energy and decreases the kinetics of the process as observed by several researchers [22]. In addition, the distance between experimental data and the models are presented in Figures 1 and 2 and can be explained by empirical data used in Models 1 and 2. Activation energy and frequency coefficients used are not specified to the studied wood. Thus, constants for the activation energies of production charcoal, gases and water vapor must be determined according to each wood type in order to have a good representation of experimental data. Also, all thermodynamic parameters, such as specific heat mass, enthalpy and thermal conductivity must be obtained after precise experiments. Using Figure 2, it is clear that thermodynamic parameters used for theoretical data must be obtained with experimental data. Figures 3 and 4 present evolutions of temperature and relative density of the cylindrical wood using inverse method. We found the values of A and E parameters to have theoretical curves near to the experimental curves given by Model 1. Thus, it was observed that when A and E are respectively equal to 1 × 10 5 s −1 and 127.4 kJ/kg, experimental and theoretical curves of the wood's relative density are close to each other. A space step of 0.05 s has been used. In effect, Di-Blasi et al. [23] show that pyrolysis reaction heat depends on the produced char and the type of studied materials. Using a cylindrical temperate sample of 3 cm diameter, Preau [24] found that temperatures of all parts of the sample increase until they reach the value of the oven (400 °C), this value is reached after 1000 s. At 1.5 cm diameter, the oven temperature (420 °C) is obtained after 400 s [25]. Influences of the heating rates experimentally considered can also increase this distance, as revealed in literature reviews [25,26].  In effect, when humidity is taken into account, water content absorbs an amount of energy and decreases the kinetics of the process as observed by several researchers [22]. In addition, the distance between experimental data and the models are presented in Figures 1 and 2 and can be explained by empirical data used in Models 1 and 2. Activation energy and frequency coefficients used are not specified to the studied wood. Thus, constants for the activation energies of production charcoal, gases and water vapor must be determined according to each wood type in order to have a good representation of experimental data. Also, all thermodynamic parameters, such as specific heat mass, enthalpy and thermal conductivity must be obtained after precise experiments. Using Figure 2, it is clear that thermodynamic parameters used for theoretical data must be obtained with experimental data. Figures 3 and 4 present evolutions of temperature and relative density of the cylindrical wood using inverse method. We found the values of A and E parameters to have theoretical curves near to the experimental curves given by Model 1. Thus, it was observed that when A and E are respectively equal to 1 × 10 5 s −1 and 127.4 kJ/kg, experimental and theoretical curves of the wood's relative density are close to each other. A space step of 0.05 s has been used. In effect, Di-Blasi et al. [23] show that pyrolysis reaction heat depends on the produced char and the type of studied materials. Using a cylindrical temperate sample of 3 cm diameter, Preau [24] found that temperatures of all parts of the sample increase until they reach the value of the oven (400 • C), this value is reached after 1000 s. At 1.5 cm diameter, the oven temperature (420 • C) is obtained after 400 s [25]. Influences of the heating rates experimentally considered can also increase this distance, as revealed in literature reviews [25,26]. Figure 5 presents the heat mass gradients when the temperatures in the oven are 673.15 K and 973.15 K. In a given time of the process and at the same distance of the center of the cylindrical wood sample, temperature is high when the oven temperature increases. From surface to the center, temperature decreases. When the process takes longer, internal parts of the cylindrical wood are hotter and the temperature increases. The thermal degradation of wood will begin from the surface to the center of samples. The average temperature increases with the duration of the process, but the limit of this evolution is near the temperature of the oven and this is because of the endothermic nature of the pyrolysis reaction, assumed in our simulation. During the first hour of the process, Figure 6 shows that the rate of temperature evolution increases with the temperature of the oven (source), and with the duration of the process.       Figure 5 presents the heat mass gradients when the temperatures in the oven are 673.15 K and 973.15 K. In a given time of the process and at the same distance of the center of the cylindrical wood sample, temperature is high when the oven temperature increases. From surface to the center, temperature decreases. When the process takes longer, internal parts of the cylindrical wood are hotter and the temperature increases. The thermal degradation of wood will begin from the surface to the center of samples. The average temperature increases with the duration of the process, but the limit of this evolution is near the temperature of the oven and this is because of the endothermic nature of the pyrolysis reaction, assumed in our simulation. During the first hour of the process, Figure 6 shows that the rate of temperature evolution increases with the temperature of the oven (source), and with the duration of the process.    Figure 5 presents the heat mass gradients when the temperatures in the oven are 673.15 K and 973.15 K. In a given time of the process and at the same distance of the center of the cylindrical wood sample, temperature is high when the oven temperature increases. From surface to the center, temperature decreases. When the process takes longer, internal parts of the cylindrical wood are hotter and the temperature increases. The thermal degradation of wood will begin from the surface to the center of samples. The average temperature increases with the duration of the process, but the limit of this evolution is near the temperature of the oven and this is because of the endothermic nature of the pyrolysis reaction, assumed in our simulation. During the first hour of the process, Figure 6 shows that the rate of temperature evolution increases with the temperature of the oven (source), and with the duration of the process.     After 600 s, only 5 mm of the wood is degraded when the temperature of the oven is equal to 673.15 K (see Figure 7a). At 973.15 K, 10 mm of wood is degraded, with 5 mm entirely degraded (see Figure 7b). A high heating temperature accelerates the degradation process. Using a low radius of wood particles with integration of the oxidation of char, Benkoussas et al. obtained the similar variations [27].  After 600 s, only 5 mm of the wood is degraded when the temperature of the oven is equal to 673.15 K (see Figure 7a). At 973.15 K, 10 mm of wood is degraded, with 5 mm entirely degraded (see Figure 7b). A high heating temperature accelerates the degradation process. Using a low radius of wood particles with integration of the oxidation of char, Benkoussas et al. obtained the similar variations [27].  Figure 8 shows that the thermal degradation of wood proceeds progressively from the surface to the center. Velocity of degradation is high with the oven temperature. The propagation of the pyrolysis front through the non-pyrolyzed solid is shown in Figure 8 where the spatial evolution of the pyrolysis reaction rate, at various times and for a heating temperature of 673.15 K and 973.15 K, is represented. The reaction rate increases rapidly, reaches a maximum and then decreases to zero. Thus, volatiles are initially produced rapidly; the rate of production slows as the pyrolysis process continues. Figure 9 shows that the average velocity of degradation increases during the start of the process but some times decreases during the process. In effect, during the degradation, wood thermal capacity decreases when wood average temperature is near the temperature of the source. Figure 10 presents the wood conversion rate versus the duration of the process. It shows that wood is progressively transformed to other components during the process. After 60 min, 57% of the wood has been transformed when the oven temperature is equal to Tf = 673.15 K. At 973.15 K, cylindrical wood is totally transformed after 40 min. The conversion should vary from 0 at the beginning of the pyrolysis to 1 at the end. As the sample is considered dry, no latency period, which would be due to moisture evaporation, is to be noted in both cases. However, at low heating temperatures, the conversion increases  Figure 8 shows that the thermal degradation of wood proceeds progressively from the surface to the center. Velocity of degradation is high with the oven temperature. The propagation of the pyrolysis front through the non-pyrolyzed solid is shown in Figure 8 where the spatial evolution of the pyrolysis reaction rate, at various times and for a heating temperature of 673.15 K and 973.15 K, is represented. The reaction rate increases rapidly, reaches a maximum and then decreases to zero. Thus, volatiles are initially produced rapidly; the rate of production slows as the pyrolysis process continues. Figure 9 shows that the average velocity of degradation increases during the start of the process but some times decreases during the process. In effect, during the degradation, wood thermal capacity decreases when wood average temperature is near the temperature of the source. Figure 10 presents the wood conversion rate versus the duration of the process. It shows that wood is progressively transformed to other components during the process. After 60 min, 57% of the wood has been transformed when the oven temperature is equal to T f = 673.15 K. At 973.15 K, cylindrical wood is totally transformed after 40 min. The conversion should vary from 0 at the beginning of the pyrolysis to 1 at the end. As the sample is considered dry, no latency period, which would be due to moisture evaporation, is to be noted in both cases. However, at low heating temperatures, the conversion increases slowly due to the fact that the temperature of the sample has not reached the temperature of the beginning of pyrolysis; during this time, the sample is only heated but no degradation takes place.           Along with Preau, Figure 11 permits the conclusion that sample temperature increases rapidly with a low diameter and an increasing temperature is influenced by the value of the diameter of samples [24]. Figure 11 presents the influence of the size of cylindrical wood on the evolution of temperature ( Figure 11a) and wood conversion rate (Figure 11b) at T f = 973.15 K. When the mass of wood increases, temperature takes time to heat all of the product. After 1000 s, all cylindrical wood of 15 mm radius is totally transformed, 2500 s are needed to totally transform a cylindrical wood with 25 mm of radius, while after 3000 s cylindrical wood with a radius higher than 35 mm is not totally transformed. The temperature increases more rapidly with small diameter samples as observed by Preau [24] in his experiments on the pyrolysis of solid wood samples. Figure 12 presents the influence of initial moisture content on the surface temperature of cylindrical wood versus the duration of the process using Model 2 (multi-reactional model) at T f = 973.15 K, R = 25 mm, RH = 0.3. When the wood is in non-hygroscopic domain (moisture content higher than 30%), free water is rapidly dried on the surface and evolutions of temperature at the surface of the cylindrical wood are almost the same. In hygroscopic domain (moisture content lower than 30%), the surface of the wood heated rapidly. However, after 400 s all curves align because the effects of the initial moisture content have been cancelled. Figure 13 shows that average moisture content in the cylindrical wood decreases exponentially and all curves align after 400 s. The samples are totally dried after 800 s. Figure 14 shows the distribution of temperature during the process at five positions of the cylindrical wood using Model 2. It is clear that wood samples will degrade progressively from the surface to the center because temperature takes time to heat internal parts of the cylindrical wood. After 500 s, the temperature of the surface is near the temperature of the source (973.15 K) but the temperature at the center has not changed. After 3500 s, the temperature of the position located at 10 mm to the axis of the cylindrical wood is equal to 750 K. Literature also shows similar experimental results [24]. Along with Preau, Figure 11 permits the conclusion that sample temperature increases rapidly with a low diameter and an increasing temperature is influenced by the value of the diameter of samples [24]. Figure 11. Influences of the wood size on wood temperature (a) and wood conversion rate (b) evolutions versus the duration using Model 1 at Tf = 973.15 K. Figure 11 presents the influence of the size of cylindrical wood on the evolution of temperature ( Figure 11a) and wood conversion rate (Figure 11b) at Tf = 973.15 K. When the mass of wood increases, temperature takes time to heat all of the product. After 1000 s, all cylindrical wood of 15 mm radius is totally transformed, 2500 s are needed to totally transform a cylindrical wood with 25 mm of radius, while after 3000 s cylindrical wood with a radius higher than 35 mm is not totally transformed. The temperature increases more rapidly with small diameter samples as observed by Preau [24] in his experiments on the pyrolysis of solid wood samples. Figure 12 presents the influence of initial moisture content on the surface temperature of cylindrical wood versus the duration of the process using Model 2 (multi-reactional model) at Tf = 973.15 K, R = 25 mm, RH = 0.3. When the wood is in non-hygroscopic domain (moisture content higher than 30%), free water is rapidly dried on the surface and evolutions of temperature at the surface of the cylindrical wood are almost the same. In hygroscopic domain (moisture content lower than 30%), the surface of the wood heated rapidly. However, after 400 s all curves align because the effects of the initial moisture  Figure 13 shows that average moisture content in the cylindrical wood decreases exponentially and all curves align after 400 s. The samples are totally dried after 800 s. Figure 14 shows the distribution of temperature during the process at five positions of the cylindrical wood using Model 2. It is clear that wood samples will degrade progressively from the surface to the center because temperature takes time to heat internal parts of the cylindrical wood. After 500 s, the temperature of the surface is near the temperature of the source (973.15 K) but the temperature at the center has not changed. After 3500 s, the temperature of the position located at 10 mm to the axis of the cylindrical wood is equal to 750 K. Literature also shows similar experimental results [24].      Figure 15 shows that the cylindrical wood surface rapidly equaled the equilibrium moisture content. Internal parts of the wood are heated rapidly when temperature of the source was high. When Tf = 673.15 K, the axial part of the wood had 0.175 kg/kg moisture content after 600 s duration (Figure 15a). After 1200 s, the axial part of the wood was equal to 0.015 kg/kg (Figure 15b). It is clear from these curves that high heating temperatures accelerate the drying process. Figure 16 shows that after 600 s, the temperature of parts located at 10 mm of the axis have the same temperature as at the start of the process. Thus, variation of moisture content near the axis presented in Figure 15 is not explained by an increase in the temperature. However, as found by Di Blasi et al. [22], the presence of moisture slows down the pyrolysis process to the point that near the center, the temperature remains constant at its initial value during the evaporation of moisture. Once the drying process is over, in which case the region near the center has reached the equilibrium moisture content, this temperature increases. Indeed, a significant amount of energy is used for the evaporation of the water contained in the wood and only a small amount is transferred to the virgin wood. When the temperature at the surface of the sample reaches about 500 K, the pyrolysis reaction starts and a layer of charcoal forms near the  Figure 15 shows that the cylindrical wood surface rapidly equaled the equilibrium moisture content. Internal parts of the wood are heated rapidly when temperature of the source was high. When T f = 673.15 K, the axial part of the wood had 0.175 kg/kg moisture content after 600 s duration (Figure 15a). After 1200 s, the axial part of the wood was equal to 0.015 kg/kg (Figure 15b). It is clear from these curves that high heating temperatures accelerate the drying process. Figure 16 shows that after 600 s, the temperature of parts located at 10 mm of the axis have the same temperature as at the start of the process. Thus, variation of moisture content near the axis presented in Figure 15 is not explained by an increase in the temperature. However, as found by Di Blasi et al. [22], the presence of moisture slows down the pyrolysis process to the point that near the center, the temperature remains constant at its initial value during the evaporation of moisture. Once the drying process is over, in which case the region near the center has reached the equilibrium moisture content, this temperature increases. Indeed, a significant amount of energy is used for the evaporation of the water contained in the wood and only a small amount is transferred to the virgin wood. When the temperature at the surface of the sample reaches about 500 K, the pyrolysis reaction starts and a layer of charcoal forms near the surface, while the temperature in the center remains constant at its initial value. It can therefore be concluded that drying and pyrolysis take place successively, at high heating temperatures, at the same positions along the radius of the pyrolyzed solid wood sample.
accelerate the drying process. Figure 16 shows that after 600 s, the temperature of parts located at 10 mm of the axis have the same temperature as at the start of the process. Thus, variation of moisture content near the axis presented in Figure 15 is not explained by an increase in the temperature. However, as found by Di Blasi et al. [22], the presence of moisture slows down the pyrolysis process to the point that near the center, the temperature remains constant at its initial value during the evaporation of moisture. Once the drying process is over, in which case the region near the center has reached the equilibrium moisture content, this temperature increases. Indeed, a significant amount of energy is used for the evaporation of the water contained in the wood and only a small amount is transferred to the virgin wood. When the temperature at the surface of the sample reaches about 500 K, the pyrolysis reaction starts and a layer of charcoal forms near the surface, while the temperature in the center remains constant at its initial value. It can therefore be concluded that drying and pyrolysis take place successively, at high heating temperatures, at the same positions along the radius of the pyrolyzed solid wood sample.   Figure 17 shows that the rate of temperature increase is high when the temperature of the source increases. This average temperature continues to increase until the temperature of all parts of the wood is near to that of the source Tf. Figure 18 shows that the density of wood decreases progressively. When the temperature of the source is equal to 673.15 K, wood density starts to decrease after 10 min. When Tf = 973.15 K, the density of the wood decreases after the first seconds of the process. In effect, Figure 16 shows that after 10 min of the process at Tf = 673.15 K, some millimeters of the wood's surface reaches a temperature equal to 500 K, the temperature from which charcoal is produced [22]. When Tf = 973.15 K, a part of the sample reaches 500 K at 150 s, thus charcoal is produced rapidly. Figures 19 and 20 show that when the wood density decreases, charcoal and gas are directly produced. Figure 19a,b show that during slow pyrolysis (Figure 19a), the density of charcoal obtained is lower than that obtained during the fast pyrolysis (Figure 19b). Values of densities found are not far from those obtained by Teixeira et al. [28] using wood pellets and wood chips at 700 °C. Figure 20 shows that the gas obtained has a density greater than 10 kg/m 3 after the first minutes of the process. Thus, it is possible to have some gas in a liquid state in order to justify the great value of gas density obtained. As showed in the literature, gas produced by the wood pyrolysis process is a mixture of many known gases such as CO, CH4, C, CO2, H2, H2O and O2 [28]. In the liquid state, it is possible that these gases give a high density value. In addition, it has been shown that fast and  Figure 17 shows that the rate of temperature increase is high when the temperature of the source increases. This average temperature continues to increase until the temperature of all parts of the wood is near to that of the source T f . Figure 18 shows that the density of wood decreases progressively. When the temperature of the source is equal to 673.15 K, wood density starts to decrease after 10 min. When T f = 973.15 K, the density of the wood decreases after the first seconds of the process. In effect, Figure 16 shows that after 10 min of the process at T f = 673.15 K, some millimeters of the wood's surface reaches a temperature equal to 500 K, the temperature from which charcoal is produced [22]. When T f = 973.15 K, a part of the sample reaches 500 K at 150 s, thus charcoal is produced rapidly. Figures 19 and 20 show that when the wood density decreases, charcoal and gas are directly produced. Figure 19a,b show that during slow pyrolysis (Figure 19a), the density of charcoal obtained is lower than that obtained during the fast pyrolysis (Figure 19b). Values of densities found are not far from those obtained by Teixeira et al. [28] using wood pellets and wood chips at 700 • C. Figure 20 shows that the gas obtained has a density greater than 10 kg/m 3 after the first minutes of the process. Thus, it is possible to have some gas in a liquid state in order to justify the great value of gas density obtained.
As showed in the literature, gas produced by the wood pyrolysis process is a mixture of many known gases such as CO, CH 4 , C, CO 2 , H 2 , H 2 O and O 2 [28]. In the liquid state, it is possible that these gases give a high density value. In addition, it has been shown that fast and slow pyrolysis give 13% and 35% of gas, respectively, near the solid products ratio given by the process [29]. Commandré et al. [30] confirm this ratio of gas density obtained during pyrolysis process and show experimentally that this ratio can reach up to 40%.
rapidly. Figures 19 and 20 show that when the wood density decreases, charcoal and gas are directly produced. Figure 19a,b show that during slow pyrolysis (Figure 19a), the density of charcoal obtained is lower than that obtained during the fast pyrolysis (Figure 19b). Values of densities found are not far from those obtained by Teixeira et al. [28] using wood pellets and wood chips at 700 °C. Figure 20 shows that the gas obtained has a density greater than 10 kg/m 3 after the first minutes of the process. Thus, it is possible to have some gas in a liquid state in order to justify the great value of gas density obtained. As showed in the literature, gas produced by the wood pyrolysis process is a mixture of many known gases such as CO, CH4, C, CO2, H2, H2O and O2 [28]. In the liquid state, it is possible that these gases give a high density value. In addition, it has been shown that fast and slow pyrolysis give 13% and 35% of gas, respectively, near the solid products ratio given by the process [29]. Commandré et al. [30] confirm this ratio of gas density obtained during pyrolysis process and show experimentally that this ratio can reach up to 40%.         [31,32].
Vapor and liquid states are highly related. In effect, Figures 21 and 22 show that the evolutions of these two pyrolysis products are directly opposite because the concentration  [31,32].
Vapor and liquid states are highly related. In effect, Figures 21 and 22 show that the evolutions of these two pyrolysis products are directly opposite because the concentration of vapor is influenced by the pressure and temperature of liquid water in the cylindrical wood. Thus, liquid density firstly decreases ( Figure 21) in order to produce the liquid vapor ( Figure 22). It has been seen that during the process, liquid in the wood is firstly transformed to gas (steam). If this gas is not evacuated to the system, it is transformed to a liquid form again (see Figure 22). of vapor is influenced by the pressure and temperature of liquid water in the cylindrical wood. Thus, liquid density firstly decreases ( Figure 21) in order to produce the liquid vapor ( Figure 22). It has been seen that during the process, liquid in the wood is firstly transformed to gas (steam). If this gas is not evacuated to the system, it is transformed to a liquid form again (see Figure 22).   Figure 23 presents the wood conversion rate versus duration of the process in the two studied temperatures. It is clear that when the temperature of the source increases, wood is degraded rapidly. Thus, charcoal, gas and water vapor are also rapidly produced of vapor is influenced by the pressure and temperature of liquid water in the cylindrical wood. Thus, liquid density firstly decreases ( Figure 21) in order to produce the liquid vapor ( Figure 22). It has been seen that during the process, liquid in the wood is firstly transformed to gas (steam). If this gas is not evacuated to the system, it is transformed to a liquid form again (see Figure 22).   Figure 23 presents the wood conversion rate versus duration of the process in the two studied temperatures. It is clear that when the temperature of the source increases, wood is degraded rapidly. Thus, charcoal, gas and water vapor are also rapidly produced (Figures 24-26). At 673.15 K, the average velocity of charcoal production increases progres-  Figure 23 presents the wood conversion rate versus duration of the process in the two studied temperatures. It is clear that when the temperature of the source increases, wood is degraded rapidly. Thus, charcoal, gas and water vapor are also rapidly produced (Figures 24-26). At 673.15 K, the average velocity of charcoal production increases progressively until 0.025 kg/(m 3 .s) and oscillates near this value. When T f = 973.15 K, a great variation is presented during the first 10 min of the process before reaching a stability near 0.1 kg/(m 3 .s). These two variations are also seen during production of gas and water vapor. At a given temperature, water vapor is more rapidly produced than charcoal, and charcoal more rapidly produced than gas (Figures 24-26).          The influence of wood pyrolysis temperature is clearly presented in Figure 23. The increasing temperature increases the wood conversion and changes the physical morphology of the charcoal (volume, form), thus the charcoal density as shown by Ananias

Conclusions
In this work, we have mathematically modeled the pyrolysis process of large pieces of wood using two kinetic models: a uni-reactionary and a multi-reactionary one step model. In the absence of Comorian wood samples, the parameters (thermo-physical and kinetic) used in this study were all extracted from the literature. The mathematical formulation of the model was used as a basis to develop a numerical model in Python.
The suggested numerical model has allowed us to follow the evolution of the pyrolysis reaction and to better understand the influence of the oven temperature, the particle size and the water content of the wood in the pyrolytic conversion process. We note that during the pyrolysis of solid wood samples, the drying and pyrolysis reactions take place successively. The water content is a factor that slows down the pyrolysis process and it is therefore recommended to use dry samples. The oven temperature has a significant influence on the temperature profile and the reaction rates; a high oven temperature reduces the drying and conversion time while it increases the rate of the production of volatile products. A wet wood sample has no obvious drying behavior at high heating temperatures. The size of the particles is also a factor influencing the process, which can longer or shorter depending on the size of the sample.
Among the two studied models in this work, the multi-reaction model better describes the process because it takes into account the effect of the water content and allows us to follow the evolution of various substances produced during the pyrolysis. This model, which is simple and does not require a long calculation time, can be used to calculate the temperature distribution in a pyrolysis bed.
However, since water content, furnace temperature and particle size are not the only factors to be controlled for an efficient and optimal pyrolysis, further studies are needed to better understand the process.
In order to specify the results of this study to Comorian wood samples, it is recommended that in future work, an extensive experimental measurement campaign be carried out to determine frequency factors, densities, thermal and mass conductivities, gas permeability, activation energies, desorption isotherms, among others, so that the numerical results of pyrolysis/combustion are more specific to Comorian woods and subsequently, The influence of wood pyrolysis temperature is clearly presented in Figure 23. The increasing temperature increases the wood conversion and changes the physical morphology of the charcoal (volume, form), thus the charcoal density as shown by Ananias et al., Hanif et al. and Arous et al. [33][34][35].

Conclusions
In this work, we have mathematically modeled the pyrolysis process of large pieces of wood using two kinetic models: a uni-reactionary and a multi-reactionary one step model. In the absence of Comorian wood samples, the parameters (thermo-physical and kinetic) used in this study were all extracted from the literature. The mathematical formulation of the model was used as a basis to develop a numerical model in Python.
The suggested numerical model has allowed us to follow the evolution of the pyrolysis reaction and to better understand the influence of the oven temperature, the particle size and the water content of the wood in the pyrolytic conversion process. We note that during the pyrolysis of solid wood samples, the drying and pyrolysis reactions take place successively. The water content is a factor that slows down the pyrolysis process and it is therefore recommended to use dry samples. The oven temperature has a significant influence on the temperature profile and the reaction rates; a high oven temperature reduces the drying and conversion time while it increases the rate of the production of volatile products. A wet wood sample has no obvious drying behavior at high heating temperatures. The size of the particles is also a factor influencing the process, which can longer or shorter depending on the size of the sample.
Among the two studied models in this work, the multi-reaction model better describes the process because it takes into account the effect of the water content and allows us to follow the evolution of various substances produced during the pyrolysis. This model, which is simple and does not require a long calculation time, can be used to calculate the temperature distribution in a pyrolysis bed.
However, since water content, furnace temperature and particle size are not the only factors to be controlled for an efficient and optimal pyrolysis, further studies are needed to better understand the process.
In order to specify the results of this study to Comorian wood samples, it is recommended that in future work, an extensive experimental measurement campaign be carried out to determine frequency factors, densities, thermal and mass conductivities, gas permeability, activation energies, desorption isotherms, among others, so that the numerical results of pyrolysis/combustion are more specific to Comorian woods and subsequently, to dimension pyrolysis/combustion/gasification equipment taking into account the specificity of these woods. Influence of the heating rates must be also studied in order to describe progressively very well all thermophysical processes.