Convergence of the MAC scheme for the compressible stationary Navier-Stokes equations

We prove in this paper the convergence of the Marker and Cell (MAC) scheme for the discretization of the steady state compressible and isentropic Navier-Stokes equations on two or three-dimensional Cartesian grids. Existence of a solution to the scheme is proven, followed by estimates on approximate solutions, which yield the convergence of the approximate solutions, up to a subsequence, and in an appropriate sense. We then prove that the limit of the approximate solutions satisfies the mass and momentum balance equations, as well as the equation of state, which is the main difficulty of this study

Discretization of the viscous dissipation term with the MAC scheme

International audienceWe propose a discretization for the MAC scheme of the viscous dissipation term τ (u) : ∇u (where τ (u) stands for the shear stress tensor associated to the velocity field u), which is suitable for the approximation of this term in a conservation equation for a scalar variable. This discretization enjoys the property that the integral over the computational domain Ω of the (discrete) dissipation term is equal to what is obtained when taking the inner product of the (discrete) momentum balance equation by u and integrating over Ω. As a consequence, it may be used as an ingredient to obtain an unconditionally stable scheme for the compressible Navier-Stokes equations. It is also shown, in some model cases, to ensure the strong convergence in L1 of the dissipation term

Error estimates for a numerical approximation to the compressible barotropic Navier-Stokes equations

We present here a general method based on the investigation of the relative energy of the system, that provides an unconditional error estimate for the approximate solution of the barotropic Navier Stokes equations obtained by time and space discretization. We use this methodology to derive an error estimate for a specific DG/finite element scheme for which the convergence was proved in [27]. This is an extended version of the paper submitted to IMAJNA

A bounded upwinding scheme for computing convection-dominated transport problems

A practical high resolution upwind differencing scheme for the numerical solution of convection-dominated transport problems is presented. The scheme is based on TVD and CBC stability criteria and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving the 1D/2D scalar advection equations, 1D inviscid Burgers’ equation, 1D scalar convection–diffusion equation, 1D/2D compressible Euler’s equations, and 2D incompressible Navier–Stokes equations. The numerical results displayed good agreement with other existing numerical and experimental data

A Hybrid Adaptive Low-Mach-Number/Compressible Method: Euler Equations

Flows in which the primary features of interest do not rely on high-frequency acoustic effects, but in which long-wavelength acoustics play a nontrivial role, present a computational challenge. Integrating the entire domain with low-Mach-number methods would remove all acoustic wave propagation, while integrating the entire domain with the fully compressible equations can in some cases be prohibitively expensive due to the CFL time step constraint. For example, simulation of thermoacoustic instabilities might require fine resolution of the fluid/chemistry interaction but not require fine resolution of acoustic effects, yet one does not want to neglect the long-wavelength wave propagation and its interaction with the larger domain. The present paper introduces a new multi-level hybrid algorithm to address these types of phenomena. In this new approach, the fully compressible Euler equations are solved on the entire domain, potentially with local refinement, while their low-Mach-number counterparts are solved on subregions of the domain with higher spatial resolution. The finest of the compressible levels communicates inhomogeneous divergence constraints to the coarsest of the low-Mach-number levels, allowing the low-Mach-number levels to retain the long-wavelength acoustics. The performance of the hybrid method is shown for a series of test cases, including results from a simulation of the aeroacoustic propagation generated from a Kelvin-Helmholtz instability in low-Mach-number mixing layers. It is demonstrated that compared to a purely compressible approach, the hybrid method allows time-steps two orders of magnitude larger at the finest level, leading to an overall reduction of the computational time by a factor of 8

Cumulative reports and publications through December 31, 1990

This document contains a complete list of ICASE reports. Since ICASE reports are intended to be preprints of articles that will appear in journals or conference proceedings, the published reference is included when it is available