Essentially Entropic Lattice Boltzmann Model

Toward a Mesoscopic Modeling Approach of Magnetohydrodynamic Blood Flow in Pathological Vessels: A Comprehensive Review

The investigation of magnetohydrodynamic (MHD) blood flow within configurations that are pertinent to the human anatomy holds significant importance in the realm of scientific inquiry because of its practical implications within the medical field. This article presents an exhaustive appraisal of the diverse applications of magnetohydrodynamics and their computational modeling in biological contexts. These applications are classified into two categories: simple flow and pulsatile flow. An alternative approach of traditional CFD methods called Lattice Boltzmann Method (LBM), a mesoscopic method based on kinetic theory, is introduced to solve complex problems, such as hemodynamics. The results show that the flow velocity reduces considerably by increasing the magnetic field intensity, and the flow separation area is minimized by the increase of magnetic field strength. The LBM with BGK collision model has shown good results in terms of precision. Finally, this literature review has revealed a number of potential avenues for further research. Suggestions for future works are proposed accordingly.

Essentially entropic lattice Boltzmann model: Theory and simulations

We present a detailed description of the essentially entropic lattice Boltzmann model. The entropic lattice Boltzmann model guarantees unconditional numerical stability by iteratively solving the nonlinear entropy evolution equation. In this paper we explain the construction of closed-form analytic solutions to this equation. We demonstrate that near equilibrium this analytic solution reduces to the standard lattice Boltzmann model. We consider a few test cases to show that the analytic solution does not exhibit any significant deviation from the iterative solution. We also extend the analytical solution for the Ellipsoidal Statistical (ES)-Bhatnagar-Gross-Krook model to remove the limitation on the Prandtl number for heat transfer problems. The simplicity of the analytic solution removes the computational overhead and algorithmic complexity associated with the entropic lattice Boltzmann models.

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.

Asymptotic method for entropic multiple relaxation time model in lattice Boltzmann method

To improve the numerical stability of the lattice Boltzmann method, Karlin et al. [Phys. Rev. E 90, 031302(R) (2014)10.1103/PhysRevE.90.031302] proposed the entropic multiple relaxation time (EMRT) collision model. The idea behind EMRT is to construct an optimal postcollision state by maximizing its local entropy value. The critical step of the EMRT model is to solve the entropy maximization problem under certain constraints, which is often computationally expensive and even not feasible. In this paper, we propose to employ perturbation theory and obtain an asymptotic solution to the maximum entropy state. With mathematical analysis of particular cases under relaxed constraints, we obtain the unperturbed form of the original problem and derive the asymptotic solution. We show that the asymptotic solution well approximates the optimal states; thus, our approach provides an efficient way to solve the constrained maximum entropy problem in the EMRT model. Also, we use the same idea of the EMRT model for the initial condition of the distribution function and propose to leave the entropy function to determine the missing information at the initial nodes. Finally, we numerically verify that the simulation results of the EMRT model obtained via the perturbation theory agree well with the exact solution to the Taylor-Green vortex problem. Furthermore, we also demonstrate that the EMRT model exhibits excellent stability performance for under-resolved simulations in the doubly periodic shear layer flow problem.

Discrete-velocity Boltzmann model: Regularization and linear stability

A discrete-velocity Boltzmann model for a nine-velocity lattice is considered. Compared to the conventional lattice Boltzmann (LB) schemes the collisions for the model are defined explicitly. Space and time discretization of the model is based on the collide and stream method; in addition, the regularization of the collision term is proposed. It is demonstrated that the regularized model can be represented as a two-relaxation-time LB model of a special type. The scheme is compared to the Onsager regularized (a specific filtered Galilean invariant model) and recursively regularized LB equations in terms of stability and dissipation properties, and linear stability analysis is performed. Several numerical experiments are carried out: double shear layer, lid-driven cavity flow, and propagation of acoustic and shear waves are considered for different grid resolutions, Mach and Reynolds numbers. It is shown that free parameters in the model corresponding to collision cross sections can be adjusted in such a way that the dissipation and stability properties can be controlled.

Onsager-regularized lattice Boltzmann method: A nonequilibrium thermodynamics-based regularized lattice Boltzmann method

The regularized class of lattice Boltzmann methods (LBMs) leverage the potency of the standard lattice-Bhatnagar-Gross-Krook method by filtering out spurious nonhydrodynamic moments from the moment space; this is achieved through evaluating regularized populations via a multiscale or a Hermite polynomial expansion approach. In this paper, we propose an alternative approach for evaluating the lattice populations. This approach is based on a kinetic theory that is consistent with nonequilibrium thermodynamics and obeys the Onsager-reciprocity principle. The proposed method is verified and validated for a number of canonical problems such as the athermal shock tube, the double periodic shear layer, the lid driven cavity, flow past square cylinder, and Poiseuille flow at nonvanishing Knudsen numbers. Additionally, the proposed method is compared to existing regularized LBM schemes and is shown to yield significant improvement in the stability and accuracy of the simulations.

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.

