H-theorem and generalized semi-detailed balance condition for lattice gas systems

For lattice gas systems obeying Fermi-Dirac statistics, an H-theorem can be proved with a more general condition that the semi-detailed balance condition. This new condition allows more flexible transition rates among states, so that it has broader applicability for various lattice gas models, including those which have multiple phase properties.     

Negative-viscosity lattice gases

A new irreversible collision rule is introduced for lattice-gas automata. The rule maximizes the flux of momentum in the direction of the local momentum gradient, yielding a negative shear viscosity. Numerical results in 2D show that the negative viscosity leads to the spontaneous ordering of the velocity field, with vorticity resolvable down to one lattice-link length. The new rule may be used in conjunction with previously proposed collision rules to yield a positive shear viscosity lower than the previous rules provide. In particular, Poiseuille flow tests demonstrate a decrease in viscosity by more than a factor of 2.     

Lattice gas with “interaction potential”

We present an extension of a simple automaton model to incorporate nonlocal interactions extending over a spatial range in lattice gases. From the viewpoint of statistical mechanics, the lattice gas with interaction range may serve as a prototype for non-ideal-gas behavior. From the density fluctuations correlation function, we obtain a quantity which is identified as a potential of mean force. Equilibrium and transport properties are computed theoretically and by numerical simulations to establish the validity of the model at macroscropic scale.     

Simulating three-dimensional hydrodynamics on a cellular automata machine

We demonstrate how three-dimensional fluid flow simulations can be carried out on the Cellular Automata Machine 8 (CAM-8), a special-purpose computer for cellular automata computations. The principal algorithmic innovation is the use of a lattice gas model with a 16-bit collision operator that is specially adapted to the machine architecture. It is shown how the collision rules can be optimized to obtain a low viscosity of the fluid. Predictions of the viscosity based on a Boltzmann approximation agree well with measurements of the viscosity made on CAM-8. Several test simulations of flows in simple geometries—channels, pipes, and a cubic array of spheres-are carried out. Measurements of average flux in these geometries compare well with theoretical predictions.     

Global invariants and equilibrium states in lattice gases

It is now well known that, in addition to the physical conserved quantities, lattice gases also have other unphysical ones related to the discretization of their phase space. From an abstract point of view a lattice gas can be considered like a full discrete Markov processL and these spurious conserved quantities yield the existence of a nonspatially homogeneous equilibrium state forL ^k. We show that a particular set of these conserved quantities is of special interest: Its elements will be called regular. These regular invariants are simply built from the local ones and their projection on each node is always a locally conserved quantity. Moreover, for most models they are one-to-one related to the Gibbs states ofL ^k which remain factorized. It turns out that all the classical known spurious invariants are regular and one can exhibit simple conditions to build models with only regular invariants. For the latter it is then justified to determine the transport coefficients of the locally conserved densities with the Green-Kubo procedure.     

Diffusion simulation with a deterministic one-dimensional lattice-gas model

A one-dimensional lattice-gas model is proposed and used to simulate diffusion processes in one dimension. Explicit forms of transport coefficients are given as a function of density and kinetic energy within the Boltzmann approximation. Without definitions of temperature and pressure, a steady nontrivial solution is given analytically in the nonconvective case when the kinetic energy is kept constant.     

A stochastic lattice gas for Burgers' equation: A practical study

We continue our investigation of stochastic lattice gases as a (highly parallel) means of simulating given PDEs, in this case Burgers' equation in one dimension. The lattice dynamics consists of stochastic unidirectional particle displacement, and our attention is turned toward the reliability of the model, i.e., its ability to reproduce the unique physical solution of Burgers' equation. Lattice gas results are discussed and compared against finite-difference calculations and exact solutions in examples which include shocks and rarefaction waves.     

Gauge-invariant lattice gas for the microcanonical Ising model

We introduce a lattice gas for particles with discrete momenta (1, 0, –1) and local deterministic microdynamics, which exactly reproduces Creutz's microcanonical algorithm for the ferromagnetic Ising model. However, because of the manifest gauge invariance of our variables, both the Ising ferromagnetic and spin-glass systems share precisely the same dynamics with different initial conditions. Additional conservation laws in the 1D Ising case result in a completely integrable system in the limit of zero or unbounded demon energy cutoff. Numerical investigations of ergodicity are presented for the pure Ising lattice gas in one and two dimensions.     

Simulation for pedestrian dynamics by real-coded cellular automata (RCA)

In this paper, we propose a new approach for pedestrian dynamics. We call it a real-coded cellular automata (RCA). The scheme is based on the real-coded lattice gas (RLG), which has been developed for fluid simulation. Similar to RLG, the position and velocity can be freely given, independent of grid points. Our strategy including the procedure for updating the position of each pedestrian is explained. It is shown that the movement of pedestrians in an oblique direction to the grid is successfully simulated by RCA, which was not taken into account in the previous CA models. Moreover, from simulations of evacuation from a room with an exit of various widths, we obtain the critical number of people beyond which the clogging appears at the exit.     

Calculation of the permeability of porous media using hydrodynamic cellular automata

