A unified gas kinetic scheme for transport and collision effects in plasma

In this study, the Vlasov-Poisson equation with or without collision term for plasma is solved by the unified gas kinetic scheme (UGKS). The Vlasov equation is a differential equation describing time evolution of the distribution function of plasma consisting of charged particles with long-range interaction. The distribution function is discretized in discrete particle velocity space. After the Vlasov equation is integrated in finite volumes of physical space, the numerical flux across a cell interface and source term for particle acceleration are computed to update the distribution function at next time step. The flux is decided by Riemann problem and variation of distribution function in discrete particle velocity space is evaluated with central difference method. A electron-ion collision model is introduced in the Vlasov equation. This finite volume method for the UGKS couples the free transport and long-range interaction between particles. The electric field induced by charged particles is controlled by the Poisson's equation. In this paper, the Poisson's equation is solved using the Green's function for two dimensional plasma system subjected to the symmetry or periodic boundary conditions. Two numerical tests of the linear Landau damping and the Gaussian beam are carried out to validate the proposed method. The linear electron plasma wave damping is simulated based on electron-ion collision operator. Compared with previous methods, it is shown that the current method is able to obtain accurate results of the Vlasov-Poisson equation with a time step much larger than the particle collision time. Highly non-equilibrium and rarefied plasma flows, such as electron flows driven by electromagnetic field, can be simulated easily.Comment: 33 pages, 13 figure

A simplified finite volume lattice Boltzmann method for simulations of fluid flows from laminar to turbulent regime, Part I: Numerical framework and its application to laminar flow simulation

In this paper, a finite volume lattice Boltzmann method (FVLBM) based on cell-center unstructured girds is presented and full studied to simulate the incompressible laminar flows, which is simple modified from the cell-vertex unstructured girds FVLBM proposed by Stiebler et al. [Computers & Fluids, 2006, 35(8): 814-819]. Compared with other complex flux reconstruct methods, the computational cost of present scheme is little and can achieve second-order spatial accuracy, the temporal accuracy is adjustable depending on the temporal discretization methods. Different boundary conditions are illustrated and easy implement to the complex geometries. Four cases are testified to validate the present method, including one plate driven Couette flow for accuracy test, flow in the square cavity, flow over the single circular cylinder and more complex double circular cylinders. Numerical experiments show that the present scheme can use relatively few grid cells to simulate relatively higher Reynolds number flow, steady and unsteady flows, demonstrate the good capability of the present method.Comment: 56 Pages; 36 figure

High-order gas-kinetic scheme with TENO class reconstruction for the Euler and Navier-Stokes equations

The high-order gas-kinetic scheme(HGKS) with WENO spatial reconstruction method has been extensively validated through many numerical experiments, demonstrating its superior accuracy efficiency, and robustness. Compared with WENO class schemes, TENO class schemes exhibit significantly improved robustness, low numerical dissipation and sharp discontinuity capturing. In this paper, two kinds of fifth-order HGKS with TENO class schemes are designed. One involves replacing WENO5 scheme with the TENO5 scheme in the conventional WENO5-GKS. WENO and TENO schemes only provide the non-equilibrium state values at the cell interface. The slopes of the non-equilibrium state along with the equilibrium values and slopes, are obtained by additional linear reconstruction. Another kind of TENO5-D GKS is similar to WENO5-AO GKS. Following a strong scale-separation procedure, a tailored novel ENO-like stencil selection strategy is proposed such that the high-order accuracy is restored in smooth regions by selecting the candidate reconstruction on the large stencil while the ENO property is enforced near discontinuities by adopting the candidate reconstruction from smooth small stencils. The such TENO schemes are TENO-AA and TENO-D scheme. The HGKS scheme based on WENO-AO or TENO-D reconstruction take advantage of the large stencil to provide point values and slopes of the non-equilibrium state. By dynamically merging the reconstructed non-equilibrium slopes, extra reconstruction of the equilibrium state at the beginning of each time step can be avoided. The simplified schemes have better robustness and efficiency than the conventional WENO5-GKS or TENO5-GKS. TENO-D GKS is also as easy to develop as WENO-AO GKS to high-order finite volume method for unstructured mesh.Comment: arXiv admin note: text overlap with arXiv:2304.05572; text overlap with arXiv:1905.08489 by other author

A global adaptive velocity space for general discrete velocity framework in predictions of rarefied and multi-scale flows

The rarefied flow and multi-scale flow are crucial for the aerodynamic design of spacecraft, ultra-low orbital vehicles and plumes. By introducing a discrete velocity space, the Boltzmann method, such as the discrete velocity method and unified methods, can capture complex and non-equilibrium velocity distribution functions (VDFs) and describe flow behaviors exactly. However, the extremely steep slope and high concentration of the gas VDFs in a local particle velocity space make it very difficult for the Boltzmann method with structured velocity space to describe high speed flow. Therefore, the adaptive velocity space (AVS) is required for the Boltzmann solvers to be practical in complex rarefied flow and multi-scale flow. This paper makes two improvements to the AVS approach, which is then incorporated into a general discrete velocity framework, such as the unified gas-kinetic scheme. Firstly, a global velocity mesh is used to prevent the interpolation of the VDFs at the physical interface during the calculation of the microscopic fluxes, maintaining the program's high level of parallelism. Secondly, rather than utilizing costly interpolation, the VDFs on a new velocity space were reconstruction using the ``consanguinity" relationship. In other words, a split child node's VDF is the same as its parent's VDF, and a merged parent's VDF is the average of its children's VDFs. Additionally, the discrete deviation of the equilibrium distribution functions is employed to maintain the proposed method's conservation. Moreover, an appropriate set of adaptive parameters is established to enhance the automation of the proposed method. Finally, a number of numerical tests are carried out to validate the proposed method