融化与凝固建模

融化是将物质状态从固体变为液体的过程。 将固体物质加热到融化温度时,会发生融化。 相反的过程称为凝固。

对于纯物质(例如水),融化与凝固通常发生在同一温度,即液相线和固相线位于相同温度。 融化与凝固瞬时发生,因此根据温度,在计算域中有液体或固体存在。

对于非纯物质(例如合金),融化与凝固会在一个温度范围内发生。 在这种情况下,液相线和固相线不在同一个温度。 在它们之间的温度间隔中,液体有部分凝固;固相和液相同时存在(如浆料)。 此温度间隔具有有限宽度,最多为几个 100 开氏度。 融化与凝固过程不是瞬时的,而是在一个时间段内发生。 在此期间,相混合物会影响流场。

凝固导致的流体阻力

融化与凝固模型仅模拟相变对能量方程的影响;它不会修改动量方程。 需要其他模型才能考虑从液态到固态的材料属性更改。

根据相对固体分数的值(指定液体的凝固程度),有三个可选模型可用:

  • 浆料粘度模型

    对于低固体分数值,混合物被视为浆料,固体颗粒几乎悬浮在熔体中,彼此之间没有相互作用。 它们对流体的影响通常建模为与相对固体分数相关的粘度。

    实现 Metzner 方法后可将粘度(浆料粘度)最高增加到临界固体分数(约为 0.3)(请参见 Eqn. (2713)),这样无需太多互动便可对随流体流动的凝固颗粒建模。 但是,使用大于两个数量级的粘度范围可能造成不稳定。

  • 两相区渗透率模型

    当固相体积分数增大时,等轴颗粒增长,开始凝聚并形成枝晶区域。 枝晶还从冷却固体表面或凝固前沿生长到熔体中。 枝晶渗透的部分凝固固定区域称为“两相区”。 与各向同性多孔介质类似,可以使用 Carman-Kozeny 两相区渗透率方法(请参见 Eqn. (2717))对两相区中的流体阻力建模。 此模型用于高于指定临界固体分数(约为 0.3)的相对固相体积分数。

  • 流体停止功能

    要完全停止正在凝固的流体,可以使用流体停止功能。 当高于指定的流体停止固体分数(约为 0.4 到 1)时,激活融化-凝固流体停止模型。 浆料粘度和两相区渗透率模型添加到移动流体的阻力取决于流体速度。 当速度为零时,这两种模型的影响也为零。 当存在压力梯度或体积力(如重力)时,蠕流的发展独立于浆料粘度或两相区渗透率的幅值。

流体阻力模型的适用性

浆料粘度模型和两相区渗透率模型仅影响部分凝固液体的流体属性。 仅当模拟非纯物质(如合金)时,才需要这些模型。 例如,当模拟冰在固定容器中的融化时不需要它们。

如果模拟纯物质的融化与凝固过程,则仅使用融化与凝固模型。 纯物质不会发生部分凝固。 要防止蠕流,还可以激活流体停止功能。

当模拟非纯物质的融化与凝固时,使用下表来确定要应用的流体阻力模型。

流体速度 固相线/液相线温度间隔 模拟时间 流体阻力模型
低 (< 1 m/s) 短 (< 1 s) 浆料粘度
低 (< 1 m/s) 长 (> 3 s) 浆料粘度

流体停止

低 (< 1 m/s) 短 (< 1 s) 浆料粘度

两相区渗透率

低 (< 1 m/s) 长 (> 3 s) 浆料粘度

两相区渗透率

流体停止

高 (> 5 m/s) 短 (< 1 s) 浆料粘度

(流体停止)

高 (> 5 m/s) 长 (> 3 s) 浆料粘度

流体停止

高 (> 5 m/s) 短 (< 1 s) 浆料粘度

两相区渗透率

(流体停止)

高 (> 5 m/s) 长 (> 3 s) 浆料粘度

两相区渗透率

流体停止