Conservative High-Order Finite-Difference Schemes for Low-Mach Number Flows

Abstract

Three finite-difference algorithms are proposed to solve a low-Mach number approximation for the Navier-Stokes equations. These algorithms exhibit fourth-order spatial and second-order temporal accuracy. They are dissipation-free, and thus well suited for DNS and LES of turbulent flows. The key ingredient common to each of the methods presented is a Poisson equation with variable coefficient that is solved for the hydrodynamic pressure. This feature ensures that the velocity field is constrained correctly. It is shown that this approach is needed to avoid violation of the conservation of kinetic energy in the inviscid limit which would otherwise arise through the pressure term in the momentum equation. An existing set of finite-difference formulae for incompressible flow is generalized to handle arbitrary large density fluctuations with no violation of conservation through the non-linear convective terms. An algorithm which conserves mass, momentum, and kinetic energy fully is obtained when an approximate equation of state is used instead of the exact one. Results from a model problem are used to show both spatial and temporal convergence rates and several test cases are presented to illustrate the performance of the algorithms.

A Cartesian grid embedded boundary method for solving the Poisson and heat equations with discontinuous coefficients in three dimensions

“…The method is conservative, a distinct advantage in certain classes of problems. For low-Mach flows with heat transfer, for instance, conservative schemes avoid unphysical results arising in marginally resolved or under-resolved situations [8,9].…”

Unpacking physics representations: Towards an appreciation of disciplinary affordance

“…For example, the clause "kinetic energy… is conserved" [41] (p. 300) can be rankshifted to form the nominal group "conservation of kinetic energy." This nominal group can then be used to fill the function of a noun in a sentence such as "…this approach is needed to avoid violation of the conservation of kinetic energy…" [42] (p. 71, italics ours). This example also illustrates how the rankshift or nominalization process can be iterated, whereby a nominal group, like "conservation of kinetic energy" mentioned above, can become embedded in an even larger nominal group such as, in the exemplar sentence just given, "violation of the conservation of kinetic energy.…”

Simplifying numerical solution of constrained PDE systems through involutive completion

“…The equations are derived from the Navier-Stokes system where assumptions are made on the state equation linking the density ρ and temperature T [16]:…”