Matania Ben-Artzi

Matania Ben-Artzi, Joseph Falcovitz, Philippe G. LeFloch

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 sphere's 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.

Matania Ben-Artzi, Guy Katriel

The biharmonic operator plays a central role in a wide array of physical models, notably in elasticity theory and the streamfunction formulation of the Navier-Stokes equations. The need for corresponding numerical simulations has led, in recent years, to the development of a discrete biharmonic calculus. The primary object of this calculus is a high-order compact discrete biharmonic operator (DBO). The numerical results have been remarkably accurate, and have been corroborated by some rigorous proofs. However, there remained the "mystery" of the "underlying reason" for this success. This paper is a contribution in this direction, expounding the strong connection between cubic spline functions (on an interval) and the DBO. It is shown in particular that the (scaled) fourth-order distributional derivative of the cubic spline is identical to the action of the DBO on grid functions. The DBO is constructed in terms of the discrete Hermitian derivative. A remarkable fact is that the kernel of the inverse of the discrete operator is (up to scaling) equal to the grid evaluation of the kernel of $\Big[\Big(\frac{d}{dx}\Big)^4\Big]^{-1} .$ Explicit expressions are presented for both kernels. The relation between the (infinite) set of eigenvalues of the fourth-order Sturm-Liouville problem and the finite set of eigenvalues of the discrete biharmonic operator is studied, and the discrete eigenvalues are proved to converge (at an "optimal" $O(h^4)$ rate) to the continuous ones. Another remarkable consequence is the validity of a comparison principle. It is well known that there is no maximum principle for the fourth-order equation. However, a positivity result is derived, both for the continuous and the discrete biharmonic equation, showing that in both cases the kernels are order preserving.

Matania Ben-Artzi, Philippe G. LeFloch

Motivated by many applications (geophysical flows, general relativity), we attempt to set the foundations for a study of entropy solutions to nonlinear 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 {\sl 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^1$ 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.

Paulo Amorim, Matania Ben-Artzi, Philippe G. LeFloch

This paper investigates some properties of entropy solutions of hyperbolic conservation laws on a Riemannian manifold. First, we generalize the Total Variation Diminishing (TVD) property to manifolds, by deriving conditions on the flux of the conservation law and a given vector field ensuring that the total variation of the solution along the integral curves of the vector field is non-increasing in time. Our results are next specialized to the important case of a flow on the 2-sphere, and examples of flux are discussed. Second, we establish the convergence of the finite volume methods based on numerical flux-functions satisfying monotonicity properties. Our proof requires detailed estimates on the entropy dissipation, and extends to general manifolds an earlier proof by Cockburn, Coquel, and LeFloch in the Euclidian case.