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.     

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.

     

A posteriori error estimates for upwind finite volume schemes for nonlinear conservation laws in multi dimensions

In this paper we shall derive a posteriori error estimates in the -norm for upwind finite volume schemes for the discretization of nonlinear conservation laws on unstructured grids in multi dimensions. This result is mainly based on some fundamental a priori error estimates published in a recent paper by C. Chainais-Hillairet. The theoretical results are confirmed by numerical experiments.

     

On the reservoir technique convergence for nonlinear hyperbolic conservation laws - Part I

This paper is devoted to the analysis of flux schemes coupled with the reservoir technique for approximating hyperbolic equations and linear hyperbolic systems of conservation laws [F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, The reservoir scheme for systems of conservation laws, in: Finite Volumes for Complex Applications, III, Porquerolles, 2002, Lab. Anal. Topol. Probab. CNRS, Marseille, 2002, pp. 247-254 (electronic); F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, Un procédé de réduction de la diffusion numérique des schémas à différence de flux d'ordre un pour les systèmes hyperboliques non linéaires, C. R. Math. Acad. Sci. Paris 335 (7) (2002) 627-632; F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, The reservoir technique: A way to make Godunov-type schemes zero or very low diffusive. Application to Colella-Glaz, Eur. J. Mech. B Fluids 27 (6) (2008)]. We prove the long time convergence of the reservoir technique and its TVD property for some specific but still general configurations. Proofs are based on a precise study of the treatment by the reservoir technique of shock and rarefaction waves.     

A composite semi-conservative scheme for hyperbolic conservation laws

In this work a first order accurate semi-conservative composite scheme is presented for hyperbolic conservation laws. The idea is to consider the non-conservative form of conservation law and utilize the explicit wave propagation direction to construct semi-conservative upwind scheme. This method captures the shock waves exactly with less numerical dissipation but generates unphysical rarefaction shocks in case of expansion waves with sonic points. It shows less dissipative nature of constructed scheme. In order to overcome it, we use the strategy of composite schemes. A very simple criteria based on wave speed direction is given to decide the iterations. The proposed method is applied to a variety of test problems and numerical results show accurate shock capturing and higher resolution for rarefaction fan.     

Hyperbolic conservation laws on manifolds. An error estimate for finite volume schemes

Following Ben-Artzi and LeFloch, we consider nonlinear hyperbolic conservation laws posed on a Riemannian manifold, and we establish an L1-error estimate for a class of finite volume schemes allowing for the approximation of entropy solutions to the initial value problem. The error in the L1 norm is of order h1/4 at most, where h represents the maximal diameter of elements in the family of geodesic triangulations. The proof relies on a suitable generalization of Cockburn, Coquel, and LeFloch＇s theory which was originally developed in the Euclidian setting. We extend the arguments to curved manifolds, by taking into account the effects to the geometry and overcoming several new technical difficulties.     

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     

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 .

     

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.     

Convergence of a level‐set algorithm in scalar conservation laws

This article is concerned with the convergence of the level‐set algorithm introduced by Aslam (J Comput Phys 167 (2001), 413–438) for tracking the discontinuities in scalar conservation laws in the case of linear or strictly convex flux function. The numerical method is deduced by the level‐set representation of the entropy solution: the zero of a level‐set function is used as an indicator of the discontinuity curves and two auxiliary states, which are assumed continuous through the discontinuities, are introduced. We rewrite the numerical level‐set algorithm as a procedure consisting of three big steps: (a) initialization, (b) evolution, and (c) reconstruction. In (a), we choose an entropy admissible level‐set representation of the initial condition. In (b), for each iteration step, we solve an uncoupled system of three equations and select the entropy admissible level‐set representation of the solution profile at the end of the time iteration. In (c), we reconstruct the entropy solution using the level‐set representation. We prove the convergence of the numerical solution to the entropy solution in for every , using ‐weak bounded variation (BV) estimates and a cell entropy inequality. In addition, some numerical examples focused on the elementary wave interaction are presented. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1310–1343, 2015     

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     

Stratified solutions for systems of conservation laws

We study a class of weak solutions to hyperbolic systems of conservation (balance) laws in one space dimension, called stratified solutions. These solutions are bounded and ``regular' in the direction of a linearly degenerate characteristic field of the system, but not in other directions. In particular, they are not required to have finite total variation. We prove some results of local existence and uniqueness.

     

非线性守恒律高阶谱粘性法的收敛性

讨论守恒型方程周期边界问题的高阶谱粘性方法逼近解的收敛性.在逼近解一致有界的假设下,通过建立其高阶导数的上界估计,证明了高阶谱粘性方法逼近解具有同二阶谱粘性方法逼近解相类似的高频衰减性质.以此为基础,用补偿列紧法证明了高阶谱粘性方法逼近解收敛于守恒型方程的物理解.     

On the design of an upwind scheme for compressible flow on general triangulations

We discuss the design features and mathematical background of an explicit upwind finite-volume method to simulate non-stationary flow of a compressible, inviscid fluid. One of the design goals was the rigorous mathematical justification of each ingredient of the method. The method itself contains elements from finite-difference methods as well as finite-element methods and is formulated in a finite volume framework. The use of well-known algorithmic ingredients in a new framework results in a robust time-accurate scheme. To be able to easily handle complex geometries as well as adaption algorithms a tringale-based formulation was chosen. Numerical tests for two-dimensional flow are presented.     

Error estimate for the approximation of nonlinear conservation laws on bounded domains by the finite volume method

In this paper we derive a priori and a posteriori error estimates for cell centered finite volume approximations of nonlinear conservation laws on polygonal bounded domains. Numerical experiments show the applicability of the a posteriori result for the derivation of local adaptive solution strategies.

     

A posteriori error estimator based on derivative recovery for the discontinuous Galerkin method for nonlinear hyperbolic conservation laws on Cartesian grids

In this article, we develop and analyze a new recovery‐based a posteriori error estimator for the discontinuous Galerkin (DG) method for nonlinear hyperbolic conservation laws on Cartesian grids, when the upwind flux is used. We prove, under some suitable initial and boundary discretizations, that the ‐norm of the solution is of order , when tensor product polynomials of degree at most are used. We further propose a very simple derivative recovery formula which gives a superconvergent approximation to the directional derivative. The order of convergence is showed to be . We use our derivative recovery result to develop a robust recovery‐type a posteriori error estimator for the directional derivative approximation which is based on an enhanced recovery technique. The proposed error estimators of the recovery‐type are easy to implement, computationally simple, asymptotically exact, and are useful in adaptive computations. Finally, we show that the proposed recovery‐type a posteriori error estimates, at a fixed time, converge to the true errors in the ‐norm under mesh refinement. The order of convergence is proved to be . Our theoretical results are valid for piecewise polynomials of degree and under the condition that each component, , of the flux function possesses a uniform positive lower bound. Several numerical examples are provided to support our theoretical results and to show the effectiveness of our recovery‐based a posteriori error estimator. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1224–1265, 2017     

Digital total variation filtering as postprocessing for Chebyshev pseudospectral methods for conservation laws

Digital total variation filtering is analyzed as a fast, robust, post-processing method for accelerating the convergence of pseudospectral approximations that have been contaminated by Gibbs oscillations. The method, which originated in image processing, can be combined with spectral filters to quickly post-process large data sets with sharp resolution of discontinuities and with exponential accuracy away from the discontinuities.     

Delta shock waves with Dirac delta function in both components for systems of conservation laws

We study a class of non-strictly and weakly hyperbolic systems of conservation laws which contain the equations of geometrical optics as a prototype. The Riemann problems are constructively solved. The Riemann solutions include two kinds of interesting structures. One involves a cavitation where both state variables tend to zero forming a singularity, the other is a delta shock wave in which both state variables contain Dirac delta function simultaneously. The generalized Rankine–Hugoniot relation and entropy condition are proposed to solve the delta shock wave. Moreover, with the limiting viscosity approach, we show all of the existence, uniqueness and stability of solution involving the delta shock wave. The generalized Rankine–Hugoniot relation is also confirmed. Then our theory is successfully applied to two typical systems including the geometric optics equations. Finally, we present the numerical results coinciding with the theoretical analysis.     

Analysis of stability and error estimates for three methods approximating a nonlinear reaction-diffusion equation

We present the error analysis of three time-stepping schemes used in the discretization of a nonlinear reaction-diffusion equation with Neumann boundary conditions, relevant in phase transition. We prove $L^\infty$ stability by maximum principle arguments, and derive error estimates using energy methods for the implicit Euler, and two implicit-explicit approaches, a linearized scheme and a fractional step method. A numerical experiment validates the theoretical results, comparing the accuracy of the methods.