A high order central-upwind scheme for hyperbolic conservation laws

A high order central-upwind scheme for approximating hyperbolic conservation laws is proposed. This construction is based on the evaluation of the local propagation speeds of the discontinuities and Peer's fourth order non-oscillatory reconstruction. The presented scheme shares the simplicity of central schemes, namely no Riemann solvers are involved. Furthermore, it avoids alternating between two staggered grids, which is particularly a challenge for problems which involve complex geometries and boundary conditions. Numerical experiments demonstrate the high resolution and non-oscillatory properties of our scheme.     

BV estimates of Lax-Friedrichs' scheme for a class of nonlinear hyperbolic conservation laws

We give uniform BV estimates and -stability of Lax-Friedrichs' scheme for a class of systems of strictly hyperbolic conservation laws whose integral curves of the eigenvector fields are straight lines, i.e., Temple class, under the assumption of small total variation. This implies that the approximate solutions generated via the Lax-Friedrichs' scheme converge to the solution given by the method of vanishing viscosity or the Godunov scheme, and then the Glimm scheme or the wave front tracking method.

     

Centred TVD schemes for hyperbolic conservation laws

New first- and high-order centred methods for conservation lawsare presented. Convenient TVD conditions for constructing centredTVD schemes are then formulated and some useful results areproved. Two families of centred TVD schemes are constructedand extended to nonlinear systems. Some numerical results arealso presented.     

Central ADER schemes for hyperbolic conservation laws

In this paper we first briefly review the very high order ADER methods for solving hyperbolic conservation laws. ADER methods use high order polynomial reconstruction of the solution and upwind fluxes as the building block. They use a first order upwind Godunov and the upwind second order weighted average (WAF) fluxes. As well known the upwind methods are more accurate than central schemes. However, the superior accuracy of the ADER upwind schemes comes at a cost, one must solve exactly or approximately the Riemann problems (RP). Conventional Riemann solvers are usually complex and are not available for many hyperbolic problems of practical interest. In this paper we propose to use two central fluxes, instead of upwind fluxes, as the building block in ADER scheme. These are the monotone first order Lax-Friedrich (LXF) and the third order TVD flux. The resulting schemes are called central ADER schemes. Accuracy of the new schemes is established. Numerical implementations of the new schemes are carried out on the scalar conservation laws with a linear flux, nonlinear convex flux and non-convex flux. The results demonstrate that the proposed scheme, with LXF flux, is comparable to those using first and second order upwind fluxes while the scheme, with third order TVD flux, is superior to those using upwind fluxes. When compared with the state of art ADER schemes, our central ADER schemes are faster, more accurate, Riemann solver free, very simple to implement and need less computer memory. A way to extend these schemes to general systems of nonlinear hyperbolic conservation laws in one and two dimensions is presented.     

Equilibrium schemes for scalar conservation laws with stiff sources

We consider a simple model case of stiff source terms in hyperbolic conservation laws, namely, the case of scalar conservation laws with a zeroth order source with low regularity. It is well known that a direct treatment of the source term by finite volume schemes gives unsatisfactory results for both the reduced CFL condition and refined meshes required because of the lack of accuracy on equilibrium states. The source term should be taken into account in the upwinding and discretized at the nodes of the grid. In order to solve numerically the problem, we introduce a so-called equilibrium schemes with the properties that (i) the maximum principle holds true; (ii) discrete entropy inequalities are satisfied; (iii) steady state solutions of the problem are maintained. One of the difficulties in studying the convergence is that there are no estimates for this problem. We therefore introduce a kinetic interpretation of upwinding taking into account the source terms. Based on the kinetic formulation we give a new convergence proof that only uses property (ii) in order to ensure desired compactness framework for a family of approximate solutions and that relies on minimal assumptions. The computational efficiency of our equilibrium schemes is demonstrated by numerical tests that show that, in comparison with an usual upwind scheme, the corresponding equilibrium version is far more accurate. Furthermore, numerical computations show that equilibrium schemes enable us to treat efficiently the sources with singularities and oscillating coefficients.

     

An adaptive mesh method for 1D hyperbolic conservation laws

