Dimension by Dimension Finite Volume HWENO Method for Hyperbolic Conservation Laws

In this paper,we propose a finite volume Hermite weighted essentially non-oscillatory(HWENO)method based on the dimension by dimension framework to solve hyperbolic conservation laws.It can maintain the high accuracy in the smooth region and obtain the high resolution solution when the discontinuity appears,and it is compact which will be good for giving the numerical boundary conditions.Furthermore,it avoids complicated least square procedure when we implement the genuine two dimensional(2D)finite volume HWENO reconstruction,and it can be regarded as a generalization of the one dimensional(1D)HWENO method.Extensive numerical tests are performed to verify the high resolution and high accuracy of the scheme.

On High-Resolution Entropy-Consistent Flux with Slope Limiter for Hyperbolic Conservation Laws

This paper proposes a new version of the high-resolution entropy-consistent(EC-Limited)flux for hyperbolic conservation laws based on a new minmod-type slope limiter.Firstly,we identify the numerical entropy production,a third-order differential term deduced from the previous work of Ismail and Roe[11].The corresponding dissipation term is added to the original Roe flux to achieve entropy consistency.The new,resultant entropy-consistent(EC)flux has a general and explicit analytical form without any corrective factor,making it easy to compute and a less-expensive method.The inequality constraints are imposed on the standard piece-wise quadratic reconstruction to enforce the pointwise values of bounded-type numerical solutions.We design the new minmod slope limiter as combining two separate limiters for left and right states.We propose the EC-Limited flux by adding this reconstruction data method to the primitive variables rather than to the conservative variables of the EC flux to preserve the equilibrium of the primitive variables.These resulting fluxes are easily applied to general hyperbolic conservation laws while having attractive features:entropy-stable,robust,and non-oscillatory.To illustrate the potential of these proposed fluxes,we show the applications to the Burgers equation and the Euler equations.

Application of a fourth-order relaxation scheme to hyperbolic systems of conservation laws

A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourthorder central weighted essentially nonoscillatory （CWENO） reconstruction for one-dimensional cases, which is generalized to two-dimensional cases by the dimension-by-dimension approach. The large stability domain Runge-Kutta-type solver ROCK4 is used for time integration. The resulting method requires neither the use of Riemann solvers nor the computation of Jacobians and therefore it enjoys the main advantage of the relaxation schemes. The high accuracy and high-resolution properties of the present method are demonstrated in one- and two-dimensional numerical experiments.

INTERACTION OF STRONG AND WEAK SINGULARITIES FOR HYPERBOLIC SYSTEM OF CONSERVATION LAWS IN MULTIDIMENSIONAL SPACE

In this paper we study the interaction of strong and weak singularities for hyperbolic system of conservation laws in multidimensional space. Under the assumption of transversal intersect of the shock front with the bicharacteristics bearing weak singularities we proved a theorem on regularity propagation across the shock front.

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.

High-order maximum-principle-preserving and positivity-preserving weighted compact nonlinear schemes for hyperbolic conservation laws

In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.

SPECTRAL/HP ELEMENT METHOD WITH HIERARCHICAL RECONSTRUCTION FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.

LARGE TIME STEP GENERALIZATION OF RANDOM CHOICE FINITE DIFFERENCE SCHEME FOR HYPERBOLIC CONSERVATION LAWS

A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.

A NOTE ON THE RIEMANN PROBLEM TO HYPERBOLIC CONSERVATION LAWS

This paper is a extension of [1], [3]. By the method in [1], the authors prove the global existence of the solutions of the Riemann problem for a class of 2 × 2 nonstrictly hyperbolic conservation laws.

Regularity of Fluxes in Nonlinear Hyperbolic Balance Laws

