Simulation of Olive Kernel Gasification in a Bubbling Fluidized Bed Pilot Scale Reactor

The main aim of this study is to develop a comprehensive process model for biomass gasification in a pilot scale bubbling fluidized bed gasifier using the ASPEN PLUS simulator. A drawback in using ASPEN PLUS is the lack of a library model to simulate fluidized bed unit operation. However, it is possible for users to input their own models, using FORTRAN codes nested within the ASPEN PLUS input file, to simulate operation of a fluidized bed. The products of homogeneous reactions are defined by Gibbs equilibrium and reaction rate kinetics are used to determine the products of char gasification. Governing hydrodynamic equations for a bubbling bed and kinetic expressions for the char combustion were adopted from the literature. Different sets of gasification results for the operation conditions (temperature and air equivalence ratio (ER)) obtained from the our pilot-scale gasifier having a capacity of 1 kg/hr of olive kernel as feeding biomass, were used to demonstrate the validation of the model. The simulation results received from the application of the model were compared with the above experimental results and showed good agreement.


Introduction
The traditional approach necessary to establish comercial plant technology is based on comprehensive experimental investigations, progressing from a laboratory scale test unit to a pilot scale plant, before building a full-scale commercial demonstration plant.For process optimisation, an extensive investigation of the plant behaviour depending on various operating parameters is required for each scale up step.To support this optimisation procedure, mathematical models are helpful to reduce the temporal and financial efforts.
Pre-condition is a reliable simulation tool, which includes the mathematical formulation of all important chemical and physical processes by describing their dependency on operating parameters and their interdependencies [1].However, only limited research has been performed to model those complex systems including also the hydrodynamics of biomass particles.
The development of numerical models for fluidizedbed gasification (FBG) documented in the literature [2] are devoted mainly to coal and less to biomass.Even though, biomass in comparison to coal is made up not only of lignin but also of cellulose and hemi cellulose, each one having its own thermal behaviour and making the modelling more difficult, the modelling approaches used for coal can serve for biomass gasification modelling as well.Ma et al. [3] developed a steady-state model which involved instantaneous devolatilization of coal at the top of the gasifier (freeboard region) and char combustion and gasification in the fluidized bed.Chejne and Hernandez [4] developed a comprehensive mathematical model to predict the behaviour of coal combustion and gasification on stacks in non-stationary operation.Sett and Bhattacharya [5] developed a mathematical model for the behaviour of a fluidised-bed charcoal gasifier, taking into account fluidised-bed hydrodynamic conditions, chemical reaction rates, mass diffusion rates, and the conservation of mass and energy.The model was solved numerically using an IBM 3083 mainframe computer.Ross et al. [6] modified their one-dimensional numerical isothermal model of a fluidised-bed coal gasifier in order to simulate the performance of a laboratory-scale gasifier.
A new numerical model based on the two-fluid model (TFM) including the kinetic theory of granular flow (KTGF) and complicated reactions has been developed to simulate coal gasification in a bubbling fluidized bed gasifier (BFBG) [7].The collision between particles was described by KTGF.The coal gasification rates were determined by combining Arrhenius rate and diffusion rate for heterogeneous reactions or turbulent mixing rate for homogeneous reactions.The flow behaviours of gas and solid phases in the bed and freeboard could be predicted, which were not easy to be measured through the experiments.Hamel and Krumm [2] have developed a mathematical model for simulation of gasification processes of solid fuels in atmospheric or pressurised bubbling fluidised beds incorporating bed and freeboard hydrodynamics, fuel drying and devolatilization, and chemical reaction kinetics is presented.The model has been used to simulate four bubbling fluidised bed gasifiers, described in literature, of different scales from atmospheric laboratory scale up to pressurised commercial scale, processing brown coal, peat and sawdust.The gasifiers have been operated within a wide range of parameters using air, air steam or oxygen steam as gasification agent, operating with or without recirculation of fines at operating pressures up to 2.5 MPa.The simulation results for overall carbon conversion, temperature and concentrations of gaseous species agreed sufficiently well with published experimental data.
Corella and Sanz [8,9], at the University of Zaragoza and Madrid (Spain) started to study the modelling of fluidized bed biomass gasifiers in the mid-1980s.More recently Corella et al. [10], discussed the reaction network existing in a CFB biomass gasifier and the problems associated with the accuracy of the kinetic equations needed for the existing complex reaction network.Furthermore, he presented a model for bubbling fluidized bed (BFB) biomass gasifiers, gasifying with pure steam [10].That model identified the four main, for modelling purposes, chemical reactions among the reaction network existing in the gasification process.With only four kinetic parameters, the model predicted quite well the BFB gasifier.More recently, Corella and Sanz [8,9] have presented a whole model for CFBBGs.Such model is 1dimensional for steady state.The model has a semi rigorous character because of the assumptions that had to be introduced by lack of accurate knowledge of some parts of the modelling.
De Souza-Santos [11] developed a comprehensive mathematical model and commercially available computer program performing a comprehensive simulation of fluidized-bed equipment (CSFB Version 3.5), to use as a tool for engineering design and operation optimisation, by predicting the behaviour of a real unit during steadystate operation.His model is regarded as complete and it includes the conservation equations for the emulsion phase and bubbles, empirical equations for hydrodynamics, and it is also includes a through mass balance which considers that both drying and volatilisation are not instantaneous.Jiang and Morey [12] developed one-dimensional, steady state, numerical model for a fluidized bed biomass gasifier.The gasifier model consisted of a fuel pyrolysis model, an oxidation model, a gasification model and a freeboard model which were validated with ex-perimental results.
Haggerty and Pulsifer [13] used the reaction model together with three different reactor models (a) plug-flow, (b) complete-mixing and (c) bubble-assemblage, where the bubble-assemblage model represents a valid alternative when modelling fluidized-bed gasifiers.Marias et al. [14] developed a mathematical model for the fluidizedbed incineration process using the waste composed of wood, cardboard and polyvinyl chloride.Aarsen's model treated an isothermal fluidized bed as two separate zones, the oxidation zone at the bottom and the gasification zone at the top [15].Fuel pyrolysis took place at such a fast rate that only pyrolysis yields, which were assumed to vary with bed temperature, were needed.Mukosei et al. [16] examined the problem of the mathematical simulation of heterogeneous processes in a fluidized-bed reactor.These works and other more recent ones such as the work by He and Rudolph [17] proposed a new approach to the modelling of gross gas-solids flow through the riser in a circulating fluidized-bed system.This approach differs from the previous ones, which are found to be theoretically incorrect based on a fundamental analysis of the riser process hydrodynamics.
Panopoulos et al. [18] used the simple way of approaching the biomass gasification modelling by predicting thermodynamic equilibrium composition through Gibbs free energy minimisation calculations for the C, H, and O atoms of the fuel and the gasification agent mixture.Samuilov et al. [19] developed a mathematical model for the gasification in carbon dioxide of a single carbon particle.The porous structure of the particle, diffusion processes, the gasification processes on the pore surface according to the Langmuir-Hinshelwood model, and reactions on active carbon sites were taken into account.Deen et al. [20] reviewed the use of discrete particle models (DPMs) for the study of the flow phenomena prevailing in fluidized beds.
The ASPEN PLUS process simulator has been used by different investigators to simulate coal conversion; examples include methanol synthesis [21,22], indirect coal liquefaction processes [23], integrated coal gasification combined cycle (IGCC) power plants [24], atmospheric fluidized bed combustor processes [25], compartmented fluidized bed coal gasifiers [26], coal hydrogasification processes [27], and coal gasification simulation [28].However, the work that has been done on biomass gasification is limited.Mansaray et al. [29] used ASPEN PLUS to simulate rice husk gasification based on material balance, energy balance, and chemical equilibrium relations.Because of the high amount of volatile material in biomass and the complexity of biomass reaction rate kinetics in fluidized beds, they ignored the char gasification and simulated the gasification process by the assumption that biomass gasification follows Gibbs equilibrium.In a typical atmospheric fluidized bed gasifier, feed, together with bed material, are fluidized by the gasifying agents, such as air and/or steam, entering at the bottom of the bed.The product gas resulting from the gasification process is fed to a gas-solid separator (i.e., cyclone) to separate solid particles carried by exhaust gas.