The permeability of two-dimensional porous media is calculated numerically as a function of porosity using the hydrodynamic cellular automata (lattice gas) approach. Results are presented for systems with up to 22 million sites (8192×2688). For randomly distributed solid obstacles whose macroscopic dimensions are much longer than the mean free path of particles in the fluid, the permeability varies with porosity as (–0.6)/(1–) for>0.7. When the solid obstacles are much smaller than the mean free path of particles in the fluid, i.e., when they form a dust of point objects, then such a relationship no longer holds and the permeability is more than an order of magnitude smaller than for the former case. The program used for the simulations is discussed and a listing is presented in the Appendix which achieved a sustained speed of 185 million sites updated per second on a single processor of the Cray-YMP. (On a Sun Sparc Workstation, the same program ran about 100 times slower.)     

Lattice Boltzmann simulations of micron-scale drop impact on dry surfaces

A lattice Boltzmann equation (LBE) method for incompressible binary fluids is proposed to model the contact line dynamics on partially wetting surfaces. Intermolecular interactions between a wall and fluids are represented by the inclusion of the cubic wall energy in the expression of the total free energy. The proposed boundary conditions eliminate the parasitic currents in the vicinity of the contact line. The LBE method is applied to micron-scale drop impact on dry surfaces, which is commonly encountered in drop-on-demand inkjet applications. For comparison with the existing experimental results [H. Dong, W.W. Carr, D.G. Bucknall, J.F. Morris, Temporally-resolved inkjet drop impaction on surfaces, AIChE J. 53 (2007) 2606–2617], computations are performed in the range of equilibrium contact angles from 31° to 107° for a fixed density ratio of 842, viscosity ratio of 51, Ohnesorge number (Oh) of 0.015, and two Weber numbers (We) of 13 and 103.     

Implementation of the FCHC lattice gas model on the Connection Machine

The 4-dimensional FCHC lattice gas model has been implemented on a Connection Machine CM-2 with 16K processors. Symmetries are used to reduce the collision table to a size that fits into local memory. This method avoids the degradation of the Reynolds coefficientR _*, but at the price of increased computing time. Bit shuffling between parallel lattices is introduced to reduce the discrepancy between measured viscosities and those predicted from the Boltzmann approximation. Thereby a model with a negative shear viscosity is obtained: a fluid having a uniform initial velocity is unstable and organized nonuniform motions develop. Because of the buildup of very strong correlations between the parallel lattices, the discrepancy with the Boltzmann values decreases only very slowly with the number of parallel lattices.     

A dilute bootstrap percolation: Lattice model for unresponsiveness in T-cell immunology

The biased majority rule of cellular automata takes a spin up if and only if at least two of its four nearest neighbors on the square lattice are up. We generalize this type of bootstrap percolation by introducing quenched site dilution as well as a random birth and decay process. Our Monte Carlo simulations then give first-order transitions qualitatively similar to our results from meanfield reaction equations describing the induction of T-cell unresponsiveness in the immune system.     

Lattice simulations for few- and many-body systems

We review the recent literature on lattice simulations for few- and many-body systems. We focus on methods that combine the framework of effective field theory with computational lattice methods. Lattice effective field theory is discussed for cold atoms as well as low-energy nucleons with and without pions. A number of different lattice formulations and computational algorithms are considered, and an effort is made to show common themes in studies of cold atoms and low-energy nuclear physics as well as common themes in work by different collaborations.     

Lattice gas simulation and experiment study of evacuation dynamics

In this paper, evacuation dynamics in an office building is studied by experiment and simulation. A lattice gas (LG) model is developed. A parameter called ‘exit bias’ is introduced into the model to describe the occupants’ familiarity with different exits in a building. The evacuation experiment, which consists of seven scenarios under various conditions, is conducted to verify the model and calibrate the model’s input parameters such as pedestrian speed and exit bias. The effect of exit width on flow rate, and the effect of occupants’ familiarity with the building on their route selections, are studied. It is found that the accuracy of simulation depends a lot on the model’s pedestrian speed. The optimal pedestrian speed is decided by not only occupant characteristics, but also flow features determined by people distribution, building structure, environment pressure, etc. LG models with proper pedestrian speed are capable of simulating the dynamic process of orderly emergency evacuations.     

Lattice gas automata for simple and complex fluids

We review some recent applications of lattice gas automata, including flow through porous media, phase transitions, thermodydrodynamics, and magnetohydrodynamics.     

Lattice gases and exactly solvable models

We detail the construction of a family of lattice gas automata based on a model of 't Hooft, proceeding by use of symmetry principles to define first the kinematics of the model and then the dynamics. A spurious conserved quantity appears; we use it to effect a radical transformation of the model into one whose spacetime configurations are equivalent to the two-dimensional states of an exactly solvable statistical mechanics model, the symmetric eight-vertex model with parameters restricted to a disorder variety. We comment on the implications of this identification for the original lattice gas.     

Bifurcations of a lattice gas flow under external forcing

We study the behavior of a Frisch-Hasslacher-Pomeau lattice gas automaton under the effect of a spatially periodic forcing. It is shown that the lattice gas dynamics reproduces the steady-state features of the bifurcation pattern predicted by a properly truncated model of the Navier-Stokes equations. In addition, we show that the dynamical evolution of the instabilities driving the bifurcation can be modeled by supplementing the truncated Navier-Stokes equation with a random force chosen on the basis of the automaton noise.