This paper addresses the issue of the formulation of weak solutions to systems of nonlinear hyperbolic conservation laws as integral balance laws.The basic idea is that the“meaningful objects”are the fluxes,evaluated across domain boundaries over time intervals.The fundamental result in this treatment is the regularity of the flux trace in the multi-dimensional setting.It implies that a weak solution indeed satisfies the balance law.In fact,it is shown that the flux is Lipschitz continuous with respect to suitable perturbations of the boundary.It should be emphasized that the weak solutions considered here need not be entropy solutions.Furthermore,the assumption imposed on the flux f(u)is quite minimal-just that it is locally bounded.

An oscillation-free Hermite WENO scheme for hyperbolic conservation laws

In this paper, the sixth-order oscillation-free Hermite weighted essentially non-oscillatory (OFHWENO) scheme is proposed for hyperbolic conservation laws on structured meshes, where the zeroth- andfirst-order moments are the variables for the governing equations. The main difference from other HWENOschemes existing in the literature is that we add high-order numerical damping terms in the first-order momentequations to control spurious oscillations for the OF-HWENO scheme. The OF-HWENO scheme not only canachieve the designed optimal numerical order, but also can be easily implemented as we use only one set ofstencils in the reconstruction procedure and the same reconstructed polynomials are applied for the zeroth- andfirst-order moment equations. In order to obtain the adaptive order resolution when facing discontinuities, atransition polynomial is added in the reconstruction, where the associated linear weights can also be any positivenumbers as long as their summation equals one. In addition, the OF-HWENO scheme still keeps compactnessas only immediate neighbor values are needed in the space discretization. Some benchmark numerical tests areperformed to illustrate the high-order accuracy, high resolution and robustness of the proposed scheme.

DECAY OF POSITIVE WAVES OF HYPERBOLIC BALANCE LAWS

