Related papers: A fully well-balanced scheme for shallow water equ…
We develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravitational source term. Here, we advocate a new paradigm based on a purely conservative reformulation of the equations using global…
The Cahn-Hilliard-Navier-Stokes (CHNS) equations represent the fundamental building blocks of hydrodynamic phase-field models for multiphase fluid flow dynamics. Due to the coupling between the Navier-Stokes equation and the Cahn-Hilliard…
A stochastic Galerkin formulation for a stochastic system of balanced or conservation laws may fail to preserve hyperbolicity of the original system. In this work, we develop hyperbolicity-preserving stochastic Galerkin formulation for the…
We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of…
We present an energy/entropy stable and high order accurate finite difference (FD) method for solving the nonlinear (rotating) shallow water equations (SWEs) in vector invariant form using the newly developed dual-pairing and…
We present a central differencing scheme for the solution of the shallow water equations with non-flat bottom topography and moving wet-dry fronts. The problem is numerically challenging due to two reasons. First, the non-flat bottom…
The shear shallow water model is an extension of the classical shallow water model to include the effects of vertical shear. It is a system of six non-linear hyperbolic PDE with non-conservative products. We develop a high-order entropy…
The discontinuous Galerkin finite element method (DGFEM) developed by Rhebergen et al. (2008) offers a robust method for solving systems of nonconservative hyperbolic partial differential equations but, as we show here, does not…
We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly…
The shallow water equations provide a useful analogue of the fully compressible Euler equations since they have similar characteristics: conservation laws, inertia-gravity and Rossby waves and a (quasi-) balanced state. In order to obtain…
This paper aims to develop and analyze a numerical scheme for solving the backward problem of semilinear subdiffusion equations. We establish the existence, uniqueness, and conditional stability of the solution to the inverse problem by…
Because of their capability to preserve steady-states, well-balanced schemes for Shallow Water equations are becoming popular. Among them, the hydrostatic reconstruction proposed in Audusse et al. (2004), coupled with a positive numerical…
We present a high-order conservative, positivity-preserving, and non-oscillatory scheme for solving the Vlasov equation. The scheme attains formal fifth-order accuracy through a convex combination of positive and non-oscillatory polynomials…
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and…
The regularisation of nonlinear hyperbolic conservation laws has been a problem of great importance for achieving uniqueness of weak solutions and also for accurate numerical simulations. In a recent work, the first two authors proposed a…
The shallow water flow model is widely used to describe water flows in rivers, lakes, and coastal areas. Accounting for uncertainty in the corresponding transport-dominated nonlinear PDE models presents theoretical and numerical challenges…
This paper develops high-order well-balanced (WB) energy stable (ES) finite difference schemes for multi-layer (the number of layers $M\geqslant 2$) shallow water equations (SWEs) on both fixed and adaptive moving meshes, extending our…
It is well known, thanks to Lax-Wendroff theorem, that the local conservation of a numerical scheme for a conservative hyperbolic system is a simple and systematic way to guarantee that, if stable, a scheme will provide a sequence of…
Our goal was to develop a robust algorithm for numerical simulation of one-dimensional shallow-water flow in a complex multiply-connected channel network with arbitrary geometry and variable topography. We apply a central-upwind scheme with…
We present a new algorithm for solving ideal relativistic hydrodynamics based on Godunov method with an exact solution of Riemann problem for an arbitrary equation of state. Standard numerical tests are executed, such as the sound wave…