Well-Balanced High-Order Discontinuous Galerkin Methods for Systems of Balance Laws

This work introduces a general strategy to develop well-balanced high-order Discontinuous Galerkin (DG) numerical schemes for systems of balance laws. The essence of our approach is a local projection step that guarantees the exactly well-balanced character of the resulting numerical method for smooth stationary solutions. The strategy can be adapted to some well-known different time marching DG discretisations. Particularly, in this article, Runge–Kutta DG and ADER DG methods are studied. Additionally, a limiting procedure based on a modified WENO approach is described to deal with the spurious oscillations generated in the presence of non-smooth solutions, keeping the well-balanced properties of the scheme intact. The resulting numerical method is then exactly well-balanced and high-order in space and time for smooth solutions. Finally, some numerical results are depicted using different systems of balance laws to show the performance of the introduced numerical strategy.

A Solution of the Junction Riemann Problem for 1D Hyperbolic Balance Laws in Networks including Supersonic Flow Conditions on Elastic Collapsible Tubes

The numerical modeling of one-dimensional (1D) domains joined by symmetric or asymmetric bifurcations or arbitrary junctions is still a challenge in the context of hyperbolic balance laws with application to flow in pipes, open channels or blood vessels, among others. The formulation of the Junction Riemann Problem (JRP) under subsonic conditions in 1D flow is clearly defined and solved by current methods, but they fail when sonic or supersonic conditions appear. Formulations coupling the 1D model for the vessels or pipes with other container-like formulations for junctions have been presented, requiring extra information such as assumed bulk mechanical properties and geometrical properties or the extension to more dimensions. To the best of our knowledge, in this work, the JRP is solved for the first time allowing solutions for all types of transitions and for any number of vessels, without requiring the definition of any extra information. The resulting JRP solver is theoretically well-founded, robust and simple, and returns the evolving state for the conserved variables in all vessels, allowing the use of any numerical method in the resolution of the inner cells used for the space-discretization of the vessels. The methodology of the proposed solver is presented in detail. The JRP solver is directly applicable if energy losses at the junctions are defined. Straightforward extension to other 1D hyperbolic flows can be performed.

Existence and Uniqueness of Isothermal, Slightly Compressible Stratified Flow

We show well-posedness for the equations describing a new model of slightly compressible fluids. This model was recently rigorously derived in Grandi and Passerini (Geophys Astrophys Fluid Dyn, 2020) from the full set of balance laws and falls in the category of anelastic Navier–Stokes fluids. In particular, we prove existence and uniqueness of global regular solutions in the two-dimensional case for initial data of arbitrary “size”, and for “small” data in three dimensions. We also show global stability of the rest state in the class of weak solutions.

Collocation Methods for High-Order Well-Balanced Methods for Systems of Balance Laws

In some previous works, two of the authors introduced a technique to design high-order numerical methods for one-dimensional balance laws that preserve all their stationary solutions. The basis of these methods is a well-balanced reconstruction operator. Moreover, they introduced a procedure to modify any standard reconstruction operator, like MUSCL, ENO, CWENO, etc., in order to be well-balanced. This strategy involves a non-linear problem at every cell at every time step that consists in finding the stationary solution whose average is the given cell value. In a recent paper, a fully well-balanced method is presented where the non-linear problems to be solved in the reconstruction procedure are interpreted as control problems. The goal of this paper is to introduce a new technique to solve these local non-linear problems based on the application of the collocation RK methods. Special care is put to analyze the effects of computing the averages and the source terms using quadrature formulas. A general technique which allows us to deal with resonant problems is also introduced. To check the efficiency of the methods and their well-balance property, they have been applied to a number of tests, ranging from easy academic systems of balance laws consisting of Burgers equation with some non-linear source terms to the shallow water equations—without and with Manning friction—or Euler equations of gas dynamics with gravity effects.

Quadratic-stretch elasticity

A nonlinear small-strain elastic theory is constructed from a systematic expansion in Biot strains, truncated at quadratic order. The primary motivation is the desire for a clean separation between stretching and bending energies for shells, which appears to arise only from reduction of a bulk energy of this type. An approximation of isotropic invariants, bypassing the solution of a quartic equation or computation of tensor square roots, allows stretches, rotations, stresses, and balance laws to be written in terms of derivatives of position. Two-field formulations are also presented. Extensions to anisotropic theories are briefly discussed.

A finite element method for modeling surface growth and resorption of deformable solids

This work explores a continuum-mechanical model for a body simultaneously undergoing finite deformation and surface growth/resorption. This is accomplished by defining the kinematics as well as the set of material points that constitute the domain of a physical body at a given time in terms of an evolving reference configuration. The implications of spatial and temporal discretization are discussed, and an extension of the Arbitrary Lagrangian–Eulerian finite element method is proposed to enforce the resulting balance laws on the grown/resorbed body in two spatial dimensions. Representative numerical examples are presented to highlight the predictive capabilities of the model and the numerical properties of the proposed solution method.

The Weyl BMS group and Einstein’s equations

We propose an extension of the BMS group, which we refer to as Weyl BMS or BMSW for short, that includes super-translations, local Weyl rescalings and arbitrary diffeomorphisms of the 2d sphere metric. After generalizing the Barnich-Troessaert bracket, we show that the Noether charges of the BMSW group provide a centerless representation of the BMSW Lie algebra at every cross section of null infinity. This result is tantamount to proving that the flux-balance laws for the Noether charges imply the validity of the asymptotic Einstein’s equations at null infinity. The extension requires a holographic renormalization procedure, which we construct without any dependence on background fields. The renormalized phase space of null infinity reveals new pairs of conjugate variables. Finally, we show that BMSW group elements label the gravitational vacua.