中国有色金属学报 2004,(11),1918-1925 DOI:10.19476/j.ysxb.1004.0609.2004.11.022
固定床中对流扩散与反应耦合问题的一个解析解
翟玉春 纪智玲
东北大学资源与土木工程学院,东北大学材料与冶金工程学院,东北大学材料与冶金工程学院 沈阳110004 ,沈阳110004 ,沈阳110004
摘 要:
考虑到冶金固定床是颗粒堆积的松散多孔介质, 首先将单一颗粒的动力学方程推广到多颗粒的微元体上, 导出了化学反应项的具体形式, 建立了冶金固定床一维对流、反应、扩散方程, 然后求出其解析解;通过实例计算分析主气流中反应气体的浓度分布。结果表明:颗粒大小和孔隙率对反应进程有很大的影响, 其影响可以用Thiele数和Peclet数定量衡量;微元体的反应类型不等于反应器的反应类型;化学反应速率对反应气体浓度的影响超过对流速度, 而扩散的影响最小;微元体的综合反应速率、对流速度以及反应器的长度本身是影响反应器反应类型的重要因素。
关键词:
颗粒动力学 ;微元体动力学 ;化学反应 ;对流 ;扩散 ;解析解 ;
中图分类号: TF01
收稿日期: 2003-11-04
基金: 国家自然科学基金重点资助项目 (50136020); 教育部重点资助项目 (01056);
An analytical solution for coupled convection-diffusion and chemical reaction in fixed-beds
Abstract:
In view of the fact that fixed bed is porous media of granule aggregation, the dynamics of a single particle is extended to the dynamics of representative-element volume including multiple particles, so that the more exact chemical reaction term can be deduced. Then one-dimensional model for instable convection, chemical reaction and diffusion in fixed bed is proposed and solved analytically. Computation cases show that the scale of pellet and the porosity of bed have great impact on the reaction process taking place in reactor, which can be measured by Thiele module and Peclet module, the reaction behavior of REV is not the same as that of reactor. For the gas concentration, chemical reaction has greater impact than convection has, because Thiele module is proportional positively to L2 /D and Peclet module is to L/D, and diffusion has the least impact. The overall rate of REV, length of reactor itself and convection all contribute remarkably to the performance of reactor.
Keyword:
dynamics of single particle; dynamics of representative-element volume; chemical reaction; convection and diffusion; analytical solution;
Received: 2003-11-04
气固反应是冶金工业中的主要反应类型, 固定床中的固相物质是人工制备的, 因此固定床可简化为均匀颗粒堆积而成的松散多孔介质。 通常关于颗粒动力学的研究是在主气流中反应气体的浓度已知的情况下进行的。 但在固定床中, 反应气体的浓度是待定的, 且气体浓度也受到反应的影响。 因此在反应器中化学反应与气体的传质过程是耦合的。
一般来说对流、 扩散和反应对反应器中的气体浓度分布都有影响。 许多文献给出了传质微分方程的通用形式
[1 ,2 ,3 ,4 ]
, 但具体的研究极少。 Loureiro等
[5 ]
研究了固定床中浓度波的传播, 其控制方程是略去了扩散项的单波方程。 Douglas等
[6 ]
利用单波方程研究了气体吸附以及对流对固定床中浓度波传播的影响。 Loureiro等
[7 ]
还研究了求解单波方程的两种数值方法。 与此同时, Ferreia等
[8 ]
研究了不存在对流作用的反应扩散方程。 Winkalman等
[9 ]
用反应扩散方程研究了稳定条件下吸附和一级化学反应同时发生的问题。 Whitaker等
[10 ,11 ]
用反应扩散方程研究无对流时固定床中的气体运动。 冶金中关于固定床等温反应的研究大多不考虑扩散
[12 ,13 ,14 ]
。 因此同时考虑对流、 反应、 扩散研究固定床式反应器中反应的工作极少见。 鞭岩利用对流反应扩散方程研究了填料床中的反应转化, 但只考虑了稳定问题
[15 ]
。
为了研究对流、 扩散与反应的耦合作用, 首先要给出化学反应项的准确形式。 冶金动力学对单颗粒的综合反应速率已作了详尽的研究。 但颗粒不能全面反映床层的性质, 如孔隙率、 渗透性等。 而这些性质对反应器中发生的反应、 对流、 扩散有很大的影响。 如果用足够多的颗粒的集合体取代单一颗粒作为反应器的表征性体积单元, 即微元体, 可以克服上述困难。 这是因为在微元体上可以定义孔隙率、 渗透性、 对流速度、 气体浓度等参数和变量。 但这首先需要将单颗粒动力学推广到微元体上。 在颗粒为等径圆球和密集堆积的条件下, 本文首先将颗粒的综合动力学推广到微元体上, 然后建立对流、 反应、 扩散模型, 求出了解析解; 并进行了实例计算与分析。 结果表明: 微元体的反应类型不等于反应器的反应类型; 反应器自身特性和外部因素对反应器的反应类型都有重要影响, 并对这些影响作了详细的分析。
1 微元体的反应动力特性
考虑等温条件下无孔隙颗粒的气固反应, 反应前后颗粒大小不变, 如式 (1) 和图1所示
A (g) +bB (s) =dD (s) (1)
式中 g表示气体, s表示固体。 将单一颗粒综合速率方程推广到微元体, 首先要给出微元体内孔隙内表面积的大小。 按晶体学
[16 ]
的研究, 等径圆球的密集堆积只有三种方式, 如图2所示。 并且在密堆积条件下的孔隙率和致密度是常数, 与球的大小无关。 因此可取一个胞体研究孔隙率、 固体致密度和孔隙内表面积, 其结果与取一个微元体相同。
由图2可知, 对于体心立方堆积, 一个胞体中有两个颗粒, 因此其致密度为
n
s
=
3
π
/
8
≈
0
.
6
8
, 孔隙率n =1-n s =0.32。 对于面心立方堆积, 致密度
n
s
=
2
π
/
6
≈
0
.
7
4
, 孔隙率n =1-n s =0.26。 对于密排立方堆积, 致密度为
n
s
=
2
π
/
6
≈
0
.
7
4
, 孔隙率n =1-n s =0.26。 同样可以取一个胞体研究颗粒表面积和孔隙内表面积。 例如, 体心立方堆积, 一个胞体中有两个颗粒, 因此颗粒表面积为f b =2×4πr 2 s , 胞体的边长
a
b
=
4
r
s
/
3
。 这样颗粒表面积与胞体体积之比就是比表面积
图1 形成固相产物的气固相反应过程 (缩核模型)
Fig.1 Gas-solid reaction process (Shrinking core model)
图2 等径圆球密集堆积
Fig.2 Idealized aggregation of pellets with uniform radius
S
b
=
2
×
4
π
r
s
2
(
4
/
3
r
s
)
=
3
(
1
-
n
)
r
s
?
?
?
(
2
)
式中 n 是孔隙率。 按相同的方法可以导出, 面心立方堆积和密排六方堆积的孔隙比表面积也可以用式 (2) 表示, 只要注意对不同的堆积方式孔隙率是不同的。 由于1-n 是固相物质的致密度, 因此上式表明单位多孔介质内的孔隙内表面积与固相物质致密度成正比, 与颗粒半径成反比。 这是易于理解的。 式 (2) 与多孔介质体积的大小无关。 这意味着与微元体内所包含的颗粒多少无关。 也与微元体的形状无关。 按照连续介质的观点, 合理的表征性体积单元取决于颗粒和所研究问题的相对尺度, 既没有绝对的大小也没有固定形状。 因此式 (2) 是孔隙比表面的合理表达式。 对非等径圆球、 非密集堆积和随机堆积, 文献
[
17 ]
利用充填机理 (filling mechanism) 和共据机理 (occupation mechanism) , 研究了颗粒堆积体的孔隙率和粒径分布的关系, 不仅可以给出定性关系, 还可以给出定量描述。 在文献
[
17 ]
的基础上, 可以研究更复杂堆积条件下单位体积多孔介质内颗粒的比表面积, 这将另文叙述。
在冶金动力学中颗粒的综合速率常以反应气体的浓度表示, 一些典型的颗粒综合速率方程, 比如外扩散控制、 化学反应控制和内扩散控制时颗粒的综合速率方程可由下式所示
[17 ]
:
-
d
G
A
d
t
=
4
π
r
s
2
k
g
C
A
=
4
π
r
c
2
k
C
A
=
4
π
D
e
r
s
r
c
r
s
-
r
c
?
?
?
(
3
)
式 (3) 的含义是单一颗粒的在单位时间因化学反应而消耗的反应气体的量。 式中G A 为反应气体A的摩尔数; r s 为颗粒半径; r c 是未反应核的半径; C A 为反应气体A的摩尔浓度; k g 为气体A的外扩散系数; k 是化学反应速率常数; D e 是固相产物层内气体的有效传质系数。 式 (3) 中4πr
s
2
和4πr
c
2
分别是颗粒表面积的大小和反应界面的大小。 为了将颗粒综合速率方程 (3) 推广到微元体上只需将颗粒的面积因素用微元体的比表面代替即可。 因此利用式 (2) , 则与式 (3) 对应的微元体综合速率为
-
d
G
A
d
t
=
3
(
1
-
n
)
r
s
k
g
C
A
=
3
(
1
-
n
)
r
s
Ζ
2
k
C
A
=
3
(
1
-
n
)
r
s
D
e
Ζ
1
-
Ζ
C
A
?
?
?
(
4
)
式 (4) 是单位体积多孔介质, 在单位时间内因化学反应所消耗的反应气体的量。 当综合速率由化学反应与外扩散联合控制时, 颗粒的综合反应速率以及在微元体上的推广为
-
d
G
A
d
t
=
4
π
r
c
2
1
/
k
+
1
/
k
g
=
3
(
1
-
n
)
r
s
Ζ
2
1
/
k
+
1
/
k
g
?
?
?
(
5
)
式中 Z =r c /r s , 是未反应核半径和颗粒半径之比, 也可以用来表征固相反应物的转化率。 由此可见利用微元体比表面的表达式 (2) 可以使我们摆脱从颗粒的综合速率方程出发研究反应器中的反应和传递过程, 同时又可利用冶金动力学关于颗粒反应动力学的现有成果。 本文作者研究外扩散控制的情况时, 微元体的综合速率方程由式 (4) 表示。
2冶金固定床的一维对流-反应-扩散方程及其解答
研究一维填充床。 若微元体的孔隙率为n , 反应是一级的, 对流速度v 是常量。 考虑到反应气体是被消耗的, 因此反应气体A的质量守恒方程为
Μ
A
?
C
?
t
=
Μ
A
D
e
?
2
C
?
x
2
-
Μ
A
v
?
C
?
x
-
R
A
?
?
?
(
6
)
式中 D e 是气体在多孔介质孔隙通道中的有效扩散系数, 即D e =D A n /τ , D A 是扩散系数, τ 是曲折因子; v 是混合气体在多孔介质孔隙通道中流动的真实速度, 与单位面积多孔介质上的平均速度的关系为vn =u ; R A 是单位时间内在单位体积微元体上因反应而消耗的反应气体的质量; M A 是反应气体A的摩尔质量。 式 (4) 的第一式即是外扩散控制条件下, 单位时间内在单位体积微元体上因反应而消耗的反应气体的质量, 以摩尔计。 若以质量表示, 则有
R
A
=
-
Μ
A
d
G
A
d
t
=
3
(
1
-
n
)
r
s
Μ
A
k
g
C
?
?
?
(
7
)
将上式代入式 (6) 可得
Μ
A
?
C
?
t
=
Μ
A
D
e
?
2
C
?
x
2
-
?
?
Μ
A
v
?
C
?
x
-
3
(
1
-
n
)
r
s
Μ
A
k
g
C
?
?
?
(
8
)
将反应气体的浓度用质量分数表示为, C * =M A C /ρ , C * 无量纲; ρ 是混合气体的密度; 则反应气体A的质量守恒为
?
C
*
?
t
=
D
e
?
2
C
*
?
x
2
-
v
?
C
*
?
x
-
3
(
1
-
n
)
n
r
s
k
g
C
*
?
?
?
(
9
)
式中 x 是空间坐标, t 是时间。 D e 是床层中反应气体的有效扩散系数。 若反应器长度用L 表示, 则可引入如下无量纲自变量
x
?
=
x
/
L
,
t
?
=
D
t
/
L
2
。 将其代入式 (9) 可得
?
C
*
?
t
?
=
D
e
?
2
C
*
?
x
?
2
-
v
?
?
C
*
?
x
?
-
μ
?
C
*
?
?
?
(
1
0
)
式中
v
?
=
v
L
/
D
是反应器的Peclet 数, 反映了对流速度与扩散速度的相对大小; 而反应器的Thiele 数为
μ
?
=
3
(
1
-
n
)
L
2
k
g
/
(
n
r
s
D
e
)
反映了化学反应速率与扩散速度的相对大小。 由化学反应计量关系可以导出以固相反应物B表示的综合速率和B的转化率
-
d
G
A
d
t
=
-
1
b
d
G
B
d
t
,
Ζ
3
=
1
-
3
b
Μ
B
ρ
B
r
s
k
g
∫
0
t
C
*
d
t
?
?
?
(
1
1
)
式中 G B 是固相物质B的摩尔数, ρ B 是其密度, M B 是其摩尔质量, b 是其化学计量系数。 将C * 对时间积分是考虑到在反应器内气体浓度是空间位置和时间的函数。 式 (11) 表明只要求出了反应气体浓度, 即可求出固相物质的反应速度和转化率。 考虑反应气体的如下初边值问题
?
C
*
?
t
?
=
?
2
C
*
?
x
?
2
-
v
?
?
C
*
?
x
?
-
μ
?
C
*
C
*
(
x
,
0
)
=
C
∞
*
C
*
(
0
,
t
)
=
C
*
(
?
C
*
?
x
+
h
?
C
*
)
|
x
=
L
=
h
?
C
∞
*
h
?
=
h
L
D
,
h
?
∞
=
h
C
∞
L
D
}
?
?
?
(
Ⅰ
)
为求问题 (Ⅰ) 的解析解, 将
C
*
(
x
?
,
t
?
)
分解为稳定和非稳定两个部分
C
*
(
x
?
,
t
?
)
=
C
1
*
(
x
?
)
+
C
2
*
(
x
?
,
t
?
)
?
?
?
(
1
2
)
令稳定部分带有非齐次边界条件, 非稳定部分带有齐次边界条件。 由于
v
?
、 D 和
μ
?
总大于零, 容易导出稳定问题的特征根λ 1 、 λ 2 是实根, 因此容易求得
C
1
*
(
x
)
=
A
1
e
λ
1
x
+
A
2
e
λ
2
x
,
λ
1
=
v
?
2
+
v
?
2
+
4
u
?
2
,
λ
2
=
v
?
2
+
v
?
2
+
4
μ
?
2
?
?
?
(
1
3
)
A
1
=
C
0
e
λ
2
(
λ
2
+
h
?
)
-
h
?
∞
e
λ
2
(
λ
2
+
h
?
)
-
e
λ
1
(
λ
1
+
h
?
)
,
A
2
=
h
?
∞
-
C
0
e
λ
1
(
λ
1
+
h
?
)
e
λ
2
(
λ
2
+
h
?
)
-
e
λ
1
(
λ
1
+
h
?
)
?
?
?
(
1
4
)
对于非稳定问题C * 2 的解, 可引入一个新函数
w
?
(
x
?
,
t
?
)
, 求得
w
?
(
x
?
,
t
?
)
和
C
2
*
(
x
?
,
t
?
)
为
w
?
(
x
?
,
t
?
)
=
∑
m
=
1
∞
1
Ν
(
β
m
)
sin
β
m
x
?
e
-
β
m
t
?
?
?
?
∫
0
1
sin
β
m
x
?
[
h
?
C
∞
-
C
1
(
x
?
)
]
e
-
v
?
2
x
?
d
x
?
?
?
?
(
1
5
)
C
2
*
(
x
?
,
t
?
)
=
w
(
x
?
,
t
?
)
exp
{
v
?
2
x
?
-
(
v
?
2
4
+
μ
?
)
t
?
}
?
?
?
(
1
6
)
1
Ν
(
β
m
)
=
2
[
β
m
2
+
(
v
?
2
+
h
?
∞
)
2
]
L
[
β
m
2
+
(
v
?
2
+
h
?
∞
2
)
]
+
(
v
?
2
+
h
?
∞
)
?
?
?
(
1
7
)
式中 β m 是下列方程的正根
(
h
?
∞
+
v
?
2
)
/
β
m
+
cot
(
h
?
∞
+
v
?
2
)
=
0
?
?
?
(
1
8
)
将式 (14) ~ (17) 代入式 (12) 即得一维固定床中反应气体A的浓度分布。
3 计算实例与分析
下面通过实例计算分析对流、 反应、 扩散对气体输运及反应器特性的影响, 本文作者主要从传输角度进行讨论。 计算中取D =0.000 3 m2 /s, C 0 =0.5, C ∞ =0, h =0.15 m/s, 反应器入口端气体压力P 0 =1.1 MPa, 反应器出口端气体压力P L =1.0 MPa, 其它参数在计算中给出。
颗粒的外扩散系数k g , 已得到较充分的研究, 可用如下关联式估算
[18 ]
S
h
=
2
.
0
+
0
.
6
R
e
1
/
2
S
c
1
/
3
,
S
h
=
k
g
r
s
D
e
,
R
e
=
2
r
s
ρ
v
μ
,
S
c
=
μ
ρ
D
e
?
?
?
(
1
9
)
式中 S h 是Sherwood数, R e 是Renold数, S c 是Schmidt数; μ 是流体的动力粘性系数; 对于颗粒堆积的多孔介质, 可从下式估算渗透率k 和渗透系数K
[20 ]
k
=
ε
3
4
5
(
1
-
ε
)
2
r
s
2
,
Κ
=
k
μ
?
?
?
(
2
0
)
若流体不可压缩, 可导出渗流速度为
v
=
-
Κ
p
L
-
p
0
L
=
c
o
n
s
t
。
从式 (19) 、 (20) 可以看出颗粒半径r s 和孔隙率n 对外扩散系数k g 和渗透系数K 有影响, 因此它们对反应进程有影响, 这种影响可用Thiele数和Peclete数衡量。 为了简便起见, 在图表中的Thiele数简记为Th , Peclet数简记为Pe 。 图3 (a) 、 (c) 所示分别是r s 对Thiele数、 Peclet数和它们的比值Th /Pe 的影响。 从图中可以看出, 当r s 增加时, Thiele数减小, Peclet数增大。 Th 的减小是因为颗粒越大则比表面积越小, 因而综合速率下降的缘故; Pe 的增大则是因为颗粒越大, 渗透系数越大, 对流加快的缘故。 图3 (b) 、 (d) 所示分别是孔隙率对Th 、 Pe 和Th /Pe 的影响。 从图中可以看出, n 增大时, Th 数减小, 这同样是因为比表面积减小的原因, Pe 数的增加是因为孔隙率增加, 渗流阻力减小, 对流加快引起的。
图3 颗粒半径、 孔隙率对Thiele数和Peclet数的影响
Fig.3 Effect of pellet radius and porosity on Thiele module and Peclet module
图4所示是r s 对气体浓度的影响。 从图中可以看出, r s 越大, 气体浓度越大。 这是因为当r s 增大时, Th /Pe 数减小, 综合速率下降, 消耗的气体减小, 同时对流加快, 由其携带进入反应器的气体增加所引起的。 图5所示是n 对气体浓度的影响, 从中可以看出, n 越大气体浓度越高, 这同样是因为n 增加时, 比值Th /Pe 下降, 对流影响上升, 而反应影响下降所引起。 图3~5表明, 颗粒半径和孔隙率是影响气体浓度和反应进程的敏感因素, 但它们的作用是通过对反应和对流的影响间接体现出来的, 因此下面直接讨论综合速率、 对流和反应器长度的影响。
图4 颗粒半径对气体浓度分布的影响
Fig.4 Effect of granule radius on concentration of reactive gas (n =0.32; L =10 m)
图6所示是反应气体浓度随时间的变化。 从图中可以看出, 随时间的发展气体浓度上升, 这主要是因为对流速度比综合反应速率大许多, Th /Pe 很小, 因此反应无法及时消耗由对流带入反应器的反应气体, 从图6还可以看出对流对浓度的非稳态分布有明显的影响, 对流较快时浓度随时间的变化更明显。
图5 孔隙率对气体浓度分布的影响
Fig.5 Effect of porosity on concentration of reactive gas (r s =0.02 m; L =10 m)
图6 气体浓度的非稳定分布
Fig.6 Instable concentration distribution of reactive gas (n =0.32; r s =0.02 m; L =10 m; k g =0.001 5 m/s)
图7所示是综合速率对浓度分布的影响。 由于其他参数不变, 外传质系数的变化决定了综合速率的变化; 从图中可以看出, 综合速率对反应转化进程有明显的影响, 综合速率越大, 反应气体浓度越小。 因此, 当综合速率很大时, 只有加快对流, 才能进一步提高反应速率。 综合速率对反应进程的影响可用Thiele数衡量, 图7给出了相应的Thiele数。
对流对气体浓度也有明显的影响, 对流越快, 由其携带的反应气体越多, 气体浓度越高, 图8示意地显示了对流的这个特性。 因此, 当对流很快时, 只有提高综合速率, 才能提高转化率。 对流对气体浓度的影响可用Pelect数衡量。
图7 综合速率对浓度分布的影响
Fig.7 Effect of overall rate on concentration distribution of reactive gas (v =5.025 2 m/s; n =0.32; L =10 m; r s =0.02)
图8 对流对气体浓度的影响
Fig.8 Effect of convection rate on gas concentration (r s =0.01 m; n =0.32; k g =0.001 518 m/s; L =10 m)
图9所示是反应器长度对气体浓度的影响, 采用无量纲长度, 可将不同长度的反应器画在同一张图上, 以方便对比。 此时图中同一坐标上的气体浓度并不是不同反应器同样位置上的气体的浓度, 而是同样比例位置上的气体浓度。 由于反应器长度不同, 此时逐点比较是无意义的, 但对反应器整体特征的比较仍是有意义的。 从图9可以看出, 反应器较短时, 对流的影响大, 气体浓度高, 此时只有提高综合速率, 便提高了反应转化率; 而当反应器长度增加时, 对流的作用下降, 气体浓度低。 由于进入反应器的反应气体的量是相同的, 反应器长度增加, 意味着气体与固体反应的机会增多。 表1给出了反应器长度与气体转化率的关系, 可以看出当L 增大时, 转化率上升。 因此, 上述现象是可以理解的。 由于Thiele数与反应器长度的平方成正比, Pelect数与长度的一次方成正比。 因此由Thiele数和Pelect数的比值定量地度量反应器长度对其自身性能的影响。 表1还给出了反应器长度对Pelect数和Thiele数的影响。
图9 反应器长度对气体浓度分布的影响
Fig.9 Effect of reactor length on concentration of reactive gas (r s =0.01 m; n =0.32; k g =0.001 5 m/s; v =5.025 2 m/s)
表1 反应器长度对Pelect数、Thiele数和气体转化率的影响
Table 1 Effect of reactor length on Peclet module, Thiele module and gas conversion
L /m
Th /103
Pe /103
Th /Pe
(C 0 -C L ) /C 0
5
1.992 2
2.093 8
0.951 5
0.613 5
10
7.968 7
4.187 7
1.902 9
0.850 7
15
17.930 0
6.281 5
2.854 4
0.942 3
20
31.875 0
8.375 3
3.805 8
0.977 7
4 结论
1) 微元体的反应类型不等于反应器的反应类型。 在等温条件下影响反应器反应特性的内部因素有微元体的综合速率、 反应器长度。 外部因素是对流速度。
2) 颗粒大小和孔隙率是影响反应进程的敏感因素。 它们的作用是通过对综合速率和对流的影响间接反应出来的, 并可以用Thiele数和Peclet数衡量。
3) 当微元体综合反应速率增加一倍时Thiele数增加倍。 而当对流速度增加一倍时, Pecelet数增加L /D 倍。 因此微元体综合反应速率的影响比对流速度更大。 这还表明反应器长度L 也是重要的影响因素。
4) 在扩散、 对流和化学反应中, 化学反应速率的影响最大, 对流的影响次之, 扩散的影响最小。 因此, 理论分析和数值分析时, 准确地给出化学反应项的具体形式至关重要。
化学反应器中, 反应与传递过程的相互作用是十分复杂的现象。 本文仅是初步的研究, 更深入的问题将在后续论文中陆续进行。
参考文献
[1] BirdRB, StewardWE, LightfoodEN.TransportPhenomena[M].NewYork:JohnWileyandSonsInc, 1960.
[2] SzekelyJ.FluidFlowPhenomenainMetalProcessing[M].NewYork:AcademicPress, 1979.
[3] 查金荣, 陈家镛.传递过程原理与应用[M].北京:冶金工业出版社, 1997.312327.ZHAJin rong, CHENJia yong.PrinciplesandApplicationsofTransportProcess[M].Beijing:MetallurgicalIndustryPress, 1997.312327.
[4] 陶文铨.数值传热学 (第二版) [M].西安:西安交通大学出版社, 2001.187237.TAOWen quan.NumericalMethodsinHeatTransfer (2rdEd) [M].Xi an:Xi anJiaotongUniversityPress, 2001.187237.
[5] LoureiroJM, CostaCA, RodriguesAE.Propagationofconcentrationwaveinfixed bedadsorptionractor[J].TheChemicalEngineeringJournal, 1983, 27:135148.
[6] LevanMD, CostaCA, RodriguesAE, etal.Fixed bedadsorptionofgases:Effectiveofvelocityontransitiontypes[J].AICHEJournal, 1988, 34 (6) :9961005.
[7] LoureiroJM, RodriguesAE.Twosolutionmethodforhyperbolicsyestemofpartialdiffentialequationsinchemicalengineering[J].ChemicalEngineeringScience, 1991, 46 (12) :32593267.
[8] RoseM, FerreiraQ, RodriguesAE.Diffusionandcstalytiezero ordorreactioninamacroreticularionexchangeresin[J].ChemicalEngineeringScience, 1993, 48 (16) :29272950.
[9] WinkelmanJGM, BeenackersAACM.Simultaneousabsorptionanddesorptionwithreversiblefirsts ordorchemicalreaction:Analyticalsolutionandnegtiveenhancementfactors[J].ChemicalEngineeringScience, 1993, 48 (16) :29512955.
[10] WhitakerS.Diffusioninparkedbedsofporousparticles[J].AICHEJournal, 1988, 34 (4) :337346.
[11] WhitakerS.Diffusionandreactioninamicroporemacroporemodelofporousmedium[J].LatAmerChemApplChem, 1983, 13:552559.
[12] BischoffKB.Conversionofsolidparticlelesincontinuousfixedbed[J].IntEngChemFundamental, 1969, 8:665.
[13] MoriyamaA.Analysisoffluid solidreactioninfixed bedreactorsunderisothermalconditions[J].ChemEng, 1970, 34:1308.
[14] 肖兴国, 谢蕴国.冶金反应工程学基础[M].北京:冶金工业出版社, 1997.281288.XIAOXing guo, XIEYun guo.BasicEngineeringPrinciplesofMetallurgicalReaction[M].Beijing:MetallurgicalIndustryPress, 1997.281288.
[15] 鞭 岩, 森山昭.冶金反应工程学[M].北京:科学出版社, 1981.7792.MuchiI, MoriyamaA.EngineeringofMetallurgicalReaction[M].Beijing:SciencePress, 1981.7792.
[16] 冯 端.金属物理学 (第一卷, 结构与缺陷) [M].北京:科学出版社, 2000.182206.FENGDuan.PhysicsofMetal (Vol.1, StructureandDefect) [M].Beijing:SciencePress, 2000.182206.
[17] AIBing yu.ThePackingofSolidParticles[D].NewSouthwalse, Australia:TheUniversityofWollongong, 1989.
[18] 莫鼎成.冶金动力学[M].长沙:中南大学出版社, 1987.131135.MODing cheng.KeneticsinMetallurgicalProcess[M].Changsha:CentralSouthIndustrialUniversityPress, 1987.131135.
[19] SohnHY, WadsworthME.提取冶金速率过程[M].北京:冶金工业出版社, 1984.24.SohnHY, WadsworthME.RateProcessesofExtractiveMetallurgyExtraction[M].Beijing:MetallurgicalIndustryPress, 1984.24.
[20] RumerRRJr.ResistancetoFlowthroughPorousMedia, FlowthroughPorousMedia[M].NewYork:AcademicPress, 1960.331376.