Historically, decay rates have been used to provide quantitative and quali- tative information on the solutions to hyperbolic conservation laws. Quantitative results include the establishment of convergence rates for approximating procedures and numer- ical schemes. Qualitative results include the establishment of results on uniqueness and regularity as well as the ability to visualize the waves and their evolution in time. This work presents two decay estimates on the positive waves for systems of hyperbolic and gen- uinely nonlinear balance laws satisfying a dissipative mechanism. The result is obtained by employing the continuity of Glimm-type functionals and the method of generalized characteristics [7, 17, 241.

Wave Interactions for the Suliciu Relaxation System

In this paper, we investigate the elementary wave interactions for the Suliciu relaxation system and construct uniquely the solution by the characteristic analysis method in the phase plane. We find that the elementary wave interactions have a much simpler structure for the Temple class than the general systems of conservation laws. It is observed that the Riemann solutions of the Suliciu relaxation system are stable under the small perturbation on the Riemann initial data.

The Formulation of Finite Difference RBFWENO Schemes for Hyperbolic Conservation Laws:An Alternative Technique

To solve conservation laws,efficient schemes such as essentially nonoscillatory(ENO)and weighted ENO(WENO)have been introduced to control the Gibbs oscillations.Based on radial basis functions(RBFs)with the classical WENO-JS weights,a new type of WENO schemes has been proposed to solve conservation laws[J.Guo et al.,J.Sci.Comput.,70(2017),pp.551–575].The purpose of this paper is to introduce a new formulation of conservative finite difference RBFWENO schemes to solve conservation laws.Unlike the usual method for reconstructing the flux functions,the flux function is generated directly with the conservative variables.Comparing with Guo and Jung(2017),the main advantage of this framework is that arbitrary monotone fluxes can be employed,while in Guo and Jung(2017)only smooth flux splitting can be used to reconstruct flux functions.Several 1D and 2D benchmark problems are prepared to demonstrate the good performance of the new scheme.

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part Ⅱ:Entropy Stability

In the hyperbolic research community,there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability.In the first part of the series[6],the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly.By applying this technique,the authors demonstrate that a pure continu-ous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way.In this work,we extend this investigation to the nonlinear case and focus on entropy conservation.By switching to entropy variables,we provide an estimation of the boundary operators also for nonlinear problems,that guarantee conservation.In numerical simulations,we verify our theoretical analysis.

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

This paper continues to study the central relaxing schemes for system of hyperbolic conservation laws, based on the local relaxation approximation. Two classes of relaxing systems with stiff source term are introduced to approximate system of conservation laws in curvilinear coordinates. Based on them, the semi-implicit relaxing schemes are con- structed as in [6, 12] without using any linear or nonlinear Riemann solvers. Numerical experiments for one-dimensional and two-dimensional problems are presented to demon- strate the performance and resolution of the current schemes.

A Finite Volume Upwind-Biased Centred Scheme for Hyperbolic Systems of Conservation Laws:Application to Shallow Water Equations

We construct a new first-order central-upwind numerical method for solving systems of hyperbolic equations in conservative form.It applies in multidimensional structured and unstructured meshes.The proposed method is an extension of the UFORCEmethod developed by Stecca,Siviglia and Toro[25],in which the upwind bias for the modification of the staggered mesh is evaluated taking into account the smallest and largest wave of the entire Riemann fan.The proposed first-order method is shown to be identical to the Godunov upwindmethod in applications to a 2×2 linear hyperbolic system.The method is then extended to non-linear systems and its performance is assessed by solving the two-dimensional inviscid shallow water equations.Extension to second-order accuracy is carried out using an ADER-WENO approach in the finite volume framework on unstructured meshes.Finally,numerical comparison with current competing numerical methods enables us to identify the salient features of the proposed method.

A NOTE ON THE INTERACTIONS OF ELEMENTARY WAVES FOR THE AR TRAFFIC FLOW MODEL WITHOUT VACUUM

In this note, we consider the interactions of elementary waves for the traffic flow model proposed by Aw and Rascle when the vacuum is not involved. The solutions are obtained constructively and globally when the initial data consist of three pieces of constant states. Furthermore, it can be found that the Riemann solutions are stable with respect to such small perturbations of the initial data in this particular situation by investigating the limits of the solutions as the perturbed parameter ε goes to zero.

Moment-Based Multi-Resolution HWENO Scheme for Hyperbolic Conservation Laws

In this paper,a high-order moment-based multi-resolution Hermite weighted essentially non-oscillatory(HWENO)scheme is designed for hyperbolic conservation laws.The main idea of this scheme is derived from our previous work[J.Comput.Phys.,446(2021)110653],in which the integral averages of the function and its first order derivative are used to reconstruct both the function and its first order derivative values at the boundaries.However,in this paper,only the function values at the Gauss-Lobatto points in the one or two dimensional case need to be reconstructed by using the information of the zeroth and first order moments.In addition,an extra modification procedure is used to modify those first order moments in the troubledcells,which leads to an improvement of stability and an enhancement of resolution near discontinuities.To obtain the same order of accuracy,the size of the stencil required by this moment-based multi-resolution HWENO scheme is still the same as the general HWENO scheme and is more compact than the generalWENO scheme.Moreover,the linear weights are not unique and are independent of the node position,and the CFL number can still be 0.6whether for the one or two dimensional case,which has to be 0.2 in the two dimensional case for other HWENO schemes.Extensive numerical examples are given to demonstrate the stability and resolution of such moment-based multi-resolution HWENO scheme.

On the Liquid-Vapor Phase-Change Interface Conditions for Numerical Simulation of Violent Separated Flows

Numerous models have been proposed in the literature to include phase change into numerical simulations of two-phase flows.This review paper presents the modeling options that have been taken in order to obtain a model for violent separated flows with application to sloshing wave impacts.A relaxation model based on linear non-equilibrium thermodynamics has been chosen to compute the rate of phase change.The integration in the system of partial differential equations is done through a non-conservative advection term.For each of these modelling choices,some alternative models from the literature are presented and discussed.The theoretical framework for all phase change model(conservation equations and entropy growth)is also summarized.