Systems Of Nonlinear Conservation Laws Research Articles

An adaptive viscosity regularization approach for the numerical solution of conservation laws: Application to finite element methods

We introduce an adaptive viscosity regularization approach for the numerical solution of systems of nonlinear conservation laws with shock waves. The approach seeks to solve a sequence of regularized problems consisting of the system of conservation laws and an additional Helmholtz equation for the artificial viscosity. We propose a homotopy continuation of the regularization parameters to minimize the amount of artificial viscosity subject to positivity-preserving and smoothness constraints on the numerical solution. The regularization methodology is combined with a mesh adaptation strategy that identifies the shock location and generates shock-aligned meshes, which allows to further reduce the amount of artificial dissipation and capture shocks with increased accuracy. We use the hybridizable discontinuous Galerkin method to numerically solve the regularized system of conservation laws and the continuous Galerkin method to solve the Helmholtz equation for the artificial viscosity. We show that the approach can produce approximate solutions that converge to the exact solution of the Burgers' equation. Finally, we demonstrate the performance of the method on inviscid transonic, supersonic, hypersonic flows in two dimensions. The approach is found to be accurate, robust and efficient, and yields very sharp yet smooth solutions in a few homotopy iterations.

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.

Riemann problems for a hyperbolic system of nonlinear conservation laws from the Liou–Steffen pressure system

This paper is devoted to a hyperbolic system of nonlinear conservation laws, that is, the pressure system independent of density and energy from the Liou–Steffen flux-splitting scheme on the compressible Euler equations. First, the one-dimensional Riemann problem is solved with eight kinds of structures. Second, the two-dimensional Riemann problem is discussed; the solution reveals a variety of geometric structures; by the generalized characteristic analysis method and studying the pointwise interactions of waves, we construct 29 kinds of structures of solution consisting of shocks, rarefaction waves and contact discontinuities; the theoretical analysis is confirmed by numerical simulations.

A Well-Balanced Runge-Kutta Discontinuous Galerkin Method for Multilayer Shallow Water Equations with Non-Flat Bottom Topography

A well-balanced Runge-Kutta discontinuous Galerkin method is presented for the numerical solution of multilayer shallow water equations with mass exchange and non-flat bottom topography. The governing equations are reformulated as a nonlinear system of conservation laws with differential source forces and reaction terms. Coupling between the flow layers is accounted for in the system using a set of exchange relations. The considered well-balanced Runge-Kutta discontinuous Galerkin method is a locally conservative finite element method whose approximate solutions are discontinuous across the inter-element boundaries. The well-balanced property is achieved using a special discretization of source terms that depends on the nature of hydrostatic solutions along with the Gauss-LobattoLegendre nodes for the quadrature used in the approximation of source terms. The method can also be viewed as a high-order version of upwind finite volume solvers and it offers attractive features for the numerical solution of conservation laws for which standard finite element methods fail. To deal with the source terms we also implement a high-order splitting operator for the time integration. The accuracy of the proposed Runge-Kutta discontinuous Galerkin method is examined for several examples of multilayer free-surface flows over both flat and non-flat beds. The performance of the method is also demonstrated by comparing the results obtained using the proposed method to those obtained using the incompressible hydrostatic Navier-Stokes equations and a well-established kinetic method. The proposed method is also applied to solve a recirculation flow problem in the Strait of Gibraltar.

On a Two-Dimensional Riemann Problem for Hyperbolic System of Nonlinear Conservation Laws

On a Two-Dimensional Riemann Problem for Hyperbolic System of Nonlinear Conservation Laws

An order-adaptive compact approximation Taylor method for systems of conservation laws

We present a new family of high-order shock-capturing finite difference numerical methods for systems of conservation laws. These methods, called Adaptive Compact Approximation Taylor (ACAT) schemes, use centered (2p+1)-point stencils, where p may take values in {1,2,…,P} according to a new family of smoothness indicators in the stencils. The methods are based on a combination of a robust first order scheme and the Compact Approximate Taylor (CAT) methods of order 2p-order, p=1,2,…,P so that they are first order accurate near discontinuities and have order 2p in smooth regions, where (2p+1) is the size of the biggest stencil in which large gradients are not detected. CAT methods, introduced in [3], are an extension to nonlinear problems of the Lax-Wendroff methods in which the Cauchy-Kovalesky (CK) procedure is circumvented following the strategy introduced in [22] that allows one to compute time derivatives in a recursive way using high-order centered differentiation formulas combined with Taylor expansions in time. The expression of ACAT methods for 1D and 2D systems of conservative laws is given and the performance is checked in a number of test cases for several linear and nonlinear systems of conservation laws, including Euler equations for gas dynamics.

Lax-Wendroff Approximate Taylor Methods with Fast and Optimized Weighted Essentially Non-oscillatory Reconstructions

