InTech Book ChapterF
InTech Book ChapterF
net/publication/221927969
CITATIONS READS
13 4,334
3 authors, including:
All content following this page was uploaded by Natalia Prieto Jiménez on 20 April 2018.
1 Introduction
Fluidization refers to the contact between a bed of solids and a flow of fluid. As a result, the
solid particles are transformed into a fluid-like behavior that can be used for different
purposes. The fluidized bed reactor is one of the most important technologies for gas-solid
heterogeneous operations chemical or petrochemical, considering catalytic or non catalytic
processes (Kunii and Levenspiel 1991). The most important industrial applications
include catalytic cracking, coal combustion and biomass combustion. One of the most
relevant type of fluidized bed reactor is the ascendant flow reactor, which is also known as
riser. The riser reactors consist of a tubular column in which both solid and gas flow
upwards. The first fluidized bed gas generator was developed in Germany by Fritz Winkler
in the 1920s. Later in the 1930s, the american petroleum industry started developing the
fluidized bed technology for oil feedstock catalytic cracking, becoming the primary
technology for such applications (Tavoulareas 1991).
Inside the riser reactor, solid particles have a wide range of residence time, which is a
disadvantage that reduces the overall conversion and the selectivity of the chemical
reactions. For that reason it has recently grown the interest in a new type of gas-solid
circulating reactor known as downer. In this reactor the gas and the solid flow cocurrently
downward, creating hydrodynamic features comparable to a plug flow reactor and allowing
a better control over the conversion, the selectivity and the catalyst deactivation. The
concept of downer reactor gas-solid appeared in the 1980s, with the first studies on the fluid
dynamics of gas-solid suspensions (Kim and Seader 1983) and with the first downer
reactors for patents developed by Texaco for the FCC process (Gross Benjamin and Ramage
Michael P 1981; Niccum Phillip K and Bunn Jr Dorrance P 1983). In these studies it is
observed that in the downer reactor has a uniform distribution of two-phase flow along the
reactor, also observed that the contact time is very low, achieving a 20% decrease in the
amounts of coke produced during the FCC process.
Fluidization occurs when a gas or liquid is forced to flow vertically through a bed of
particles at such a rate that the buoyed weight of the particles is completely supported by
the drag force imposed by the fluid.
As the superficial gas velocity, U, is increased stepwise beyond the minimum fluidization
velocity, it is observed different types of flow regimes. The principal ones are schematically
shown in Figure 1. The flow regimes are listed by increasing value of U as follows:
Bubbling fluidization
Slug flow
Turbulent fluidization
The behavior of solids fluidized by gases fall into four clearly recognizable groups,
characterized by density difference (ρs – ρf) and mean particle size. The features of the
groups are: powders in group A exhibit dense phase expansion after minimum fluidization
and prior to the commencement of bubbling; those in group B bubble at the minimum
fluidization velocity; those in group C are difficult to fluidize at all and those in group D can
form stable spouted beds (Geldart 1973; González Silva et al. 2010). Desirable properties
of particles and gas for fluidized bed are delineated in Table 1.
For better understanding of these phenomena and to facilitate the solution of mathematical
models is necessary to make an analysis of experimental data. This experimental analysis
requires specialized measurement techniques are able to explain the flow field must also be
automated to minimize human involvement in the process of collecting data.
NMT Ref for more details
PIV measures whole velocity fields by taking (van Buijtenen et al. 2011;
two images shortly after each other and Fu et al. 2011; He et al.
Particle
calculating the distance individual particles 2009; Hernández-Jiménez
Image
travelled within this time. The displacement et al. n.d.; Kashyap and
Velocimetry
of the particle images is measured in the Gidaspow 2011; Laverman
(PIV)
plane of the image and used to determine the et al. 2008; Sathe et al.
displacement of the particles 2010)
Table 2. Non-intrusive measurement techniques.
The measurement techniques, to capture the important fluids dynamic behavior of the two-
phase flow, can be classified as non-intrusive (NMT) and intrusive (IMT) techniques. The
intrusive techniques are generally probes used to study local basic flow phenomena. Some
of these are intended only as research instruments. The most common parameters that are
measured with such probes are solids mass flows, radial and axial solids concentration,
solids velocities, and distribution.
The particles can be deposited in the measuring device reducing its performance or causing
malfunction. Besides this, the flow area reduction makes of the intrusive devices not the best
solution. Non-intrusive techniques to characterize the flow within a fluidized bed are more
desirable because it does not disturb the flow behavior. In the Table 2 and Table 3
classification techniques are included and recent successes have been achieved.
IMT References
Governing equations for Eulerian-Eulerian model are here presented in tensor notation.
t
g g g g vg 0 (1)
s s s s vs 0 (2)
t
Where α, 𝜌 and v are volume fraction, density and the vector velocity, respectively. No
mass transfer is allowed between phases.
t
g g vg g g vg vg gp g g g g vs vg (3)
p and are fluid pressure and gravity acceleration. β is the drag coefficient between the
phases g and s. The stress tensor is given by:
s s vs s s vs vs sGs s s s g vg vs (5)
t
2
s s s vs vs s svs
T
(6)
3
G is the modulus of elasticity given by:
Where αs,max is the maximum solid volume fraction and β is the interface momentum
transfer proposed by Gidaspow, (1994):
s 1 g g s g vs vg
150 1.75 | g 0.8
g d p2
dp
(8)
3 s g g vs vg 2.65
CD g | g 0.8
4 dp
Where dp and CD are the particle diameter and the drag coefficient, based in the relative
Reynolds number (Res)
24 1 0.15Re0.687 | Re 1000
s
CD s (9)
Res
0.44 | Res 1000
g vs vg
Res (10)
g
H g g vg H g g gTg Ts Tg g g Hr Ctr
t g g g
(11)
r
s s Hs s s vs Hs ssTs Tg Ts (12)
t
Where
H = Specific enthalpy
T = Temperature
γ = Interface heat transfer coefficient: Nu / d p
λ = Thermal conductivity
Turbulence is that state of fluid motion which is characterized by random and chaotic three-
dimensional vorticity. When turbulence is present, it usually dominates all other flow
phenomena and results in increased energy dissipation, mixing, heat transfer, and drag. The
physical turbulence models provide the solution the closure problem in solving Navier –
Stokes equations. While there are ten unknown variables (mean pressure, three velocity
components, and six Reynolds stress components), there are only four equations (mass
balance equation and three velocity component momentum balance equations). This
disparity in number between unknowns and equations make a direct solution of any
turbulent flow problem impossible in this formulation. The fundamental problem of
turbulence modeling is to relate the six Reynolds stress components to the mean flow
quantities and their gradients in some physically plausible manner.
The most important numerical methods used to approximate the partial differential
equations by a system of algebraic equations in terms of the variables at some discrete
locations in space and time (called “discretization method”) are the Finite Volume (FV), the
Finite Difference (FD) and the Finite Element (FE) methods. In this book, the finite volume
method and the commercial software CFX® 12.0 were chosen; the solution domain is
discretized in a computational mesh that can be structured or unstructured (G. González et
al. 2012). (G. González et al. 2012)
v S (13)
t
II III IV
I
I. Transient term
The transport equations are integrated in each computational cell using the divergence
theorem over a given time interval ∆t:
t t
t
V t
dV v dA dA S dV dt
v
(14)
Figure 2. Gas flow over a flat solid surface (left to right) experimental picture, refined mesh near the wall and contrast
between experiment and discretization.
Linearization and interpolation techniques can be clarified considering the finite volume P
shown in Figure 3
Aw
dA P W Dw P W (15)
hw
Where S P is the value of source term in the center of the cell P and VP is the volume of
computational cell centered on node P. The method to represent S P was suggested by
Patankar, 1980
*
dS
S P S P P P*
*
P (17)
d
This type of linearization is recommended since the source term decreases with increasing
Φ. The source term coefficients are represented by:
* dS P * *
P VP
SC S P (18)
d
*
dS
SP P VP
(19)
d
4.3.2 Spatial discretization
The most widely used in CFD is first and second order Upwind methods. In the first order
one, quantities at cell faces are determined by assuming that the cell-center values of any
field variable represent a cell-average value and hold throughout the entire cell. The face
value (Φw) are equal to the cell-center value of Φ in the upstream cell.
v dA v A
w w W CwW (20)
In the second order one, quantities at cell faces are computed using a multidimensional
linear reconstruction approach (Jespersen and Barth 1989). In this approach, higher-order
accuracy is achieved at cell faces through a Taylor series expansion of the cell-centered
solution about the cell centroid. Thus, the face value Φw is computed using the following
expression:
3 1 1
w W WW W W WW (22)
2 2 2
The east face coefficient and matrix coefficient are shown below
3 1
e P W (23)
2 2
1
Aw MAX Cw ,0 MAX Ce ,0 Dw (24)
2
4.3.3 Temporal discretization
Temporal discretization involves the integration of every term in the differential equations
over a time step ∆t. A generic expression for the time evolution of a variable Φ is given by
F (25)
t
Where the function F incorporates any spatial discretization. The first-order accurate
temporal discretization is given by
n1 n
F (26)
t
And the second-order discretization is given by
3 n1 4 n n1
F (27)
2t
5 Case Studies
In order to give a better introduction with regards to the simulation of fluidized beds, in this
chapter there are presented three case studies that were carried out by using a CFD software
package. The case studies were carried out using simulations in dynamic state. These
simulations were set up taking into account the average value of the Courant number,
which is recommended to be near 1. Besides this, it was used a constant step time, in this
way was possible to have numerical stability during the execution of each of the
simulations.
Lab scale riser reactor (Samuelsberg and Hjertager 1996; V Mathiesen 2000). Riser height, 1
m; riser diameter, 0.032 m. Experimental data and LES - Smagorinsky simulations were
compared for three velocities with initial particle bed, 5cm.
∆x = 2 mm
The boundary conditions for both cases are shown in Table 5 and Table 6.
Gas velocity = 0.36; 1.42 m/s
In
Particle mass flow equal to the output
Out Opening = atmospheric pressure
Particles = free slip and No slip
Wall
Gas = no slip
Initial height Bed height = 0,05 m
Particles 60 μm; 1600 kg/m3
Table 5. Boundary conditions for the Case 1.
In addition, tests were made with a 500.000 control volume mesh with same block
distribution (the description of volume distribution in the meshes, are presented in Table 7).
Obtaining similar results with the 100.000 control volume mesh. Both meshes are shown in
Figure 4.
Volumes
Mesh dxdp dx
Number
I 15 99900 0.05
II 10 467313 0.08
Table 7. Volume discretization of the meshes.
Numeric calculations performed (Vreman, Geurts, and Kuerten 1997; Chow and Moin 2003)
showed that the required values to obtain an accurate numerical solution, it is necessary to
use a ratio dx 0.25 for the second order spatial scheme, and a ratio dx < 0.5 for the sixth
order scheme. The values of dx presented in Table 7 are within the range recommended in
the literature (Chow and Moin 2003; Agrawal et al. 2001; van Wachem 2000; Ahmed
and Elghobashi 2000; Vreman, Geurts, and Kuerten 1997).
Figure 4. Schematic diagram of the Table 7 meshes. Up: Mesh I. Down: Mesh II
Figure 5 presents the solid volume fraction time evolution for the mesh II with superficial
velocity 1 m/s. At the beginning, the solids present in the riser are forced to flow in the
upward direction, similar to a plug flow. When the bed of solids starts to expand, it is
observed high solid particle concentration at the center of the tube and near the walls
(Figure 5). This reordering of solid particles is a counteraction in order to offer a lower
resistance to the gas flow. This type of flow regime is known as pre-fluidized bed.
It is important to mention that one of most relevant characteristics of the fluidization is the
high contact area between the solid particles and the fluid. In this way, a cubic meter of
particles of 100 micron contains a superficial area of around 30000 m2. The advantage of this
high surface area is reflected in a high mass and heat transfer rates between the solid and
the fluid.
Figure 5. Evolution of the volume fraction field in a fluidized bed at 0, 11, 35, 70, 90, 132, 165,
185, 198, 220, 242, 264, 275, 290, and 317 ms.
Figure 6 shows the similarity between results presented by Miller and Gidaspow (1992).
Here it is represented the regions of high and low solid concentration. Near the walls
velocity is negative and near the center velocity is positive.
The annular-core behavior is something that detrimental in the units of Fluid Catalytic
Cracking (FCC), since big fraction of the oil is converted in a region where the catalyst
works less efficient. In addition to this, the particles that flow at center core are expose to
bigger concentrations of oil compounds, which is something that produces faster
deactivation of the catalyst. One the strategies to solve this issue is to inject pressurized gas
in perpendicular direction to the flow in the reaction zone. Another solution is to include
rings connected to walls, with the purpose of redirecting the solids from the wall towards
the center.
2.5
H= 0.32 m
1.5
Solids Velocity , m/s
0.5
-0.5
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
r/R
Figure 6. Comparison of solid phase velocity profile presented by Miller and Gidaspow
(1992) with the CFD simulations (-▲-) and experimental data performed by Samuelsberg
and B. H. Hjertager (1996) (●).
To get an impression regarding the flow behavior inside the column, the time averaged
solid volume fraction is plotted at different column heights, 0.16 m, 0.32 m and 0.48 m
(Figure 7). Here it can be observed the strong tendency of the solid particles to be near the
wall.
Figure 7. Axial profile of the solid phase volume fraction fields in the center (left) and radial
profiles at 0.48 m, 0.32 m, 0.16 m (right up to down). Superficial velocity 0.36 m s-1
5.2 Case 3
Pilot plant scale riser reactor (Bader, R., Findlay, J. and Knowlton, TM 1988). Riser height:
13 m, riser diameter 0.3 m. Entrance with angle 60°, gas superficial velocity 3.7 m and solids
flux 98 kg/(s.m^2) as shown in Figure 8.
Figure 8. Solids volumetric fraction in the center of the riser. Simulation time 15 sec. Left to
right: LES Smagorinsky, LES WALE, LES Dynamic model, Detached Eddy Simulation
(DES).
In the
Figure 8 can be observed that the solid particles enter to the reactor uniformly distributed,
after a short distance these particles start falling due to the gravity and they start flowing
over the wall of the inclined pipe. After this, the solids fall into a turbulent zone where they
get mixed. Some of the particles will continue falling over the vertical wall opposite to the
entrance. The core-annular zone is formed at some height in the middle of the column.
6 Conclusions
Computational fluids dynamics is a very powerful tool understanding the behavior of multi
phase in engineering applications.
Large eddy simulation (LES) turbulence method provides a very detailed description of two
phase flow, which makes it suitable for simulation models that are validated with
experimental data. By applying the LES method, it is possible to characterize different
regions of a fluidized bed (core-annulus). LES can be considered as a valuable method for
development and validation of closure models that include additional phenomena like heat
exchange, mass transfer and chemical reactions.
It is important to constantly monitor the simulation, using parameters such as the Courant
number, creating a function that calculates the maximum and average number of the control
volume courant. The average value is recommended that is near or less than unity.
7 Acknowledgements
The author G. Gonzalez is grateful to PETROBRAS and the National Council for Scientific
and Technological Development (CNPq) for the financial support to this research.
8 References