Well-balanced Scheme Research Articles

A New Approach for Designing Well-Balanced Schemes for the Shallow Water Equations: A Combination of Conservative and Primitive Formulations

A New Approach for Designing Well-Balanced Schemes for the Shallow Water Equations: A Combination of Conservative and Primitive Formulations

A Godunov-type scheme for the Saint-Venant-Exner equations with a moving steady states

The Saint-Venant-Exner system is widely used in industrial codes to model the transport of bed sediments. In practice, most of the methods used for its numerical resolution suffer from significant stability problems due to the fact that they are mainly guaranteed by separation techniques which allow weak coupling between hydraulic and morphodynamic software. The search for an adequate alternative method to this problem is a focus towards which several approaches converge. Recently, many authors have proposed acceptable methods but they are not so trivial to implement in an industrial context and some do not take into account all aspects such as mobile steady states. In this work, we derive a well-balanced scheme that takes into account all stable states, to approach weak solutions of the sediment transport model in shallow waters. We demonstrate that it captures exactly all regular steady solutions with non-evanescent velocities, retains water level positivity and is able to degenerate to a classical shallow water model when the sedimentary transport flow is null. A number of numerical test cases are also presented to illustrate these properties.

An explicit well-balanced scheme on staggered grids for barotropic Euler equations

In this paper, we introduce a specific modification of the numerical fluxes in order to insure the well-balanced property of schemes on staggered grids for the Euler equations. This property is crucial for the numerical representation of equilibrium solutions of balance laws with source terms, like when describing flows subjected to gravity and a complex topography. We propose first and second order versions of the well-balanced scheme. The performances of the method are evaluated through a series of 1D and 2D benchmarks.

A fully well-balanced hydrodynamic reconstruction

Abstract The present work focuses on the numerical approximation of the weak solutions of the shallow water model over a non-flat topography. In particular, we pay close attention to steady solutions with nonzero velocity. The goal of this work is to derive a scheme that exactly preserves these stationary solutions, as well as the commonly preserved lake at rest steady solution. These moving steady states are solution to a nonlinear equation. We emphasize that the method proposed here never requires solving this nonlinear equation; instead, a suitable linearization is derived. To address this issue, we propose an extension of the well-known hydrostatic reconstruction. By appropriately defining the reconstructed states at the interfaces, any numerical flux function, combined with a relevant source term discretization, produces a well-balanced scheme that preserves both moving and non-moving steady solutions. This eliminates the need to construct specific numerical fluxes. Additionally, we prove that the resulting scheme is consistent with the homogeneous system on flat topographies, and that it reduces to the hydrostatic reconstruction when the velocity vanishes. To increase the accuracy of the simulations, we propose a well-balanced high-order procedure, which still does not require solving any nonlinear equation. Several numerical experiments demonstrate the effectiveness of the numerical scheme.

Improved well-balanced AWENO schemes with hydrostatic reconstruction for the Euler equations under gravitational fields

The Euler equations under gravitational fields allow the hydrostatic equilibrium states, which requires that the numerical scheme of the system should also have this characteristic. In our previous work, a well-balanced finite difference conservative AWENO scheme has been constructed to preserve the isothermal equilibrium state accurately (Fu et al., Appl. Numer. Math., 180:1-15, 2022). However, the scheme fails to maintain more complicated hydrostatic equilibrium states, such as the isentropic equilibrium state. In this study, the improved well-balanced AWENO schemes with hydrostatic reconstruction are designed to maintain the general equilibrium states numerically. The key technologies to achieve the well-balanced property are rewriting the source term into an appropriate form based on the equilibrium equation, acting on the conserved variables with an affine-invariant WENO or linearized WENO interpolation operator, employing hydrostatic reconstruction methods to create auxiliary conserved variables, and modifying the numerical flux with auxiliary conserved variables and higher-order conserved terms based on pressure. The proposed well-balanced schemes maintain the equilibrium solution near the machine precision and capture small perturbations of the equilibrium state without spurious oscillations even at the coarse mesh. Furthermore, they can be applicable to the general hydrostatic equilibrium more than the isothermal equilibrium state. Several numerical experiments are tested to show that the developed well-balanced schemes perform well in maintaining equilibrium states and capturing small perturbations.

A Lifshitz-Slyozov type model for adipocyte size dynamics: limit from Becker-Döring system and numerical simulation.

