Search

Semi-Lagrangian (SL) schemes are highly efficient for simulating transport equations and are widely used across various applications. Despite their success, designing genuinely multi-dimensional and conservative SL schemes remains a significant challenge. Building on our previous work [Chen et al., J. Comput. Phys., V490 112329, (2023)], we introduce a conservative machine-learning-based SL finite difference (FD) method that allows for extra-large time step evolution. At the core of our approach is a novel dynamical graph neural network designed to handle the complexities associated with tracking accurately upstream points along characteristics. This proposed neural transport solver learns the conservative SL FD discretization directly from data, improving accuracy and efficiency compared to traditional numerical schemes, while significantly simplifying algorithm implementation. We validate the method' s effectiveness and efficiency through numerical tests on benchmark transport equations in both one and two dimensions, as well as the nonlinear Vlasov-Poisson system., Comment: arXiv admin note: text overlap with arXiv:2309.04943

Machine-learning (ML) based discretization has been developed to simulate complex partial differential equations (PDEs) with tremendous success across various fields. These learned PDE solvers can effectively resolve the underlying solution structures of interest and achieve a level of accuracy which often requires an order-of-magnitude finer grid for a conventional numerical method using polynomial-based approximations. In a previous work in [13], we introduced a learned finite volume discretization that further incorporates the semi-Lagrangian (SL) mechanism, enabling larger CFL numbers for stability. However, the efficiency and effectiveness of such methodology heavily rely on the availability of abundant high-resolution training data, which can be prohibitively expensive to obtain. To address this challenge, in this paper, we propose a novel multi-fidelity ML-based SL method for transport equations. This method leverages a combination of a small amount of high-fidelity data and sufficient but cheaper low-fidelity data. The approach is designed based on a composite convolutional neural network architecture that explore the inherent correlation between high-fidelity and low-fidelity data. The proposed method demonstrates the capability to achieve a reasonable level of accuracy, particularly in scenarios where a single-fidelity model fails to generalize effectively. We further extend the method to the nonlinear Vlasov-Poisson system by employing high order Runge-Kutta exponential integrators. A collection of numerical tests are provided to validate the efficiency and accuracy of the proposed method.

Semi-Lagrangian (SL) schemes are known as a major numerical tool for solving transport equations with many advantages and have been widely deployed in the fields of computational fluid dynamics, plasma physics modeling, numerical weather prediction, among others. In this work, we develop a novel machine learning-assisted approach to accelerate the conventional SL finite volume (FV) schemes. The proposed scheme avoids the expensive tracking of upstream cells but attempts to learn the SL discretization from the data by incorporating specific inductive biases in the neural network, significantly simplifying the algorithm implementation and leading to improved efficiency. In addition, the method delivers sharp shock transitions and a level of accuracy that would typically require a much finer grid with traditional transport solvers. Numerical tests demonstrate the effectiveness and efficiency of the proposed method., Comment: 24 pages

We present an energy-conserving numerical scheme to solve the Vlasov-Maxwell (VM) system based on the regularized moment method proposed in [Z. Cai, Y. Fan, and R. Li. CPAM, 2014]. The globally hyperbolic moment system is deduced for the multi-dimensional VM system under the framework of the Hermite expansions, where the expansion center and the scaling factor are set as the macroscopic velocity and local temperature, respectively. Thus, the effect of the Lorentz force term could be reduced into several ODEs about the macroscopic velocity and the moment coefficients of higher order, which could significantly reduce the computational cost of the whole system. An energy-conserving numerical scheme is proposed to solve the moment equations and the Maxwell equations, where only a linear equation system needs to be solved. Several numerical examples such as the two-stream instability, Weibel instability, and the two-dimensional Orszag Tang vortex problem are studied to validate the efficiency and excellent energy-preserving property of the numerical scheme.

