Particles on Demand for Kinetic Theory

Probing double-distribution-function models in discrete-velocity Boltzmann methods for highly compressible flows: Particles-on-demand realization

The double distribution function approach is an efficient route toward an extension of kinetic solvers to compressible flows. With a number of realizations available, an overview and comparative study in the context of high-speed compressible flows is presented. We discuss the different variants of the energy partition, analyses of hydrodynamic limits, and a numerical study of accuracy and performance with the particles on demand realization. Out of three considered energy partition strategies, it is shown that the nontranslational energy split requires a higher-order quadrature for proper recovery of the Navier-Stokes-Fourier equations. The internal energy split, on the other hand, while recovering the correct hydrodynamic limit with fourth-order quadrature, comes with a nonlocal-both in space and time-source term that contributes to higher computational cost and memory overhead. Based on our analysis, the total energy split demonstrates the optimal overall performance.

Particles on demand method: Theoretical analysis, simplification techniques, and model extensions

The particles on demand method [Phys. Rev. Lett. 121, 130602 (2018)0031-900710.1103/PhysRevLett.121.130602] was recently formulated with a conservative finite-volume discretization and validated against challenging benchmarks. In this work, we focus on the properties of the reference frame transformation and its implications on the accuracy of the model. Based on these considerations, we propose strategies that simplify the scheme and generalize it to include a tunable Prandtl number via quasi-equilibrium relaxation. Finally, we adapt concepts from the multiscale semi-Lagrangian lattice Boltzmann formulation to the proposed framework, further improving the potential and the operating range of the kinetic model. Numerical simulations of high Mach compressible flows demonstrate excellent accuracy and stability of the model over a wide range of conditions.

Discrete unified gas kinetic scheme for continuum compressible flows

In this paper, a discrete unified gas kinetic scheme (DUGKS) is proposed for continuum compressible gas flows based on the total energy kinetic model [Guo et al., Phys. Rev. E 75, 036704 (2007)1539-375510.1103/PhysRevE.75.036704]. The proposed DUGKS can be viewed as a special finite-volume lattice Boltzmann method for the compressible Navier-Stokes equations in the double distribution function formulation, in which the mass and momentum transport are described by the kinetic equation for a density distribution function (g), and the energy transport is described by the other one for an energy distribution function (h). To recover the full compressible Navier-Stokes equations exactly, the corresponding equilibrium distribution functions g^{eq} and h^{eq} are expanded as Hermite polynomials up to third and second orders, respectively. The velocity spaces for the kinetic equations are discretized according to the seventh and fifth Gauss-Hermite quadratures. Consequently, the computational efficiency of the present DUGKS can be much improved in comparison with previous versions using more discrete velocities required by the ninth Gauss-Hermite quadrature.

Particles on demand for flows with strong discontinuities

Particles-on-demand formulation of kinetic theory [B. Dorschner, F. Bösch and I. V. Karlin, Phys. Rev. Lett. 121, 130602 (2018)0031-900710.1103/PhysRevLett.121.130602] is used to simulate a variety of compressible flows with strong discontinuities in density, pressure, and velocity. Two modifications are applied to the original formulation of the particles-on-demand method. First, a regularization by Grad's projection of particles populations is combined with the reference frame transformations in order to enhance stability and accuracy. Second, a finite-volume scheme is implemented which allows tight control of mass, momentum, and energy conservation. The proposed model is validated with an array of challenging one- and two-dimensional benchmarks of compressible flows, including hypersonic and near-vacuum situations, Richtmyer-Meshkov instability, double Mach reflection, and astrophysical jet. Excellent performance of the modified particles-on-demand method is demonstrated beyond the limitations of other lattice Boltzmann-like approaches to compressible flows.

Multiscale analysis of the particles on demand kinetic model

We present a thorough investigation of the particles on demand kinetic model. After a brief introduction of the method, an appropriate multiscale analysis is carried out to derive the hydrodynamic limit of the model. In this analysis, the effect of the time-space dependent comoving reference frames are taken into account. This could be regarded as a generalization of the conventional Chapman-Enskog analysis applied to the lattice Boltzmann models which feature global constant reference frames. Further simulations of target benchmarks provide numerical evidence confirming the theoretical predictions.