Modelling Approach
Mathematical modelling is an alternative means to study fluidisation and fluidised bed biomass gasification.Modelling can not only account for the fundamental hydrodynamic behaviour of fluidisation and the coal gasification process, but also serve as a predictive tool to assist with the design, optimisation and scale-up of fluidised bed gasifiers.Therefore, in recent years mathematical modelling of fluidised beds has increasingly attracted more attention.

Assumptions
The following assumptions were considered in modelling the gasification process:

Combustion Zone
The combustible substance of a solid fuel is usually composed of elements carbon, hydrogen and oxygen.In complete combustion carbon dioxide is obtained from carbon in fuel and water is obtained from the hydrogen, usually as steam.The combustion reaction is exothermic and yields a theoretical oxidation temperature of 1450˚C.The main reactions, therefore, are:

Reduction Zone
The products of partial combustion (water, carbon dioxide and uncombusted partially cracked pyrolysis products) now pass through a red-hot charcoal bed where the following reduction reactions take place.
Reactions ( 3) and ( 4) are main reduction reactions and being endothermic have the capability of reducing gas temperature.Consequently the temperatures in the reduction zone are normally 800˚C -1000˚C.Lower the reduction zone temperature (~ 700˚C -800˚C), lower is the calorific value of gas.

Pyrolysis Zone
Wood pyrolysis is an intricate process that is still not completely understood 14.The products depend upon temperature, pressure, residence time and heat losses.However following general remarks can be made about them.Up to the temperature of 200˚C only water is driven off.Between 200˚C to 280˚C carbon dioxide, acetic acid and water are given off.The real pyrolysis, which takes place between 280˚C to 500˚C, produces large quantities of tar and gases containing carbon dioxide.Besides light tars, some methyl alcohol is also formed.Between 500˚C to 700˚C the gas production is small and contains hydrogen.Thus it is easy to see that updraft gasifier will produce much more tar than downdraft one.In downdraft gasifier the tars have to go through combustion and reduction zone and are partially broken down.

