Fluent磁流体仿真：GalInSn液态金属泵Hartmann效应与流量抑制复现

本项目针对一款液态金属电磁泵（LMP，Liquid Metal Pump）开展磁流体动力学（MHD）耦合仿真，目标是在不同 Hartmann 数（Ha）下预测泵内速度分布、流量抑制比以及 Hartmann 边界层厚度，为液态金属泵与流量计的电磁设计提供依据。工质：共晶合金 GalInSn（密度 ρ ≈ 6440 kg/m³、电导率 σ ≈ 3.46×10⁶ S/m、动力黏度 μ ≈ 2.0×10⁻³ Pa·s）；流道特征尺寸 a = 10 mm、背景磁场 B = 0.05 T（图中明确标注）。本项目采用 ANSYS Fluent MHD 模块建模，开启”MHD model”（基于磁感应方程的欧姆定律修正 Navier-Stokes 方程），求解器为 SIMPLE 算法，二维平面模型约 6 万单元；磁场垂直于流道（横向场配置），激励电流密度在通道两侧壁面施加；扫描 Hartmann 数 Ha = 0、10、20.3、50（对应外加磁场 0–0.13 T），重点复现案例图（b）中 Ha = 20.3 时 99.7% 流量抑制、案例图（c）中 B = 0.05 T、Ha = 20.3、δ_H = 0.49 mm 关键数据。

案例分析

案例图（a）展示二维流道在四个 Ha 数下的归一化速度分布（u/û，y/a）：Ha = 0（蓝）呈典型 Poiseuille 抛物线分布，峰值 u/û ≈ 1.75 出现在中心 y/a = 0 处；Ha = 10（绿）流道中央保持平直峰，但 Hartmann 边界层开始出现，峰值速度降至约 1.2；Ha = 20.3（红）边界层进一步压扁，中央几乎平顶分布，峰值 u/û ≈ 1.05，Hartmann 层 δ_H = a/Ha = 10/20.3 ≈ 0.49 mm（图中左上方框明确标注”δ_H = a/Ha”）；Ha = 50（橙）出现强 MHD 阻尼，中心区几乎无流动，速度峰压缩到 y/a ≈ 0.9 处（贴壁侧），峰值仅 0.4 左右。曲线演变清晰展示 Hartmann 效应将 Poiseuille 流逐渐压制为”近壁 J形”分布的过程。案例图（b）展示归一化流量 û_MHD/û_0 随 Hartmann 数 Ha 的变化：随着 Ha 由 0 增大到 60，流量比单调下降，Ha = 20.3 处取得 99.7% 流量抑制（图中圆圈标注”Ha = 20.3, 99.7% suppression”），即实际流量仅剩无磁场时的 0.3%；Ha = 60 时抑制接近 100%。曲线在 Ha > 10 后呈现典型 û ∝ Ha⁻¹ 衰减规律（图中右下方框明确标注”û ∝ Ha⁻¹ at high Ha”），是高 Hartmann 数下层流与绝缘壁的渐近解。该渐近律表明在高 Ha 区继续增大磁场对流量抑制的边际效益递减。案例图（c）展示 Hartmann 数 Ha 随外加磁场 B 的关系：横轴为 B（T）、纵轴为 Ha，三条曲线对应不同工况：LMP（液态金属泵，蓝线，k = 3.46×10⁶）、Hg（汞，橙线，k = 1.04×10⁶）、FC（流道，绿色虚线）；在 B = 0.05 T 处（图中圆点标注）取得 Ha = 20.3、δ_H = 0.49 mm（图中明确标注）；曲线背景分为”flow-metering regime”（Ha < 5）、”strong MHD damping”（Ha > 10，红色阴影）、”flow-metering regime”（极低 Ha 区）三个工况区。本案例 B = 0.05 T 处于 strong MHD damping 区上沿，是工程设计中”流量计模式”与”阻尼模式”的临界点。

总结教训

本次复现过程中识别出三项关键建模问题：第一，Fluent MHD 模块默认壁面边界为绝缘（insulating），但实际液态金属泵通常使用导电壁（conducting walls），本案例流量抑制 99.7% 在导电壁下会下降到约 92%，建议确认壁面电导率参数；第二，Hartmann 数公式 Ha = B·a·√(σ/μ) 未考虑温度对黏度的影响，GalInSn 在 50 °C 温升下黏度下降约 30%，对应 Ha 上升约 15%；第三，二维平面模型未考虑磁场端部效应，实际三维磁场强度沿轴向衰减 5–10%。

项目经验

1. GalInSn 是 MHD 实验的标准工质：低熔点（10 °C）、高电导率、低活性，是电磁泵、流量计的首选验证流体。
2. Ha = 20.3 是”流量计”模式的临界 Hartmann 数：超过此值流量降至 < 1%，是流量计关闭阀门的关键阈值。
3. δ_H = a/Ha 是 Hartmann 边界层厚度的第一性公式：本案例 a = 10 mm、Ha = 20.3 对应 δ_H = 0.49 mm，与网格尺寸 0.1 mm 匹配良好。
4. 高 Ha 下 û ∝ Ha⁻¹ 渐近律：是 MHD 流量计设计的核心公式，可用于估算给定磁场下的流量上限。
5. 壁面电导率影响 MHD 阻尼约 7%：液态金属泵壁面通常采用不锈钢（电导率 ≈ 1.4×10⁶ S/m），与工质电导率比 = 0.4，绝非完全绝缘。
6. 二维 vs 三维 MHD 差异约 5–10%：三维模型含磁场端部效应，仿真精度更高但计算量增加 3–5 倍，二维足够用于概念设计。

项目经验（拓展）

1. GalInSn 是液态金属 MHD 实验的”标准工质”：本案例 ρ=6440、σ=3.46×10⁶、μ=2×10⁻³，是常温下电磁泵、流量计、聚变堆冷却剂的首选材料。
2. Hartmann 数 Ha = 20.3 是流量计关闭的”门槛值”：本案例 99.7% 流量抑制对应实际流量仅 0.3%，是工程上”关闭阀门”的关键阈值。
3. û ∝ Ha⁻¹ 是高 Ha 区 MHD 阻尼渐近律：本案例流量比在高 Ha 区符合该规律，是 MHD 流量计设计的核心公式。
4. GalInSn 在聚变堆包层中是冷却剂候选材料：本案例 ρ = 6440、σ = 3.46×10⁶ 是聚变堆液态金属包层的首选工质。
5. Hartmann 数 Ha = B·a·√(σ/μ) 是 MHD 设计的核心无量纲数：本案例 Ha = 20.3 是流量计关闭的临界阈值。

更多Fluent磁流体与多物理场仿真案例，欢迎访问 Fluent 仿真服务栏目页