98 research outputs found
High Order Maximum Principle Preserving Semi-Lagrangian Finite Difference WENO schemes for the Vlasov Equation
In this paper, we propose the parametrized maximum principle preserving (MPP)
flux limiter, originally developed in [Z. Xu, Math. Comp., (2013), in press],
to the semi- Lagrangian finite difference weighted essentially non-oscillatory
scheme for solving the Vlasov equation. The MPP flux limiter is proved to
maintain up to fourth order accuracy for the semi-Lagrangian finite difference
scheme without any time step restriction. Numerical studies on the
Vlasov-Poisson system demonstrate the performance of the proposed method and
its ability in preserving the positivity of the probability distribution
function while maintaining the high order accuracy
Positivity-Preserving Finite Difference WENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations
In this paper, we utilize the maximum-principle-preserving flux limiting
technique, originally designed for high order weighted essentially
non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to
develop a class of high order positivity-preserving finite difference WENO
methods for the ideal magnetohydrodynamic (MHD) equations. Our schemes, under
the constrained transport (CT) framework, can achieve high order accuracy, a
discrete divergence-free condition and positivity of the numerical solution
simultaneously. Numerical examples in 1D, 2D and 3D are provided to demonstrate
the performance of the proposed method.Comment: 21 pages, 28 figure
Maximum principle preserving high order schemes for convection-dominated diffusion equations
The maximum principle is an important property of solutions to PDE. Correspondingly, it\u27s of great interest for people to design a high order numerical scheme solving PDE with this property maintained. In this thesis, our particular interest is solving convection-dominated diffusion equation. We first review a nonconventional maximum principle preserving(MPP) high order finite volume(FV) WENO scheme, and then propose a new parametrized MPP high order finite difference(FD) WENO framework, which is generalized from the one solving hyperbolic conservation laws. A formal analysis is presented to show that a third order finite difference scheme with this parametrized MPP flux limiters maintains the third order accuracy without extra CFL constraint when the low order monotone flux is chosen appropriately. Numerical tests in both one and two dimensional cases are performed on the simulation of the incompressible Navier-Stokes equations in vorticity stream-function formulation and several other problems to show the effectiveness of the proposed method
An explicit high-order single-stage single-step positivity-preserving finite difference WENO method for the compressible Euler equations
In this work we construct a high-order, single-stage, single-step
positivity-preserving method for the compressible Euler equations. Space is
discretized with the finite difference weighted essentially non-oscillatory
(WENO) method. Time is discretized through a Lax-Wendroff procedure that is
constructed from the Picard integral formulation (PIF) of the partial
differential equation. The method can be viewed as a modified flux approach,
where a linear combination of a low- and high-order flux defines the numerical
flux used for a single-step update. The coefficients of the linear combination
are constructed by solving a simple optimization problem at each time step. The
high-order flux itself is constructed through the use of Taylor series and the
Cauchy-Kowalewski procedure that incorporates higher-order terms. Numerical
results in one- and two-dimensions are presented
High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics
The paper develops high-order accurate physical-constraints-preserving finite
difference WENO schemes for special relativistic hydrodynamical (RHD)
equations, built on the local Lax-Friedrich splitting, the WENO reconstruction,
the physical-constraints-preserving flux limiter, and the high-order strong
stability preserving time discretization. They are extensions of the
positivity-preserving finite difference WENO schemes for the non-relativistic
Euler equations. However, developing physical-constraints-preserving methods
for the RHD system becomes much more difficult than the non-relativistic case
because of the strongly coupling between the RHD equations, no explicit
expressions of the primitive variables and the flux vectors, in terms of the
conservative vector, and one more physical constraint for the fluid velocity in
addition to the positivity of the rest-mass density and the pressure. The key
is to prove the convexity and other properties of the admissible state set and
discover a concave function with respect to the conservative vector replacing
the pressure which is an important ingredient to enforce the
positivity-preserving property for the non-relativistic case. Several one- and
two-dimensional numerical examples are used to demonstrate accuracy,
robustness, and effectiveness of the proposed physical-constraints-preserving
schemes in solving RHD problems with large Lorentz factor, or strong
discontinuities, or low rest-mass density or pressure etc.Comment: 39 pages, 13 figure
- …