磁流体力学 (MHD)

通过 Simcenter STAR-CCM+,可对导电流体(如熔融金属、电解质和等离子体)与电磁场之间的相互作用进行建模。

导电流体相对于磁场的运动会产生电流密度:

1. EQUATION_DISPLAY
JL=σv×B
(4369)

其中, σ 为导电率(请参见 Eqn. (4228)), v 为流体速度, B 为磁通量密度。Eqn. (4369) 影响总电流密度 J (请参见 Eqn. (4242))。

JL 进而产生磁通量密度,该密度会影响 Eqn. (4369) 中的总磁通量密度。Simcenter STAR-CCM+ 使用单向或双向耦合法计算 JL 。在单向耦合法中,Simcenter STAR-CCM+ 相对于指定的磁通量密度 Bp 来计算 JL 。在双向耦合法中,Simcenter STAR-CCM+ 计算总磁通量密度 B ,该总密度也将由 JL 产生的磁通量密度纳入其中。

单向耦合磁流体动力

产生的磁通量密度和指定的磁通量密度之间的比率由无量纲数(称为磁雷诺数)来定义:

2. EQUATION_DISPLAY
Remag=μ0σUL
(4370)

其中, μ 0 为真空渗透率, U 为特征流体速度, L 为特征长度。

单向耦合磁流体动力方法适用于小的 Remag (其中产生的磁通量密度明显低于指定的磁通量密度)。

在这种情况下,Simcenter STAR-CCM+ 以解耦的方式求解 Eqn. (4242),忽略随时间变化的磁场所产生的涡流(请参见Low-Frequency Electromagnetics in Conducting Media)。在没有涡流的情况下,将 JL (Eqn. (4369)) 添加到总电流密度 J (Eqn. (4234)) 会得出:

3. EQUATION_DISPLAY
J=σϕ+σv×B
(4371)
其中, B 为指定的磁通量密度:
4. EQUATION_DISPLAY
B=Bp
(4372)
由于 J=σE(请参见 Eqn. (4228)),因此电场定义如下:
5. EQUATION_DISPLAY
E=ϕ+v×B
(4373)
双向耦合磁流体动力
使用此方法,Simcenter STAR-CCM+ 将对 Eqn. (4241) 求解(对电势建模时,还将对 Eqn. (4242) 求解)。在这种情况下,电场计算如下:
6. EQUATION_DISPLAY
E=ϕAt+v×B
(4374)
总电流密度为:
7. EQUATION_DISPLAY
J=σE=σ(ϕAt+v×B)
(4375)

其中:

8. EQUATION_DISPLAY
B=∇×A
(4376)
在准非稳态分析中,会忽略 Eqn. (4374)Eqn. (4375) 中的 A/t 项。当 Simcenter STAR-CCM+ 不求解电势时,也会忽略 Eqn. (4374)Eqn. (4375) 中的 ϕ 项。

为了简单起见,Eqn. (4371)Eqn. (4375) 假设 JL 是对 Eqn. (4234) 中的 Jex 的唯一贡献。相关章节中将讨论对 Jex 的额外贡献。

在断路器模拟(等离子体)中,来自电势 Jϕ=σϕ 梯度的对电流密度的贡献通常优于其他贡献。在电磁搅拌应用中,主要贡献通常来自涡流 JA=σA/t 。当磁通量密度和流体速度很重要时, JL 变为主项。这通常是湍流控制应用中的情况,即使用 Lorentz 力来调节或抑制流体湍流。

Lorentz 力

导电流体所受的每单位体积的体积力(称为 Lorentz 力)为:

9. EQUATION_DISPLAY
fL=J×B
(4377)

其中, J Eqn. (4371)Eqn. (4375) 中定义的总电流密度。在磁流体动力应用中,Simcenter STAR-CCM+ 包括动量方程 (Eqn. (655)) 中的 Lorentz 力。

对于正弦场(请参见谐波时间依赖性),Simcenter STAR-CCM+ 将计算循环平均 Lorentz 力,如下所示:

10. EQUATION_DISPLAY
fL=12Re(J^×B^*)
(4378)

其中,在 Eqn. (4269) 中定义 J^B^*B^ 的共轭复数(请参见 Eqn. (4268))。在谐波平衡应用中,Eqn. (4369) 不会影响总电流密度。