Physics:HPP model

Small scale demonstration of the square lattice HPP model (click the picture for animated version).

The Hardy–Pomeau–Pazzis (HPP) model is a fundamental lattice gas automaton for the simulation of gases and liquids. It was a precursor to the lattice Boltzmann methods. From lattice gas automata, it is possible to derive the macroscopic Navier-Stokes equations. Interest in lattice gas automaton methods levelled off in the early 1990s, due to rising interest in the lattice Boltzmann methods.

It was first introduced in papers published in 1973 and 1976 by Jean Hardy, Yves Pomeau and Olivier de Pazzis, whose initials give the model its name. The model can be used as a simple model for both the movement of gases and fluid.

Model

Visual examples of the four rules governing the HPP Model.

In this model, the lattice takes the form of a two-dimensional square grid, with particles capable of moving to any of the four adjacent grid points which share a common edge, and particles cannot move diagonally. This means each grid point can only have one of sixteen possible interactions.

• Particles exist only on the grid points, never on the edges or surface of the lattice.
• Each particle has an associated direction (from one grid point to another immediately adjacent grid point).
• Each lattice grid cell can only contain a maximum of one particle for each direction, i.e., contain a total of between zero and four particles.

An implementation of the HPP model. A 300x200 pixel grid was used with a torus like geometry, meaning particles leaving the grid from the left enter form the right again, particles leaving from the top enter from the bottom, etc. Particle density is displayed by a white pixel being densest with 4 particles and black the least dense with 0 particles. The simulation starts with a high pressure 100x100 square in the middle.

The following rules also govern the model:

1. A single particle moves in a fixed direction until it experiences a collision.
2. Two particles experiencing a head-on collision are deflected perpendicularly.
3. Two particles experience a collision which isn't head-on simply pass through each other and continue in the same direction.
4. Optionally, when a particles collides with the edges of a lattice it can rebound.

The HPP models follows a two-stage update process.

Collision step

In this step, the above rules 2., 3., and 4. are checked and applied if any collisions have occurred. This results in head-on collision particles changing direction, pass-through collisions continuing unchanged, or non-colliding particles simple remaining the same.

Transport step

The second step consists of each particle moving one lattice step in the direction they are currently travelling, which could have been changed by the above Collision Step.

Formal definition

The model operates on an infinite two-dimensional square lattice, where four unit vectors are associated with the following numbers: [math]\displaystyle{ 1 \mapsto \left( \begin{array}{c} 1 \\ 0 \end{array} \right), 2 \mapsto \left( \begin{array}{c} 0 \\ 1 \end{array} \right), 3 \mapsto \left( \begin{array}{c} -1 \\ 0 \end{array} \right), 4 \mapsto \left( \begin{array}{c} 0 \\ -1 \end{array} \right) }[/math].

Let [math]\displaystyle{ X }[/math] be an allowed configuration. The function [math]\displaystyle{ \sigma_X }[/math] checks for the existence of a particle with a certain velocity, while [math]\displaystyle{ \bar{\sigma}_X }[/math] does the opposite:

[math]\displaystyle{ \sigma_X(p,q,i) = \begin{cases} 1 \text{ if the site } (p,q) \text{ is occupied by a particle of } X \text{ of velocity } i, \\ 0 \text{ otherwise.} \end{cases} }[/math]

[math]\displaystyle{ \bar{\sigma}_X(p,q,i) = 1 - \sigma_X(p,q,i) }[/math]

The successor [math]\displaystyle{ T(X) }[/math] of the configuration [math]\displaystyle{ X }[/math] can be calculated using the formulas form the original paper:[math]\displaystyle{ \begin{array}{lcl} \sigma_{T(X)}(p,q,1) & = & \sigma_X(p-1,q,1)-\Psi_X(p,q) \\ \sigma_{T(X)}(p,q,2) & = & \sigma_X(p,q-1,2)+\Psi_X(p,q) \\ \sigma_{T(X)}(p,q,3) & = & \sigma_X(p+1,q,3)-\Psi_X(p,q) \\ \sigma_{T(X)}(p,q,4) & = & \sigma_X(p,q+1,4)+\Psi_X(p,q) \end{array} }[/math][math]\displaystyle{ \begin{array}{lcl} \Psi_X(p,q) & = & \psi^1_X(p,q) - \psi^2_X(p,q) \\ \psi^1_X(p,q) & = & \sigma_X(p-1,q,1)\sigma_X(p+1,q,3)\bar{\sigma}_X(p,q-1,2)\bar{\sigma}_X(p,q+1,4)\\ \psi^2_X(p,q) & = & \bar{\sigma}_X(p-1,q,1)\bar{\sigma}_X(p+1,q,3)\sigma_X(p,q-1,2)\sigma_X(p,q+1,4) \end{array} }[/math]

Shortcomings

The model is badly flawed, as momentum is always conserved in both the horizontal and vertical lanes. No energy is ever removed from the model, either by collisions or movement, so it will continue indefinitely.

The HPP model lacked rotational invariance, which made the model highly anisotropic. This means for example, that the vortices produced by the HPP model are square-shaped.

References

• Sauro Succi (2001). The Lattice Boltzmann Equation, for fluid dynamics and beyond. Oxford Science Publications. ISBN 0-19-850398-9.  (Chapter 2 on Lattice gas Cellular Automata)
• Neil Gershenfeld (1998). The Nature of Mathematical Modeling. Cambridge University Press. ISBN 978-0521570954.