American Institute of Mathematical Sciences

Advanced
March  2021, 14(3): 769-783. doi: 10.3934/dcdss.2020232

Numerical simulation of fluidization for application in oxyfuel combustion

Michal Beneš , Pavel Eichler , Jakub Klinkovský , Miroslav Kolář , Jakub Solovský , Pavel Strachota , and  Alexandr Žák

Department of Mathematics, FNSPE CTU in Prague, Trojanova 13, 120 00 Praha 2, Czech Republic

* Corresponding author: Pavel Strachota

Received  January 2019 Revised  October 2019 Published  December 2019

This paper is concerned with the simulation of multiphase flow hydrodynamics in an experimental oxyfuel fluidized bed combustor designed for biomass fuels. The aim is to perform cross-validation between several models and solvers that differ in the description of some phenomena in question. We focus on the influence of turbulence modeling, inter-phase drag force models, the presence of biomass in the mixture. Also the possibility to simplify the full 3D description to a quasi-1D model is tested. However, the results indicate that such simplification is not suitable for chaotic phenomena in considered scenarios. The models were developed using ANSYS Fluent and OpenFOAM CFD software packages as well as our in-house CFD code CFBSim. The quantities relevant for comparison (the densities of the dispersed solid phases and the phase velocities) are presented in the form of cross-section averaged vertical profiles.

Keywords: Granular flow, fluidized bed, oxyfuel combustion, inter-phase drag force, ANSYS Fluent, OpenFOAM.
Mathematics Subject Classification: Primary: 76T25, 76N15; Secondary: 70-08.
Citation: Michal Beneš, Pavel Eichler, Jakub Klinkovský, Miroslav Kolář, Jakub Solovský, Pavel Strachota, Alexandr Žák. Numerical simulation of fluidization for application in oxyfuel combustion. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 769-783. doi: 10.3934/dcdss.2020232
References:
[1]

ANSYS, Inc., ANSYS Fluent Theory Guide, Release 15.0, Canonsburg, Pennsylvania, 2013. Google Scholar

[2] P. Basu, Combustion and Gasification in Fluidized Beds, CRC Press, 2006.  doi: 10.1201/9781420005158.  Google Scholar
[3]

P. Bauer, M. Beneš, R. Fučík, H. D. Hoang, V. Klement, R. Máca, J. Mach, T. Oberhuber, P. Strachota, V. Žabka and V. Havlena, Numerical simulation of flow in fluidized beds, Discrete. Cont. Dyn. S. S, 8 (2015), 833-846, URL http://www.aimsciences.org/journals/displayArticlesnew.jsp?paperID=11376. Google Scholar

[4]

M. Beneš, P. Strachota, R. Máca, V. Havlena and J. Mach, A quasi-1D model of biomass co-firing in a circulating fluidized bed boiler, in Finite Volumes for Complex Applications VII. Elliptic, Parabolic, and Hyperbolic Problems, Springer Proc. Math. Stat., Springer, Cham, 78 (2014), 791-799. doi: 10.1007/978-3-319-05591-6_79.  Google Scholar

[5] C. T. Bowman and D. J. Seery, Emissions from Continuous Combustion Systems, Plenum Press, 1972.   Google Scholar
[6]

J. Christiansen, Handbook Series Numerical Integration: Numerical solution of ordinary simultaneous differential equations of the 1st order using a method for automatic step change, Numer. Math., 14 (1970), 317-324.  doi: 10.1007/BF02165587.  Google Scholar

[7] D. Gidaspow, Multiphase Flow and Fluidization: Continuum and Kinetic Theory Description, Academic Press, Inc., Boston, MA, 1994.   Google Scholar
[8]

B. E. Launder and D. B. Spalding, The numerical computation of turbulent flows, Numerical Prediction of Flow, Heat Transfer, Turbulence and Combustion, (1983), 96-116. doi: 10.1016/B978-0-08-030937-8.50016-7.  Google Scholar

[9]

Y. F. Liu and O. Hinrichsen, CFD modeling of bubbling fluidized beds using OpenFOAM: Model validation and comparison of TVD differencing schemes, Computers and Chemical Engineering, 69 (2014), 75-88.  doi: 10.1016/j.compchemeng.2014.07.002.  Google Scholar

[10]

M. Mercedes Maroto-Valer (ed.), Developments and Innovation in Carbon Dioxide (CO2) Capture and Storage Technology, Woodhead Publishing, 2010. Google Scholar

[11]

