第２２卷第６期

能源环境保护

・１５・

知％‘＞毒¨锦）

＝一以ｐ绋一一月ｔ酝

（３）液滴控制方程表达式：

液雾颗粒相连续方程：

昙阱＋兰◇。％）＝一氆镌

瓦阱＋丽妒々％户一氆镌

液雾相数密度连续性方程：

警＋昙＠％）＝ｏ

液雾颗粒褶动量方程：ｗ

（８）

（９）

（１０）

昙饥≯毒‰）＝尝ｔ瑚

＋ｐ湛ｉ—ｖ—Ｉｍｉ

２臻憋白五一％乏嚣，憋

（１２）

将以上方程组进行＆舭展开，本模拟剃

用了Ｂｏｕｓ８妇ｓｑ假说，将雷诺应力项与时均速度

梯度联系起来，使得守恒方程数与变量数相等。

３

输运方程的初始边界条件

本模拟中涉及到三个边界条件：人翻初始条

件、蹙面边界条件和出口条件。

（１）瞻已知模拟条件及经验公式给定入曩边

界上的Ⅱ、舻、ｌ‘，、灿五、ｅ分布，其中ｋｑ，如＝ｏ．００３ｕ２ｑ

￡蚋：船

（１３）

￡一向２瓦丽面尹

（１３）

（２）固体壁面边界条件：采用无滑移条件，即

珏蒯＝影ｗ＊＝球Ｍｌ森＿盛＝￡柚

由于壁面附近速度迅速下降并趋于零。粘性

力将在瀛动中怒主导侔蔫，逶用予离＆数流动懿

一方程不适用予壁面附近的流动，本模拟采用壁

面函数法进行处理：

（３）轴对称线（ｒ＝Ｏ）处边界条件：

芒：要：嬖：要：粤：ｏ

（。４）

ｉ２ｉ２＿２●２÷２ｕ

（１４）

仪

承

暾

礅

＆

（４）出口边界条件：在反应器的出口处烟气分

布不均，放出瞄处采掰篷力边界条件。

４数值求解方法

采用ｓ啪算法，对流相采用二阶迎风差

分格式。模拟计算中所有变量收敛条件定为：￡≤

ｌ×１０一气

符号说明：１

ｅｌ，ｅ厂一湍流模登参数；

ｃｏ，Ｃ扩颗粒与流体间拽力系数；

Ｃｒ气体定篮吃热，Ｊ傀．Ｋ；

Ｃ，厂一液漓颗粒比热，Ｊ傀．Ｋ；

Ｃ广脱硫剂颗粒比热，Ｊ傀．Ｋ；

卜扩散系数；

ＤⅨ——液雾颗粒直径，ｍ；
舀。——脱硫颗粒壹径，ｍ；

肘广一液雾相单个颗粒质鬣，ｋｇ
ｍ，——单个液雾颗粒质量变化率，ｌ【ｇ／ｓ；

ｒ壁蕊法线方向；

ｎｒ液雾颗粒数密度；

珏∥一脱硫剡颗粒数密度；

Ｐ——静压，Ｐａ；

口广一滚匿相阉换热，ｊ／ｓ；

Ｑ，一反应吸热，Ｊ傀；

豫——Ｒ础准数；

ｒ——径向坐标，ｍ；

ｒ——温度，Ｋ；

于广气化湿度，Ｋ；

ｒｂ——大气温度，Ｋ；

乳广一液雾齐温速率，Ｋ／ｓ；

彬——切向速度，ｍ／ｓ；

埘，——化学反应速率；

茗——轴肉坐标，ｍ；

ｌ，∥一气相组分ｓ的质量浓度分数；

卜通用扩散系数；

￡——气相动能耗散率。

５

结果与分析

（下转第３８页）

万方数据