No Arabic abstract
We consider entropy solutions to the initial value problem associated with scalar nonlinear hyperbolic conservation laws posed on the two-dimensional sphere. We propose a finite volume scheme which relies on a web-like mesh made of segments of longitude and latitude lines. The structure of the mesh allows for a discrete version of a natural geometric compatibility condition, which arose earlier in the well-posedness theory established by Ben-Artzi and LeFloch. We study here several classes of flux vectors which define the conservation law under consideration. They are based on prescribing a suitable vector field in the Euclidean three-dimensional space and then suitably projecting it on the spheres tangent plane; even when the flux vector in the ambient space is constant, the corresponding flux vector is a non-trivial vector field on the sphere. In particular, we construct here equatorial periodic solutions, analogous to one-dimensional periodic solutions to one-dimensional conservation laws, as well as a wide variety of stationary (steady state) solutions. We also construct confined solutions, which are time-dependent solutions supported in an arbitrarily specified subdomain of the sphere. Finally, representative numerical examples and test-cases are presented.
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.
In this paper, we propose a hybrid finite volume Hermite weighted essentially non-oscillatory (HWENO) scheme for solving one and two dimensional hyperbolic conservation laws. The zeroth-order and the first-order moments are used in the spatial reconstruction, with total variation diminishing Runge-Kutta time discretization. The main idea of the hybrid HWENO scheme is that we first use a shock-detection technique to identify the troubled cell, then, if the cell is identified as a troubled cell, we would modify the first order moment in the troubled cell and employ HWENO reconstruction in spatial discretization; otherwise, we directly use high order linear reconstruction. Unlike other HWENO schemes, we borrow the thought of limiter for discontinuous Galerkin (DG) method to control the spurious oscillations, after this procedure, the scheme would avoid the oscillations by using HWENO reconstruction nearby discontinuities and have higher efficiency for using linear approximation straightforwardly in the smooth regions. In addition, the hybrid HWENO scheme still keeps the compactness. A collection of benchmark numerical tests for one and two dimensional cases are performed to demonstrate the numerical accuracy, high resolution and robustness of the proposed scheme.
In this paper, a fifth-order Hermite weighted essentially non-oscillatory (HWENO) scheme with artificial linear weights is proposed for one and two dimensional hyperbolic conservation laws, where the zeroth-order and the first-order moments are used in the spatial reconstruction. We construct the HWENO methodology using a nonlinear convex combination of a high degree polynomial with several low degree polynomials, and the associated linear weights can be any artificial positive numbers with only requirement that their summation equals one. The one advantage of the HWENO scheme is its simplicity and easy extension to multi-dimension in engineering applications for we can use any artificial linear weights which are independent on geometry of mesh. The another advantage is its higher order numerical accuracy using less candidate stencils for two dimensional problems. In addition, the HWENO scheme still keeps the compactness as only immediate neighbor information is needed in the reconstruction and has high efficiency for directly using linear approximation in the smooth regions. In order to avoid nonphysical oscillations nearby strong shocks or contact discontinuities, we adopt the thought of limiter for discontinuous Galerkin method to control the spurious oscillations. Some benchmark numerical tests are performed to demonstrate the capability of the proposed scheme.
The paper proposes a scheme by combining the Runge-Kutta discontinuous Galerkin method with a {delta}-mapping algorithm for solving hyperbolic conservation laws with discontinuous fluxes. This hybrid scheme is particularly applied to nonlinear elasticity in heterogeneous media and multi-class traffic flow with inhomogeneous road conditions. Numerical examples indicate the schemes efficiency in resolving complex waves of the two systems. Moreover, the discussion implies that the so-called {delta}-mapping algorithm can also be combined with any other classical methods for solving similar problems in general.
This paper addresses the three concepts of textit{ consistency, stability and convergence } in the context of compact finite volume schemes for systems of nonlinear hyperbolic conservation laws. The treatment utilizes the framework of balance laws. Such laws express the relevant physical conservation laws in the presence of discontinuities. Finite volume approximations employ this viewpoint, and the present paper can be regarded as being in this category. It is first shown that under very mild conditions a weak solution is indeed a solution to the balance law. The schemes considered here allow the computation of several quantities per mesh cell (e.g., slopes) and the notion of consistency must be extended to this framework. Then a suitable convergence theorem is established, generalizing the classical convergence theorem of Lax and Wendroff. Finally, the limit functions are shown to be entropy solutions by using a notion of Godunov compatibility, which serves as a substitute to the entropy condition.