Related papers: High order asymptotic preserving scheme for linear…
Symmetry preserving difference schemes approximating second and third order ordinary differential equations are presented. They have the same three or four-dimensional symmetry groups as the original differential equations. The new…
In this paper, we propose a general framework to design asymptotic preserving schemes for the Boltzmann kinetic kinetic and related equations. Numerically solving these equations are challenging due to the nonlinear stiff collision (source)…
The aim of this paper is to apply a high-order discontinuous-in-time scheme to second-order hyperbolic partial differential equations (PDEs). We first discretize the PDEs in time while keeping the spatial differential operators…
We investigate time complexities of finite difference methods for solving the high-dimensional linear heat equation, the high-dimensional linear hyperbolic equation and the multiscale hyperbolic heat system with quantum algorithms (hence…
In magnetized plasma, the magnetic field confines the particles around the field lines. The anisotropy intensity in the viscosity and heat conduction may reach the order of $10^{12}$. When the boundary conditions are periodic or Neumann,…
Convergence properties of model-free two-timescale asymptotic simulations of singularly perturbed hybrid inclusions are developed. A hybrid inclusion combines constrained differential and difference inclusions to capture continuous (flow)…
A fully discrete finite element method, based on a new weak formulation and a new time-stepping scheme, is proposed for the surface diffusion flow of closed curves in the two-dimensional plane. It is proved that the proposed method can…
In contrast with the diffusion equation which smoothens the initial data to $C^\infty$ for $t>0$ (away from the corners/edges of the domain), the subdiffusion equation only exhibits limited spatial regularity. As a result, one generally…
We propose an efficient numerical strategy for simulating fluid flow through porous media with highly oscillatory characteristics. Specifically, we consider non-linear diffusion models. This scheme is based on the classical homogenization…
An efficient, iterative semi-implicit (SI) numerical method for the time integration of stiff wave systems is presented. Physics-based assumptions are used to derive a convergent iterative formulation of the SI scheme which enables the…
The convective Allen-Cahn equation has been widely used to simulate multi-phase flows in many phase-field models. As a generalized form of the classic Allen-Cahn equation, the convective Allen-Cahn equation still preserves the maximum bound…
Discontinuous Galerkin (DG) methods are widely adopted to discretize the radiation transport equation (RTE) with diffusive scalings. One of the most important advantages of the DG methods for RTE is their asymptotic preserving (AP)…
Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems.Particular attention is paid to the case when first order…
In this article we numerically revisit the long-time solution behavior of the Camassa-Holm equation. The finite difference solution of this integrable equation is sought subject to the newly derived initial condition with Delta-function…
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…
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…
High-speed turbulent flows are encountered in most space-related applications (including exploration, tourism and defense fields) and represent a subject of growing interest in the last decades. A major challenge in performing high-fidelity…
In this paper the analysis of an asymptotic preserving (AP) IMEX-RK finite volume scheme for the wave equation system in the zero Mach number limit is presented. The accuracy of a numerical scheme at low Mach numbers is its ability to…
Anisotropic diffusion is a well recognized tool in digital image processing, including edge detection and denoising. We present here a particular nonlinear time-dependent operator together with an appropriate high-order discretization for…
High order schemes are known to be unstable in the presence of shock discontinuities or under-resolved solution features for nonlinear conservation laws. Entropy stable schemes address this instability by ensuring that physically relevant…