颗粒温度

通过求解颗粒能量守恒方程，可以估计颗粒温度 $θ_{p}$ 。

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

\frac{3}{2} [\frac{\partial}{\partial t} \int_{V} α_{p} ρ_{p} θ_{p} d V + \oint_{A} α_{p} ρ_{p} θ_{p} v_{p} \cdot d a] = \int_{V} τ_{p} : \nabla v_{p} d V + \oint_{A} κ_{e f f} \nabla θ_{p} \cdot d a - \int_{V} γ d V - \int_{V} J_{p} d V + \int_{V} α_{p} ρ_{p} ϵ_{p} d V

(2375)

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

有效颗粒扩散系数 $κ_{e f f}$ 为：

κ_{e f f} = κ + \frac{3 μ_{p t}}{2 σ_{p t}}

(2376)

其中：

$κ$ 为颗粒扩散系数
$μ_{p t}$ 为颗粒相的湍流粘度
$σ_{p t}$ 为湍流颗粒扩散普朗特数。默认值为 1。

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

Gidaspow 形式：

κ = \frac{150 ρ_{p} d_{p} \sqrt{π θ_{p}}}{384 (1 + e) g_{0}} [\frac{α_{p}}{α_{p, s u m}} + \frac{12}{5} α_{p} (1 + e) g_{0} + \frac{36}{25} α_{p} α_{p, s u m} {(1 + e)}^{2} {g_{0}}^{2}] + 2 α_{p} α_{p, s u m} ρ_{p} d_{p} (1 + e) g_{0} \sqrt{\frac{θ_{p}}{π}}

(2377)

Syamlal 形式：

κ = \frac{15 α_{p} ρ_{p} d_{p} \sqrt{π θ_{p}}}{4 (41 - 33 η)} [1 + \frac{12}{5} η^{2} (4 η - 3) α_{p, s u m} g_{0} + \frac{16}{15 π} (41 - 33 η) η α_{p, s u m} g_{0}]

(2378)

其中：

$ρ_{p}$ 为颗粒密度。
$d_{p}$ 为颗粒直径。相互作用长度尺度设为与颗粒直径相等。
$η = \frac{1}{2} (1 + e)$

颗粒能量的耗散 $γ$ 为：

γ = \frac{12 (1 - e^{2}) g_{0}}{d_{p} \sqrt{π}} α_{p} α_{p, sum} ρ_{p} {θ_{p}}^{3 / 2}

(2379)

由于相间曳力导致的颗粒能量耗散 $J_{p}$ 为：

J_{p} = A_{D} [3 θ_{p} - \frac{A_{D} d_{p} {(v_{g} - v_{p})}^{2}}{4 α_{p} ρ_{p} g_{0} \sqrt{π θ_{p}}}]

(2380)

其中， $A_{D}$ 为相间动量传递系数。

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

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

颗粒温度的代数模型

通过假设局部平衡（其中能量耗散和结果相互平衡），颗粒能量方程 Eqn. (2375) 可以简化为直接针对 $θ_{p}$ 求解的二次方程：

α_{p} K_{4 p} θ_{p} + α_{p} K_{1 p} t r (S) \sqrt{θ_{p}} - {K_{2 p} {[t r (S)]}^{2} + 2 K_{3 p} (S : S)} = 0

(2381)

其中：

K_{1 p} = 2 g_{0} ρ_{p} (1 + e) + 3 (e^{2} - 1) ρ_{p} g_{0}

(2382)

K_{2 p} = \frac{4}{3 \sqrt{π}} α_{p} ρ_{p} d_{p} g_{0} (1 + e) - \frac{2}{3} K_{3 p}

(2383)

K_{3 p} = \frac{4}{5 \sqrt{π}} α_{p} ρ_{p} d_{p} g_{0} (1 + e) + \frac{μ_{p}^{K}}{\sqrt{θ_{p}}}

(2384)

K_{4 p} = \frac{12 (1 - e^{2})}{d_{p} \sqrt{π}} ρ_{p} g_{0}

(2385)

其中 $S$ 为应变率张量， $t r (S)$ 为应变率张量轨迹， $ρ_{p}$ 为颗粒密度， $d_{p}$ 为颗粒直径。相互作用长度尺度与颗粒直径相等。对有效颗粒粘度 $μ_{p}^{K}$ 的动能贡献使用采用 Eqn. (2395) 或 Eqn. (2396) 的某一形式进行定义。

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

K_{3 p} = \frac{4}{5 \sqrt{π}} α_{p} ρ_{p} d_{p} g_{0} (1 + e)

(2386)

默认情况下，域中的最小颗粒温度设为 $1 \times 10^{- 10} m^{2} / s^{2}$ ，但用户可以根据需要更改其值。要估计边界处的颗粒温度，可推算面网格单元处的颗粒温度，或者使用用户自定义的值。指定的颗粒温度值用于计算边界处的固相压力，它对结果具有显著影响。

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

P_{p} = P_{p}^{C} + P_{p}^{K}

(2387)

P_{p}^{C} = 2 g_{0} ρ_{p} α_{p}^{2} θ_{p} (1 + e)

(2388)

P_{p}^{K} = ρ_{p} α_{p} θ_{p}

(2389)

总体粘度为 [546]：

ξ_{p} = \frac{4}{3} {α_{p}^{2} ρ_{p} d_{p} g}_{0} (1 + e) \sqrt{\frac{θ_{p}}{π}}

(2390)

颗粒相的固体粘度由碰撞、动能和摩擦效应组成：

μ_{p, s} = μ_{p}^{C} + μ_{p}^{K} + μ_{p}^{f}

(2391)

层流的有效颗粒粘度定义为：

μ_{p}

μ_{p} = μ_{p, s}

(2392)

对于湍流，计算有效颗粒粘度时还要考虑湍流粘度的影响：

μ_{p} = μ_{p}^{C} + μ_{p}^{K} + μ_{p}^{f} + α_{p} μ_{t}

(2393)

用户自定义的最大固体粘度限制应用于 $μ_{p}$ 。

μ_{p}^{C} = \frac{4}{5} {α_{p}^{2} ρ_{p} d_{p} g}_{0} (1 + e) \sqrt{\frac{θ_{p}}{π}}

(2394)

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

Gidaspow [464] 形式为：

$μ_{p}^{K} = \frac{10 ρ_{p} d_{p} \sqrt{{π θ}_{p}}}{96 (1 + e) g_{0}} {[1 + \frac{4}{5} g_{0} α_{p} (1 + e)]}^{2}$

(2395)

Syamlal [429] 形式为：

$μ_{p}^{K} = \frac{{α_{p} ρ}_{p} d_{p} \sqrt{{π θ}_{p}}}{6 (3 - e)} [0.5 (3 e + 1) + \frac{2}{5} (1 + e) (3 e - 1) g_{0} α_{p}]$

(2396)

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