fbpx
维基百科

格子气自动机

十月 19, 2023

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

气体流仿真的HPP模型。像素灰度正比于对应格点的气体粒子密度（0～4）。气体被一层以黄色标注的反射元胞包围以形成封闭流场。

基本思想 编辑

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

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

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

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

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

HPP模型 编辑

 
HPP模型的小规模演示

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

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

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

FHP模型 编辑

 
三种FHP模型的碰撞规则

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

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

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

三维模型 编辑

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

宏观动力学 编辑

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

优势与不足 编辑

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

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

备注 编辑

  1. ^ Succi, section 2.3 describes the process
  2. ^ Succi, section 2.6
  3. ^ Buick, section 3.4
  4. ^ Wolfram, Stephen, A New Kind of Science, Wolfram Media: 459–464, 2002, ISBN 1-57955-008-8 .
  5. ^ Buick, section 3.2.1
  6. ^ Succi, footnote p. 22
  7. ^ Buick, section 3.2.2
  8. ^ Wolf-Gladrow 3.2.6, figure 3.2.3
  9. ^ Wolf-Gladrow 3.2.1
  10. ^ Succi, footnote p. 23
  11. ^ 11.0 11.1 Wolf-Gladrow, sections 3.4 - 3.5
  12. ^ Buick, section 3.5.1
  13. ^ Buick, section 3.8
  14. ^ Succi, section 2.4
  15. ^ Succi, section 2.5

参考文献 编辑

  • Sauro Succi. The Lattice Boltzmann Equation, for fluid dynamics and beyond. Oxford Science Publications. 2001. ISBN 0-19-850398-9.  (Chapter 2 is about lattice gas Cellular Automata)
  • James Maxwell Buick (1997). Lattice Boltzmann Methods in Interfacial Wave Modelling. PhD Thesis, University of Edinburgh. (Chapter 3 is about the lattice gas model.) () 2008-11-13
  • Dieter A. Wolf-Gladrow. Lattice-Gas Cellular Automata and Lattice Boltzmann Models. Springer. 2000. ISBN 3-540-66973-6. 

