Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow: Two-Dimensional Case

Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from $0$ to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax--Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.     

Lattice Boltzmann Simulations of Blood Flow: Non-Newtonian Rheology and Clotting Processes

The numerical simulation of thrombosis in stented aneurysms is an important issue to estimate the efficiency of a stent. In this paper, we consider a Lattice Boltzmann (LB) approach to bloodflow modeling and we implement a non-Newtonian correction in order to reproduce more realistic flow profiles. We obtain a good agreement between simulations and Casson’s model of blood rheology in a simple geometry. Finally we discuss how, by using a passive scalar suspension model with aggregation on top of the LB dynamics, we can describe the clotting processes in the aneurysm     

复杂流动的LBGK模拟

ＬＢＧＫ方法在保留了格子气模型的一些优点的同时，克服了格子气模型的不足之处．应用ＬＢＧＫ方法，可以对更复杂的流体进行模拟．本文讨论了一种九点ＬＢＧＫ模型，并用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开方法和多尺度技术证明其在二阶精度上表现为标准的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程．用该模型对圆柱绕流和管排流动进行模拟的结果显示，该模型能较好地模拟复杂流动现象，并具有一定工程应用背景．     

Two-dimensional MRT LB model for compressible and incompressible flows

In the paper we extend the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model pro- posed in [Europhys. Lctt., 2010, 90： 54003] so that it is suitable also for incompressible flows. To decrease tile artificial oscillations, the convection term is discretized by the flux linfiter scheme with splitting technique. A new model is validated by some well-known benchmark tests, including Rie- mann problem and Couette flow, and satisfying agreements are obtained between the sinmlation results and ana.lytical ones. In order to show the merit of LB model over traditional methods, the non-equilibrium characteristics of system are solved. The simulation results are consistent with the physical analysis.     

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional (3D) Lattice Boltzmann (LB) model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 (2004) 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative (NND) scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model isvalidated by well-known benchmarks, (i) Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; (ii) reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.     

Lattice Boltzmann modeling with discontinuous collision components: Hydrodynamic and Advection-Diffusion Equations

Irrespective of the nature of the modeled conservation laws, we establish first the microscopic interface continuity conditions for Lattice Boltzmann (LB) multiple-relaxation time, link-wise collision operators with discontinuous components (equilibrium functions and/or relaxation parameters). Effective macroscopic continuity conditions are derived for a planar implicit interface between two immiscible fluids, described by the simple two phase hydrodynamic model, and for an implicit interface boundary between two heterogeneous and anisotropic, variably saturated soils, described by Richard’s equation. Comparing the effective macroscopic conditions to the physical ones, we show that the range of the accessible parameters is restricted, e.g. a variation of fluid densities or a heterogeneity of the anisotropic soil properties. When the interface is explicitly tracked, the interface collision components are derived from the leading order continuity conditions. Among particular interface solutions, a harmonic mean value is found to be an exact LB solution, both for the interface kinematic viscosity and for the interface vertical hydraulic conductivity function. We construct simple problems with the explicit and implicit interfaces, matched exactly by the LB hydrodynamic and/or advection-diffusion schemes with the aid of special solutions for free collision parameters.     

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer-Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. 90 (2010) 54003]. We discuss the reason why the Gram-Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 (2004) 035701(R)] is only valid in subsonic flows.The von Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer-Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.     

Simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

Fresh cement mortar is a type of workable paste,which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering.In this paper,Papanastasiou’s model for Bingham fluids is solved by using the multiplerelaxation-time lattice Boltzmann model(MRT-LB).Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou’s model provides a good approximation of realistic Bingham plastics for values of m108.For lower values of m,Papanastasiou’s model is valid for fluids between Bingham and Newtonian fluids.The MRT-LB model is validated by two benchmark problems:2D steady Poiseuille flows and lid-driven cavity flows.Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability.We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle.Besides the rich flow structures obtained in this work,the dynamics fluid force on the round particle is calculated.Results show that both the Reynolds number Re and the Bingham number Bn afect the drag coefcients CD,and a drag coefcient with Re and Bn being taken into account is proposed.The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed.Finally,the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields.These results help simulation of fresh concrete flowing in porous media.     

Investigation of cavitation bubble collapse in hydrophobic concave using the pseudopotential multi-relaxation-time lattice Boltzmann method

