Fluent的融化凝固模型介绍

Fluent的融化凝固模型（来源于帮助文档）

一、 简介

ANSYS Fluent可用于解决涉及在一个温度（例如，纯金属）或一个温度范围（例如，二元合金）下发生的凝固和/或熔化的流体流动问题。ANSYS Fluent使用焓-孔隙度公式，而不是显式跟踪液固界面。将液-固糊状区视为多孔区，其孔隙率等于液体分数，并在动量方程中添加适当的动量汇项，以解释固体材料存在引起的压降。还将会添加到湍流方程中，以说明固体区域中孔隙度的降低。
ANSYS Fluent为凝固和熔化建模提供了以下功能：

• 纯金属和二元合金液固凝固/熔化的计算
• 连铸过程建模（即将固体材料“拉出”领域）
• 固化材料和墙壁之间接触热阻的建模（例如，由于存在气隙）
• 凝固/熔融过程中相变化的模拟
• 凝固/熔化相关量的后处理（即液体分数和拉动速度）
  这些建模功能使ANSYS Fluent能够模拟各种凝固/熔化问题，包括熔化、冻结、晶体生长和连铸。《用户指南建模凝固和熔化》中描述了用于解决凝固/熔化问题计算的物理方程以及设置方法。
  ANSYS Fluent中的公式可用于模拟纯材料以及合金的凝固/熔化。对于合金，ANSYS Fluent提供了两种确定液体分数与温度关系的规则，即线性杠杆规则和非线性Scheil规则。

二、 限制

以下限制适用于ANSYS Fluent中的凝固/熔化模型：

• 凝固/熔化模型只能与基于压力的解算器一起使用；它不适用于基于密度的解算器。
• 凝固/熔化模型不能用于可压缩流。
• 在一般的多相模型（VOF、混合物和欧拉模型）中，只有VOF模型可以与凝固/熔化模型一起使用。
• 除物质扩散率外，不能通过用户界面为固体和液体材质指定单独的材质属性。但是，如果需要，可以使用DEFINE PROPERTY宏使用用户定义的函数指定单独的固体和液体特性。有关使用DEFINE PROPERTY宏的详细信息，请参阅Fluent自定义手册中的DEFINE\U PROPERTY UDF。
• 当将凝固/熔融模型与反应的物质迁移模型结合使用时，没有任何机制将反应仅限于液体区域；也就是说，反应在任何地方都得到了解决。

三、 介绍

ANSYS Fluent中使用焓-孔隙度技术[612]（第979页）、[614]（第979页）、[615]（第979页）对凝固/熔化过程进行建模。在这种技术中，熔体界面没有被明确跟踪。相反，一个称为液体分数的量（表示液态体积的分数）与域中的每个单元相关联。根据焓平衡，在每次迭代时计算液体分数。糊状区是液体分数介于0和1之间的区域。糊状区被模拟为一种“伪”多孔介质，在这种介质中，随着材料的固化，孔隙度从1降低到0。当材料在电池中完全固化时，孔隙率变为零，因此速度也降至零。本节概述了凝固/熔化理论。有关焓孔隙度法的详细信息，请参阅Voller和Prakash【615】（第979页）。

四、 能量方程

温度的解本质上是能量方程（方程20.5（第741页））和液体分数方程（方程20.3（第741页））之间的迭代。直接使用方程20.3（第741页）[616]（第979页）来更新液体分数通常会导致能量方程的收敛性差。在ANSYS Fluent中，使用Voller和Swaminathan建议的方法更新液体分数。对于和相等的纯金属，使用Voller和Prakash[615]（第979页）给出的基于比热的方法。

五、 动量方程

焓孔隙度技术将糊状区域（部分凝固区域）视为多孔介质。每个单元中的孔隙度设置为该单元中的液体分数。在完全凝固区域，孔隙率等于零，从而使这些区域的速度消失。由于糊状区域中孔隙度降低，动量汇呈现以下形式：

式中，是液体体积分数，是一个小数字（0.001），以防止被零除，是糊状区常数，是固化材料拉出畴时的固体速度（也称为拉动速度）。
糊状区常数测量阻尼的振幅；该值越高，材料固化时速度向零的过渡越陡。非常大的值可能会导致解振荡。
当凝固材料在连铸过程中不断从该区域中抽出时，拉速被包括在内，以说明凝固材料的移动。方程式20.6（第741页）中存在该项，允许新固化材料以拉伸速度移动。如果固化材料未从域中拉出。连铸拉速（第744页）中提供了有关拉速的更多详细信息。

六、 湍流方程

在糊状和凝固区的所有湍流方程中加入汇，以说明固体物质的存在。下沉项与动量下沉项非常相似（方程式20.6（第741页））