自适应直角网格下模拟复杂几何界面的MHD流动算法发展.pdfVIP

中国工程热物理学会 传热与传质 学术会议论文 编号：123186 自适应直角网格下模拟复杂几何界面的 MHD流动算法发展 张杰，倪明玖* 中国科学院研究生院物理科学学院，北京，100049 * Tel : 010, E-mail: mjni@ 摘要：本文将相容守恒格式扩展到自适应（AMR）网格下进行磁流体流动模拟，保证了不同层网格间的 相容性和守恒性，并利用直角网格处理复杂几何界面，大大节省了计算机资源。文章结合切割网格法 和VOF方法，使得边界条件能精确施加在直角网格上，实现了复杂构型下MHD流的精确模拟。计算结 果表明，切割网格法能在直角网格下精确模拟有复杂边界的磁流体流动。 关键词：MHD流动，复杂边界，直角网格，VOF方法，相容守恒性 0引言 磁流体力学（magnetohydrodynamics ，MHD ）是流体力学的一个分支，它主要研究 在磁场作用下流体呈现出的不同于普通流动的状态[1,2]，比如高压降，变化的速度分布， 传热效应的增强或者削弱等。而在核聚变堆液态自冷包层中，由于存在环形流动或者突 扩流动等更复杂流动状态，因而能否精确模拟复杂几何形状的MHD 流动成为数值模拟 研究的关键[3]。 在这种情况下，处理复杂几何构型流动的一个主要挑战就是如何模拟复杂边界。在 传统方法中，非结构网格或者贴体网格是处理这类问题的通常手段，因为边界条件加载 到网格上很容易[4,5,6]。Ni 等[4]曾用他们发展的非结构网格下的相容守恒格式计算了圆 管中的MHD 流动，Vantieghem[5] 也曾借鉴这种方法计算了导电壁面的圆管流动， Mutschke[6]在模拟MHD 效应下的圆柱绕流时也采用了相同的网格生成方式。 最近在流体力学数值模拟中，笛卡尔直角网格越来越多的被用到处理复杂几何界面 流动上，这得益于它简单的网格生成方式，不需要任何坐标变换，大大节省了计算机资 源。现在比较常见的基于直角网格的方法包括嵌入边界法（IBM ）[7], 虚拟网格法 [GCM][8]，切割网格法[CCM][9]等。这种类处理技术中，整个N-S 控制方程的离散仍然 基于直角坐标系，而并不是贴合住复杂几何界面的网格，因而造成这种方法的主要难点 就是如何精确加载边界条件。 另外在MHD 流动中，边界附近会形成很薄的“边界层”，根据与磁场方向的位置关 *基金项目：国家自然科学基金项科技部ITER 专项2009GB104001 1 1 系，分别称之为Ha 层和侧层，其无量纲厚度分别为 和 ，其中Ha 数是MHD Ha Ha 流动中表征磁场力和粘性力比值的无量纲参数。在极薄的Ha 层内有大量电流流过，需 要进行精确模拟，因而要布置多个网格，而在主流区域则只需要少量网格。本文采用自 适应网格分布，在Ha 层和侧层处自动加密，大大提高了计算效率。 在本文中，首次提出用切割网格技术在直角网格下处理复杂几何界面的MHD 流动， 并用自适应网格对壁面处和边界层处进行精确捕捉。在均匀不带壁面的情况下，本文和 文献[10]中提出的格式完全一致，具有相容守恒性质，并在此基础上推广到自适应网格， 能精确计算电流和洛伦兹力。 1数值方法 1.1控制方程 在有外加磁场时，粘性不可压缩流体的控制方程主要包括质量守恒和动量守恒方 程： ∇⋅u 0 (1) ∂u ρ( ) p