Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme

TVD differencing on three-dimensional unstructured meshes with monotonicity-preserving correction of mesh skewness

This data set contains the data accompanying the article F. Denner and B. van Wachem, TVD differencing on three-dimensional unstructured meshes with monotonicity-preserving correction of mesh skewness, Journal of Computational Physics (2015), http://dx.doi.org/10.1016/j.jcp.2015.06.008.This data set contains the data accompanying the article F. Denner and B. van Wachem, TVD differencing on three-dimensional unstructured meshes with monotonicity-preserving correction of mesh skewness, Journal of Computational Physics (2015), http://dx.doi.org/10.1016/j.jcp.2015.06.008

A Moving Frame Algorithm for High Mach Number Hydrodynamics

We present a new approach to Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in astrophysical hydrodynamic simulations. The Eulerian fluid conservation equations are solved in an adaptive frame moving with the fluid where Mach numbers are minimized. The moving frame approach uses a velocity decomposition technique to define local kinetic variables while storing the bulk kinetic components in a smoothed background velocity field that is associated with the grid velocity. Gravitationally induced accelerations are added to the grid, thereby minimizing the spurious heating problem encountered in cold gas flows. Separately tracking local and bulk flow components allows thermodynamic variables to be accurately calculated in both subsonic and supersonic regions. A main feature of the algorithm, that is not possible in previous Eulerian implementations, is the ability to resolve shocks and prevent spurious heating where both the preshock and postshock Mach numbers are high. The hybrid algorithm combines the high resolution shock capturing ability of the second-order accurate Eulerian TVD scheme with a low-diffusion Lagrangian advection scheme. We have implemented a cosmological code where the hydrodynamic evolution of the baryons is captured using the moving frame algorithm while the gravitational evolution of the collisionless dark matter is tracked using a particle-mesh N-body algorithm. The MACH code is highly suited for simulating the evolution of the IGM where accurate thermodynamic evolution is needed for studies of the Lyman alpha forest, the Sunyaev-Zeldovich effect, and the X-ray background. Hydrodynamic and cosmological tests are described and results presented. The current code is fast, memory-friendly, and parallelized for shared-memory machines.Comment: 19 pages, 5 figure

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities

A rarefaction-tracking method for hyperbolic conservation laws

We present a numerical method for scalar conservation laws in one space dimension. The solution is approximated by local similarity solutions. While many commonly used approaches are based on shocks, the presented method uses rarefaction and compression waves. The solution is represented by particles that carry function values and move according to the method of characteristics. Between two neighboring particles, an interpolation is defined by an analytical similarity solution of the conservation law. An interaction of particles represents a collision of characteristics. The resulting shock is resolved by merging particles so that the total area under the function is conserved. The method is variation diminishing, nevertheless, it has no numerical dissipation away from shocks. Although shocks are not explicitly tracked, they can be located accurately. We present numerical examples, and outline specific applications and extensions of the approach.Comment: 21 pages, 7 figures. Similarity 2008 conference proceeding

Relations between WENO3 and Third-order Limiting in Finite Volume Methods

Weighted essentially non-oscillatory (WENO) and finite volume (FV) methods employ different philosophies in their way to perform limiting. We show that a generalized view on limiter functions, which considers a two-dimensional, rather than a one-dimensional dependence on the slopes in neighboring cells, allows to write WENO3 and $3^\text{rd}$-order FV schemes in the same fashion. Within this framework, it becomes apparent that the classical approach of FV limiters to only consider ratios of the slopes in neighboring cells, is overly restrictive. The hope of this new perspective is to establish new connections between WENO3 and FV limiter functions, which may give rise to improvements for the limiting behavior in both approaches.Comment: 22 page

Relativistic MHD and black hole excision: Formulation and initial tests

A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations with AMR. The fluid equations are solved using a finite difference Convex ENO method. Excision is implemented using overlapping grids. Elliptic and hyperbolic divergence cleaning techniques allow for maximum flexibility in choosing coordinate systems, and we compare both methods for a standard problem. Numerical results of standard test problems are presented in two-dimensional flat space using excision, overlapping grids, and elliptic and hyperbolic divergence cleaning.Comment: 22 pages, 8 figure