Conservative Semi-Lagrangian Finite Difference WENO Formulations with Applications to the Vlasov Equation

被引:46
作者
Qiu, Jing-Mei [1 ]
Shu, Chi-Wang [2 ]
机构
[1] Colorado Sch Mines, Dept Math & Comp Sci, Golden, CO 80401 USA
[2] Brown Univ, Div Appl Math, Providence, RI 02912 USA
基金
美国国家科学基金会;
关键词
Semi-Lagrangian methods; finite difference/finite volume scheme; conservative scheme; WENO reconstruction; Vlasov equation; Landau damping; two-stream instability; SCHEMES; INTERPOLATION; LAWS;
D O I
10.4208/cicp.180210.251110a
中图分类号
O4 [物理学];
学科分类号
0702 ;
摘要
In this paper, we propose a new conservative semi-Lagrangian (SL) finite difference (FD) WENO scheme for linear advection equations, which can serve as a base scheme for the Vlasov equation by Strang splitting [4]. The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume (FV) WENO scheme [3]. However, instead of inputting cell averages and approximate the integral forth of the equation in a FV scheme, we input point values and approximate the differential form of equation in a FD spirit, yet retaining very high order (fifth order in our experiment) spatial accuracy. The advantage of using point values, rather than cell averages, is to avoid the second order spatial error, due to the shearing in velocity (v) and electrical field (E) over a cell when performing the Strang splitting to the Vlasov equation. As a result, the proposed scheme has very high spatial accuracy, compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson (VP) system. We perform numerical experiments on linear advection, rigid body rotation problem; and on the Landau damping and two-stream instabilities by solving the VP system. For comparison, we also apply (1) the conservative SL FD WENO scheme, proposed in [22] for incompressible advection problem, (2) the conservative SL FD WENO scheme proposed in [21] and (3) the non-conservative version of the SL FD WENO scheme in [3] to the same test problems. The performances of different schemes are compared by the error table, solution resolution of sharp interface, and by tracking the conservation of physical norms, energies and entropies, which should be physically preserved.
引用
收藏
页码:979 / 1000
页数:22
相关论文
共 29 条
[1]  
Boyd JohnP, 2001, CHEBYSHEV FOURIER SP
[2]   A weighted essentially nonoscillatory, large time-step scheme for Hamilton -Jacobi equations [J].
Carlini, E ;
Ferretti, R ;
Russo, G .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2005, 27 (03) :1071-1091
[3]   Nonoscillatory interpolation methods applied to Vlasov-based models [J].
Carrillo, J. A. ;
Vecil, F. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2007, 29 (03) :1179-1206
[4]   INTEGRATION OF VLASOV EQUATION IN CONFIGURATION SPACE [J].
CHENG, CZ ;
KNORR, G .
JOURNAL OF COMPUTATIONAL PHYSICS, 1976, 22 (03) :330-351
[5]   CHARACTERISTIC GALERKIN METHODS FOR SCALAR CONSERVATION-LAWS IN ONE DIMENSION [J].
CHILDS, PN ;
MORTON, KW .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1990, 27 (03) :553-594
[6]  
Cockburn Bernardo, 1998, Advanced Numerical Approximation of Nonlinear Hyperbolic Equations
[7]   THE PIECEWISE PARABOLIC METHOD (PPM) FOR GAS-DYNAMICAL SIMULATIONS [J].
COLELLA, P ;
WOODWARD, PR .
JOURNAL OF COMPUTATIONAL PHYSICS, 1984, 54 (01) :174-201
[8]   Hermite spline interpolation on patches for parallelly solving the Vlasov-Poisson equation [J].
Crouseilles, Nicolas ;
Latu, Guillaume ;
Sonnendruecker, Eric .
INTERNATIONAL JOURNAL OF APPLIED MATHEMATICS AND COMPUTER SCIENCE, 2007, 17 (03) :335-349
[9]   Conservative semi-Lagrangian schemes for Vlasov equations [J].
Crouseilles, Nicolas ;
Mehrenberger, Michel ;
Sonnendruecker, Eric .
JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 229 (06) :1927-1953
[10]   Comparison of Eulerian Vlasov solvers [J].
Filbet, F ;
Sonnendrücker, E .
COMPUTER PHYSICS COMMUNICATIONS, 2003, 150 (03) :247-266