Modelo HPP

Summary

El modelo Hardy-Pomeau-Pazzis (HPP) es un autómata de gas en celosía fundamental para la simulación de gases y líquidos. Fue precursor de los métodos de Boltzmann reticular. A partir de los autómatas de gas en red, es posible derivar las ecuaciones de Navier-Stokes macroscópicas.[1]​ El interés por los métodos de autómatas de gas en red se estabilizó a principios de la década de 1990, debido al creciente interés por los métodos de Boltzmann en red.[2]

Demostración a pequeña escala del modelo HPP de red cuadrada (haga clic en la imagen para ver la versión animada).

Se introdujo por primera vez en artículos publicados en 1973 y 1976 por Jean Hardy, Yves Pomeau y Olivier de Pazzis,[3]​ cuyas iniciales dan nombre al modelo. El modelo puede utilizarse como modelo simple tanto para el movimiento de gases como de fluidos.[4]

Modelo

editar
 
Ejemplos visuales de las cuatro reglas que rigen el modelo HPP.

En este modelo, la red adopta la forma de una cuadrícula bidimensional, con partículas capaces de moverse a cualquiera de los cuatro puntos adyacentes de la cuadrícula que comparten un borde común, y las partículas no pueden moverse en diagonal. Esto significa que cada punto de la cuadrícula solo puede tener una de las dieciséis interacciones posibles.

  • Las partículas solo existen en los puntos de la red, nunca en los bordes o en la superficie de la red.
  • Cada partícula tiene una dirección asociada (de un punto de la red a otro punto inmediatamente adyacente).
  • Cada células de red de la red solo puede contener un máximo de una partícula por cada dirección, es decir, contiene un total de entre cero y cuatro partículas.
 
Una implementación del modelo HPP. Se utilizó una red de 300 x 200 píxeles con una geometría similar a un toro, lo que significa que las partículas que salen de la red por la izquierda vuelven a entrar por la derecha, las que salen por la parte superior entran por la parte inferior, etc. La densidad de las partículas se muestra mediante píxeles blancos, que son los más densos con 4 partículas, y negros, que son los menos densos con 0 partículas. La simulación comienza con un cuadrado de alta presión de 100 x 100 en el centro.

El modelo también se rige por las siguientes reglas:

  1. Una sola partícula se mueve en una dirección fija hasta que experimenta una colisión.
  2. Dos partículas que experimentan una colisión frontal se desvían perpendicularmente.
  3. Dos partículas que experimentan una colisión que no es frontal simplemente se atraviesan y continúan en la misma dirección.
  4. Opcionalmente, cuando una partícula colisiona con los bordes de una red, puede rebotar.

El modelo HPP sigue un proceso de actualización en dos etapas.

Paso de colisión

editar

En este paso, se comprueban las reglas 2., 3. y 4. anteriores y se aplican si se ha producido alguna colisión. Esto da como resultado que las partículas en colisión frontal cambien de dirección, las colisiones de paso continúen sin cambios o las partículas que no colisionan simplemente permanezcan igual.

Paso de transporte

editar

El segundo paso consiste en que cada partícula se mueva un paso de la red en la dirección en la que se está desplazando, que podría haber cambiado en el paso de colisión anterior.

Definición formal

editar

El modelo opera en una red cuadrada bidimensional infinita, en la que cuatro vectores unitarios están asociados a los siguientes números:  .[5]

Sea   una configuración permitida. La función   comprueba la existencia de una partícula con una velocidad determinada, mientras que   hace lo contrario:  [5] [5]

El sucesor   de la configuración   se puede calcular utilizando las fórmulas del artículo original[5]

 

 

Deficiencias

editar

El modelo presenta graves deficiencias, ya que el momento siempre se conserva tanto en los carriles horizontales como en los verticales. Nunca se elimina energía del modelo, ni por colisiones ni por movimiento, por lo que continuará indefinidamente.

El modelo HPP carecía de invarianza rotacional, lo que lo hacía altamente anisotrópico. Esto significa, por ejemplo, que los vórtices producidos por el modelo HPP son de forma cuadrada.[6]

Referencias

editar
  1. Succi, sección 2.3 describe el proceso
  2. Succi, sección 2.6
  3. Rothman, Daniel H.; Zaleski, Stiphane, eds. (1997), «A simple model of fluid mechanics», Lattice-Gas Cellular Automata: Simple Models of Complex Hydrodynamics, Collection Alea-Saclay: Monographs and Texts in Statistical Physics (Cambridge: Cambridge University Press): 1-11, ISBN 978-0-521-60760-5, doi:10.1017/CBO9780511524714.002, consultado el 24 de agosto de 2022 .
  4. Gershenfeld, pp. 103
  5. a b c d Hardy, J.; Pomeau, Y.; de Pazzis, O. (30 de julio de 1973). «Time Evolution of a Two-Dimensional Classical Lattice System». Physical Review Letters 31 (5): 276-279. doi:10.1103/PhysRevLett.31.276. 
  6. Succi, nota al pie de página 22

Bibliografía

editar
  • 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.