Hyperbolic Systems Of Conservation Laws Research Articles

KFVM-WENO: A High-order Accurate Kernel-based Finite Volume Method for Compressible Hydrodynamics

This paper presents a fully multidimensional kernel-based reconstruction scheme for finite volume methods applied to systems of hyperbolic conservation laws, with a particular emphasis on the compressible Euler equations. Nonoscillatory reconstruction is achieved through an adaptive-order weighted essentially nonoscillatory (WENO) method cast into a form suited to multidimensional reconstruction. A kernel-based approach inspired by radial basis functions and Gaussian process modeling, which we call kernel-based finite volume method with WENO, is presented here. This approach allows the creation of a scheme of arbitrary order of accuracy with simply defined multidimensional stencils and substencils. Furthermore, the fully multidimensional nature of the reconstruction allows for a more straightforward extension to higher spatial dimensions and removes the need for complicated boundary conditions on intermediate quantities in modified dimension-by-dimension methods. In addition, a new simple yet effective set of reconstruction variables is introduced, which could be useful in existing schemes with little modification. The proposed scheme is applied to a suite of stringent and informative benchmark problems to demonstrate its efficacy and utility. A highly parallel multi-GPU implementation using Kokkos and the message-passing interface is also provided.

Existence of a global smooth solution for a degenerate problem of the hyperbolic conservation law systems

In this paper, we consider the existence of a global smooth solution for a degenerate problem of the hyperbolic conservation law systems. We use the characteristic decomposition to establish the global existence of smooth solutions to the degenerate Goursat problem up to a boundary and prove that the vacuum bubble is a point. This problem comes from the expansion of a wedge of gas with constant velocity into vacuum in gas dynamics.

Erratum to : “A family of second-order dissipative finite volume schemes for hyperbolic systems of conservation laws”

Erratum to : “A family of second-order dissipative finite volume schemes for hyperbolic systems of conservation laws”

An almost fail-safe a-posteriori limited high-order CAT scheme

In this paper we blend the high order Compact Approximate Taylor (CAT) numerical schemes with an a posteriori Multi-dimensional Optimal Order Detection (MOOD) paradigm to solve hyperbolic systems of conservation laws in 2D. The resulting scheme presents high accuracy on smooth solutions, essentially non-oscillatory behavior on irregular ones, and almost fail-safe property concerning positivity issues. The numerical results on a set of sanity test cases and demanding ones are presented assessing the appropriate behavior of the CAT-MOOD scheme.

Riemann problem for a 2times 2 hyperbolic system with time-gradually-degenerate damping

This paper is focused on the Riemann problem for a 2times 2 hyperbolic system of conservation laws with a time-gradually-degenerate damping. Two kinds of non-self-similar solutions involving the delta-shocks and vacuum are obtained using the variable substitution method. The generalized Rankine-Hugoniot relation and entropy condition are clarified for the delta-shock. Furthermore, the vanishing viscosity method proves the existence, uniqueness, and stability of non-self-similar solutions.

The mathematical model and analysis of the nanoparticle-stabilized foam displacement

This work proposes a mathematical model to study the foam displacement in porous media stabilized by nanoparticles. We consider a simplification of the Stochastic Bubble Population balance model in local equilibrium, with nanoparticle dependence inspired by the experimental data from the literature. It consists of a non-strictly hyperbolic system of conservation laws, which is solved for the generic initial and injection conditions. We investigate the existence of a global solution as a sequence of waves following the Conservation Laws Theory. When the solution is composed of two or more waves, we present necessary and sufficient conditions to guarantee the compatibility of these wave sequences. The analytical solution for the nanoparticle-stabilized foam displacement in porous media allowed us to quantify the effect of nanoparticles on foam displacement, focusing on the breakthrough time and cumulative water production. In agreement with the literature, when only gas is injected, the breakthrough time and the water production increase with the nanoparticle concentration. Although, we also observe that the effect of nanoparticles is less pronounced for high nanoparticle concentration. Counterintuitively, during gas-water co-injection for a certain parameter range, adding nanoparticles changes the mathematical solution qualitatively, yielding a negligible effect on water production. We discuss the most favorable conditions to observe the action of nanoparticles in laboratory experiments.

The Exact Boundary Controllability of Nodal Profile for Entropy Solutions to 1-D Quasilinear Hyperbolic Systems of Conservation Laws with Linearly Degenerate Negative (Resp., Positive) Characteristic Fields

The Exact Boundary Controllability of Nodal Profile for Entropy Solutions to 1-D Quasilinear Hyperbolic Systems of Conservation Laws with Linearly Degenerate Negative (Resp., Positive) Characteristic Fields

