The cellular automata are models that discretize the space time and define automata in each point (cell) of the network that act in function of what their neighbors do (automata because they have a defined form of reacting). An example is a hexagonal structure:

image

Model D2Q7 (two dimensions and 7 elements per cell - 6 sides and 1 center)

In the case that it is applied to a particulate gas, each node may or may not contain (states 0 and 1) a particle that can only have velocities with the directions that links between cells.

In the simulation with models like cellular automata there are two phases:

- cell acts on the others
- cell processes actions of the environment

In the special case of modeling a gas the first step corresponds to streaming while the second to collision (collision).

The mathematical description is performed by the particle distribution function f(\vec{x},\vec{v},t) where \vec{x} is Position, \vec{v} the speed and t the time. As in this case there are only discrete speeds \vec{e}_i tend to indicate the distribution function as a set of functions f_i such that

equation=8466

In the case of the discretization in the LBM models we work not with functions of the speed if not with discrete components. In this way the i component is defined by:

equation

where w_i is the relative weight.

The dispersion is associated with the relaxation time \tau which in the hydrodynamic approximation is reflected by the viscosity

equation

with \rho the density, \Delta x the length of the cell and \Delta t simulation time step.

By discretizing we assume that the particles move with a mean velocity c that is associated to the mesh of the network of cells of long \Delta x and time \Delta t using

equation

El modelo D2Q9 es un modelo bidimensional (D2) en que se se conecta el nodo (punto central) en nodos a lo largo de los ejes cartesianos\\n\\nen el origen\\n\\n

$\vec{e}_0=(0,0)$

\\n\\nen las esquinas\\n\\n

$\vec{e}_1=(1,0)$

(E),\\n

$\vec{e}_2=(0,1)$

(N), \\n

$\vec{e}_3=(-1,0)$

(W) y \\n

$\vec{e}_4=(0,-1)$

(S)\\n\\ny en las diagonales\\n\\n

$\vec{e}_5=(1,1)$

(NE), \\n

$\vec{e}_6=(-1,1)$

(SE), \\n

$\vec{e}_7=(-1,-1)$

(SW) y \\n

$\vec{e}_8=(1,-1)$

(NW)

lo que se representa en la siguiente gr fica:

image

The D3Q15 model is a two-dimensional model (D3) in which the node (central point) is connected in nodes along Cartesian axes\\n\\n

$(1,0,0), (-1,0,0), (0,1,0), (0,-1,0), (0,0,1) (0,0,-1)$

\\n\\nand in the corners of the bucket\\n\\n

$(1,0,1), (-1,0,1), (0,1,1) , (0,-1,1), (1,0,-1), (-1,0,-1), (0,1,-1) , (0,-1,-1)$

which is represented in the following graph:

image

It is easy to build models of type D3Q19 (including side edge halves) or D3Q27 (all possible points).