Efficient time stepping for reactive turbulent simulations with stiff chemistry

A combination of a steady-state preserving operator splitting method and a semi-implicit integration scheme is proposed for efficient time stepping in simulations of unsteady reacting flows, such as turbulent flames, using detailed chemical kinetic mechanisms. The operator splitting is based on the Simpler balanced splitting method, which is constructed with improved stability properties and reduced computational cost. The method is shown to be capable of stable and accurate prediction of ignition and extinction for reaction-diffusion systems near critical conditions. The ROK4E scheme is designed for semi-implicit integration of spatially independent chemically reacting systems. Being a Rosenbrock-Krylov method, ROK4E utilizes the low-rank approximation of the Jacobian to reduce the cost for integrating the system of ODEs that have relative few stiff components. The efficiency of the scheme is further improved via the careful choice of coefficients to require three right-hand-side evaluations over four stages. Combing these two methods, efficient calculation is achieved for large-scale parallel simulations of turbulent flames.

[1]  C. Law,et al.  Toward accommodating realistic fuel chemistry in large-scale computations , 2009 .

[2]  M. Ihme,et al.  A Flamelet Model with Heat-Loss Effects for Predicting Wall-Heat Transfer in Rocket Engines , 2017 .

[3]  G. Strang On the Construction and Comparison of Difference Schemes , 1968 .

[4]  C. W. Gear,et al.  The automatic integration of ordinary differential equations , 1971, Commun. ACM.

[5]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[6]  B. A. Schmitt,et al.  ROWMAP—a ROW-code with Krylov techniques for large stiff ODEs , 1997 .

[7]  Hao Wu,et al.  A general probabilistic approach for the quantitative assessment of LES combustion models , 2017, 1702.05539.

[8]  M. Ihme,et al.  Application of Pareto-efficient combustion modeling framework to large eddy simulations of turbulent reacting flows , 2017, 1712.00966.

[9]  Comparison of flamelet/progress-variable and finite-rate chemistry LES models in a preconditioning scheme , 2017 .

[10]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[11]  Wenting Sun,et al.  Global Pathway Analysis of the Extinction and Re-ignition of a Turbulent Non-Premixed Flame , 2017 .

[12]  Matthias Ihme,et al.  A Pareto-efficient combustion framework with submodel assignment for predicting complex flame configurations , 2015 .

[13]  Ghenadie Bulat,et al.  Large Eddy Simulation of an industrial gas turbine combustor using reduced chemistry with accurate pollutant prediction , 2017 .

[14]  Matthias Ihme,et al.  An entropy-stable hybrid scheme for simulations of transcritical real-fluid flows , 2017, J. Comput. Phys..

[15]  L. Trefethen,et al.  Numerical linear algebra , 1997 .

[16]  D. Goodwin,et al.  Cantera: An Object-oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes. Version 2.2.0 , 2015 .

[17]  Matthias Ihme,et al.  A Fidelity Adaptive Modeling Framework for Combustion Systems Based on Model Trust-Region , 2015 .

[18]  Jonathan F. MacArt,et al.  Semi-implicit iterative methods for low Mach number turbulent reacting flows: Operator splitting versus approximate factorization , 2015, J. Comput. Phys..

[19]  Matthias Ihme,et al.  Compliance of combustion models for turbulent reacting flow simulations , 2016 .

[20]  Adrian Sandu,et al.  Rosenbrock-Krylov Methods for Large Systems of Differential Equations , 2013, SIAM J. Sci. Comput..

[21]  P. Rentrop,et al.  Generalized Runge-Kutta methods of order four with stepsize control for stiff ordinary differential equations , 1979 .

[22]  C. Law,et al.  The effect of flame structure on soot formation and transport in turbulent nonpremixed flames using direct numerical simulation , 2007 .

[23]  Yuan Xuan,et al.  A computationally-efficient, semi-implicit, iterative method for the time-integration of reacting flows with stiff chemistry , 2015, J. Comput. Phys..

[24]  Tianfeng Lu,et al.  Systematic approach to obtain analytic solutions of quasi steady state species in reduced mechanisms. , 2006, The journal of physical chemistry. A.

[25]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[26]  P. Moin,et al.  Unstructured Large Eddy Simulation for Prediction of Noise Issued from Turbulent Jets in Various Configurations , 2011 .

[27]  Gurij Ivanovich Marchuk,et al.  Some application of splitting-up methods to the solution of mathematical physics problems , 1968 .

[28]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[29]  S. Menon,et al.  Sensitivity of predictions to chemical kinetics models in a temporally evolving turbulent non-premixed flame , 2017 .

[30]  Vigor Yang,et al.  A Preconditioned Flux-Differencing Scheme for Chemically Reacting Flows at all Mach Numbers , 1997 .

[31]  Vigor Yang,et al.  Parallel on-the-fly adaptive kinetics in direct numerical simulation of turbulent premixed flame , 2017 .

[32]  Jacqueline H. Chen,et al.  Numerical and experimental investigation of turbulent DME jet flames , 2015 .

[33]  L. Petzold A description of dassl: a differential/algebraic system solver , 1982 .

[34]  R. Alexander Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[35]  A. Hindmarsh,et al.  CVODE, a stiff/nonstiff ODE solver in C , 1996 .

[36]  Kjell Gustafsson,et al.  Control theoretic techniques for stepsize selection in explicit Runge-Kutta methods , 1991, TOMS.

[37]  Zhuyin Ren,et al.  A dynamic adaptive method for hybrid integration of stiff chemistry , 2015 .

[38]  Zhuyin Ren,et al.  Analysis of operator splitting errors for near-limit flame simulations , 2017, J. Comput. Phys..

[39]  William H. Green,et al.  Balanced Splitting and Rebalanced Splitting , 2013, SIAM J. Numer. Anal..

[40]  Hao Wu,et al.  MVP-Workshop Contribution: Modeling of Volvo bluff body flame experiment , 2017 .

[41]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[42]  B. A. Schmitt,et al.  Matrix-free W-methods using a multiple Arnoldi iteration , 1995 .