We show that an enslaved phase-separation front moving with diffusive speeds can leave alternating domains of increasing size in their wake. We find the size and spacing of these domains is identical to Liesegang patterns. For equal composition of the components we are able to predict the exact form of the pattern analytically. To our knowledge this is the first fully analytical derivation of the Liesegang laws. We also show that there is a critical value for C below which only two domains are formed. Our analytical predictions are verified by numerical simulations using a lattice Boltzmann method.

The dynamic behaviors of continuous and discrete flows in superhydrophobic microchannels are investigated with a lattice Boltzmann model. Typical characters of the superhydrophobic phenomenon are well observed from our simulations, including air trapped in the surface microstructures, high contact angles, low contact angle hysteresis, and reduced friction to fluid motions. Increasing the roughness of a hydrophobic surface can produce a large flow rate through the channel due to the trapped air, implying less friction or large apparent slip. The apparent slip length appears to be independent to the channel width and could be considered as a surface property. For a moving droplet, its behavior is affected by the surface roughness from two aspects: the contact angle difference between its two ends and the surface-liquid interfacial friction. As a consequence, the resulting droplet velocity changes with the surface roughness as firstly decreasing and then increasing. Simulation results are also compared with experimental observations and better agreement has been obtained than that from other numerical method. The information from this study could be valuable for microfluidic systems.

This work proposes an extension to Boltzmann BGK equation for dense gases. The present model has an H-theorem and it allows choice of the Prandtl number as an independent parameter. I show that similar to Enskog equation this equation can reproduce dynamics of dense gases.

We show the methodology and advantages of asymptotic analysis when applied to lattice Boltzmann outflow treatments. On the one hand, one can analyze outflow algorithms formulated directly in terms of the lattice Boltzmann variables, like the extrapolation method, to find the induced outflow conditions in terms of the Navier-Stokes variables. On the other hand, one can check the consistency and accuracy of lattice Boltzmann outflow treatments to given hydrodynamic outflow conditions like the Neumann or average pressure condition. As example how the gained insight can be used, we propose an improvement of the well known extrapolation method.

The exact solution to the hierarchy of nonlinear lattice Boltzmann kinetic equations, for the stationary planar Couette flow for any Knudsen number was presented by S. Ansumali et al. [Phys. Rev. Lett., 98 (2007), 124502]. In this paper, simulation results at a non-vanishing value of the Knudsen number are compared with the closed-form solutions for the higher-order moments. The order of convergence to the exact solution is also studied. The lattice Boltzmann simulations are in excellent agreement with the exact solution.

Formal analogies between the Car-Parrinello (CP) ab-initio molecular dynamics for quantum many-body systems, and the Lattice Boltzmann (LB) method for classical and quantum fluids, are pointed out. A theoretical scenario, whereby the quantum LB would be coupled to the CP framework to speed-up many-body quantum simulations, is also discussed, together with accompanying considerations on the computational efficiency of the prospective CP-LB scheme.

Homogeneous isotropic turbulence has been playing a key role in the research of turbulence theory. And the pseudo-spectral method is the most popular numerical method to simulate this type of flow fields in a periodic box, where fast Fourier transform (FFT) is mostly effective. However, the bottle-neck in this method is the memory of computer, which motivates us to construct a memory-saving algorithm for spectral method in present paper. Inevitably, more times of FFT are needed as compensation. In the most memory-saving situation, only 6 three-dimension arrays are employed in the code. The cost of computation is increased by a factor of 4, and that 38 FFTs are needed per time step instead of the previous 9 FFTs. A simulation of isotropic turbulence on 20483 grid can be implemented on a 256G distributed memory clusters through this method.

The lattice Boltzmann method (LBM) for multicomponent immiscible fluids is applied to the simulations of solid-fluid mixture flows including spherical or non-spherical particles in a square pipe at Reynolds numbers of about 100. A spherical solid particle is modeled by a droplet with strong interfacial tension and large viscosity, and consequently there is no need to track the moving solid-liquid boundary explicitly. Nonspherical (discoid, flat discoid, and biconcave discoid) solid particles are made by applying artificial forces to the spherical droplet. It is found that the spherical particle moves straightly along a stable position between the wall and the center of the pipe (the Segré-Silberberg effect). On the other hand, the biconcave discoid particle moves along a periodic helical path around the center of the pipe with changing its orientation, and the radius of the helical path and the polar angle of the orientation increase as the hollow of the concave becomes large.

In this paper we present a generalized lattice Boltzmann based approach for sliding-mesh local reference frame. This scheme exactly conserves hydrodynamic fluxes across local reference frame interface. The accuracy and robustness of our scheme are demonstrated by benchmark validations.

Water management is a key to ensuring high performance and durability of polymer electrolyte fuel cell (PEFC), and it is important to understand the behavior of liquid water in PEFC. In this study, the two-phase lattice Boltzmann method is applied to the simulations of water discharge from gas diffusion layers (GDL) to gas channels. The GDL is porous media composed of carbon fibers with hydrophobic treatment, and the gas channels are hydrophilic micro-scale ducts. In the simulations, arbitrarily generated porous materials are used as the structures of the GDL. We investigate the effects of solid surface wettabilities on water distribution in the gas channels and the GDL. Moreover, the results of X-ray computed tomography images in the operating PEFC are compared with the numerical simulations, and the mechanism of the water transport in PEFC is considered.