十月 19, 2023
格子气自动机, 或格子气细胞自动机是一种用来模拟流体流动的细胞自动机, 它是格子玻尔兹曼方法, 英语, lattice, boltzmann, methods, 的前身, 通过的部分模型可以推导得出宏观的纳维, 斯托克斯方程, 20世纪90年代初, 格子玻尔兹曼方法, 英语, lattice, boltzmann, methods, 开始兴起, 对方法的研究则趋缓了, 气体流仿真的hpp模型, 像素灰度正比于对应格点的气体粒子密度, 气体被一层以黄色标注的反射元胞包围以形成封闭流场, 目录, 基本思想, hpp模型. 格子气自动机或格子气细胞自动机是一种用来模拟流体流动的细胞自动机 它是格子玻尔兹曼方法 英语 Lattice Boltzmann methods 的前身 通过格子气自动机的部分模型可以推导得出宏观的纳维 斯托克斯方程 1 20世纪90年代初 格子玻尔兹曼方法 英语 Lattice Boltzmann methods 开始兴起 2 对格子气自动机方法的研究则趋缓了 气体流仿真的HPP模型 像素灰度正比于对应格点的气体粒子密度 0 4 气体被一层以黄色标注的反射元胞包围以形成封闭流场 目录 1 基本思想 2 HPP模型 3 FHP模型 4 三维模型 5 宏观动力学 6 优势与不足 7 备注 8 参考文献基本思想 编辑细胞自动机模型中 晶格格点具有有限的离散状态 在格子气自动机模型中 这些离散状态用以表示具有不同速度的粒子 流体模拟在离散的时间序列中演进 在每个时间步长后 格点的状态可由该时间步长前此格点与相邻格点的状态确定 每个格点以一系列布尔值来描述给定位置上各速度方向是否有粒子存在 每个时间步长的状态演进都分为两个阶段 即迁移和碰撞 3 迁移 除非发生碰撞 每个粒子将沿它的速度方向向距离最近的格点移动 并维持该速度 泡利不相容原理还要求同一格点同一速度方向的粒子最多只能有一个 碰撞 碰撞规则用于确定多个粒子到达同一格点时粒子的行为方式 并应满足质量守恒与动量守恒 4 请注意 泡利不相容原理并不影响同一格点的两个粒子以相反方向迁移 HPP模型 编辑 nbsp HPP模型的小规模演示在1973年和1976年发表的论文中 Hardy Pomeau和de Pazzis提出了首个格子气自动机模型 后根据三位作者的名字命名为HPP模型 HPP模型是一个二维的流体粒子相互作用模型 其晶格为正方形 粒子以单位速度独立运动 粒子可以移动至相邻的四个晶格而不能沿对角方向移动 发生对头碰撞时 如一个向左移动的粒子与一个向右移动的粒子到达同一个格点 则两个粒子的速度将分别旋转90 离开该格点 5 HPP模型缺乏旋转不变性 从而使模型具有高度各向异性 例如 HPP模型产生的涡流是正方形的 6 FHP模型 编辑 nbsp 三种FHP模型的碰撞规则1986年 乌里尔 弗里施 英语 Uriel Frisch 布罗斯 哈斯拉彻和伊夫 波莫 英语 Yves Pomeau 在论文中首次提出了以规则六边形为晶格的二维格子气自动机模型 即FHP模型 在基本FHP模型 又称FHP I模型 中 流体粒子具有指向相邻的6个晶格的6种离散速度 碰撞规则包括对头二体碰撞和对称三体碰撞 在此基础上 FHP II模型引入了代表粒子静止的第七种速度 其碰撞规则相应增加了含静止粒子的碰撞及非对称三体碰撞 FHP III模型更进一步引入了对称四体碰撞 至此包括了能满足质量和动量守恒的所有碰撞规则 7 碰撞次数的增加会增大雷诺数 因此与FHP I模型相比 FHP II和FHP III模型可以模拟粘性更小的流体 8 与HPP模型不同 FHP模型中满足质量 动量守恒的碰撞规则具有不确定性 部分碰撞类型会产生两种结果 此情况下会随机选择其中一种 相应粒子也将依照两种碰撞结果之一迁移 由于纯计算方法不能生成随机数 因此通常采用伪随机算法 9 FHP模型比HPP模型具有更好的对称性 但这一结论并不显然 10 三维模型 编辑对于三维晶格 唯一能密铺三维空间的规则多面体是立方体 而唯一能密铺三维空间 具有足够大对称性 否则模型将具有和HPP模型类似的缺陷 的多胞形组合是十二面体和二十面体 因此 模拟三维流场需要增加晶格的维数 例如1986年D Humieres Lallemand和Frisch提出的面心超立方体模型 11 宏观动力学 编辑可以通过每个格点的粒子数来计算格点的密度 以对应单位速度加权求和即可获得该位置的动量 12 此计算过程会受到大量噪声的影响 在模拟中应在较大区域内取均值以获得更合理的结果 通常使用系综平均 英语 Ensemble average 进一步降低统计噪声 13 优势与不足 编辑格子气自动机的主要优势是只涉及布尔运算因此避免了舍入误差 且作为细胞自动机具有天然的并行性 14 其不足之处在于缺乏伽利略不变性 英语 Galilean invariance 和随机因素造成的统计噪声 15 另外 格子气自动机难以直接扩展至三维情形 需要引入更高维度来维持足够的对称性 11 备注 编辑 Succi section 2 3 describes the process Succi section 2 6 Buick section 3 4 Wolfram Stephen A New Kind of Science Wolfram Media 459 464 2002 ISBN 1 57955 008 8 Buick section 3 2 1 Succi footnote p 22 Buick section 3 2 2 Wolf Gladrow 3 2 6 figure 3 2 3 Wolf Gladrow 3 2 1 Succi footnote p 23 11 0 11 1 Wolf Gladrow sections 3 4 3 5 Buick section 3 5 1 Buick section 3 8 Succi section 2 4 Succi section 2 5参考文献 编辑Sauro Succi The Lattice Boltzmann Equation for fluid dynamics and beyond Oxford Science Publications 2001 ISBN 0 19 850398 9 Chapter 2 is about lattice gas Cellular Automata James Maxwell Buick 1997 Lattice Boltzmann Methods in Interfacial Wave Modelling PhD Thesis University of Edinburgh Chapter 3 is about the lattice gas model archive org 2008 11 13 Dieter A Wolf Gladrow Lattice Gas Cellular Automata and Lattice Boltzmann Models Springer 2000 ISBN 3 540 66973 6 取自 https zh wikipedia org w index php title 格子气自动机 amp oldid 75886004, 维基百科，wiki，书籍，书籍，图书馆，

文章

，阅读，下载，免费，免费下载，mp3，视频，mp4，3gp， jpg，jpeg，gif，png，图片，音乐，歌曲，电影，书籍，游戏，游戏。