Related papers: A free energy satisfying finite difference method …
We present a discretization-free scalable framework for solving a large class of mass-conserving partial differential equations (PDEs), including the time-dependent Fokker-Planck equation and the Wasserstein gradient flow. The main…
We present a finite-volume based numerical scheme for a nonlocal Cahn-Hilliard equation which combines ideas from recent numerical schemes for gradient flow equations and nonlocal Cahn-Hilliard equations. The equation of interest is a…
For a general class of nonlinear port-Hamiltonian systems we develop a high-order time discretization scheme with certain structure preservation properties. The finite or infinite-dimensional system under consideration possesses a…
We propose a new Monte Carlo algorithm for the free energy calculation based on configuration space sampling. We implement this algorithm for Ising model. Comparison with the exact free energy shows an excellent agreement. We analyse the…
Well balanced and free energy dissipative first- and second-order accurate finite volume schemes are proposed for a general class of hydrodynamic systems with linear and nonlinear damping. The natural Liapunov functional of the system,…
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…
With a small parameter $\epsilon$, Poisson-Nernst-Planck (PNP) systems over a finite one-dimensional (1D) spatial domain have steady state solutions, called 1D boundary layer solutions, which profiles form boundary layers near boundary…
In this paper, we present a novel second order in time mixed finite element scheme for the Cahn-Hilliard-Navier-Stokes equations with matched densities. The scheme combines a standard second order Crank-Nicholson method for the…
In this paper, energy-preserving methods are formulated and studied for solving charged-particle dynamics. We first formulate the scheme of energy-preserving methods and analyze its basic properties including algebraic order and symmetry.…
We propose finite-volume schemes for general continuity equations which preserve positivity and global bounds that arise from saturation effects in the mobility function. In the case of gradient flows, the schemes dissipate the free energy…
We propose a new numerical method to solve the Cahn-Hilliard equation coupled with non-linear wetting boundary conditions. We show that the method is mass-conservative and that the discrete solution satisfies a discrete energy law similar…
This paper is concerned with finite element approximations of $W^{2,p}$ strong solutions of second-order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. A nonstandard (primal)…
We analyze a new framework for expressing finite element methods on arbitrarily many intersecting meshes: multimesh finite element methods. The multimesh finite element method, first presented in [40], enables the use of separate meshes to…
This paper detailedly discusses the locally one-dimensional numerical methods for efficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional…
In this paper, a linearized semi-implicit finite difference scheme is proposed for solving the two-dimensional (2D) space fractional nonlinear Schr\"{o}dinger equation (SFNSE).The scheme has the property of mass and energy conservation on…
In this paper, we study the biharmonic equation with the Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the 4th-order problem into a system of Poisson equations. Different…
We consider the linear Schr\"odinger equation and its discretization by split-step methods where the part corresponding to the Laplace operator is approximated by the midpoint rule. We show that the numerical solution coincides with the…
This paper concerns the construction and analysis of a numerical scheme for a mixed discrete-continuous fragmentation equation. A finite volume scheme is developed, based on a conservative formulation of a truncated version of the…
The parabolic-elliptic Keller-Segel equation with sensitivity saturation, because of its pattern formation ability, is a challenge for numerical simulations. We provide two finite-volume schemes whose goals are to preserve, at the discrete…
We study the low-energy solutions to the 3D compressible Navier-Stokes-Poisson equations. We first obtain the existence of smooth solutions with small $L^2$-norm and essentially bounded densities. No smallness assumption is imposed on the…