交错网格下MHD相容守恒格式的发展.PDFVIP

第 31 卷 第 3 期 核 聚 变 与 等 离 子 体 物 理 Vol.31, No.3 2 0 1 1 年 9 月 Nuclear Fusion and Plasma Physics Sep. 2011 文章编号：0254−6086(2011)03−0219−05 交错网格下 MHD 相容守恒格式的发展 李俊峰，倪明玖 (中国科学院研究生院物理科学学院，北京 100049) 摘 要：在低磁场雷诺数条件下，基于电势泊松方程，发展了交错网格下可以精确计算电流和洛伦兹力( 电磁 力)的相容守恒格式。采用压力为变量的原始变量法求解不可压缩 Navier-Stokes 方程，所计算的电流满足电荷守恒 定律，所计算的电磁力满足动量守恒定律。对金属流体在 Hartmann 数 50~5000 范围内验证了格式的精确性。交错 网格下相容守恒格式的发展为后续 MHD 稳定性分析、湍流的大涡模拟及直接数值模拟提供很好的选择。 关键词：磁流体动力学；相容守恒格式；交错网格；洛伦兹力 中图分类号：O361.3 文献标识码：A 1 引言 小，基于 Ohm 定律和电荷守恒定律可以精确地描 强磁场作用下的金属流体流动已经成为热核 述其 MHD 效应。与一般流体流动相比，金属流体 聚变研究领域非常热门的课题，是液态锂铅包层和 MHD 模拟的最大困难是如何保证计算过程中电荷 偏滤器研发的核心科学问题之一[1] 。研究强磁场作 守恒及电磁力计算的精确性。Sterl 在直角均分交错 用下金属流体的流动与传热性能是保证热核聚变 网格上对金属流体流动进行了数值模拟[4] ，但在很 液态包层设计成功与否的关键所在。 薄的 Hartmann 层内需要至少 2 ～3 个网格，极大地 金属流体流动过程中受磁场作用，会产生感应 增加了计算的网格量，在现有的计算机条件下，需 电流，其流动会受到洛伦兹力的影响，与一般流体 要较长的计算时间。由于 Hartmann 数越大，要求 流动会有较大差别。二维直通道内充分发展的磁流 网格越细，所以 Sterl 在计算中，对于二维问题 体力学问题已经得到很好的研究，获得了一些特定 Hartmann 数限制在 1000 以内，对于三维问题 情况下的解析解[2, 3] ，这对分析金属流体的磁流体 Hartmann 数限制在 50 以内，实际包层设计中， 4 5 动力学(MHD)效应很有价值。无量纲参数 Hartmann Hartmann 数高达 10 ~10 ，需要采用非均分网格来 数的平方表示劳伦兹力和粘性力的比值。在三维直 提高计算效率。完全交错网格下 Leboucher [5]发展 通道金属流体流动中，与磁场垂直方向的边界层叫 的格式可以实现强磁场下 MHD 的精确计算，但是 −1 仅限于单方向应用磁场及简单几何构型。倪明玖等 Hartmann 层，其无量纲厚度用 Ha 来度量；与磁 场平行方向的边界层叫侧层，其无量纲厚度用