A multiple-relaxation-time collision model by Hermite expansion

The Bhatnagar-Gross-Krook (BGK) single-relaxation-time collision model for the Boltzmann equation serves as the foundation of the lattice BGK (LBGK) method developed in recent years. The description of the collision as a uniform relaxation process of the distribution function towards its equilibrium is, in many scenarios, simplistic. Based on a previous series of papers, we present a collision model formulated as independent relaxations of the irreducible components of the Hermite coefficients in the reference frame moving with the fluid. These components, corresponding to the irreducible representation of the rotation group, are the minimum tensor components that can be separately relaxed without violating rotation symmetry. For the 2nd, 3rd and 4th moments, respectively, two, two and three independent relaxation rates can exist, giving rise to the shear and bulk viscosity, thermal diffusivity and some high-order relaxation process not explicitly manifested in the Navier-Stokes-Fourier equations. Using the binomial transform, the Hermite coefficients are evaluated in the absolute frame to avoid the numerical dissipation introduced by interpolation. Extensive numerical verification is also provided. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

High-order semi-Lagrangian kinetic scheme for compressible turbulence

Turbulent compressible flows are traditionally simulated using explicit time integrators applied to discretized versions of the Navier-Stokes equations. However, the associated Courant-Friedrichs-Lewy condition severely restricts the maximum time-step size. Exploiting the Lagrangian nature of the Boltzmann equation's material derivative, we now introduce a feasible three-dimensional semi-Lagrangian lattice Boltzmann method (SLLBM), which circumvents this restriction. While many lattice Boltzmann methods for compressible flows were restricted to two dimensions due to the enormous number of discrete velocities in three dimensions, the SLLBM uses only 45 discrete velocities. Based on compressible Taylor-Green vortex simulations we show that the new method accurately captures shocks or shocklets as well as turbulence in 3D without utilizing additional filtering or stabilizing techniques other than the filtering introduced by the interpolation, even when the time-step sizes are up to two orders of magnitude larger compared to simulations in the literature. Our new method therefore enables researchers to study compressible turbulent flows by a fully explicit scheme, whose range of admissible time-step sizes is dictated by physics rather than spatial discretization.

Multiscale semi-Lagrangian lattice Boltzmann method

We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles' velocity space. Different velocity sets of lower and higher order are consistently and efficiently coupled, allowing us to use the higher-order model only when and where needed. This includes regions of high Mach or high Knudsen numbers. The coupling procedure of discrete velocity sets consists of either a projection of the higher-order populations onto the lower-order lattice or lifting of the lower-order populations to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme is formulated for both a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is validated with an advection of an athermal vortex and in a jet flow setup. The performance of the proposed scheme is further investigated in the shock structure problem and a high-Knudsen-number Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.

Mesoscopic kinetic approach for studying nonequilibrium hydrodynamic and thermodynamic effects of shock wave, contact discontinuity, and rarefaction wave in the unsteady shock tube

This paper presents a detailed description of a molecular velocity distribution-based mesoscopic kinetic approach that enables a better understanding of various nonequilibrium hydrodynamic and thermodynamic effects in shock waves, contact discontinuities, and rarefaction waves. This builds on the mesoscopic kinetic approach in a previous investigation into regular reflection shocks by further addressing the mesoscopic physical meaning of kinetic moments from the view of kinetics and the implications of the magnitude and sign of nonequilibrium kinetic moments. To deepen understanding of nonequilibrium effects, this work focuses on the one-dimensional unsteady shock tube problem, which contains the typical and essential features of the discontinuous flows, and has no interference of two-dimensional flow direction. The approach uses a lattice Boltzmann method to solve the flow field, and describes nonequilibrium effects through the nonequilibrium kinetic moments of molecular velocity distribution functions. The mechanism of nonequilibrium effect in discontinuous flows is further probed. This work develops the mesoscopic kinetic approach and clarifies the mesoscopic physics of shock waves, contact discontinuities, and rarefaction waves.

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.

