Lattice+Boltzmann方法及其在枝晶生长模拟中的应用研究.pdfVIP

Lattice Boltzmann 方法及其在枝晶生长模拟中的应用 苏 敏 孙东科 朱鸣芳 东南大学 材料科学与工程学院 摘要：本文介绍了Lattice Boltzmann 方法（LBM ）的基本原理，模拟流场和浓度场的LBM 控制方 程，讨论了不同情况下边界条件的处理方法。将LBM 和一个改进 cellular automaton(MCA)模型相耦 合，模拟了合金中单枝晶在纯扩散和存在对流的两种条件下的生长行为。结果表明在纯扩散条件下， 枝晶呈现对称生长形貌。但存在对流时，枝晶的生长在上游方向得到促进，而在下游方向受到抑制， 形成了非对称的枝晶生长形貌。 关键词：Lattice Boltzmann 方法；流场；浓度场；边界条件；凝固；枝晶生长 Lattice Boltzmann 方法(LBM)是上世纪80 年代中期提出研究流体流动的一种新的数值计算方法， 其基本原理来源于格子气自动机方法(Lattice gas automata) [1~4] 。LBM 不仅适合于大规模的并行计算， 而且可以处理任意复杂的边界，其计算速度同空间复杂程度无关。由于 LBM 具有上述优点，在近 20 年内得到了迅速发展，现已成功地应用于流体力学诸多方面的数值计算，如湍流、多相流、多孔 介质流、空气动力学、磁流体力学、燃烧问题以及一些化学反应和具有复杂边界的流动等[1~7] 。近年 来LBM 在凝固过程数值模拟领域如对流作用下的枝晶生长、金属熔化和泡沫金属成型过程等方面也 得到了初步应用，显示出很好的应用潜力[1] 。 1 Lattice Boltzmann 方法模型 1.1 计算流场的LBM 模型和控制方程 LBM 用f (x ,t) 表示t 时刻在位置x 方向为i (具有速度e ) 的粒子分布函数。宏观密度ρ ，动量 i i ρu 可用如下式子表示： ρ ∑f i (x ,t) （1） ρu ∑e f (x ,t) （2 ） i i 这样分布函数就满足方程 ： f (x +e ∆t,t +∆t) −f (x ,t) ∆tΩ （3 ） i i i i 式中，Ω 为碰撞项。式（3 ）即为Lattice Boltzmann 方程。陈十一等人基于单一松弛时间模型对碰 i 撞函数进行简化[3,7] ： 1 Ω − (f (x ,t) −f eq (x ,t)) （4 ） i i i τ 基金项目: 国家自然科学基金资助项目 第一作者简介：苏敏，女，1979 年出生，硕士研究生，东南大学材料科学与工程学院，南京(210096) 通讯作者：朱鸣芳，教授，东南大学材料科学与工程学院，Email: zhumf@seu.edu.cn.