Hydrodynamic Parameters
The following assumptions were made in simulating the hydrodynamics:

•
Fluidized bed reactors divided into two regions bed and freeboard

•
The fluidization state in the bed is maintained in the bubbling regime The volume fraction of solids decreases as height increases, corresponding to the coalescence of bubbles in the bed and the returning of solid particles to the bed in the TDH zone • Volumetric flow rate of gas increases along with height, corresponding to the production of gaseous products • The mixing of solid particles, consisting of ash, char particles, and bed material, is perfect The reactor is divided into a finite number of equal elements with constant hydrodynamic parameters • The fluidized bed is one-dimensional; any variations in conditions are considered to occur only in the axial direction.

Global Model of Fluidized Bed Gasifier
The proposed gasification dynamic model represents two stages: instantaneous devolatilization of straw and combustion of the char at the bottom of the gasifier and the gasification in the fluidized bed.Thus, a two-phase representation of the fluidized bed incorporates the phenomena of jetting, bubbling, slugging and dynamic mass and energy balances.Figure 1 shows the fluidized bed domain divided into three regions: jet, bubbles and slugs.
The following assumptions were made regarding to the gas low divisions and the low patterns of gas and solids: 1) the fluidized bed consists of a dilute phase (jets, bubbles and = or slugs) and an emulsion phase, 2) the emulsion phase is divided into an interstitial gas phase and a solid phase 3) mass and heat exchange take place between the dilute phase and interstitial gas and between the interstitial gas and the solids, 4) the fluidizing gas enters the bed through nozzles in a jet form.The jets degenerate into bubbles, which rise through the bed and grow by coalescence with other bubbles to form slugs, 5) slugging occurs if the bubble diameter becomes larger than one third of the reactor diameter, (Slugging occurs in improper fluidization 6) the dilute gas and the interstitial gas are in plug low and the dilute gas is free of solids, 7) the gas behaves ideally, and 8) the produced gas consist of CO, CO 2 , H 2 , H 2 O, CH 4 and N 2 .