Estimating the herd immunity threshold by accounting for the hidden asymptomatics using a COVID-19 specific model

A quantitative COVID-19 model that incorporates hidden asymptomatic patients is developed, and an analytic solution in parametric form is given. The model incorporates the impact of lock-down and resulting spatial migration of population due to announcement of lock-down. A method is presented for estimating the model parameters from real-world data, and it is shown that the various phases in the observed epidemiological data are captured well. It is shown that increase of infections slows down and herd immunity is achieved when active symptomatic patients are 10-25% of the population for the four countries we studied. Finally, a method for estimating the number of asymptomatic patients, who have been the key hidden link in the spread of the infections, is presented.

Linear stability and isotropy properties of athermal regularized lattice Boltzmann methods

The present work proposes a general methodology to study stability and isotropy properties of lattice Boltzmann (LB) schemes. As a first investigation, such a methodology is applied to better understand these properties in the context of regularized approaches. To this extent, linear stability analyses of two-dimensional models are proposed: the standard Bhatnagar-Gross-Krook collision model, the original precollision regularization, and the recursive regularized model, where off-equilibrium distributions are partially computed thanks to a recursive formula. A systematic identification of the physical content carried by each LB mode is done by analyzing the eigenvectors of the linear systems. Stability results are then numerically confirmed by performing simulations of shear and acoustic waves. This work allows drawing fair conclusions on the stability properties of each model. In particular, among the aforementioned models, recursive regularization turns out to be the most stable one for the D2Q9 lattice, especially in the zero-viscosity limit. Two major properties shared by every regularized model are highlighted: (1) a mode filtering property and (2) an incorrect, and broadly anisotropic, dissipation rate of the modes carrying physical waves in under-resolved conditions. The first property is the main source of increased stability, especially for the recursive regularization. It is a direct consequence of the reconstruction of off-equilibrium populations before each collision process, decreasing the rank of the system of discrete equations. The second property seems to be related to numerical errors directly induced by the equilibration of high-order moments. In such a case, this property is likely to occur with any collision model that follows such a stabilization methodology.

Enhanced multi-relaxation-time lattice Boltzmann model by entropic stabilizers

The difficulty of choice of relaxation rates in multi-relaxation-time lattice Boltzmann model (MRT-LBM) is surmounted by solution of least-square problem of entropic stabilizer equations. Relaxation rates in the enhanced MRT-LBM are evolving with time rather than remain constants. To derive entropic stabilizer equations, nonequilibrium population is split into different modes in terms of column vectors in the inverse transform matrix. The entropic stabilizer equations are achieved by minimization of H-function. Different moment representations in MRT-LBM, such as Gram-Schmidt orthogonal moment, natural moment, and central moment, are tested for double periodic shear flow, shock tube problem, and lid-driven cavity flow, which demonstrates the potential of enhanced MRT-LBM.

Impact of collision models on the physical properties and the stability of lattice Boltzmann methods

The lattice Boltzmann method (LBM) is known to suffer from stability issues when the collision model relies on the BGK approximation, especially in the zero viscosity limit and for non-vanishing Mach numbers. To tackle this problem, two kinds of solutions were proposed in the literature. They consist in changing either the numerical discretization (finite-volume, finite-difference, spectral-element, etc.) of the discrete velocity Boltzmann equation (DVBE), or the collision model. In this work, the latter solution is investigated in detail. More precisely, we propose a comprehensive comparison of (static relaxation time based) collision models, in terms of stability, and with preliminary results on their accuracy, for the simulation of isothermal high-Reynolds number flows in the (weakly) compressible regime. It starts by investigating the possible impact of collision models on the macroscopic behaviour of stream-and-collide based D2Q9-LBMs, which clarifies the exact physical properties of collision models on LBMs. It is followed by extensive linear and numerical stability analyses, supplemented with an accuracy study based on the transport of vortical structures over long distances. In order to draw conclusions as generally as possible, the most common moment spaces (raw, central, Hermite, central Hermite and cumulant), as well as regularized approaches, are considered for the comparative studies. LBMs based on dynamic collision mechanisms (entropic collision, subgrid-scale models, explicit filtering, etc.) are also briefly discussed. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

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.

Lattice Boltzmann model for weakly compressible flows

We present an energy conserving lattice Boltzmann model based on a crystallographic lattice for simulation of weakly compressible flows. The theoretical requirements and the methodology to construct such a model are discussed. We demonstrate that the model recovers the isentropic sound speed in addition to the effects of viscous heating and heat flux dynamics. Several test cases for acoustics and thermal and thermoacoustic flows are simulated to show the accuracy of the proposed model.

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.

Comprehensive comparison of collision models in the lattice Boltzmann framework: Theoretical investigations