An adaptive method is developed for solving one-dimensional systems of hyperbolic conservation laws, which combines the rezoning approach with the finite volume weighted essentially non-oscillatory (WENO) scheme. An a posteriori error estimate, used to equidistribute the mesh, is obtained from the differences between respective numerical solutions of 5th-order WENO (WENO5) and 3rd-order ENO (ENO3) schemes. The number of grids can be adaptively readjusted based on the solution structure. For higher efficiency, mesh readjustment is performed every few time steps rather than every time step. In addition, a high order conservative interpolation is used to compute the physical solutions on the new mesh from old mesh based on the finite volume ENO reconstruction. Extensive examples suggest that this adaptive method exhibits more accurate resolution of discontinuities for a similar level of computational time comparing with that on a uniform mesh.     

一种新的求解双曲型守恒律的三阶中心差分格式

本文提出了一种求解双曲型守恒律新的三阶中心差分格式，主要是引入了一种推广的三阶重构，并证明了这种重构在网格边界无振荡．所提的格式保持了中心差分格式简单的优点，不需用Riemann解算器，避免了进行特征解耦．数值试验结果表明本文格式是高精度、高分辨率的。     

BV-estimates of Lax–Friedrichs' scheme for hyperbolic conservation laws with relaxation

In this paper, we will give BV-estimates of Lax–Friedrichs' scheme for a simple hyperbolic system of conservation laws with relaxation and get the global existence and uniqueness of BV-solution by the BV-estimates above. Furthermore, our results show that the solution converge towards the solution of an equilibrium model as the relaxation time ε>0 tends to zero provided sub-characteristic condition holds. Copyright © 2007 John Wiley & Sons, Ltd.     

An adaptive wavelet viscosity method for hyperbolic conservation laws

We extend the multiscale finite element viscosity method for hyperbolic conservation laws developed in terms of hierarchical finite element bases to a (pre‐orthogonal spline‐)wavelet basis. Depending on an appropriate error criterion, the multiscale framework allows for a controlled adaptive resolution of discontinuities of the solution. The nonlinearity in the weak form is treated by solving a least‐squares data fitting problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Systems of hyperbolic conservation laws with a resonant moving source

In this paper, we study the stability of a single transonic shock wave solution to the hyperbolic conservation laws with a resonant moving source. Compared with the previous results [W.-C. Lien, Hyperbolic conservation laws with a moving source, Comm. Pure Appl. Math. 52 (9) (1999) 1075-1098; T.P. Liu, Nonlinear stability and instability of transonic flows through a nozzle, Comm. Math. Phys. 83 (2) (1982) 243-260] on this stability problem, in this paper, the transonic ith shock is assumed to be relatively strong and stable in the sense of Majda. Then the framework of [M. Lewicka, L¹ stability of patterns of non-interacting large shock waves, Indiana Univ. Math. J. 49 (4) (2000) 1515-1537; M. Lewicka, Stability conditions for patterns of noninteracting large shock waves, SIAM J. Math. Anal. 32 (5) (2001) 1094-1116 (electronic)] can be applied. A new criterion is obtained to test whether such a shock is time asymptotically stable or not. And by constructing the Liu-Yang functional, one can prove the L¹ stability of the shock under the stability condition. This is an extension of the result [S.-Y. Ha, T. Yang, L¹ stability for systems of hyperbolic conservation laws with a resonant moving source, SIAM J. Math. Anal. 34 (5) (2003) 1226-1251 (electronic); W.-C. Lien, Hyperbolic conservation laws with a moving source, Comm. Pure Appl. Math. 52 (9) (1999) 1075-1098] to a more general case.     

High-order accurate monotone compact running scheme for multidimensional hyperbolic equations

Monotone absolutely stable conservative difference schemes intended for solving quasilinear multidimensional hyperbolic equations are described. For sufficiently smooth solutions, the schemes are fourth-order accurate in each spatial direction and can be used in a wide range of local Courant numbers. The order of accuracy in time varies from the third for the smooth parts of the solution to the first near discontinuities. This is achieved by choosing special weighting coefficients that depend locally on the solution. The presented schemes are numerically efficient thanks to the simple two-diagonal (or block two-diagonal) structure of the matrix to be inverted. First the schemes are applied to system of nonlinear multidimensional conservation laws. The choice of optimal weighting coefficients for the schemes of variable order of accuracy in time and flux splitting is discussed in detail. The capabilities of the schemes are demonstrated by computing well-known two-dimensional Riemann problems for gasdynamic equations with a complex shock wave structure.     

