AMUSIG 破碎和聚结
对于层流,聚结和破碎由连续相的剪切速率驱动。对于湍流,聚结和破碎由连续相的湍流尺度驱动。破碎率和聚结率对为连续相和颗粒引发的湍流选择的湍流模型具有一定敏感度。要使用不同的湍流模型获得可比较的结果,可能需要重新校准模型。
Simcenter STAR-CCM+ AMUSIG 破碎和聚结模型不适用于颗粒流体。但是,可以使用 AMUSIG 对颗粒流体中的固体颗粒相尺寸分布建模:对于相关应用,运用对聚结和破碎现象的理解和知识,根据需要校准幂次定律聚结和破碎模型。
有一组模型可用于描述不同应用的破碎和聚结过程:
- Luo ([513]) 和 O’Rourke ([524]) 模型适用于湍流聚结。
- Tsouris 和 Tavlarides ([515]) 模型、Martinez-Bazan ([515]) 模型、Kocamustafaogullari ([515]) 模型和 Coulaloglou 和 Eskin ([456]) 模型适用于湍流液体/液体和气泡流破碎。
- 剪切感应破碎和聚结模型适用于层流。
- 此外,对于测试、演示并提供一些用户指定的建模,也可以实施常数和幂次定律模型。
聚结率模型包含碰撞率和聚结效率子模型。破碎率模型包含破碎率和碎片颗粒尺寸分布子模型,可定义给定尺寸颗粒的移除项。所有比率模型都有一个整体标定常数,用户可以根据自身需求调整该常数。临界韦伯数为可调参数,与破碎率模型相关。
单个尺寸组足以对深舱或深塔中上升的稀释泡沫的静水膨胀建模。但是,要使碰撞/破碎研究有意义,建议至少使用三个组。AMUSIG 模型具有自适应性,可以使用少量组研究颗粒尺寸的显著和快速变化。增加组数可能会产生与分辨率无关的结果,作用类似于网格加密。计算成本确定了尺寸组数的实际上限。
层流和湍流破碎模型
破碎过程建模分为两个部分:破碎率和碎片颗粒尺寸分布。破碎率具有时间-1 维度,因此单个颗粒在时间跨度
- 幂次定律破碎率
-
此模型适用于湍流破碎和层流破碎。幂次定律破碎率模型描述给定尺寸颗粒的破碎速度。该模型是一个通用模型,针对特定颗粒尺寸
(按特征直径 进行比例缩放)下数密度的破碎率乘数 提供了可调参数。(2244)由于
和数密度的乘积是一个比率(以事件/秒/ 为单位),因此标定常数 的维度是秒。此破碎率如同任何其他破碎率一样输入守恒方程:
(2245)
其中:
-
为连续相的湍流耗散率。 -
为颗粒直径。 -
为韦伯数。 -
为临界韦伯数。它确定可以承受给定湍流强度的颗粒尺寸。较大的 值表示破碎程度较弱。 -
为可用于调整破碎时间尺度的标定常数。
对于湍流,韦伯数由以下公式给出:
校正因子
- Martinez-Bazan
- 对于 Eqn. (2246) 中的函数
,Martinez-Bazan 模型假设平方根形状。Martinez-Bazan 模型表示低于最小直径时不发生破碎。 - Tsouris 和 Tavlarides
-
对于 Eqn. (2246) 中的函数
,Tsouris-Tavlarides 模型使用指数。Tsouris-Tavlarides 模型预测任何液滴都有可能破碎(无最小直径),但破碎概率会随液滴直径呈指数下降。 - Coulagloglou-Eskin
- Coulaloglou-Eskin 模型预测的尺寸分布范围大于另外两个模型。
- Kocamustafaogullari
- Kocamustafaogullari 模型适用于对连续气体中的液滴破碎进行建模。该模型考虑液滴内表面张力和粘性力抵抗的变形。此行为的特征度量为韦伯数和奥内佐格数。
对于层流,液滴破碎由毛细管数控制。毛细管数是尝试使液滴变形和破碎的粘性剪切应力与保持液滴球形的毛细管压力之比,由以下公式给出:
其中:
-
为连续相的动力粘度。 -
为平均剪切速率。 -
为表面张力。
层流破碎模型假设存在临界值
其中:
-
为低粘度比率指数前因子。 -
为高粘度比率指数前因子。 -
为低粘度比率指数。 -
为高粘度比率指数。 -
为最大粘度比率。
对于层流破碎模型,破碎率具有以下结构:
其中:
层流破碎模型之间的区别在于 Eqn. (2250) 中的函数
- 克里斯蒂尼 (Cristini)
- 克里斯蒂尼模型将
定义为 ([445]): - Tsouris 和 Tavlarides
- 对于 Tsouris 和 Tavlarides 剪切模型,破碎概率随液滴直径呈指数级下降。函数
由以下公式给出: - Coulaloglou 和 Eskin
- Coulaloglou 和 Eskin 剪切模型将
定义为:
- Kocamustafaogullari
- Kocamustafaogullari 模型基于随机二次液滴 (SSD) 模型,其中假设
以下液滴瞬时破碎率:
(2255)其中:
-
为临界韦伯数。默认值为 12。 -
is a parameter regulating the speed of breakup. Large values signify that breakup occurs later. The default value is . -
为液体的密度。 -
为液滴直径。 -
为气体的密度。 -
为液滴通过连续相的湍流滑移速度。
Kocamustafaogullari ([488]) 通过将 Levich [501] 公式代入湍流液滴的运动方程,关联环形流中液滴尺寸分布的实验数据:(2256)其中, 和 分别为液体和气体的波动速度。假设曳力系数为 ,Eqn. (2256) 随后重新排列如下:(2257)替换 Eqn. (2256) 中的 Eqn. (2258) 得到:
(2260)Eqn. (2260) 在以下条件下达到最大值:(2261)无量纲滑移速度为
,其中, and 为无量纲平均和波动滑移速度。归一化分布由以下公式给出:(2262)(2263)(2264)(2265) -
聚结模型
在碰撞率和聚结效率这两个零部件中对聚结过程建模。
- 湍流碰撞率
-
仅适用于湍流聚结。聚结率是碰撞率和聚结效率的乘积。湍流碰撞率计算为两个球体(直径分别为
和 ,并以随机相对速度 移动)的碰撞概率。该模型没有拟合参数,但聚结效率子模型除外。
- Luo 聚结效率
-
Luo 聚结效率模型 ([513]) 假设当接触时间
大于 时发生聚结,需要该时间才能使分离两个气泡的液膜破碎:(2269)其中,
是该模型的唯一标定常数。 较大表示无聚结;默认值为 。 - O'Rourke 聚结效率
O’Rourke 聚结效率模型 ([524]) 用于检测连续气相中的液滴碰撞。碰撞的液滴会因连续相的湍流波动发生混乱运动。碰撞结果由碰撞韦伯数控制:
(2270)其中:-
为液滴半径。 -
为液滴 1 和 2 的平均表面张力。 -
为平均密度。 -
为两个液滴的相对速度,仅受在 Eqn. (2261) 中计算的湍流波动影响。
-
- 液膜导流聚结效率
-
液膜导流聚结效率模型是与任何碰撞率模型均兼容的通用模型。
(2279)其中,
为碰撞韦伯数:(2280) 为碰撞雷诺数:(2281)碰撞速度
由碰撞率模型提供。 为碰撞液滴或气泡的有效直径。
- 层流碰撞率
- 对于层流碰撞率,两个球体的相对速度是由于平均剪切速率所致: (2282)
- 增强的碰撞率幂次定律
-
此碰撞率模型适用于湍流聚结和层流聚结。
增强的碰撞率幂次定律及其可调参数可用于研究两个尺寸分别为
和 (按特征直径 进行比例缩放)的颗粒之间的替代凝结内核 。(2283)由于内核与两个碰撞尺寸的数密度之积是一个比率(以事件/秒/
为单位),因此标定常数 的维度是 /秒。