颗粒温度

通过求解颗粒能量守恒方程,可以估计颗粒温度 θp

波动颗粒动能的守恒方程给定为:

1. EQUATION_DISPLAY
32[tVαpρpθpdV+Aαpρpθpvpda]=Vτp:vpdV+Aκeffθpda-VγdVVJpdV+VαpρpϵpdV
(2366)

上面的 Eqn. (2366) 中的第一项是颗粒能量结果。壁面处理用于计算与壁面相邻的第一个网格单元中的颗粒能量结果的已修正值。

有效颗粒扩散系数 κeff 为:

2. EQUATION_DISPLAY
κeff=κ+3μpt2σpt
(2367)

其中:

  • κ 为颗粒扩散系数
  • μpt 为颗粒相的湍流粘度
  • σpt 为湍流颗粒扩散普朗特数。默认值为 1。

颗粒扩散系数 κ 有两种形式:

  • Gidaspow 形式:
3. EQUATION_DISPLAY
κ=150ρpdpπθp384(1+e)g0[αpαp,sum+125αp(1+e)g0+3625αpαp,sum(1+e)2g02]+2αpαp,sumρpdp(1+e)g0θpπ
(2368)
  • Syamlal 形式:
4. EQUATION_DISPLAY
κ=15αpρpdpπθp4(41-33η)[1+125η2(4η-3)αp,sumg0+1615π(41-33η)ηαp,sumg0]
(2369)

其中:

  • ρp 为颗粒密度。
  • dp 为颗粒直径。相互作用长度尺度设为与颗粒直径相等。
  • η=12(1+e)

颗粒能量的耗散 γ 为:

5. EQUATION_DISPLAY
γ=12(1-e2)g0dpπαpαp,sumρpθp3/2
(2370)

由于相间曳力导致的颗粒能量耗散 Jp 为:

6. EQUATION_DISPLAY
Jp=AD[3θpADdp(vgvp)24αpρpg0πθp]
(2371)

其中,AD 为相间动量传递系数。

第一项归功于 Gidaspow 的 [466],将始终包括在内。第二项归功于 Louge 等人的 [510],仅当颗粒能量传递模型的交叉相关项属性设为 Louge 时才包括在内。

Eqn. (2366) 中的最后一项 αpρpϵp 是颗粒相湍动能耗散所产生的源项。仅当使用湍流粘滞态对颗粒相建模时,此项才存在。

颗粒温度的代数模型

通过假设局部平衡(其中能量耗散和结果相互平衡),颗粒能量方程 Eqn. (2366) 可以简化为直接针对 θp 求解的二次方程:

7. EQUATION_DISPLAY
αpK4pθp+αpK1ptr(S)θp{K2p[tr(S)]2+2K3p(S:S)}=0
(2372)

其中:

8. EQUATION_DISPLAY
K1p=2g0ρp(1+e)+3(e21)ρpg0
(2373)
9. EQUATION_DISPLAY
K2p=43παpρpdpg0(1+e)-23K3p
(2374)
10. EQUATION_DISPLAY
K3p=45παpρpdpg0(1+e)+μpKθp
(2375)
11. EQUATION_DISPLAY
K4p=12(1-e2)dpπρpg0
(2376)

其中 S 为应变率张量, tr(S) 为应变率张量轨迹, ρp 为颗粒密度, dp 为颗粒直径。相互作用长度尺度与颗粒直径相等。对有效颗粒粘度 μpK 的动能贡献使用采用 Eqn. (2384)Eqn. (2385) 的某一形式进行定义。

当使用颗粒温度的代数模型时,将忽略 Eqn. (2375) 中的最后一项:

12. EQUATION_DISPLAY
K3p=45παpρpdpg0(1+e)
(2377)

默认情况下,域中的最小颗粒温度设为 1×1010m2/s2 ,但用户可以根据需要更改其值。要估计边界处的颗粒温度,可推算面网格单元处的颗粒温度,或者使用用户自定义的值。指定的颗粒温度值用于计算边界处的固相压力,它对结果具有显著影响。

固相压力、总体粘度和有效颗粒粘度的表达式作为颗粒温度的函数获得, θp 。固相压力包含碰撞和运动效应 [564]

13. EQUATION_DISPLAY
Pp=PpC+PpK
(2378)
14. EQUATION_DISPLAY
PpC=2g0ρpαp2θp(1+e)
(2379)
15. EQUATION_DISPLAY
PpK=ρpαpθp
(2380)

总体粘度为 [542]

16. EQUATION_DISPLAY
ξp=43αp2ρpdpg0(1+e)θpπ
(2381)

有效颗粒粘度 μp 包含碰撞、运动和摩擦效应:

17. EQUATION_DISPLAY
μp=μpC+μpK+μpf
(2382)

用户自定义的最大固体粘度限制应用于 μp

18. EQUATION_DISPLAY
μpC=45αp2ρpdpg0(1+e)θpπ
(2383)

可以在两种形式的颗粒动力学粘度之间选择:

  • Gidaspow [466] 形式为:
    19. EQUATION_DISPLAY
    μpK=10ρpdpπθp96(1+e)g0[1+45g0αp(1+e)]2
    (2384)
  • Syamlal [433] 形式为:
    20. EQUATION_DISPLAY
    μpK=αpρpdpπθp6(3-e)[0.5(3e+1)+25(1+e)(3e-1)g0αp]
    (2385)

Syamlal 版本预测颗粒体积分数低于 0.3 时的较低粘度 [564]。建议根据所选的颗粒动力学粘度版本使用相应的线性拖曳法(GidaspowSyamlal)。