A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy

doi:https://doi.org/10.1016/j.jcp.2016.02.059

Title: A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy

Abstract

A novel and efficient algorithm is presented in this paper to deal with DNS of turbulent reacting flows under the low-Mach-number assumption, with detailed chemistry and a quasi-spectral accuracy. The temporal integration of the equations relies on an operating-split strategy, where chemical reactions are solved implicitly with a stiff solver and the convection-diffusion operators are solved with a Runge-Kutta-Chebyshev method. The spatial discretisation is performed with high-order compact schemes, and a FFT based constant-coefficient spectral solver is employed to solve a variable-coefficient Poisson equation. The numerical implementation takes advantage of the 2DECOMP&FFT libraries developed by, which are based on a pencil decomposition method of the domain and are proven to be computationally very efficient. An enhanced pressure-correction method is proposed to speed up the achievement of machine precision accuracy. It is demonstrated that a second-order accuracy is reached in time, while the spatial accuracy ranges from fourth-order to sixth-order depending on the set of imposed boundary conditions. In conclusion, the software developed to implement the present algorithm is called HOLOMAC, and its numerical efficiency opens the way to deal with DNS of reacting flows to understand complex turbulent and chemical phenomena in flames.

Authors:

^[1]; Abraham, J. ^[2]

Univ. of Adelaide, SA (Australia). School of Mechanical Engineering
Univ. of Adelaide, SA (Australia). School of Mechanical Engineering; Purdue Univ., West Lafayette, IN (United States). School of Mechanical Engineering

Publication Date:: 2016-02-27

Research Org.:: Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1525148

Grant/Contract Number:: AC02-05CH11231

Resource Type:: Journal Article: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 313; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; DNS; Low-Mach-number; Detailed chemistry; Turbulent reacting flow; High-order methods; Spectral accuracy; Operator splitting

Citation Formats


                    Motheau, E., and Abraham, J. A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy.  United States: N. p., 2016. 
        Web.  doi:10.1016/j.jcp.2016.02.059.

Copy to clipboard


                    Motheau, E., & Abraham, J. A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy.  United States.  https://doi.org/10.1016/j.jcp.2016.02.059

Copy to clipboard


                    Motheau, E., and Abraham, J. Sat .  
        "A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy".  United States.  https://doi.org/10.1016/j.jcp.2016.02.059.  https://www.osti.gov/servlets/purl/1525148.

