gas if ication modelling

3
1st ERCOFTAC Conference on Simula tion of Mul tiphase Flows in Gasifica tion and Combustion , 18–21 September 2011, Dresden, Germany MODELLING OF BIOMASS PYROL YSIS AND GASIFICATION I N INDUSTRIAL-SCALE GASIFIER Kamil Kwiatkowski 12 , Jakub Korotko 3 , Pawe l  ˙ Zuk 1 , Konrad Bajer 12 1 Faculty of Physics, University of Warsaw 2 Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw 3 Faculty of Power and Aeronautical Engineering, Warsaw University of Technology email:  ka [email protected]. pl, korotkoj@gmail. com, [email protected] , konrad.bajer@f uw.edu.pl and info@biomassgasi fication.eu The solution to the problem of dealing with waste, both biomass and municipal, is its utilisation via the gasicatio n process. The direct product of biomas s gasicatio n, the bioma ss syngas, is consi dered as an attrac- tive, versatile and fully renewable source of energy. It contains mainly nitrogen, hydrogen and carbon monoxide and residual tars and other pollutant. Almost unlimited diversity in shapes and properties of wastes, complex homo- and heterogeneous chem- istry and complex multiphase ows makes gasication a dicult phenomenon for mathematical and numerical modelling. There are two main approaches to pyrolysis and gasication modelling, the rst focuses on the detail analysis of the ow, reactions and heat transfer occurring in one representative biomass particle, the second focuses on modelli ng the whole gasic ation chamber but use several simplic ation s and subscale models. Driv en by the immediate demand from industry we have been developing the latter approach, focusing on the microscopic phenomena mainly to parametrise them. In this paper we present the numerical model and its results for industrial-scale biomass gasier obtained par- allel from our ANSYS Fluent  biomass gasication  UDF-extention (more details in  [ 1])and in-house OpenFOAM porousBiomassGasicationFoam  solver (see also [2]). The model takes into cons idera tion air and syngas ow and reactions within porous layers, water evaporation from biomass, pyrolysis, gasication and heterogeneous combustion of carbon remnants  [ 3]. Figure 1: I ns ta ll at io n of wood chi ps ga si ca- tion and sawdut drying to pellets production, Szepi- etowo. Fig ure 2: Dry and wet wood ch ip s wait in g to be gasied. The com- po sition of the wood chips is complex, contains fr ag me nte d pl yw ood, chipboard, straw, boards, bark, sa wdust and wood shaving. The biomass bed formed in the gasication chamber is treated as a porous zone, where Darcy’s law is applicable. This approximation is sucient, since the Reynolds number is less than 300. Still both Fluent and OpenFoam solve full Navier-Stokes equation with the additional momentum source term. All parameters of the porous lay er are anisotropic, non-un iform and transient. Gene rally porosit y is increased in the pyrolysi s zone then the bed rapidly crash es and p orosi ty is decrease d. T o initialise the model we made observations in the real-life installation and developed supplementary zonal and one-dimensional model of gasier in Matlab (more details in  [4]). In the zonal appro ach each individual zone is art ic ial ly separa ted and modelled, while in the 1D and CFD approaches there are no articial division for zones, only the sequence of the process is clarify. In order to model the four mentioned above processes it is crucial to properly mode l the heat tra nsfer b et we en sol id (bi oma ss) and gas (air and synga s). We assumed that the heat exchange in the porous medium is purely convective and we the used well established heat transfe r parameters . Then we introduce the global heat transfer coecient based on experimental results. In order to make our approach versatile, and to deal with various kind of waste and biomass, we implemented the data structures containing the elemental com- position of solid, its moisture content, volatile particles, xed carbon and ash con- tent [5]. Additionally vol atiles are quantied by fract ional compositi on of lignin, cellulose and hemicellulose for vegetable biomass  [ 6]. We include the homogeneous chemical reactions of water-gas shift, combustion of CO, CH4 ,H2 and heterogeneous reactions such as gasication and char combus- tion. The reaction rates are dened by the chemical kinetics based on the Arrhenius equat ion (for homog eneou s react ions) and the Langmuir-Hinshe lwood (for hetero- gene ous ones). In the near future we inten d to repla ce the coeci ents tak en from literature [ 5]  by the data obtained from the thermogravimetric measurements that have already been scheduled. The results, which we present below are prepared for the industrial updraft, xed- bed gasier fuelled by wood chi ps. It is located in Szepieto wo and achiev es 3.5MW of heat power used for sawdust drying (g.  1  and  2). The samp le results obtained from both OpenFOAM and Fluent are presented in gures  3  and  4 .  The shape of the feed, marked in the gures in red, was assumed drawin g from the experience of the installation operator. This shape can v ary depending mainly on the sizes of the chips. Our model will be experimentally veried for this installation. 1

Upload: jean-carlos-gonzalez-hernandez

Post on 04-Jun-2018

216 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Gas if Ication Modelling

8/13/2019 Gas if Ication Modelling

http://slidepdf.com/reader/full/gas-if-ication-modelling 1/2

1st ERCOFTAC Conference onSimulation of Multiphase Flows in Gasification and Combustion,18–21 September 2011, Dresden, Germany

MODELLING OF BIOMASS PYROLYSIS AND GASIFICATION IN

INDUSTRIAL-SCALE GASIFIER

Kamil Kwiatkowski12, Jakub Korotko3, Pawel  Zuk1, Konrad Bajer12

1Faculty of Physics, University of Warsaw 2

Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw 3Faculty of Power and Aeronautical Engineering, Warsaw University of Technology email:   [email protected], [email protected], [email protected], [email protected]

and [email protected]

