格子气模型
格子玻尔兹曼方法及其在大气湍流研究中的应用_图文(精)
第22卷第3期2007年3月地球科学进展A DVAN CE S I N E AR T H S C I E N C EV o l.22N o.3M a r.,2007文章编号:1001-8166(200703--12格子玻尔兹曼方法及其在大气湍流研究中的应用*程雪玲,胡非,赵松年,姜金华(中国科学院大气物理研究所大气边界层物理和大气化学国家重点实验室,北京100029摘要:文章的目的是对格子玻尔兹曼方法进行系统的介绍,格子玻尔兹曼方法(L a tti ce B o lt z m a nn M e t hod的出现直接来源于20世纪60年代的元胞自动机(C e ll u l a r A u t om a t a思想,而这一方法用于解决流动现象时,又可以追溯到19世纪的分子运动论,求解的是B o lt z m a nn提出的玻尔兹曼输运方程,因此将这一方法称为格子玻尔兹曼方法,之前也被称为格子气自动机(L a tti c e G a s A u t om a-t o n。
该方法多用于研究复杂现象,如材料晶体凝聚时的生长过程、城市土地利用的演化等方面。
在20世纪70年代由H a r dy、P om e au和P a zz i s建立了第一个用于研究流体运动的格子气自动机,此后,这一方法被广泛用来模拟各种流动问题,诸如二相流、孔隙介质中的渗流等,并根据这一方法开发了相应的商业软件P ow e r F l o w。
同时,格子玻尔兹曼方法由于其在微观水平描述运动的特点,成为研究湍流的一个很好的数值计算工具,特别是用其进行直接数值模拟(D N S计算,成为继传统的差分法、有限体积法和谱方法之后的又一有力的手段。
而作为大气运动的一个主要现象的大气湍流,比普通湍流更加复杂,在这里着重介绍了大气湍流的特点和应用格子玻尔兹曼方法模拟湍流的发展过程。
关键词:格子玻尔兹曼;元胞自动机;格子气;直接数值模拟;大气湍流中图分类号:P425.2文献标识码:A1引言元胞自动机是空间和时间都离散、物理参量只取有限数值集的物理系统的理想化模型。
偏微分方程求解的一种新颖方法——格子Boltzmann模型
７ ６ 邻 节 点
大 学 数 学
第２ ７卷
） 撞 ， 一 个 节 点 上 从 相 邻 节 点 运 动 来 的 粒 子 发 生 碰 撞 ， 据 质 量 、 量 和 能 量 守 恒 规 则 改 碰 在 根 动
其 中 ｒ 松弛 时间 尺度 ， 制达 到平衡 的速 度 （ 是 控 可根据 需 要 进行 设 置 ） 由于稳 定 性 的原 因 ， 过 实 际测 ， 经 算 ｒ必须 大于 １ｅ ／．
事 实 上 不 同 的 网 格 剖 分 有 着 不 同 的平 衡 分 布 函数 ， Ｂ 建 立 模 型 的 核 心 问 题 就 是 根 据 不 同 的 网 格 Ｌ Ｍ
［ 键 词 ］ 格 子 Ｂ ｌ ｍａ ｎ方 法 ； 衡 态 分 布 函 数 ； Ｑ 关 ｏｔ ｎ ｚ 平 Ｄ２ ９模 型 ； ａ ｉ — ｔｋ ｓ 程 ； 流一 扩 散 方 程 Ｎ ｖｅ Ｓｏ ｅ 方 ｒ 对 ［ 图 分 类 号 ］ Ｏ２ １ ８ 中 ４ ．２ ［ 献标识码］Ａ 文 ［ 章 编 号 ］ １７ —４ ４ ２ １ ） ３０ ７ —８ 文 ６ ２１ ５ （０ １ ０ —０ ５０
在 低 Ｍａｈ 马赫 ） 的假 设下 （ ｌ ｃ）其 中粒子平 衡态 分布 函数 为 ｃ（ 数 Ｉ ， Ｕ《
～ Ｐ
［ ＋ 一 ］ ＋
且
Ｃ ＝ｃ ４ ／ ￣，
。 ／， 一４ ９
１ Ｕ 一 ３ 一（２ 一 ４ １ ９， ∞ 一 ６ ７ ８ １ ３ ， — ／ ５ 一∞ 一∞ — ／ ６
第２ ７卷 第 ３期
基于元胞自动机模型的人员疏散行为模拟
基于元胞自动机模型的人员疏散行为模拟郭良杰;赵云胜【摘要】通过引入静态场、动态场及其他参数,建立了二维元胞自动机模型,并利用MATLAB编写模拟软件,可实现对人员疏散过程中的环境熟悉度、从众行为、摩擦阻碍作用、惯性行为、拥挤跌倒行为和竞争行为的模拟.通过实例仿真模拟结果表明:不同的期望速度下会出现“欲速则不达”的现象,同时出口处会产生不同程度的人员聚集;出口处障碍物相对出口纵向放置比横向放置更有利于人员疏散;在疏散环境陌生,或紧急情况下对环境判断能力降低时,适当的从众行为利于最优疏散路径信息的传递,从而有利于人员疏散,但是从众行为过于严重则易造成出口利用率降低或利用不平衡.【期刊名称】《安全与环境工程》【年(卷),期】2014(021)004【总页数】6页(P101-106)【关键词】元胞自动机;地面场模型;人员疏散行为;仿真模拟;从众行为【作者】郭良杰;赵云胜【作者单位】中国地质大学工程学院,湖北武汉430074;中国地质大学工程学院,湖北武汉430074【正文语种】中文【中图分类】X913.4近几年我国灾害频发，如近期发生的雅安地震、吉林宝源丰禽业公司氨气泄漏爆炸火灾以及厦门BRT公交车火灾等。
灾害发生后，人员聚集场所的紧急安全疏散再次成为人们关心的问题和研究的热点。
国内外针对此方面的研究主要集中在对紧急条件下疏散时间和人员逃生行为特征的研究，但由于实际数据的缺乏和实验准确性的不足，建立人员疏散模型并利用计算机进行仿真模拟已成为研究人员疏散行为的主要手段之一［1］。
目前较常用的人员疏散模型有连续型（社会力模型等）和离散型（格子气模型、元胞自动机模型等）之分。
本文基于元胞自动机模型对人员疏散过程中的现象和疏散行为进行了定性研究。
1 元胞自动机模型建立1.1 元胞自动机简介元胞自动机（Cellular Automata，CA）是由大量简单一致的个体通过局部联系组成的离散、分散及空间可扩展系统，它是在均匀一致的网格上由有限状态的变量（元胞）构成的动力系统［1］。
基于格子Boltzrnann方法的回采工作面瓦斯渗流模拟
基于格子Boltzrnann方法的回采工作面瓦斯渗流模拟摘要：工作面回采会导致其前方煤体中的应力重新分布，在工作面前方形成减压区、增压区和稳压区，致使回采工作面前方煤体中的渗透率沿煤层走向呈现不均匀分布，因此，不能再将煤体视为均匀介质。
鉴于这一工程事实，采用格子B oltz m an n 方法，结合算例，研究了煤体渗透率沿走向呈现不均匀分布情形下的瓦斯渗流规律，为回采工作面瓦斯灾害防治提供依据。
研究表明: 采用格子B oltz m an n 方法可以有效模拟回采工作面瓦斯渗流，能够生动再现瓦斯流向回采工作面这一过程。
关键词：回采工作面瓦斯渗流渗透率格子B oltz m an n 压力分布瓦斯涌出量研究回采工作面瓦斯渗流规律对防治煤与瓦斯突出、预测瓦斯涌出量以及合理设计通风具有重要指导意义。
关于瓦斯在采场中运移的研究，目前的观点认为瓦斯在工作面的运移为紊流运动，瓦斯在采空区的运移为气体在多孔介质中的渗流运动，通过构建偏微分方程组，应用传统的计算流体力学、单元法原理和网络解算技术求解[1]。
格子玻尔兹曼方法(lattice Boltzmannmethod，LBM)，在对物理系统细观运动描述的基础上，建立了直接模拟物理系统演化的格子气模型。
物理系统的宏观物理量直接由模型的演化方程来计算，从而建立了宏观物理量与细观结构运动演化的相互联系。
本文采用L B M 方法，对回采工作面前方煤体中瓦斯渗流进行了研究，考虑到实际工程中回采工作会使煤体应力重新分布，导致工作面前方煤体的渗透率沿煤层走向分布不均匀，研究过程中将煤体渗透率视为随地应力变化的参数，研究结果为回采工作面瓦斯灾害防治提供依据。
1 多孔介质中流体渗流的LB E2 回采工作面前方煤体中瓦斯流场的LB M 模型及算例模拟利用L B M 方法对回采工作面瓦斯渗流进行模拟，回采工作面前方煤体的竖直剖面图如图 1 所示。
由于顶底板的透气性远小于煤体的透气性，因此在研究中可以将顶底板视为不透气的挡板。
元胞自动机(CellularAutomata),简称CA,也有人译为细胞
元胞自动机(Cellular Automata)，简称CA，也有人译为细胞自动机、点格自动机、分子自动机或单元自动机)。
是一时间和空间都离散的动力系统。
散布在规则格网 (Lattice Grid)中的每一元胞(Cell)取有限的离散状态，遵循同样的作用规则，依据确定的局部规则作同步更新。
大量元胞通过简单的相互作用而构成动态系统的演化。
不同于一般的动力学模型，元胞自动机不是由严格定义的物理方程或函数确定，而是用一系列模型构造的规则构成。
凡是满足这些规则的模型都可以算作是元胞自动机模型。
因此，元胞自动机是一类模型的总称，或者说是一个方法框架。
其特点是时间、空间、状态都离散，每个变量只取有限多个状态，且其状态改变的规则在时间和空间上都是局部的。
元胞自动机的构建没有固定的数学公式，构成方式繁杂，变种很多，行为复杂。
故其分类难度也较大，自元胞自动机产生以来，对于元胞自动机分类的研究就是元胞自动机的一个重要的研究课题和核心理论，在基于不同的出发点，元胞自动机可有多种分类，其中，最具影响力的当属S. Wolfram在80年代初做的基于动力学行为的元胞自动机分类，而基于维数的元胞自动机分类也是最简单和最常用的划分。
除此之外，在1990年， Howard A.Gutowitz提出了基于元胞自动机行为的马尔科夫概率量测的层次化、参量化的分类体系(Gutowitz, H.A. ,1990)。
下面就上述的前两种分类作进一步的介绍。
同时就几种特殊类型的元胞自动机进行介绍和探讨S. Wolfrarm在详细分忻研究了一维元胞自动机的演化行为，并在大量的计算机实验的基础上，将所有元胞自动机的动力学行为归纳为四大类 (Wolfram. S.，1986):(1)平稳型:自任何初始状态开始,经过一定时间运行后，元胞空间趋于一个空间平稳的构形，这里空间平稳即指每一个元胞处于固定状态。
不随时间变化而变化。
(2)周期型:经过一定时间运行后，元胞空间趋于一系列简单的固定结构(Stable Paterns)或周期结构(Perlodical Patterns)。
人流模拟总结
现有的人群仿真研宄将其研宄重点大体可以分成三个方面,即路径规划问题、行为决策问题、动画表现问题;这三点也被称为是人群仿真研宄关注的三要素。
一．人员疏散模拟：宏观模型：（讲人群整体作为建模对象，不考虑个体间相互作用）流体力学模型。
微观模型：（对行人的详细行为进行描述）（1）连续型：(a)磁场力模型(疏散方式单一，参数由实验得出，没有计算公式)(b)社会力模型（变量所代表的物理意义是可计算的。
结果显示了现实中人群运动自组织现象，可以模拟人员的冲撞，挤压，恐慌等特征）。
（2）离散型：（时间，空间离散化）(a)元胞自动机（根据行为规则建模，行为模型根据出现的结果改变状态）(b)格子气模型（元胞自动机的特殊形式）(c)排队网络模型（基于离散事件的蒙特卡洛仿真）(d)智能体模型（agent-based对行人决策和动作方式建模, 更符合人类的行为表现;能模拟人类的思维过程，不能够真实模拟行人所处的道路网）。
模拟步骤：疏散模型相关因素：（1）楼宇空间结构：走廊过道的长度与宽度、出入口数量及宽度等（2）人流因素：人群速度，人群密度，人体所占空间，人群构成与状态等现有人群疏散研究方向：（1）单一、多个出口的房间、大厅；单向、双向，多向通道的人流交汇场景。
（2）人群流动特性分析（拥挤人群的基本特性, 如密度、速度与激波的关系）（3）房间，通道出口对人群疏散的影响（单出口，多出口，宽度，角度等因素）二．路径规划模拟以楼宇内部各节点网络的构建作为切入点，分析楼宇节点网络模型的建立、空间数据的拓扑分析和最优路径规划。
建立楼宇网络模型，主要包括点和线两种类型，即拐点、中心点、障碍点、弧段等。
进行拓扑分析和网络分析，进行基于距离最短、时间最少和有障碍物情况下等不同情况的最短路径规划。
优势是将最短路径分析的内容从距离延伸到了时间、花费、障碍物等复杂要素，从而提高该算法的适用性和可靠性。
不足是时间复杂度较高，需要人为设定要素的不同权重。
格子Boltzmann方法原理及其应用
格子Boltzmann方法原理及其应用摘要在上世纪八十年代后期提出的格子Boltzamnn方法克服了格子气方法的缺点，其本身也在不断的发展之中.格子Boltzamnn方法在流体运动计算方面展现了非凡的风采，成功地模拟了包括均相不可压缩湍流和多孔介质中的多相流动在内的流体动力学问题.但和成熟的流体动力学计算方法相比，特别在工程实际应用上，该方法还有许多值得研究的地方.本文主要介绍工程实际应用时，具体模型的选择问题.首先从理论上对应用最为广泛的几种基本模型进行了详尽的分析和比较.选择了Poiseuille流动，然后从计算精度、数值稳定性和收敛速度这几个方面进行了细致的比较.从理论和实验两个角度验证了D2G9模型的优越性，为工程实际应用上模型的具体选择提供了一定的参考依据.通过研究二阶精确的格子Boltzamnn模型，提出了非牛顿流体.非牛顿流动性是使用幂法则模型实现的.它可以估算出模型的精确程度，同时不会限制这个模型.二阶精度由剪切变稀和剪切增稠液体的幂法则模型参数范围给出.这些结果与Gabbanelli等人的结果相比，精确度更高，并且得到了更快的计算效率.结果表明了格子Boltzamnn方法适用于非牛顿流体模拟.对于实际流动模拟，本文应用二维9速度模型模拟了四种情况的方柱绕流问题.在第一种情况中，单个方柱位于流场中央，给出了流线图，等涡线图，模拟了卡门涡街现象，并计算了升、阻力系数，Strouhal数等参数；在第二种情况中，计算细长矩板截面柱绕流问题，得到了Strouhal数随着矩形长宽不同的比值下的变化情况；在第三种情况中，两个方柱并列位于流场中央，考察了方柱间距对于流场的影响；在第四种情况中，计算了水平来流为剪切流的方柱绕流问题，比较了速度梯度取不同值下流场的变化情况.所有有关力的求解均采用动量转换法.所得结果，包括流线、等涡线、升/阻力系数曲线等均与已有文献的实验或数值结果基本一致，显示LBM方法及其力的求解方法——动量转换法是有效的，能够精确的模拟各流场.其次，我们还引入一种两相耦合机制对D2G9模型进行了修正，从而使之可以正确处理气固两相流中输运相和颗粒相之间的相互作用.随后，我们模拟了后台阶流动，并和传统CFD方法的模拟结果以及修正其他模型的模拟结果进行了验证，得到了令人满意的结论.从一定程度上验证了两相耦合机制的可行性.通过软件模拟获得了水包油、过渡流型和油包水三种流型的典型模拟图.经分析发现：由软件模拟的流型特点和由探针获得的流型特点具有较好的一致性.在本文最后，我们介绍了以经典算例一方腔流为例，对格子Boltzamnn方法的核心代码进行了优化的方法，主要讲述对时间和空间上的优化，优化的程序使计算效率提高数倍.在并行的框架下，核心演化的代码换为优化后的程序，计算效率有大幅度的提高.关键词：格子方法；格子Boltzamnn 方法；格子气自动机；格子Boltzamnn模型.AbstractIn the latter of 80’s，the Lattice Boltzamnn Method（LBM）was introduced mainlyto cope with major drawbacks of its ancestor，the Lattice Gas Automata（LGA）．Eversince，it has undergone a number of refinements and extensions which have taken it tothe point where it can successfully compute a number of non trivial flows，raging fromhomogeneous incompressible turbulence to multiphase flows in porous geometries．Yet，when compared with conventional computational fluids dynamics methods，such as finiteelement，finite difference，it is apparent that there is still a way to go before LBM canachieve full engineering status．In this paper，we mostly focus on the choice of the basic LB models in theengineering application fields．Firstly，we expatiate the basic LB models in theory．Then，we simulate the Poiseuille flow with those basic LB models．And wecompare the simulation results from the computation precision、the numerical stabilityand the convergence rate．Finally，we draw a conclusion that the D2G9 model is the bestchoice in the engineering application fields．Simulation of Flow past square cylinder with LB Method.For the simulation of actual flow，we use D2Q9 investigate fourcases of flow past square cylinders in this paper.For case 1，one singlesquare cylinder is located at the center of the channel，we describe thestreamline contour，vortices contours，simulate the Karman vortex，then compute the lift coefficient，drag coefficient，Strouhal numbersetc.For the case 2，simulate the flow past a cylinder of rectangularcross-section;compute the change of Strouhal numbers varying withthe side ratio.For case 3:two square cylinders arranged side by side inthe center of the channel，the flow features at different spacing ratiosare studied.For case 4:we compute the linear shear flow over a squarecylinder，compare the evolution of flow with different velocitygradient.The results of thesimulation including the streamlines，vorticity contours，lift and drag coefficients etc.are agreed with thoseof available literatures，and show that LB method and itsmomentum-exchange method can achieve accurate results and obtainthe reasonable flow in detail.we employ a two-way coupling mechanisms to modify theD2G9 model．With the modified D2G9 model，we can handle with the interactionsbetween carrier phase and dispersed phase in the model．Then，we simulate abackward-facing step model，and the results are compared qualitatively with the result ofthe traditional CFD method and the other modified LB models．Though the comparison，we can see that the two-way coupling mechanisms can handle with the gas-solid twophases flows successfully．Three kinds of flow pattern，which are oil-in-water flow，transitional flow andwater-in-oil flow，have been got by simulation.According to the result of simulation，theoil-water two-phase flow pattern transition boundary model has been got by.By the analysisof simulation，the characteristic of three kinds of flow pattern of vertical oil which has beengot by analysis of the signals is consistent with results by simulation.We take the classical problem-cavity flow as an example and optimize the kerne codes of the LBM. The optimization include two aspects :time and space .The efficiency of the optimized code increased much more .In the parallel frame，the efficiency also increased if the kernel code is taken the optimized code.Key word：1atrice method；1atrice bohzmann method；lattice gas automata;LBM目录第1章概述 11.1研究格子 Boltzamnn方法的意义 11.2 格子 Boltzamnn方法的发展历程 31.2.1孕育阶段 31.2.2 萌芽到成长阶段 31.3 格子 Boltzamnn方法应用概况及优缺点 51.3.1格子Boltzamnn方法应用概况 51.3.2格子Boltzamnn的优缺点 61.4本论文的研究目的 81.5 相关研究的综述与专注情况 8第2章格子Boltzamnn方法介绍 102.1 Boltzamnn方程的产生 102.2细胞自动机（CA） 112.3格子气自动机（LGA） 122.4格子Boltzamnn方法（LBM） 132.5 格子Boltzamnn的基本结构 162.6本章小结 17第3章格子Boltzamnn方法的基本模型比较 183.1 格子 Boltzamnn 方法基本模型概述 183.2 进行常压力梯度驱动的Poiseuille流动模拟比较几种基本模型 23 3.3本章小结 27第4章格子Boltzamnn方法的算法设计 284.1格子Boltzamnn方法的算法实现 284.2格子Boltzamnn方法的高效算法设计 304.2.1优化算法 304.2.2优化实验 324.3 本章小结 34第5章格子Boltzamnn方法的实际应用 355.1二阶精确格子Boltzamnn非牛顿流体的流动模拟 35 5.1.1理论背景 355.1.2方法和计算结果分析 385.1.3 本节小结 405.2 格子Boltzamnn方法的方柱绕流模拟 405.2.1 单个方柱位于流场中央的绕流问题 405.2.2 细长矩形截面住绕流问题 425.2.3 两个并列方柱的绕流问题 445.2.4来流为剪切流的绕流问题 495.3格子Boltzamnn方法模拟气固两相流 515.3.1对气固两相流的模拟模拟对象简介 515.3.2 计算结果分析 545.3.3本节小结 565.4 格子Boltzamnn方法模拟油水两相流软件设计 565.4.1 LBM油水两相流的关键因素选取 575.4.2 软件的设计 605.4.3 本节小结 635.5 简述格子Boltzamnn方法在其他领域中的应用 645.5.1 颗粒悬浮问题的模拟 645.5.2 热导和对流—扩散问题的模拟 645.5.3 偏微分方程的模拟 655.5.4 多相流和多元流的模拟 65结论及展望 67参考文献 68第1章概述1.1研究格子Boltzamnn方法的意义自从二十世纪四十年代出现了第一台电子计算机以来，人们开始进入了电子信息时代.随着高存储、高速度计算机的出现，人们所能解决的问题也越来越广泛，同时所面临的问题也越来越复杂.在对流动现象的研究中，以往人们大部分依靠的是解析方法，但所解决的问题非常有限.而现实生活中所面临的流动问题往往十分复杂，如航空航天器的亚跨超音速飞行、舰船的航行等等，依靠解析的方法来解决这些复杂的流动现象是不可能的.到现今为止，人们对流体运动的研究主要靠实验方法和数值计算方法.实验方法具有直观、结果基本可靠的特点.但也存在较大的缺点：耗费大、周期长，并且结果受实验条件的影响也较大，尤其是如今的航空航天飞行，速度高、飞行条件复杂，用风洞来模拟困难是相当大的.而流体的运动可以由一组偏微分方程描述.在大多数情况下，这些方程（如N-S方程）都是高度非线性的，采用解析的求解方法是不实际也是不可行的.随着大型计算机的出现，使人们可以借助于计算机用数值计算方法来解决复杂的流动问题.因此，在二十世纪六十年代，用数值方法分析求解流动问题的学科——计算流体力学（CFD）逐渐发展起来.伴随着电子计算机的飞速发展以及各种新颖算法的不断出现，CFD已经形成了一门独立的学科，并且在航空航天、船舶、大型能源装置（如核电站）、新型交通工具、海洋工程、环境保护等众多工程技术部门和领域都得到了广泛的应用.随着计算技术的发展、巨型计算机的出现、计算方法的不断改进，计算流体力学在解决流动的理论和工程实际问题中愈加显示出它的巨大作用.目前，计算流体力学已经成为现代计算科学的最有力的推动力之一.在计算流体力学中，传统的数值模拟方法可以分为两大类：（1）从宏观角度出发，基于连续介质假设，采用数值计算方法，求解全位势方程或Euler方程或N-S方程；（2）从微观角度出发，采用分子动力学的方法，对流动进行数值模拟.其中，格子Boltzamnn方法就是典型的一种.格子Boltzamnn方法（Lattice Boltzamnn Method，LBM）1.1.2格子Boltzamnn法（lattice Boltzamnn method）起源于格子气自动机（Lattice Gas Automata，LGA）.LGA方法是元胞自动机（Cellular Automata，CA）在流体力学中的具体应用，是空间、时间和速度空间都离散的一个虚拟微观模型，与以连续微分方程为基础的宏观计算流体力学方法有着本质的不同.LGA的微观特性使得它的边界条件非常容易实现，并且计算也很简单.因此，LGA方法非常适于处理边界复杂的问题.更为重要的是，LGA的计算具有局部性和并行性，非常容易在并行机上实现.LGA的出现不但为并行计算提供了许多新思想，而且对并行计算机制造技术产生了重要的影响.但是，LGA方法也有许多不足之处.例如，由于含有随机因素，LGA的计算结果往往包含很大的统计噪声，LGA的宏观方程也不是标准的流体运动宏观方程.格子Boltzamnn方法是为克服LGA方法的一些内在不足而发展起来的一种新方法.LBM不但克服了LGA的缺点，继承了LGA的主要优点，而且还有许多新的优点，如计算量小、计算效率高、编程简单等.LBM的产生与发展，不仅在计算流体力学领域中产生了深远的影响，它所使用的处理方法和观点对其他许多学科也是富有启发性的.格子Boltzamnn法是一种应用非连续介质思想研究宏观物理现象，并可平行运行，求解流体力学问题的新方法.它是由格子气自动机（lattice gas automata，简称LGA）方法发展而来的.该法把流体及其存在的时间、空间完全离散，把流体看成由许多只有质量没有体积的微小粒子组成，所有这些粒子同步地随着离散的时间步长，根据给定碰撞规则在网格点上相互碰撞，并沿网格线在节点之间运动.碰撞规则遵循质量、动量和能量守恒定律.流体运动的宏观特征是由微观流体格子相互碰撞并在整体上表现出来的统计规律.该法是直接从微观模型出发，经过Boole化处理后进行计算，可认为是N-S差分法逼近的一种无限稳定的格式.被广泛应用于复杂几何边界流体流动、多孔介质流、多相流及反应流等.格子气自动机的基本思想是，把计算区域分成许多均匀的正三角形（或正方形）的网格，而那些只有质量无体积的粒子只能在网格点上存在，并沿着网格线在网格间运动.当某一个粒子从某一网格点到邻近的网格点时，有可能和从其他网格点到达该点的粒子相碰撞.根据Pauli不相容原理，在同一时刻同一点上，沿着每一网格线运动方向最多只有一个粒子，流场中的粒子速度不是0（静止）就是1（设格子边长及时间间隔都为1）.以三角形网格为例，每一个网格上在某一时刻，其周围的6个网格上粒子沿着网格线聚集到该点，加上该点可能还有一个静止粒子，这样，可能有7个粒子在该点发生碰撞，然后根据碰撞规则再散射出去，演化为新的运动粒子流向各节点的邻居，形成格子气自动机.1986年MeNamaxa和Zaneltti，提出把格子气自动机中的整数运算变成实数运算，建立了格子Boltzamnn 模型，克服了格子气自动机的数值噪声的缺点.后来陈十一和钱跃宏采用了单一时间松弛方法，满足了各项同性，GalIean不变性，并得到了独立于速度的压力项．使格子Boltzamnn模型保留了格子气自动机的优点，克服了其不足，并在理论分析和数值模拟方面都具有很大灵活性，而且程序编制简单，计算效率较高.从格子Boltzamnn方法诞生至今天已有20年，20年间，其在理论和应用研究等方面都取得了迅速发展，并逐渐成为在相关领域研究的国际热点之一，受到国内外众多学者关注.与之传统模拟方法不同，格子Boltzamnn方法基于分子动理论，具有清晰的物理背景.该方法在宏观上是离散方法，微观上是连续方法，因而被称为介观模拟方法.在许多传统模拟方法难以胜任的领域，入微尺度流动与换热、多孔介质、生物流动、磁流体、晶体生长等，格子Boltzamnn方法都可以进行有效的模拟，因此它被用于多种复杂现象的机理研究，推动了相关学科的发展.可以说，格子Boltzamnn方法不仅仅是一种数值模拟方法，而且是一项重要的科学研究手段.此外，格子Boltzamnn方法还具有天生的并行特性，以及边界条件处理简单、程序易于实施等优点.可以预计，随着计算机技术的进一步发展，以及计算方法的逐渐丰富，格子Boltzamnn方法将会取得更多成果，并为科技发展发挥更重要的作用.1.2 格子Boltzamnn方法的发展历程格子Boltzamnn方法自诞生至今年已取得了长足发展，被誉为现代流体力学的一场变革.1.2.1孕育阶段：对格子Boltzamnn方法发展使得了解，得先从格子自动机说起.格子气自动机使更广泛的元胞自动机在流体学中的应用.元胞自动机是一个时间和空间离散的数学模型.20世纪60年代，Broadwell等人首先提出了离散速度模型，用以研究流体中的激波结构.20世纪70年代，为了研究流体的运输性质，法国的Hardy、Pomeau和Pazzis提出了第一个完全离散模型，该模型命名HPP模型.这是历史上的第一个格子气自动机模型.1986年，法国的Frisch、Pomeau和美国的Hasslacher提出具有足够对称的二维正六变形格子气自动机模型，，命名为FHP模型.由于这些方法在还处在一些缺点：（1）有格子气自动机演化方程推导出来的动量方程不满足Gaililei不变形；（2）流体状态方程不仅仅依赖于密度和温度，还与宏观流速有关；（3）破装蒜子具有指数复杂性，对计算量和存储量也有较大要求.因而，我们将这一段格子气自动机的发展过程称作格子Boltzamnn方法的孕育期.1.2.2 萌芽到成长阶段：自1988年底一篇关于格子Boltzamnn方法的论文出现至今，格子Boltzamnn方法从萌芽逐渐成长壮大，并成为目前一大国际研究热点，受到越来越多学者的关注.1988年，McNamra和Zanetti提出把格子气自动机中的Bool运算变成时数运算，格子点上的粒子数不是用整数0或1来表征，而是用实数f来表示系综平均后的局部粒子分布函数，用Boltzamnn方程代替格子气自动机的演化方程，并将该模型用于流体的数值计算.这是最早的格子Boltzamnn模型，从此开启了格子Boltzamnn方法的历史大门.1989年，Higuera和Jimenez提出了一种简化模型：通过引入平衡分布函数，将碰撞算子线性化.该模型不需要碰撞模型，并忽略各自粒子间的碰撞细节，相比于多粒子碰撞模型，容易构造.同年，Higuera等进一步提出了强化碰撞算子方法，以增加模型的数值稳定性.这两模型统成为矩阵模型.经历了上述两类模型，格子Boltzamnn方法消除了统计噪声，克服了碰撞算子指数复杂性，但是由于依然使用Fermi-Dirac平衡态分布函数，格子气自动机的其他缺点仍然存在.1991年，Chen等提出了单松弛时间法，用同一个时间松弛系数来控制不同例子靠近各自平衡态的快慢，进一步简化了碰撞算子；Qian等人在1992年也提出了类似的方法，称之为格子BGK（LBGK）模型.LBGK模型与矩阵模型类似，但与前面两种模型不同的是，当粒子种类数增加时，碰撞算子本身发生生变化，不会变得复杂.至此，格子Boltzamnn方法完全克服了格子气自动机的一系列缺点，并逐渐成熟，成为国际研究的热点.早期的格子Boltzamnn模型只能用于等温不可压缩流动的模拟.但因为存在可压缩效应，会引起一定的误差.为了消除或强敌有可压缩效应引起的误差，许多学者致力于新的格子Boltzamnn模型的研究，并提出了多种等温不可压模型.而后，一些不可压缩热模型成功实现了对有效范围温度变化的热力学和传热学问题的模型.其中，最成功的要数双分布函数模型.他是在密度分布函数的基础上引入了温度分度函数、或内能分布函数、或总能分布函数，并用密度分布函数演化得到速度场，这类模型具有与等温不可压模型相同的数值稳定性，而且可以从根本上解决压缩功和耗热问题.边界处理方面，经历了20年的发展，格子Boltzamnn方法已逐渐发展出适合不同边界条件、不同模型的边界处理格式.网格划分方面，最初的格子Boltzamnn方法是基于正六边形或正四边形的均匀对称网格.由于均匀网格在计算效率、计算精度等方面的不足，从而促进了非均匀网格、多快以及多重网格、无网格等多技术出现.总的来说，这些网格技术延展了格子Boltzamnn方法的应用范围，使得格子Boltzamnn方法主机去年从理论的神殿走向更可能多的实际应用领域.1.3 格子boltzamnn方法应用概况及优缺点1.3.1格子boltzamnn方法应用概况与传统的宏观数值方法相比，具有介观特性的格子Boltzamnn方法其主要优点是物理图像清晰、便捷容易处理以及并行性能好等.因而自诞生之日起，格子Boltzamnn方法就得到了国内外学术界的广泛关注，并寄希望该方法能再注入为尺度流体、多相流、多孔介质内流动与换热、化学反应流等传统法就延受限的领域取得开拓性进展.事实上，在20年的发展过程中，格子Boltzamnn方法的确也已成一个十分活跃极具发展前景的模拟手段.并迅速在微/纳米尺度流、多孔介质流、多相多质流、非牛顿流体、粒子悬隔i浮流、湍流、化学反应流、燃烧问题、磁流体、晶体生长等许多领域得到应用.下面分别以多孔介质流、多相流和非牛顿流体三个方面为例，做较详细说明.由于格子Boltzamnn方法边界条件易于实施，在模拟具有复杂几何构型的问题具有较大的优势，因而这个方向的发展非常迅速.目前，采用格子Boltzamnn方法对多孔介质流进行模拟主要在空隙尺度和代表单元尺度上进行.在孔隙尺度上，可以直接使用格子Boltzamnn方法描述孔隙内的流体流动，多孔介质则当做固体壁面，流体与介质相互作用使用边界处理格式来描述.在多相流方面，由于真实的流动问题常常是多相的，因而对其开展研究具有重要的现实意义.由于格子Boltzamnn方法的介质特性，它可以方便地描述数流动中不同相之间的相互作用，因而在多相流领域具有较好的应用前景.按照设计方法的不用，现有模拟多相流的格子Boltzamnn模型可分为四大类：着色模型、伪势模型、自由模型和其他模型.格子Boltzamnn方法在非牛顿流体领域的应用刚刚起步，主要研究对象是非牛顿幂律流体.Aharonov等最早提出使用矩阵碰撞该算子来计算幂律流问题，即在每一个时步内，调整碰撞算自来该表局部的动力学黏性系数.Boek用该模型模拟了幂律流体在简化多孔介质中模型的流动，模拟结果与达西定律符合良好.最近，Gabbanelli又对上述模型进行了改进，引入分段幂律方程描述剪切率和表现黏度的关系.以上可看出，到目前为止，格子Boltzamnn方法的研究者主要局限在科学界.尽管如此，随着格子Boltzamnn 方法理论体系逐渐完善，以及计算机技术的进一步发展，格子Boltzamnn方法也会走向更加广泛的工业实际应用中.1.3.2格子Boltzamnn的优缺点流体力学的理论描述通常建立在纳维--斯托克斯方程的基础上，作为流体力学的基石，它已处在了一个多世纪.在通常尺度下，|人们对此方程的物理可靠性即准确性并不抱异议.理论上人们一般通过求纳维--斯托克斯方程及其各种简化形式的途径来处理复杂的流体力学问题，现行的计算流体力学研究也主要是围绕着纳维--斯托克斯方程的计算方法展开的.然而，基于其本质上的非线性以及边界条件处理的困难，除少数简单问题外，解析和数值求解纳维--斯托克斯方程都是极具挑战性的任务.除了求解的困难外，作为一种对流体物理的描述，与描述经典力学运动的牛顿运动方程，或与描述量子力学运动的薛定谔方程等原理方程不同，纳维--斯托克斯方程是从更根本的原理性方程出发，在合理地假定某些物理机制可以忽略后，经过统计平均得到的.本质上纳维--斯托克斯方程当然不可能描述那些被忽略了的物理机制带来的宏观现象，比如流体系统中的相变、非牛顿的本构关系以及在分子运动自由程尺度上的物理现象，在这些领域，纳维--斯托克斯方程明显的显示出了他的局限性.从20世纪80年代末开始，一种对于流体力学的全新的理论表相及有效的计算方法初步形成，这就是现在人们通常所谓的格子Boltzamnn方法.关于格子Boltzamnn方法的早期发展，上文已有较全面的综述，在此仅作简单介绍.从历史角度来讲，格子Boltzamnn方法最初是从所谓的格子气模型演化而来的，而后者是一种抽象简化的分子运动数学模型.格子Boltzamnn方法最初的引入有两个主要原因：一是为了降低模型导致的数值噪音；而是能够克服格子气模型里处在的非物理缺陷.可以证明，格子Boltzamnn系统的宏观表象基本满足纳维--斯托克斯方程.从而，人们可以模拟格子Boltzamnn系统地方法来间接地解纳维--斯托克斯方程.标准格子Boltzamnn方程一般用一下的数学表达式描述：式中——粒子分布函数；——碰撞项.用格子玻尔兹曼模型进行流体的数值模拟有一些明显的优越性．如，它的对流（advection）过程是通过常数值速度实现的．这相应的计算是一项极其简单的操作步骤．当适当的格子网格选定后，该过程通常可以用完全平移的方式实现．用计算数学里的常规有限插值语言来讲，它对应于上风插值．但所不同的是其对应的柯郎数（Courant Number）等于1．相比之下，纳维——斯托克思方程的对流项是一个随时空变化的非线性函数．众所周知，对于它的计算不是一项简单的事，并且，数值稳定性的要求迫使人们在实际问题的计算中只能使用比1小得多的柯朗数．在给定空间分辨度的情况下，小柯朗数意味着小时间步长，从而大大延长了计算时间：同时，小柯朗数也增大了数值扩散误差，迫使人们采用更高精度格式或隐式格式．其后果是，或者算法变得极为复杂，并行效率大大降低；或者计算只限制在处理定常流的情况下．事实上，定常流是对流动情况的极大限制．许多重要的流体力学问题，如分离流，即使我们只关心它的时间平均的结果，也是不能用定常流假设来近似的．在此我们也要提一下格子玻尔兹曼方程的另一个本质特性：所有非线性效应在格子玻尔兹曼方法里都包含在碰撞项中，并且是以纯粹局部信息的方式体现的．这进一步发挥了并行计算的长处．所有这些理由意味着格子玻尔兹曼方法是对非定常流动实行大规模并行模拟计算的一种比较优越的方法．相比之下，以流体力学方程（纳维一斯托克思方程或Burnett类型方程）宏观描述为基础的传统计算方法对许多这类问题存存基本困难．除边界条件之外，利用各种封闭性假设推导出的超越纳维一斯托克思的宏观方程直至现今仍存在对其数学规范性的疑问和争议，多相流的计算也存存同样问题．众所周知，流体系统中存在多相的物理机制是分子问的长程作用力，这种机制早已超出了流体力学方程所能描述的物理现象范围．以流体力学方程为基础的多相流计算方法必须依赖额外的模型来模拟流体力学方程本身所不包含的物理现象．除了实际数值结果显示的问题之外，这种方法本质上隐含着严重的基本物理缺陷，这种缺陷集中表现在对相交界面的准确描述上面，即在十分尖锐的相界面附近，纳维一斯托克思方程之类近平衡态的近似表象是有相当疑问的．这也反映在相界面和兀滑动（no—slip）固体边界条件的互斥性上面，为了修补这一缺憾，人们不得不引入各种滑动经验模型.反之，以细观（mesoscopic）为表象基础的格子玻尔兹曼方法可容忍更大的非平衡态程度及更广义的严格边界条件．另外，压力的状态方程在细观表象中是由粒子的相互作用自然得出的，而不用直接输入和处理．在相变情况下，物体的宏观特性将产生不连续性，而对应的微观和细观力学机制并无改变．格子玻尔兹曼方法在模拟多相流上有着广泛的使用.然而，这种为大多数人所熟悉的格子玻尔兹曼方法的理论框架存在本质上的缺陷．由于它运用逆向切普曼一安斯柯格展开的途径来适定平衡态分布函数中的关键参数，以达到复建宏观物理体系的目的，这就使其。
【国家自然科学基金】_格子气模型_基金支持热词逐年推荐_【万方软件创新助手】_20140802
推荐指数 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2013年 序号 1 2 3 4
科研热词 相分离过程 生长律 格子气模型 monte carlo模拟推荐指数 1 1 1 1
2014年 序号 1 2 3 4
2014年 科研热词 预辐照 氦泡生长 抑制作用 monte carlo模拟 推荐指数 1 1 1 1
2008年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
科研热词 人员疏散 群体疏散 疏散模型 生物实验视角 格子气模型 格子气 微观模型 建筑空间 宏观模型 多速度人群 多格子模型 博弈论 区间估计 创新性成果 代表性成果 人流疏散
推荐指数 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2009年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13
科研热词 计算效率 行人交通仿真 虚拟实验 自组织 群体疏散 离散 社会力模型 格子气 微观仿真 双向行人流 元胞自动机模型 个体差异 agent模型
推荐指数 1 1 1 1 1 1 1 1 1 1 1 1 1
2011年 科研热词 非阿尔奇现象 物理模型 格子气自动机 数字岩心 效率 并行 含水饱和度 加密系统 推荐指数 1 1 1 1 1 1 1 1
2012年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
科研热词 气固两相流 视距 自激振荡 背景场 研究进展 瓶颈 热声装置 热声效应 潜意识 格子波尔兹曼 格子气模型 格子气 拦截 惯性碰撞 布袋除尘 布朗扩散 大涡模拟 后台阶流动 双向行人流 元胞自动机模型 openfoam
行人交通模型与微观仿真
方法融合了解析方法和实验方法的优势 , 在采集行人 运动行为特性数据的基础上 , 建立行人运动的模型 , 通过离散事件或连续事件仿真来直观地刻画行人交 通 。相比单纯的解析方法和实验方法 , 仿真更适宜描 述大规模、复杂的行人交通 , 并且能够进行重复性的、 现实中难以开展的实验 (如发生紧急情况 ) , 非常适合 与行人交通有关的设计以及运营方案的评估和检验。 本文回顾行人交通仿真的发展概况 , 探讨和比较 宏观模型 、中观模型 、微观模型的原理和算法 , 重点 选择了 Legion、 STEPS、 SimW alk 和 AnyLogic 四种适 用于正常模式和疏散模式的微观仿真软件 , 研究了如 何利用上述软件进行快速建模 , 并比较分析了其技术 性能 。
212 中观模型
中观模型以格子气模型为代表 , 融合了宏观模型 和微观模型 , 将平面划分为小格子或三角形 , 行人位 于交点处 , 其运动方向为前 、左 、右 3 个方向 , 依照 这 3 个方向的确定值来决定下一步运动方向 。格子气 模型从个体行人的角度建模 , 大多不考虑行人间的相 互作用 , 较微观模型略显粗糙 。
→作用力的影响 : ( 1 ) 驱动 力 , 主观意识对个体行为的影响可化为个体所受自己 施加的 “ 社会力 ” , 体现了行人以渴望的速度移动到 目的地的动机 ; ( 2 ) 人与人之间的作用力 , 指试图与 其他行人保持一定距离所施加的 “ 力” ; ( 3 ) 人与边界 之间的作用力 , 边界和障碍对人的影响类似于人与人 之间的作用 。模型用如下方程组表示 。
在微观仿真中每位行人被当作独立个体进行描述和计算不仅能从宏观和整体上把握行人交通流而且能从局部和细节上刻画行人交通的复杂行为已被成功用于交通枢纽和地铁车站的乘客集散分析体育比赛或集会活动的人群疏散研2行人交通模型行人交通仿真的困难主要体现在行人运动随意性大步行行为复杂影响因素复杂等诸多方面
REV尺度多孔介质格子Boltzmann方法的数学模型及应用的研究进展
CHEMICAL INDUSTRY AND ENGINEERING PROGRESS 2016年第35卷第6期·1698·化 工 进 展REV 尺度多孔介质格子Boltzmann 方法的数学模型及应用的研究进展张潇丹1，2，雍玉梅2，李文军3，赵元生4，李媛媛2，杨巧文1，杨超2（1中国矿业大学（北京）化学与环境工程学院，北京 100083；2中国科学院过程工程研究所绿色过程与工程重点实验室，北京 100190；3华北科技学院环境工程学院，河北 廊坊 065201；4中国石油化工研究院渣油加氢实验室，北京 102200）摘要：综述了多孔介质表征体元尺度（REV ）格子Boltzmann 模型的研究进展，根据对多孔介质处理方式主要分为部分反弹模型和阻力模型两类，分析归纳了各类模型的优缺点。
由于阻力模型中渗流的广义格子Boltzmann 方程（GLBE ）的作用力是基于GUO 等的作用力模型，可以准确得到宏观方程，不存在离散误差，且模型的平衡分布函数和作用力项中都包含反应介质特性的孔隙率，因而应用最为广泛。
本文还重点介绍了REV 尺度多孔介质LBE 模型在流动、传热、传质、化学反应及相变等过程中的具体应用，认为REV 尺度多孔介质内的三传一反数学模型中需要加入孔隙尺度因素，在更大工程尺度上应该考虑过程参数的各向异性，展望了REV 尺度多孔介质LBE 模型的发展和应用前景。
关键词：多孔介质；表征体元尺度；格子Boltzmann 方法；流动；传热；传质中图分类号：TQ021.9 文献标志码：A 文章编号：1000–6613（2016）06–1698–15 DOI ：10.16085/j.issn.1000-6613.2016.06.010Models and application of lattice Boltzmann method at REV-scalein porous mediaZHANG Xiaodan 1，2，YONG Yumei 2，LI Wenjun 3，ZHAO Yuansheng 4，LI Yuanyuan 2，YANG Qiaowen 1，YANG Chao 2（1School of Chemical & Environmental Engineering ，China University of Mining & Technology (Beijing)，Beijing100083，China ；2Key Laboratory of Green Process and Engineering ，Institute of Process Engineering ，Chinese Academy of Sciences ，Beijing 100190，China ；3School of Environmental Engineering ，North China Institute of Science and Technology ，Langfang 065201，Hebei ，China ；4Laboratory of Residue Hydrotreating ，Research Institute of PetroleumProcessing ，PetroChina ，Beijing 102200，China ）Abstract ：This paper discusses the lattice Boltzmann model at representative elementary volume (REV)scale for porous media. According to different treatments of porous media ，the lattice Boltzmann model at REV-scale for porous media can be classified into two categories ，the partially bouncing-back model and the resistance model. The advantages and disadvantages of various models are analyzed. The Generalized lattice Boltzmann equation (GLBE model) in the resistance model is most widely used. Firstly ，the force item of the GLBM model is based on the method proposed by Guo et al ，which can be第一作者：张潇丹（1989—），女，硕士研究生，主要从事化学工程数值模拟。
计算流体力学中的LBGK方法的理论研究
武汉理工大学硕士学位论文计算流体力学中的LBGK方法的理论研究姓名：胡娟申请学位级别：硕士专业：流体力学指导教师：熊鳌魁2003.4.1摘要ＬａｔｔｉｃｅＢＧＫ，简称ＬＢＧＫ，是一种非常新颖的数值计算方法，它在格子气（ＬａｔｔｉｃｅＧａｓＡｕｔｏｍａｔｏｎ，简称ＬＧＡ）的基础上发展而来，并己在计算流体力学中得到广泛的应用，它可以用来模拟各种流动现象。
作为一种新的计算方法，对它的研究是非常有意义的。
本文总结了各种ＬＢＧＫ模型，特别是对ＬＢＧＫ模拟二相流和多相流模型进行了认真地研究，同时对ＬＢＧＫ进行了理论研究，用全新的方法探讨了ＬＢＧＫ和Ｎ．ｓ方程的本质联系，为推广ＬＢＧＫ方法进行了基础理论研究工作。
ＬＢＧＫ方法有许多优点，如边界条件容易处理，并行度高。
当然，它同时也有一些缺点，如数值不稳定。
因此，本论文对ＬＢＧＫ的稳定性进行了研究，提出了一种新的解决方法。
本论文的研究，旨在全面剖析ＬＢＧＫ方法，希望通过这些基础研究工作，促进ＬＢＧＫ方法在计算流体力学中的广泛应用。
本论文的主要工作是：对ＬＢＧＫ的各种模型进行认真地总结分析，对ＬＢＧＫ方法作了基础研究工作；通过对ＬＢＧＫ方法的理论研究，指出了ＬＢＧＫ方法和Ｎ．ｓ方程的本质内在联系：研究ＬＢＧＫ方法的稳定性，提出了新的提高稳定性的方法。
即对平衡分布函数取对数，以保证平衡分布函数恒为正数，来改善ＬＢＧＫ方法的稳定性。
关键词：ＬＢＧＫ方法，模型，稳定性ＡｂｓｔｒａｃｔＬＢＧＫｍｅｔｈｏｄｉＳａｖｅｒｙｈｏｖｅｌｎｕｍｅｒｉｃａｌｍｅｔｈｏｄ．Ｉｔｄｅｖｅｌｏｐｓｏｎｔｈｅｂａｓｉｓｏｆｔｈｅｌａｔｔｉｃｅｇａｓａｕｔｏｍａｔｏｎ．ＴｈｅｓｕｃｃｅｓｓｆｕｌｓｉｍｕｌａｔｉｏｎｏｆｆｌｕｉｄｆｌｏｗｂｙｔｈｅＬＢＧＫｍｅｔｈｏｄｈａｓｄｅｍｏｎｓｔｒａｔｅｄａｗｉｄｅａｐｐｌｉｃａｔｉｏｎｉｎｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｆｌｕｉｄｄｙｎａｍｉｃｓ．ＩｔｉｓｖｅｒｙｗｏｒｔｈｗｈｉｌｅｔｏｒｅｓｅａｒｃｈＬＢＧＫｍｅｔｈｏｄｓｗｈｉｃｈｉｓａｎｅｗｎｕｍｅｒｉｃａｌｍｅｔｈｏｄ．ＴｈｉｓｔｈｅｓｉｓｓｕｍｍａｒｉｚｅｓｕｐａＵｋｉｎｄｓｏｆｔｈｅＬＢＧＫｍｅｔｈｏｄ，ｅｓｐｅｃｉａｌｌｙｔｈｅｍｏｄｅｌｏｆｔｈｅｔｗｏ—ｐｈａｓｅａｎｄｍｕｌｔｉｃｏｍｐｏｎｅｎｔｆｌｕｉｄ．Ａｔｔｈｅｓａｍｅｔｉｍｅ，ｔｈｉｓｔｈｅｓｉｓｆｏｃｕｓｅｓｏｎｔｈｅｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｎｔｈｅＬＢＧＫｍｅｔｈｏｄａｎｄｄｉｓｃｕｓｓｅｓｔｈｅｉｎｎｅｒｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｈｅＬＢＧＫａｎｄｔｈｅＮ－Ｓｅｑｕａｔｉｏｎｓ．ＬＢＧＫｍｅｔｈｏｄｈａｓｍａｎｙａｄｖａｎｔａｇｅｓ，ｓｕｃｈａｓｅａｓｙｄｅａｌｉｎｇｗｉｔｈｔｈｅｃｏｍｐｌｅｘｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎｓａｎｄｈｉｇｈａｍｅｎａｂｉｌｉｔｙｔｏｐａｒａｌｌｅｌｃｏｍｐｕｔｉｎｇ．ＯｆＣＯＵｒＳｅ，ｉｔｈａｓｓｏｍｅｄｉｓａｄｖａｎｔａｇｅ，ｅｓｐｅｃｉａｌｌｙｔｈｅｉｎｓｔａｂｉｌｉｔｙ．ＳｏＩｐｕｔｆｏｒｔｈａｎｅｗｍｅｔｈｏｄｔｏｉｍｐｒｏｖｅｔｈｅｓｔａｂｉｌｉｔｙｏｆＬＢＧＫｍｅｔｈｏｄ．Ｔｈｉｓｐａｐｅｒ’ＳａｉｍｉｓｔｏａｎａｌｙｚｅＬＢＧＫｍｅｔｈｏｄｉｎｅｖｅｒｙａｓｐｅｃｔａｎｄｔｏｄｏｓｏｍｅｆｕｎｄａｍｅｎｔａｌｒｅｓｅａｒｃｈｏｎｔｈｉｓｍｅｔｈｏｄｆｏｒｔｈｅｉｍｐｒｏｖｅｍｅｎｔｏｆｉｔｓａｐｐｌｉｃａｔｉｏｎｉｎＣＦＤ．Ｔｈｅｍａｉｎｗｏｒｋｏｆｔｈｅｐａｐｅｒｉｓｆｏｌｌｏｗｅｄ：Ｆｉｒｓｔ，ｔｈｅｂａｓｉｃｔｈｅｏｒｙａｎｄｍａｎｙｍｏｄｅｌｓｏｆＬＢＧＫｍｅｔｈｏｄａｒｅｉｎｔｒｏｄｕｃｅｄａｎｄｄｉｓｃｕｓｓｅｄ．Ｓｅｃｏｎｄ，ｔｈｅｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｈｅＬＢＧＫｍｅｔｈｏｄａｎｄｔｈｅＮ－Ｓｅｑｕａｔｉｏｎｓｉｓｒｅａｖｅｌｅｄｔｈｒｏｕｇｈａｎｅｗｍｅｔｈｏｄ．ＩｃｌｅａｒｌｙｅｘｐｌａｉｎｅｄｔｈｅｒｅａｓｏｎａｂｌｅａｓｐｅｃｔａｎｄｔｈｅｐｒｏｂｌｅｍｏｆＬＢＧＫｍｅｔｈｏｄ．Ｆｉｎａｌｌｙ，ｔｈｅｓｔａｂｉｌｉｔｙｏｆＬＢＧＫｍｅｔｈｏｄｉｓａｎａｌｙｚｅｄａｎｄｔｈｅｎｅｗｍｅｔｈｏｄｔｏｉｍｐｒｏｖｅｔｈｅｓｔａｂｉｌｉｔｙｏｆＬＢＧＫｍｅｔｈｏｄｉｓｐｒｏｐｏｓｅｄ．Ｔｈａｔｉｓ，Ｔｈｅｅｑｕｉｌｉｂｒｉｕｍｄｉｓｔｒｉｂｕｔｉｏｎｆｕｎｃｔｉｏｎｔａｋｅｓｔｈｅｆｏｒｍｏｆｌｏｇａｒｉｔｈｍｉｎｏｒｄｅｒｔｏｉｎｓｕｒｅｔｈｅｅｑｕｉｌｉｂｒｉｕｍｄｉｓｔｒｉｂｕｔｉｏｎｆｕｎｃｔｉｏｎｂｅｉｎｇｔｈｅｐｏｓｉｔｉｖｅｎｕｍｂｅｒ，ｗｈｉｃｈＣａｎｉｍｐｒｏｖｅｔｈｅｓｔａｂｉｌｉｔｙｏｆＬＢＧＫｍｅｔｈｏｄ．Ｋｅｙｗｏｒｄｓ：ＬＢＧＫｍｅｔｈｏｄ，ｍｏｄｅｌ，ｓｔａｂｉｌｉｔｙⅡ武汉理工大学硕士学位论文第１章绪论１．１ＬＢＧＫ方法的发展史纵观古今，任何新生事物总是由旧事物发展变化而来，体现着唯物辨证法中亘古不变的发展与普遍联系原理。
常用元胞自动机
常用元胞自动机在元胞自动机是由空间上各项同性的一系列元胞所组成，是在有限元胞自动机基础上发展起来的，用于模拟和分析几何空间内的各种现象。
典型的元胞自动机在元胞自动机的发展过程中，科学家们构造了各种各样的元胞自动机模型。
其中，以下几个典型模型对元胞自动机的理论方法的研究起到了极大的推动作用，因此，它们又被认为是元胞自动机发展历程中的几个里程碑。
l. S. Wolfram和初等元胞自动机初等元胞自动机(Elementary Cellular Automata，简称ECA)是状态集S只有两个元素{s1，s2}，即状态个数k=2，邻居半径r=l的一维元胞自动机(谢惠民，1994、李才伟，1997、Wolfram，S，1986)。
它几乎是最简单的元胞自动机模型。
由于在S中具体采用什么符号并不重要，它可取{0，1}，{-l，1}，{静止，运动}，{黑，白}，{生，死}等等，这里重要的是S所含的符号个数，通常我们将其记为{0，1}。
此时，邻居集N的个数2r=2，局部映射f:S3→S可记为:其中变量有三个，每个变量取两个状态值，那么就有2×2×2=8种组合，只要给出在这八个自变量组合上的值，f就完全确定了。
例如以下映射便是其中的一个规则:通常这种规则也可表示为以下图形方式(黑色方块代表l，白色方块代表0):这样，对于任何一个一维的0，1序列，应用以上规则，可以产生下一时刻的相应的序列。
以下序列就是应用以上规则产生的:t: 010111110101011100010t+1:1010001010101010001以上八种组合分别对应0或1，因而这样的组合共有28=256种，即初等元胞自动机只可能有256种不同规则。
S. Wolfram定义由上述八种构形产生的八个结果组成一个二进制(注意高低位顺序)，如上可得01001100，然后计算它的十进制值R:R在[0，255]内，S. Wolfram定义R为初等元胞自动机的标号，则上面的元胞自动机模型就是76号初等元胞自动机(谢惠民，1994；李才伟，1997)。
格子玻尔兹曼方法
格子玻尔兹曼方法
格子玻尔兹曼方法是一种用于模拟气体动力学的计算方法。
它通过将模拟区域划分为若干个小的格子，然后在每个格子中模拟气体分子的运动状态，来预测整个气体系统的宏观性质。
在格子玻尔兹曼方法中，每个格子内部的气体分子相互作用被表示为碰撞模型。
根据碰撞模型的不同，格子玻尔兹曼方法可以分为两类：老化型和多尺度型。
在老化型格子玻尔兹曼方法中，每个格子内的气体分子与周围格子内的分子发生碰撞，从而实现能量和动量的传递。
通过不断迭代，系统中的能量和动量逐渐趋于平衡状态。
而在多尺度型格子玻尔兹曼方法中，每个格子内部的气体分子被分为多个不同尺度的子格子，每个子格子内的分子与周围子格子内的分子进行碰撞。
这种方法可以更加准确地模拟气体分子的运动，但计算复杂度也更高。
格子玻尔兹曼方法的一个重要应用是在计算流体力学中。
它可以用来研究气体的流动、传热等问题。
由于该方法可以充分考虑气体分子间的相互作用，因此在研究微观尺度上的气体行为时具有一定的优势。
总的来说，格子玻尔兹曼方法是一种能够模拟气体动力学的计算方法，它通过将模拟区域划分为小的格子，模拟分子之间的碰撞来预测气体系统的宏观性质。
这种方法在计算流体力学中有广泛的应用。
格子气模型在机场出港大厅旅客疏散过程应用
而且由于采
既保证了各行人具有同等的选择权
� 重制约经济发展的 "瓶颈 " 由于公共场所行人拥挤造成的各种事 � 利 又不用解决行人之间相互冲突的问题 所以使得编程十分简 故数量近年来一直呈上升的趋势 研究人流在公共场所的流动问 单 而且用唯一的参数来定义规则 使得模型具有灵活 易于调整 � � � 题已经越来越受到重视 近年来研究人流的方法主要有分子动力 的特点 很容易就可以将其应用到各种不同的环境中 � �
建立了偏向随机行走格子
将其应用到机场候机楼出港大厅 � � 并对出港大厅内的出港旅客流的形成 平衡 缓解和消散四个阶段进行
� 近年来 随着社会的飞速发展和人口的急剧 � 增长 对于人群 仿真 得到了一些出港旅客流运动性能指标的规律 和资源相对集中的公共场所而言 � 日益加剧的人流问题已成为严 用了随机的 串行更新规则
"
系统在每一个时步的状态由它上一个时 在每一时步中 所有的行人按 行人在每一次更
步的状态以及演化的规则来决定
照事 先安排好的随机次序轮流进行位置更新
新位 置时最多只能移动一个格点的距离 至于移动到哪一个格 � 点则由演化规则来决定 为了使各行人具有同等的选择权利 在 � 每一 个时步的开始 首先对处于系统中的所有行人进行随机编 号 以确定各行人在该时步的行走次序 � 则如下 � 关 于人 群在 公共 � 场 所的 运 动规 � 律研 究 已经 有许 多 成果 � 1)在每一时步的开始 如不受阻挡 旅客从各入口进入大厅
每人每时
型通道等 但将格子气模型运用到候机楼旅客出港大厅的旅客疏 步最多只能 走一格 �
不允许两个或两个以上的旅客同时占据一
散过程的研究还较少 本文基于� 个格点 旅客向着出口的方向运动 并可以以一定的概率转到旁 � -. /01023 . 等人提出的偏向随机 行走格子气模型中方向选择概率的计算方法和该模型的运行规 � � 边的两个格点中 但是不允许后退 图 2 列举了在偏向随机行走 则 建立串行更新的一速偏向随机行走格子气行人流模型 � 将该 � 规则下旅客所有可能的运动状态 其中黑色实心圆圈代表旅客 箭头表示旅客可能的运 模型应用到机场候机楼出港大厅的出港旅客流 � 并对候机楼出港 打叉处代表 该格点正被其他行人占据
一种考虑人数分布特性的人员疏散格子气模型
一
，
如 何针 对人 群特 征 建 立 准 确完 善 的 模 型 和规 则
是 研究 的关 键． 截 至 ２ ０ ０ ５年 初 ，已有 超 过 ５ ０种 人 员 疏散 模 型在各 类 文 献 中被 介 绍 ［ ． 目前 人 员 疏 散 １ ］ 模 型研究 大致 可 以分 为 两种 ：一 种是 将 人 看 作 连续
维普资讯
自美科手遗展 第１卷 第５ ２ ８ 月 ８ 期 ０ 年５ ０
一
种考 虑 人 数 分布 特 性 的 人 员 疏 散 格 子 气 模 型 ＊
宋卫 国 张 俊 胥 旋 刘 轩 于彦 飞
中 国科 学 技 术 大 学 火 灾 科 学 国家 重 点 实 验 室 ，合 肥 ２ ０ ２ ３０ ６
员 的个体 差异 以及 人 员 的行 为． 微 观模 型又 包 括连
续 性模 型和 离散模 型 ，连续 性 模 型 中最 为著 名 的是 Ｈｅｂｎ ｌｉｇ等 提 出 的社 会 力 模 型 ［ ，离 散 模 型 以 格 ２ ］ 子 气模 型 （ Ｇ）７１ 和 元胞 自动 机 （ Ａ）１ ｚ 使 用 较 Ｌ ［ ７ ］ Ｃ ［ １ ８ ］
摘要
紧 急情况 下 （ 火灾等 ） 如 的人 员疏 散 研 究 可 以 为减 少 人 员伤 亡 提供 科 学指 导 ， 同 时也 可 以 为
建 筑设 计提供 实用 的依据 和方 法． 但 是 在 实 际 的疏 散 分 析 中 ，往 往 不 知 道准 确 的待 疏散 人 数 ，而 只
格子气模型
由于人员有偏随机运动，且有8个移动方向，因 此需要将优先方向D投影到各个方向上
人员之间的相互作用的强度用F 表示，不 同方向上大小不同
按照移动规则以及影响因素，建立如下的数学模型
✓ 为了保证各个方向上的概率之和为1，设置归一化因子N ✓ i, j 表示该方向是否被其他个体占据，取0、1两个值 ✓ Ii, j 表示是否沿上一时间步的方向走，>1表示是，=1表示没 有这 一性质，每一步均重新选择方向
规则：
1、每个个体占据多个格子 2、每个时间步长移动一格 3、允许相邻个体之间的压盖 4、规定8个移动方向
场景：方形的大厅和一个 宽度为W的出口，人员有 偏随机运动
在人群疏散过程中不可避 免地存在着人员之间的相 互作用力以及人和建筑物 之间的相互作用力。将个 体当作有质量无体积的质 点则无法准确描述这种作 用力。
数学模型的建立方法 从前面提到的几种模型的建立方法可以看 出，基本上是从场景着手，定义个体的移 动规则，以及更新规则。按照规则建立相 应的数学模型。
如果模拟场景有多个出口，怎么建模？ 如果考虑到人员的从众效应，怎么建模？
✓ Di, j 表示优先方向D沿(i,j)方向的分量
✓ fi, j 、 fw 分别表示人与人之间的作用力以及人与墙之间的
P
w
作用力
规则：每个个体占据两个 格子，有四个可能的移动 方向
每人可能会处在如右图所 示的9种状况下
场景：在某些紧急情况下， 例如地震、在某些顶比较低 的通道里疏散等等，人们不 大可能正常地行走。考虑人 们爬行的状态，需要重新建 立模型来模拟。
基本格子气模型（Lattice-Gas Automata） 是元胞自动机的具体化，利用元胞自动机 的动态特征来模拟流体粒子的运动。
应用d2q9格子玻尔兹曼模型
对高雷诺数流动的模拟
该模型在模拟高雷诺数流动时可能存 在局限性，需要采用更精细的离散方 法。
05
模型优化与改进建议
算法优化
采用快速格子玻尔兹曼方法
通过采用快速格子玻尔兹曼方法，可以降低算法的时间复杂度， 提高模型的计算效率。
优化边界条件处理
针对边界条件的处理进行优化，例如采用周期性边界条件或者采用 更加高效的边界处理方法，可以减少计算量和内存占用。
意义
通过研究d2q9格子玻尔兹曼模型，可以更好地理解流体流动和传热过程的微观 机制，为实际工程应用提供理论支持和指导。
研究现状与发展
现状
目前，d2q9格子玻尔兹曼模型已经在多个领域得到广泛应用，如航空航天、能 源、材料科学等。然而，该模型在某些复杂流动和传热环境下仍存在一定的局限 性。
发展
随着计算机技术的进步和数值模拟方法的不断完善，d2q9格子玻尔兹曼模型将 有望在更多领域发挥其优势，同时其应用范围也将不断扩大。
误差，提高模型的精度。
并行计算实现
01
02
03
分布式内存并行化
通过分布式内存并行化， 可以利用多台计算机的内 存资源，提高模型的计算 速度和效率。
并行化碰撞处理
通过并行化碰撞处理，可 以利用多核CPU的计算资 源，提高模型的计算速度 和效率。
并行化边界处理
通过并行化边界处理，可 以利用多核CPU的计算资 源，提高模型的计算速度 和效率。
D2Q9模型的数值实现
D2Q9模型的数值实现主要包括以下几个步骤
碰撞过程：根据分子速度和位置进行碰撞，更新速度分 布函数。
时间步进：根据设定的时间步长重复进行碰撞和迁移过 程，直到达到稳定状态或达到预设的计算时间。
·
路网的警力布控作者：周婧琳来源：《法制博览》2013年第09期【摘要】本文基于C#语言编写的计算机软件和既得数据（历史经验数据和即时路况数据），基于格子气模型对数据进行分析处理，从而实现对道路上需要追踪的运动车辆轨迹进行研判，得出对其实施追捕或拦截的路线，指导应急条件下的警力指挥和调度。
研究成果将为公安实战部门提供有力的决策支持，并可以应用到犯罪追捕、消防、医疗、各类应急抢险等各个应急行业中。
在此成果的基础上加以改进后可应用到智能交通领域，建立数学模型，以指导城市道路规划和犯罪辑逃，解决当前日益突出的城市交通犯罪问题。
【关键词】路网；警力；格子气；仿真一、引言随着我国城市化进程的加快，城市的道路建设获得了突飞猛进的发展。
截止2010年底，全国公路通车总里程超过400万公里，在“十一五”计划的五年里新增公路65万公里[1]而且，现阶段，全国的车辆保有量已超过7000万辆，20个城市的汽车保有量超过百万，且部分二三线城市的车辆保有量还会出现快速增长[2]。
这些条件都直接导致了犯罪分子驾车沿路网逃逸的可能性大大增加。
在目前的研究中，警力指挥调度一直是一个研究热点。
在国外，警力指挥调度一般以成熟的业务机制为基础，通过情报分析研判来辅助有经验的管理人员来做出决策。
例如最先开始于美国、现在在西方发达国家广为应用的CompStat模式就是一种警力指挥调度模式[3]。
该模式以计算机信息分析为基础，通过对大量的情报信息进行综合研判，提出当前犯罪问题的主要热点区域分布，并对一段时间以内的犯罪问题进行预测，随后决策人员根据现有的警力分布来做出针对性的部署。
这种指挥调度模式主要以专家经验为基础，优点是反应快，一旦发现问题可以立即作出响应，但也存在着可靠性不高、不利于警员的广泛参与等问题[4]。
本文基于格子气模型和软件仿真，对城市路网中车辆在路段上所用时间和对在路口处对不同行驶方向的选择的分析讨论，结合实际路网运行状况和路网拓扑结构图，对城市道路上车辆的运行路线进行初步的研究，为城市路网上的警力调度提供科学依据。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
i, j i, j i, j
f 、 f 分别表示人与人之间的作用力以及人与墙之间的 作用力
i, j w P w
场景：在某些紧急情况下， 例如地震、在某些顶比较低 的通道里疏散等等，人们不 大可能正常地行走。考虑人 们爬行的状态，需要重新建 立模型来模拟。
规则：每个个体占据两个 格子，有四个可能的移动 方向 每人可能会处在如右图所 示的9种状况下
由于人员有偏随机运动，且有8个移动方向，因 此需要将优先方向D投影到各个方向上
人员之间的相互作用的强度用F 表示，不 同方向上大小不同
按照移动规则以及影响因素，建立如下的数学模型
为了保证各个方向上的概率之和为1，设置归一化因子N 表示该方向是否被其他个体占据，取0、1两个值 I 表示是否沿上一时间步的方向走，>1表示是，=1表示没 有这 一性质，每一步均重新选择方向 D 表示优先方向D沿(i,j)方向的分量
基本格子气模型（Lattice-Gas Automata） 是元胞自动机的具体化，利用元胞自动机 的动态特征来模拟流体粒子的运动。
个体有质量、无体积 只能在网格点上存在， 并沿着网格线在网格间 运动 在同一时刻同一点上， 沿着每一网格线运动方 向最多只有一个粒子
建模思路 1、移动规则 2、更新规则
作为离散模型，二者均有很高的计算效率 相对于元胞自动机，格子气模型也是从微观上进 行模拟群体的疏散过程，但是在模拟疏散时更关 注群体疏散中表现出来的物理上的流体特性 二者均关注整个疏散过程所耗费的时间，而不关 注单独的个体的疏散行为。 在具体的建模思想上，二者有很大的共性。但是 元胞自动机模型中，每一个元胞的状态是由周围 邻居的上一时间步长状态决定的。格子气模型中， 每一个粒子下一刻的位置是由上一时间步长的各 个可移动方向的方向概率决定的。
Dx = D|x − x0|/(|x − x0| + |y − y0|) Dy = D|y − y0|/(|x − x0| + |y − y0|)
其中，DX+DY=D，保证了每种情况下向各个方向移动的总概率是1。
更新规则： 计算好每个个体向各个方向移动的概率后 ，按照上述移动规则，依照随机 顺序对每个个体的位置进行更新
8种可能存在的情况 由于场景以及逃生规则，采用了 偏向随机且无退后的模型，因此 前进的方向只能是三个
移动规则：对于该图中在出口左边的个体,对于8种可能的情况有如下移动规 则
A、pt;y = Dy + (1 − D)/3; pt; x = Dx + (1 − D)/3; pt;−x = (1 − D)/3 B、pt;y = Dy + (1 − D)/2; pt; x = Dx + (1 − D)/2; pt;−x = 0 C、pt;y = D + (1 − D)/2; pt; x = 0; pt;−x = (1 − D)/2 D、pt;y =0; pt; x = D + (1 − D)/2; pt;−x = (1 − D)/2 E、pt;y =1; pt; x = 0; pt;−x = 0 F、pt;y =0; pt; x = 0; pt;−x = 1 G、pt;y =0; pt; x = 1; pt;−x = 0 H、pt;y = pt; x = pt;−x = 0 D是优先方向偏移量,指向出口处，其值不变，但是方向随着运动物体的位置不断改变。 DX、DY分别是D在X方向上的偏移分量和在Y方向上的偏移分量。
改进模型（3种） 按照格网划分方式不一样区分 1、按照人所占的空间范围划分，即一个人占一小格，每个时间步 长，每个人移动一个格子的距离 2、为了研究人群之间的相互作用，需要划分精细格网 3、为了模拟人的匍匐前进，将个体看做是长条形，即每个个体占 两个格网空间
场景：方形的大厅和一个宽度 为W的出口。在这种场景中人 的运动规律一般是从各个角落 涌入出口。因此在这种情况下 人们会有一个主方向，即无论 从任何位置均朝向出口，并且 不会后退。
规则： 1、每个个体占据多个格子 2、每个时间步长移动一格 3、允许相邻个体之间的压盖 4、规定8个移动方向
场景：方形的大厅和一个 宽度为W的出口，人员有 偏随机运动 在人群疏散过程中不可避 免地存在着人员之间的相 互作用力以及人和建筑物 之间的相互作用力。将个 体当作有质量无体积的质 点则无法准确描述这种作 用力。