The goal of this work is to introduce new families of shock-capturing high-order numerical methods for systems of conservation laws that combine Fast WENO (FWENO) and Optimal WENO (OWENO) reconstructions with Approximate Taylor methods for the time discretization. FWENO reconstructions are based on smoothness indicators that require a lower number of calculations than the standard ones. OWENO reconstructions are based on a definition of the nonlinear weights that allows one to unconditionally attain the optimal order of accuracy regardless of the order of critical points. Approximate Taylor methods update the numerical solutions by using a Taylor expansion in time in which, instead of using the Cauchy–Kovalevskaya procedure, the time derivatives are computed by combining spatial and temporal numerical differentiation with Taylor expansions in a recursive way. These new methods are compared between them and against methods based on standard WENO implementations and/or SSP-RK time discretization. A number of test cases are considered ranging from scalar linear 1d problems to nonlinear systems of conservation laws in 2d.

On a new centered strategy to control the accuracy of weighted essentially non oscillatory algorithm for conservation laws close to discontinuities

AbstractIn this paper we try to redesign an algorithm analyzed in previous works by Amat et al. in order to obtain a centered strategy aimed to approximate the solution of linear and nonlinear systems of conservation laws. We perform the analysis of the new algorithm in order to explain why the new strategy works close to shocks in the solution of conservation laws. The objective is to force the use of the most centered stencil of those used whenever the data are affected by the presence of a high gradient.

A Source of Uncertainty in Computed Discontinuous Flows

It is speculated that some discontinuous weak solutions of boundary-value problems for nonlinear systems of conservation laws are computed, however routinely, with prescribed boundary data insufficient to uniquely determine such a solution. Stationary, transonic fluid flow exemplifies applications of present concern. A supplemental, a posteriori computation is described, which can potentially resolve this issue in any specific case.

Elementary waves, Riemann problem, Riemann invariants and new conservation laws for the pressure gradient model

We study the pressure gradient system of nonlinear conservation laws. We present the parametrization of shocks and rarefaction waves, then the solution of Riemann problem is given. Relying on the Riemann solution, we establish the uniqueness of the Riemann problem for the pressure gradient system. We also formulate the Riemann invariants corresponding to the pressure gradient system. We will present some interesting applications for the Riemann invariants. The first application is to represent the pressure gradient model in a diagonal representation. Actually this diagonal form admits the existence of global smooth solution for the pressure gradient system. The second application is to give new conservation laws corresponding to the pressure gradient model.

Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts

A class of wavelet-based methods, which combine adaptive grids together with second-order central and central-upwind high-resolution schemes, is introduced for accurate solution of first order hyperbolic systems of conservation laws and related equations. It is shown that the proposed class recovers stability which otherwise may be lost when the underlying central schemes are utilized over non-uniform grids. For simulations on irregular grids, both full-discretized and semi-discretized forms are derived; in particular, the effect of using certain shifts of the center of the computational allow to use standard slope limiters on non-uniform cells. Thereafter, the questions of numerical entropy production, local truncation errors, and the Total Variation Diminishing (TVD) criterion for scalar equations, are investigated. Central schemes are sensitives for irregular grids. To cure this drawback in the present context of multiresolution, the adapted grid is locally modified around high-gradient zones by adding new neighboring points. For an adapted point of resolution scale j, the new points are locally added in both resolution j and all successive coarser resolutions. It is shown that this simple grid modification improves stability of numerical solutions; this, however, comes at the expense that cell-centered central and central-upwind high resolution schemes, may not satisfy the TVD property. Finally, we present numerical simulations for both scalar and systems of non-linear conservation laws which confirm the simplicity and efficiency of the proposed method. These simulations demonstrate the high accuracy, the entropy production of wavelet-based algorithms which can effectively detect high gradient zones of shock waves, rarefaction regions, and contact discontinuities.

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case

This is the second part to our companion paper [18]. Herein, we generalize to two space dimensions the C-method developed in [20,18] for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts.We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm.To test the methodology, we use a highly simplified WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 1: The 1-D case

In this first part of two papers, we extend the C-method developed in [40] for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities in one space dimension. For gas dynamics, the C-method couples the Euler equations to a scalar reaction-diffusion equation, whose solution C serves as a space-time smooth artificial viscosity indicator.The purpose of this paper is the development of a high-order numerical algorithm for shock-wall collision and bounce-back. Specifically, we generalize the original C-method by adding a new collision indicator, which naturally activates during shock-wall collision. Additionally, we implement a new high-frequency wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We show that our scheme improves the order of accuracy of our WENO algorithm, handles extremely strong discontinuities (ranging up to nine orders of magnitude), allows for shock collision and bounce back, and removes high frequency noise. The causes of the well-known “wall heating” phenomenon are discussed, and we demonstrate that this particular pathology can be effectively treated in the framework of the C-method. This method is generalized to two space dimensions in the second part of this work [41].

High Order Edge Sensors with $\ell^1$ Regularization for Enhanced Discontinuous Galerkin Methods