Biological data show that the size distribution of adipose cells follows a bimodal distribution. In this work, we introduce a Lifshitz-Slyozov type model, based on a transport partial differential equation, for the dynamics of the size distribution of adipose cells. We prove a new convergence result from the related Becker-Döring model, a system composed of several ordinary differential equations, toward mild solutions of the Lifshitz-Slyozov model using distribution tail techniques. Then, this result allows us to propose a new advective-diffusive model, the second-order diffusive Lifshitz-Slyozov model, which is expected to better fit the experimental data. Numerical simulations of the solutions to the diffusive Lifshitz-Slyozov model are performed using a well-balanced scheme and compared to solutions to the transport model. Those simulations show that both bimodal and unimodal profiles can be reached asymptotically depending on several parameters. We put in evidence that the asymptotic profile for the second-order system does not depend on initial conditions, unlike for the transport Lifshitz-Slyozov model.

High-order well-balanced numerical schemes for one-dimensional shallow-water systems with Coriolis terms

The goal of this work is to develop high-order well-balanced schemes for the one-dimensional shallow-water equations with Coriolis terms. The main contribution is the development of general numerical methods that allow the achievement of arbitrary high-order for the shallow-water equations with Coriolis terms while preserving all the stationary solutions.

A scalable well-balanced numerical scheme for the simulation of fast landslides with efficient time stepping

We consider a single-phase depth-averaged model for the numerical simulation of fast-moving landslides with the goal of constructing a well-balanced, yet scalable and efficient, second-order time-stepping algorithm. We apply a Strang splitting approach to distinguish between parabolic and hyperbolic problems. For the parabolic contribution, we adopt a second-order Implicit-Explicit Runge-Kutta-Chebyshev scheme, while we use a two-stage Taylor discretization combined with a path-conservative strategy, to deal with the purely hyperbolic contribution. The proposed strategy allows to decouple hyperbolic from parabolic-reaction stiff contributions resulting in an overall well-balanced scheme subject just to stability restrictions of the hyperbolic term. The spatial discretization we adopt is based on a standard finite element method, associated with a hierarchically refined Cartesian grid. After providing numerical evidence of the well-balancing property, we demonstrate the capability of the proposed approach to select time steps larger than the ones adopted by a classical Taylor-Galerkin scheme. Finally, we provide some meaningful scaling results on ideal and realistic scenarios.

A well-balanced scheme using exact solutions to the two species Vlasov-Poisson system

In this work, we consider the numerical approximation of the two-species Vlasov-Poisson system using Eulerian methods. A family of exact non homogeneous stationary solutions are constructed using elliptic functions. Then, specific numerical schemes are proposed to compute solutions which remain close to a given stationary solution since standard schemes fail to capture such a dynamics on coarse meshes. The strategy is based on a suitable decomposition of the solution (in the spirit of δ f approaches) which can be easily combined with any classical Vlasov solvers. For unstable dynamics, a projection technique is proposed in order to dynamically change the equilibrium. Finally, numerical tests are proposed to illustrate the good behavior of the proposed strategy.

High Order Asymptotic Preserving and Well-Balanced Schemes for the Shallow Water Equations with Source Terms

High Order Asymptotic Preserving and Well-Balanced Schemes for the Shallow Water Equations with Source Terms

High-order compact gas-kinetic scheme for two-layer shallow water equations on unstructured mesh

For the two-layer shallow water equations, a high-order compact gas-kinetic scheme (GKS) on triangular mesh is proposed. The two-layer shallow water equations have complex source terms in comparison with the single layer equations. The main focus of this study is to construct a time-accurate evolution solution at a cell interface and to design a well-balanced scheme. The evolution model at a cell interface provides not only the numerical fluxes, but also the flow variables. In the development of the well-balanced scheme, the time-dependent flow variables at the cell interfaces surrounding the cell can be utilized to update the cell-averaged gradients. These time-evolving gradients are subsequently applied in the discretization of source terms inside each control volume. Based on the cell-averaged flow variables and their gradients, high-order initial data reconstruction can be achieved with compact stencils. The compact high-order GKS has advantages to simulate the flow evolution in complex domain covered by unstructured mesh. Many test cases are used to validate the accuracy and robustness of the scheme for the two-layer shallow water equations.

An operator-splitting approach with a hybrid finite volume/finite difference scheme for extended Boussinesq models

