4.7 Article

Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: Theoretical analysis and application to the Vlasov-Poisson system

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 230, Issue 23, Pages 8386-8409

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2011.07.018

Keywords

Vlasov simulation; Discontinuous Galerkin method; Semi-Lagrangian method; Positivity preserving method; Characteristic/Lagrangian Galerkin method

Funding

  1. NSF [DMS-0914852, DMS-0809086]
  2. AFOSR [FA9550-09-1-0126]
  3. Division Of Mathematical Sciences
  4. Direct For Mathematical & Physical Scien [914852] Funding Source: National Science Foundation

Ask authors/readers for more resources

Semi-Lagrangian (SL) methods have been very popular in the Vlasov simulation community [29,3,4,18,32,5,24,26]. In this paper, we propose a new Strang split SL discontinuous Galerkin (DG) method for solving the Vlasov equation. Specifically, we apply the Strang splitting for the Vlasov equation [6], as a way to decouple the nonlinear Vlasov system into a sequence of 1-D advection equations, each of which has an advection velocity that only depends on coordinates that are transverse to the direction of propagation. To evolve the decoupled linear equations, we propose to couple the SL framework with the semi-discrete DG formulation. The proposed SL DC method is free of time step restriction compared with the Runge-Kutta DG method, which is known to suffer from numerical time step limitation with relatively small CFL numbers according to linear stability analysis. We apply the recently developed positivity preserving (PP) limiter [37], which is a low-cost black box procedure, to our scheme to ensure the positivity of the unknown probability density function without affecting the high order accuracy of the base SL DG scheme. We analyze the stability and accuracy properties of the SL DG scheme by establishing its connection with the direct and weak formulations of the characteristics/Lagrangian Galerkin method [23]. The quality of the proposed method is demonstrated via basic test problems, such as linear advection and rigid body rotation, and via classical plasma problems, such as Landau damping and the two stream instability. Published by Elsevier Inc.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

On a new centered strategy to control the accuracy of weighted essentially non oscillatory algorithm for conservation laws close to discontinuities

Sergio Amat, Antonio Baeza, Juan Ruiz, Chi-Wang Shu

Summary: This paper aims to redesign an algorithm previously analyzed by Amat et al. to achieve a centered strategy for approximating the solutions of linear and nonlinear systems of conservation laws. The analysis of the new algorithm explains its effectiveness near shocks in conservation law solutions and emphasizes the goal of using the most centered stencil in high gradient conditions.

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS (2021)

Article Engineering, Multidisciplinary

A high order positivity-preserving conservative WENO remapping method on 2D quadrilateral meshes

Nuo Lei, Juan Cheng, Chi-Wang Shu

Summary: This paper presents a high-order accurate positivity-preserving conservative remapping algorithm based on WENO reconstruction, which has been demonstrated through a series of numerical experiments to be conservative, positivity-preserving, and efficient.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Computer Science, Interdisciplinary Applications

Weighted ghost fluid discontinuous Galerkin method for two-medium problems

Yun-Long Liu, Chi-Wang Shu, A-Man Zhang

Summary: A new interface treatment method is proposed for simulating compressible two-medium problems using the RKDG method. The method ensures a smooth transition of the interface while minimizing overshoots or undershoots with the adoption of entropy-fix technique. It demonstrates high accuracy and compactness in handling interfaces with large entropy ratios.

JOURNAL OF COMPUTATIONAL PHYSICS (2021)

Article Mathematics, Applied

A discontinuous Galerkin method and its error estimate for nonlinear fourth-order wave equations

Qi Tao, Yan Xu, Chi-Wang Shu

Summary: This paper presents an ultra-weak local discontinuous Galerkin (UWLDG) method for a class of nonlinear fourth-order wave equations, designed and analyzed. The method demonstrates energy conserving properties and optimal error estimates, which are confirmed through numerical experiments. Compatible high order energy conserving time integrators are also proposed.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2021)

Article Mathematics, Applied

Stability analysis of inverse Lax-Wendroff boundary treatment of high order compact difference schemes for parabolic equations

Tingting Li, Jianfang Lu, Chi-Wang Shu

Summary: This paper investigates the stability of a numerical boundary treatment of high order compact finite difference methods for parabolic equations. The study utilizes the simplified inverse Lax-Wendroff procedure and third order TVD Runge-Kutta method, along with two analysis techniques to check algorithm stability, yielding consistent results in both semi-discrete and fully-discrete cases. Numerical experimental results are presented to validate the theoretical findings.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2022)

Article Mathematics, Applied

AN ESSENTIALLY OSCILLATION-FREE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC SYSTEMS

