Advanced

An improved high-order scheme for DNS of low Mach number turbulent reacting flows based on stiff chemistry solver

Yu, Rixin LU ; Yu, Jiangfei LU and Bai, Xue-Song LU (2012) In Journal of Computational Physics 231(16). p.5504-5521
Abstract
We present an improved numerical scheme for numerical simulations of low Mach number turbulent reacting flows with detailed chemistry and transport. The method is based on a semi-implicit operator-splitting scheme with a stiff solver for integration of the chemical kinetic rates, developed by Knio et al. [O.M. Knio, H.N. Najm, P.S. Wyckoff, A semi-implicit numerical scheme for reacting flow II. Stiff, operator-split formulation, Journal of Computational Physics 154 (2) (1999) 428-467]. Using the material derivative form of continuity equation, we enhance the scheme to allow for large density ratio in the flow field. The scheme is developed for direct numerical simulation of turbulent reacting flow by employing high-order discretization for... (More)
We present an improved numerical scheme for numerical simulations of low Mach number turbulent reacting flows with detailed chemistry and transport. The method is based on a semi-implicit operator-splitting scheme with a stiff solver for integration of the chemical kinetic rates, developed by Knio et al. [O.M. Knio, H.N. Najm, P.S. Wyckoff, A semi-implicit numerical scheme for reacting flow II. Stiff, operator-split formulation, Journal of Computational Physics 154 (2) (1999) 428-467]. Using the material derivative form of continuity equation, we enhance the scheme to allow for large density ratio in the flow field. The scheme is developed for direct numerical simulation of turbulent reacting flow by employing high-order discretization for the spatial terms. The accuracy of the scheme in space and time is verified by examining the grid/time-step dependency on one-dimensional benchmark cases: a freely propagating premixed flame in an open environment and in an enclosure related to spark-ignition engines. The scheme is then examined in simulations of a two-dimensional laminar flame/vortex-pair interaction. Furthermore, we apply the scheme to direct numerical simulation of a homogeneous charge compression ignition (HCCI) process in an enclosure studied previously in the literature. Satisfactory agreement is found in terms of the overall ignition behavior, local reaction zone structures and statistical quantities. Finally, the scheme is used to study the development of intrinsic flame instabilities in a lean H-2/air premixed flame, where it is shown that the spatial and temporary accuracies of numerical schemes can have great impact on the prediction of the sensitive nonlinear evolution process of flame instability. Crown Copyright (C) 2012 Published by Elsevier Inc. All rights reserved. (Less)
Please use this url to cite or link to this publication:
author
organization
publishing date
type
Contribution to journal
publication status
published
subject
keywords
Low Mach number, Stiff, Detailed chemistry, Operator splitting, High, order scheme, Variable density, DNS, Turbulent reacting flow
in
Journal of Computational Physics
volume
231
issue
16
pages
5504 - 5521
publisher
Elsevier
external identifiers
  • wos:000305340900016
  • scopus:84862262046
ISSN
0021-9991
DOI
10.1016/j.jcp.2012.05.006
language
English
LU publication?
yes
id
bac99ce3-2f3c-43d0-b36b-b2dc5365d572 (old id 2887437)
alternative location
http://www.sciencedirect.com/science/article/pii/S0021999112002306
date added to LUP
2012-07-25 11:03:23
date last changed
2017-07-23 03:06:23
@article{bac99ce3-2f3c-43d0-b36b-b2dc5365d572,
  abstract     = {We present an improved numerical scheme for numerical simulations of low Mach number turbulent reacting flows with detailed chemistry and transport. The method is based on a semi-implicit operator-splitting scheme with a stiff solver for integration of the chemical kinetic rates, developed by Knio et al. [O.M. Knio, H.N. Najm, P.S. Wyckoff, A semi-implicit numerical scheme for reacting flow II. Stiff, operator-split formulation, Journal of Computational Physics 154 (2) (1999) 428-467]. Using the material derivative form of continuity equation, we enhance the scheme to allow for large density ratio in the flow field. The scheme is developed for direct numerical simulation of turbulent reacting flow by employing high-order discretization for the spatial terms. The accuracy of the scheme in space and time is verified by examining the grid/time-step dependency on one-dimensional benchmark cases: a freely propagating premixed flame in an open environment and in an enclosure related to spark-ignition engines. The scheme is then examined in simulations of a two-dimensional laminar flame/vortex-pair interaction. Furthermore, we apply the scheme to direct numerical simulation of a homogeneous charge compression ignition (HCCI) process in an enclosure studied previously in the literature. Satisfactory agreement is found in terms of the overall ignition behavior, local reaction zone structures and statistical quantities. Finally, the scheme is used to study the development of intrinsic flame instabilities in a lean H-2/air premixed flame, where it is shown that the spatial and temporary accuracies of numerical schemes can have great impact on the prediction of the sensitive nonlinear evolution process of flame instability. Crown Copyright (C) 2012 Published by Elsevier Inc. All rights reserved.},
  author       = {Yu, Rixin and Yu, Jiangfei and Bai, Xue-Song},
  issn         = {0021-9991},
  keyword      = {Low Mach number,Stiff,Detailed chemistry,Operator splitting,High,order scheme,Variable density,DNS,Turbulent reacting flow},
  language     = {eng},
  number       = {16},
  pages        = {5504--5521},
  publisher    = {Elsevier},
  series       = {Journal of Computational Physics},
  title        = {An improved high-order scheme for DNS of low Mach number turbulent reacting flows based on stiff chemistry solver},
  url          = {http://dx.doi.org/10.1016/j.jcp.2012.05.006},
  volume       = {231},
  year         = {2012},
}