Kinetic Simulations of Compressible Non-Ideal Fluids: From Supercritical Flows to Phase-Change and Exotic Behavior

We investigate a kinetic model for compressible non-ideal fluids. The model imposes the local thermodynamic pressure through a rescaling of the particle’s velocities, which accounts for both long- and short-range effects and hence full thermodynamic consistency. The model is fully Galilean invariant and treats mass, momentum, and energy as local conservation laws. The analysis and derivation of the hydrodynamic limit is followed by the assessment of accuracy and robustness through benchmark simulations ranging from the Joule–Thompson effect to a phase-change and high-speed flows. In particular, we show the direct simulation of the inversion line of a van der Waals gas followed by simulations of phase-change such as the one-dimensional evaporation of a saturated liquid, nucleate, and film boiling and eventually, we investigate the stability of a perturbed strong shock front in two different fluid mediums. In all of the cases, we find excellent agreement with the corresponding theoretical analysis and experimental correlations. We show that our model can operate in the entire phase diagram, including super- as well as sub-critical regimes and inherently captures phase-change phenomena.

Thermokinetic lattice Boltzmann model of nonideal fluids

We present a kinetic model for nonideal fluids, where the local thermodynamic pressure is imposed through appropriate rescaling of the particle's velocities, accounting for both long- and short-range effects and hence full thermodynamic consistency. The model features full Galilean invariance together with mass, momentum, and energy conservation and enables simulations ranging from subcritical to supercritical flows, which is illustrated on various benchmark flows such as anomalous shock waves or shock droplet interaction.

Semi-Lagrangian implicit Bhatnagar-Gross-Krook collision model for the finite-volume discrete Boltzmann method

In order to increase the accuracy of temporal solutions, reduce the computational cost of time marching, and improve the stability associated with collisions for the finite-volume discrete Boltzmann method, an advanced implicit Bhatnagar-Gross-Krook (BGK) collision model using a semi-Lagrangian approach is proposed in this paper. Unlike existing models, in which the implicit BGK collision is resolved either by a temporal extrapolation or by a variable transformation, the proposed model removes the implicitness by tracing the particle distribution functions (PDFs) back in time along their characteristic paths during the collision process. An interpolation scheme is needed to evaluate the PDFs at the traced-back locations. By using the first-order interpolation, the resulting model allows for the straightforward replacement of f_{α}^{eq,n+1} by f_{α}^{eq,n} no matter where it appears. After comparing the proposed model with the existing models under different numerical conditions (e.g., different flux schemes and time-marching schemes) and using the proposed model to successfully modify the variable transformation technique, three conclusions can be drawn. First, the proposed model can improve the accuracy by almost an order of magnitude. Second, it can slightly reduce the computational cost. Therefore, the proposed scheme improves accuracy without extra cost. Finally, the proposed model can significantly improve the Δt/τ limit compared to the temporal interpolation model while having the same Δt/τ limit as the variable transformation approach. The proposed scheme with a second-order interpolation is also developed and tested; however, that technique displays no advantage over the simple first-order interpolation approach. Both numerical and theoretical analyses are also provided to explain why the developed implicit scheme with simple first-order interpolation can outperform the same scheme with second-order interpolation, as well as the existing temporal extrapolation and variable transformation schemes.

Efficient supersonic flow simulations using lattice Boltzmann methods based on numerical equilibria

A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on a numerical equilibrium that has been extensively used by discrete velocity methods since the late 1990s. Here, it is extended to reproduce an arbitrary number of moments of the Maxwell-Boltzmann distribution. These extensions to the standard 5-constraint (mass, momentum and energy) approach lead to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high Reynolds-number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the supersonic flow past a NACA0012 aerofoil confirms the excellent behaviour of this model in a low-viscosity and supersonic regime. The flow past a sphere is further simulated to investigate the 3D behaviour of our model in the low-viscosity supersonic regime. The proposed model is shown to be substantially more efficient than the previous 5-moment D3Q343 DDF-LBM for both CPU and GPU architectures. It then opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Compressibility in lattice Boltzmann on standard stencils: effects of deviation from reference temperature

