格子Boltzmann方法模拟融化相变过程研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093084 格子Boltzmann 方法模拟融化相变过程 1,2 1 1 杲东彦 ，陈振乾 ，施明恒 （1.东南大学能源与环境学院，江苏 南京 210096 ； 2. 南京工程学院能源与动力工程学院，江苏 南京 211167） （Tel: E-mail:yzgdyg96@） 摘 要 ：建立基于Boussinesq 假设耦合的格子Boltzmann 双分布函数模型，采用焓法迭代求解相变 非线性源项，并通过变松弛时间参数方法处理固液两相变热物性问题。热传导融化过程的计算值与分 析解的对比分析说明该模型能够准确地模拟融化过程。自然对流条件下融化过程的模拟结果表明自然 对流在一定程度上影响了融化界面形状特征及温度分布等，体现出与热传导融化不同的物理机制。 关键词：格子-Boltzmann；融化；数值模拟；自然对流 0 前言 融化是自然界和工程技术领域中的一种常见现象，例如冰的融化、晶体的熔融生长、 相变储能等。由于融化过程存在移动界面的特点及固液界面处伴随有潜热的吸收过程， 因此融化传热问题表现出强非线性的特点，求解十分困难。而液相自然对流等因素存在 使得融化问题更加复杂。融化传热问题的求解除少数一维理想化问题可以有准确的解析 解之外，大多数问题需要数值求解。由于移动界面问题，数值求解需要特殊方法，如前 沿追踪方法、自适应网格方法、相场方法、VOF 方法等。 格子-Boltzmann 方法（简称LB）是近十几年来发展起来一种不同于传统数值计算方 法的流体计算和建模新方法。该方法在许多传统计算方法难以胜任的领域，如微尺度流 动和换热、多孔介质、多相流等，都取得了成功的应用[1-3] 。近年来也开展了应用格子 -Boltzmann 方法研究固液相变传热问题，由于问题复杂性，研究初期大多是针对热传导 型的融化问题，Jiaung 等应用格子-Boltzmann 方法计算了固液相变的热传导型单区域、 双区域、三区域问题[4] ；钱吉裕等应用D2Q5 热传导模型分析了多孔介质内固液相变问题 [5] 。Christian 等采用D2Q5 温度模型耦合D2Q9 速度场模型，分析了自然对流情况下融 化传热问题[6] 。Christian 在计算中采用固液两相热扩散系数相同的假设，然而固液两 相热物性通常有所不同，如热扩散系数冰是水的三倍左右等。针对自然对流条件下的融 化相变传热问题，本文建立基于 Boussinesq 假设的耦合温度场与速度场的双分布函数 格子-Boltzmann 模型，并结合Jiaung 基于焓的方法迭代求解相变非线性项[4]及通过变 松弛时间参数方法处理固液两相变热物性问题。 1 计算模型 1.1 融化相变传热的宏观方程 针对不可压缩流体，忽略流动中黏性热耗散，采用Boussinesq 假设的自然对流条件 下的融化相变流动及传热方程组为： ∇⋅u 0 （1） ∂u 2 (u )u P υ u gβ(T T ) （2 ） + ⋅∇ =−∇ + ∇ − − 0 ∂t 基金项目：国家自然科学基金资助项目（No） ∂T 2 L ∂(f l )