地下矿巷道火灾风流紊乱现象仿真方法研究-采矿工程专业论文.docx

西安建筑科技大学硕士学位论文 西安建筑科技大学硕士学位论文 PAGE PAGE 3 进了矿井火灾风流的计算机模拟的发展 ，为创建高效、可靠、适应性强的蔓延模 型提供了理论支持。格子 Boltzmann 方法!Il (Lattice Boltzmann Method ，LBM) 作 为一种新的流体计算方法，因为具有物理定义清晰、算法简单、边界条件易处理、 压力直接计算 、本质并行等显著优点，引起了众多领域的学者对其理论和应用进 行研究。目前， LBM 在热流动、端流、多相流、化学反应扩散流、生物流场、掺 流、粒子悬浑流等众多领域中都得到了比较成功的应用。 运用格子 Boltzmann 方法研究地下矿井火灾风流流动规律与能量传递过程的优 势主要表现在以下几方面 : (1)格子 Boltzmann 方法从离散运动论层次描述物理规律，克服传统的计算 流体力学无法解算多数微分方程组的难题: (2)传统的计算流体力学对矿井火灾的研究常常假设火灾巷道为水 平圆管巷 道，格子 Boltzmann 方法则可以通过对边界的分析在任意形状的巷道进行 : (3) 格子 Boltzmann 方程中压力以状态方程表示，因其解法 无需反复计算而 使其适应于大规模并行的计算: (4)精度要求较低时，模型的边界条件容易设定。 1.3 格子 Boltzm ann 方法 格子 Boltzmann 方法历史上是由格 子气 自动机 (LGA) 发展而来的，因此与 LGA 相似，格子 Boltzmann 方法也可以看作 是一种流体运动的微观模型:另 一方 面，格子 Bol刷刷方法可以看作是连续 Boltzmann 方程的一个特殊的离散格式， 因此格子 Boltzmann 方法又可以看作是 -个介观的动理论格式。 1.3.1 计算流体力学研究的两种方法 在流体力学中，流体运动论层次的描述所遵循的基本原则是物理守恒定律，它 要求所有微观过程保持质量守恒、?;)j量守恒以及能量守恒。离散运动论层次想象 流体是由大量的微观粒 子所组成的，这些粒子的运动遵守力学 定律，并且服从统 计学规律。 ( 1) 基于宏观连续模型的自顶向下方法 该方法的描述是以场论作为工 具，假设流体连续的分布于流体所占的整个区 域，其他物理量均是空间 与时间的足够光滑的函数 ，从基本的质量守恒、?;)J量守 忡。恒和z能量d守(恒 以及空间不变性导出流体的一组偏微分方程，即: 忡。 θ(ρu) 一百一+ V.(ρuu)=v. σ+ρg (1.1) 一-δE(pe) +V.(ρue)= σ:如-V.q 一 虽然这些方程形式简单，但只有一些很简单的问题才可以求得解析解 I21，实际 应用起来十分困难 。传统的计算流体力学方法就是在宏观连续模型 N-S 方程的基 础土，运用一些数值离散方法进行求解。 (2) 基于微观离散模型的白底向上方法 该方法的描述是以非平衡统计力学为出发点，认为流体是由大量的微观粒子所 组成，并且这些粒子的运动遵守力学定律以及统计学规律。得到的基本方程是 Boltzmann 方程，即: 宇+V Vf=(」)(f -f 呵) (1.2 ) OI ? 该方程描述的是稀薄气体分子的动力学规律。值得注意的是模型中的流体是由 大量离散的分子组成，而不再是连续介质。自底向上方法从分子的观点对流体进 行离散化处理，但是空间和时间依然连续。这种方法不需要求解复杂的 偏微分方 程组，完全应用简单的局部动力行为展现，打破了传统的数学建模的观念，为模 拟大规模复杂现象提供可行途径，也为可视化提供了新的方法。 1.3.2 格子 Boltzmann 方法的发展 (1)格子气自动机 格子气自动机基于场论对流体进行空间和时间的离散，将流体存在的空间划分 为离散网格，将时间离散为整时间步，流体假设为只有质量而没有体积的微粒。 粒子存在于网格节点之上并沿网格线运动，流场所有的粒子依据所制定的碰撞规 则随时间步长同步地相互碰撞和迁移。同步意味着某一时刻所有粒 子在网格节点 处发生碰撞后，下 一时刻所有粒子正好运动到临近的网格节点并在此与到达该节 点的其它粒子碰撞。 根据以上规定，就可以建立一个格子气自动机，即: N;(x+c;l!J，1+!li)= Nj(x，t)+0.(n) (1.3) 式中 N，衍，t)表示 t 时刻Jx 节点上以 Ci 速度运动的粒子数，是个整数 ;[Ji(的表示 碰撞函数。 格子气自动机求解流体力学问题步骤 z 第二 ，建立格子气自动机。包括流体离散，空间离散，时间离散和碰撞函数: 第二，运行格子气自动机。流体柿子在离散网格上根据模型建立的碰撞规则随 时间步长演化; 第三，应用统计学规律求数值结果。 (2) 格子 Boltzmann 方法 格子 Boltzm