:

*
[email protected]

:
.
.
.
.
.
Shooting .
:

Unreacted shrinking core

.
.
.
] .[

.

Shaft Furnace

.


.
)Fe(s) + HO (g
)(
)FeO (s) + H (g
)FeO (s)+ CO(g
)Fe(s) + CO (g
)(

].[
.
.
.
.

) (

) (
.
:
.

](kg>>D//r) [
.
.
.
:
)(

Cj j ij

i j .

) ( :

i
j Xi

C j = C 0j +

:
)(2

X 1 = C H0 2 C H 2

)(3

0
X 2 = CCO
CCO

:
0
X 2 = 3(C Fe
2O 3

)(
] [ ) (rc
C Fe 2O3 ) = X 3

X rs = X 1 +

) ( .
X 3M W 1 / 3
)
n p 4

)(

rc = (r03

r np Mw
.
:

)(6
)(7
)(8

)(9

)(10

=0

dX 3
+ n p ( R1 ( X 1 , X 3 ) + R 2 ( X 2 , X 3 )) = 0
dz
np
dTsol

dz
) Gmsol ( X 3 )Cpsol ( X 3 , Tsol

u sol

H (T

) sol ) Ri ( X i , X 3 , Tsol

dX 1
+ n p R1 ( X 1 , X 3 ) = 0
dz
dX 2
+ n p R 2 ( X 2 , X 3 ) = 0
u
dz
dTg
) n p Ap h(Tsol Tg

=0
) dz Gmg Cp g ( X 1 , X 2 , Tg
u

Ap h(Tsol Tg )

:
X 2 ( z = L) = 0

X 1 ( z = L) = 0

Tg ( z = L) = Tgin

X 3 ( z = 0) = 0
Ts ( z = 0) = Tatm

)(11

X X Tg ) z= ( )(
) ( z=L .
:
)(12

4rc2Cn
) (1 / kn + rc / Dn rc2 / r0 Dn

= R

n= H CO . :
)(13

R = n p R

.
X :
)(14

+ (r03 rc3 ) ps / M wps

/ M wps

rc3 ) ps

:
)(15

wrs

[1 / 2r / M
[r / M

3
wrs + ( r0

3
c rs

3
c rs

L
Gmsol (rc ) = Gmsol

rs
ps
3
3
3
) p rs rc + C p ps ( r0 rc
) rs rc3 + ps (r03 rc3

[C

= C Ps

.
) (h ) ( :
)(16

1
Pr 2 / 3
C p Gmg

=h

Pr .
R [] E a1 / R g = 179.14 R
[] Ea 2 / Rg = 342.43.
z
MATLAB .

].[
) (
.
) :(

Gas
Gas
Nm/h

rate

flow

,
,
,
,
,
C

)Inlet composition(at z=L

,
,
,
,

length

Solid
Production
t/h
Mineral pellet density
g/cm
Sponge iron density
g/cm
Pellet ratio
cm
np
pellets/cm

zone

Reactor
Reaction
cm
Diameter
cm

:
) ( ].[
) :(

(,exp
(,exp
, T, cm/s
, T, cm/s
cal/cm/K

K
/,/T) cm/s
K
/,/T) cm/s
D
D
h

) ( Shooting
.
) :(

data
(Dry

,
,
--,
,
,

siderca
model data
(Dry base) %
(Dry base) %

model data
base) %

orthogonal collocation
Shooting
)outlet gas composition(z=
H
,
,
CO
,

HO
----, ,
CO
,
CH+N
, ,
,
metallization
,
,
)CPU Time(min
---

:
.( )
:()

Gas
Gas flow rate
Nm/h
Inlet composition(at z=L)
H
,%
CO
,%
HO
,%
CO
,%
CH+N
,%
pressure
atm
Solid
Production
t/h
Mineral pellet density
g/cm
Sponge iron density
g/cm
Pellet ratio
cm
np
pellets/cm
Reactor
Reaction zone length
cm
Diameter
cm

,
,
,
,
,

:
) ( ].[
) :(
(,exp
(,exp
, T, cm/s
, T, cm/s
cal/cm/K

K
/,/T) cm/s
K
/,/T) cm/s
D
D
h

) ( .
) :(
model
)(Dry base

gilmore data
model data
(Dry base) %
(Dry base) %

data
%

Orthogonal

,
,
,
,
,
,

collocation
Shooting
)outlet gas composition(z=
H
,
,
CO
,
,
,
HO
,
CO
,
,
,
CH+N
,
metallization

,
)CPU Time(min
--

.
.

[] W.D Munro, N.R. Amundson, Solid-fluid heat exchange in moving beds, Ind. Eng. Chem, ()
.
[] N.R. Amundson, Solid-fluid interaction in fixed and moving beds, Ind. Eng. Chem, () .
[] C.W. Siegmund, W.D. Munro, N.R. Amundson, Solid-fluid interaction in fixed and moving beds. Two
problems on moving beds. Ind. Eng. Chem, () .
[] R.J. Schaefer, D. Vortmeyer, C.C. Watson, Steady state behavior of moving bed reactors, Chem. Eng.
Sci. () .
[] H. Yoon, J. Wei, M.M. Denn, A model for moving bed coal gasification reactors. AIChE J. ()
.
[] N.R. Amundson, L.E. Arri, Char gasification in counter current reactor, AIChE J. () .
[] Takenaka Y., Kimura Y., Narita K. and Kaneko D., Mathematical model of direct reduction shaft
furnace and its application to actual operation of a model plant, Comp. and Chem. Engng. (),.
[] Towhidi N. , Direct Reduction, V., Tehran University press,
[] Patisson F. and Ablitzer, Physicochemical and thermal modeling of the reaction between a porous
pellet and a gas, Powder Technology () .

[]

Daniel R. Parisi, Miguel A. Laborde , Modeling of counter current moving bed gas-solid reactor
used in direct reduction of iron ore, Chemical Engineering journal () .