In this paper, we present a hybrid numerical scheme combining finite volume and finite difference methods (FV/FD), for solving two extended Boussinesq-type (BT) models. We propose two new reformulations specially designed to allow a separation between the hyperbolic part and dispersive part of the equations. Taking advantage of the proposed reformulations, we employ a splitting scheme specifically tailored to comprehensively capture the broad spectrum of natural events that are observed in the study of coastal oceanography. We employ a high-order FV well-balanced scheme to effectively handle the hyperbolic portion of the equations, guaranteeing non-negativity of water depth and accommodating dry regions. We solve the remaining dispersive component by employing a traditional FD approach. Several test cases in one horizontal dimension are validated showing that the proposed approach is able to accurately replicate coastal wave dynamics. In comparison to the non-linear shallow water equations (NSWE), the BT models exhibit significantly enhanced accuracy in simulating highly dispersive waves across varying water depths.

A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume / finite element scheme for the incompressible MHD equations

We present a new exactly divergence-free and well-balanced hybrid finite volume/finite element scheme for the numerical solution of the incompressible viscous and resistive magnetohydrodynamics (MHD) equations on staggered unstructured mixed-element meshes in two and three space dimensions. The equations are split into several subsystems, each of which is then discretized with a particular scheme that allows to preserve some fundamental structural features of the underlying governing PDE system also at the discrete level. The pressure is defined on the vertices of the primary mesh, while the velocity field and the normal components of the magnetic field are defined on an edge-based/face-based dual mesh in two and three space dimensions, respectively. This allows to account for the divergence-free conditions of the velocity field and of the magnetic field in a rather natural manner. The non-linear convective and the viscous terms in the momentum equation are solved at the aid of an explicit finite volume scheme, while the magnetic field is evolved in an exactly divergence-free manner via an explicit finite volume method based on a discrete form of the Stokes law in the edges/faces of each primary element. The latter method is stabilized by the proper choice of the numerical resistivity in the computation of the electric field in the vertices/edges of the 2D/3D elements. To achieve higher order of accuracy, a piecewise linear polynomial is reconstructed for the magnetic field, which is guaranteed to be exactly divergence-free via a constrained L2 projection. Finally, the pressure subsystem is solved implicitly at the aid of a classical continuous finite element method in the vertices of the primary mesh and making use of the staggered arrangement of the velocity, which is typical for incompressible Navier-Stokes solvers. In order to maintain non-trivial stationary equilibrium solutions of the governing PDE system exactly, which are assumed to be known a priori, each step of the new algorithm takes the known equilibrium solution explicitly into account so that the method becomes exactly well-balanced. We show numerous test cases in two and three space dimensions in order to validate our new method carefully against known exact and numerical reference solutions. In particular, this paper includes a very thorough study of the lid-driven MHD cavity problem in the presence of different magnetic fields and the obtained numerical solutions are provided as free supplementary electronic material to allow other research groups to reproduce our results and to compare with our data. We finally present long-time simulations of Soloviev equilibrium solutions in several simplified 3D tokamak configurations, showing that the new well-balanced scheme introduced in this paper is able to maintain stationary equilibria exactly over very long integration times even on very coarse unstructured meshes that, in general, do not need to be aligned with the magnetic field.

On high order positivity-preserving well-balanced finite volume methods for the Euler equations with gravitation

This paper studies three high-order structure-preserving finite volume weighted essentially non-oscillatory (WENO) methods, which are not only well balanced (WB) for a general known hydrostatic equilibrium state but also preserve the positivity of density and pressure, for the compressible Euler equations under gravitational fields. These methods are built on a simple local scaling positivity-preserving (PP) limiter and a modified WENO-ZQ reconstruction exactly preserving the cell average value and scaling invariance. The WB properties of these three methods are achieved based on suitable numerical fluxes and approximation to the gravitational source terms. Based on some convex decomposition techniques as well as several critical properties of the admissible states and numerical flux, we carry out rigorous positivity-preserving analyses for these three WB schemes. We rigorously prove that the three WB methods, coupled with the PP limiter and a strong-stability-preserving time discretization, are always PP under suitable Courant-Friedrichs-Lewy conditions. Extensive numerical examples are provided to confirm WB and PP properties of three methods.

A Well-Balanced Scheme for Euler Equations with Singular Sources

A Well-Balanced Scheme for Euler Equations with Singular Sources

Numerical modelling of dense snow avalanches with a well-balanced scheme based on the 2D shallow water equations

Abstract A common technique for simulating non–Newtonian fluid dynamics, such as snow avalanches, is to solve the Shallow Water Equations (SWE), together with a rheological model describing the momentum dissipation by shear stresses. Friction and cohesion terms are commonly modelled using the Voellmy friction model and, recently, the Bartelt cohesion model. Here, an adaptation of the Roe scheme that ensures the balance between the flux and pressure gradients and the friction source term is presented. An upwind scheme was used for the discretisation of the SWE numerical fluxes and the non–velocity-dependent terms of the friction–cohesion model, whereas a centred scheme was used for the velocity-dependent source terms. The model was tested in analytically solvable settings, laboratory experiments and real cases. In all cases, the model performed well, avoiding numerical instabilities and achieving stable and consistent solution even for an avalanche stopping on a sloping terrain.