1 F. Sher, M. A. Pans, C. Sun, C. Snape and H. Liu, Oxy-fuel combustion study of biomass fuels in a 20 kw fluidized bed combustor, Fuel, 215 (2018), 778-786.   Google Scholar

[12] L. D. Smoot and P. J. Smith, Coal Combustion and Gasification, Plenum Press, 1985.  doi: 10.1007/978-1-4757-9721-3.  Google Scholar
[13]

M. Syamlal and T. J. O'Brien, Computer simulation of bubbles in a fluidized bed, AIChE Symp. Ser., 85 (1989), 22-31.   Google Scholar

[14]

L. G. Zheng, Oxy-Fuel Combustion for Power Generation and Carbon Dioxide (CO2) Capture, Woodhead Publishing, 2011. Google Scholar

show all references

References:
[1]

ANSYS, Inc., ANSYS Fluent Theory Guide, Release 15.0, Canonsburg, Pennsylvania, 2013. Google Scholar

[2] P. Basu, Combustion and Gasification in Fluidized Beds, CRC Press, 2006.  doi: 10.1201/9781420005158.  Google Scholar
[3]

P. Bauer, M. Beneš, R. Fučík, H. D. Hoang, V. Klement, R. Máca, J. Mach, T. Oberhuber, P. Strachota, V. Žabka and V. Havlena, Numerical simulation of flow in fluidized beds, Discrete. Cont. Dyn. S. S, 8 (2015), 833-846, URL http://www.aimsciences.org/journals/displayArticlesnew.jsp?paperID=11376. Google Scholar

[4]

M. Beneš, P. Strachota, R. Máca, V. Havlena and J. Mach, A quasi-1D model of biomass co-firing in a circulating fluidized bed boiler, in Finite Volumes for Complex Applications VII. Elliptic, Parabolic, and Hyperbolic Problems, Springer Proc. Math. Stat., Springer, Cham, 78 (2014), 791-799. doi: 10.1007/978-3-319-05591-6_79.  Google Scholar

[5] C. T. Bowman and D. J. Seery, Emissions from Continuous Combustion Systems, Plenum Press, 1972.   Google Scholar
[6]

J. Christiansen, Handbook Series Numerical Integration: Numerical solution of ordinary simultaneous differential equations of the 1st order using a method for automatic step change, Numer. Math., 14 (1970), 317-324.  doi: 10.1007/BF02165587.  Google Scholar

[7] D. Gidaspow, Multiphase Flow and Fluidization: Continuum and Kinetic Theory Description, Academic Press, Inc., Boston, MA, 1994.   Google Scholar
[8]

B. E. Launder and D. B. Spalding, The numerical computation of turbulent flows, Numerical Prediction of Flow, Heat Transfer, Turbulence and Combustion, (1983), 96-116. doi: 10.1016/B978-0-08-030937-8.50016-7.  Google Scholar

[9]

Y. F. Liu and O. Hinrichsen, CFD modeling of bubbling fluidized beds using OpenFOAM: Model validation and comparison of TVD differencing schemes, Computers and Chemical Engineering, 69 (2014), 75-88.  doi: 10.1016/j.compchemeng.2014.07.002.  Google Scholar

[10]

M. Mercedes Maroto-Valer (ed.), Developments and Innovation in Carbon Dioxide (CO2) Capture and Storage Technology, Woodhead Publishing, 2010. Google Scholar

[11]

1 F. Sher, M. A. Pans, C. Sun, C. Snape and H. Liu, Oxy-fuel combustion study of biomass fuels in a 20 kw fluidized bed combustor, Fuel, 215 (2018), 778-786.   Google Scholar

[12] L. D. Smoot and P. J. Smith, Coal Combustion and Gasification, Plenum Press, 1985.  doi: 10.1007/978-1-4757-9721-3.  Google Scholar
[13]

M. Syamlal and T. J. O'Brien, Computer simulation of bubbles in a fluidized bed, AIChE Symp. Ser., 85 (1989), 22-31.   Google Scholar

[14]

L. G. Zheng, Oxy-Fuel Combustion for Power Generation and Carbon Dioxide (CO2) Capture, Woodhead Publishing, 2011. Google Scholar