Yong Liu, Jianfang Lu, Chi-Wang Shu

Summary: In this paper, an essentially oscillation-free discontinuous Galerkin method is developed for systems of hyperbolic conservation laws. The method introduces numerical damping terms to control spurious oscillations. Both classical Runge-Kutta method and modified exponential Runge-Kutta method are used in time discretization. Extensive numerical experiments demonstrate the robustness and effectiveness of the algorithm.

SIAM JOURNAL ON SCIENTIFIC COMPUTING (2022)

Article Mathematics, Applied

Development and analysis of two new finite element schemes for a time-domain carpet cloak model

Jichun Li, Chi-Wang Shu, Wei Yang

Summary: This paper focuses on the time-domain carpet cloak model and proposes two new finite element schemes to address the numerical stability issue of previous schemes. The unconditional stability of the Crank-Nicolson scheme and the conditional stability of the leap-frog scheme are proved, both inheriting the exact form of the continuous stability.

ADVANCES IN COMPUTATIONAL MATHEMATICS (2022)

Article Engineering, Multidisciplinary

A high order positivity-preserving conservative WENO remapping method on 3D tetrahedral meshes

Nuo Lei, Juan Cheng, Chi-Wang Shu

Summary: We propose a high-order positivity-preserving conservative remapping method on 3D tetrahedral meshes based on the weighted essentially non-oscillatory reconstruction method. By accurately calculating the overlaps between meshes, our method allows for wider range of mesh movements and simplifies the remapping process. Utilizing the third order multi-resolution WENO reconstruction procedure, we distribute nonlinear weights based on the smoothness of polynomials to achieve optimal accuracy and avoid numerical oscillations. Our method also incorporates efficient local limiting to preserve positivity without compromising high-order accuracy and conservation. Numerical tests confirm the properties of our remapping algorithm, such as high-order accuracy, conservation, non-oscillatory performance, positivity-preserving, and efficiency.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2022)

Article Mathematics, Applied

Local discontinuous Galerkin methods for diffusive-viscous wave equations

Dan Ling, Chi-Wang Shu, Wenjing Yan

Summary: The paper focuses on the design of numerical methods for the diffusive-viscous wave equations with variable coefficients and develops a local discontinuous Galerkin (LDG) method. Numerical experiments are provided to demonstrate the optimal convergence rate and effectiveness of the proposed LDG method.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2023)

Article Computer Science, Interdisciplinary Applications

An improved simple WENO limiter for discontinuous Galerkin methods solving hyperbolic systems on unstructured meshes

Jie Du, Chi-Wang Shu, Xinghui Zhong

Summary: In this paper, we propose an improved simple WENO limiter for the Runge-Kutta discontinuous Galerkin method in solving two-dimensional hyperbolic systems on unstructured meshes. The major improvement is reducing the number of polynomials transformed to the characteristic fields for each direction, resulting in reduced computational cost and improved efficiency. The improved limiter provides a simpler and more practical way for the characteristic-wise limiting procedure, while maintaining uniform high-order accuracy in smooth regions and controlling nonphysical oscillations near discontinuities. Numerical results show that the improved limiter outperforms the original one in terms of accuracy and resolution.

JOURNAL OF COMPUTATIONAL PHYSICS (2022)

Article Mathematics, Applied

On the approximation of derivative values using a WENO algorithm with progressive order of accuracy close to discontinuities

Sergio Amat, Juan Ruiz-Alvarez, Chi-Wang Shu, Dionisio F. Yanez

Summary: This article introduces a new WENO algorithm for approximating derivative values of a function on a non-regular grid. The algorithm adapts ideas from a previous study to design nonlinear weights that maximize accuracy near discontinuities. The article provides proofs, discusses stencil selection, and presents explicit formulas for the weights and smoothness indicators. Numerical experiments are also conducted to validate the theoretical results.

COMPUTATIONAL & APPLIED MATHEMATICS (2022)

Article Computer Science, Interdisciplinary Applications

A primal-dual approach for solving conservation laws with implicit in time approximations

Siting Liu, Stanley Osher, Wuchen Li, Chi -Wang Shu

Summary: In this work, a novel framework for numerically solving time-dependent conservation laws with implicit schemes is proposed. The approach involves casting the initial value problem as a saddle point of a min-max problem and using iterative optimization methods to find the saddle point. The flexibility in the choice of time and spatial discretization schemes, as well as the large regions of stability gained from the implicit structure, make this approach advantageous. It is highly parallelizable and easy to implement, while avoiding non-linear inversions.

JOURNAL OF COMPUTATIONAL PHYSICS (2023)