Existence of multi‐dimensional contact discontinuities for the ideal compressible magnetohydrodynamics

AbstractWe establish the local existence and uniqueness of multi‐dimensional contact discontinuities for the ideal compressible magnetohydrodynamics (MHD) in Sobolev spaces, which are most typical interfacial waves for astrophysical plasmas and prototypical fundamental waves for hyperbolic systems of conservation laws. Such waves are characteristic discontinuities for which there is no flow across the discontinuity surface while the magnetic field crosses transversely, which lead to a two‐phase free boundary problem where the pressure, velocity and magnetic field are continuous across the interface whereas the entropy and density may have jumps. To overcome the difficulties of possible nonlinear Rayleigh–Taylor instability and loss of derivatives, here we use crucially the Lagrangian formulation and Cauchy's celebrated integral (1815) for the magnetic field. These motivate us to define two special good unknowns; one enables us to capture the boundary regularizing effect of the transversal magnetic field on the flow map, and the other one allows us to get around the troublesome boundary integrals due to the transversality of the magnetic field. In particular, our result removes the additional assumption of the Rayleigh–Taylor sign condition required by Morando, Trakhinin and Trebeschi (J. Differ. Equ. 258 (2015), no. 7, 2531–2571; Arch. Ration. Mech. Anal. 228 (2018), no. 2, 697–742) and holds for both 2D and 3D and hence gives a complete answer to the two open questions raised therein. Moreover, there is no loss of derivatives in our well‐posedness theory. The solution is constructed as the inviscid limit of solutions to some suitably‐chosen nonlinear approximate problems for the two‐phase compressible viscous non‐resistive MHD.

A triangle-based positive semi-discrete Lagrangian–Eulerian scheme via the weak asymptotic method for scalar equations and systems of hyperbolic conservation laws

In this paper, we show how to construct a positive semi-discrete Lagrangian–Eulerian scheme on triangular grids that asymptotically satisfies multidimensional scalar hyperbolic problems and first-order nonlinear systems of conservation laws. The construction of such numerical scheme is based on the improved concept of space–time no-flow curves. For scalar equations, through a sequence of computational grid functions that tend to satisfy the underlying hyperbolic conservation law(s) in a weak sense in the space variable and in a strong sense in the time variable, we prove that the proposed scheme converges to the unique weak solution satisfying a type of Kruzhkov entropy condition. For nonlinear systems of first-order hyperbolic conservation laws, we show that the proposed scheme also satisfies a (weak) positivity principle. This provides an effective numerical analysis tool (based on weak asymptotic methods) to construct a family of continuous functions that asymptotically satisfies scalar conservation laws with discontinuous nonlinearity and systems that have irregular solutions. In order to assess the robustness of the proposed approach, we solve a series of numerical tests on several hyperbolic problems (scalar and systems) using the positive semi-discrete Lagrangian–Eulerian scheme on triangular grids. The results confirm that the proposed method is efficient (since it is Riemann-solver-free) and provides good resolution for scalar models and for the solution of systems of conservation laws.

Experimental convergence rate study for three shock‐capturing schemes and development of highly accurate combined schemes

AbstractWe study experimental convergence rates of three shock‐capturing schemes for hyperbolic systems of conservation laws: the second‐order central‐upwind (CU) scheme, the third‐order Rusanov‐Burstein‐Mirin (RBM), and the fifth‐order alternative weighted essentially non‐oscillatory (A‐WENO) scheme. We use three imbedded grids to define the experimental pointwise, integral, and convergence rates. We apply the studied schemes to the shallow water equations and conduct their comprehensive numerical convergence study. We verify that while the studied schemes achieve their formal orders of accuracy on smooth solutions, after the shock formation, a part of the computed solutions is affected by shock propagation and both the pointwise and integral convergence rates reduce there. Moreover, while the convergence rates for the CU and A‐WENO schemes, which rely on nonlinear stabilization mechanisms, reduce to the first order, the RBM scheme, which utilizes a linear stabilization, is clearly second‐order accurate. Finally, relying on the conducted experimental convergence rate study, we develop two new combined schemes based on the RBM and either the CU or A‐WENO scheme. The obtained combined schemes can achieve the same high order of accuracy as the RBM scheme in the smooth areas while being non‐oscillatory near the shocks.

A New Entropy Stable Finite Difference Scheme for Hyperbolic Systems of Conservation Laws

