重力

重力产生的力通常在流体和固体的动量平衡中起着重要作用。 重力产生的体积力作为源项添加到动量方程中。

1. EQUATION_DISPLAY
fg=ρg
(860)

其中,g 为重力矢量。

当流体密度 ρ 随温度变化时,重力矢量作用于密度,导致自然对流。

当格拉斯霍夫 (Gr) 数和雷诺 (Re) 数的比率接近或超过 1 时,非均匀密度场导致的浮力视为有效:
2. EQUATION_DISPLAY
GrRe2=ρβΔT|g|Lρv2
(861)

其中,β 为热膨胀系数。

对于自然对流,当温度差很小时,可以使用 Boussinesq 模型来模拟浮力效应,而不是使用与温度相关的密度。 Boussinesq 模型中使用的密度在所有求解方程中均为常数,只有浮力源项发生变化。

工作压力

当重力模型处于活动状态时,Simcenter STAR-CCM+ 工作压力会变为测压压力,这时会考虑参考高度(由用户设置)。 测压压力如下:

3. EQUATION_DISPLAY
ppiezo=pstatic-ρrefg(x-x0)
(862)

对于可变密度流和恒密度流,都需要设置参考高度 x0,它用于计算测压压力,如上所示。 参考高度指测压压力与静压被视为相等的点。

对于重力模型处于活动状态的模拟,设置初始或边界压力时实际设置的是测压压力值。 Simcenter STAR-CCM+ 会计算对应的静压。

可变密度流

对于可变密度流,浮力源项为:

4. EQUATION_DISPLAY
fg=(ρ-ρref)g
(863)

其中, ρ ref 为参考密度。

对于 VOF 多相流体,密度变化可能非常显著(好几个数量级),导致数值不稳定。 要在这种情况下获得可靠且稳定的数值解,需要将参考密度设为零或最轻相的密度。

Boussinesq 模型

可以通过使用 Boussinesq 模型增强自然对流模拟的收敛。 如果通过 Boussinesq 近似 ρ=ρref(1βΔT) 从浮力项中去除了 ρ,控制方程的所有项中都使用恒密度 ρref,则浮力源项近似为:
5. EQUATION_DISPLAY
fg=ρrefgβ(Tref-T)
(864)

其中, T ref 为工作温度。

此近似仅对小的温度差(因此仅对小的密度变化)有效。 如果 β(TT0)1,则 Boussinesq 近似有效。