Semi-implicit iterative methods for low Mach number turbulent reacting flows: Operator splitting versus approximate factorization

Jonathan F. MacArt, Michael Edward Mueller

Research output: Contribution to journalArticlepeer-review

55 Scopus citations

Abstract

Two formally second-order accurate, semi-implicit, iterative methods for the solution of scalar transport–reaction equations are developed for Direct Numerical Simulation (DNS) of low Mach number turbulent reacting flows. The first is a monolithic scheme based on a linearly implicit midpoint method utilizing an approximately factorized exact Jacobian of the transport and reaction operators. The second is an operator splitting scheme based on the Strang splitting approach. The accuracy properties of these schemes, as well as their stability, cost, and the effect of chemical mechanism size on relative performance, are assessed in two one-dimensional test configurations comprising an unsteady premixed flame and an unsteady nonpremixed ignition, which have substantially different Damköhler numbers and relative stiffness of transport to chemistry. All schemes demonstrate their formal order of accuracy in the fully-coupled convergence tests. Compared to a (non-)factorized scheme with a diagonal approximation to the chemical Jacobian, the monolithic, factorized scheme using the exact chemical Jacobian is shown to be both more stable and more economical. This is due to an improved convergence rate of the iterative procedure, and the difference between the two schemes in convergence rate grows as the time step increases. The stability properties of the Strang splitting scheme are demonstrated to outpace those of Lie splitting and monolithic schemes in simulations at high Damköhler number; however, in this regime, the monolithic scheme using the approximately factorized exact Jacobian is found to be the most economical at practical CFL numbers. The performance of the schemes is further evaluated in a simulation of a three-dimensional, spatially evolving, turbulent nonpremixed planar jet flame.

Original languageEnglish (US)
Pages (from-to)569-595
Number of pages27
JournalJournal of Computational Physics
Volume326
DOIs
StatePublished - Dec 1 2016

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Keywords

  • Chemical Jacobian
  • Direct numerical simulation
  • Operator splitting
  • Strang splitting
  • Turbulent reacting flows

Fingerprint

Dive into the research topics of 'Semi-implicit iterative methods for low Mach number turbulent reacting flows: Operator splitting versus approximate factorization'. Together they form a unique fingerprint.

Cite this