With growing interest in the simulation of compressible flows using the lattice Boltzmann (LB) method, a number of different approaches have been developed. These methods can be classified as pertaining to one of two major categories: (i) solvers relying on high-order stencils recovering the Navier-Stokes-Fourier equations, and (ii) approaches relying on classical first-neighbour stencils for the compressible Navier-Stokes equations coupled to an additional (LB-based or classical) solver for the energy balance equation. In most cases, the latter relies on a thermal Hermite expansion of the continuous equilibrium distribution function (EDF) to allow for compressibility. Even though recovering the correct equation of state at the Euler level, it has been observed that deviations of local flow temperature from the reference can result in instabilities and/or over-dissipation. The aim of the present study is to evaluate the stability domain of different EDFs, different collision models, with and without the correction terms for the third-order moments. The study is first based on a linear von Neumann analysis. The correction term for the space- and time-discretized equations is derived via a Chapman-Enskog analysis and further corroborated through spectral dispersion-dissipation curves. Finally, a number of numerical simulations are performed to illustrate the proposed theoretical study. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Semi-Lagrangian lattice Boltzmann method for compressible flows

This work thoroughly investigates a semi-Lagrangian lattice Boltzmann (SLLBM) solver for compressible flows. In contrast to other LBM for compressible flows, the vertices are organized in cells, and interpolation polynomials up to fourth order are used to attain the off-vertex distribution function values. Differing from the recently introduced Particles on Demand (PoD) method [Dorschner, Bösch, and Karlin, Phys. Rev. Lett. 121, 130602 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.130602], the method operates in a static, nonmoving reference frame. Yet the SLLBM in the present formulation grants supersonic flows and exhibits a high degree of Galilean invariance. The SLLBM solver allows for an independent time step size due to the integration along characteristics and for the use of unusual velocity sets, like the D2Q25, which is constructed by the roots of the fifth-order Hermite polynomial. The properties of the present model are shown in diverse example simulations of a two-dimensional Taylor-Green vortex, a Sod shock tube, a two-dimensional Riemann problem, and a shock-vortex interaction. It is shown that the cell-based interpolation and the use of Gauss-Lobatto-Chebyshev support points allow for spatially high-order solutions and minimize the mass loss caused by the interpolation. Transformed grids in the shock-vortex interaction show the general applicability to nonuniform grids.

Consistent vortex initialization for the athermal lattice Boltzmann method

A barotropic counterpart of the well-known convected vortex test case is rigorously derived from the Euler equations along with an athermal equation of state. Starting from a given velocity distribution corresponding to an intended flow recirculation, the athermal counterpart of the Euler equations are solved to obtain a consistent density field. The present initialization is assessed on a standard lattice Boltzmann solver based on the D2Q9 lattice. Compared to the usual isentropic initialization, a much lower spurious relaxation toward the targeted solution is observed, which is due to the spatial resolution rather than approximated macroscopic quantities. The amplitude of the spurious waves can be further reduced by including an off-equilibrium part in the initial distribution functions.

Extensive analysis of the lattice Boltzmann method on shifted stencils