Relaxation WENO schemes for multidimensional hyperbolic systems of conservation laws

We present a class of high‐order weighted essentially nonoscillatory (WENO) reconstructions based on relaxation approximation of hyperbolic systems of conservation laws. The main advantage of combining the WENO schemes with relaxation approximation is the fact that the presented schemes avoid solution of the Riemann problems due to the relaxation approach and high‐resolution is obtained by applying the WENO approach. The emphasis is on a fifth‐order scheme and its performance for solving a wide class of systems of conservation laws. To show the effectiveness of these methods, we present numerical results for different test problems on multidimensional hyperbolic systems of conservation laws. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

An efficient TVD‐WENO method for conservation laws

In this article, we present a high‐resolution hybrid scheme for solving hyperbolic conservation laws in one and two dimensions. In this scheme, we use a cheap fourth order total variation diminishing (TVD) scheme for smooth region and expensive seventh order weighted nonoscillatory (WENO) scheme near discontinuities. To distinguish between the smooth parts and discontinuities, we use an efficient adaptive multiresolution technique. For time integration, we use the third order TVD Runge‐Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

An efficient WENO scheme for solving hyperbolic conservation laws

In this paper we propose a new WENO scheme, in which we use a central WENO [G. Capdeville, J. Comput. Phys. 227 (2008) 2977-3014] (CWENO) reconstruction combined with the smoothness indicators introduced in [R. Borges, M. Carmona, B. Costa, W. Sun Don, J. Comput. Phys. 227 (2008) 3191-3211] (IWENO). We use the central-upwind flux [A. Kurganov, S. Noelle, G. Petrova, SIAM J. Sci. Comp. 23 (2001) 707-740] which is simple, universal and efficient. For time integration we use the third order TVD Runge-Kutta scheme. The resulting scheme improves the convergence order at critical points of smooth parts of solution as well as decrease the dissipation near discontinuities. Numerical experiments of the new scheme for one and two-dimensional problems are reported. The results demonstrates that the proposed scheme is superior to the original CWENO and IWENO schemes.     

An optimal error estimate for upwind Finite Volume methods for nonlinear hyperbolic conservation laws

The purpose of this paper is to show that the cell-centered upwind Finite Volume scheme applied to general hyperbolic systems of m conservation laws approximates smooth solutions to the continuous problem at order one in space and time. As it is now well understood, there is a lack of consistency for order one upwind Finite Volume schemes: the truncation error does not tend to zero as the time step and the grid size tend to zero. Here, following our previous papers on scalar equations, we construct a corrector that allows us to prove the expected error estimate for nonlinear systems of equations in one dimension.     

Viscosity methods for piecewise smooth solutions to scalar conservation laws

It is proved that for scalar conservation laws, if the flux function is strictly convex, and if the entropy solution is piecewise smooth with finitely many discontinuities (which includes initial central rarefaction waves, initial shocks, possible spontaneous formation of shocks in a future time and interactions of all these patterns), then the error of viscosity solution to the inviscid solution is bounded by in the -norm, which is an improvement of the upper bound. If neither central rarefaction waves nor spontaneous shocks occur, the error bound is improved to .

     

On the optimization of flux limiter schemes for hyperbolic conservation laws

A classic strategy to obtain high‐quality discretizations of hyperbolic partial differential equations is to use flux limiter (FL) functions for blending two types of approximations: a monotone first‐order scheme that deals with discontinuous solution features and a higher order method for approximating smooth solution parts. In this article, we study a new approach to FL methods. Relying on a classification of input data with respect to smoothness, we associate specific basis functions with the individual smoothness notions. Then, we construct a limiter as a linear combination of the members of parameter‐dependent families of basis functions, and we explore the possibility to optimize the parameters in interesting model situations to find a corresponding optimal limiter. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

BV solutions for hyperbolic systems of balance laws with relaxation

For hyperbolic systems of balance laws with source manifesting relaxation, it is shown that the Kawashima condition, which yields global classical solutions with smooth initial values near equilibrium, is also instrumental in inducing the existence of global admissible BV solutions, accommodating shocks.