Over the last decades, several types of collision models have been proposed to extend the validity domain of the lattice Boltzmann method (LBM), each of them being introduced in its own formalism. This article proposes a formalism that describes all these methods within a common mathematical framework, and in this way allows us to draw direct links between them. Here, the focus is put on single and multirelaxation time collision models in either their raw moment, central moment, cumulant, or regularized form. In parallel with that, several bases (nonorthogonal, orthogonal, Hermite) are considered for the polynomial expansion of populations. General relationships between moments are first derived to understand how moment spaces are related to each other. In addition, a review of collision models further sheds light on collision models that can be rewritten in a linear matrix form. More quantitative mathematical studies are then carried out by comparing explicit expressions for the post-collision populations. Thanks to this, it is possible to deduce the impact of both the polynomial basis (raw, Hermite, central, central Hermite, cumulant) and the inclusion of regularization steps on isothermal LBMs. Extensive results are provided for the D1Q3, D2Q9, and D3Q27 lattices, the latter being further extended to the D3Q19 velocity discretization. Links with the most common two and multirelaxation time collision models are also provided for the sake of completeness. This work ends by emphasizing the importance of an accurate representation of the equilibrium state, independently of the choice of moment space. As an addition to the theoretical purpose of this article, general instructions are provided to help the reader with the implementation of the most complicated collision models.

Stability of the lattice kinetic scheme and choice of the free relaxation parameter

The lattice kinetic scheme (LKS), a modified version of the classical single relaxation time (SRT) lattice Boltzmann method, was initially developed as a suitable numerical approach for non-Newtonian flow simulations and a way to reduce memory consumption of the original SRT approach. The better performances observed for non-Newtonian flows are mainly due to the additional degree of freedom allowing an independent control over the relaxation of higher-order moments, independently from the fluid viscosity. Although widely applied to fluid flow simulations, no theoretical analysis of LKS has been performed. The present work focuses on a systematic von Neumann analysis of the linearized collision operator. Thanks to this analysis, the effects of the modified collision operator on the stability domain and spectral behavior of the scheme are clarified. Results obtained in this study show that correct choices of the "second relaxation coefficient" lead, to a certain extent, to a more consistent dispersion and dissipation for large values of the first relaxation coefficient. Furthermore, appropriate values of this parameter can lead to a larger linear stability domain. At moderate and low values of viscosity, larger values of the free parameter are observed to increase dissipation of kinetic modes, while leaving the acoustic modes untouched and having a less pronounced effect on the convective mode. This increased dissipation leads in general to less pronounced sources of nonlinear instability, thus improving the stability of the LKS.

Comparison of the lattice-Boltzmann model with the finite-difference time-domain method for electrodynamics

A three-dimensional lattice-Boltzmann model (LBM) for the simulation of the Maxwell equations is presented. The inclusion of media follows an extension of a special limit described in the literature which is applicable to this LBM and does not harm the stability of simulations. The focus of the present study lies on the properties of numerical accuracy and stability of the LBM in comparison to the standard finite-difference time-domain (FDTD) method based on Yee's method. Typical examples, often investigated in the context of numerical simulations, are considered. These include the propagation of electrodynamic (EM) fields in one- and three-dimensional systems. Results of this simulations are compared to the ones of their theoretical predictions. Further on, long-time simulations are done in systems with periodic boundary conditions to check if the total energy is conserved. To investigate the effect of the numeric impedance, the propagation of an EM pulse is monitored spatially and temporarily in a two-dimensional system. The simulation results indicate, in contrast to the one obtained from the FDTD method, that the presented LBM does fulfill the expected energy conservation and is not effected by the numerical impedance. This LBM therefore represents a valuable alternative for the simulation of EM problems like long-time simulations by avoiding intrinsic properties the FDTD method suffers from.

Particles on Demand for Kinetic Theory

A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes restrictions on flow speed and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. In the new kinetic theory, discrete particles are determined by a rigorous limit process which avoids ad hoc assumptions about their velocities. Classical benchmarks for incompressible and compressible flows demonstrate that the proposed discrete-particles kinetic theory opens up an unprecedented wide domain of applications for computational fluid dynamics.

Discrete Boltzmann trans-scale modeling of high-speed compressible flows

We present a general framework for constructing trans-scale discrete Boltzmann models (DBMs) for high-speed compressible flows ranging from continuum to transition regime. This is achieved by designing a higher-order discrete equilibrium distribution function that satisfies additional nonhydrodynamic kinetic moments. To characterize the thermodynamic nonequilibrium (TNE) effects and estimate the condition under which the DBMs at various levels should be used, two measures are presented: (i) the relative TNE strength, describing the relative strength of the (N+1)th order TNE effects to the Nth order one; (ii) the TNE discrepancy between DBM simulation and relevant theoretical analysis. Whether or not the higher-order TNE effects should be taken into account in the modeling and which level of DBM should be adopted is best described by the relative TNE intensity and/or the discrepancy rather than by the value of the Knudsen number. As a model example, a two-dimensional DBM with 26 discrete velocities at Burnett level is formulated, verified, and validated.