Standard lattice Boltzmann methods (LBMs) are based on a symmetric discretization of the phase space, which amounts to study the evolution of particle distribution functions (PDFs) in a reference frame at rest. This choice induces a number of limitations when the simulated flow speed gets closer to the sound speed, such as velocity-dependent transport coefficients. The latter issue is usually referred to as a Galilean invariance defect. To restore the Galilean invariance of LBMs, it was proposed to study the evolution of PDFs in a comoving reference frame by relying on asymmetric shifted lattices [N. Frapolli, S. S. Chikatamarla, and I. V. Karlin, Phys. Rev. Lett. 117, 010604 (2016)].PRLTAO0031-900710.1103/PhysRevLett.117.010604 From the numerical viewpoint, this corresponds to overcoming the rather restrictive Courant-Friedrichs-Lewy conditions on standard LBMs and modeling compressible flows while keeping memory consumption and processing costs to a minimum (therefore using the standard first-neighbor stencils). In the present work systematic physical error evaluations and stability analyses are conducted for different discrete equilibrium distribution functions (EDFs) and collision models. Thanks to them, it is possible to (1) better understand the effect of this solution on both physics and stability, (2) assess its viability as a way to extend the validity range of LBMs, and (3) quantify the importance of the reference state as compared to other parameters such as the equilibrium state and equilibration path. The results clearly show that, in theory, the concept of shifted lattices allows the scheme to deal with arbitrarily high values of the nondimensional velocity. Furthermore, just like the zero-Mach flow for the standard stencils, it is observed that setting the shift velocity to the fluid velocity results in optimal physical and numerical properties. In addition, a detailed analysis of the obtained results shows that the properties of different collision models and EDFs remain unchanged under the shift of stencil. In other words, by introducing a velocity shift in the stencil, the optimal operating point, in terms of physics and numerics, will also be shifted by the same vector regardless of the EDF or collision model considered. Eventually, while limited to the D2Q9 stencil with the nine possible first-neighbor shifts, the present study and corresponding conclusions can be extended to other stencils and velocity shifts in a straightforward manner.

Pseudoentropic derivation of the regularized lattice Boltzmann method

The lattice Boltzmann method (LBM) facilitates efficient simulations of fluid turbulence based on advection and collision of local particle distribution functions. To ensure stable simulations on underresolved grids, the collision operator must prevent drastic deviations from local equilibrium. This can be achieved by various methods, such as the multirelaxation time, entropic, quasiequilibrium, regularized, and cumulant schemes. Complementing a part of a unified theoretical framework of these schemes, the present work presents a derivation of the regularized lattice Boltzmann method (RLBM), which follows a recently introduced entropic multirelaxation time LBM by Karlin, Bösch, and Chikatamarla (KBC). It is shown that both methods can be derived by locally maximizing a quadratic Taylor expansion of the entropy function. While KBC expands around the local equilibrium distribution, the RLBM is recovered by expanding entropy around a global equilibrium. Numerical tests were performed to elucidate the role of pseudoentropy maximization in these models. Simulations of a two-dimensional shear layer show that the RLBM successfully reproduces the largest eddies even on a 16×16 grid, while the conventional LBM becomes unstable for grid resolutions of 128×128 and lower. The RLBM suppresses spurious vortices more effectively than KBC. In contrast, simulations of the three-dimensional Taylor-Green and Kida vortices show that KBC performs better in resolving small scale vortices, outperforming the RLBM by a factor of 1.8 in terms of the effective Reynolds number.

Temperature-scaled collision process for the high-order lattice Boltzmann model

We postulate that the relaxations of the distribution function in the lattice Boltzmann model should be self-similar under temperature scaling. Based on this postulation, a multiple-relaxation-time collision model in the relative, temperature-scaled reference frame is devised with Hermite expansion. Resorting to the relation between the Hermite basis with the temperature-scaled relative velocity and the Hermite basis with the raw velocity, the relaxations in the temperature-scaled reference frame can be converted to those in the raw reference frame with some correction terms to eliminate the cross-talk effects among the relaxations of different orders. The highest-order nonequilibrium relative central moment is filtered due to the insufficient discrerization in the velocity space. The highest-order collision term can be recursively obtained from the lower-order collision terms. The improved performance is validated by the double shear layer flow, shock tube flow, and the Taylor-Green vortex flow.

Lattice Boltzmann model for compressible flows on standard lattices: Variable Prandtl number and adiabatic exponent

A lattice Boltzmann model for compressible flows on standard lattices is developed and analyzed. A consistent two-population thermal lattice Boltzmann is used which allows a variable Prandtl number and a variable adiabatic exponent, and appropriate correction terms are introduced into the kinetic equations to compensate for deviations in the hydrodynamic limit. Using the concept of a shifted lattice, the model is extended to supersonic flows involving shock waves, and the shock-vortex interaction problem is simulated to show the accuracy of the proposed model. Numerical results demonstrate that the proposed model is a viable candidate for compressible flow simulations.