Figure 1.  Dimensions of the combustion chamber and meshes considered. (A) Dimensions of the cross-section of combustion chamber. (B) Coarse mesh generated by ANSYS Meshing. (C) Cross-section of the mesh in ANSYS. (D) Cross-section of the mesh in OpenFOAM generated by BlockMesh
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 3.  Profiles of the gas velocity averaged across the cross-section of the chamber obtained from the simulations for the circulating (upper curves) and bubbling (lower curves) fluidization regimes
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 2.  ANSYS Fluent simulation of the sand particles dynamics in the lower part of the combustion chamber in the bubbling regime shown by means of the sand volume fraction. Time sequence of snapshots of the longitudinal slice through the center of the 3D domain is displayed. The sequence shows a two-dimensional image of penetration dynamics of the gas phase (the lighter regions) through sand particles (their distribution is indicated by darker regions)
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 4.  Linear density of sand particles along the vertical axis of the combustion chamber. Circulating regime
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 5.  Velocity of sand particles averaged across the cross-section of the chamber. Circulating regime
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 6.  Comparison of biomass linear densities between 1D CFBSim code and ANSYS Fluent
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 7.  Comparison between the profiles of the biomass velocity averaged over the cross-section of the chamber obtained by 1D CFBSim code and ANSYS Fluent
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 8.  Velocity of sand particles averaged across the cross-section of the chamber in circulating regime for Gidaspow [7] and Syamlal [13] drag models
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 9.  Linear density of sand particles in circulating regime in dependence on drag models according to Gidaspow [7] and Syamlal [13]
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 10.  Comparison of linear densities of sand particles with RAS turbulence model and laminar flow available in the OpenFOAM package (solver twoPhaseEulerFoam)
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 11.  Comparison of the profiles of the velocity of sand particles averaged over the cross-section of the chamber obtained using the RAS turbulence model and laminar flow model available in the OpenFOAM package (solver twoPhaseEulerFoam). The profiles are also compared with the velocity profile obtained by the 1D CFBSim model
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 12.  Comparison of linear densities in the two-phase laminar model (OpenFOAM solver twoPhaseEulerFoam) and three-phase laminar model (OpenFOAM solver reactingMultiphaseEulerFoam)
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 13.  Comparison of particle velocities averaged over the cross-section of the chamber in the two-phase laminar model (OpenFOAM solver twoPhaseEulerFoam) and three-phase laminar model (OpenFOAM solver reactingMultiphaseEulerFoam)
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 14.  Comparison of the gas velocity averaged over the cross-section of the chamber under the oxyfuel setting and the air setting. Computed in ANSYS Fluent
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 15.  Comparison of the velocity of sand particles averaged over the cross-section of the chamber under the oxyfuel setting and the air setting. Computed in ANSYS Fluent
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 16.  Comparison of the linear density of sand particles under the oxyfuel setting and the air setting. Computed in ANSYS Fluent
Figure Options

Download full-size image

Download as PowerPoint slide

Table 1.  Phase properties corresponding to given temperature and gas composition (air, oxyfuel). Diameter $ d_b $ is calculated so that the volume of a particle corresponds to the average volume of biomass pellets used in [11]. Note that the kinematic viscosity values of all the phases have the same order of magnitude ($ \propto 10^{-4}\:\mathrm{m}^2\;\mathrm{s}^{-1} $)
Temperature $ T = 1000 \: \mathrm{K} $
Air - density $ \rho_g = 0.3564 \: \mathrm{kg}\;\mathrm{m}^{-3} $
Air - dynamic viscosity $ \mu_g = 4.1923 \cdot 10^{-5} \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Oxyfuel gas - density $ \rho_g = 0.495 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Oxyfuel gas - dynamic viscosity $ \mu_g = 5.8227 \cdot 10^{-5} \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Sand - particle diameter $ d_s = 0.78 \: \mathrm{mm} $
Sand - density $ \rho_s = 2655 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Sand - dynamic viscosity [7] $ \mu_s = 0.5 \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Sand - packing limit $ 0.59 $
Biomass (Miscanthus) - particle diameter $ d_b = 10.36 \: \mathrm{mm} $
Biomass (Miscanthus) - density $ \rho_b = 603 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Biomass (Miscanthus) - dynamic viscosity [7] $ \mu_b = 0.5 \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Biomass (Miscanthus) - packing limit $ 0.5 $
Table Options

Download as excel

Temperature $ T = 1000 \: \mathrm{K} $
Air - density $ \rho_g = 0.3564 \: \mathrm{kg}\;\mathrm{m}^{-3} $
Air - dynamic viscosity $ \mu_g = 4.1923 \cdot 10^{-5} \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Oxyfuel gas - density $ \rho_g = 0.495 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Oxyfuel gas - dynamic viscosity $ \mu_g = 5.8227 \cdot 10^{-5} \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Sand - particle diameter $ d_s = 0.78 \: \mathrm{mm} $
Sand - density $ \rho_s = 2655 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Sand - dynamic viscosity [7] $ \mu_s = 0.5 \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Sand - packing limit $ 0.59 $
Biomass (Miscanthus) - particle diameter $ d_b = 10.36 \: \mathrm{mm} $
Biomass (Miscanthus) - density $ \rho_b = 603 \: \mathrm{kg} \; \mathrm{m}^{-3} $
Biomass (Miscanthus) - dynamic viscosity [7] $ \mu_b = 0.5 \: \mathrm{kg} \; \mathrm{m}^{-1} \; \mathrm{s}^{-1} $
Biomass (Miscanthus) - packing limit $ 0.5 $
[1]