Article Transportation

On the existence and uniqueness properties of the Hoogendoorn-Bovy pedestrian flow model

Liangze Yang, Chi-Wang Shu, S. C. Wong, Mengping Zhang, Jie Du

Summary: This study investigates the existence and uniqueness of solutions to the Hoogendoorn-Bovy (HB) pedestrian flow model. The results show that the HB model can be transformed into a forward conservation law equation and a backward Hamilton-Jacobi equation, ensuring the existence and uniqueness of solutions for both equations when suitable parameters are chosen.

TRANSPORTMETRICA B-TRANSPORT DYNAMICS (2023)

Article Mathematics, Applied

A Fixed-Point Fast Sweeping WENO Method with Inverse Lax-Wendroff Boundary Treatment for Steady State of Hyperbolic Conservation Laws

Liang Li, Jun Zhu, Chi-Wang Shu, Yong-Tao Zhang

Summary: Fixed-point fast sweeping WENO methods are efficient high-order numerical methods used to solve steady-state solutions of hyperbolic PDEs. They have high-order accuracy and good performance.

COMMUNICATIONS ON APPLIED MATHEMATICS AND COMPUTATION (2023)

Article Physics, Mathematical

Moment-Based Multi-Resolution HWENO Scheme for Hyperbolic Conservation Laws

Jiayin Li, Chi-Wang Shu, Jianxian Qiu

Summary: In this paper, a high-order moment-based multi-resolution HWENO scheme is proposed for hyperbolic conservation laws. The scheme reconstructs the function values at Gauss-Lobatto points using the information of the zeroth and first order moments, leading to improved stability and resolution. Compared to general HWENO and WENO schemes, this moment-based scheme has a more compact stencil size and a higher CFL number restriction.

COMMUNICATIONS IN COMPUTATIONAL PHYSICS (2022)

Article Computer Science, Interdisciplinary Applications

A new type of non-polynomial based TENO scheme for hyperbolic conservation laws

Tian Liang, Lin Fu

Summary: In this work, a new shock-capturing framework is proposed based on a new candidate stencil arrangement and the combination of infinitely differentiable non-polynomial RBF-based reconstruction in smooth regions with jump-like non-polynomial interpolation for genuine discontinuities. The resulting scheme achieves high order accuracy and resolves genuine discontinuities with sub-cell resolution.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A high-order residual-based viscosity finite element method for incompressible variable density flow

Lukas Lundgren, Murtazo Nazarov

Summary: In this paper, a high-order accurate finite element method for incompressible variable density flow is introduced. The method addresses the issues of saddle point system and stability problem through Schur complement preconditioning and artificial compressibility approaches, and it is validated to have high-order accuracy for smooth problems and accurately resolve discontinuities.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Convergence analysis and optimization of a Robin Schwarz waveform relaxation method for time-periodic parabolic optimal control problems

Gabriele Ciaramella, Laurence Halpern, Luca Mechelli

Summary: This paper presents a novel convergence analysis of the optimized Schwarz waveform relaxation method for solving optimal control problems governed by periodic parabolic PDEs. The analysis is based on a Fourier-type technique applied to a semidiscrete-in-time form of the optimality condition, which enables a precise characterization of the convergence factor at the semidiscrete level. The behavior of the optimal transmission condition parameter is also analyzed in detail as the time discretization approaches zero.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Data-driven Whitney forms for structure-preserving control volume analysis

Jonas A. Actor, Xiaozhe Hu, Andy Huang, Scott A. Roberts, Nathaniel Trask

Summary: This article introduces a scientific machine learning framework that uses a partition of unity architecture to model physics through control volume analysis. The framework can extract reduced models from full field data while preserving the physics. It is applicable to manifolds in arbitrary dimension and has been demonstrated effective in specific problems.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Higher-continuity s-version of finite element method with B-spline functions

Nozomi Magome, Naoki Morita, Shigeki Kaneko, Naoto Mitsume

Summary: This paper proposes a novel strategy called B-spline based SFEM to fundamentally solve the problems of the conventional SFEM. It uses different basis functions and cubic B-spline basis functions with C-2-continuity to improve the accuracy of numerical integration and avoid matrix singularity. Numerical results show that the proposed method is superior to conventional methods in terms of accuracy and convergence.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A cell-centred Eulerian volume-of-fluid method for compressible multi-material flows

Timothy R. Law, Philip T. Barton