Mass and Energy Balances
The olive kernel feed flow is calculated by the following equation [30]: For the calculation of the air that is demanded for the stoichiometric combustion of biomass is used [30]: where the coefficients a, b, c, d represent the molecular fractions H/C, O/C, N/C and S/C .Therefore we can calculate the air flow for several values of equivalence ratio (ER) according to previous equation [30] With the assumption that gasification char is consisted of fixed carbon and ash, is calculated the solid mass flow from the equation [30]: Then the overall mass balance is formed for the gasification reactor and the mass flow for the gas gasification is estimated as following [30]: Correspondingly the overall energy balance is described as below [30] Es represents the thermal content of biomass feed flow which is estimated from the following equation [30]: For the air it is assumed that enters the reactor at environmental temperature, so its enthalpy content is considered zero (Eair, gasif = 0) .
Eth represents the thermal energy and the enthalpy of gas product and it is given by the following equation [30]: where The following equation estimates the enthalpy of each component of gas gasification in gasification temperature Finally Eloss represents the energy losses from the system due to ash and char removal and the losses to environment from the walls of the gasifier.Therefore: It is assumed that the energy losses from the gasifier walls constitute the 10% of the overall losses.The energy losses due to char gasification removal is calculated as below [30]:

Air Equivalence Ratio
tes a significant factor of Equivalence ratio (ER) constitu biomass gasification.The wise choice of it discourages the propulsion of the oxidation reactions.It is defined by the following equation [31]: A π represents the mass ratio between air and bi th pr r Gas LHV nificantly affect the heating omass at has been used in laboratory and A θ represents the stoichiometrically demanded mass ratio of air for the total combustion of the same quantity of biomass.
The alteration of air flow effects the process since it leads it to the two terminal situations: 1) High air equivalence ratio (ER > 1) approaches the ocess of combustion giving as gasification product mainly CO 2 2) Low air equivalence ratio (ER < 0.25) benefits the process of pyrolysis giving primarily as product syngas (CO + H 2 ).

Produce
Temperature and ER also sig value of producer gas.Main research target is the production of a producer gas enriched in CO, H 2 , and CH 4 .The presence of these combustible species contributes to the production of a medium to high heating value gas, suitable for further exploitation in internal combustion engines (ICE) and turbines for power production.LHV calculation was made using the following equation [32]: where, CO, H 2 and CH 4 are the molar ratios of the specteristics used as biomass feed-heating or in t fluidized bed gasifier red for biogenic mah a GC system for quantitative dete The FBG unit is equipped with a feeding hopper suitable g takes place cies in the producer gas as measured in our gasification experiments.

Biomass Chara
In the laboratory olive kernel was ing and its characteristics are illustrated in the Table 1.
Olive kernel is the final solid residue emerging from olive oil production industries.Olives are processed in an oil extracting plant to recover the oil content.After a first residue drying, residual oil is hexane extracted generating a residual solid product with a moisture content of around 10 -12 wt% called olive kernel.Traditionally, such kind of residue is sold as fuel for combustion in small boilers and specially designed furnaces due to its significant calorific value (HHV ~21 MJ/kg).Nowadays, there are several boilers fuelled with olive kernels destined for residential and greenhouse some cases to generate steam for electricity production in conventional steam turbines.However, a more efficient and environmentally friendly alternative with respect to CO 2 and CH 4 mitigation to the atmosphere, is olive kernel gasification; aim is to produce a high calorific value gas to fuel a gas engine or even a gas turbine in a combined cycle of heat/power production.

The Bench Scale FBG Uni
A 5 kW , bench scale, bubbling th (FBG) was designed and manufactu terial gasification.It is an atmospheric air blown FBG unit that consists of four sections: 1) the biomass feeding section, 2) the air supply, control and preheating section, 3) the gasification facility, 4) the gas sampling and off line analysis section.
Producer gas was sampled using airtight gas sampling bags and analyzed wit rmination of CO, CO 2 , H 2 , CH 4 , C 2 H 4 and C 2 H 6 .

