Gary D. Doolen

Multicomponent lattice-Boltzmann model with interparticle interaction

A lattice Boltzmann model for simulating fluids with multiple components and interparticle forces proposed by Shan and Chen is described in detail. Macroscopic equations governing the motion of each component are derived by using the Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirment by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity, so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confiremoed numerically.

Lattice-Boltzmann Simulations of Fluid Flows in MEMS

The lattice Boltzmann model is a simplified kinetic method based on the particle distribution function. We use this method to simulate problems in MEMS, in which the velocity slip near the wall plays an important role. It is demonstrated that the lattice Boltzmann method can capture the fundamental behaviors in micro-channel flow, including velocity slip, nonlinear pressure drop along the channel and mass flow rate variation with Knudsen number. The Knudsen number dependence of the position of the vortex center and the pressure contour in micro-cavity flows is also demonstrated.

Lattice Boltzmann Computational Fluid Dynamics in Three Dimensions

The recent development of the lattice gas method and its extension to the lattice Boltzmann method have provided new computational schemes for fluid dynamics. Both methods are fully paralleled and can easily model many different physical problems, including flows with complicated boundary conditions. In this paper, basic principles of a lattice Boltzmann computational method are described and applied to several three-dimensional benchmark problems. In most previous lattice gas and lattice Boltzmann methods, a face-centered-hyper-cubic lattice in four-dimensional space was used to obtain an isotropic stress tensor. To conserve computer memory, we develop a model which requires 14 moving directions instead of the usual 24 directions. Lattice Boltzmann models, describing two-phase fluid flows and magnetohydrodynamics, can be developed based on this simpler 14-directional lattice. Comparisons between three-dimensional spectral code results and results using our method are given for simple periodic geometries. An important property of the lattice Boltzmann method is that simulations for flow in simple and complex geometries have the same speed and efficiency, while all other methods, including the spectral method, are unable to model complicated geometries efficiently.

Thermodynamic Foundations of Kinetic Theory and Lattice Boltzmann Models for Multiphase Flows

This paper demonstrates that thermodynamically consistent lattice Boltzmann models for single-component multiphase flows can be derived from a kinetic equation using both Enskogs theory for dense fluids and mean-field theory for long-range molecular interaction. The lattice Boltzmann models derived this way satisfy the correct mass, momentum, and energy conservation equations. All the thermodynamic variables in these LBM models are consistent. The strengths and weaknesses of previous lattice Boltzmann multiphase models are analyzed.

On statistical correlations between velocity increments and locally averaged dissipation in homogeneous turbulence

Kolmogorov postulated in 1962 [J. Fluid Mech. 13, 82 (1962)] that the magnitude of velocity increments δur across an inertial range distance r in high Reynolds number flows is typically (rer)1/3, where er is the locally averaged dissipation rate. This refined similarity hypothesis has been widely used in discussions of anomalous exponents of velocity structure functions in connection with the scaling exponents of er. Recently Hosokawa and Yamamoto [Phys. Fluids A 4, 457 (1992)] have presented numerical evidence from turbulence simulations that δur is uncorrelated with er in moderate Reynolds number flows. In the present paper, results of similar measurements are offered for flow fields with a wide range of Reynolds numbers obtained from high‐resolution numerical simulations of both forced and decaying isotropic turbulence. The present results show clear evidence of correlations between δur and er, irrespective of the Reynolds number. Kolmogorov’s hypothesis is verified for r somewhat larger than the visco...

A improved incompressible lattice Boltzmann model for time-independent flows

It is well known that the lattice Boltzmann equation method (LBE) can model the incompressible Navier-Stokes (NS) equations in the limit where density goes to a constant. In a LBE simulation, however, the density cannot be constant because pressure is equal to density times the square of sound speed, hence a compressibility error seems inevitable for the LBE to model incompressible flows. This work uses a modified equilibrium distribution and a modified velocity to construct an LBE which models time-independent (steady) incompressible flows with significantly reduced compressibility error. Computational results in 2D cavity flow and in a 2D flow with an exact solution are reported.

Lattice methods and their applications to reacting systems

Abstract The recent development of the lattice gas automata method and its extension to the lattice Boltzmann method have provided new computational schemes for solving a variety of partial differential equations and modeling chemically reacting systems. The lattice gas method, regarded as the simplest microscopic and kinetic approach which generates meaningful macroscopic dynamics, is fully parallel and can, as a result, be easily programmed on parallel machines. In this paper, we introduce the basic principles of the lattice gas method and the lattice Boltzmann method, their numerical implementations and applications to chemically reacting systems. Comparisons of the lattice Boltzmann method with the lattice gas technique and other traditional numerical schemes, including the finite difference scheme and the pseudo-spectral method, for solving the Navier-Stokes hydrodynamic fluid flows will be discussed. Recent developments of the lattice gas and the lattice Boltzmann method and their applications to pattern formation in chemical reaction-diffusion systems, multiphase fluid flows and polymeric dynamics will be presented.

On the three-dimensional Rayleigh–Taylor instability

The three-dimensional Rayleigh–Taylor instability is studied using a lattice Boltzmann scheme for multiphase flow in the nearly incompressible limit. This study focuses on the evolution of the three-dimensional structure of the interface. In addition to the bubble and spike fronts, a saddle point is found to be another important landmark on the interface. Two layers of heavy-fluid roll-ups, one at the spike tip and the other at the saddle point, were observed. The secondary instability in the horizontal planes entangles the already complicated structure of the interface. Parallel computations are utilized to accommodate the massive computational requirements of the simulations.

Lattice gas automata for flow through porous media

Abstract Lattice gas hydrodynamic models for flows through porous media in two and three dimensions are described. The computational method easily handles arbitrary boundaries and a large range of Reynolds numbers. Darcys law is confirmed for Poiseuille flow and for complicated boundary flows. Multiply connected pore structures similar to actual sandstone with fixed fractal dimension and porosity are generated. Permeability as a function of fractal dimension and porosity is calculated and compared with results of other methods and experiments.

Simulation of Combustion Field with Lattice Boltzmann Method

Turbulent combustion is ubiquitously used in practical combustion devices. However, even chemically non-reacting turbulent flows are complex phenomena, and chemical reactions make the problem even more complicated. Due to the limitation of the computational costs, conventional numerical methods are impractical in carrying out direct 3D numerical simulations at high Reynolds numbers with detailed chemistry. Recently, the lattice Boltzmann method has emerged as an efficient alternative for numerical simulation of complex flows. Compared with conventional methods, the lattice Boltzmann scheme is simple and easy for parallel computing. In this study, we present a lattice Boltzmann model for simulation of combustion, which includes reaction, diffusion, and convection. We assume the chemical reaction does not affect the flow field. Flow, temperature, and concentration fields are decoupled and solved separately. As a preliminary simulation, we study the so-called “counter-flow” laminar flame. The particular flow geometry has two opposed uniform combustible jets which form a stagnation flow. The results are compared with those obtained from solving Navier–Stokes equations.