This paper investigates the use of $\ell^1$ regularization for solving hyperbolic conservation laws based on high order discontinuous Galerkin (DG) approximations. We first use the polynomial annihilation method to construct a high order edge sensor which enables us to flag troubled elements. The DG approximation is enhanced in these troubled regions by activating $\ell^1$ regularization to promote sparsity in the corresponding jump function of the numerical solution. The resulting $\ell^1$ optimization problem is efficiently implemented using the alternating direction method of multipliers. By enacting $\ell^1$ regularization only in troubled cells, our method remains accurate and efficient, as no additional regularization or expensive iterative procedures are needed in smooth regions. We present results for the inviscid Burgers' equation as well as a nonlinear system of conservation laws using a nodal collocation-type DG method as a solver.

Continuum Dynamic Traffic Assignment Model for Autonomous Vehicles in a Polycentric Urban City with Environmental Consideration

This paper proposed a continuum dynamic model for autonomous vehicles in a polycentric urban city by considering the environment impact of traffic emission. The model assumes that homogeneous autonomous vehicles are continuously distributed over the urban areas which tend to choose a path to minimize their total travel cost from origin to destination. To describe the path choice behavior of travelers, we presented the continuum dynamic traffic assignment model which consists of a two-dimensional hyperbolic system of nonlinear conservation laws with source terms and an Eikonal-type equation. The elastic demand is considered using a function which associating each copy of flow with its total instantaneous travel cost. For the environmental impacts, here we consider the influence of CO emission and include the cost of emission into the actual transportation cost. A solution algorithm for the model is designed as a cell-centered finite volume method for conservation law equations and a fast sweeping method for Eikonal-type equations on unstructured grids. Numerical examples are given to demonstrate the model and the proposed solution algorithm. Further, the results of the travel cost considering CO emissions and not considering CO emissions are compared.

An Implicit Eulerian–Lagrangian WENO3 Scheme for Nonlinear Conservation Laws

We present a new, formally third order, implicit Weighted Essentially Non-Oscillatory (iWENO3) finite volume scheme for solving systems of nonlinear conservation laws. We then generalize it to define an implicit Eulerian–Lagrangian WENO (iEL-WENO) scheme. Implicitness comes from the use of an implicit Runge–Kutta (RK) time integrator. A specially chosen two-stage RK method allows us to drastically simplify the computation of the intermediate RK fluxes, leading to a computationally tractable scheme. The iEL-WENO3 scheme has two main steps. The first accounts for particles being transported within a grid element in a Lagrangian sense along the particle paths. Since this particle velocity is unknown (in a nonlinear problem), a fixed trace velocity v is used. The second step of the scheme accounts for the inaccuracy of the trace velocity v by computing the flux of particles crossing the incorrect tracelines. The CFL condition is relaxed when v is chosen to approximate the characteristic velocity. A new Roe solver for the Euler system is developed to account for the Lagrangian tracings, which could be useful even for explicit EL-WENO schemes. Numerical results show that iEL-WENO3 is both less numerically diffusive and can take on the order of about 2–3 times longer time steps than standard WENO3 for challenging nonlinear problems. An extension is made to the advection–diffusion equation. When advection dominates, the scheme retains its third order accuracy.

An Entropy Stable h / p Non-Conforming Discontinuous Galerkin Method with the Summation-by-Parts Property

This work presents an entropy stable discontinuous Galerkin (DG) spectral element approximation for systems of non-linear conservation laws with general geometric (h) and polynomial order (p) non-c ...

Divided difference estimates and accuracy enhancement of discontinuous Galerkin methods for nonlinear symmetric systems of hyperbolic conservation laws

In this paper, we investigate the accuracy-enhancement for the discontinuous Galerkin (DG) method for solving one-dimensional nonlinear symmetric systems of hyperbolic conservation laws. For nonlinear equations, the divided difference estimate is an important tool that allows for superconvergence of the post-processed solutions in the local L2-norm. Therefore, we first prove that the L2-norm of the α-th order (1≤ α≤ k+1) divided difference of the DG error with upwind fluxes is of order k+(3-α)/2, provided that the flux Jacobian matrix, f'(u), is symmetric positive definite. Furthermore, using the duality argument, we are able to derive superconvergence estimates of order 2k+(3-α)/2 for the negative-order norm, indicating that some particular compact kernels can be used to extract at least (3k/2+1)-th order superconvergence for nonlinear systems of conservation laws. Numerical experiments are shown to demonstrate the theoretical results.

An overview on the approximation of boundary Riemann problems through physical viscosity

This note aims at providing an overview of some recent results concerning the viscous approximation of so-called boundary Riemann problems for nonlinear systems of conservation laws in small total variation regimes.

Crowd dynamics through non-local conservation laws

We present a Lax-Friedrichs type scheme to compute the solutions of a class of non-local and non-linear systems of conservation laws in several space dimensions. The convergence of the approximate solutions is proved by providing suitable L 1 , L ∞ and BV uniform bounds. To illustrate the performances of the scheme, we consider an application to crowd dynamics. Numerical integrations show the formation of lanes in groups moving in opposite directions. This is joint work with R.M. Colombo (INDAM Unit, University of Brescia).