The Olive Kernel Feeding Section
for feeding of low density biomass.Feedin through a system of a constant speed screw feeder and a variable speed screw feeder (motorized via inverter), separated by a rotary valve in order to avoid the hot gas backflow that could pyrolyse olive kernels and cause serious tar fouling problems.Additionally, a small proportion of fluidization air (~1 l/min) is supplied for backflow prevention.The motorized screw feeder is able to control the biomass feeding, based on a series of calibration curves and the maximum mass flow rate of olive kernel could reach 120 g/min.Olive kernel pass through the metal to metal rotary valve and is introduced into the reactor by the second screw feeder, into the hot zone of the reactor, 120 mm above the air distributor.s it plays a dual role acting as heat transfer met an inert material, as an internal e Air Supply, Contr heatin ry's air supply ary and second tion and gasification process, always in hypo-stoichiometric amounts obtaining the desired ER ratio.It consists of a central compressor, a main valve to control the air flow rate and an air flowmeter with maximum flow capacity of 200 l/min.The primary air for gasification process is introduced from the bottom of the FBG and a uniformly perforated plate (151 holes of 1mm diameter) is used in order to retain bed fluidization.Primary air is preheated up to 350˚C, by passing through a stainless steel tube that is kept in touch with the hot zone of the reactor inside the furnace.In the beginning of each operation, while the reactor is still cold, air preheating is obtained via an electric resistance placed at the bottom of the FBG.Secondary air contributes in producer gas temperature increase, which results in further tar cracking improving carbon conversion efficiency.

The 5 kW th FBG Test Unit
Figure 2 shows the bench scale FBG reactor with nominal capacity of 5 kW th made of stainless steel with diamet mm respectively; the total reactor height is 1400 mm.A more detailed analysis about the reactor design and construction can be found in a previous published work [33].FBG design enables the easy dismantling and thorough cleaning, in order to avoid not only memory effects, but also possible tar fouling under the hot flow experimental conditions.The gas cleaning system consists of a cyclone of 10 μm cut point.Limited experimental work has done using high temperature ceramic filters and various types of metallic foam for particulate abatement and tar elimination.

Bed Material
The bed material is of great importance in a fluidised bed reactor a dium and, if it is no catalyst for extended tar cracking.
A major problem of gasification applicability is tar formation and its minimisation is seen as one of the greatest technical challenges to overcome for the successful development of commercially attractive gasification technologies.According to the literature, tar production in an air steam gasification process could be as low as ~20 g/Nm 3 [31].
Tar treatment technologies are categorised either as upstream, when the tar problem is tackled inside the gasifier (named also primary or internal methods) or downstream, when tar is cleaned in separate hot gas devices (named also secondary or external methods).Although the secondary methods are proven to be effective and fully controlled, primary treatment technologies are gain-ing more interest due to their economic benefits.In primary methods operating parameters such as temperature, gasifying agent, bed material, equivalence ratio and residence time can play an important role on the formation and decomposition of tars.Cao et al. [31] concluded that the utilization of some catalysts in the gasifier and the concepts of two stage gasification are of prime importance.They worked on a laboratory scale, two-region, combined fluidised bed reactor for air gasification of sawdust producing a H 2 rich gas, with an LHV of 5 MJ/Nm 3 and tar concentration below 10 mg/Nm 3 .It is stated that several in-gasifier materials have catalytic action and eliminate in-bed tar by activating reforming reactions.Such materials modify the kinetic constants of most catalytic reactions involved in gasification.The in-bed additives are classified into two categories: 1) laboratory scale catalysts, which are very active but expensive; those catalysts are deactivated by coking when used at high (above 1000 kg/hm 2 cross sectional area) biomass throughputs (realistic for commercial applications), 2) natural occurring solids such as dolomite, limestone, olivine, high iron-content solids, which have been proved to have only a relatively small tar-elimination activity.
Nordgreen et al. [34] used metallic iron and iron oxides to catalytically crack tars in a secondary reactor, downstream the atmospheric fluidized bed gasifier fuelled with Swedish birch.They worked at a temperature range of 700˚C -900˚C and low equivalence ratio almost near to pyrolysis conditions (0 < ER < 0.2), concluding that catalytic breakdown of the tar reached 100%.
The first set of experiments was carried out with quartz sand (Table 2) as bed material, in the size range: 500 -425 μm.Although quartz sand was an easy and cheap to find material, it caused severe defluidization due to inevitable agglomeration and tar formation tendency at relatively low gasification temperature (T < 800˚C).Thus, and taking into account that an iron containing bed material might act as a catalyst for tar destruction, quartz sand was replaced with olivine (Table 2).A 500 -425 μm size fraction of olivine was also used and it was found to be durable for a prolongued operating period and in temperature range 800˚C -950˚C.

