Related papers: Structure preserving schemes for a class of Wasser…
In this article, we formulate topology optimization problems concerning the mass distribution as minimization problems for functionals on the Wasserstein space. We relax optimization problems regarding non-convex objective functions on the…
In this work we consider an extension of a recently proposed structure preserving numerical scheme for nonlinear Fokker-Planck-type equations to the case of nonconstant full diffusion matrices. While in existing works the schemes are…
We study a Lagrangian numerical scheme for solution of a nonlinear drift diffusion equation on an interval. The discretization is based on the equation's gradient flow structure with respect to the Wasserstein distance. The scheme inherits…
We propose and analyze two novel decoupled numerical schemes for solving the Cahn-Hilliard-Stokes-Darcy (CHSD) model for two-phase flows in karstic geometry. In the first numerical scheme, we explore a fractional step method (operator…
We propose a novel class of temporal high-order parametric finite element methods for solving a wide range of geometric flows of curves and surfaces. By incorporating the backward differentiation formulae (BDF) for time discretization into…
In this paper we present a numerical scheme for nonlinear continuity equations, which is based on the gradient flow formulation of an energy functional with respect to the quadratic transportation distance. It can be applied to a large…
We develop a set of numerical schemes for the Poisson--Nernst--Planck equations. We prove that our schemes are mass conservative, uniquely solvable and keep positivity unconditionally. Furthermore, the first-order scheme is proven to be…
This paper studies two hybrid discontinuous Galerkin (HDG) discretizations for the velocity-density formulation of the compressible Stokes equations with respect to several desired structural properties, namely provable convergence, the…
Structure-preserving finite-difference schemes for general nonlinear fourth-order parabolic equations on the one-dimensional torus are derived. Examples include the thin-film and the Derrida-Lebowitz-Speer-Spohn equations. The schemes…
Variational problems that involve Wasserstein distances have been recently proposed to summarize and learn from probability measures. Despite being conceptually simple, such problems are computationally challenging because they involve…
An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived, by approximating the Hamiltonian formulation, based on the Poisson brackets and the vorticity-divergence variables, of the inviscid shallow water…
We develop a new finite volume method using unstructured mesh-vertex grids for coupled systems modeling shallow water flows and solute transport over complex bottom topography. Novel well-balanced positivity preserving discretization…
Wasserstein gradient flows have become a central tool for optimization problems over probability measures. A natural numerical approach is forward-Euler time discretization. We show, however, that even in the simple case where the energy…
Two-fluid relativistic plasma flow equations combine the equations of relativistic hydrodynamics with Maxwell's equations for electromagnetic fields, which involve divergence constraints for the magnetic and electric fields. When developing…
In this paper, we start from a very natural system of cross-diffusion equations, which can be seen as a gradient flow for the Wasserstein distance of a certain functional. Unfortunately, the cross-diffusion system is not well-posed, as a…
We present a positivity-preserving method for multi-resolution simulations of compressible flows involving extreme conditions such as near vacuum and strong discontinuities. The novelty of this work is due to two aspects. First we extend…
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux…
Split form schemes for Euler and Navier-Stokes equations are useful for computation of turbulent flows due to their better robustness. This is because they satisfy additional conservation properties of the governing equations like kinetic…
This report presents a low computational and cognitive complexity, stable, time accurate and adaptive method for the Navier-Stokes equations. The improved method requires a minimally intrusive modification to an existing program based on…
For gradient flows, the existing structure-preserving schemes are difficult to achieve arbitrary high-order accuracy in time while preserving maximum-principle (MBP) and energy dissipating simultaneously. In this paper, we develop a new…