Abstract:The computation of compressible flows becomes more challenging when the Mach number has different orders of magnitude. When the Mach number is of order one, modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions. However, if the Mach number is small, the acoustic waves lead to stiffness in time and excessively large numerical viscosity, thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation. In this paper, we develop an all-speed asymptotic preserving (AP) numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uniformly stable and accurate for all Mach numbers. Our idea is to split the system into two parts: one involves a slow, nonlinear and conservative hyperbolic system adequate for the use of modern shock capturing methods and the other a linear hyperbolic system which contains the stiff acoustic dynamics, to be solved implicitly. This implicit part is reformulated into a standard pressure Poisson projection system and thus possesses sufficient structure for efficient fast Fourier transform solution techniques. In the zero Mach number limit, the scheme automatically becomes a projection method-like incompressible solver. We present numerical results in one and two dimensions in both compressible and incompressible regimes.

Analysis of the Picard-Newton iteration for the incompressible Boussinesq equations

Stability Analysis of Picard Iteration for Coupled Neutronics/Thermal-Hydraulics Simulations

Improving Newton's method performance by parametrization: the case of Richards equation

Assessing the benefits of approximately exact step sizes for Picard and Newton solver in simulating ice flow (FEniCS-full-Stokes v.1.3.2)

A study on iterative methods for solving Richards’ equation

Multilevel Picard algorithm for general semilinear parabolic PDEs with gradient-dependent nonlinearities

Accelerating and enabling convergence of nonlinear solvers for Navier-Stokes equations by continuous data assimilation

Global convergence of Newton's method for the regularized $p$-Stokes equations

Two-Level Newton Iterative Method for the 2D/3D Stationary Incompressible Magnetohydrodynamics

Global convergence of Newton’s method for the regularized p-Stokes equations

Numerical study on stability and accuracy of the fractional step algorithm for the incompressible N-S equations

Newton iterative method based on finite element discretization for the stationary Darcy–Brinkman equations

On Two-Level Oseen Iterative Methods for the 2D/3D Steady Navier–Stokes Equations

A decoupled, parallel, iterative finite element method for solving the steady boussinesq equations

A New Higher Order Fractional-Step Method for the Incompressible Navier-Stokes Equations

An iterative stabilized scheme for unsteady incompressible viscous flow

An All-Speed Asymptotic-Preserving Method for the Isentropic Euler and Navier-Stokes Equations

An Efficient Two-Step Algorithm for the Incompressible Flow Problem.

The Carleman-Newton method to globally reconstruct the initial condition for nonlinear parabolic equations

Convergence Analysis of an Implicit Fractional-Step Method for the Incompressible Navier–Stokes Equations

Computational orders of convergence of iterative methods for Richards' equation