skip to main content

DOE PAGESDOE PAGES

Title: Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.
Authors:
 [1] ;  [2] ;  [2]
  1. North Carolina State Univ., Raleigh, NC (United States)
  2. Univ. of Minnesota, Minneapolis, MN (United States)
Publication Date:
Grant/Contract Number:
NA0002382; FA9550-12-1-0461; FA9550-12-1-0064
Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 348; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Research Org:
California Inst. of Technology (CalTech), Pasadena, CA (United States)
Sponsoring Org:
USDOE National Nuclear Security Administration (NNSA), Office of Defense Programs (DP) (NA-10); US Air Force Office of Scientific Research (AFOSR)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Numerical methods; Compressible flow; Boundedness; Reacting flow; Scalar excursions; Multispecies
OSTI Identifier:
1434335

Subbareddy, Pramod K., Kartha, Anand, and Candler, Graham V.. Scalar conservation and boundedness in simulations of compressible flow. United States: N. p., Web. doi:10.1016/j.jcp.2017.08.001.
Subbareddy, Pramod K., Kartha, Anand, & Candler, Graham V.. Scalar conservation and boundedness in simulations of compressible flow. United States. doi:10.1016/j.jcp.2017.08.001.
Subbareddy, Pramod K., Kartha, Anand, and Candler, Graham V.. 2017. "Scalar conservation and boundedness in simulations of compressible flow". United States. doi:10.1016/j.jcp.2017.08.001. https://www.osti.gov/servlets/purl/1434335.
@article{osti_1434335,
title = {Scalar conservation and boundedness in simulations of compressible flow},
author = {Subbareddy, Pramod K. and Kartha, Anand and Candler, Graham V.},
abstractNote = {With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.},
doi = {10.1016/j.jcp.2017.08.001},
journal = {Journal of Computational Physics},
number = C,
volume = 348,
place = {United States},
year = {2017},
month = {8}
}