格子気体法

格子気体法について

格子気体法（Lattice Gas Cellular Automata）は、流体の動きを数値的にシミュレーションするための手法です。この方法はセル・オートマトンに基づいており、流体や計算空間を格子と粒子に分けることで離散化されています。具体的には、粒子に衝突と並進という2つの主要な運算を行い、それを繰り返すことで流体の時間的な変化を模擬します。この手法は、様々なモデルやバージョンが存在し、それぞれの格子の形状や衝突ルールが異なる点が特徴です。

モデルの概要

HPPモデル

1973年、Hardy、Pomeau、de Pazzisによって開発されたHPPモデルは、格子気体法の基盤となるモデルです。このモデルの主な特徴は、計算空間を正方形の格子で離散化し、各格子に最大4つの粒子を配置可能とすることです。

各粒子は、ノイマン近傍の方向へ向かう限られた速度を持ち、同一の格子内では異なる速度が設定されています。衝突ルールは、同じ格子に存在する2つの粒子が平行な速度を持つ場合にのみ、その速度が90度回転するというものです。また、並進ルールでは、全ての粒子が1時間あたりその速度に従って隣接する格子に移動します。これによって、運動量と質量の保存則が守られます。しかし、HPPモデルには流体の応力テンソルの等方性が満たされないなどの問題が存在します。

FHPモデル

1986年にFrisch、ハシュラヒャー、Pomeauによって開発されたFHPモデルは、HPPモデルの改良版として知られています。このモデルでは、HPPモデルの持つテンソルの等方性に関する問題が解消され、現在では格子気体法の主流となっています。

FHPモデルの特徴として、格子が正六角形で構成されている点が挙げられます。各格子には最大で7つの粒子が存在でき、粒子の速度は六角形の中心から各頂点へ向かう6種類に加え、静止している状態も考慮されています。また、衝突ルールは運動量保存則を基にしており、2体衝突だけでなく、3～5体の衝突も規定されています。このモデルによるシミュレーションは、整数演算や真偽値演算を使用することで、高速かつ安定した計算が可能となっています。

まとめ

格子気体法は、流体力学のシミュレーションにおいて多様なアプローチを提供する重要な手法です。HPPモデルからFHPモデルへの進化を通じて、より高度な計算が実現され、多くの分野でその応用が期待されています。セル・オートマトンや格子ボルツマン法と並んで、流体の動きを理解し、模擬するための基盤が築かれています。

もう一度検索