The hyperbolic problem has a unique entropy solution, which maintains the entropy inequality. As such, people hope that the numerical results should maintain the discrete entropy inequalities accordingly. In view of this, people tend to construct entropy stable (ES) schemes. However, traditional numerical schemes cannot directly maintain discrete entropy inequalities. To address this, we here construct an ES finite difference scheme for the nonlinear hyperbolic systems of conservation laws. The proposed scheme can not only maintain the discrete entropy inequality, but also enjoy high-order accuracy. Firstly, we construct the second-order accurate semi-discrete entropy conservative (EC) schemes and ensure that the schemes meet the entropy identity when an entropy pair is given. Then, the second-order EC schemes are used as a building block to achieve the high-order accurate semi-discrete EC schemes. Thirdly, we add a dissipation term to the above schemes to obtain the high-order ES schemes. The term is based on the Weighted Essentially Non-Oscillatory (WENO) reconstruction. Finally, we integrate the scheme using the third-order Runge–Kutta (RK) approach in time. In the end, plentiful one- and two-dimensional examples are implemented to validate the capability of the scheme. In summary, the current scheme has sharp discontinuity transitions and keeps the genuine high-order accuracy for smooth solutions. Compared to the standard WENO schemes, the current scheme can achieve higher resolution.

Tightly coupled hyperbolic treatment of buoyant two-phase flow and transport in porous media

Buoyant incompressible multiphase flow and transport in porous media typically is simulated by coupling an elliptic flow equation with a hyperbolic transport description. The tight coupling between flow and transport either requires a fully implicit solution algorithm or very small time steps, if solved sequentially. In any case, however, a large linear system has to be solved each time step.Here a new solution approach is devised, which relies on solving a coupled hyperbolic system of conservation laws with an explicit finite volume method. Consequently, only local operations have to be performed, which is a great advantage in case of massive parallel simulations and if GPUs are employed.The devised method is based on the isothermal Euler equations with momentum source terms accounting for resistance due to the porous medium and buoyancy. In this system the pressure is proportional to the density and in case of very small Mach numbers and relative density variation, the computed velocity field is divergence free and converges towards the mean total volumetric flow in the porous domain. To account for saturation transport, the system was augmented by an additional hyperbolic equation. In order to obtain the numerical fluxes, a characteristic based approximative Riemann solver was developed and 2nd order accuracy in space and time is achieved by piecewise linear reconstruction.Two-phase flow test cases with buoyant plumes and a lock exchange problem demonstrate that the devised hyperbolic approach recovers the correct incompressible solutions. Numerical experiments also demonstrate that the method is accurate and efficient.Thus the devised method is very promising, and conceptually it also can be applied for flow and transport in heterogeneous media; but it remains to be investigated how well it performs for such cases.

On the Consistency of Arnoldi‐Based Krylov Methods for Conservation Laws

AbstractConservation and consistency are fundamental properties of discretizations of systems of hyperbolic conservation laws. Re‐ cently, these concepts have been extended to the realm of iterative methods by defining locally conservative and flux consistent iterations. In this note, the current status of such iterative methods is summarized. In particular, it has been shown that Krylov subspace methods are locally conservative, but that they are not flux consistent. Here, we approach the problem of quantifying the flux inconsistency of Krylov subspace methods. Krylov methods introduce a time retardation factor into discretizations of linear conservation laws. It has thusfar been unknown how to compute the precise value of this factor. This issue is resolved herein for Arnoldi‐based Krylov subspace methods.

Global entropy stability for a class of unlimited second-order schemes for 2D hyperbolic systems of conservation laws on unstructured meshes

We introduce a class of second-order accurate finite volume schemes defined on two dimensional unstructured meshes. They satisfy a global discrete entropy stability for an arbitrary strictly convex entropy provided the second-order discretization verifies a dissipative inequality. These schemes do not require usual limitation techniques. Numerical experiments are performed to confirm the accuracy and the shock capturing efficiency of the schemes.

Stability of Riemann solutions to a class of non-strictly hyperbolic systems of conservation laws

This paper is concerned with Riemann problem for a class of non-strictly hyperbolic systems of conservation laws. The two kinds of Riemann solutions including vacuum and delta shock wave are constructed. Under the generalized Rankine-Hugoniot relation and entropy condition, we establish the existence and uniqueness of delta shock wave solutions. Furthermore, by studying the interactions among of the delta shock wave and vacuum as well as contact discontinuity, the Riemann solutions with four kinds of different structures are obtained. Additionally, the stability of the Riemann solutions is obtained under certain perturbation of the initial data.

Divergence-Preserving MHD Solutions Using a Noniterative Hyperbolic Conservation Law Scheme

