Abstract:Context. There exists a zoo of different astrophysical fluid dynamics solvers, most of which are based on an explicit formulation and hence stability-limited to small time steps dictated by the Courant number expressing the local speed of sound. With this limitation, the modeling of low-Mach-number flows requires small time steps that introduce significant numerical diffusion, and a large amount of computational resources are needed. On the other hand, implicit methods are often developed to exclusively model the fully incompressible or 1D case since they require the construction and solution of one or more large (non)linear systems per time step, for which direct matrix inversion procedures become unacceptably slow in two or more dimensions. Aims. In this work, we present a globally implicit 3D axisymmetric Eulerian solver for the compressible Navier–Stokes equations including the energy equation using conservative formulation and a fully simultaneous approach. We use the second-order-in-time backward differentiation formula for temporal discretization as well as the κ scheme for spatial discretization. We implement different limiter functions to prohibit the occurrence of spurious oscillations in the vicinity of discontinuities. Our method resembles the well-known monotone upwind scheme for conservation laws (MUSCL). We briefly present efficient solution methods for the arising sparse and nonlinear system of equations. Methods. To deal with the nonlinearity of the Navier–Stokes equations we used a Newton iteration procedure in which the required Jacobian matrix-vector product was reconstructed with a first-order finite difference approximation to machine precision in a matrix-free way. The resulting linear system was solved either completely matrix-free with a combination of a sufficient Krylov solver and an approximate Jacobian preconditioner or semi-matrix-free with the incomplete lower upper factorization technique as a preconditioner. The latter was dependent on a standalone approximation of the Jacobian matrix, which was optionally calculated and needed solely for the purpose of preconditioning. Results. We show our method to be capable of damping sound waves and resolving shocks even at Courant numbers larger than one. Furthermore, we prove the method's ability to solve boundary value problems like the cylindrical Taylor-Couette flow (TC), including viscosity, and to model transition flows. To show the latter, we recover predicted growth rates for the vertical shear instability, while choosing a time step orders of magnitude larger than the explicit one. Finally, we verify that our method is second order in space by simulating a simplistic, stationary solar wind.

On the Parallel Implementation and Performance Study of High-Order Rosenbrock-type Implicit Runge-Kutta Methods for the FR/CPR Solutions of the Navier-Stokes Equations

Comparison of ROW, ESDIRK, and BDF2 for Unsteady Flows with the High-Order Flux Reconstruction Formulation

Parallel multilevel restricted Schwarz preconditioners for implicit simulation of subsurface flows with Peng-Robinson equation of state

High-order space–time parallel computing of the Navier–Stokes equations

Achieving high performance and portable parallel GMRES algorithm for compressible flow simulations on unstructured grids

A High-Order Flux Reconstruction/correction Procedure Via Reconstruction Formulation for Unsteady Incompressible Flow on Unstructured Moving Grids

Investigation of 3D Internal Flow Using New Flux-Reconstruction High Order Method

Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

On the accuracy of stiff-accurate diagonal implicit Runge-Kutta methods for finite volume based Navier-Stokes equations

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Implicit Block Data-Parallel Relaxation Scheme of Navier–Stokes Equations Using Graphics Processing Units

Explicit two‐step Runge‐Kutta methods for computational fluid dynamics solvers

Immersed boundary method for high-order flux reconstruction based on volume penalization

Fast and high-order approximation of parabolic equations using hierarchical direct solvers and implicit Runge-Kutta methods

A Rational High-Order Compact Difference Method for the Steady-State Stream Function-Vorticity Formulation of the Navier-Stokes Equations

Numerical Analysis of a 4Th-Order Time Parallel Algorithm for the Time-Dependent Navier-Stokes Equations

On the implementation of a robust and efficient finite element-based parallel solver for the compressible Navier-Stokes equations

Comparative study of inner-outer Krylov solvers for linear systems in structured and high-order unstructured CFD problems

A two-stage fourth-order gas-kinetic CPR method for the Navier-Stokes equations on triangular meshes

MATRICS: The implicit matrix-free Eulerian hydrodynamics solver