Producer Gas Sampling and Analysis
The sampling line is kept at a ~300˚C in order to avoid tar condensation and line fouling.Tar condensation and particle removal takes place in the gas cleaning section (Figure 3) consisting of a water washing, moisture trap, impinger bottles with isopropanole, a fiber filter and a silica gel filter.
When the appropriate gasification conditions are achieved producer gas is sampled, using a membrane pum rtight gas sampling bags and analyzed in laboratory's GC system (Model 6890N, Agilent Technologies equipped with FID and TCD detectors connected in series).The Gas Chromatograph is fitted with two columns HP-Plot Q with helium as carrier gas.GC's temperature profile was an isothermal at 50˚C and the retention time of the analysis process was 38 min.The standard gas mixture used.

Simulation Software
ASPEN Plus was selected for modelling the gasifier.It is a steady state chemical process simulator, which was developed to evaluate synthetic fuel technologies.It uses unit operation blocks, which are models of specific process operations (reactors heaters, pumps etc.).The user places these blocks on a flowsheet, specifying material and energy streams.An extensive built in physical properties database is used for the simulation calculations.ASPEN Plus has the capability to incorporate Fortran code into the model.This feature is utilised for the definition of non-conventional fuels, e.g.biomass, specific coals and for ensuring the system operates within user defined limits and constraints.The development of a model in ASPEN Plus involves the follo lowsheet (unit operation blocks s).e, composition, ).
wing steps: 1) Stream class specification and property method selection.
3) Defining the process f and connecting material and energy stream 4) Specifying feed streams (flow rat and thermodynamic condition 5) Specifying unit operation blocks (thermodynamic condition, chemical reactions etc.).

Flow Sheet Development
ASPEN PLUS blocks that have been used in order to simulate the gasification process are described in Table 3. Figure 4 shows the flowsheet for the gasification process.

Drying
The ASPEN PLUS stoichiometric reactor, RSTOIC, was used to simulate the process of biomass drying.In this step biomass is separated from moisture.A SEPARA-TION COLUMN model was used before the RYIELD reactor to separate the moisture from dried biomass.Within the ASPEN PLUS environment, the separation column is the most appropriate unit operation to achieve this goal.The amount of volatile material can be specified from the biomass approximate analysis.

Biomass Decomposition
The ASPEN PLUS yield reactor, RYIELD, was used to simulate the decomposition of the feed.In this step, biomass is converted into its constituting components including carbon, hydrogen, oxygen, sulfur, nitrogen, and ash, by specifying the yield distribution according to the biomass ultimate analysis.

Volatile Reactions
The ASPEN PLUS Gibbs reactor, RGIBBS, was used for volatile combustion, in conformity with the assumption that volatile reactions follow the Gibbs equilibrium.Biomass consists of mainly C, H, N, O, S, and ash.Carbon will partly constitute the gas phase, which takes part in devolatilization, and the remaining carbon comprises part of the solid phase (char) and subsequently results in char gasification.A COMPONENT SEPARATOR model was used before the RGIBBS reactor to separate S and N within the ASPEN PLUS environment, the component separator is the most appropriate unit operation to achieve this goal.The amount of volatile material can be specified from the biomass approximate analysis.The MIXER block was used before RGIBBS so as to add air to gas mixture.