The interaction between cavitation bubble and solid surface is a fundamental topic which is deeply concerned for the utilization or avoidance of cavitation effect.The complexity of this topic is that the cavitation bubble collapse includes many extreme physical phenomena and variability of different solid surface properties.In the present work,the cavitation bubble collapse in hydrophobic concave is studied using the pseudopotential multi-relaxation-time lattice Boltzmann model(MRT-LB).The model is modified by involving the piecewise linear equation of state and improved forcing scheme.The fluid-solid interaction in the model is employed to adjust the wettability of solid surface.Moreover,the validity of the model is verified by comparison with experimental results and grid-independence verification.Finally,the cavitation bubble collapse in a hydrophobic concave is studied by investigating density field,pressure field,collapse time,and jet velocity.The superimposed effect of the surface hydrophobicity and concave geometry is analyzed and explained in the framework of the pseudopotential LBM.The study shows that the hydrophobic concave can enhance cavitation effect by decreasing cavitation threshold,accelerating collapse and increasing jet velocity.     

Multi-relaxation-time lattice Boltzmann model for incompressible flow

In this Letter an incompressible MRT-LB model has been proposed. The equilibria in momentum space are derived from an earlier incompressible LBGK model by Guo et al. Through the Chapman–Enskog expansion the incompressible Navier–Stokes equations can be recovered without artificial compressible effects. The steady Poiseuille flow, the driven cavity flow and the double shear flow have been carried on by the incompressible MRT-LB model. The numerical simulation results agree well with the analytical solutions or the existing results. It is found that the incompressible MRT-LB model shows better numerical stability.     

A note on equilibrium boundary conditions in lattice Boltzmann fluid dynamic simulations

It is shown that Lattice Boltzmann (LB) simulations using simple equilibrium boundary conditions at solid walls, provide quantitatively accurate results for backward-facing step flows at moderate Reynolds numbers. The basic reason for such favorable behavior is that well-resolved LB simulations operate in the so-called weak non-equilibrium regime, in which shear effects at the scale of a single lattice spacing are weak, meaning by this that the cell-shear time scale is much longer than the molecular time scales, so that the LB collisional relaxation takes place in a quasi-homogeneous velocity field. Due to their simplicity, it is suggested that equilibrium boundary conditions may represent a viable option for the LB simulation of complex flows with solid boundaries at moderate Reynolds numbers.     

Simulation of blood flow using extended Boltzmann kinetic approach

Lattice Boltzmann (LB) simulations are conducted to obtain the detailed hydrodynamics in a variety of blood vessel setups, including a prototype stented channel and four human coronary artery geometries based on the images obtained from real patients. For a model of stented flow involving an S-shape stent, a pulsatile flow rate is applied as the inlet boundary condition, and the time- and space-dependent flow field is computed. The LB simulation is found to reproduce the analytical solutions for the velocity profiles and wall shear stress distributions for the pulsatile channel flow. For the coronary arteries, the distributions of wall shear stress, which is important for clinical diagnostic purposes, are in good agreement with the conventional CFD predictions.     

Numerical Stability Analysis of FDLBM

We analyze the numerical stability of Finite Difference Lattice Boltzmann Method (FDLBM) by means of von Neumann stability analysis. The stability boundary of the FDLBM depends on the BGK relaxation time, the CFL number, the mean flow velocity, and the wavenumber. As the BGK relaxation time is increased at constant CFL number, the stability of the central difference LB scheme may not be ensured. The limits of maximum stable velocity are obtained around 0.39, 0.43, and 0.43 for the central difference, for the explicit upwind difference, and for the semi-implicit upwind difference schemes, respectively. We derive artificial viscosities for every difference scheme and investigate their influence on numerical stability. The requirements for artificial viscosity is consistent with the conditions derived from von Neumann stability analysis. This analysis elucidates that the upwind difference schemes are suitable for simulation of high Reynolds number flows.     

Modeling and simulation of thermocapillary flows using lattice Boltzmann method

To understand how thermocapillary forces manipulate droplet motion in microfluidic channels, we develop a lattice Boltzmann (LB) multiphase model to simulate thermocapillary flows. The complex hydrodynamic interactions are described by an improved color-fluid LB model, in which the interfacial tension forces and the Marangoni stresses are modeled in a consistent manner using the concept of a continuum surface force. An additional convection–diffusion equation is solved in the LB framework to obtain the temperature field, which is coupled to the interfacial tension through an equation of state. A stress-free boundary condition is also introduced to treat outflow boundary, which can conserve the total mass of an incompressible system, thus improving the numerical stability for creeping flows.The model is firstly validated against the analytical solutions for the thermocapillary driven convection in two superimposed fluids at negligibly small Reynolds and Marangoni numbers. It is then used to simulate thermocapillary migration of three-dimensional deformable droplet at various Marangoni numbers, and its accuracy is once again verified against the theoretical prediction in the limit of zero Marangoni number. Finally, we numerically investigate how the localized heating from a laser can block the microfluidic droplet motion through the induced thermocapillary forces. The droplet motion can be completely blocked provided that the intensity of laser exceeds a threshold value, below which the droplet motion successively undergoes four stages: constant velocity, deceleration, acceleration, and constant velocity. When the droplet motion is completely blocked, four steady vortices are clearly visible, and the droplet is fully filled by two internal vortices. The external vortices diminish when the intensity of laser increases.     

