Related papers: Positivity preserving high order schemes for angio…
We introduce the first-order system proximal Galerkin (FOSPG) method, a locally mass-conserving, hybridizable finite element method for solving heterogeneous anisotropic diffusion and obstacle problems. Like other proximal Galerkin methods,…
In this paper, we propose a class of explicit positivity preserving numerical methods for general stochastic differential equations which have positive solutions. Namely, all the numerical solutions are positive. Under some reasonable…
We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is…
Godunov type numerical schemes for the class of hyperbolic systems, admitting non-classical $\delta-$ shocks are proposed. It is shown that the numerical approximations converge to the solution and preserve the physical properties of the…
We present an asymptotic preserving method for the radiative transfer equations in the framework of PN method. An implicit and explicit method is proposed to solve the P N system based on the order analysis of the expansion coefficients of…
We aim to solve a topology optimization problem where the distribution of material in the design domain is represented by a density function. To obtain candidates for local minima, we want to solve the first order optimality system via…
We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general…
The first-order linear positivity preserving schemes in time are available for the time dependent Poisson-Nernst-Planck (PNP) equations, second-order linear ones are still challenging. In this paper, we propose the first- and second-order…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics, via solving the high-dimensional Fokker-Planck equation in terms of low-order cluster moments or…
Here, we introduce a numerical approach for a class of Fokker-Planck (FP) equations. These equations are the adjoint of the linearization of Hamilton-Jacobi (HJ) equations. Using this structure, we show how to transfer the properties of…
We analyze an advection-diffusion-reaction problem with non-homogeneous boundary conditions that models the chromatography process, a vital stage in bioseparation. We prove stability and error estimates for both constant and affine…
We present high-order variational Lagrangian finite element methods for compressible fluids using a discrete energetic variational approach. Our spatial discretization is mass/momentum/energy conserving and entropy stable. Fully implicit…
The thermal radiative transfer equations model temperature evolution through a background medium as a result of radiation. When a large number of particles are absorbed in a short time scale, the dynamics tend to a non-linear diffusion-type…
This paper presents a general positivity-preserving algorithm for implicit high-order finite volume schemes solving Euler and Navier-Stokes equations. Previous positivity-preserving algorithms are mainly based on mathematical analyses,…
In this work, we develop novel structure-preserving numerical schemes for a class of nonlinear Fokker--Planck equations with nonlocal interactions. Such equations can cover many cases of importance, such as porous medium equations with…
In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method keeps the original scheme unchanged and detects critical numerical fluxes which may lead to…
The main aim of this paper is to document the performance of $p$-refinement with respect to maximum principles and the non-negative constraint. The model problem is (steady-state) anisotropic diffusion with decay (which is a second-order…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
We study the numerical approximation of advection-diffusion equations with highly oscillatory coefficients and possibly dominant advection terms by means of the Multiscale Finite Element Method. The latter method is a now classical, finite…