A Godunov--type scheme for a scalar conservation law with space-time flux discontinuity

Kwame Atta Gyamfi

2024-02-06

Abstract:We present and analyze a new finite volume scheme of Gudonov-type for a nonlinear scalar conservation law whose flux function has a discontinuous coefficient due to time-dependent changes in its sign along a Lipschitz continuous curve.

Numerical Analysis

What problem does this paper attempt to address?

Convergence of a Finite Volume Scheme for Compactly Heterogeneous Scalar Conservation Laws

Abraham Sylla

2024-07-21

Abstract:We build a finite volume scheme for the scalar conservation law $\partial_t u + \partial_x (H(x, u)) = 0$ with bounded initial condition for a wide class of flux function $H$, convex with respect to the second variable. The main idea for the construction of the scheme is to use the theory of discontinuous flux. We prove that the resulting approximating sequence converges boundedly almost everywhere on $\mathopen]0, +\infty\mathclose[$ to the entropy solution.

Numerical Analysis
Numerical Study on Several Godunov-type Schemes with ALE Formulation

Baolin TIAN,Weidong SHEN,Yan LIU,Junbo CHENG,Shuanghu WANG

DOI: https://doi.org/10.3969/j.issn.1001-246X.2007.05.006

2007-01-01

Abstract:A Godunov-type ALE (Arbitrary Lagrangian-Eulerian) method based on finite volume technique is developed to simulate compressible multimaterial flows with large deformation. The method adapts well to Lagrangian, ALE or Eulerian formulations due to its arbitrary moving speed of meshes. Numerical features are studied in the formulations. Numerical traits of several Godunov-type schemes are studied. Resolutions of scheme for shock and contact are compared.
Discontinuous Galerkin Methods for Nonlinear Scalar Conservation Laws: Generalized Local Lax-Friedrichs Numerical Fluxes.

Jia Li,Dazhi Zhang,Xiong Meng,Boying Wu,Qiang Zhang

DOI: https://doi.org/10.1137/19m1243798

IF: 3.039

2020-01-01

SIAM Journal on Numerical Analysis

Abstract:In this paper, we study the discontinuous Galerkin (DG) method with a class of generalized numerical fluxes for one-dimensional scalar nonlinear conservation laws. The generalized local Lax-Friedrichs (GLLF) fluxes with two weights, which may not be monotone, are proposed and analyzed. Under a condition for the weights, we first show the monotonicity for the flux and thus the L-2 stability of the scheme. Then, by constructing and analyzing a special piecewise global projection which commutes with the time derivative operator, we are able to show optimal error estimates for the DG scheme with GLLF fluxes. The result is sharp for monotone numerical fluxes, for which only suboptimal estimates can be proved in previous work. Moreover, optimal error estimates are still valid for fluxes that are not monotone, allowing us to choose some suitable weights to achieve less numerical dissipation and thus to better resolve shocks. Numerical experiments are provided to show the sharpness of theoretical results.
Scalar conservation laws with stochastic discontinuous flux function

Lukas Brencher,Andrea Barth

DOI: https://doi.org/10.48550/arXiv.2107.00549

2021-07-01

Abstract:A variety of real-world applications are modeled via hyperbolic conservation laws. To account for uncertainties or insufficient measurements, random coefficients may be incorporated. These random fields may depend discontinuously on the state space, e.g., to represent permeability in a heterogeneous or fractured medium. We introduce a suitable admissibility criterion for the resulting stochastic discontinuous-flux conservation law and prove its well-posedness. Therefore, we ensure the pathwise existence and uniqueness of the corresponding deterministic setting and present a novel proof for the measurability of the solution, since classical approaches fail in the discontinuous-flux case. As an example of the developed theory, we present a specific advection coefficient, which is modeled as a sum of a continuous random field and a pure jump field. This random field is employed in the stochastic conservation law, in particular a stochastic Burgers' equation, for numerical experiments. We approximate the solution to this problem via the Finite Volume method and introduce a new meshing strategy that accounts for the resulting standing wave profiles caused by the flux-discontinuities. The ability of this new meshing method to reduce the sample-wise variance is demonstrated in numerous numerical investigations.

Numerical Analysis,Analysis of PDEs,Probability
Compactness estimates for difference schemes for conservation laws with discontinuous flux

Kenneth H. Karlsen,John D. Towers

2023-10-28

Abstract:We establish quantitative compactness estimates for finite difference schemes used to solve nonlinear conservation laws. These equations involve a flux function $f(k(x,t),u)$, where the coefficient $k(x,t$ is $BV$-regular and may exhibit discontinuities along curves in the $(x,t)$ plane. Our approach, which is technically elementary, relies on a discrete interaction estimate and the existence of one entropy function. While the details are specifically outlined for the Lax-Friedrichs scheme, the same framework can be applied to other difference schemes. Notably, our compactness estimates are new even in the homogeneous case ($k\equiv 1$).

Numerical Analysis
Conservation Laws with Discontinuous Gradient-Dependent Flux: the Unstable Case

Debora Amadori,Alberto Bressan,Wen Shen

2024-11-21

Abstract:The paper is concerned with a scalar conservation law with discontinuous gradient-dependent flux. Namely, the flux is described by two different functions $f(u)$ or $g(u)$, when the gradient $u_x$ of the solution is positive or negative, respectively. We study here the unstable case where $f(u)>g(u)$ for all $u\in {\mathbb R}$. Assuming that both $f$ and $g$ are strictly convex, solutions to the Riemann problem are constructed. Even for smooth initial data, examples show that infinitely many solutions can occur. For an initial data which is piecewise monotone, i.e., increasing or decreasing on a finite number of intervals, a global solution of the Cauchy problem is obtained. It is proved that such solution is unique under the additional requirement that the number of interfaces, where the flux switches between $f$ and $g$, remains as small as possible.

Analysis of PDEs
A Monotone Finite Volume Scheme for Diffusion Equations on General Non-Conforming Meshes

Qi Zhang,Zhiqiang Sheng,Guangwei Yuan

DOI: https://doi.org/10.1016/j.amc.2017.05.041

IF: 4.397

2017-01-01

Applied Mathematics and Computation

Abstract:A nonlinear monotone finite volume scheme on general non-conforming meshes for diffusion equations is introduced, which deals with discontinuous tensor coefficients rigorously. Since the expression of normal flux depends on auxiliary unknowns defined at cell-vertex including hanging nodes, we propose a new method to eliminate vertex-unknown by using primary unknowns at the centers of the cells sharing the vertex. Especially the unknowns defined on hanging nodes are eliminated by flux continuous conditions. The resulting scheme is monotone and preserves positivity of analytical solutions for strongly anisotropic and heterogeneous full tensor coefficient problems. Numerical results show that the convergent order of the monotone scheme by different methods of eliminating vertex unknowns will vary remarkably, and our new method can assure that it has almost second order accuracy and more accurate than some existing methods.
Superconvergence of Discontinuous Galerkin Methods for Scalar Nonlinear Conservation Laws in One Space Dimension.

Xiong Meng,Chi-Wang Shu,Qiang Zhang,Boying Wu

DOI: https://doi.org/10.1137/110857635

IF: 3.039

2012-01-01

SIAM Journal on Numerical Analysis

Abstract:In this paper, an analysis of the superconvergence property of the semidiscretediscontinuous Galerkin (DG) method applied to one-dimensionaltime-dependent nonlinear scalar conservation laws is carried out. Weprove that the error between the DG solution and a particularprojection of the exact solution achieves $\left(k + \frac32\right)$th ordersuperconvergence when upwind fluxes are used. The results hold truefor arbitrary nonuniform regular meshes and for piecewisepolynomials of degree $k$ ($k \ge 1$), under the condition that$|f'(u)|$ possesses a uniform positive lower bound. Numericalexperiments are provided to show that the superconvergence propertyactually holds true for nonlinear conservation laws with generalflux functions, indicating that the restriction on $f(u)$ isartificial.
Gas Evolution Dynamics in Godunov-Type Schemes

K. Xu

1998-01-01

Abstract:As a continuous eeort to understand the Godunov-type schemes, following the paper \Projection in this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers. More speciically, the underlying dynamics of Flux Vector Splitting (FVS) and Flux Diierence Splitting (FDS) schemes will be analyzed. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and numerical formulations, based on the governing equation of the discretized KFVS scheme, the weakness and advantages of FVS scheme are clearly observed. Also, in this paper, the implicit equilibrium assumption in the Godunov ux will be analyzed. Due to the numerical shock thickness related to the cell size, the numerical scheme should be able to capture both equilibrium and non-equilibrium ow behavior in smooth and discontinuous regions. The Godunov ux basically lacks the mechanism to capture nonequilibrium eeects in the artiicially enlarged numerical shock region and to stabilize the numerical shock structure. Consequently, the Godunov method is exposed to the possible spurious solutions, such as the carbuncle phenomena and odd-even decoupling, once the dissipation provided in the projection stage is not enough.
Conservation Laws with Discontinuous Gradient-Dependent Flux: the Stable Case

Debora Amadori,Alberto Bressan,Wen Shen

2024-11-16

Abstract:The paper is concerned with a scalar conservation law with discontinuous gradient-dependent flux. Namely, the flux is described by two different functions $f(u)$ or $g(u)$, when the gradient $u_x$ of the solution is positive or negative, respectively. We study here the stable case where $f(u)<g(u)$ for all $u\in {\mathbb R}$, with $f,g$ smooth but possibly not convex. A front tracking algorithm is introduced, proving that piecewise constant approximations converge to the trajectories of a contractive semigroup on $\mathbf{L}^1({\mathbb R})$. In the spatially periodic case, we prove that semigroup trajectories coincide with the unique limits of a suitable class of vanishing viscosity approximations.

Analysis of PDEs
Hyperbolic conservation laws on spacetimes. A finite volume scheme based on differential forms

Philippe G. LeFloch,Baver Okutmustur

DOI: https://doi.org/10.48550/arXiv.0810.0255

2008-10-02

Abstract:We consider nonlinear hyperbolic conservation laws, posed on a differential (n+1)-manifold with boundary referred to as a spacetime, and in which the "flux" is defined as a flux field of n-forms depending on a parameter (the unknown variable). We introduce a formulation of the initial and boundary value problem which is geometric in nature and is more natural than the vector field approach recently developed for Riemannian manifolds. Our main assumption on the manifold and the flux field is a global hyperbolicity condition, which provides a global time-orientation as is standard in Lorentzian geometry and general relativity. Assuming that the manifold admits a foliation by compact slices, we establish the existence of a semi-group of entropy solutions. Moreover, given any two hypersurfaces with one lying in the future of the other, we establish a "contraction" property which compares two entropy solutions, in a (geometrically natural) distance equivalent to the L1 distance. To carry out the proofs, we rely on a new version of the finite volume method, which only requires the knowledge of the given n-volume form structure on the (n+1)-manifold and involves the {\sl total flux} across faces of the elements of the triangulations, only, rather than the product of a numerical flux times the measure of that face.

Analysis of PDEs,Numerical Analysis
A personal discussion on conservation, and how to formulate it

Remi Abgrall

2023-04-18

Abstract:Since the celebrated theorem of Lax and Wendroff, we know a necessary condition that any numerical scheme for hyperbolic problem should satisfy: it should be written in flux form. A variant can also be formulated for the entropy. Even though some schemes, as for example those using continuous finite element, do not formally cast into this framework, it is a very convenient one. In this paper, we revisit this, introduce a different notion of local conservation which contains the previous one in one space dimension, and explore its consequences. This gives a more flexible framework that allows to get, systematically, entropy stable schemes, entropy dissipative ones, or accomodate more constraints. In particular, we can show that continuous finite element method can be rewritten in the finite volume framework, and all the quantities involved are explicitly computable. We end by presenting the only counter example we are aware of, i.e a scheme that seems not to be rewritten as a finite volume scheme.

Numerical Analysis
A Scalar Conservation Law for Plume Migration in Carbon Sequestration

Elisabeth Brown,Michael Shearer

DOI: https://doi.org/10.1137/17m1127089

IF: 1.9

2018-01-01

SIAM Journal on Applied Mathematics

Abstract:A quasi-linear hyperbolic partial differential equation with a discontinuous flux models geologic carbon dioxide (CO$_2$) migration and storage [M. Hesse, F. Orr, and H. Tchelepi, Fluid Mech., 611 (2008), pp. 35--60]. Two flux functions characterize the model, giving rise to flux discontinuities. One convex flux describes the invasion of the plume into pore space, and the other captures the flow as the plume leaves CO$_2$ bubbles behind, which are then trapped in the pore space. We investigate the method of characteristics, the structure of shock and rarefaction waves, and the result of binary wave interactions. The dual flux property introduces unexpected differences between the structure of these solutions and those of a scalar conservation law with a convex flux. During our analysis, we introduce a new construction of cross-hatch characteristics in regions of the space-time plane where the solution is constant, and there are two characteristic speeds. This construction is used to generalize the notion of the Lax entropy condition for admissible shocks and is crucial to continuing the propagation of a shock wave if its speed becomes characteristic.

mathematics, applied
Dissipative Mechanism in Godunov‐type Schemes

K Xu,ZW Li

DOI: https://doi.org/10.1002/fld.160

IF: 1.847

2001-01-01

International Journal for Numerical Methods in Fluids

Abstract:This paper concerns the dissipative mechanism in the shock capturing schemes, where exact or approximate Riemann solvers are used in the flux evaluation. More specifically, we are going to analyze the dissipation in the flux vector splitting (FVS) scheme and the Godunov method, from which some pathological phenomena from the FVS scheme and the Godunov method will be explained, such as the artificial dissipation and the shock instability. Copyright © 2001 John Wiley & Sons, Ltd.
A Discontinuous Galerkin Method for Stochastic Conservation Laws.

Yunzhang Li,Chi-Wang Sho,Shanjian Tang

DOI: https://doi.org/10.1137/19m125710x

IF: 2.968

2020-01-01

SIAM Journal on Scientific Computing

Abstract:In this paper we present a discontinuous Galerkin (DG) method to approximate stochastic conservation laws, which is an efficient high-order scheme. We study the stability for the semidiscrete DG methods for fully nonlinear stochastic equations. Error estimates are obtained for smooth solutions of semilinear stochastic equations with variable coefficients. We also establish a derivative-free second-order time discretization scheme for matrix-valued stochastic ordinary differential equations. Numerical experiments are performed to confirm the analytical results.
Analysis of the Local Discontinuous Galerkin Method with Generalized Fluxes for One-Dimensional Nonlinear Convection-Diffusion Systems

Hongjuan Zhang,Boying Wu,Xiong Meng

DOI: https://doi.org/10.1007/s11425-022-2035-y

2022-01-01

Science China Mathematics

Abstract:In this paper, we present optimal error estimates of the local discontinuous Galerkin method with generalized numerical fluxes for one-dimensional nonlinear convection-diffusion systems. The upwind-biased flux with the adjustable numerical viscosity for the convective term is chosen based on the local characteristic decomposition, which is helpful in resolving discontinuities of degenerate parabolic equations without enforcing any limiting procedure. For the diffusive term, a pair of generalized alternating fluxes are considered. By constructing and analyzing generalized Gauss-Radau projections with respect to different convective or diffusive terms, we derive optimal error estimates for nonlinear convection-diffusion systems with the symmetrizable flux Jacobian and fully nonlinear diffusive problems. Numerical experiments including long time simulations, different boundary conditions and degenerate equations with discontinuous initial data are provided to demonstrate the sharpness of theoretical results.
On Large Time Step Godunov Scheme for Hyperbolic Conservation Laws

Jinghua Wang,Hairui Wen,Tie Zhou

DOI: https://doi.org/10.4310/cms.2004.v2.n3.a7

2004-01-01

Communications in Mathematical Sciences

Abstract:In this paper we study the large time step (LTS) Godunov scheme for scalar hyperbolic conservation laws proposed by LeVeque.We show that for an arbitrary Courant number, all the possible wave interactions in each time step occur only in a finite number of cells, and the number of cells is bounded by a constant depending on the Courant number for a given initial value problem.As an application of the result mentioned above, we show that for any given Courant number, if the total variation of the initial value satisfies some conditions, then the numerical solutions of the LTS Godunov scheme will converge to the entropy solutions of the convex scalar conservation laws.
A positive meshless finite difference scheme for scalar conservation laws with adaptive artificial viscosity driven by fault detection

Cesare Bracco,Oleg Davydov,Carlotta Giannelli,Alessandra Sestini

2024-09-24

Abstract:We present a meshless finite difference method for multivariate scalar conservation laws that generates positive schemes satisfying a local maximum principle on irregular nodes and relies on artificial viscosity for shock capturing. Coupling two different numerical differentiation formulas and adaptive selection of the sets of influence allows to meet a local CFL condition without any a priori time step restriction. Artificial viscosity term is chosen in an adaptive way by applying it only in the vicinity of the sharp features of the solution identified by an algorithm for fault detection on scattered data. Numerical tests demonstrate a robust performance of the method on irregular nodes and advantages of adaptive artificial viscosity. The accuracy of the obtained solutions is comparable to that for standard monotone methods available only on Cartesian grids.

Numerical Analysis
Error analysis of an Algebraic Flux Correction Scheme for a nonlinear Scalar Conservation Law Using SSP-RK2

Christos Pervolianakis

2024-09-27

Abstract:We consider a scalar conservation law with linear and nonlinear flux function on a bounded domain $\Omega\subset{\R}^2$ with Lipschitz boundary $\partial\Omega.$ We discretize the spatial variable with the standard finite element method where we use a local extremum diminishing flux limiter which is linearity preserving. For temporal discretization, we use the second order explicit strong stability preserving Runge--Kutta method. It is known that the resulting fully-discrete scheme satisfies the discrete maximum principle. Under the sufficiently regularity of the weak solution and the CFL condition $k = \mathcal{O}(h^2)$, we derive error estimates in $L^{2}-$ norm for the algebraic flux correction scheme in space and in $\ell^\infty$ in time. We also present numerical experiments that validate that the fully-discrete scheme satisfies the temporal order of convergence of the fully-discrete scheme that we proved in the theoretical analysis.

Numerical Analysis

A Godunov--type scheme for a scalar conservation law with space-time flux discontinuity

Convergence of a Finite Volume Scheme for Compactly Heterogeneous Scalar Conservation Laws

Numerical Study on Several Godunov-type Schemes with ALE Formulation

Discontinuous Galerkin Methods for Nonlinear Scalar Conservation Laws: Generalized Local Lax-Friedrichs Numerical Fluxes.

Scalar conservation laws with stochastic discontinuous flux function

Compactness estimates for difference schemes for conservation laws with discontinuous flux

Conservation Laws with Discontinuous Gradient-Dependent Flux: the Unstable Case

A Monotone Finite Volume Scheme for Diffusion Equations on General Non-Conforming Meshes

Superconvergence of Discontinuous Galerkin Methods for Scalar Nonlinear Conservation Laws in One Space Dimension.

Gas Evolution Dynamics in Godunov-Type Schemes

Conservation Laws with Discontinuous Gradient-Dependent Flux: the Stable Case

Hyperbolic conservation laws on spacetimes. A finite volume scheme based on differential forms

A personal discussion on conservation, and how to formulate it

A Scalar Conservation Law for Plume Migration in Carbon Sequestration

Dissipative Mechanism in Godunov‐type Schemes

A Discontinuous Galerkin Method for Stochastic Conservation Laws.

Analysis of the Local Discontinuous Galerkin Method with Generalized Fluxes for One-Dimensional Nonlinear Convection-Diffusion Systems

On Large Time Step Godunov Scheme for Hyperbolic Conservation Laws

A positive meshless finite difference scheme for scalar conservation laws with adaptive artificial viscosity driven by fault detection

Error analysis of an Algebraic Flux Correction Scheme for a nonlinear Scalar Conservation Law Using SSP-RK2