The solution to the problem of dealing with waste, both biomass and municipal, is its utilisation via thegasification process. The direct product of biomass gasification, the biomass syngas, is considered as an attrac-tive, versatile and fully renewable source of energy. It contains mainly nitrogen, hydrogen and carbon monoxideand residual tars and other pollutant.

Almost unlimited diversity in shapes and properties of wastes, complex homo- and heterogeneous chem-istry and complex multiphase flows makes gasification a difficult phenomenon for mathematical and numericalmodelling.

There are two main approaches to pyrolysis and gasification modelling, the first focuses on the detail analysis

of the flow, reactions and heat transfer occurring in one representative biomass particle, the second focuses onmodelling the whole gasification chamber but use several simplifications and subscale models. Driven by theimmediate demand from industry we have been developing the latter approach, focusing on the microscopicphenomena mainly to parametrise them.

In this paper we present the numerical model and its results for industrial-scale biomass gasifier obtained par-allel from our ANSYS Fluent  biomass gasification  UDF-extention (more details in [1])and in-house OpenFOAMporousBiomassGasificationFoam   solver (see also [2]). The model takes into consideration air and syngas flowand reactions within porous layers, water evaporation from biomass, pyrolysis, gasification and heterogeneouscombustion of carbon remnants [3].

Figure 1: Installationof wood chips gasifica-tion and sawdut drying topellets production, Szepi-etowo.

Figure 2: Dry and wetwood chips waiting tobe gasified. The com-position of the woodchips is complex, containsfragmented plywood,chipboard, straw, boards,

bark, sawdust and woodshaving.

The biomass bed formed in the gasification chamber is treated as a porous zone,where Darcy’s law is applicable. This approximation is sufficient, since the Reynoldsnumber is less than 300. Still both Fluent and OpenFoam solve full Navier-Stokes

equation with the additional momentum source term. All parameters of the porouslayer are anisotropic, non-uniform and transient. Generally porosity is increasedin the pyrolysis zone then the bed rapidly crashes and porosity is decreased. Toinitialise the model we made observations in the real-life installation and developedsupplementary zonal and one-dimensional model of gasifier in Matlab (more detailsin   [4]). In the zonal approach each individual zone is artificially separated andmodelled, while in the 1D and CFD approaches there are no artificial division forzones, only the sequence of the process is clarify.

In order to model the four mentioned above processes it is crucial to properlymodel the heat transfer b etween solid (biomass) and gas (air and syngas). Weassumed that the heat exchange in the porous medium is purely convective and wethe used well established heat transfer parameters. Then we introduce the globalheat transfer coefficient based on experimental results.

In order to make our approach versatile, and to deal with various kind of wasteand biomass, we implemented the data structures containing the elemental com-position of solid, its moisture content, volatile particles, fixed carbon and ash con-tent [5]. Additionally volatiles are quantified by fractional composition of lignin,cellulose and hemicellulose for vegetable biomass [6].

We include the homogeneous chemical reactions of water-gas shift, combustionof CO, CH4 ,H2 and heterogeneous reactions such as gasification and char combus-tion. The reaction rates are defined by the chemical kinetics based on the Arrheniusequation (for homogeneous reactions) and the Langmuir-Hinshelwood (for hetero-geneous ones). In the near future we intend to replace the coefficients taken fromliterature [5] by the data obtained from the thermogravimetric measurements thathave already been scheduled.

The results, which we present below are prepared for the industrial updraft,fixed-bed gasifier fuelled by wood chips. It is located in Szepietowo and achieves3.5MW of heat power used for sawdust drying (fig.   1 and 2). The sample resultsobtained from both OpenFOAM and Fluent are presented in figures 3  and  4.  The shape of the feed, marked inthe figures in red, was assumed drawing from the experience of the installation operator. This shape can varydepending mainly on the sizes of the chips. Our model will be experimentally verified for this installation.

1

Page 2: Gas if Ication Modelling

8/13/2019 Gas if Ication Modelling

http://slidepdf.com/reader/full/gas-if-ication-modelling 2/2

(a) OpenFOAM (b) ANSYS Fluent

Figure 3: Results of our simulations of the flow in a porous layer of biomass (the shape of the porous zone shownin by red) in the simplified biomass gasifier geometry. The velocity fields simulated by the in-house OpenFoamsolver are compared with the result from Fluent (with additional UDF). The air inlets are symmetrically locatedat the bottom. Outlets are located in the top, right-hand side of the gasifier.

Figure 4: Comparison of the pressure distributions obtained with OpenFoam and Fluent. Small disagreementis visible in the lower part of the gasifier. The decrease of pressure with height is in good agreement for bothcodes.

References

[1] Korotko, J., Kwiatkowski, K., Bajer, K. Numerical modelling of thermodynamic processes in biomass bed,IX Workshop Modelling multiphase flows in thermochemical systems, Wiezyca, Poland 2011.

[2] Kwiatkowski,K.,  Zuk, P., Wedolowski, K., Bajer, K. Flow, reactions and production of syngas in the porousbiomass layer. 6th OpenFOAM Workshop, PennState University, USA, 2011.

[3] Souza-Santoz de,M.L.,  Solid Fuels Combustion and Gasification , MARCEL DEKKER Inc., 2004.

[4] Gorecki, B., Kwiatkowski, K., Bajer, K. Numerical simulation of biomass gasification - computational codeand modelling, IX Workshop Modelling multiphase flows in thermochemical systems, Wiezyca, Poland 2011.

[5] Basu, P., Biomass gasification and pyrolysis, Practical design , Elsevier, 2010.

[6] Di Blasi, C., Modelling chemical and physical processes of wood and biomass pyrolysis, Progress in Energyand Combustion Science, 34, 47–90, 2008.

2