Julian Tugaut. Captivity of the solution to the granular media equation. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2021002
[2]

Yancong Xu, Lijun Wei, Xiaoyu Jiang, Zirui Zhu. Complex dynamics of a SIRS epidemic model with the influence of hospital bed number. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021016
[3]

Xiuli Xu, Xueke Pu. Optimal convergence rates of the magnetohydrodynamic model for quantum plasmas with potential force. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 987-1010. doi: 10.3934/dcdsb.2020150
[4]

Taige Wang, Bing-Yu Zhang. Forced oscillation of viscous Burgers' equation with a time-periodic force. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1205-1221. doi: 10.3934/dcdsb.2020160
[5]

Bernard Bonnard, Jérémy Rouot. Geometric optimal techniques to control the muscular force response to functional electrical stimulation using a non-isometric force-fatigue model. Journal of Geometric Mechanics, 2020  doi: 10.3934/jgm.2020032
[6]

Kuo-Chih Hung, Shin-Hwa Wang. Classification and evolution of bifurcation curves for a porous-medium combustion problem with large activation energy. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020281
[7]

Helmut Abels, Andreas Marquardt. On a linearized Mullins-Sekerka/Stokes system for two-phase flows. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020467
[8]

Helmut Abels, Johannes Kampmann. Existence of weak solutions for a sharp interface model for phase separation on biological membranes. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 331-351. doi: 10.3934/dcdss.2020325
[9]

Tian Ma, Shouhong Wang. Topological phase transition III: Solar surface eruptions and sunspots. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 501-514. doi: 10.3934/dcdsb.2020350
[10]

Lin Shi, Dingshi Li, Kening Lu. Limiting behavior of unstable manifolds for spdes in varying phase spaces. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021020
[11]

Tomáš Smejkal, Jiří Mikyška, Jaromír Kukal. Comparison of modern heuristics on solving the phase stability testing problem. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1161-1180. doi: 10.3934/dcdss.2020227
[12]

Shuxing Chen, Jianzhong Min, Yongqian Zhang. Weak shock solution in supersonic flow past a wedge. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 115-132. doi: 10.3934/dcds.2009.23.115
[13]

Liupeng Wang, Yunqing Huang. Error estimates for second-order SAV finite element method to phase field crystal model. Electronic Research Archive, 2021, 29 (1) : 1735-1752. doi: 10.3934/era.2020089
[14]

Caterina Balzotti, Simone Göttlich. A two-dimensional multi-class traffic flow model. Networks & Heterogeneous Media, 2020  doi: 10.3934/nhm.2020034
[15]

Shuang Liu, Yuan Lou. A functional approach towards eigenvalue problems associated with incompressible flow. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3715-3736. doi: 10.3934/dcds.2020028
[16]

Pablo D. Carrasco, Túlio Vales. A symmetric Random Walk defined by the time-one map of a geodesic flow. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020390
[17]

Joan Carles Tatjer, Arturo Vieiro. Dynamics of the QR-flow for upper Hessenberg real matrices. Discrete & Continuous Dynamical Systems - B, 2021, 26 (3) : 1359-1403. doi: 10.3934/dcdsb.2020166
[18]

Petr Pauš, Shigetoshi Yazaki. Segmentation of color images using mean curvature flow and parametric curves. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1123-1132. doi: 10.3934/dcdss.2020389
[19]

Gui-Qiang Chen, Beixiang Fang. Stability of transonic shock-fronts in three-dimensional conical steady potential flow past a perturbed cone. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 85-114. doi: 10.3934/dcds.2009.23.85
[20]

Peter Frolkovič, Viera Kleinová. A new numerical method for level set motion in normal direction used in optical flow estimation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 851-863. doi: 10.3934/dcdss.2020347

2019 Impact Factor: 1.233

Tools

Metrics

  • PDF downloads (341)
  • HTML views (447)
  • Cited by (0)

Other articles
by authors

[Back to Top]

Copyright © 2020 American Institute of Mathematical Sciences