In this paper, we develop two finite difference weighted essentially non-oscillatory (WENO) schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi (DP) and $\mu$-Degasperis-Procesi ($\mu$DP) equations, which contain nonlinear high order derivatives, and possibly peakon solutions or shock waves. By introducing auxiliary variable(s), we rewrite the DP equation as a hyperbolic-elliptic system, and the \mdp equation as a first order system. Then we choose a linear finite difference scheme with suitable order of accuracy for the auxiliary variable(s), and two finite difference WENO schemes with unequal-sized sub-stencils for the primal variable. One WENO scheme uses one large stencil and several smaller stencils, and the other WENO scheme is based on the multi-resolution framework which uses a series of unequal-sized hierarchical central stencils. Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil, both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights. Another advantage is that the final reconstructed polynomial on the target cell is a polynomial of the same degree as the polynomial over the big stencil, while the classical finite difference WENO reconstruction can only be obtained for specific points inside the target interval. Numerical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.

We develop a numerical framework to implement the cumulative density function (CDF) method for obtaining the probability distribution of the system state described by a kinematic wave model. The approach relies on Monte Carlo Simulations (MCS) of the fine-grained CDF equation of system state, as derived by the CDF method. This fine-grained CDF equation is solved via the method of characteristics. Each method of characteristics solution is far more computationally efficient than the direct solution of the kinematic wave model, and the MCS estimator of the CDF converges relatively quickly. We verify the accuracy and robustness of our procedure via comparison with direct MCS of a particular kinematic wave system, the Saint-Venant equation., Comment: 19 pages, 6 figures, 3 tables

In this paper, we propose energy-conserving Eulerian solvers for the two-species Vlasov-Amp\`{e}re (VA) system and apply the methods to simulate current-driven ion-acoustic instability. The algorithm is generalized from our previous work for the single-species VA system and Vlasov-Maxwell (VM) system. The main feature of the schemes is their ability to preserve the total particle number and total energy on the fully discrete level regardless of mesh size. Those are desired properties of numerical schemes especially for long time simulations with under-resolved mesh. The conservation is realized by explicit and implicit energy-conserving temporal discretizations, and the discontinuous Galerkin (DG) spatial discretizations. We benchmarked our algorithms on a test example to check the one-species limit, and the current-driven ion-acoustic instability. To simulate the current-driven ion-acoustic instability, a slight modification for the implicit method is necessary to fully decouple the split equations. This is achieved by a Gauss-Seidel type iteration technique. Numerical results verified the conservation and performance of our methods.

In this paper, we generalize the idea in our previous work for the Vlasov-Amp\`{e}re (VA) system \cite{cheng_va} and develop energy-conserving discontinuous Galerkin (DG) methods for the Vlasov-Maxwell (VM) system. The VM system is a fundamental model in the simulation of collisionless magnetized plasmas. Compared to \cite{cheng_va}, additional care needs to be taken for both the temporal and spatial discretizations to achieve similar type of conservation when the magnetic field is no longer negligible. Our proposed schemes conserve the total particle number and the total energy at the same time, and therefore can obtain accurate, yet physically relevant solutions. The main components of our methods include second order and above, explicit or implicit energy-conserving temporal discretizations, and DG methods for Vlasov and Maxwell's equations with carefully chosen numerical fluxes. Benchmark numerical tests such as the streaming Weibel instability are provided to validate the accuracy and conservation of the schemes., Comment: arXiv admin note: substantial text overlap with arXiv:1306.0931 Add more numerical results on page 28

In this paper, we propose energy-conserving numerical schemes for the Vlasov-Amp\`{e}re (VA) systems. The VA system is a model used to describe the evolution of probability density function of charged particles under self consistent electric field in plasmas. It conserves many physical quantities, including the total energy which is comprised of the kinetic and electric energy. Unlike the total particle number conservation, the total energy conservation is challenging to achieve. For simulations in longer time ranges, negligence of this fact could cause unphysical results, such as plasma self heating or cooling. In this paper, we develop the first Eulerian solvers that can preserve fully discrete total energy conservation. The main components of our solvers include explicit or implicit energy-conserving temporal discretizations, an energy-conserving operator splitting for the VA equation and discontinuous Galerkin finite element methods for the spatial discretizations. We validate our schemes by rigorous derivations and benchmark numerical examples such as Landau damping, two-stream instability and bump-on-tail instability.