Summary: This paper presents a practical cell-centred volume-of-fluid method for simulating compressible solid-fluid problems within a pure Eulerian setting. The method incorporates a mixed-cell update to maintain sharp interfaces, and can be easily extended to include other coupled physics. Various challenging test problems are used to validate the method, and its robustness and application in a multi-physics context are demonstrated.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Two-step multi-resolution reconstruction-based compact gas-kinetic scheme on tetrahedral mesh

Xing Ji, Fengxiang Zhao, Wei Shyy, Kun Xu

Summary: This paper presents the development of a third-order compact gas-kinetic scheme for compressible Euler and Navier-Stokes solutions, constructed particularly for an unstructured tetrahedral mesh. The scheme demonstrates robustness in high-speed flow computation and exhibits excellent adaptability to meshes with complex geometrical configurations.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Multiscale sampling for the inverse modeling of partial differential equations

Alsadig Ali, Abdullah Al-Mamun, Felipe Pereira, Arunasalam Rahunanthan

Summary: This paper presents a novel Bayesian statistical framework for the characterization of natural subsurface formations, and introduces the concept of multiscale sampling to localize the search in the stochastic space. The results show that the proposed framework performs well in solving inverse problems related to porous media flows.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Constrained optimized dynamic mode decomposition with control for physically stable systems with exogeneous inputs

Jacob Rains, Yi Wang, Alec House, Andrew L. Kaminsky, Nathan A. Tison, Vamshi M. Korivi

Summary: This paper presents a novel method called constrained optimized DMD with Control (cOptDMDc), which extends the optimized DMD method to systems with exogenous inputs and can enforce the stability of the resulting reduced order model (ROM). The proposed method optimally places eigenvalues within the stable region, thus mitigating spurious eigenvalue issues. Comparative studies show that cOptDMDc achieves high accuracy and robustness.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A hybridizable discontinuous Galerkin formulation for the Euler-Maxwell plasma model

Andrea La Spina, Jacob Fish

Summary: This work introduces a hybridizable discontinuous Galerkin formulation for simulating ideal plasmas. The proposed method couples the fluid and electromagnetic subproblems monolithically based on source and employs a fully implicit time integration scheme. The approach also utilizes a projection-based divergence correction method to enforce the Gauss laws in challenging scenarios. Numerical examples demonstrate the high-order accuracy, efficiency, and robustness of the proposed formulation.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Numerical solution of the cavity scattering problem for flexural waves on thin plates: Linear finite element methods

Junhong Yue, Peijun Li

Summary: This paper proposes two numerical methods (IP-FEM and BP-FEM) to study the flexural wave scattering problem of an arbitrary-shaped cavity on an infinite thin plate. These methods successfully decompose the fourth-order plate wave equation into the Helmholtz and modified Helmholtz equations with coupled conditions on the cavity boundary, providing an effective solution to this challenging problem.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Fast and scalable computation of shape-morphing nonlinear solutions with application to evolutional neural networks

William Anderson, Mohammad Farazmand

Summary: We develop fast and scalable methods, called RONS, for computing reduced-order nonlinear solutions. These methods have been proven to be highly effective in tackling challenging problems, but become computationally prohibitive as the number of parameters grows. To address this issue, three separate methods are proposed and their efficacy is demonstrated through examples. The application of RONS to neural networks is also discussed.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A second order directional split exponential integrator for systems of advection-diffusion-reaction equations

Marco Caliari, Fabio Cassini

Summary: In this paper, a second order exponential scheme for stiff evolutionary advection-diffusion-reaction equations is proposed. The scheme is based on a directional splitting approach and uses computation of small sized exponential-like functions and tensor-matrix products for efficient implementation. Numerical examples demonstrate the advantage of the proposed approach over state-of-the-art techniques.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A conservative semi-Lagrangian method for inhomogeneous Boltzmann equation

Sebastiano Boscarino, Seung Yeon Cho, Giovanni Russo

Summary: This work proposes a high order conservative semi-Lagrangian method for the inhomogeneous Boltzmann equation of rarefied gas dynamics. The method combines a semi-Lagrangian scheme for the convection term, a fast spectral method for computation of the collision operator, and a high order conservative reconstruction and a weighted optimization technique to preserve conservative quantities. Numerical tests demonstrate the accuracy and efficiency of the proposed method.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Uniqueness and numerical scheme for spherical shell-structured sources from the far field patterns with at most two frequencies

Jialei Li, Xiaodong Liu, Qingxiang Shi

Summary: This study shows that the number, centers, scattering strengths, inner and outer diameters of spherical shell-structured sources can be uniquely determined from the far field patterns. A numerical scheme is proposed for reconstructing the spherical shell-structured sources, which includes a migration series method for locating the centers and an iterative method for computing the inner and outer diameters without computing derivatives.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)