Realization of the Shan-Chen multiphase flow lattice Boltzmann model is considered in the framework of the higher-order Galilean invariant lattices. The present multiphase lattice Boltzmann model is used in two-dimensional simulation of droplet collisions at high Weber numbers. Results are found to be in a good agreement with experimental findings.

Using the lattice Boltzmann multiphase model, numerical simulations have been performed to understand the dynamics of droplet formation in a microfluidic cross-junction. The influence of capillary number, flow rate ratio, viscosity ratio, and viscosity of the continuous phase on droplet formation has been systematically studied over a wide range of capillary numbers. Two different regimes, namely the squeezinglike regime and the dripping regime, are clearly identified with the transition occurring at a critical capillary number Cacr. Generally, large flow rate ratio is expected to produce big droplets, while increasing capillary number will reduce droplet size. In the squeezing-like regime (Ca ≤ Cacr), droplet breakup process is dominated by the squeezing pressure and the viscous force; while in the dripping regime (Ca ≤ Cacr), the viscous force is dominant and the droplet size becomes independent of the flow rate ratio as the capillary number increases. In addition, the droplet size weakly depends on the viscosity ratio in both regimes and decreases when the viscosity of the continuous phase increases. Finally, a scaling law is established to predict the droplet size.

The objective of this paper is to seek an alternative to the numerical simulation of the Navier-Stokes equations by a method similar to solving the BGK-type modeled lattice Boltzmann equation. The proposed method is valid for both gas and liquid flows. A discrete flux scheme (DFS) is used to derive the governing equations for two distribution functions; one for mass and another for thermal energy. These equations are derived by considering an infinitesimally small control volume with a velocity lattice representation for the distribution functions. The zero-order moment equation of the mass distribution function is used to recover the continuity equation, while the first-order moment equation recovers the linear momentum equation. The recovered equations are correct to the first order of the Knudsen number (Kn); thus, satisfying the continuum assumption. Similarly, the zero-order moment equation of the thermal energy distribution function is used to recover the thermal energy equation. For aerodynamic flows, it is shown that the finite difference solution of the DFS is equivalent to solving the lattice Boltzmann equation (LBE) with a BGK-type model and a specified equation of state. Thus formulated, the DFS can be used to simulate a variety of aerodynamic and hydrodynamic flows. Examples of classical aeroacoustics, compressible flow with shocks, incompressible isothermal and non-isothermal Couette flows, stratified flow in a cavity, and double diffusive flow inside a rectangle are used to demonstrate the validity and extent of the DFS. Very good to excellent agreement with known analytical and/or numerical solutions is obtained; thus lending evidence to the DFS approach as an alternative to solving the Navier-Stokes equations for fluid flow simulations.

A new boundary condition, aimed at inhibiting near-wall condensation effects in lattice Boltzmann simulations of capillary flows in micro-corrugated channels, is introduced. The new boundary condition is validated against analytical solutions for smooth channels and demonstrated for the case of three-dimensional microflows over randomly corrugated walls.

Recent studies on applications of the lattice Boltzmann method (LBM) and the finite-difference lattice Boltzmann method (FDLBM) to velocity slip simulations are mostly on one-dimensional (1D) problems such as a shear flow between parallel plates. Applications to a 2D problem may raise new issues. The author performed numerical simulations of rotational slip flow in coaxial cylinders as an example of 2D problem. Two types of 2D models were used. The first were multi-speed FDLBM models proposed by the author. The second was a standard LBM, the D2Q9 model. The simulations were performed applying a finite difference scheme to both the models. The study had two objectives. The first was to investigate the accuracies of LBM and FDLBM on applications to rotational slip flow. The second was to obtain an experience on application of the cylindrical coordinate system. The FDLBM model with 8 directions and the D2Q9 model showed an anisotropic flow pattern when the relaxation time constant or the Knudsen number was large. The FDLBM model with 24 directions showed accurate results even at large Knudsen numbers.

We analyze the Lattice Boltzmann method for the simulation of fluctuating hydrodynamics by Adhikari et al. [Europhys. Lett., 71 (2005), 473-479] and find that it shows excellent agreement with theory even for small wavelengths as long as a stationary system is considered. This is in contrast to other finite difference and older lattice Boltzmann implementations that show convergence only in the limit of large wavelengths. In particular cross correlators vanish to less than 0.5%. For larger mean velocities, however, Galilean invariance violations manifest themselves.

In order to model the dispersal of volcanic particles in the atmosphere and their deposition on the ground, one has to simulate an advection-diffusion-sedimentation process on a large spatial area. Here we compare a Lattice Boltzmann and a Cellular Automata approach. Our results show that for high Peclet regimes, the cellular automata model produce results that are as accurate as the lattice Boltzmann model and is computationally more effective.

The lattice Boltzmann method is applied to the investigations of the diffusivity and the permeability in the gas diffusion layer (GDL) of the polymer electrolyte fuel cell (PEFC). The effects of the configuration of water droplets, the porosity of the GDL, the viscosity ratio of water to air, and the surface wettability of the GDL are investigated. From the simulations under the PEFC operating conditions, it is found that the heterogeneous water network and the high porosity improve the diffusivity and the permeability, and the hydrophobic surface decreases the permeability.