α D (C sat − C H 2 O )Se ∂u g 1 ∂ ∂vg ∂P ∂ + −α g + µe r ∂r rµ e ∂z + α g ρ g g z − ∂z ∂z ∂ z ∂u g ∂vg v g ∂ 2 2 + + α g ρ gk + µe ∂z ∂r r 3 ∂z 3 ∂u g ∂ ∂v g 2µ ev g ∂P ∂ + + µ e r∂r rµ e ∂r − r 2 + ∂r ∂z ∂ r ∂u g ∂v g vg ∂ 2 2 α g ρ g g r − α g ρ g k + µe + + 3 r ∂r ∂r ∂z 3 −α g Gk − α g ρ g ε ε Cε1Gk − Cε 2α g ρ gε k
2
数理模型
基金项目 国家自然科学基金重点项目 50136020 教育部博 士点基金项目 2000000304 Project Supported by National Natural Science Foundation of China (50136020).
2.1 数学模型 循环流化床烟气脱硫设备装置及过程描述参 见文 [2] [3] 将循环流化床烟气脱硫设备反应
第 11 期
颜
岩等
循环流化床内烟气脱硫模拟分析
175
2.2 相关参数 颗粒动力学参数可参见文[6]~[8] 1 颗粒湍动特性分析 存在着颗粒团 还存在颗粒团 循环流化床烟气脱硫设备中 的生成和解体 除单颗粒脉动外 脉动 对颗粒湍动特性描述
液相中 SO 2 传质为 NSO 2 = φα D ,SO 2 ,l CSO 2 ,i 系, 即 PSO 2 ,i = HCSO 2 , i (10) 联立式 (8)~(10) 得到 NSO 2 = φαSO 2 , l CSO 2 / φα SO 2 ,l / α SO 2 , g + H / RT g (10) 总的 SO2 传质系数为 K SO 2 = φα SO 2 ,l /(φ α SO 2 ,l α SO 2 , g + H / RT g )
dp 颗粒直径

β (vg − v p )
湍能 k p 方程 湍流耗散率 εp 方程 能量方程 拟温度方程 SO2 浓度方程 H2 O 浓度方程 kg(m2 ⋅s)−1 P
kp εp cpp Tp
µ s + (µ s , t / ξ ε ) µ s + (µ s , t / ξ ε ) λp µ s, t 2 Γ Θ+ 3 ξΘ D SO 2 , e D H 2O ,e
2 −1
m D 扩散系数 m⋅s−1
Gs 固体质量通量
径向位置; m
Sc Schmidt 准则数 T 温度 K u 轴向速度 m⋅s−1 ; v 运动粘度, m2 ⋅s−1 ; w 床壁面
径向速度, m⋅s
τ 剪切应力 N⋅(m )
ρ 密度 kg⋅(m )
µ 粘度 Pa⋅s φ 增强因子 v
万方数据
Tab. 1
方程 连续方程 ϕ 1 Γϕ 0
间相互作用 下标 i 取 g 和 p 分别代表气相和颗 粒相 控制方程的具体表达式见表 1
φ=0.3m h=4.5m
图 1 循环流化床烟气脱硫设备反应器示意图 Fig. 1 Diagrammatic sketch of CFB-FGD reactor
表 1 控制方程中各项具体表达式 The expression for the governing equations
Sϕ Sϕgp 0
动量方程(z)
ug
µe
− β (u g − u p )
气 相 动量方程(r) vg µe
− β (v g − v p )
湍能 k 方程 湍流耗散率 ε方程 能量方程 连续方程
k ε
µ g + µ g ,t / ξ k µ g + µ g ,t / ξε λe
0
Gp
(
)
(ε / k )Cε 3G p −α H (Tg − T p )Se
0
c pg T g
1
µφ g + rα D (C sat − C H 2 O )Se −α D (Csat − C H 2 O )Se ∂u p 1 ∂ ∂v p ∂P ∂ ∂Ps + + µ p r ∂r rµ p ∂z + α p ρ p g z − ∂z − ∂z ∂z ∂ z ∂u p ∂v p v p ∂ 2 2 + + α p ρ p k p + µ p − ζ s r ∂z ∂r 3 ∂z 3
第 23 卷 第 11 期 2003 年 11 月 文章编号 0258-8013 (2003) 11-0173-05
中 国 电 机 工 程 学 报 Proceedings of the CSEE 中图分类号 TK221 文献标识码 A
Vol.23 No.11 Nov. 2003 ©2003 Chin.Soc.for Elec.Eng. 学科分类号 470⋅2040
−α p
动量方程( z ) up µp
β (ug − u p )
颗 粒 相
动量方程( r )
vp
µp
∂P ∂ ∂uv ∂ ∂v 2µ v + µ v + rµ v v − v v + ∂r ∂z ∂r r∂r ∂r r2 ∂P ∂ 2 ∂v v 2 ∂u α v ρ v g r − 2 − α v ρ v k v + µv − ς 2 v + v + v ∂r ∂r 3 ∂r r 3 ∂z − αv G kp − α p ρ p ε p ε p / k p Cε 1 Gkp − C ε 2 α p ρ p ε p µφ g − rα D C sat − C H 2 O Se G pp − (2 / 3)γ / 3 + (2 / 3)α p ρ pε p − K SO2 C SO2 Se α D (Csat − C H 2 O ) Se
汽化潜热 kJ(kg) −1 ; Sh
3 −1
G gp
(
)
(ε p / k p )Cε 3G gp −α H (T g − T p ) Se
0 0 0
(
)
Θ C SO2
C H 2O
注 aD 对流传质系数 m⋅s−1 ; aH 压力 Pa; r
−1
对流换热系数, kJ(m2 ⋅K)−1 cp 定压比热 kJ⋅K−1 Sherwood 准则数
ABSTRACT A model was developed from two-fluid model including the effect of evaporation and desulfurization to simulate the flue gas desulfurization process in a circulating fluidized bed. The influence of collision between cluster ans particles and wall was considered by introducing the coefficient of restitution, and the particle turbulent energy and dissipation rate of turbulent energy were introduced to descript the cluster turbulence. These considerations greatly improve the consistence of theoretical predictions with experiments and better understand some phenomena such as local solids fraction distribution. The numerical simulation indicated that the cluster could decrease the slurry evaporation rate and enhance the desulfurization. KEY WORDS: Two-fluid model; Circulating fluidization bed; Flue gas desulfurization; Cluster 摘要: 基于气固两相流双流体模型及床内浆滴蒸发和 SO2 吸收过程的分析 团 特性 过程 该文发展了循环流化床烟气脱硫过程 采用恢复系数修正颗粒与颗粒 合理解释固相含率沿径向分布 的双流体 / 蒸发脱硫模型 壁面之间碰撞影响
循环流化床内烟气脱硫模拟分析
颜 岩 彭晓峰 王补宣
(清华大学热能工程系 北京 100084)
INVESTIGATION ON FLUE GAS DESULFURIZATION IN A CIRCULATING FLUIDIZED BED
YAN Yan, PENG Xiao-feng, WANG Bu-xuan (Thermal Engineering Department, Tsinghua University, Beijing 100084, China)
引入颗粒相湍能和湍能耗散率描述颗粒团的脉动 分析不同表观气速和颗粒循环量时床内颗粒相轴 及颗粒团形成和脉动对浆滴蒸 分析指出颗粒团的形成可减小浆
向速度沿径向变化规律 发和脱硫效果的影响 滴蒸发速度
增强脱硫效果
关键词: 双流体模型 循环流化床 烟气脱硫 颗粒团