Numerical Simulation of Sediment Transport in Unsteady Open Channel Flow

This paper presented a two-dimensional, well-balanced hydrodynamic and sediment transport model based on the solutions of variable-density shallow-water equations (VDSWEs) and the Exner equation for bed change for simulating sediment transport in unsteady flows. Those equations are solved in a coupled way by the first-order Godunov-type finite volume method. The Harten–Lax–van Leer–Contact (HLLC) Riemann solver is extended to find the local Riemann fluxes to maintain the exact balance between the momentum term and the bed slope term. A well-balanced scheme is superior to an unbalanced scheme to minimize numerical dispersion as demonstrated by the synthetic standing contact-discontinuity test case. Following this, the model is employed to simulate two laboratory experiments and a field case, the 1996 Lake Ha! Ha! flood event in Canada. The results of water surface elevations and bed surface profiles agree well with the measurements. The accuracy and robustness of the numerical schemes make the model a good candidate for practical engineering applications.

A family of well-balanced WENO and TENO schemes for atmospheric flows

We herein present a novel methodology to construct very high order well-balanced schemes for the computation of the Euler equations with gravitational source term, with application to numerical weather prediction (NWP). The proposed method is based on augmented Riemann solvers, which allow preserving the exact equilibrium between fluxes and source terms at cell interfaces. In particular, the augmented HLL solver (HLLS) is considered. Different spatial reconstruction methods can be used to ensure a high order of accuracy in space (e.g. WENO, TENO, linear reconstruction), being the TENO reconstruction the preferred method in this work. To the knowledge of the authors, the TENO method has not been applied to NWP before, although it has been extensively used by the computational fluid dynamics community in recent years. Therefore, we offer a thorough assessment of the TENO method to evidence its suitability for NWP considering some benchmark cases which involve inertia and gravity waves as well as convective processes. The TENO method offers an enhanced behavior when dealing with turbulent flows and underresolved solutions, where the traditional WENO scheme proves to be more diffusive. The proposed methodology, based on the HLLS solver in combination with a very high-order discretization, allows carrying out the simulation of meso- and micro-scale atmospheric flows in an implicit Large Eddy Simulation manner. Due to the HLLS solver, the isothermal, adiabatic and constant Brunt-Väisälä frequency hydrostatic equilibrium states are preserved with machine accuracy.

Aiolos – A multipurpose 1D hydrodynamics code for planetary atmospheres

ABSTRACT We present a new 1D multiphysics simulation code with use cases intended for, but not limited to, hydrodynamic escape problems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species as separated hydrodynamic fields, couples them via friction laws, allows for a multiband flux-limited radiation transport, and tracks ionization fronts in high-energy irradiation bands. Besides coupling various known numerical solution techniques together, we improve on the numerical stability of deep hydrostatic atmospheres by using a well-balanced scheme, hence preventing unphysical driving of atmospheric in- or outflow. We demonstrate the correct physical behaviour of the individual code modules and present a few simple, new applications, such as a proof-of-concept simulations of combined core-powered mass-loss and UV-driven atmospheric escape, along with a fully time-dependent core-collapse giant planet simulation. The multispecies nature of the code opens up the area of exploring simulations that are agnostic towards the dominant atmospheric species and can lead to implementations of advanced planetary evolution schemes.

A Well-Balanced FVC Scheme for 2D Shallow Water Flows on Unstructured Triangular Meshes

We consider in this work the numerical resolution of a 2D shallow water system with a Coriolis effect and bottom friction stresses on unstructured meshes by a new Finite Volume Characteristics (FVC) scheme, which has been introduced in the preliminary works that will be cited below. Our main goal is to extend this approach to 2D unstructured formalism while preserving the physical and mathematical properties of the system, including the C-property. First, we present our extension by preserving the advantages of the finite volume discretization such as conservation property and the method of characteristics such as elimination of Riemann solvers. Afterward, an approach was applied to the topography source term that leads to a well-balanced scheme satisfying the steady-state condition of still water. A semi-implicit treatment will also be presented in this study to avoid stability problems for the other source terms. Finally, the proposed finite volume method is verified on several benchmark tests and shows good agreement with analytical solutions and experimental results; moreover, it gives a noticeable accuracy and rapidity improvement compared to the original approaches.