Numerical and Experimental Investigation of Equivalence Ratio ( ER ) and Feedstock Particle Size on Birchwood Gasification

This paper discusses the characteristics of Birchwood gasification using the simulated results of a Computational Fluid Dynamics (CFD) model. The CFD model is developed and validated with the experimental results obtained with the fixed bed downdraft gasifier available at the University of Agder (UIA), Norway. In this work, several parameters are examined and given importance, such as producer gas yield, syngas composition, lower heating value (LHV), and cold gas efficiency (CGE) of the syngas. The behavior of the parameters mentioned above is examined by varying the biomass particle size. The diameters of the two biomass particles are 11.5 mm and 9.18 mm. All the parameters investigate within the Equivalences Ratio (ER) range from 0.2 to 0.5. In the simulations, a variable air inflow rate is used to achieve different ER values. For the different biomass particle sizes, CO, CO2, CH4, and H2 mass fractions of the syngas are analyzed along with syngas yield, LHV, and CGE. At an ER value of 0.35, 9.18 mm diameter particle shows average maximum values of 60% of CGE and 2.79 Nm3/h of syngas yield, in turn showing 3.4% and 0.09 Nm3/h improvement in the respective parameters over the 11.5 mm diameter biomass particle.


Introduction
Among the available energy sources, biomass is envisaged to play a major role in the future energy supplement.It produces no net carbon emission, while being the fourth largest energy source available in the world [1][2][3][4].Therefore, biomass has high potential in contributing to satisfy the future energy demands of the world.Further, it is seen that for countries where the economy is mostly based on agriculture, they can utilize the potential of biomass efficiently.To recover energy from biomass, either thermochemical or biochemical process can be used [5][6][7].The gasification process is a thermochemical process which gives a set of gases as output, consisting of CO, CO 2 , H 2 , CH 4, and N 2 by converting organic or carbonaceous materials like coal or biomass [8][9][10].In gasification, several types of reactors used, such as entrained flow, fixed bed, fluidized bed, and moving bed gasifiers [7,[11][12][13].
The producer gas is a result of a set of endothermic and exothermic reactions.The required heat for the endothermic reactions is provided by the exothermic reactions [11,14].Once the steady state is achieved, the gasifier could work in a certain temperature range.Then, the producer gas can be obtained, as long as the fuel is being fed [14,15].
It was proven that producer gas was more versatile and useful than the biomass [16].The quality of the producer gas is established to be one of the most important aspects of the gasification which should be enhanced, as it is used to generate power [11,16].Chemical reactions inside the gasifier consist of both homogeneous and heterogeneous reactions.There are some factors which are important for the improvement of the gasification process such as reactor temperature, reactor pressure, the flow Both heterogeneous and homogeneous reactions are considered and variation of syngas composition, syngas yield, LHV of the syngas, and CGE are measured for two sizes of feedstock, as well as for ER.Fixed bed downdraft gasifiers have demonstrated to perform well in the ER range from 0.25 to 0.43.Therefore, in the present work, the variations of the gasification parameters are considered in the ER range from 0.2 to 0.5 [14,15].For the experimental work, the downdraft gasifier is used, which is available at the University of Agder (UIA), Norway.In addition, the CFD model is developed as the same geometry of the gasifier used for the experiments.For the validation purposes, experiments are done with the same gasifier with one of the wood particle sizes (11.5 mm woodchip) as it is the only biomass particle size available for the experiments.Then, with the developed simulation model, both biomass particle sizes are simulated.

Feedstock Characteristics
For the study, Birchwood is used as the feedstock.The ultimate and proximate analysis results for the Birchwood are listed down in Table 1.The data is extracted from a previous study [33], which is done on Birchwood in the same lab and under the same conditions.Table 1.Proximate and ultimate analysis results of Birchwood [33].According to the analysis above, Birchwood has a lower amount of moisture in comparison to some of the other wood types [33][34][35], which is an important aspect for gasification.Apart from the moisture content, Birchwood has shown a lower amount of ash, Cl, and S, while higher values of calorific value and volatile matter [34] has been indicated.It is obvious that Birchwood has a higher calorific value because of the higher amount of carbon content [33].In addition, when a biomass consists of less ash, it can lead to a higher conversion process which has lesser slag [14].Hence, Birchwood has some positive characteristics in the gasification perspective.According to some of the literature, K, Si, and Ca can have a small effect in determining gasification characteristics.However, in this work, those are not considered as they are not measured in the ultimate analysis [36].

Methods
In this section, first, the experimental setup and the experimental process is described.Then, a brief description of the equations which use for the calculations is presented.It is then followed by the numerical study and simulation process.

Experimental Study
For the experiment, initially, about 6 kg of biomass is fed into the gasifier along with 1 kg of charcoal.Charcoal is added to the system to help the ignition process.The simulations are also carried out following the same procedure.The average size of the Birchwood chip which was used for the experiments is approximately 11.5 mm × 11.5 mm.A schematic diagram of the system is shown in Figure 1.The hopper is a conically shaped stainless-steel structure with a volume of 0.13 m 3 and doublewalled thick walls.It has a radius of 0.5 m.Biomass is fed through this hopper, and it is connected to an electrical shaker which can provide 15 shakes per minute (shakes/min).This shaking process helps to reduce or avoid some common problems, such as channeling and bridging.To drain out the condensed moisture, a drain valve is also attached to the bottom of the hopper.
The primary gasification reactor has a diameter of 0.26 m.It is 0.6 m in height.The air inlet system consists of six same size nozzles with equal spacing among each of a 5 mm diameter, and they are located from about 0.4 m from the bottom of the reactor.At the bottom of the reactor part, a reciprocally oscillating grate is connected and takes out the ash to the lower part of the ash collector.The grate is oscillating 30 s/min.To measure the weight loss during the gasification process, a mass scale locates at the bottom of the gasifier system.In addition, two pressure sensors and five K-type thermocouples are connected to the gasifier assembly to monitor and measure the parameters.
All the necessary experimental data is measured via a Lab view program which is specially designed to obtain data from the gasifier system.Further, to measure the composition of the syngas, a gas sample is sent to a gas analyzer, and throughout the experiment, changes in the composition of the syngas output is produced by the gas analyzer.

Theoretical and Numerical Study
In this study, the effect of biomass particle size and ER on the behavior of gasification parameters is examined using Birchwood as the feedstock.Thus, the variance of producer gas composition, producer gas yield, LHV of the syngas, and CGE is mainly assessed on two different values of biomass particle sizes.
Although downdraft gasifier demonstrates better performances at the ER ratio of 0.25, acceptable performances showed even close to the ER value of 0.43 [33,37].Therefore, ER value range from 0.2 to 0.5 is selected in this study to investigate the best performance.ER is calculated using the following expression Equation (1) [8,15].The hopper is a conically shaped stainless-steel structure with a volume of 0.13 m 3 and double-walled thick walls.It has a radius of 0.5 m.Biomass is fed through this hopper, and it is connected to an electrical shaker which can provide 15 shakes per minute (shakes/min).This shaking process helps to reduce or avoid some common problems, such as channeling and bridging.To drain out the condensed moisture, a drain valve is also attached to the bottom of the hopper.
The primary gasification reactor has a diameter of 0.26 m.It is 0.6 m in height.The air inlet system consists of six same size nozzles with equal spacing among each of a 5 mm diameter, and they are located from about 0.4 m from the bottom of the reactor.At the bottom of the reactor part, a reciprocally oscillating grate is connected and takes out the ash to the lower part of the ash collector.The grate is oscillating 30 s/min.To measure the weight loss during the gasification process, a mass scale locates at the bottom of the gasifier system.In addition, two pressure sensors and five K-type thermocouples are connected to the gasifier assembly to monitor and measure the parameters.
All the necessary experimental data is measured via a Lab view program which is specially designed to obtain data from the gasifier system.Further, to measure the composition of the syngas, a gas sample is sent to a gas analyzer, and throughout the experiment, changes in the composition of the syngas output is produced by the gas analyzer.

Theoretical and Numerical Study
In this study, the effect of biomass particle size and ER on the behavior of gasification parameters is examined using Birchwood as the feedstock.Thus, the variance of producer gas composition, producer gas yield, LHV of the syngas, and CGE is mainly assessed on two different values of biomass particle sizes.
Although downdraft gasifier demonstrates better performances at the ER ratio of 0.25, acceptable performances showed even close to the ER value of 0.43 [33,37].Therefore, ER value range from 0.2 to 0.5 is selected in this study to investigate the best performance.ER is calculated using the following expression Equation (1) [8,15].

ER =
Actual air fuel ratio Stoichiometric air fuel ratio (1) here, the actual air-fuel ratio is calculated from the actual values measured from the experiment, while the stoichiometric air-fuel ratio is determined using Equation (2).
where a and b can be determined using the results of the ultimate analysis of the biomass (Birchwood).
Here, y and z represent the stoichiometric coefficients of H 2 O and CO 2 , respectively, while x represent the stoichiometric coefficient of O 2 and 3.25 times of the N 2 amount.Further, the other parameters are calculated using the following Equation (3) to Equation ( 5) [33,[38][39][40].
LHV syngas = ((25.7 Since the Hydro-Carbons (HC) higher than CH 4 are not considered for this study, they have not been used to calculate LHV syngas .
Producer gas yield Cold gas efficiency (CGE) = LHV syngas × V syngas LHV biomass × 100 (5) where LHV syngas represents the dry base lower heating value of the syngas in MJ/Nm 3 , V syngas represents the syngas yield Nm 3 /kg feed stock, Q air represents the air input in Nm 3 /h, m b represents the biomass input in kg/h, LHV biomass represents the lower heating value of the feedstock in MJ/kg, and N 2 represents nitrogen mass fraction of the output gas.
In the simulation model, a spherical particle is used for modeling.Hence, spherical particles with a diameter of 11.5 mm (available wood chip size in the lab) and 9.18 mm are utilized for the simulations (a spherical particle which has the same surface area as an actual wood chip, with a diameter of 9.18 mm consider here).
In the presented work, the STAR CCM+ software package is used for the modeling and simulation.A 3D geometric model of the actual gasifier is designed to create the computational model.For the gasifier model, 0.01m mesh is used, while in the air inlets and air-fuel mixing region, a 0.001 m mesh has been used.Figure 2 shows the customized meshfor the gasifier (0.01 m) and for air inlets (0.001 m).

ER = Actual air fuel ratio Stoichiometric air fuel ratio
(1) here, the actual air-fuel ratio is calculated from the actual values measured from the experiment, while the stoichiometric air-fuel ratio is determined using Equation (2).
where  and  can be determined using the results of the ultimate analysis of the biomass (Birchwood).Here, y and z represent the stoichiometric coefficients of H 2 O and CO 2 , respectively, while x represent the stoichiometric coefficient of O 2 and 3.25 times of the N 2 amount.Further, the other parameters are calculated using the following Equation (3) to Equation ( 5) [33,[38][39][40].
Since the Hydro-Carbons (HC) higher than CH4 are not considered for this study, they have not been used to calculate   .

Producer gas yield (𝑉
Cold gas efficiency (CGE) =   ×     × 100 (5) where   represents the dry base lower heating value of the syngas in MJ Nm 3 ⁄ ,   represents the syngas yield Nm 3 kg feed stock ⁄ ,   represents the air input in Nm 3 h ⁄ ,   represents the biomass input in kg h ⁄ ,   represents the lower heating value of the feedstock in MJ kg ⁄ , and N2 represents nitrogen mass fraction of the output gas.
In the simulation model, a spherical particle is used for modeling.Hence, spherical particles with a diameter of 11.5 mm (available wood chip size in the lab) and 9.18 mm are utilized for the simulations (a spherical particle which has the same surface area as an actual wood chip, with a diameter of 9.18 mm consider here).
In the presented work, the STAR CCM+ software package is used for the modeling and simulation.A 3D geometric model of the actual gasifier is designed to create the computational model.For the gasifier model, 0.01m mesh is used, while in the air inlets and air-fuel mixing region, a 0.001 m mesh has been used.Figure 2 shows the customized meshfor the gasifier (0.01 m) and for air inlets (0.001 m).

Figure 2. Meshed gasifier geometry and customized mesh at air inlets.
A segregated flow model is implemented to solve the transport equations.In addition, the Lagrangian particle modeling method is used to model the particle behavior (solid phase).A segregated flow model is implemented to solve the transport equations.In addition, the Lagrangian particle modeling method is used to model the particle behavior (solid phase).

Assumptions
The underlying assumptions of the modeling process are summarized as follows:

•
Steady flow is considered inside the gasifier.

•
The flow inside the domain is considered as incompressible and turbulent.

•
Spherical particles are used.

•
Evenly distributed particle regime is utilized.

•
The No-slip condition is imposed on inside wall surfaces.

Eulerian-Lagrangian Method
Reynolds averaged Navier-Stokes (RANS) equations are solved using a Eulerian-Lagrangian reference frame to solve the numerical scheme.The gas phase is considered as a reacting gas phase, and Lagrangian particle models are used to address the solid (particle) phase.Standard k-ε model is used for turbulence modeling.In addition, both finite rates of chemical reactions and eddy dissipation rates are used in the computational model.Moreover, for the modeling work, the conservation equations for transport of energy, momentum, and mass are used (Appendix A).Also, the default sub models available for moisture evaporation, devolatilization and char oxidation are employed (Appendix A).In the Eulerian-Eulerian approach, gas phase, as well as the solid phase, is considered as continuum [41].In the Lagrangian approach, usually a large number of particles are tracked transiently.The method is started by solving the transient momentum equation for each particle.In order to calculate the trajectory of a particle, Newton's second law is used and written in a Lagrangian reference frame [41][42][43][44].Momentum balance for particles is described in Equation ( 6), where V p represents the velocity of the particle, m p is the mass of the particle, and t is time.F s and F b represent the surface force and body force, respectively.When these forces are observed in depth, each of them is a composition of few forces available in the system.The following Equations ( 7) and ( 8) describe the two forces.
where F d , F p , and F vm represent the drag force on the particle, pressure force as a gradient in the static pressure, and the extra virtual mass added while the acceleration in the phase, respectively.F g and F u represent the gravitational force and the user defined force added to the system, respectively.In addition, the drag force can be defined as follows in Equation ( 9).
where C D is the drag coefficient, ρ is the density, A p represents the surface area of the particle, and V s is the slip velocity of the particle.
Using the Nusselt number (Nu p ), the heat transfer coefficient can be found for the Lagrangian particle, which is the energy model for the Lagrangian particle.This is shown in Equation (10), where k has represented the thermal conductivity of the phase and D p is the particle diameter.Nusselt number (Nu p ) is expressed in the following Equation (11), where it is calculated by the Ranz-Marshall correlation and presented below in Equation (12).
p Pr 1 3 (11) where Re represents the Reynolds number, while the expression Pr expresses the Prandtl number [45,46].
where h is the heat coefficient, d p is the particle diameter, k represents the thermal conductivity of the gas, Re is the Reynolds number, and Pr represents the prandtl number.The letters a, c, m, and n are numerical constants which are determined by the flow field.According to the Ranz-Marshall correlation, at a fluid velocity (V) equal to zero, the heat transfer can only be affected by conduction.therefore, the values of a, c, m, and n become 2, 0.6, 0.5, and 0.33, respectively [47].

Reactions Modeling
In this work, the following approach is used to simulate the gasification process.This method has been used by [8,48].In this study, it is assumed that the moisture is evaporated completely in the drying phase [49].
In the pyrolysis phase, the Biomass (dry) is converted into char, tar, ash, and volatile matter [49,50].
Biomass (dry) where, α is the stoichiometric coefficient of tar.
Here, the volatile matter is described as follows.It is assumed that the volatile matter consists of only CO, CO 2 , H 2 , and CH 4 .
The following method is used to calculate the mass distribution in volatile matter and mass fractions in dry biomass.For the calculations, the data available in Table 1 is used.These values can be used to combine with the ultimate analysis.In this method, a set of linear equations is used to solve the C, O, and H components balance in the system [8,48] In the above equations, the coefficients are the ratios of molar mass of each atom to the molar mass of each species.(As an example, in Equation ( 15), 4  16 ).Here, with the molar masses, the numerical values are rounded up for the simplification of the calculations.For finding the mass fractions Y C,vol , Y O,vol , and Y H,vol , the following method is used along with the data from the ultimate analysis of the biomass.
According to the mass conservation, total mass before devolatilization − mass of fixed carbon = total mass after devolatilization Thus, according to the above criteria, The ratio between CO and CO 2 is temperature dependent.In this study, the ratio has been determined as 2.43 according to [8].For the calculations, the bed temperature is used as 907 K, which was obtained during the experiments.Below Equation ( 21) is used to determine the connection between CO and CO 2 at 907 K.
x y = 2.43 (22) The linear equations from Equation (9) to Equation (15) are solved using the MATLAB program, and the following values are obtained from the program.
x = 0.5605 According to the above results, 1 kg of volatile matter produce devolatilized products according to the following expression.
Using the following molar mass values, the above Equation ( 23) can be used as a molar based equation.
Therefore, according to the following (R4), the devolatilization of the volatiles is modeled.Based on the above devolatilization reaction (R4), the gasification process is modeled.In addition, for the modeling work, the following homogeneous and heterogeneous reactions are considered [8,28,51].
Water-gas shift reaction:

Results and Discussions
In this section, the computational model validation is conducted by comparing the available experimental data and numerical simulation results.Moreover, the comparison of the two sizes of biomass particle sizes is carried out using the gasification characteristics.

Development of Mass Fractions of Considered Gases in the Syngas
A specially designed lab view program was used to measure all the gas component mass fractions throughout the experiment.From the obtained data, the variation of the mass fractions of the gas components during the experiment is presented below in Figure 3.For the tests, a biomass woodchip with an average size of 11.5 mm × 11.5 mm is used, which is available in the laboratory.
In the simulation model, for two iterations, one data point has been recorded.600 iterations are used in the simulations, and 300 data points are recorded throughout the simulation.Using the obtained data points gas component mass fractions are graphed in Figure 4.For the simulations, two types of wood particle sizes are used.One of them is 11.5 mm in diameter, while the other is 9.18 mm in diameter.Below, Figure 4 was obtained using 11.5 mm diameter wood particle.
biomass particle sizes is carried out using the gasification characteristics.

Development of Mass Fractions of Considered Gases in the Syngas
A specially designed lab view program was used to measure all the gas component mass fractions throughout the experiment.From the obtained data, the variation of the mass fractions of the gas components during the experiment is presented below in Figure 3.For the tests, a biomass woodchip with an average size of 11.5 mm × 11.5 mm is used, which is available in the laboratory.By comparing Figures 3 and 4, it shows the fluctuation in gas components mass friction in the simulation values, whereas in the experiment, it shows smooth variation.It can be due to the fact that the interval between two data points are bigger.However, an average value is used to plot graphs.
The data obtained from the model is also compared with some of the data collected from the literature (Study 1 [52], and Study 2 [53]), as well as with the data achieved during the experiment of the present work.Figure 5 shows the comparison of the predicted data along with the experimental data.According to the graph, predicted data shows good agreement with the literature and experimental data.

Effect of ER and Biomass Particle Size on Quality of Producer Gas
The variation of gas components in producer gas in both experiments and simulation models is illustrated in the following Figure 6.For the simulations, a biomass particle with a diameter of 11.5 mm and a biomass particle with a diameter of 9.18 mm were used.(The change of gas component contours in various ER values are illustrated in Appendix A (Figures A1-A3)).By comparing Figures 3 and 4, it shows the fluctuation in gas components mass friction in the simulation values, whereas in the experiment, it shows smooth variation.It can be due to the fact that the interval between two data points are bigger.However, an average value is used to plot graphs.
The data obtained from the model is also compared with some of the data collected from the literature (Study 1 [52], and Study 2 [53]), as well as with the data achieved during the experiment of the present work.Figure 5 shows the comparison of the predicted data along with the experimental data.According to the graph, predicted data shows good agreement with the literature and experimental data.By comparing Figures 3 and 4, it shows the fluctuation in gas components mass friction in the simulation values, whereas in the experiment, it shows smooth variation.It can be due to the fact that the interval between two data points are bigger.However, an average value is used to plot graphs.
The data obtained from the model is also compared with some of the data collected from the literature (Study 1 [52], and Study 2 [53]), as well as with the data achieved during the experiment of the present work.Figure 5 shows the comparison of the predicted data along with the experimental data.According to the graph, predicted data shows good agreement with the literature and experimental data.

Effect of ER and Biomass Particle Size on Quality of Producer Gas
The variation of gas components in producer gas in both experiments and simulation models is illustrated in the following Figure 6.For the simulations, a biomass particle with a diameter of 11.5 mm and a biomass particle with a diameter of 9.18 mm were used.(The change of gas component

Effect of ER and Biomass Particle Size on Quality of Producer Gas
The variation of gas components in producer gas in both experiments and simulation models is illustrated in the following Figure 6.For the simulations, a biomass particle with a diameter of 11.5 mm and a biomass particle with a diameter of 9.18 mm were used.(The change of gas component contours in various ER values are illustrated in Appendix A (Figures A1-A3)).According to Figure 6, mass fractions of CO and CO2 are increased when the ER is increased up to a certain level (ER = 0.2 to 0.35).Then, ER contributes to a reduction of the mass fraction of CO, while further increasing CO2.This can be a result of the combined effect of Combustion reaction (R9), Boudouard reaction (R10), and Carbon partial reaction (R13).Therefore, in lower ER values, (R10) tends to produce more CO to the system at the expense of CO2, while in high ER values, (R9) and (R13) convert more CO into CO2 [2].
Rgearding CH4 and H2, increasing ER is effected on a decrease in both gas species.Nevertheless, H2 is varied in a higher mass fraction range than of CH4.This is mainly due to the effects of dry reforming (R7) and steam reforming (R8) reactions, which produce more H2 to the system at the expense of CH4.
Since a variable air inflow rate is used in the experiments, an average value of the air inflow rate is used for the calculation of ER values.However, a fixed air inflow rate is used in the simulations.Therefore, the impact of variable air inflow may have caused the deviation in the graphs of the experimental and simulation data.In the simulations, spherical shape wood particles are used, while wood chips with different shapes are used for the experiments.The difference of behavior of spherical wood particle and divergement shaped wood chips can also create the inequalities in the graphs.In addition, during the experiments, channeling and bridging can also have taken place [7].The effect of that phenomenon may have also affected the experimental results.
According to Figure 6, all the considered gas component mass fractions, except H2, are increased for the 9.18 mm wood particle over 11.5 mm wood particle.This pattern can be found in some of the literature [54,55].When the particle size is getting smaller, the area to volume ratio of the wood According to Figure 6, mass fractions of CO and CO 2 are increased when the ER is increased up to a certain level (ER = 0.2 to 0.35).Then, ER contributes to a reduction of the mass fraction of CO, while further increasing CO 2 .This can be a result of the combined effect of Combustion reaction (R9), Boudouard reaction (R10), and Carbon partial reaction (R13).Therefore, in lower ER values, (R10) tends to produce more CO to the system at the expense of CO 2 , while in high ER values, (R9) and (R13) convert more CO into CO 2 [2].
Rgearding CH 4 and H 2 , increasing ER is effected on a decrease in both gas species.Nevertheless, H 2 is varied in a higher mass fraction range than of CH 4 .This is mainly due to the effects of dry reforming (R7) and steam reforming (R8) reactions, which produce more H 2 to the system at the expense of CH 4 .
Since a variable air inflow rate is used in the experiments, an average value of the air inflow rate is used for the calculation of ER values.However, a fixed air inflow rate is used in the simulations.Therefore, the impact of variable air inflow may have caused the deviation in the graphs of the experimental and simulation data.In the simulations, spherical shape wood particles are used, while wood chips with different shapes are used for the experiments.The difference of behavior of spherical wood particle and divergement shaped wood chips can also create the inequalities in the graphs.In addition, during the experiments, channeling and bridging can also have taken place [7].The effect of that phenomenon may have also affected the experimental results.
According to Figure 6, all the considered gas component mass fractions, except H 2 , are increased for the 9.18 mm wood particle over 11.5 mm wood particle.This pattern can be found in some of the literature [54,55].When the particle size is getting smaller, the area to volume ratio of the wood particle is increased.This can result in an increased contact with the reactants [54][55][56].Hence, the reaction rates can be increased and gasification reactions can also be improved.Further, the enhancement of hydrogasification (R12) is crucial for the behavior of the gas component mass fractions.When the biomass particle is getting smaller, it can enhance the reactivity and produce more CH 4 with R12. Figure 7 shows the variation of syngas yield, LHV and CGE of the syngas as a function of ER.According to Figure 7, the behavior of syngas yield could be an effect of higher volatilization [28].
Energies 2017, 10, 1232 12 of 20 particle is increased.This can result in an increased contact with the reactants [54][55][56].Hence, the reaction rates can be increased and gasification reactions can also be improved.Further, the enhancement of hydrogasification (R12) is crucial for the behavior of the gas component mass fractions.When the biomass particle is getting smaller, it can enhance the reactivity and produce more CH4 with R12.

Effect of Equivalences Ratio (ER) and Biomass Particle Size on Syngas Yield, Lower Heating Value (LHV), and Cold Gas Efficiency (CGE)
Figure 7 shows the variation of syngas yield, LHV and CGE of the syngas as a function of ER.According to Figure 7, the behavior of syngas yield could be an effect of higher volatilization [28].During the simulations, both biomass particles are evaluated in the same ER range.However, the gas yields with two separate sizes of biomass particles have shown a difference.This difference of syngas yield should be due to a different conversion.For the context of LHV of syngas, increasing ER has shown a positive effect within the range from 0.2 to 0.35.It is evident that, as the gas species mass fractions increase, LHV of syngas increases as well.
According to Figure 7, the smaller particle size has resulted in a significant increment of the LHV of the gas, although it has followed the same variation along the ER range as the bigger biomass particle.The increment of LHV can be mainly due to the increment of CH4 with the smaller biomass particle.As mentioned before, LHV is a direct function of combustible gas species, and CH4 has the highest contribution for the LHV among the gas species.
For CGE, a positive effect of increasing ER is observed from 0.2 to 0.35 of ER.Within this ER range, CGE could have been improved by increasing airflow.Increasing airflow (increasing ER) can During the simulations, both biomass particles are evaluated in the same ER range.However, the gas yields with two separate sizes of biomass particles have shown a difference.This difference of syngas yield should be due to a different conversion.For the context of LHV of syngas, increasing ER has shown a positive effect within the range from 0.2 to 0.35.It is evident that, as the gas species mass fractions increase, LHV of syngas increases as well.
According to Figure 7, the smaller particle size has resulted in a significant increment of the LHV of the gas, although it has followed the same variation along the ER range as the bigger biomass particle.The increment of LHV can be mainly due to the increment of CH 4 with the smaller biomass particle.As mentioned before, LHV is a direct function of combustible gas species, and CH 4 has the highest contribution for the LHV among the gas species.
For CGE, a positive effect of increasing ER is observed from 0.2 to 0.35 of ER.Within this ER range, CGE could have been improved by increasing airflow.Increasing airflow (increasing ER) can contribute to creating a high composition of combustible gas species in syngas by enhancing the reaction chemistry of homogeneous and heterogeneous reactions.Therefore, as the gas species are increased up to 0.35 of ER, so does the CGE.As the gas species have shown a decrement in the mass fractions in the ER range from 0.35 to 0.5, CGE has also shown a decrement in its value.With more O 2 input to the system, the system can move towards combustion rather than gasification.This can be a result of the reduction of CGE after the ER value of 0.35.According to the graph in Figure 7, with the 9.18 mm biomass particle, CGE has shown an average maximum value of 60% at ER value of 0.35, compared to an average maximum value of 56.63 % with 11.5 mm diameter biomass particle.
For both particle sizes, same ER range is considered.Hence, an increase of LHV, syngas yield, and combustible gas spices in syngas with smaller biomass particle could be the reason for the increase of the CGE.

Temperature Profile
Temperature profile inside the gasifier in the simulation model is illustrated in Figure 8.In Figure 9, the temperature values recorded in the experiments and the simulation model are demonstrated in a graph.Figure 8 is mainly used to help understand the actual placement of the thermocouples.In the simulation model, seven temperature measurements are taken along the mid-axis, which is shown in Figure 9.However, in the experimental setup, only four thermocouples are available.Hence, to develop the curve, those four temperature measurement values are used.According to Figure 9, the maximum temperature is shown close to the air inlets in both experimental and simulation setups.In the experimental setup, the thermo couples are placed on the gasifier wall.Hence, the actual temperature can be slightly higher than what it is shown here.
Energies 2017, 10, 1232 13 of 20 contribute to creating a high composition of combustible gas species in syngas by enhancing the reaction chemistry of homogeneous and heterogeneous reactions.Therefore, as the gas species are increased up to 0.35 of ER, so does the CGE.As the gas species have shown a decrement in the mass fractions in the ER range from 0.35 to 0.5, CGE has also shown a decrement in its value.With more O2 input to the system, the system can move towards combustion rather than gasification.This can be a result of the reduction of CGE after the ER value of 0.35.According to the graph in Figure 7, with the 9.18 mm biomass particle, CGE has shown an average maximum value of 60% at ER value of 0.35, compared to an average maximum value of 56.63 % with 11.5 mm diameter biomass particle.
For both particle sizes, same ER range is considered.Hence, an increase of LHV, syngas yield, and combustible gas spices in syngas with smaller biomass particle could be the reason for the increase of the CGE.

Temperature Profile
Temperature profile inside the gasifier in the simulation model is illustrated in Figure 8.In Figure 9, the temperature values recorded in the experiments and the simulation model are demonstrated in a graph.Figure 8 is mainly used to help understand the actual placement of the thermocouples.In the simulation model, seven temperature measurements are taken along the midaxis, which is shown in Figure 9.However, in the experimental setup, only four thermocouples are available.Hence, to develop the curve, those four temperature measurement values are used.According to Figure 9, the maximum temperature is shown close to the air inlets in both experimental and simulation setups.In the experimental setup, the thermo couples are placed on the gasifier wall.Hence, the actual temperature can be slightly higher than what it is shown here.

Conclusions
Based on the comparative study, it is concluded that most of the measured and calculated parameters (including syngas yield, LHV of syngas, and CGE) have shown an increasing trend as ER is increased up to 0.35, and when biomass particle size is reduced.Furthermore, the average value of CGE, for instance, reached a maximum value of 60% with the 9.18 mm diameter wood particle, and with a corresponding ER about 0.35.Moreover, 11.5 mm diameter wood particle has shown an average maximum value of 56.6% at the same ER value, which has been both an improvement and a significant output of this work.Reaching an average value of 2.79 Nm 3 /h with a corresponding ER of 0.35, the average value of Syngas yield has shown an improvement of 0.09 Nm 3 /h with the 9.18 mm diameter wood particle over 11.5 mm wood particle, which has been an important aspect in energy harvesting from the syngas.Therefore, as the result shows, an ER value close to 0.35 shows the best results in CGE, syngas yield, and LHV of the syngas.Furthermore, smaller biomass particles show improvements in the results.Therefore, according to our study, in order to obtain better outputs, ER should be kept close to 0.35.In the future, more simulations and experiments could be done with finer biomass particles.Further, more user defined functions could be used to improve the gasification process and to represent more realistic biomass particle shapes.

Figure 2 .
Figure 2. Meshed gasifier geometry and customized mesh at air inlets.

Figure 3 .Figure 3 .
Figure 3. Mass fraction variation with time in the experimental study.In the simulation model, for two iterations, one data point has been recorded.600 iterations are used in the simulations, and 300 data points are recorded throughout the simulation.Using the

Figure 4 .
Figure 4. Mass fraction variation with iteration in the simulation model.

Figure 5 .
Figure 5. Composition of the syngas predicted by model compared with experimental results.

Figure 4 .
Figure 4. Mass fraction variation with iteration in the simulation model.

Energies 2017 ,
10, 1232 10 of 20obtained data points gas component mass fractions are graphed in Figure4.For the simulations, two types of wood particle sizes are used.One of them is 11.5 mm in diameter, while the other is 9.18 mm in diameter.Below, Figure4was obtained using 11.5 mm diameter wood particle.

Figure 4 .
Figure 4. Mass fraction variation with iteration in the simulation model.

Figure 5 .
Figure 5. Composition of the syngas predicted by model compared with experimental results.

4 Figure 5 .
Figure 5. Composition of the syngas predicted by model compared with experimental results.

Energies 2017, 10 , 1232 11 of 20 Figure 6 .
Figure 6.Variation of gas component mass fractions as a function of ER.

Figure 6 .
Figure 6.Variation of gas component mass fractions as a function of ER.

4. 3 .
Effect of Equivalences Ratio (ER) and Biomass Particle Size on Syngas Yield, Lower Heating Value (LHV), and Cold Gas Efficiency (CGE)

Figure 1 .
Figure 1.Variation of Syngas yield lower heating value (LHV) and cold gas efficiency (CGE) of syngas as a function of Equivalences Ratio (ER).

Figure 8 .
Figure 8. Temperature profile inside the gasifier.Figure 8. Temperature profile inside the gasifier.

Figure 8 .
Figure 8. Temperature profile inside the gasifier.Figure 8. Temperature profile inside the gasifier.

Figure 9 .
Figure 9. Temperature variation along the mid-axis of the gasifier simulation values vs. experimental values.

Physical or Chemical Property Value
. and w are unknown mass fractions of CO, CO 2 , CH 4 , and H 2 , respectively in volatile species after devolatilization.Y C,vol , Y O,vol , and Y H,vol are, respectively, the mass fractions of C, O, and H in biomass volatile.