Superdiffusivity of Two Dimensional Lattice Gas Models

It was proved [Navier–Stokes Equations for Stochastic Particle System on the Lattice. Comm. Math. Phys. (1996) 182, 395; Lattice gases, large deviations, and the incompressible Navier–Stokes equations. Ann. Math. (1998) 148, 51] that stochastic lattice gas dynamics converge to the Navier–Stokes equations in dimension d=3 in the incompressible limits. In particular, the viscosity is finite. We proved that, on the other hand, the viscosity for a two dimensional lattice gas model diverges faster than (log t)^1/2. Our argument indicates that the correct divergence rate is (log t)^1/2. This problem is closely related to the logarithmic correction of the time decay rate for the velocity auto-correlation function of a tagged particle.     

Accuracy analysis of high-order lattice Boltzmann models for rarefied gas flows

In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is shown to be able to reduce to the linearized Bhatnagar–Gross–Krook (BGK) equation. Therefore, when the same Gauss–Hermite quadrature is used, LB method closely resembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be directly correlated with the approximation order in terms of the Knudsen number to the BGK equation for incompressible flows. Meanwhile, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role in capturing non-equilibrium effect for low-speed flows. By contrast, appropriate Gauss–Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. Our simulation results, where the effect of wall/gas interactions is excluded, can lead to conclusion on the LB modeling capability that the models with higher-order quadratures provide more accurate results. For the same order Gauss–Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss–Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at the Navier–Stokes level and rarefied gas flows at the linearized Boltzmann model equation level.     

A simple enthalpy-based lattice Boltzmann scheme for complicated thermal systems

To extend the lattice Boltzmann (LB) method to describe the applicable energy systems, the first key step is to build a suitable thermal LB model and corresponding boundary treatments. There are two main shortcomings in the existing related works: either some additional energy source terms are inconvenient to be naturally incorporated or the implementation of non-Dirichlet-type thermal boundary conditions is extremely difficult and sometimes impossible in them for complicated thermal systems, which restrict their applicability to only a few special classes of problems. In order to overcome these drawbacks by a simple way, in this paper a thermal LB model and corresponding boundary treatments are constructed based on the total enthalpy. The specific benefits due to the introduction of the total enthalpy are analyzed and it is found that the numerical results obtained by the present scheme agree well with the analytical solutions and/or the data reported in previous studies.     

Error analysis and correction for Lattice Boltzmann simulated flow conductance in capillaries of different shapes and alignments

We study the relative error in conductance calculations, for simulated flow of a single component single phase fluid through a capillary in three dimensions, by the Lattice Boltzmann (LB) method with bounce-back boundary conditions. The relative error with respect to analytical results for capillary cross-sections of circular, triangular and square shapes are calculated as a function of the cross-section diameter, a, and for different alignment of the cross-section relative to the underlying lattice grid. It is shown, when the shapes are not aligned perfectly to the lattice, that the relative error decreases systematically with the size, a, as ～1/a when a is evaluated by mapping the computed cross-sectional area, in terms of the enclosed number of grid points, to the respective geometrical shapes concerned. For perfectly aligned geometries, viz. the square capillary aligned to the LB lattice grid or rotated with its side along the diagonal of the LB grid, the relative error decreases as ～1/a². A simple method is suggested to locate the boundary wall depending on its orientation relative to the grid, such that the exact conductance of the new shape matches the LB computed conductance.     

A velocity-estimation subgrid model constrained by subgrid scale dissipation

Purely dissipative eddy-viscosity subgrid models have proven very successful in large-eddy simulations (LES) at moderate resolution. Simulations at coarse resolutions where the underlying assumption of small-scale universality is not valid, warrant more advanced models. However, non-eddy viscosity models are often unstable due to the lack of sufficient dissipation. This paper proposes a simple modeling approach which incorporates the dissipative nature of existing eddy viscosity models into more physically appealing non-eddy viscosity SGS models. The key idea is to impose the SGS dissipation of the eddy viscosity model as a constraint on the non-eddy viscosity model when determining the coefficients in the non-eddy viscosity model. We propose a new subgrid scale model (RSEM), which is based on estimation of the unresolved velocity field. RSEM is developed in physical space and does not require the use of finer grids to estimate the subgrid velocity field. The model coefficient is determined such that total SGS dissipation matches that from a target SGS model in the mean or least-squares sense. The dynamic Smagorinsky model is used to provide the target dissipation. Results are shown for LES of decaying isotropic turbulence and turbulent channel flow. For isotropic turbulence, RSEM displays some level of backward dissipation, while yielding as good results as the dynamic Smagorinsky model. For channel flow, the results from RSEM are better than those from the dynamic Smagorinsky model for both statistics and instantaneous flow structures.