强磁场下基于热格子Boltzmann方法的自然对流数值模拟.pdfVIP

第六届两北地区计算物理学术会议论文集 强磁场下基于热格子Boltzmann方法的自然对流数值模拟 钟诚文·解建飞卓丛山 周小旭熊生伟 (两北T业大学哭型、叶栅空气动力学国防科技重点实验室，陕西西安710072) 格式的对流扩散方程进行了改进，形成了可以模拟非等温流动的热格子Boltzmann方法。本文用这种方法对强磁 场环境下低导电溶液的方腔自然对流进行丫数值模拟。在强磁场环境下，方腔中的低导电溶液的流动除了受到浮 力外，还受到磁化力和洛仑兹力，自然对流会受到很人的抑制。我们对Gr--4．64×104和Pr=7．0的低导电溶液 (11wt％KCl溶液)柚i强磁场环境下的自然对流进行了数值模拟。详细分析了自然对流在强磁场环境下受到的影响， 并与没有J,bDii磁场的数据进行了定量比较。本文的工作提供r模拟强磁场环境下自然对流的一种新的数值方法。 关键字：热latticeBoltzmann方法，磁场，自然对流，方腔 引言 分方程数值解法【ll，并且可以利用并行计算机模拟输运现象。格子气(LG)米源于连续力学的 微观观点。在微观方面，粒子在不同的方向以不同的速率移动，这些粒子总的效应就是压强和温 度，在宏观的角度可以被测量。因此，通过模拟每个粒子的运动，我们就可以模拟宏观的流体行 为。在格子气中，粒子在每个时间步只能移动到与其相邻的格点上，所移动到的邻点由速度的方 向来决定。在这个过程中，它们可以和其它的粒子碰撞，并且交换它们的动量。格子Boltzmann 方法是建立在格子气元胞自动机的基础上的，所以它们用相同的格子和应用相同的碰撞过程。与 格子气元胞自动机不同的是，格子Boltzmann方法用连续分布函数代替粒子来处理碰撞，这些分 布函数在当地互相作用(只是在某一点分布函数被引入时)，并且在碰撞后迁移到相邻格点。格 子Boltzmann方法后续的发展是对碰撞算子的简化和选取不同类型的分布函数。 格子Boltzmann方法是由格子气自动机发展而来的，因此从格子气自动机那里继承了一些优 于传统方法的优判21。由于粒子之间碰撞的当地性和粒子信息的输运，格子Boltzmann方法具有 天然并行性，很适合大型的并行计算。格子Boltzmann方法的数值精确度，数值鲁棒性，对于复 杂边界的适应性，计算效率，也证明其可以作为一种新型的计算流体力学(CFD)数值方法。 本文在J．A．Somers提出的格子Boltzmann方法pJ的基础上，加入格子Boltzmann离散格式 的温度扩散方程141，提出一种研究温度场演化过样的热格子Boltzmann方法和一个比较实用的边 界条件，用来模拟1F等温流动。并给出了强磁场环境下模拟盐溶液的自然对流的详细过程和定量 结果，并将结果和没有外加磁场的情况作了定量的比较。 1．数值方法 1．1格子Boltzmann方程 的演化包括两步：迁移步，拖曳所有变量以使在位置x的质量密度M移动到x+c，；碰撞步，在 每个当地格点的速度方向重新分布质量密度。我们用这个格式求解卜．列耦合的偏微分方程组： ‘本文受如下项日支持：国家自然科学攮金项IE!C0 工业人学基础科学研究基金重点项目。 第六届两J匕地区计算物理学术会议论文集 2008 ＼／ ／＼ 图1D2Q9多速度模型 色M+c，·叫=Q。(N) (1) 其中，碰撞算子Q，(N)遵守质量和动餐守恒定律：∑。Q。(N)=o，∑，c。Q，(N)--f。 力f(x,tJ是体积力，包含重力和电磁力，会在后面详细讨论。质量密度平衡分布Ⅳj叫【x，f)： V·u]}(2) ％：{4，4，1，4，1，4，1，4，1}，f：0，．．．，8，y是运动粘性系数。 为了确定碰撞算子Qr(N)，我们考虑格子Bolt砷ann方程的交错格式： Ⅳl(x+扣，+三)=M(x一扣f一三)崛(N)(3) 把M(x+三c∥+吾)和弼(x一三c；,t--去)在Ⅳj扛，f)进行泰勒级数展开得到： q(N)=M(x+知，+三)一