格子气自动机

✍ dations ◷ 2025-04-03 12:38:14 #计算流体力学,计算物理学,细胞自动机

格子气自动机或格子气细胞自动机是一种用来模拟流体流动的细胞自动机。它是格子玻尔兹曼方法(英语:Lattice Boltzmann methods)的前身。通过格子气自动机的部分模型可以推导得出宏观的纳维-斯托克斯方程。20世纪90年代初,格子玻尔兹曼方法(英语:Lattice Boltzmann methods)开始兴起,对格子气自动机方法的研究则趋缓了。

细胞自动机模型中,晶格格点具有有限的离散状态。在格子气自动机模型中,这些离散状态用以表示具有不同速度的粒子。流体模拟在离散的时间序列中演进,在每个时间步长后,格点的状态可由该时间步长前此格点与相邻格点的状态确定。

每个格点以一系列布尔值来描述给定位置上各速度方向是否有粒子存在。

每个时间步长的状态演进都分为两个阶段,即迁移和碰撞。

迁移:除非发生碰撞,每个粒子将沿它的速度方向向距离最近的格点移动,并维持该速度。泡利不相容原理还要求同一格点同一速度方向的粒子最多只能有一个。

碰撞:碰撞规则用于确定多个粒子到达同一格点时粒子的行为方式,并应满足质量守恒与动量守恒。 请注意,泡利不相容原理并不影响同一格点的两个粒子以相反方向迁移。

在1973年和1976年发表的论文中,Hardy,Pomeau和de Pazzis提出了首个格子气自动机模型,后根据三位作者的名字命名为HPP模型。HPP模型是一个二维的流体粒子相互作用模型,其晶格为正方形,粒子以单位速度独立运动。粒子可以移动至相邻的四个晶格而不能沿对角方向移动。

发生对头碰撞时,如一个向左移动的粒子与一个向右移动的粒子到达同一个格点,则两个粒子的速度将分别旋转90°离开该格点。

HPP模型缺乏旋转不变性,从而使模型具有高度各向异性。 例如,HPP模型产生的涡流是正方形的。

1986年,乌里尔·弗里施(英语:Uriel Frisch),布罗斯·哈斯拉彻(英语:Brosl Hasslacher)和伊夫·波莫(英语:Yves Pomeau)在论文中首次提出了以规则六边形为晶格的二维格子气自动机模型,即FHP模型。在基本FHP模型(又称FHP-I模型)中,流体粒子具有指向相邻的6个晶格的6种离散速度,碰撞规则包括对头二体碰撞和对称三体碰撞;在此基础上,FHP-II模型引入了代表粒子静止的第七种速度,其碰撞规则相应增加了含静止粒子的碰撞及非对称三体碰撞;FHP-III模型更进一步引入了对称四体碰撞,至此包括了能满足质量和动量守恒的所有碰撞规则。 碰撞次数的增加会增大雷诺数,因此与FHP-I模型相比,FHP-II和FHP-III模型可以模拟粘性更小的流体。

与HPP模型不同,FHP模型中满足质量、动量守恒的碰撞规则具有不确定性,部分碰撞类型会产生两种结果,此情况下会随机选择其中一种,相应粒子也将依照两种碰撞结果之一迁移。由于纯计算方法不能生成随机数,因此通常采用伪随机算法。

FHP模型比HPP模型具有更好的对称性,但这一结论并不显然。

对于三维晶格,唯一能密铺三维空间的规则多面体是立方体,而唯一能密铺三维空间、具有足够大对称性(否则模型将具有和HPP模型类似的缺陷)的多胞形组合是十二面体和二十面体。因此,模拟三维流场需要增加晶格的维数,例如1986年D'Humières,Lallemand和Frisch提出的面心超立方体模型。

可以通过每个格点的粒子数来计算格点的密度,以对应单位速度加权求和即可获得该位置的动量。此计算过程会受到大量噪声的影响,在模拟中应在较大区域内取均值以获得更合理的结果。通常使用系综平均(英语:Ensemble average)进一步降低统计噪声。

格子气自动机的主要优势是只涉及布尔运算因此避免了舍入误差,且作为细胞自动机具有天然的并行性。

其不足之处在于缺乏伽利略不变性(英语:Galilean invariance)和随机因素造成的统计噪声。另外,格子气自动机难以直接扩展至三维情形,需要引入更高维度来维持足够的对称性。

相关

  • 什切青旧城什切青旧城(波兰语:Stargard Szczeciński;卡舒比语、波美拉尼亚语:Stôrgard;德语:Stargard in Pommern)是波兰西北部波美拉尼亚地区的一个城镇,斯塔加德县首府,位于伊纳河畔,2005年
  • 手帕手绢,亦称手帕、手捏子、手巾、面帕,是一种生活用品。
  • 大兴校区北京大学软件与微电子学院正式成立于2002年,是中华人民共和国教育部和中华人民共和国国家发展和改革委员会以教高6号文正式批准北京大学等35所高等学校试办示范性软件学院之
  • 约瑟夫·泰勒约瑟夫·胡顿·泰勒(英语:Joseph Hooton Taylor,1941年3月29日-),美国物理学家。他和拉塞尔·赫尔斯共同发现史上第一个位于双星系统脉冲星PSR B1913+16,并通过对其深入研究首次发
  • 卡方分布卡方分布(chi-square distribution, ²-distribution,或写作²分布)是概率论与统计学中常用的一种概率分布。k个独立的标准正态分布变量的平方和服从自由度为k的卡方分布。卡方
  • 蕾蜜·海德利蕾蜜·海德利是福克斯电视网出品的电视剧《豪斯医生》中的重要角色之一,扮演者为奥利维亚·王尔德。她是医学博士,曾是格里高利·豪斯医生的第二批诊断小组的成员。该人物的绰
  • 西格丽德·温塞特西格丽德·温塞特(挪威语:Sigrid Undset,1882年5月20日-1949年6月10日)是一位挪威语小说家,1928年诺贝尔文学奖获得者。值得一提的是,温塞特获奖后,她将奖金全部捐赠给社会福利机关
  • 马尼乌斯·库里乌斯·登塔图斯马尼乌斯·库里乌斯·登塔图斯(Manius Curius Dentatus,?~前270年)古罗马军事活动家和政治家,曾4次当选为执政官(前290年,前284年,前275年和前274年),一次当选为监察官(前272年)。他因领
  • 北寨墓群北寨墓群位于山东省临沂市沂南县县城西北3公里的北寨村内,是中国现已发现的、保存最完整的一座东汉晚期(2世纪)的画像石墓。有M1、M2号两座墓。一号墓葬有夫妻二人。左妻右夫,根
  • 认证互连设计者认证互连设计者(Certified Interconnect Designer)简称CID,是国际电子工业联接协会(英语:IPC (electronics))(IPC)针对有经验印刷电路板设计专家提出的认证。CID+是CID的进阶版本。