A quasilinear operator for wave–particle interactions in the electron cyclotron range is derived using the cold plasma dispersion relation. The finite width of the beam implies a broadening of the resonance layer in the velocity space which allows the use of a numerically efficient treatment of the operator. The specific case of a pure heating injected at the outer mid-plane is then treated. It allows to treat the beam propagation and absorption easily, simplifying the implementation of the source term in the gyrokinetic code. This specific case is implemented in the gyrokinetic particle-in-cell code ORB5 using a Langevin approach. The derivation and verification of the source operator are presented.

No abstract

While accurate simulations of dense gas flows far from the equilibrium can be achieved by direct simulation adapted to the Enskog equation, the significant computational demand required for collisions appears as a major constraint. In order to cope with that, an efficient yet accurate solution algorithm based on the Fokker-Planck approximation of the Enskog equation is devised in this paper; the approximation is very much associated with the Fokker-Planck model derived from the Boltzmann equation by Jenny et al. [“A solution algorithm for the fluid dynamic equations based on a stochastic model for molecular motion,” J. Comput. Phys. 229, 1077–1098 (2010)] and Gorji et al. [“Fokker–Planck model for computational studies of monatomic rarefied gas flows,” J. Fluid Mech. 680, 574–601 (2011)]. The idea behind these Fokker-Planck descriptions is to project the dynamics of discrete collisions implied by the molecular encounters into a set of continuous Markovian processes subject to the drift and diffusion. Thereby, the evolution of particles representing the governing stochastic process becomes independent from each other and thus very efficient numerical schemes can be constructed. By close inspection of the Enskog operator, it is observed that the dense gas effects contribute further to the advection of molecular quantities. That motivates a modelling approach where the dense gas corrections can be cast in the extra advection of particles. Therefore, the corresponding Fokker-Planck approximation is derived such that the evolution in the physical space accounts for the dense effects present in the pressure, stress tensor, and heat fluxes. Hence the consistency between the devised Fokker-Planck approximation and the Enskog operator is shown for the velocity moments up to the heat fluxes. For validation studies, a homogeneous gas inside a box besides Fourier, Couette, and lid-driven cavity flow setups is considered. The results based on the Fokker-Planck model are compared with respect to benchmark simulations, where good agreement is found for the flow field along with the transport properties.

Maximum-Entropy Distributions offer an attractive family of probability densities suitable for moment closure problems. Yet finding the Lagrange multipliers which parametrize these distributions, turns out to be a computational bottleneck for practical closure settings. Motivated by recent success of Gaussian processes, we investigate the suitability of Gaussian priors to approximate the Lagrange multipliers as a map of a given set of moments. Examining various kernel functions, the hyperparameters are optimized by maximizing the log-likelihood. The performance of the devised data-driven Maximum-Entropy closure is studied for couple of test cases including relaxation of non-equilibrium distributions governed by Bhatnagar-Gross-Krook and Boltzmann kinetic equations.

The kinetic theory of rarefied gases and numerical schemes based on the Boltzmann equation, have evolved to the cornerstone of non-equilibrium gas dynamics. However, their counterparts in the dense regime remain rather exotic for practical non-continuum scenarios. This problem is partly due to the fact that long-range interactions arising from the attractive tail of molecular potentials, lead to a computationally demanding Vlasov integral. This study focuses on numerical remedies for efficient stochastic particle simulations based on the Enskog-Vlasov kinetic equation. In particular, we devise a Poisson type elliptic equation which governs the underlying long-range interactions. The idea comes through fitting a Green function to the molecular potential, and hence deriving an elliptic equation for the associated fundamental solution. Through this transformation of the Vlasov integral, efficient Poisson type solvers can be readily employed in order to compute the mean field forces. Besides the technical aspects of different numerical schemes for treatment of the Vlasov integral, simulation results for evaporation of a liquid slab into the vacuum are presented. It is shown that the proposed formulation leads to accurate predictions with a reasonable computational cost.