Char Gasification
The ASPEN PLUS PFR reactor, RPFR, performs char gasification by including gasification reactions and chemical kinetics, as mentioned previously.The hydrodynamic and kinetic parameters, such as superficial velocity and voidage are constant in these reactors.Also considering the assumption that char contains only carbon and ash, the amount of carbon in the volatile portion can be calculated by deducting the total amount of carbon in char from the total carbon in biomass.

Solid Separation
A CYCLONE SEPARATOR block was used in this model in order to distinguish carbon solid from the gas mixture that it is received as final product.A CYCLONE SEPARATOR block in ASPEN PLUS environm izing cyclones.
A model was developed for the gasification of biomass sifier using the ASPEN with air.The simulation results for the product gas composition versus temperature ere compared with experimenum conditions is shown in Tables 4 and 5. Carbon monoxide and methane show the highest and lowest ts of experiments.ent separates solids from gas using gas vortex in a cyclone and it is used for rating and s

Model Validation
In order to validate the simulation results, experimental data from gasification of olive kernel in a pilot-scale fluidized bed gasifier was used.Biomass produces more tar and unburned hydrocarbon in lower temperatures, which decreases hydrogen production.The error related to the prediction of hydrogen is the result of ignoring tar production in the simulation. in a bubbling fluidized bed ga PLUS simulator.To provide the model, several ASPEN PLUS unit operation blocks were combined and, where necessary, kinetic expressions and hydrodynamic models were developed using data and models from the literature.The model was used to predict the results of pilot-scale gasification of olive kernel and air equivalence ratio w tal results.
Sensitivity analysis is performed to monitor the dependence of different parameters on the composition of the produced gas from the reactor.In this work, the temperature and air equivalence ratio (ER) are analysed in order to evaluate their effects on the composition of the produced gas and on the LHV of gasification gas.The optimum operating conditions according to simulation results is ER = 0.2 and T = 750˚C.During the sensitivity analysis the model input data was kept the same as for model validation with one parameter being varied at any given time Simulation results were compared with all sets of experimental data.The analysis of data for product gases in the optim error, respectively, in all se

Modelling Predictions
Since the reliability of proposed model has been proved it ca gas composition and LHV for several values of temperature behaviour n be now used in order to predict gasification

H 2 /CO Ratio in Producer Gas
High yields of syngas (H 2 and CO) were obtained from olive kernel gasification with air and this behaviour could be attributed to the increased content of cellulose and hemicellulose in the raw material.Such information was known from previous experimentation in a laboratory scale, fixed bed gasification unit [35].Uses of syngas could be distinguished according to the H 2 /CO molar ratio.Producer gas with H 2 /CO ratio between 1 to 2 can be used for the production of 2 nd generation biofuels via Fischer Tropsch synthesis (H 2 /CO = 2) and high added value chemical products like methanol.The present study indicates that producer gas of olive kernels gasification at 750˚C and ER = 0.2 with a molar ratio of H 2 /CO = 1.75 could be useful for Fisher Tropsch    synthesis.The effect of temperature on H 2 /CO molar ratio