Copy to clipboard


                    
@article{osti_1525148,

  title        = {A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy},

  author       = {Motheau, E. and Abraham, J.},

  abstractNote = {A novel and efficient algorithm is presented in this paper to deal with DNS of turbulent reacting flows under the low-Mach-number assumption, with detailed chemistry and a quasi-spectral accuracy. The temporal integration of the equations relies on an operating-split strategy, where chemical reactions are solved implicitly with a stiff solver and the convection-diffusion operators are solved with a Runge-Kutta-Chebyshev method. The spatial discretisation is performed with high-order compact schemes, and a FFT based constant-coefficient spectral solver is employed to solve a variable-coefficient Poisson equation. The numerical implementation takes advantage of the 2DECOMP&FFT libraries developed by, which are based on a pencil decomposition method of the domain and are proven to be computationally very efficient. An enhanced pressure-correction method is proposed to speed up the achievement of machine precision accuracy. It is demonstrated that a second-order accuracy is reached in time, while the spatial accuracy ranges from fourth-order to sixth-order depending on the set of imposed boundary conditions. In conclusion, the software developed to implement the present algorithm is called HOLOMAC, and its numerical efficiency opens the way to deal with DNS of reacting flows to understand complex turbulent and chemical phenomena in flames.},

  doi          = {10.1016/j.jcp.2016.02.059},

  url          = {https://www.osti.gov/biblio/1525148},
  journal      = {Journal of Computational Physics},

  issn         = {0021-9991},

  number       = C,

  volume       = 313,

  place        = {United States},

  year         = {2016},

  month        = {2}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2016.02.059

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 20 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referencing / citing this record:

A fourth-order adaptive mesh refinement algorithm for the multicomponent, reacting compressible Navier–Stokes equations
journal, January 2019

Emmett, Matthew; Motheau, Emmanuel; Zhang, Weiqun
Combustion Theory and Modelling, Vol. 23, Issue 4
https://doi.org/10.1080/13647830.2019.1566574

Similar records in OSTI.GOV collections:

PIROCK: A swiss-knife partitioned implicit–explicit orthogonal Runge–Kutta Chebyshev integrator for stiff diffusion–advection–reaction problems with or without noise

Journal Article Abdulle, Assyr ; Vilmart, Gilles - Journal of Computational Physics

A partitioned implicit–explicit orthogonal Runge–Kutta method (PIROCK) is proposed for the time integration of diffusion–advection–reaction problems with possibly severely stiff reaction terms and stiff stochastic terms. The diffusion terms are solved by the explicit second order orthogonal Chebyshev method (ROCK2), while the stiff reaction terms (solved implicitly) and the advection and noise terms (solved explicitly) are integrated in the algorithm as finishing procedures. It is shown that the various coupling (between diffusion, reaction, advection and noise) can be stabilized in the PIROCK method. The method, implemented in a single black-box code that is fully adaptive, provides error estimators for themore »« less
https://doi.org/10.1016/J.JCP.2013.02.009
Recovery Discontinuous Galerkin Jacobian-Free Newton-Krylov Method for All-Speed Flows

Conference Park, HyeongKae ; Nourgaliev, Robert ; Mousseau, Vincent ; ...

A novel numerical algorithm (rDG-JFNK) for all-speed fluid flows with heat conduction and viscosity is introduced. The rDG-JFNK combines the Discontinuous Galerkin spatial discretization with the implicit Runge-Kutta time integration under the Jacobian-free Newton-Krylov framework. We solve fully-compressible Navier-Stokes equations without operator-splitting of hyperbolic, diffusion and reaction terms, which enables fully-coupled high-order temporal discretization. The stability constraint is removed due to the L-stable Explicit, Singly Diagonal Implicit Runge-Kutta (ESDIRK) scheme. The governing equations are solved in the conservative form, which allows one to accurately compute shock dynamics, as well as low-speed flows. For spatial discretization, we develop a “recovery” familymore »« less
Full Text Available
Fast and accurate implementation of Fourier spectral approximations of nonlocal diffusion operators and its applications

Journal Article Du, Qiang ; Yang, Jiang - Journal of Computational Physics

This work is concerned with the Fourier spectral approximation of various integral differential equations associated with some linear nonlocal diffusion and peridynamic operators under periodic boundary conditions. For radially symmetric kernels, the nonlocal operators under consideration are diagonalizable in the Fourier space so that the main computational challenge is on the accurate and fast evaluation of their eigenvalues or Fourier symbols consisting of possibly singular and highly oscillatory integrals. For a large class of fractional power-like kernels, we propose a new approach based on reformulating the Fourier symbols both as coefficients of a series expansion and solutions of some simplemore »« less
https://doi.org/10.1016/J.JCP.2016.11.028
Numerical study of pseudospectral methods in shock wave applications

Journal Article Don, Wai Sun - Journal of Computational Physics; (United States)

In this paper the authors use shock capturing spectral methods to simulate compressible flows in the presence of shock waves. Three applications are considered. The first is an interaction of a one-dimensional shock with an entropy wave. The second case deals with interactions between shock wave and an entropy wave in two space dimensions and the third case is two-dimensional shock-vortex interactions. The first two applications are found in the study of turbulence in high speed flow. The last application is a key element in understanding the acoustic dynamics in aeroacoustic and in design of supersonic jet in aerodynamics. Themore »« less
https://doi.org/10.1006/jcph.1994.1008
An efficient implementation of the Chebyshev Rational Approximation Method (CRAM) for solving the burnup equations

Conference Pusa, M. ; Leppaenen, J.

The Chebyshev Rational Approximation Method (CRAM) has been recently introduced by the authors for solving the burnup equations with excellent results. This method has been shown to be capable of simultaneously solving an entire burnup system with thousands of nuclides both accurately and efficiently. The method was prompted by an analysis of the spectral properties of burnup matrices and it can be characterized as the best rational approximation on the negative real axis. The coefficients of the rational approximation are fixed and have been reported for various approximation orders. In addition to these coefficients, implementing the method only requires amore »« less

Similar Records

Title: A high-order numerical algorithm for DNS of low-Mach-number reactive flows with detailed chemistry and quasi-spectral accuracy

Abstract

Citation Formats

A fourth-order adaptive mesh refinement algorithm for the multicomponent, reacting compressible Navier–Stokes equations journal, January 2019

A fourth-order adaptive mesh refinement algorithm for the multicomponent, reacting compressible Navier–Stokes equations
journal, January 2019