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 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.     

High order positivity-preserving finite volume WENO schemes for a hierarchical size-structured population model

In this paper we develop high order positivity-preserving finite volume weighted essentially non-oscillatory (WENO) schemes for solving a hierarchical size-structured population model with nonlinear growth, mortality and reproduction rates. We carefully treat the technical complications in boundary conditions and global integration terms to ensure high order accuracy and the positivity-preserving property. Comparing with the previous high order difference WENO scheme for this model, the positivity-preserving finite volume WENO scheme has a comparable computational cost and accuracy, with the added advantages of being positivity-preserving and having L¹ stability. Numerical examples, including that of the evolution of the population of Gambusia affinis, are presented to illustrate the good performance of the scheme.     

Efficient implementation of WENO schemes to nonuniform meshes

Most of the standard papers about the WENO schemes consider their implementation to uniform meshes only. In that case the WENO reconstruction is performed efficiently by using the algebraic expressions for evaluating the reconstruction values and the smoothness indicators from cell averages. The coefficients appearing in these expressions are constant, dependent just on the scheme order, not on the mesh size or the reconstruction function values, and can be found, for example, in Jiang and Shu (J Comp Phys 126:202–228, 1996). In problems where the geometrical properties must be taken into account or the solution has localized fine scale structure that must be resolved, it is computationally efficient to do local grid refinement. Therefore, it is also desirable to have numerical schemes, which can be applied to nonuniform meshes. Finite volume WENO schemes extend naturally to nonuniform meshes although the reconstruction becomes quite complicated, depending on the complexity of the grid structure. In this paper we propose an efficient implementation of finite volume WENO schemes to nonuniform meshes. In order to save the computational cost in the nonuniform case, we suggest the way for precomputing the coefficients and linear weights for different orders of WENO schemes. Furthermore, for the smoothness indicators that are defined in an integral form we present the corresponding algebraic expressions in which the coefficients obtained as a linear combination of divided differences arise. In order to validate the new implementation, resulting schemes are applied in different test examples.      

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.     

High order Hybrid central—WENO finite difference scheme for conservation laws

In this article we present a high resolution hybrid central finite difference—WENO scheme for the solution of conservation laws, in particular, those related to shock–turbulence interaction problems. A sixth order central finite difference scheme is conjugated with a fifth order weighted essentially non-oscillatory WENO scheme in a grid-based adaptive way. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme while the smooth regions are computed with the more efficient and accurate central finite difference scheme. The application of high order filtering to mitigate the dispersion error of central finite difference schemes is also discussed. Numerical experiments with the 1D compressible Euler equations are shown.     

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.     

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.     

Balanced semi-implicit Roe scheme for open channel flow equations

In this paper, we present an extension of balanced upwind explicit finite volume schemes designed in Vuković and Sopta (SIAM J Sci Comput 24(5):1630–1649, 2003) to semi-implicit ones. Particularly, those schemes are applied to open-channel flows with general geometries and we verify the exact conservation property (C-property). We present the algorithm, the proof of exact C-property and results for several test cases. Also, we test and compare balanced semi-implicit and explicit schemes on standard test cases and on cases involving friction, non-uniform bed slopes and strong channel width variations.      

High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems

This paper is concerned with the development of high order methods for the numerical approximation of one-dimensional nonconservative hyperbolic systems. In particular, we are interested in high order extensions of the generalized Roe methods introduced by I. Toumi in 1992, based on WENO reconstruction of states. We also investigate the well-balanced properties of the resulting schemes. Finally, we will focus on applications to shallow-water systems.

     

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     

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.     

Shock reflection for a system of hyperbolic balance laws

This paper concerns shock reflection for a system of hyperbolic balance laws in one space dimension. It is shown that the generalized nonlinear initial-boundary Riemann problem for a system of hyperbolic balance laws with nonlinear boundary conditions in the half space admits a unique global piecewise C¹ solution u=u(t,x) containing only shocks with small amplitude and this solution possesses a global structure similar to that of self-similar solution of the corresponding homogeneous Riemann problem, if each characteristic field with positive velocity is genuinely nonlinear and the corresponding homogeneous Riemann problem has only shocks but no centered rarefaction waves and contact discontinuities. This result is also applied to shock reflection for the flow equations of a model class of fluids with viscosity induced by fading memory.     

