The Solvability of Mixed Value Problem for the First and Second Approximations of One-Dimensional Nonlinear System of Moment Equations with Microscopic Boundary Conditions

The paper gives a derivation of a new one-dimensional non-stationary nonlinear system of moment equations, that depend on the flight velocity and the surface temperature of an aircraft. Maxwell microscopic condition is approximated for the distribution function on moving boundary, when one fraction of molecules reflected from the surface specular and another fraction diffusely with Maxwell distribution. Moreover, macroscopic boundary conditions for the moment system of equations depend on evenness or oddness of approximation $${f}_{k}(t,x,c)$$ f k ( t , x , c ) , where $${f}_{k}(t,x,c)$$ f k ( t , x , c ) is partial expansion sum of the molecules distribution function over eigenfunctions of linearized collision operator around local Maxwell distribution. The formulation of initial and boundary value problem for the system of moment equations in the first and second approximations is described. Existence and uniqueness of the solution for the above-mentioned problem using macroscopic boundary conditions in the space of functions $$C\left(\left[0,T\right];{L}^{2}\left[-a,a\right]\right)$$ C 0 , T ; L 2 - a , a are proved.

From kinetic to fluid models of liquid crystals by the moment method

<p style='text-indent:20px;'>This paper deals with the convergence of the Doi-Navier-Stokes model of liquid crystals to the Ericksen-Leslie model in the limit of the Deborah number tending to zero. While the literature has investigated this problem by means of the Hilbert expansion method, we develop the moment method, i.e. a method that exploits conservation relations obeyed by the collision operator. These are non-classical conservation relations which are associated with a new concept, that of Generalized Collision Invariant (GCI). In this paper, we develop the GCI concept and relate it to geometrical and analytical structures of the collision operator. Then, the derivation of the limit model using the GCI is performed in an arbitrary number of spatial dimensions and with non-constant and non-uniform polymer density. This non-uniformity generates new terms in the Ericksen-Leslie model.</p>

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N-body simulations as well as observations of gravitational systems.

Kinetic simulations of collision-less plasmas in open magnetic geometries†

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (1) Scrape-Of-Layer (SOL) and divertor regions in toroidal confinement fusion devices , (2) linear divertor simulators, (3) plasma-based thrusters and (4) magnetic mirrors. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic Particle-In-Cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (1) the effect of non-uniform magnetic fields, (2) finite fully-absorbing boundaries for the particles and (3) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. Results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, inclusion of a Monte-Carlo based Fokker-Planck collision operator is discussed in the context of future work.

Analytical Solution for Boltzmann Collision Operator for the1-D Diffusion equation

In this paper, we have presented the analytical solution of the collision operator for the Boltzmann equation of onedimensional diffusion equation using the analytical solution of the one-dimensional Navier Stoke diffusion equation. Keywords: Boltzmann equation; analytical collision operator; one-dimensional diffusion equation.

Development of advanced linearized gyrokinetic collision operators using a moment approach

The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge et al., J. Plasma Phys., vol. 85, 2019, 905850604) and of the widely used linearized GK Sugama collision operator (Sugama et al., Phys. Plasmas, vol. 16, 2009, 112503) is reported. An approach based on a Hermite–Laguerre moment expansion of the perturbed gyrocentre distribution function is used, referred to as gyromoment expansion. This approach allows the considering of arbitrary perpendicular wavenumber and expressing the two linearized GK operators as a linear combination of gyromoments where the expansion coefficients are given by closed analytical expressions that depend on the perpendicular wavenumber and on the temperature and mass ratios of the colliding species. The drift-kinetic (DK) limits of the GK linearized Coulomb and Sugama operators are also obtained. Comparisons between the gyromoment approach and the DK Coulomb and GK Sugama operators in the continuum GK code GENE are reported, focusing on the ion-temperature-gradient instability and zonal flow damping, finding an excellent agreement. It is confirmed that stronger collisional damping of the zonal flow residual by the Sugama GK model compared with the GK linearized Coulomb (Pan et al., Phys. Plasmas, vol. 27, 2020, 042307) persists at higher collisionality. Finally, we show that the numerical efficiency of the gyromoment approach increases with collisionality, a desired property for boundary plasma applications.

A unified lattice Boltzmann model and application to multiphase flows

In this work, we develop a unified lattice Boltzmann model (ULBM) framework that can seamlessly integrate the widely used lattice Boltzmann collision operators, including the Bhatnagar–Gross–Krook or single-relation-time, multiple-relaxation-time, central-moment or cascaded lattice Boltzmann method and multiple entropic operators (KBC). Such a framework clarifies the relations among the existing collision operators and greatly facilitates model comparison and development as well as coding. Importantly, any LB model or treatment constructed for a specific collision operator could be easily adopted by other operators. We demonstrate the flexibility and power of the ULBM framework through three multiphase flow problems: the rheology of an emulsion, splashing of a droplet on a liquid film and dynamics of pool boiling. Further exploration of ULBM for a wide variety of phenomena would be both realistic and beneficial, making the LBM more accessible to non-specialists. This article is part of the theme issue ‘Progress in mesoscale methods for fluid dynamics simulation’.