Related papers: Strang splitting structure-preserving high-order c…
Many multiscale problems have a high contrast, which is expressed as a very large ratio between the media properties. The contrast is known to introduce many challenges in the design of multiscale methods and domain decomposition…
In this paper, we propose a high order residual distribution conservative finite difference scheme for solving steady state conservation laws. A new type of WENO (weighted essentially non-oscillatory) termed as WENO-ZQ integration is used…
In the past decades, the finite difference methods for space fractional operators develop rapidly; to the best of our knowledge, all the existing finite difference schemes, including the first and high order ones, just work on uniform…
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment…
We propose and analyze a second-order Strang splitting method for a class of stiff matrix differential equations with Sylvester-type structure. The method splits the dynamics into a stiff linear part, treated exactly via matrix…
An efficient numerical scheme for solving transport equations for tokamak plasmas within an integrated modelling framework is presented. The plasma transport equations are formulated as diffusion-advection equations in two coordinates (a…
In this paper, based on the developed nonlinear fourth-order operator and method of order reduction, a novel fourth-order compact difference scheme is constructed for the mixed-type time-fractional Burgers' equation, from which…
In this work we propose a nonlinear stabilization technique for convection-diffusion-reaction and pure transport problems discretized with space-time isogeometric analysis. The stabilization is based on a graph-theoretic artificial…
Our goal is to develop a flux limiter of the Flux-Corrected Transport method for a nonconservative convection-diffusion equation. For this, we consider a hybrid difference scheme that is a linear combination of a monotone scheme and a…
This paper is devoted to the analysis of a numerical scheme for the coagulation and fragmentation equation with diffusion in space. A finite volume scheme is developed, based on a conservative formulation of the space nonhomogeneous…
A series of third- and fifth-order hybrid compact least-squares central weighted essentially non-oscillatory schemes are proposed and applied to curvilinear structured grids for the finite volume method. In smooth regions, compact…
In this paper, we introduce a Lagrange multiplier approach to construct linearly implicit energy-preserving schemes of arbitrary order for general Hamiltonian PDEs. Unlike the widely used auxiliary variable methods, this novel approach does…
We extend to multi-dimensions the work of [1], where new fully explicit kinetic methods were built for the approximation of linear and non-linear convection-diffusion problems. The fundamental principles from the earlier work are retained:…
Simulating infiltration in porous media using Richards' equation remains computationally challenging due to its parabolic structure and nonlinear coefficients. While a wide range of numerical methods for differential equations have been…
A fourth-order compact scheme is proposed for a fourth-order subdiffusion equation with the first Dirichlet boundary conditions. The fourth-order problem is firstly reduced into a couple of spatially second-order system and we use an…
We extend to the N-level Bloch model the splitting scheme which use exact numerical solutions of sub-equations. These exact solutions involve matrix exponentials which we want to avoid to calculate at each time step. We use Newton…
In this paper, we present a class of nonuniform time-stepping, high-order linear stabilized schemes that can preserve both the discrete energy stability and maximum-bound principle (MBP) for the time-fractional Allen-Cahn equation. To this…
A finite difference scheme is presented for the Dirac equation in (1+1)D. It can handle space- and time-dependent mass and potential terms and utilizes exact discrete transparent boundary conditions (DTBCs). Based on a space- and…
We consider an initial-boundary value problem for a 2D time-dependent Schr\"odinger equation on a semi-infinite strip. For the Numerov-Crank-Nicolson finite-difference scheme with discrete transparent boundary conditions, the Strang-type…
The isentropic compressible Cahn-Hilliard-Navier-Stokes equations is a system of fourth-order partial differential equations that model the evolution of some binary fluids under convection. The purpose of this paper is the design of…