Numerical smoothness and error analysis on WENO for the nonlinear conservation laws

In this study, we give an a posteriori error analysis on the weighted essentially nonoscillatory schemes for the nonlinear scalar conservation laws. This analysis is based on the new concept of numerical smoothness, with some new error analysis mechanisms developed for the finite difference and finite volume discretizations. The local error estimate is of optimal order in space and time. The global error estimate grows linearly in time, because of the direct application of the L₁ ‐contraction between entropy solutions in the error propagation analysis. As a beginning, we only deal with smooth solutions in this article. Within the same error propagation framework, when we deal with piecewise smooth solutions later, we only need to work on estimating the local error where smoothness is lost. The smoothness indicators not only serve the purpose of local error estimation, but also serve as a monitor on both the possible numerical instability and the expected solution shapening. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Some techniques for improving the resolution of finite difference component-wise WENO schemes for polydisperse sedimentation models

Polydisperse sedimentation models can be described by a system of conservation laws for the concentration of each species of solids. Some of these models, as the Masliyah–Locket–Bassoon model, can be proven to be hyperbolic, but its full characteristic structure cannot be computed in closed form. Component-wise finite difference WENO schemes may be used in these cases, but these schemes suffer from an excessive diffusion and may present spurious oscillations near shocks. In this work we propose to use a flux-splitting that prescribes less numerical viscosity for component-wise finite difference WENO schemes. We compare this technique with others to alleviate the diffusion and oscillatory behavior of the solutions obtained with component-wise finite difference WENO methods.     

Well-posedness theory for geometry-compatible hyperbolic conservation laws on manifolds

Motivated by many applications (geophysical flows, general relativity), we attempt to set the foundations for a study of entropy solutions to non-linear hyperbolic conservation laws posed on a (Riemannian or Lorentzian) manifold. The flux of the conservation laws is viewed as a vector-field on the manifold and depends on the unknown function as a parameter. We introduce notions of entropy solutions in the class of bounded measurable functions and in the class of measure-valued mappings. We establish the well-posedness theory for conservation laws on a manifold, by generalizing both Kruzkov's and DiPerna's theories originally developed in the Euclidian setting. The class of geometry-compatible (as we call it) conservation laws is singled out as an important case of interest, which leads to robust L^p estimates independent of the geometry of the manifold. On the other hand, general conservation laws solely enjoy the L¹ contraction property and leads to a unique contractive semi-group of entropy solutions. Our framework allows us to construct entropy solutions on a manifold via the vanishing diffusion method or the finite volume method.     

Shock waves for the hyperbolic balance laws with discontinuous sources

In this paper, we are concerned with the behavior of shock waves in a 2 × 2 balance law with discontinuous source terms. We obtain the existence of a local shock wave solution of this problem and deduce that the discontinuous source terms create a weak discontinuity in this solution. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

The Riemann problem for a class of resonant hyperbolic systems of balance laws

We solve the Riemann problem for a class of resonant hyperbolic systems of balance laws. The systems are not strictly hyperbolic and the solutions take their values in a neighborhood of a state where two characteristic speeds coincide. Our construction generalizes the ones given earlier by Isaacson and Temple for scalar equations and for conservative systems. The class of systems under consideration here includes, in particular, a model from continuum physics that describes the evolution of a fluid flow in a nozzle with discontinuous cross-section.

Résumé

Nous résolvons le problème de Riemann pour une classe de systèmes hyperboliques non-conservatifs et résonants. Ces systèmes ne sont pas strictement hyperboliques et les solutions considérées prennent leurs valeurs au voisinage d'un état constant où deux des vitesses caractéristiques coincident. Notre construction généralise celle donnée précédemment par Isaacson et Temple pour les équations scalaires et les systèmes conservatifs. La classe générale de systèmes étudiée ici comprend, en particulier, un modèle important de la dynamique des milieux continus qui décrit l'évolution d'un fluide dans une tuyère dont la section est discontinue.     

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.