磁气流体力学（MHD）
磁气流体力学（MHD）的理论基础
磁气流体力学的基础
教授，磁气流体力学（MHD）是流体力学和电磁学的结合吗？
没错。MHD是处理导电性流体（液态金属、等离子体、电解质溶液等）中电磁场与流体运动相互作用的学科。电流流经流体时，磁场作用会产生洛伦兹力，流动随之变化。反过来，导电性流体的运动会诱导磁场。
实际应用案例：
- 连续铸造：溶钢的电磁制动、电磁搅拌
- 铝电解精炼：Hall-Herault槽内熔融铝的流动控制
- 核聚变炉：等离子体的磁场约束
- MHD泵：无可动部件地传送液态金属
- 太空推进：MPD推进器
MHD的支配方程
MHD是Navier-Stokes方程与Maxwell方程的耦合体系。
修正Navier-Stokes方程（添加洛伦兹力）：
其中 $\mathbf{J} \times \mathbf{B}$ 是洛伦兹力（体积力）。
欧姆定律（运动导电流体中）：
$\sigma$ 是电导率，$\mathbf{E}$ 是电场，$\mathbf{u} \times \mathbf{B}$ 是流体运动引起的感应电动势。
磁场诱导方程：
右侧第一项是磁场的平流（冻结），第二项是磁场的扩散。
流体运动方程中添加了电磁力，同时磁场方程中也包含流速。这是完全的双向耦合。
MHD的无次元数
让我们整理MHD特有的无次元数。
| 无次元数 | 定义 | 物理意义 |
|---|---|---|
| 哈特曼数 Ha | $BL\sqrt{\sigma/\mu}$ | 电磁力/粘性力 |
| 磁雷诺数 Rm | $\mu_0 \sigma U L$ | 磁场平流/磁场扩散 |
| Stuart数（相互作用参数）N | $\sigma B^2 L / (\rho U) = \text{Ha}^2/\text{Re}$ | 电磁力/惯性力 |
| 磁Prandtl数 Pm | $\mu_0 \sigma \nu$ | 运动量扩散/磁场扩散 |
工业液态金属（溶钢、熔融铝等）中 $\text{Pm} \sim 10^{-6}$ 极小。这意味着磁场扩散远快于流体的运动量扩散，$\text{Rm} \ll 1$。此时磁场由外部印加场基本决定，只需考虑其对流体的影响（低Rm近似）。
哈特曼数越大，磁场影响越强，对吗？
是的。Ha = 0时是普通流体力学，Ha → ∞时流动完全受磁场约束。连续铸造中Ha通常在100-1000范围内。
核聚变炉的等离子体约束——MHD在人类梦想中的角色
ITER（国际核聚变实验炉）中，1亿度的等离子体通过托卡马克型磁场约束。等离子体稳定性分析的基础就是MHD理论。等离子体是否会"逃逸"磁场？是否会产生"扭结不稳定性"或"气球不稳定性"？——这些问题通过求解MHD方程的特征值问题来解决，从而优化线圈设计和电流分布。核聚变虽然声称"30年后实现"已多年，这正反映了MHD稳定性分析的复杂性。
磁气流体力学（MHD）的数值计算方法
MHD的数值解法
MHD耦合方程在CFD中如何求解？
低Rm近似（大多数工业液态金属）和高Rm问题（等离子体、天体物理）的方法不同。
低Rm近似的解法
当 $\text{Rm} \ll 1$ 时，诱导磁场相对于印加磁场可忽略，磁场可视为已知的外部场 $\mathbf{B}_0$。电场的支配方程为：
$\phi$ 是电位。求解此泊松方程得到 $\mathbf{E} = -\nabla\phi$，再计算 $\mathbf{J} = \sigma(-\nabla\phi + \mathbf{u} \times \mathbf{B}_0)$、$\mathbf{F}_L = \mathbf{J} \times \mathbf{B}_0$，并作为源项加入N-S方程。
计算步骤：
1. 从流场计算 $\mathbf{u} \times \mathbf{B}_0$
2. 求解电位的泊松方程
3. 计算电流密度 $\mathbf{J}$
4. 将洛伦兹力 $\mathbf{F}_L$ 加入N-S方程以更新流场
5. 迭代到收敛
只是增加了一个电位方程，基本上是N-S的扩展。
高Rm问题的解法
当 $\text{Rm} \gg 1$（等离子体物理、天体物理）时，必须完整求解磁场诱导方程。
主要数值方法：
- 约束传输法（Constrained Transport）：通过单元面上的磁通量守恒来严格保证 $\nabla \cdot \mathbf{B} = 0$
- 发散清理法（Divergence Cleaning）：加入补正方程以阻尼 $\nabla \cdot \mathbf{B}$ 误差
- 矢量位法：用 $\mathbf{B} = \nabla \times \mathbf{A}$ 自动保证 $\nabla \cdot \mathbf{B} = 0$
维持 $\nabla \cdot \mathbf{B} = 0$ 是数值MHD的核心问题。违反此条件会产生非物理的磁单极力，导致计算失败。
哈特曼流 — MHD的基本验证问题
MHD最基本的解析解是哈特曼流。在平行平板间的流动中，对垂直于板的方向施加磁场 $B_0$。
当 Ha = 0时为抛物线分布（Poiseuille流），Ha → ∞时除了哈特曼层（厚度 $\delta_H \sim H/\text{Ha}$）外流动分布趋于平坦。
| Ha | 速度分布 | 哈特曼层厚度 |
|---|---|---|
| 0 | 抛物线 | 无 |
| 10 | 略平坦化 | $H/10$ |
| 100 | 中部近似平坦 | $H/100$ |
| 1000 | 仅在哈特曼层有梯度 | $H/1000$ |
Ha = 1000时哈特曼层只有通道宽度的1/1000。网格要求看起来很严格。
没错。为解析哈特曼层，需要沿壁面方向配置与Ha数成正比的网格层。这是MHD计算的主要成本之一。
MHD数值方法的难关——维持"磁场发散为零"的困境
MHD分析中最令人困扰的是 ∇·B=0（磁通连续性）的保持。在有限体积法中离散磁场会积累数值误差，导致 ∇·B≠0，产生非物理的"磁单极"力，解会严重扭曲。为此开发了"发散清理"和"约束传输法"等技术。特别在太阳风大规模MHD模拟中，发散误差超过1%时会产生虚假的"磁单极"，等离子体被错误地偏转。磁场发散为零既是物理约束，也是数值方法的战场。
磁气流体力学（MHD）的实务应用
MHD解析的实务流程
请教我实际进行MHD的CFD分析的步骤。
以连续铸造的电磁制动（EMBr）为例说明。
第1步：磁场计算
- 先用电磁场分析软件（ANSYS Maxwell、COMSOL AC/DC、JMAG等）从磁石和线圈几何计算磁场分布 $\mathbf{B}_0(\mathbf{x})$
- 将结果导入CFD求解器（映射）
第2步：流体域网格划分
- 哈特曼层解析：第一层厚度 ≈ $H/\text{Ha}$ 的1/5以下
- 典型例：$H = 0.1\,\text{m}$、$\text{Ha} = 500$ → $\delta_H = 0.2\,\text{mm}$ → 第一层 = 0.04 mm
- 从壁面用1.1-1.2的增长率配置棱柱层
第3步：物性值设置
| 物性 | 溶钢 (1550°C) | 熔融铝 (700°C) |
|---|---|---|
| 密度 $\rho$ | 7000 kg/m$^3$ | 2385 kg/m$^3$ |
| 粘度 $\mu$ | 0.006 Pa·s | 0.0013 Pa·s |
| 电导率 $\sigma$ | $7.14 \times 10^5$ S/m | $3.5 \times 10^6$ S/m |
| 热导率 $k$ | 40 W/(m·K) | 100 W/(m·K) |
第4步：边界条件
- 入口：喷嘴流入速度（通常 1-3 m/s）
- 出口：压力指定
- 墙面：不滑移 + 电绝缘或导电
- 电边界条件很重要：$\mathbf{J} \cdot \mathbf{n} = 0$（绝缘墙）或 $\phi = 0$（完全导电墙）
墙面的电边界条件是MHD特有的。导电墙和绝缘墙的结果会不同吗？
差异很大。导电墙允许电流通过，制动效应增加。哈特曼层厚度也不同。绝缘墙时 $\delta_H \sim H/\text{Ha}$，导电墙时 $\delta_H \sim H/\sqrt{\text{Ha}}$。
验证的要点
MHD分析的验证检查清单：
- 哈特曼流再现：与解析解比较。速度分布和压力降应一致
- 压力损失验证：MHD附加压力损失：$\Delta p_{\text{MHD}} \sim \sigma U B^2 L$（与Ha²成正比）
- 电流守恒：检查 $\nabla \cdot \mathbf{J} = 0$ 的满足程度
- 磁场无散度：监测 $\nabla \cdot \mathbf{B}$ 最大值
- 实验比较：与UDV（超声多普勒速度计）测得的液态金属速度数据比较
液态金属不透明无法用PIV，所以用超声波。
没错。液态金属速度测量用UDV或CIFT（非接触式感应流量断层扫描）。实验数据有限，所以解析解验证特别重要。
电磁搅拌（EMS）改善连续铸造品质——钢铁业的MHD应用
在钢铁工业的连续铸造工艺中，"电磁搅拌（EMS：Electromagnetic Stirring）"已成为标准质量改善技术。在铸型或二冷区安装的线圈产生旋转磁场，在溶钢中诱导洛伦兹力驱动的搅拌流，实现：①中心偏析抑制、②夹杂物浮起促进、③等轴晶率提高。CFD-MHD耦合分析用于同步模拟溶钢流动、温度场和凝固前沿，优化EMS电流、频率和配置。国内钢铁厂实例显示，CFD优化EMS条件后中心偏析指标改善0.05以上（0→C从0.20→0.15），高强度钢产率提升5%，已在钢铁学会刊物上发表。
磁气流体力学（MHD）的软件比较
MHD支持工具比较
请介绍可用于MHD模拟的软件。
MHD需要电磁场与CFD耦合，软件选型很关键。
ANSYS Fluent（MHD模块）
Fluent集成了MHD模块。
- 支持：低Rm近似（电位法）。读入外部磁场，自动计算洛伦兹力
- 设置：Models > MHD 启用。磁场由用户定义或配置文件输入
- 耦合：Ansys Maxwell或Ansys EMF的磁场数据通过Workbench映射
- 湍流：MHD对湍流的影响建模（Sommeria-Moreau、Knaepen-Moreau）
- 限制：高Rm问题、自诱导磁场不支持
COMSOL Multiphysics
COMSOL在MHD方面很强。
- AC/DC模块 + CFD模块耦合求解MHD
- 低Rm/高Rm均支持：既支持电位法也支持完整诱导方程
- 洛伦兹力自动计算：电磁场与流体场双向耦合是标准功能
- 交流磁场支持：电磁搅拌等交流磁场建模
- 缺点：大规模3D问题（>1000万单元）内存和计算时间困难
OpenFOAM
OpenFOAM包含MHD求解器。
- mhdFoam：完整求解磁场诱导方程的求解器（高Rm支持）
- electrokinetics：电场驱动流（EOF、电化学）
- 自定义实现：许多研究者自开发低Rm近似的电位法求解器
- 社区：大量哈特曼流和MHD湍流教程
专业工具
| 工具 | 用途 | 特点 |
|---|---|---|
| ANSYS Maxwell | 电磁场计算 | 磁场分布计算。用于CFD耦合 |
| JMAG | 电磁场计算 | 日本制造。电动机和变压器设计强项 |
| Plato/Athena | MHD专用（天体物理） | 基于Godunov法的高Rm求解器 |
| FLASH | 天体MHD | AMR支持，核聚变和天体物理 |
| SFEMaNS | MHD专用 | 谱有限元法。精度高 |
耦合工作流
若用不同软件分别计算电磁场和CFD，工作流是怎样的？
典型耦合工作流（连续铸造情况）：
```
Maxwell/JMAG → 磁场分布B(x,y,z) → 导出
↓
Fluent/STAR-CCM+ 将磁场
映射到CFD网格
↓
MHD耦合计算（定常或非定常）
↓
```
磁场映射精度很重要。电磁计算网格与CFD网格通常不同，需注意插值误差。ANSYS Workbench的System Coupling可自动执行网格间映射。
选型指南
| 应用 | 推荐 | 理由 |
|---|---|---|
| 连续铸造EMBr | Fluent MHD + Maxwell | 应用案例丰富 |
| 铝电解槽 | COMSOL | 交流磁场支持 |
| 核聚变毛细血管 | OpenFOAM | 高Ha数定制实现 |
| 天体物理 | Athena++、FLASH | 高Rm、可压缩MHD |
| 研究（通用） | COMSOL | 学习成本低 |
工业用Fluent+Maxwell，学术研究用COMSOL或OpenFOAM。
MHD求解器——开源走在商业前面的罕见领域
MHD分析是少数开源软件领先于商用工具的流体学科之一。OpenFOAM的mhdFoam和GPU支持的Athena++等由核聚变和宇宙物理社区长期主导开发，新算法的实现速度快于商业代码。而ANSYS Fluent的MHD模块则在工业用途（电磁搅拌、铝铸造）积累了验证数据。形成了"研究用开源、工业用ANSYS"的明确分工，这是MHD领域的特色。
磁气流体力学（MHD）的先进研究
MHD湍流
磁场对湍流有什么影响？
MHD湍流与通常湍流在质的方面有所不同。外部磁场引起湍流的各向异性。
形成"雪茄形"涡结构，磁场方向的尺度大，磁场垂直方向的尺度小。这是因为焦耳散逸会衰减垂直于磁场的速度分量。
MHD湍流的能量谱偏离通常的Kolmogorov理论（$E(k) \propto k^{-5/3}$）。
- 弱MHD（$N \ll 1$）：接近Kolmogorov谱
- 中等（$N \sim 1$）：各向异性的级联
- 强MHD（$N \gg 1$）：Iroshnikov-Kraichnan谱 $E(k) \propto k^{-3/2}$（Alfven波相互作用）
强磁场会完全改变能量级联的机制。
为用RANS模型表达这种情况，提出了带MHD修正的湍流模型。Knaepen & Moreau (2008) 的修正 $k$-$\varepsilon$ 模型中增加了焦耳散逸项：
这一项表示湍流能量的电磁散逸。
液态金属电池
近年引人注目的MHD应用是液态金属电池。由熔融金属正负极和熔融盐电解质组成，是大型储能电池，用于可再生能源的电网存储。
MHD的课题：
- 充放电电流的自诱导磁场搅拌液态金属，引起层间混合（短路）
- 泰勒不稳定性：轴向电流超过临界值时产生涡
- 晃动不稳定性：界面波动与电磁力的耦合
CFD用VOF或Level-Set追踪3层（正极/电解质/负极）界面，同时耦合电磁力求解。OpenFOAM派生的专用求解器由多个研究小组开发。
三层液态金属被磁场搅动，设计上要抑制这种扰动。
核聚变炉的MHD
核聚变炉的毛细血管采用PbLi（铅锂合金）或FLiBe（氟化锂铍）等液态金属作冷却剂兼钍增殖剂。
磁场强度达5-10 T，哈特曼数达 $10^4 \sim 10^5$。如此极高的Ha数，常规CFD网格无法解析哈特曼层和侧层。
应对方法：
- 壁函数型方法：将哈特曼层的解析解作为壁面边界条件
- 准2D（Sommeria-Moreau）模型：假设沿磁场方向变化小，化为2D
- 谱法：用高阶多项式精确捕捉哈特曼层的急剧变化
最新研究动向
MHD的最前沿在哪里？
从核聚变到天体物理，MHD的应用范围真的很广。
核聚变炉的MHD——液态金属冷却剂的磁场下流动与"MHD制动效应"
核聚变炉（ITER/DEMO）第一壁冷却采用液态金属（锂或共晶合金PbLi）。但在强磁场（4〜12T）中流动的液态金属受洛伦兹力作用，流动被拉伸至磁场方向，横断面被抑制——出现"哈特曼效应（Hartmann Effect）"。这引发"MHD制动效应"，导热性大幅下降，1MW/m²的设计热流量难以移除。CFD（MHD方程与N-S耦合）用于毛细血管设计，是核聚变工程的前沿课题。ITER机构和德国KIT合作发布了标准计算。Ha>1000的核聚变条件下常规RANS模型失效，需专用MHD求解器。
磁气流体力学（MHD）故障排除
MHD计算的典型故障
MHD计算中最容易踩的坑是什么？
按模式整理MHD特有的故障。
1. 哈特曼层未解析导致压力损失过小
现象：MHD引起的压力损失远低于理论值
原因：哈特曼层（厚度 $\delta_H \sim H/\text{Ha}$）网格解析度不足。壁面附近的速度梯度和电流无法精确计算。
对策：
- 哈特曼层内至少配5个单元
- Ha = 100、$H$ = 0.1 m 的情况：$\delta_H$ = 1 mm、第一层 = 0.2 mm
- Ha = 1000 的情况：$\delta_H$ = 0.1 mm、第一层 = 0.02 mm
- 壁面棱柱层成长率控制在1.1以下
哈特曼数越大网格要求越严格。
2. 电流守恒破裂
现象：$\nabla \cdot \mathbf{J}$ 不为零，非物理的充电现象
原因：电位泊松方程解不充分，或电流密度计算格式不为守恒型
对策：
- 电位方程收敛基准足够严格（残差 $< 10^{-8}$）
- 电流密度面值用守恒格式计算（梯度的面心估值）
- 检查壁面电边界条件（绝缘壁：$\partial\phi/\partial n = (\mathbf{u} \times \mathbf{B}) \cdot \mathbf{n}$）
3. 磁场div B误差（高Rm问题）
现象：非物理的磁单极力出现，等离子体行为异常
原因：诱导方程离散化没有保证 $\nabla \cdot \mathbf{B} = 0$ 的数值满足
对策：
- 约束传输：磁通在面上定义，用精确的离散Stokes定理更新
- 发散清理：引入 $\psi$，$\partial\mathbf{B}/\partial t + \nabla\psi = ...$ 以阻尼磁单极
- Powell 8-波法：在源项中加 $-\nabla \cdot \mathbf{B}$ 比例项稳定
div B = 0 的维持是MHD计算的生命线。
4. 湍流模型不当
现象：本应层流化的区域湍流粘度过大
原因：标准 $k$-$\varepsilon$ 模型未考虑MHD对湍流的抑制
对策：
- 用带MHD修正的湍流模型（Fluent MHD 模块内含）
- 采用LES（MHD效应通过尺度分解自然考虑）
- DNS（可行时）
- Stuart数 $N = \text{Ha}^2/\text{Re}$ > 1 时，忽略湍流先做层流计算
调试步骤
MHD计算的系统调试流程：
1. 无磁场（Ha = 0）收敛流场
2. 弱磁场（Ha ≈ 10）施加，与哈特曼流比较验证实现
3. 逐步升高磁场至目标值
4. 各阶段比较压力损失：$\Delta p_{\text{MHD}}/\Delta p_0 \propto \text{Ha}^2$ （高Ha渐近）
5. 可视化电流密度分布，检查物理合理性
又是分阶段的策略。从零逐步增加哈特曼数。
MHD耦合强，贸然高Ha数开始易发散。哈特曼流是确立的验证问题，用上它是明智之举。
MHD CFD发散——洛伦兹力项的隐式/显式处理选择失误
MHD CFD的常见发散问题：洛伦兹力项J×B用显式（Explicit）加入N-S方程时，电流密度J与速度u的耦合变得不稳定，CFL条件变极小——Ha>100强磁场下必须使CFL小于Ha²/Re数倍，实际计算时间变得不可行。解决方案是采用"隐式洛伦兹力修正"格式，或改用"完整MHD"求解器同时求解诱导方程。低Rm（磁雷诺数<1）的液态金属条件可用"低Rm近似"，忽略诱导磁场大幅简化计算。
帮助
更详细
错误