The divergence constraint, <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\nabla$</tex-math> </inline-formula> <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\cdot$</tex-math> </inline-formula> B <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$=$</tex-math> </inline-formula> 0, presents a central challenging numerical feature of electromagnetic and magnetohydrodynamic solvers. In general, this constraint is not strictly enforced in hyperbolic solvers and tends to accumulate error over the duration of a simulation, polluting the solution with unphysical behavior. In this article, we present a numerical scheme that constructs a system of hyperbolic conservation laws that already satisfy the divergence constraint, starting from the potential form of Maxwell’s equations and basic derivative relationships. The new system consists only of the first-order spatial differentiation, in contrast to the potential equations which are second order; this allows them to be directly integrated into the existing hyperbolic solvers with no additional implementation changes needed. Furthermore, since the system of conservation laws constructed automatically satisfies the divergence criteria, no subiteration is required to resolve the divergence-free solution. We present the results of integrating this new conservation law approach in an existing finite volume plasmadynamic solver that discretizes the coupled Maxwell and Navier–Stokes equations. Numerical tests indicate that the new approach can strongly resolve physical solutions. The divergence error is also measured and quantified for these problems, showing that the new system suppresses divergence error by many orders of magnitude.

A family of second-order dissipative finite volume schemes for hyperbolic systems of conservation laws

We propose and study a family of formally second-order accurate schemes to approximate weak solutions of hyperbolic systems of conservation laws. Theses schemes are based on a dissipative property satisfied by the second-order discretization in space. They are proven to satisfy a global entropy inequality for a generic strictly convex entropy. These schemes do not involve limitation techniques. Numerical results are provided to illustrate their accuracy and stability.

T5 configurations and hyperbolic systems

In this paper, we study the rank-one convex hull of a differential inclusion associated to entropy solutions of a hyperbolic system of conservation laws. This was introduced in [B. Kirchheim, S. Müller and V. Šverák, Studying Nonlinear PDE by Geometry in Matrix Space (Springer, 2003), Sec. 7], and many of its properties have already been shown in [A. Lorent and G. Peng, Null Lagrangian measures in subspaces, compensated compactness and conservation laws, Arch. Ration. Mech. Anal. 234(2) (2019) 857–910; A. Lorent and G. Peng, On the Rank-1 convex hull of a set arising from a hyperbolic system of Lagrangian elasticity, Calc. Var. Partial Differential Equations 59(5) (2020) 156]. In particular, in [A. Lorent and G. Peng, On the Rank-1 convex hull of a set arising from a hyperbolic system of Lagrangian elasticity, Calc. Var. Partial Differential Equations 59(5) (2020) 156], it is shown that the differential inclusion does not contain any [Formula: see text] configurations. Here, we continue that study by showing that the differential inclusion does not contain [Formula: see text] configurations.

“A posteriori” limited high order and robust schemes for transient simulations of fluid flows in gas dynamics

In this paper, we propose a novel approximation strategy for time-dependent hyperbolic systems of conservation laws for the Euler system of gas dynamics that aims to represent the dynamics of strong interacting discontinuities. The goal of our method is to allow an approximation with a high-order of accuracy in smooth regions of the flow, while ensuring robustness and a non-oscillatory behaviour in the regions of steep gradients, in particular across shocks.Following the Multidimensional Optimal Order Detection (MOOD) ([15,17]) approach, a candidate solution is computed at a next time level via a high-order accurate explicit scheme ([3,5]). A so-called detector determines if the candidate solution reveals any spurious oscillations or numerical issue and, if so, only the troubled cells are locally recomputed via a more dissipative scheme. This allows to design a family of “a posteriori” limited, robust and positivity preserving, as well as high accurate, non-oscillatory and effective scheme. Among the detecting criteria of the novel MOOD strategy, two different approaches from literature, based on the work of [15,17] and on [35], are investigated. Numerical examples in 1D and 2D, on structured and unstructured meshes, are proposed to assess the effective order of accuracy for smooth flows, the non-oscillatory behaviour on shocked flows, the robustness and positivity preservation on more extreme flows.

A MOOD-like compact high order finite volume scheme with adaptive mesh refinement

In this paper, a novel Finite Volume (FV) scheme for obtaining high order approximations of solutions of multi-dimensional hyperbolic systems of conservation laws within an Adaptive Mesh Refinement framework is proposed. It is based on a point-wise polynomial reconstruction that avoids the recalculation of reconstruction stencils and matrices whenever a mesh is refined or coarsened. It also couples both the limiting of the FV scheme and the refinement procedure, taking advantage of the Multi-dimensional Optimal Order Detection (MOOD) detection criteria. The resulting computational procedure is employed to simulate test cases of increasing difficulty using two models of Partial Differential Equations: the Euler system and the radiative M1 model, thus demonstrating its efficiency.