Conclusions
The objective of this study was to develop simulation capable of predicting the steady-state performance of an atmospheric fluidized bed gasifier by considering the hydrodynamic and reaction rate kinetics simultaneously.The products of homogeneous reactions were defined by Gibbs equilibrium, and reaction rate kinetics were used to determine the products of char gasification.
A BFB biomass gasifier model was developed using ASPEN Plus.The results obtained from the sensitivity analyses are in good agreement with experimental work asifier performance over a wide range of operating conditions.
The influence of ER and temperature were investigated and the main conclusion is that the optimum operating conditions are ER = 0.2 and temperature 750˚C something that agrees with experimental data.At these conditions the concentration of syngas (CO + H 2 ) in gasification gas was maximized (41.7 v/v%) and LHV took its highest value (6.77MJ/Nm 3 ).
By increasing the ER the concentrations of CO and H 2 were decreasing in the favour of CO 2 production.As a result of sults Based on the results obtained in the present study, the following conclusions are drawn: 1) The proposed combined transport and kinetic model is successfully validated with the experimental data reported in the literature.
2) The model predicted composition of producer gas matches very well with the experimental data.
3) The molar fractions of CO and CH 4 increase and that of water vapor decreases with time in the pyrolysis zone of the gasifier due to a subsequent increase in the temperature of pyrolysis zone.
4) The molar fraction of water vapor in the gaseou phase is the highest in the pyrolysis zone and decreases s e gasifier.
5) The molar fractions of CO and CH 4 in the gaseous phase increase in the pyrolysis zone as gas travels through it.However, the molar fractions of CO, H 2 and CH 4 decrease in the oxidation zone due to high temperature oxidation and hence the molar fraction of CO 2 increases.
Good qualitative agreement between model prediction and experimental data was achieved .However, to improve the simulation results, some modifications should be considered.Those modifications concern tar formation and mass transfer inside the particles: 1) In this study simulation results of parametric study of temperature and air ratio on gas composition (H 2 , CO, formation consideration in the ve the predicted results.For this rearate depe rticles.This is another feature that should be in -910(01)00356-4 .Therefore, the model is capable of predicting accurately g this the LHV of gasification gas was also decreasing which are in accordance with the experimental re s a th it travels through the oxidation and reduction zones of CO 2 , and CH 4 ) without considering tar formation were achieved.However, tar model might impro son, detailed experimental data about the influence of operating conditions on the formation of tar along with the kinetics studies is needed in order to obtain a thorough evaluation.
2) For the future experimental work it is proposed to use the chemical formula C x H y O z as that of tar.The parameters (x, y, z) are temperature and heating ndent.Tar production can be taking in account by defining non-equilibrium products in the RGIBBS.
3) Mass transfer inside solid particles is also an important parameter in gas-solid reactions and heat transfer inside particles.Heat transfer is also important between phases, pa cluded in the model in order to achieve better simulation prediction.Radial dispersion inside the reactor will also help to see wall effects on the hydrodynamics of the fluidized bed reactor.
by specifying reaction yields of each component.This model is useful when reaction stoichiometry and kinetics are unknown and yield distribution data or correlations are available RG simultaneous phase and chemical equilibrium by s is known.This model is ctor with specified rsion.Reactors where nown or unimportant but eaction are known.IBBS minimizing Gibbs free energy, subject to atom Models single-phase chemical equilibrium, or balance constraints.This model is useful when temperature and pressure are known and reaction stoichiometry is unknown RPFR Models a plug flow reactor.This model is useful when reaction kinetic useful when solids, such as char, are participating in the reactions Models stoichiometric rea RSTOIC reaction extent or conve reaction kinetics are unk stoichiometry and extent of r

Figure 4 .
Figure 4. Comprehensive simulation diagram for the fluidized bed gasification process.

5. 1 . 1 .
Effect of Temperature Figures 5-10 show the simulation results compared with experimental data for product gas composition versus three different temperatures in the range of 750˚C -9850˚C.
versus three different air equivalence ratios in the ra 2 -0.4 hown in Figures 11-16 .
air equivalence ratio.Figures 17-19 presents model predictions versus temperature and Figures 20-22 the predictions versus air equivalence ratio while Figure 23 shows LHV versus temperature and air equivalence ratio.

Figure 20 .
Figure 20.Effect of air ratio on gas composition (T = 750˚C).

Figure 21 .
Figure 21.Effect of air ratio on gas composition (T = 800˚C).

Figure 22 .
Figure 22.Effect of air ratio on gas composition (T = 850˚C).

Figure 23 .
Figure 23.Effect of air ratio and temperature on LHV.

Figure 24 .
Figure 24.Effect of temperature and ER on H /CO molar 2

Table 1 . Olive kernel characteristics.
. Bed and freeboard sections are ers of 60 mm and 90

Table co
GC calibration composed from CO, CO 2 , H 2 , CH 4 , C 2 H 4 and C 2 H 6 1 vv% of each one of the above components balanced in helium.
Figure 3. Producer gas sampling section.