CFD simulations of steam cracking furnaces using detailed combustion mechanisms

Abstract A three-dimensional mathematical model has been developed for the simulation of flow, temperature and concentration fields in the radiation section of industrial scale steam cracking units. The model takes into account turbulence–chemistry interactions through the Eddy Dissipation Concept (EDC) model and makes use of Detailed Reaction Kinetics (DRK), which allows the detailed investigation of the flame structure. Furthermore, simulation results obtained with the EDC-DRK model are compared with simulation results obtained with a simplified model combining the Eddy Break Up (EBU)/finite rate formulation with Simplified Reaction Kinetics (SRK). When the EBU-SRK model is used, much faster fuel oxidation and products formation is predicted. The location of the peak temperature is shifted towards the burner, resulting in a smaller flame and the confinement of the combustion process into a smaller area. This is most likely because of the inherent deficiency of the simplified model to correctly predict the overall (effective) burning rate when the turbulent mixing rate and the reaction rate are comparable. It is shown that when neither the “fast-chemistry” nor the “slow-chemistry” approximation is satisfied, the overall burning rate is overpredicted. The smaller flame volumes obtained with the EBU-SRK model have important effects on the predicted temperature distribution in the furnace as well as on other significant design parameters like the refractory wall and tube skin temperatures. It is suggested that more sophisticated turbulence–chemistry interaction models like the EDC model and more Detailed Reaction Kinetics should be used for combustion modeling in steam cracking furnaces under normal firing conditions.

[1]  Gilbert F. Froment,et al.  Kinetic modeling of the thermal cracking of hydrocarbons. 2. Calculation of activation energies , 1988 .

[2]  S. Orszag,et al.  Development of turbulence models for shear flows by a double expansion technique , 1992 .

[3]  J. Steelant,et al.  Coupled solution of the steady compressible Navier-Stokes equations and the k -e turbulence equations with a multigrid method , 1997 .

[4]  A. Brink,et al.  Possibilities and limitations of the Eddy Break-Up Model , 2000 .

[5]  Ahmad H. Haidari,et al.  Design reactors via CFD , 2001 .

[6]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[7]  Bjørn F. Magnussen,et al.  A Numerical Study of a Bluff-Body Stabilized Diffusion Flame. Part 2. Influence of Combustion Modeling And Finite-Rate Chemistry , 1996 .

[8]  Guy Marin,et al.  Three-Dimensional Asymmetric Flow and Temperature Fields in Cracking Furnaces† , 2001 .

[9]  Rodney O. Fox,et al.  Computational Methods for Turbulent Reacting Flows in the Chemical Process Industry , 1996 .

[10]  C. Westbrook,et al.  Simplified Reaction Mechanisms for the Oxidation of Hydrocarbon Fuels in Flames , 1981 .

[11]  G. Froment,et al.  Coke Formation in the Thermal Cracking of Hydrocarbons. 4. Modeling of Coke Formation in Naphtha Cracking , 1994 .

[12]  Debasis Chakraborty,et al.  Evaluation of combustion models for high speed H2/air confined mixing layer using DNS data , 2000 .

[13]  Tobias S. Christensen,et al.  Post-processing of detailed chemical kinetic mechanisms onto CFD simulations , 2004, Comput. Chem. Eng..

[14]  Grégoire Mallard,et al.  The NIST Chemical Kinetics Database , 1998 .

[15]  S. Pope.  PDF methods for turbulent reactive flows , 1985 .

[16]  W. Jones,et al.  Calculation methods for reacting turbulent flows: A review , 1982 .

[17]  Gilbert F. Froment,et al.  Simulation of the run length of an ethane cracking furnace , 1990 .

[18]  F. Williams,et al.  Turbulent Reacting Flows , 1981 .

[19]  I. S. Ertesvåg,et al.  The Eddy Dissipation Turbulence Energy Cascade Model , 2000 .

[20]  D. K. Edwards,et al.  Radiation Heat Transfer in Nonisothermal Nongray Gases , 1967 .

[21]  Guy Marin,et al.  Effect of radial temperature profiles on yields in steam cracking , 2004 .

[22]  R. W. Dibble,et al.  Application of reduced chemical mechanisms for prediction of turbulent nonpremixed methane jet flames , 1991 .

[23]  Geraldine Heynderickx,et al.  Calculation of three-dimensional flow and pressure fields in cracking furnaces , 2003 .

[24]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[25]  Uwe Schnell,et al.  Detailed modeling of hybrid reburn/SNCR processes for NOX reduction in coal-fired furnaces , 2003 .

[26]  Erik Dick,et al.  Three‐dimensional flow patterns in cracking furnaces with long‐flame burners , 2001 .

[27]  Maria A. Founti,et al.  Radiative heat transfer in natural gas-fired furnaces , 2000 .

[28]  Claudio Bruno,et al.  The coupling of turbulence and chemistry in a premixed bluff-body flame as studied by LES , 2004 .

[29]  Erik Dick,et al.  Multigrid formulation of polynomial flux-difference splitting for steady Euler equations , 1990 .

[30]  Jack Dongarra,et al.  LAPACK Users' Guide, 3rd ed. , 1999 .

[31]  Gilbert F. Froment,et al.  Kinetic modeling of the thermal cracking of hydrocarbons. 1. Calculation of frequency factors , 1988 .

[32]  B. Hjertager,et al.  On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion , 1977 .