Abstract:Active Flux is a modified Finite Volume method that evolves additional Degrees of Freedom for each cell that are located on the interface by a non-conservative method to compute high-order approximations to the numerical fluxes through the respective interface to evolve the cell-average in a conservative way. In this paper, we apply the method to the Vlasov-Poisson system describing the time evolution of the time-dependent distribution function of a collisionless plasma. In particular, we consider the evaluation of the flux integrals in higher dimensions. We propose a dimensional splitting and three types of formulations of the flux integral: a one-dimensional reconstruction of second order, a third-order reconstruction based on information along each dimension, and a third-order reconstruction based on a discrepancy formulation of the Active Flux method. Numerical results in 1D1V phase-space compare the properties of the various methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an efficient and accurate numerical method for solving the Vlasov - Poisson system that describes the time evolution of collisionless plasmas. Specifically, the paper introduces a method called split - step Active Flux, aiming at performing highly efficient and high - precision numerical simulations in high - dimensional phase space.
### Core of the problem
1. **High - dimensionality**: The Vlasov equation is a six - dimensional problem (three - dimensional position and three - dimensional velocity space plus time), so the numerical simulation has a very high computational cost.
2. **Large time - step limitation**: Although the Vlasov equation is formally a relatively simple transport equation, since the components of the velocity \(v\) can be very large, this will lead to severe time - step limitations.
3. **Preservation of physical properties**: The numerical scheme should preserve the positivity and conservation of the distribution function \(f_s\) to ensure the correctness of the physical results.
4. **Energy conservation**: During the numerical discretization process, it is necessary to ensure that the total energy (including particle kinetic energy and electromagnetic field energy) remains constant over time.
5. **Constraints**: The numerical solution should satisfy the constraint conditions of the electromagnetic field, such as \(\nabla\cdot E=\rho / \epsilon_0\) and \(\nabla\cdot B = 0\).
### Solution
The method proposed in the paper combines the Active Flux method and the splitting technique, specifically including:
- **Active Flux method**: This is an improved finite - volume method. By introducing additional degrees of freedom (DOF) on the cell interface, the numerical flux can be approximated with high - order accuracy, thus evolving the cell averages in a conservative manner.
- **Splitting technique**: To deal with high - dimensional problems, the paper adopts a directional splitting strategy to separate the convection in the configuration space and the velocity space. This can avoid the huge computational burden brought by directly dealing with six - dimensional problems and allows the use of one - dimensional update steps to gradually solve multi - dimensional problems.
### Main contributions
1. **High - order accuracy**: By introducing the splitting technique and the Active Flux method, high - order convergence can be achieved in time and space.
2. **Low energy consumption**: The Active Flux method can reduce numerical dissipation, thereby reducing the non - physical heating of the system.
3. **Efficient computation**: The splitting technique makes it only necessary to deal with one - dimensional problems at each step, greatly reducing the computational complexity and making large - scale simulations possible.
### Mathematical formulas
The Vlasov equation describes the time evolution of collisionless plasmas:
\[
\frac{\partial f_s}{\partial t}+\nabla_x\cdot (v f_s)+\nabla_v\cdot\left(\frac{q_s}{m_s}(E + v\times B)f_s\right)=0
\]
where \(f_s = f_s(x, v, t)\) is the probability density function of species \(s\), depending on time \(t\), macroscopic position \(x\in\Omega\subset\mathbb{R}^3\) and particle velocity \(v\in\mathbb{R}^3\). \(B = B(x, t)\in\mathbb{R}^3\) is the magnetic field, \(E = E(x, t)\in\mathbb{R}^3\) is the electric field, and \(q_s\) and \(m_s\) are the charge and mass of species \(s\), respectively.
In the electrostatic limit, the Maxwell equations are simplified to:
\[
\nabla^2\phi=-\frac{\rho}{\epsilon_0},\quad E = -\nabla\phi
\]
where \(\rho\) is the total charge density, calculated from the moments of the distribution function \(f_s\):
\[
\rho:=\sum_s q_s n_s,\quad j:=\sum_s q_s n_s