3D numerical simulation of gas-solid hydrodynamics and coal combustion in an industrial scale circulating fluidized bed Massoud Massoudi Farid Hyo Jae Jeong Jong-Min Lee (KEPCO) Dong-Won Kim (KEPCO) Jungho Hwang Department of Mechanical Engineering Yonsei University
47
Embed
3D numerical simulation of gas -solid hydrodynamics and …. M... · · 2017-02-213D numerical simulation of gas -solid hydrodynamics and coal combustion in an ... Solid and fluid
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
3D numerical simulation of gas-solid hydrodynamics and coal combustion in an industrial scale circulating fluidized bed
Massoud Massoudi Farid Hyo Jae Jeong
Jong-Min Lee (KEPCO) Dong-Won Kim (KEPCO)
Jungho Hwang
Department of Mechanical Engineering Yonsei University
Introduction
Circulating Fluidized Bed
After 10 seconds no combustion occurred ( it took around 10 days)
Cyclone
Return leg
Furnace (Riser)
Air inlet
Wing Wall tubes
Steam inlets
Coal Feeder
1) Air enters the furnace through nozzles located at the bottom of the furnace and fluidizes sand particles which are packed initially in the furnace.
2) Air reacts with coal particles which are injected to the furnace by coal feeders located around the furnace.
3) Sand particles carried by mixture of air and combustion products to the top area of the furnace, enter the cyclones and are separated from the combustion products.
4) Combustion products then leave the cyclones through exit ports and sand particles return to the furnace through return legs.
5) Water steam flow into wing wall tubes, receive heat from gas-solid mixture inside the furnace.
Eulerian-Eulerian
In Eulerian frame
In Eulerian frame
Assumptions 1) It is a multiphase flow method and N.S equations are solved for both phases. 2) Solid and fluid phases are treated in the same way as interpenetrating continua. 3) granular properties are added to solid phase by solving an additional equation for solids
fluctuating energy or granular temperature. 4) There is no limitation for volume fraction. Hence effects of volume fraction should be considered in N.S equations
advantages: 1) It predicts well main features
of the flow. 2) lower computational time
compared to other methods
Multiphase Eulerian-Granular Method(MEGM) (For interaction of Gas and Sand particles)
Discrete Phase Method(DPM) (For interaction of Gas and Coal)
Assumptions 1) Solid particles have low Volume fraction
and they can not affect fluid motion. Hence Volume fraction is not considered in N.S equations.
2) There is no particles interaction. 3) Trajectory of a discrete phase particle is
predicted by integrating the force balance equation of the particle.
Eulerian-Lagrangian
In Lagrangian frame
In Eulerian frame
Particle force Balance for Solid Phase
xp
ppD
p Fg
uuFt
u +
−+−=
∂∂
ρρρ )(
)( 24Re18
2D
ppD
Cd
Fρ
µ=
Momentum of Fluid Phase
)()()( otherDPM FFgpvvvt
++ρ+τ⋅∇+−∇=ρ⋅∇+ρ
∂∂
Boiler Geometry
After 10 seconds no combustion occurred ( it took around 10 days)
5 Super heater wing walls
1 Evaporator wing wall Due to
computational barriers
Real CFB
340 MWe CFB boiler located in Yeosu, Korea.
Boundary Conditions and Coal properties
After 10 seconds no combustion occurred ( it took around 10 days)
Once air was introduced to the bottom of the furnace, the bed began to expand (around t=0.3s) , with some particles hitting the top wall and falling down and some particles close to the cyclone were dragged by gas into the cyclone. At t=35s bed height reached to the highest level and from then onward the bed height remained almost constant which shows the quasi-steady state situation where the input and output rates of particles for each of the flowing units involved are equal.
Bed expansion by time steps
0
2
4
6
8
10
12
14
16
18
20
22
24
26
28
30
32
34
36
0 0,02 0,04 0,06 0,08 0,1 0,12 0,14 0,16 0,18 0,2
Sand Volume Fraction
H
Sand-Volume fraction
In the riser, solid volume fraction is high at the bottom and low at the top (the so called S-Shaped Profile)
1) There is a dense bottom coexisting with dilute top in both the furnace chamber and the return leg.
2) A clear boundary between solid phase and gas phase in the return leg can be discerned which is a bubbling fluidized bed, and packed at the beginning of the simulation to serve as solid seal to prevent riser gas leakage via the return leg.
3) The solids volume fraction was normally high near the walls and low in the center of the furnace, except from some clusters,
1
2
Grace Et. al
Sand-Velocity-Vectors (colored by sand Z velocity value )
Largest Company
1) Wherever there is a velocity vector it means that there is a particle at that point.
2) Solids vertical velocity is mainly positive in those areas where volume fraction is
low and negative those areas where volume fraction is high or on the other words negative velocity owes to clustering of particles (points 1 and 2)
1
2
Core-annulus
Largest Company
Z=15m Z=20m Z=30m
1) The volume fraction curves are typical to the so-called core–annulus structure, which causes solid back-mixing.
2) The core–annulus structure can be confirmed by Z velocity curves showing a downward movement of solids near-wall region, while rising particles in the center.
3) At Z=15 m the cross section is under the wing wall tubes and core-annulus structure is more obvious.
4) At Z=20 m the cross section is at the bottom part of wing wall tubes. At this point the core-annulus structure is affected little by wing wall tubes.
5) There are some local core annulus between wing wall tubes.
6) At Z=30 solid concentration is low and solid concentration is higher at the entrance of the cyclone. Hence core-annulus structure is weaker compared to previous pictures but still near the walls and wing wall tubes solid concentration is high and solid Z velocity is negative.
Gas-Velocity
ATmV total
cal ×ρ=
)(
Z=20m Varea-av-num=7.24m/s
Vcal =5.38m/s based on T=1356.6 K ρ=0.255936
Z=30m Varea-av-num=6m/s
Vcal =5.23m/s based on T=1314 K ρ=0.262421
Z=10m Varea-av-num=7.2m/s
Vcal =4.91m/s based on T=1291 K ρ=0.280238
1) In some regions due to existence of solid particles gas velocity is higher than other regions and in general gas velocity is higher at the middle of the furnace and lower near the walls due to no slip boundary condition. At the entrances and the exit of the cyclone gas velocity is high.
2) Theoretically velocity can be calculate by knowing inlet mass flow rate, furnace cross section and averaged density.
3) Differences between theoretical values and numerical ones could be due to existence of solid particles, coal combustion, and gas turbulence flow.
Sand-Velocity-Cyclone (colored by sand Z velocity value )
Largest Company
Y=-3.1 Y=-2.5 Y=0 Z=35
Most particles congregate at the wall immediately after the inlet and then descend in strands in a cyclone. also there is particle accumulation in the apex region where particle–particle interaction is intensive. These phenomena have been reported before by other authors.
Gas velocity vectors in cyclone (colored by gas Z velocity)
Largest Company
%96.99)/(157)/(06.01100)1( =−=×−=
skgskg
mm
efficiencyCycloneinletletSand
outletSand
Largest Company
Pressure Loop and solid circulation rate
a
b
c d
f
g
e
inlet Cyclone inlet Cyclone outlet
numerical real numerical real numerical real
P(pa) 9778.273 8002 -1791.72 0 -1304 -1304
Gs(kg/s-m2) 3.31
h
Pcyc-en=0 pa
Pnuz-out=8002 pa
a
b
c
e
f
g
d
0
1
2
3
4
51 56 61G
s t
0
1
2
3
4
5
0 20 40 60
Gs
t
Combustion Results
Largest Company
Gas Temperature.
1) Gas Temperature is high at the bottom of the furnace due to combustion reactions.
2) Due to heat transfer with wing walls and furnace wall, temperature falls and reaches at 1095K at the cyclone outlet.
3) The outlet Temperature in real CFB is around 1127K; therefore the numerical error is only around 2.84%.
Z=10m
Z=20m
Z=30m Coal Feeder No.2 Y=3.024
Coal feeder No.1 Y=0.56
CO Mole fraction
Largest Company
1) Due to complete combustion and conversion of CO to CO2 the amount of CO in the domain is generally low.
2) The majority of CO exists at the bottom of the furnace and near the coal feeders due to large number of carbon particles and de-volatilization.
Z=10m
Z=20m
Z=30m
Largest Company
CO2 Mole fraction.
1) The amount of CO2 in the domain is high especially at the upper side of the furnace due to conversion of CO to CO2 by water gas shift reaction and CO and O2 reaction, and also volatile combustion.
2) The reason that at X=0 cross-section CO2 mole fraction is higher at the right side is because of the position of the coal feeders.
Coal Feeder No.2 Y=3.024
Coal feeder No.1 Y=0.56
Z=10m
Z=20m
Z=30m
Largest Company
Volatiles Mole fraction
Similar to CO mole fraction the amount of coal volatiles mole fraction is low in the domain due to burning and converting to combustion products.
Z=10m
Z=20m
Z=30m
Largest Company
Char Concentration
the concentration of char is high at the bottom of the furnace and near the coal feeders.
Largest Company
Outlet dry based spices
N2 O2 CO2 H2 CO CH4 volatiles Mol % 79.673 8.261 12.066 0 0 0 0
0
20
40
60
80
N2O2
CO2
Conclusion
• Eulerian-Eulerian method was used (multiphase granular model) for modeling sand and gas hydrodynamics
• Hydrodynamic simulation results such as sand circulating rate, pressure loop, and sand distribution in furnace were reasonable.
• 3D numerical simulation of air and sand particle hydrodynamics along with coal combustion in an industrial scale circulating fluidized bed was done.
• Eulerian-Lagrangian method was used (discreet phase model (DPM)) for modeling coal particle trajectories
• Combustion simulation results such as combustion products distribution in the furnace were reasonable.
• Exit gas temperature was in good agreement with real data obtained from a 340 MWe CFB boiler located in Yeosu, Korea.