An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation

Su, Yunde; Kim, Seung Hyun

doi:10.1016/j.jcp.2018.06.021

Title: An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation

Abstract

A consistent, conservative and truly high order finite difference scheme for low Mach number turbulent variable density flow simulation is presented. In the present scheme, the scalar transport equation is discretized with the explicit weighted compact nonlinear scheme (WCNS) while the continuity and momentum equations are discretized with a conservative centered high order finite difference scheme on a staggered grid. Bounded hybrid high order upwind interpolation is proposed for the evaluation of density on the cell faces of the staggered grid. The use of explicit WCNS along with the bounded density interpolation resolves the issue of degraded order of accuracy in the existing high order finite difference schemes where the centered finite difference scheme for mass and momentum conservation is combined with a total variation stable scheme for scalar transport. With the same difference operator used for the continuity and scalar transport equations, the present scheme is consistent naturally, and resolves the inconsistency problem reported in Trisjono et al. (J. Comp. Phys., 327, 612-628, 2016) without employing any modification to a scheme. Such improvements are demonstrated for several test cases.

Authors:

Su, Yunde ^[1]; Kim, Seung Hyun ^[1]

The Ohio State Univ., Columbus, OH (United States)

Publication Date:: Fri Jun 15 00:00:00 EDT 2018

Research Org.:: The Ohio State Univ., Columbus, OH (United States)

Sponsoring Org.:: USDOE Office of Energy Efficiency and Renewable Energy (EERE), Vehicle Technologies Office (EE-3V); USDOE Office of Energy Efficiency and Renewable Energy (EERE)

OSTI Identifier:: 1513501

Alternate Identifier(s):: OSTI ID: 1702368

Grant/Contract Number:: EE0007334

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

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

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 42 ENGINEERING; 97 MATHEMATICS AND COMPUTING; 02 PETROLEUM; High order finite difference scheme; Low Mach number variable density flow; Conservative scheme; Scalar transport; Weighted compact nonlinear scheme (WCNS); Turbulent reacting flow

Citation Formats


                    Su, Yunde, and Kim, Seung Hyun. An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation.  United States: N. p., 2018. 
Web.  doi:10.1016/j.jcp.2018.06.021.

Copy to clipboard


                    Su, Yunde, & Kim, Seung Hyun. An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation.  United States.  https://doi.org/10.1016/j.jcp.2018.06.021

Copy to clipboard


                    Su, Yunde, and Kim, Seung Hyun. Fri .  
"An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation".  United States.  https://doi.org/10.1016/j.jcp.2018.06.021.  https://www.osti.gov/servlets/purl/1513501.

Copy to clipboard


                    
@article{osti_1513501,

  title        = {An improved consistent, conservative, non-oscillatory and high order finite difference scheme for variable density low Mach number turbulent flow simulation},

  author       = {Su, Yunde and Kim, Seung Hyun},

  abstractNote = {A consistent, conservative and truly high order finite difference scheme for low Mach number turbulent variable density flow simulation is presented. In the present scheme, the scalar transport equation is discretized with the explicit weighted compact nonlinear scheme (WCNS) while the continuity and momentum equations are discretized with a conservative centered high order finite difference scheme on a staggered grid. Bounded hybrid high order upwind interpolation is proposed for the evaluation of density on the cell faces of the staggered grid. The use of explicit WCNS along with the bounded density interpolation resolves the issue of degraded order of accuracy in the existing high order finite difference schemes where the centered finite difference scheme for mass and momentum conservation is combined with a total variation stable scheme for scalar transport. With the same difference operator used for the continuity and scalar transport equations, the present scheme is consistent naturally, and resolves the inconsistency problem reported in Trisjono et al. (J. Comp. Phys., 327, 612-628, 2016) without employing any modification to a scheme. Such improvements are demonstrated for several test cases.},

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

  journal      = {Journal of Computational Physics},

  number       = C,

  volume       = 372,

  place        = {United States},

  year         = {Fri Jun 15 00:00:00 EDT 2018},

  month        = {Fri Jun 15 00:00:00 EDT 2018}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

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

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 4 works

Citation information provided by
Web of Science

Figures / Tables:

Figure 1: 2D Cartesian staggered grid.

All figures and tables (14 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar Records in DOE PAGES and OSTI.GOV collections:

An explicit staggered-grid method for numerical simulation of large-scale natural gas pipeline networks

Journal Article Gyrya, Vitaliy ; Zlotnik, Anatoly - Applied Mathematical Modelling

Here, we present an explicit second order staggered finite difference (FD) discretization scheme for forward simulation of natural gas transport in pipeline networks. By construction, this discretization approach guarantees that the conservation of mass condition is satisfied exactly. The mathematical model is formulated in terms of density, pressure, and mass flux variables, and as a result permits the use of a general equation of state to define the relation between the gas density and pressure for a given temperature. In a single pipe, the model represents the dynamics of the density by propagation of a non-linear wave according to amore »« less
Cited by 14
https://doi.org/10.1016/j.apm.2018.07.051
Second order multidimensional sign-preserving remapping for ALE methods

Conference Hill, Ryan N ; Szmelter, J

A second-order conservative sign-preserving remapping scheme for Arbitrary Lagrangian-Eulerian (ALE) methods is developed utilising concepts of the Multidimensional Positive Definite Advection Transport Algorithm (MPDATA). The algorithm is inherently multidimensional, and so does not introduce splitting errors. The remapping is implemented in a two-dimensional, finite element ALE solver employing staggered quadrilateral meshes. The MPDATA remapping uses a finite volume discretization developed for volume coordinates. It is applied for the remapping of density and internal energy arranged as cell centered, and velocity as nodal, dependent variables. In the paper, the advection of scalar fields is examined first for test cases with prescribedmore »« less
Full Text Available
Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

Technical Report Wang, Dean ; Downar, Thomas ; Xu, Yunlin ; ...

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transportmore »« less
https://doi.org/10.2172/1511575

Full Text Available
Segregated prediction of 3-D compressible subsonic fluid flows using collocated grids

Journal Article Zhang, G ; Assanis, D N ; Tamamidis, P - Numerical Heat Transfer. Part A, Applications

A segregated approach for the prediction of three-dimensional, compressible, subsonic flows is presented. The method uses a collocated finite volume scheme in body-fitted coordinates. For computational expediency and ease of implementation of high-order schemes, the continuity, momentum, and state equations are decoupled and solved sequentially. In order to eliminate potential ``checkerboard`` pressure distributions due to the decoupling of the pressure-velocity terms in conjunction with a collocated grid, a novel method of nonlinear cell-face interpolation is adopted. The high-order QUICK scheme is implemented for discretizing convective terms, while the central difference scheme is used for discretizing diffusion terms. A four-stage Runge-Kuttamore »« less
https://doi.org/10.1080/10407789608913818
A mass and momentum conserving unsplit semi-Lagrangian framework for simulating multiphase flows

Journal Article Owkes, Mark ; Desjardins, Olivier - Journal of Computational Physics

In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas–liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the three-dimensional, unsplit, second-order semi-Lagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum defined on a staggered grid and discrete conservation of mass and momentum, evenmore »« less
https://doi.org/10.1016/J.JCP.2016.11.046

Similar Records