Accepted Manuscript Title: Computational Fluid Dynamic Modelling of FCC Riser: A review Author: Milinkumar T. Shah Ranjeet P. Utikar Vishnu K. Pareek Geoffrey M. Evans Jyeshtharaj B. Joshi PII: DOI: Reference:
S0263-8762(16)30077-6 http://dx.doi.org/doi:10.1016/j.cherd.2016.04.017 CHERD 2263
To appear in: Received date: Revised date: Accepted date:
26-9-2015 19-4-2016 25-4-2016
Please cite this article as: Shah, M.T., Utikar, R.P., Pareek, V.K., Evans, G.M., Joshi, J.B.,Computational Fluid Dynamic Modelling of FCC Riser: A review, Chemical Engineering Research and Design (2016), http://dx.doi.org/10.1016/j.cherd.2016.04.017 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
Computational Fluid Dynamic Modelling of FCC Riser: A review
2
Milinkumar T Shah1, Ranjeet P Utikar1, Vishnu K Pareek1, Geoffrey M Evans2 and Jyeshtharaj B
3
Joshi3 1
ip t
Department of Chemical Engineering, Curtin University, Western Australia, Australia Department of Chemical Engineering, University of Newcastle, New South Wales, Australia 3 Homi Bhabha National Institute, India *Corresponding author:
[email protected]
2
cr
4 5 6 7
Highlights
9
(1) CFD models and experimental data of FCC riser are critically reviewed.
us
8
(2) Governing equations and constitutive correlations of CFD models are explained.
11
(3) Effect of important closure models on predictions is analysed.
12
(4) Shortcomings of CFD models are identified and suggestions for future work are made.
13
Abstract
14
Design and scale-up of fluid catalytic cracking (FCC) riser is still largely empirical, owing to limited
15
understanding of inherent multiphase flow in this equipment. The multiphase flow of FCC riser has
16
therefore been extensively investigated both experimentally and computationally. The experiments
17
have provided significant insight into gas-solid flow patterns inside cold-flow risers, but simultaneous
18
observations on flow and performance parameters (conversion and yields) in FCC riser are rarely
19
found in literature. Consequently, computational fluid dynamic (CFD) models of FCC riser that can
20
simultaneously account for flow, interphase interactions, droplet vaporization and cracking kinetics
21
have been developed. The CFD modelling of FCC riser, despite several efforts, has still remained a
22
challenge as it requires careful consideration of governing equations and closure models. This review
23
presents state-of-the-art in CFD modelling and experimental analysis of gas-solid hydrodynamics and
24
reactive flow of FCC riser. The CFD models are explained in greater detail with governing equations,
25
constitutive relations, and physical significance of all the terms. A brief review of DNS studies on
26
cluster formation, gas-solid drag, and turbulent interactions is also presented. Impact of important
27
closure models such as drag models, viscous stress models, boundary conditions, droplet vaporization
28
models, and kinetic models on predictions is critically examined. The review identifies major
29
shortcomings of current CFD models and makes detailed recommendations for future work.
Ac ce
pt
ed
M
an
10
Page 1 of 1
Page 1 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Keywords: Fluid catalytic cracking, riser, hydrodynamics, computation, simulation
2
Content
3
1. Introduction
4
2. Gas-solid flow in FCC riser 2.1. “Core-annulus” radial profile
6
2.2. “S-shape” axial profile
7
2.3. Clusters
8
2.4. Residence time distribution (RTD)
9
2.5. Entry and exit effects
3.1. Eulerian-Eulerian model 3.1.1.Governing equations
13
3.1.2.Constitutive equations
Ac ce
18
3.2.4.Effect of boundary conditions
16
us
3.2. Effect of closure models on flow predictions
17
15
an
pt
12
14
M
11
3. CFD models of gas-solid flow in riser
ed
10
cr
5
ip t
1
3.2.1.Effect of drag models 3.2.2.Effect of gas phase stress models 3.2.3.Effect of KTGF models
19
3.3. Shortcomings
20
3.4. Eulerian-Lagrangian model
Page 2 of 2
Page 2 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
3.5. Drift flux analysis
2
4. Direct numerical simulations
3
5. CFD models of a reactive flow in FCC riser 5.1. FCC riser performance
5
5.2. Coupling of gas-solid flow model with droplet flow and vaporization models
6
5.3. Coupling of gas-solid flow model and cracking kinetic model
7
5.4. Predictions
8
5.5. Shortcomings
cr
us an
9
ip t
4
6. Recommendation for future work
Ac ce
pt
ed
M
10
Page 3 of 3
Page 3 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1. Introduction
2
In petroleum refinery, an FCC unit converts low-value heavy residuals into valuable light products such as
3
gas (C1-C4), gasoline (C5-C12), and diesel (C10-C15). In an FCC, a riser (Figure-1) is a long vertical pipe
4
acting as a reactor where cracking reactions take place. In the FCC riser, steam, as an inert fluidizing
5
medium, enters from the bottom and fluidizes regenerated catalyst. A few meters above the catalyst and
6
steam inlets, hydrocarbon feedstock (vacuum gas oil - VGO, residue of vacuum distillation column) enters as
7
a liquid spray through a set of atomizing nozzles. The fine droplets of the feedstock vaporize due to rapid
8
mixing and heat transfer, and endothermic reactions then take place in the gas phase. The catalyst, droplets,
9
and gas (hydrocarbon vapor and steam) concurrently flow upward along with reactions and droplet
10
vaporization. This multiphase flow governs distribution of phases, temperature and concentrations, and in
11
turn it dictates overall performance of the FCC riser. Consequently, the multiphase flow of FCC riser has
12
been investigated both experimentally and computationally.
us
cr
ip t
1
an
Dimension and flow conditions (Buchanan, 1994; Berruti et al., 1995; Grace et al., 1997)
ed
M
Dimensions Height Diameter Operating conditions Inlet temperature of feedstock
30-40 m 1-2 m
150-300 ⁰C
Ac ce
pt
Inlet temperature of catalyst
Riser top temperature 500-550 ⁰C
Hydrocarbon Feedstock FCC Catalyst
`
Steam
675-750 ⁰C
Solids circulation rate Catalyst to oil ratio Dispersion steam Pressure Solids residence time Catalyst particles Average size Particle density Bulk density Typical Umf Typical VT Geldart classification Feed droplets Droplet diameter
>250 kg/m2s 4-10 wt% 0-5 wt% 150-300 KPa 3-15 s 70 μm 1200-1700 kg/m3 750-1000kg/m3 0.001 m/s 0.1 m/s A up to 2000 μm
Page 4 of 4 Page 4 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review Number of feed nozzles
4-12
Figure-1: Schematic of FCC riser with dimensions and flow conditions Several experiments have been conducted to investigate gas-solid flow in riser. These experiments have
2
shown (i) large variations in the volume fraction of the solids phase in different sections of riser, (ii) particle-
3
scale heterogeneity due to formation of particle aggregates called 'clusters', (iii) considerable back mixing of
4
the solids phase, and (iv) significant effects of entry and exit boundaries. While the experimental
5
observations have provided significant understanding of the gas-solid flow patterns inside riser, they do not
6
provide detailed information on hydrodynamics. In addition, experiments that provide observations of both
7
flow and performance parameters (conversion and yields) in FCC riser are rarely found in literature. As a
8
result, CFD models of riser have been developed. The CFD modelling of the riser flow involves careful
9
selection and/or customization of the governing equations and closure models. Mostly two-fluid Eulerian-
10
Eulerian gas-solid flow models have been used to model the gas-solid flow, while the two-fluid models
11
coupled with cracking kinetics models have been used to model reactive flow in FCC riser. The two-fluid
12
model is often extended by introducing feed droplet phase, which is modelled as either Eulerian or discrete
13
Lagrangian phase. The model with the third (droplet) phase is coupled with appropriate droplet vaporization
14
model. Predictions from such models are largely dependent on closure models. Multiple models are proposed
15
in literature for solids phase properties, gas-solid drag force, viscous stresses, boundary conditions, droplet
16
vaporization, cracking kinetics etc. In last few years, new closure models have been proposed based on
17
multiscale approaches, cluster-based hypothesis, and direct numerical simulations. The closure models have
18
been investigated in previous studies, who have recommended different closure models after comparing their
19
predictions with experimental data. A systematic review of previous studies is necessary to identify most
20
suitable set of closure models.
21
Only a few reviews on CFD modelling of risers are available in literature. Berutti et al. (1995) has presented
22
a detail review of flow patterns and models of gas-solid flow of the riser. Berruti et al. (1995) have
23
thoroughly explained models that predict only axial variation of solids suspension density, and those that
24
predict only radial variation of solids suspension density. Berruti et al. (1995) have briefly summarized
25
previous works on CFD models of the gas-solid flow in riser. Godfrey et al. (1999) have evaluated the
26
models that predict either axial or radial profiles by comparing predictions with experimental data. Wang
27
(2009) has reviewed CFD models of a gas-solid flow focusing mainly on drag models available to capture
28
sub-grid scale flow structures such as clusters. Wang et al. (2010) have critically reviewed applicability of a
29
drag model based on the energy minimization multiscale (EMMS) approach in CFD simulations of riser.
30
Gupta et al. (2010) have reviewed modelling of FCC riser by covering literature on models for cracking
31
kinetics, droplet vaporization and hydrodynamics. The hydrodynamic modelling section of Gupta et al.
32
(2010) doesn’t focus exclusively on CFD models. In last one decade, a large number of CFD studies of the
33
riser flow have been published with them being mostly focused on evaluation and modification of different
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 5 of 5 Page 5 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
closure models. During this period, several studies have also been conducted using coupled CFD models to
2
investigate alternate design and flow conditions in FCC riser. A comprehensive review of CFD models of
3
FCC riser, particularly with focus on recent findings, is not available in literature.
4
The current study presents a critical review of CFD models of FCC risers accompanying with analysis of
5
experimental data. The review divides the available literature on this subject in two groups, those on gas-
6
solid flow and those on reactive flow in FCC riser. The CFD models and experimental data published under
7
both groups are comprehensively analysed. The CFD models are explained in detail with governing
8
equations, constitutive relations, and physical significance of all the terms. A brief review on direct
9
numerical simulation (DNS) studies on cluster formation, gas-solid drag, and turbulent interactions is also
10
included. The impact of important closure models such as drag models, viscous stress models, boundary
11
conditions, droplet vaporization models, and kinetic models on predictions is critically analysed. More
12
specifically, conventional drag models vs. multiscale drag models, laminar vs. turbulent models for the gas
13
phase, 2D vs. 3D boundary conditions, and no-slip vs. partial-slip wall boundary conditions are analysed for
14
predictions of gas-solid flow in riser. Furthermore, uncertainty associated with the selection of droplet
15
vaporization models and cracking kinetic models is also critically examined. Based on the thorough analysis,
16
major shortcomings of current CFD models are explicitly brought out and detailed suggestions for future
17
work are given.
18
2. Gas-solid flow in FCC riser
19
The gas-solid flow in FCC riser has been extensively investigated by cold-flow pilot-scale CFB setups. In
20
these experiments, controlled quantities of the solids and gas (typically air) are circulated through a riser. As
21
summarised in Table-1, these experiments are conducted for different riser dimensions, particle properties,
22
and flow conditions (such as mass flux and fluidizing gas velocity). Measurements of the pressure, velocity,
23
and volume fraction are taken at multiple locations along height and radial positions by using various
24
invasive (pressure taps, optical probes and capacitance probe) or non-invasive (high speed videography,
25
tomography, particle image velocimetry, electrical capacitance volume tomography, laser Doppler
26
velocimetry, radioactive particle tracking and positron emission particle tracking) techniques. Though the
27
operating conditions in these experiments do not resemble those in the actual FCC riser, the cold-flow
28
experiments can still provide useful observations of inherent gas-solid hydrodynamics in FCC riser. Some of
29
the important flow characteristics can be summarised as:
30
Table-1
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 6 of 6 Page 6 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Table-1: Cold-flow FCC riser experiments
Ac ce
pt
cr
us
an M
ed
Yerushalmi et al. (1976) Rhodes and Geldart (1986) Ishii et al. (1989) Bader et al. (1988) Li et al. (1988) Pressure taps, Optical probes Ambler et al. (1990) Radioactive tracer detection Louge and Chang (1990) Pressure taps, Capacitance probes Bai et al. (1992) Pressure taps Yang et al. (1992) LDA Miller and Gidaspow (1992) X-ray densitometer, Extraction probe Brereton and Grace (1993) Needle capacitance probe Patience and Chaouki(1993) Radioactive tracer detection Horio and Kuroki (1994) High speed camera Wei et al. (1995) Optical fibre sensors Knowlton et al. (1995) Pressure taps Pugsley and Berruti (1996) Pressure taps Nieuwland et al. (1996b) Optical probes, Pressure taps Pugsley et al. (1997) Pressure taps Wei et al. (1998) LDA Cheng et al. (1998) LDA Noymer and Glicksman (1998) TIV Issangya et al. (2000) Pressure taps Bai et al. (1999) Optical fibre probes Mathiesen et al. (2000) LDA/PDA Smolders and Baeyens (2000) Electric conductivity of tracer Sharma et al. (2000) Pressure taps, Needle capacitance probe Pärssinen and Zhu (2001) Fibre optic probes Ibsen et al. (2002) LDV/PIV Harris et al. (2003) Light sensitive photomultiplier tube Zhang et al. (2003) Optical fibre probe, LDV Pandey et al. (2004) LDV Bhusarapu et al. (2005) CARPT and CT Mabrouk et al. (2007) RPT He et al. (2006) PIV Van de Velden et al. (2007) PEPT Kim et al. (2008) Pressure taps Yan et al. (2009) Phosphor tracer technique Chang et al. (2010) PEPT Wang et al. (2010) Pressure sensors Xu and Zhu (2011) High speed camera, Optical fibre probe Heynderickx et al. (2011) LDA Gao et al. (2012) Optical probes Pantzali et al. (2013) LDA Zhang et al. (2013) Pressure sensors Rodrigues et al. (2015) Pressure taps Monazam et al. (2016) Pressure taps Experimental observations: (a) Axial profiles of volume fraction of solids phase (b) Axial profiles of pressure drops (c) Axial profiles of velocity of solids phase (d) Radial profiles of volume fraction of solids phase (e) Radial profiles of velocity of solids phase (f) Observations on clusters (g) RTD (h) Effect of entries or exit configurations
Experimental setup and operating conditions Dt (m) H (m) dp (μm) Gs (kg/m2s) ug (m/s) 0.0762 7.3 60 52.90 3.65 0.152 6 64 8.5-107 2.5-4.5 0.05 2.79 60 10.7 1.29 0.305 12.2 76 147 9.1 0.09 10.5 54 14.3 1.52 0.05 3 106 124.305 4.5-7.1 0.197 7 61 17 2 0.140; 0.186 8 54, 280, 165, 31 30-180 2-8 0.140 11 54 26.56-119.43 1.5-6.5 0.075 6.58 75 12 2.89 0.152 9.3 148 42, 48, 62 6.5 0.083 5 277 20-140 4-8 0.2 1.6 61.3 0.016-0.6 0.15-0.6 0.2 1.6 61.3 0.016-0.6 0.15-0.6 0.20 14.2 76 485 5.2 0.05 5 208 51.3-700 8.5 0.054 8 129 100-300 10 0.1 and 0.2 6 and 12 220, 230, 71, 80 10-45; 10-85 4-6 0.186 8 54 200 2.3-6.2 0.186 8 54 200 2.3-6.2 0.159 (square) 2.44 69 25, 18 3, 6 0.762 6.1 70 18-325 4-8 0.0762; 0.102 6.1, 8.2 70 425 8 0.032 1 120, 185 0.8, 1, 1.2 0.1 6.5 90 5.19-34.1 2.82-4.92 0.15 11 70, 120 75 4, 4.9, 6 0.076 10 67 100, 300 8 0.02 2 250 0.45 0.14 (square) 5.8 25 1.3-4.5, 2-25.5 0.418 18 77 19-180 1.8-8 0.305 15.2 812 3.4-17.1 3.75-5.4 0.152 7.9 150 33.7, 36.9 3.9 and 4.5 0.52 1 250, 170 2, 12 2D, 0.05 x 1.5 1 335 5,10, 20 2-2.3 0.10 6.5 90 25, 89.7 3.6 and 9 0.05 4.5 70 7-300 0.2-9 0.186 10 78 40.8, 229.4 3.156, 5.989 0.16 7.9 120 5, 622 1, 10 0.27 (square) 10 330 113-165 15.5 0.019 x 0.114 7.6 67-288 50-200 3-8 0.1 8.765 77 3.0 2.48–7.43 0.095 and 0.14 2.2 139 0.94-1.25 0.1 8.7 77 1 3.5 and 5.3 0.316×0.08, 0.1×0.1 4.5, 4.3 350,175 12.7, 13.9, 16.1 4.4, 5.3, 5 0.3 18 300,250, 280 60 8-10 0.3 16.65 60, 230, 812 6.89-281 5.35-7.7 Abbreviations: LDA = Laser Doppler Anemometry PDA = Phase Doppler Anemometry PIV = Particle Image Velocimetry CARPT = Computer Aided Radioactive Particle Tracking CT = Computed Tomography RPT = Radioactive Particle Tracking TIV = Thermal Image Velocimetry PEPT = Positron Emission Particle Tracking
ip t
Measurement technique Pressure taps, High-speed camera Pressure taps Optical fibre probes
ρs (kg/m3) 510.98 1800 1000 1714 1654 Sand 1545, 706, 750, 2660 1545 930 2650 2630 1780 1780 1712 2580 2540 2500, 2200, 2600, 1500 1398 1398 6970 1600 1600 2400 2200 1500 2400 4060 1398 189 2550 2500; 3400 2500 1740 1225 2260 2630 1877-2498 1550 2400 1550 2600 2650, 2600, 3300 2550, 1250, 189
Observations (f) (a) (a), (b), (f) (a), (d) (a), (e), (d) (g) (a), (b) (a) (d) € (b), (d), (e) (a) (a), (e), (g) (a), (f) (f) (b), (d), (e) (b) (d), (e) (b), (h) (a), (d), (e) (d), (h) (f) (a) (d) (d), (e) (g) (b), (e), (f) (a), (d), (e) (a), (b), (d) (g) (d), (e), (g) (f), (e) (d), (e) (a), (g) (e) (g) (a), (h) (e), (g) (g) (e) (f) (h) (a), (c), (d) (e) (h) (b) (a), (b), (h)
31 Page 7 of 7 Page 7 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
2.1. “Core-annulus” radial profile
2
Bader et al. (1988), Li et al. (1988), Miller and Gidaspow (1992), Yang et al. (1992), Nieuwland et al.
3
(1996b), Knowlton et al. (1995), Mathiesen et al. (2000), Bhusarapu et al. (2005), and Wang et al. (2010)
4
have reported radial profiles of the volume fraction and velocities of the solids phase. A lower volume
5
fraction and higher velocity of the solids phase are observed at centre of riser (Figure-2). This region is
6
referred as a “core” region, where the solids are carried upward due to flow of the gas. Contrary to the core
7
region, a lower velocity or downward flow of the solids with a higher volume fraction of the solids are
8
observed near wall. The region near wall is defined as an “annular” region. The heterogeneous radial
9
distribution of velocity and volume fraction is widely known as the “core-annular” flow pattern. The higher
10
solids volume fraction in the annular region is largely attributed to friction between flow and wall and energy
11
dissipation due to particle-particle and particle-wall collisions. The solid segregation is also caused by
12
formation of clusters (Noymer and Glicksman, 2000 and Pandey et al., 2004). Bader et al., (1988), Yang et
13
al. (1992), Miller and Gidaspow (1992), Noymer and Glicksman (2000) and Pandey et al. (2004) have
14
observed downward sliding motion of segregated particles near wall. Yang et al. (1992) have observed
15
concurrent movement of phases in the core region, while they have observed counter current movement of
16
the solids and gas near wall. The counter current movement of the phases lead to higher slip velocity near
17
wall than that in core region (Yang et al., 1992). Experiments for various flow conditions and riser systems
18
have revealed that the core region occupies almost 60-80% of cross-section with the annular region
19
occupying rest 20-40%. Miller and Gidaspow (1992) have observed (Figure-2a and b) that an increase in the
20
solids flux at a constant gas velocity results in higher solids volume fraction in both the core and annular
21
regions; while an increase in the solids flux had a minor impact on velocity of the solids in the core region.
22
Nieuwland et al. (1996b) have observed (Figure-2c and d) that an increase in velocity of the gas phase at a
23
constant solids flux results in a lower solid volume fraction in both the core and annular regions. Perssinen
24
and Zhu (2001) have reported the solids volume fraction at different heights, and observe that the values in
25
both the core and annular regions decrease with an increase in the height.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 8 of 8 Page 8 of 114
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
us
(a)
pt
ed
M
an
(b)
Ac ce
(c) (d) Figure-2: Core-annulus radial profiles, (a) and (b) effect of solids flux, (b) and (c) effect of gas velocity. 1
2.2. “S-shape” axial profile
2
Li et al. (1988), Louge and Chang (1990), Pärssinen and Zhu (2001), Rhodes and Geldart (1986), Wei et al.
3
(1998), Issangya et al. (2000) and Zhang et al. (2003) have reported axial variation in the solids volume
4
fraction (Figure-3a). The axial profiles have shown an “S-shape” (Figure-3a and b) with three distinct
5
regions (i) a dense bottom region, (ii) a dilute top section, and (iii) middle section presenting a transition
6
from the dense to dilute flow. Experimental under various flow conditions in different riser systems have
7
showed several similarities, i.e. the volume fraction of the solids in the dense bottom section is in a range of
8
0.3-0.15, and that in the dilute top section is in a range of 0-0.05. The experimental data shows significant
9
variations in the height of the transition from the dense to dilute flow (20-60% of the total height). Pärssinen
10
and Zhu (2001) have found significant impact of the inlet gas velocity at a constant solids flux on axial
Page 9 of 9 Page 9 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
profile the solid volume fraction (Figure-3c). A lower gas velocity gives a clear distinction between the
2
dense and dilute region; while a higher gas velocity results in a smooth transition. Rhodes and Geldart (1986)
3
and Issangya et al. (2000) have investigated the effect of the solids flux on axial profiles (Figure-3b). Rhodes
4
and Geldart (1986) have found that a higher solids flux gives dense flow with higher values of the solids
5
volume fraction than those from lower solids flux condition. Issangya et al. (2000) have found that a lower
6
solids flux give a longer developing section than that resulted at a higher solids flux. Li et al. (1988) have
7
investigated the effect of the solids inventory in CFB on axial profiles (Figure-3d). Notably, in Li et al.’s
8
experiment, there is no mechanical valve that controls quantity of the solids circulated through the riser. Li et
9
al. (1988) have found that a higher solids inventory in the CFB loop gave a longer dense phase and transition
(a)
pt
ed
M
an
us
cr
regions than those from a lower solids inventory.
Ac ce
10
ip t
1
(b)
Page 10 of 10 Page 10 of 114
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
an
(d) (c) Figure-3: (a) “S” shape axial profiles of solids volume fraction, dimensionless height Vs. solid volume fraction; effect of (b) solids flux, (c) gas velocity, and (d) solid inventory on axial profiles of the solids volume fraction Miller and Gidaspow (1992); Knowlton et al. (1995); Sharma et al. (2000) and Ibsen et al. (2002) have
2
reported variations in pressure drop along height, where it is observed that the bottom section of riser has
3
higher pressure drop with lower values at the top (Figure-4). Miller and Gidaspow (1992) have found that the
4
pressure drop in riser increases with increase in the solids flux (Figure-4a); while Sharma et al. (2000) have
5
found that an increase in the gas velocity reduces the pressure drop in riser (Figure-4b).
Ac ce
pt
ed
M
1
(a) (b) Figure-4: (a) Effect of circulating flux and (b) Effect of velocity of the gas phase on pressure drop profile Page 11 of 11 Page 11 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
2.3. Clusters
2
Yerushalmi et al. (1976) have observed segregation of particles and movement of the segregated solids
3
within a dilute phase flow. They have characterized shape of the segregated solids as “strands”, “steamer”
4
and “ribbons”. Later, Ishii et al. (1989) and Horio and Kuroki (1994) have reconfirmed presence of the
5
particle segregations and defined them as “clusters”. The authors have observed that the clusters are present
6
not only near wall but also in dilute phase at centre of riser. Ishii et al. (1989) and Horio and Kuroki (1994)
7
have also observed that the shape of the clusters continuously changing with flow conditions. Soong et al.
8
(1994) have measured instantaneous concentrations of the solids at various locations in riser, and give three
9
criteria to identify the clusters, (1) solids concentration inside cluster must be significantly higher than local
10
time–mean solids concentration; (2) perturbation caused by cluster must exist for significantly longer time
11
than sampling time interval; (3) this perturbation must also be sensed by a sampling volume which has a
12
characteristic length scale greater than one to two orders of particle diameter. Based on Soong et al.’s
13
criteria, Noymer and Glicksman (1998); Sharma et al. (2000); Manyele et al. (2002); Liu et al. (2005); and
14
He et al. (2006) have recorded measurements of cluster properties. More recently, Shaffer et al. (2010; 2013)
15
have used high speed PIV to capture clusters and their movements along wall of riser. They have found large
16
variations in cluster structures from several particles to several riser diameters. Furthermore, Shaffer et al.’s
17
studies have also revealed that a presence of stream of the gas phase influences make-up and break-up of
18
clusters near wall. Formation of clusters has two profound effects, (i) apparent solids diameter increases
19
which leads to higher terminal velocity and (ii) gas tends to flow around clusters without penetrating inside
20
the clusters (Li and Kwauk, 1994). Consequently, formation of clusters significantly impacts gas-solid drag
21
in riser (discussed in section-3.2.1).
22
2.4. Residence time distribution (RTD)
23
Ambler et al. (1990), Patience and Chaouki (1993), Smolders and Baeyens (2000), Harris et al. (2003),
24
Mabrouk et al. (2007), Van de Velden et al. (2007) and Chan et al. (2010), have reported RTDs from
25
stimulus response experiments at various flow conditions. In the stimulus response experiments, a tracer is
26
injected at bottom of riser and its concentration at outlet is measured at different times after injection. Harris
27
et al. (2003) have measured RTD in a square riser by using a small portion of the solids particles as tracers,
28
and investigate the effect of three types of exit geometries, namely smooth, abrupt and highly refluxing, on
29
RTDs at various flow conditions. Figure-5(a) shows variations in the RTDs caused by variations in the gas
30
phase velocity for the smooth exit riser. Harris et al. (2003) have observed that an increase in the superficial
31
gas velocity decreases both mean residence time and variance of distribution over a range of investigated
32
solids fluxes. At a higher gas velocity, the RTD curves are narrow, relatively unskewed and having a
33
significant long tail. In contrast, at a lower gas velocity, the curves are broad with a long tail. Van de Valden
34
et al. (2007) have also reported RTDs for two operating gas velocities at a fixed solids flux (Figure-5b),
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 12 of 12 Page 12 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
where the RTD from a higher gas velocity is closer to that of ideal plug flow reactor; while the RTD from a
2
lower gas velocity is closer to that of the ideal mixed flow reactor.
us
cr
ip t
1
an
(a) Figure-5: Effect of gas velocity on residence time distribution (a) E(t) Vs. time and (b) F(t) Vs. t/ τ
(b)
M
; τ = half the injection time, approximately 1 s (Van de Valden et al., 2007) 2.5. Entry and exit effects
4
Bai et al. (1992), Brereton and Grace (1994), Zheng and Zhang (1994), Pugsley et al. (1997), Cheng et al.
5
(1998), Kim et al. (2008) and Mabrouk et al. (2008) have investigated the effect of exit configuration on riser
6
hydrodynamics. Cheng et al. (1998) and Kim et al. (2008) have also studied the effect on inlet
7
configurations. Figure-6(a) shows the effect of two types of exit configuration (abrupt exit – L or T-type exit,
8
and smooth exit – long radius bend type exit) on the axial profiles of pressure drop as observed by Pugsley et
9
al. (1997). It can be seen that the effect of the abrupt exit causes higher pressure drop in upper section of
10
riser, while the pressure drop in the bottom section with both types of the exit configurations are the same.
11
Pugsley et al. (1997) have also observed that the effect of the abrupt exit is only restricted to the top section,
12
and this effect is even higher for a lower solid circulated mass flux. Cheng et al. (1998) have also compared
13
axial profiles of the gas volume fraction from abrupt and smooth exits. Cheng et al. (1998) have found that
14
the abrupt exit gives a lower gas volume fraction near the exit, while values in the bottom section from the
15
two types of exits are in reasonable agreement. Cheng et al. (1998) have studied the effect of three types of
16
air entries in an internal nozzle type air inlet arrangement; While Kim et al. (2008) have compared three
17
different types (L-valve, J-valve and Loop-seal type) solid entries. Figure-6(b) shows the effect these three
18
types of solid entries on solid circulation rate vs. solid inventory plots (Kim et al., 2008). An increase in the
19
solid inventory increases the solid circulation rate in riser. At a higher value of the solid inventory, the values
20
of circulating rates from the L-valve and J-valve are in reasonable agreement. For an entire range of
Ac ce
pt
ed
3
Page 13 of 13 Page 13 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
experimental values of the solid inventory (12-36 kg), the loop-seal type solid entry results in significantly
2
lower solid circulation rate than that in the other two types of entries.
(a)
an
us
cr
ip t
1
M
(b) Figure-6: (a) Effect of exit configuration, dimensionless height Vs. pressure drop per unit length (Gs = 85 kg/m2s and ug = 4.2 m/s); (b) Effect of solid phase entry, solid circulation rate Vs. solid inventory (ug = 7 m/s) Experimental observations show large flow structures such as radial and axial distributions extended across
4
length and width of riser. These structures are directly influenced by small scale flow phenomena such as
5
formation of clusters, particle-wall collisions and boundary effects. Capturing these multiple flow structures
6
that have disparate time and length scales in CFD predictions is a challenging task. Higher spatial resolution
7
using fine grid discretization of flow domain and lower time step for transient simulations are essential
8
simulation parameters for riser simulations. However, these parameters are generally constrained by
9
available computational resources. Hence careful selection of modelling approach and governing equations is
10
necessary. Each modelling approach needs closure models such as models for interphase interactions,
11
turbulent forces and boundary conditions. These closure models have been derived empirically or by
12
conducting high resolution sub-grid simulations or DNS. For each closure model, several models are
13
available in literature. Realistic prediction of riser flow critically depends on the selection of closure models.
14
For example, various drag models have been proposed to cater different flow conditions. But for industrial
15
scale riser simulation, a drag model that can capture clustering phenomena even with coarse grid spatial
16
discretization is essential. Similarly, proper inlet and outlet configurations of flow domain are required to
17
capture entry and exit effect. Moreover, capturing dissipation due to particle-wall collisions is critical in
18
configuring wall boundary conditions. Various modelling approaches, model equations and closure models
19
are discussed in detail in the next section.
Ac ce
pt
ed
3
Page 14 of 14 Page 14 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
3. CFD models of gas-solid flow in riser
2
Different types of models of the gas-solid flow in riser have been proposed over last three decades. Harris
3
and Davidson (1994) have classified the models published before early 90’s into three categories, i.e. (i)
4
those that predict axial variation of suspension density, but not radial variation; (ii) those that predict radial
5
variation, but not axial variation; (iii) those that employ fundamental fluid dynamic equations for two phase
6
gas-solid flow. Among these three types of models, the category-(iii) models are most useful to investigate
7
the effect of local flow structures and design configurations (Harris and Davidson, 1994; and Berruti et al.,
8
1995). The models based on fundamental fluid dynamic equations can be further classified in two broad
9
categories, namely the Eulerian-Eulerian (E-E) or the Eulerian-Langragian (E-L) models. The E-E model
10
considers both the gas and solids phases to be continuous and fully interpenetrating. In the E-E model,
11
averaged Navier-Stokes equations of the mass, momentum and energy conservations are solved for both the
12
phases. Closure laws based on the kinetic theory of granular flow (KTGF) (Jenkins and Savage, 1983; Lun
13
and Savage, 1984; Lun et al., 1984) are used to calculate stresses in the solids phase. Coupling between the
14
momentum balance equations is modelled by using correlations for interphase forces; while coupling
15
between the energy balance equations is modelled by using correlations for heat transfer between the two
16
phases. The E-L model represents solids phase as a collection of discrete particles. The motion of each
17
individual particle is solved by Newton’s law, accounting for the effect of particle collisions and forces
18
acting on the particle by the gas. Solution of the force balance equations gives updated velocity and position
19
of each particle. In a further advanced E-L model, discrete element model (DEM), contact force due to
20
collisions between particles are also captured at a much shorter timescale by means of the empirical
21
coefficient of restitution and friction (hard sphere approach) or empirical spring stiffness and a friction
22
coefficient (soft-sphere approach) (Cundall and Strack, 1979; Tsuji et al., 1993; Deen et al., 2007). The
23
interactions between two phases are still modelled using interphase drag and heat transfer correlations
24
similar to those used in the E-E model. Use of the E-L approach requires tracking of billions of solid
25
particles by solving the force equation for each particle and then calculating contact forces between nearby
26
pairs of the solid particles. Hence, the E-L approach becomes computationally impractical for simulations of
27
industrial- or laboratory-scale risers. On the other hand, the E-E approach is computationally less expensive,
28
and hence, it has been widely used for the gas-solid flows in riser. Many of important simulation studies
29
using the E-E gas-solid models for riser simulations are summarised in Table-2.
30
Table-2
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 15 of 15 Page 15 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Drag: (a) (b) (c) (d) (e) (f) (g) (h)
Wen and Yu (1966) Gidaspow et al. (1994) Syamlal – O’Brien (1987) EMMS (Li and Kwauk, 1994) SGS – Filtered drag Richardson and Zaki (1954) Arastoopour et al. (1990) Other
(b), (d) (c) (b), (e) (c) (a) (b), (d) (b), (c), (g) (b), (e) (d) (b), (d), (e) (b), (d) (a) (d) (d) (b)
NC NC (f) NC (e) NC (e)
(a) (a) (a) (a) (a) (a)
(c) (e) (c) (e) (c) (c)
(f) NC NC
(a) (a) (a)
(b) (b) NC NC (b)
(d) (b), (d) (a), (b), (c), (h)
(b),(c), (d) (a)
KTGF closures Collision Shear dissipation rate viscosity
Radial distribution function
Granular energy exchange
Investigations
Boundary condition Solids Walls (for inlets solids) (f) (d) (b) (c)
(b) (c) (a) (b) (a) (a) (a) (b)
(d) (j) (e) (c) (d) (j) (c) (c) (c) (f)
(b) (b) (d) (c) (b) (b) (b) (b) (c) (b)
(a) (a) (a) (a) (a) ((b)
(b) ((c) (e) (e) (d) (b) (d)
(b) (a) (b) (a) (a) (a) (b) (b) (b) ((c)
(a) (b) (a) (b) (a) (b)
(k) (c) (f) (c) (f),(g),(h) (c) (c)+(i)
(e) (b) (b) (b) (b) (b)
(b) (b),(c), (d) (a) (a)
(b) (d) (d) (d) (d) (b) (b)
(b) (b) (a), (c) (b) (b) (b)
(e) (c) (e)
(b) (a) (b)
(f) (c) (f)+(i)
(b) (b) (f)
(b) (a) (b)
(d) (b) (a), (b)
(b) (b)
(a), (b)
(b) (b) (a) (a) (a), (b)
(c) (e) (c)
(b) (b) (a) (b) (b)
(c)+(i) (b) (c) (f)+(i) (c)+(i)
(a) (e) (b) (f) (a)
(a) (a) (f) (b) (a)
(b) (d) (b) (b) (b), (c)
(a) (b) (b) (b) (a),(b)
(a), (f) (a), (b), (e), (f)
(a) (b) (a)
(c)
(b) (b) (a)
(c) (c)+(i) (c)
(b) (a) (c)
(a) (a) (a)
(b) (b) (c)
(b) (b) (b)
ip t
(a) (d) (c) (a) (c) (c) (c) (e)
M
an
us
(a) (a) (a) (a) (a) (a) (a) (a)
ed
Shah et al. (2011b,c) Naren and Ranade (2011) Shuai et al. (2011) Benyahia (2012) Shah et al. (2011a) Chalermsinsuwan et al. (2013) Zhou and Wang (2014) Shah et al. (2015) Baharanchi et al. (2015)
Granular thermal conductivity
NC (b) NC (f) (b) (b) (f) (b) NC NC
pt
Igci et al. (2008) Lu et al. (2009) Benyahia, (2009)
Granular temperature formulation
(a) (f) (b) (b) (a) (f) (b) (b) (g) (b)
Ac ce
Arastoopour et al. (1990) Pita and Sundaresan (1993) Dasgupta et al. (1994) Nieuwland et al. (1996a) Samuelsberg and Hjertager (1996) Hrenya and Sinclair (1997) Dasgupta et al. (1998) Mathiesen et al. (2000) Neri and Gidaspow (2000) Benyahia et al. (2001) Agrawal et al. (2001) Huilin et al. (2003) Yang et al. (2003b) Huilin et al. (2005) Andrews IV et al. (2005) Jiradilok et al. (2006) Benyahia et al. (2007) Wang et al. (2007) Almuttahar and Taghipour (2008)
Gas phase turbulent
cr
Table-2: CFD models of cold-flow gas-solid flow in riser Drag
(b) (f) (d), (e) (d) (a) (a), (f) (a), (b) (b), (c) (a), (b), (c), (d)
(a)
(d), (e) (a), (c), (e), (f) (a)
Gas phase:
Granular temperature:
Granular shear viscosity:
Solids Inlets:
Investigations:
(a) NC
(a) Partial differential equation (Lun et al., 1984) (b) Algebraic equation (Syamlal et al., 1993)
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k)
(a) (b) (c) (d) (e)
Effect of (a) drag models (b) KTGF models (c) viscous models (d) inlet and/or outlet boundary conditions (e) wall boundary conditions (f) flow conditions
(b) Standard k- mixture model
(c) Standard k- dispersed model
Granular thermal conductivity: (a) (b) (c) (d) (e)
Lun et al. (1984) Syamlal et al. (1993) Gidaspow (1994) Nieuwland et al. (1996) Agrawal et al. (2001)
Lun et al. (1984) Syamlal et al. (1993) Gidaspow (1994) Hrenya and Sinclair (1997) Nieuwland et al. (1996) Agrawal et al. (2001) Blazer et al. (1996) Cao and Ahmadi (1995) Frictional viscosity (Schaeffer, 1987) Dasgupta et al. (1994) Yang et al. (2004)
Real 3D 2D- two-sided solids inlet 2D-single solids inlet Periodic One common inlet for gas and solids phases
Walls (for solids): (a) No-slip (b) Partial-slip (c) Free-slip
Page 16 of 16 Page 16 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(d) Standard k- per-phase model
Collision dissipation:
Radial distribution function:
(a) Jinkin and Sevage (1983) (b) Lun et al. (1984) (c) Nieuwland et al. (1996)
(a) (b) (c) (d) (e) (f)
(e) Modified k- model (f) Other
Lun and Sevage (1986) Sinclair and Jackson (1989) Gidaspow (1984) Ma and Ahmadi (1986) Iddir and Arastoopour (2005) Carnahan and Starling (1995)
Note: (1) All the studies are for cold flow in the riser. Consequently, energy balance is ignored. (2) All the studies use governing equations by Anderson and Jackson (1967). (3) All the studies have ignored interphase forces other than drag.
Granular energy exchange: Gidaspow (1994) Koch and Sangani (1999) Cao and Ahmadi (1995) Blazer et al. (1996)
cr
ip t
(a) (b) (c) (d)
Ac ce
pt
ed
M
an
us
31
Page 17 of 17 Page 17 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
3.1. Eulerian-Eulerian model
2
In the formulation of the E-E model, local instantaneous equations are written for mass, momentum and
3
energy balances for each phase in each control volume. These equations can be solved by direct simulation
4
using size of control volume finer than the smallest length scales of the flow and a time step shorter than the
5
time scale of the fastest fluctuations. However, direct simulation would be computationally highly intensive.
6
Thus averaging of local instantaneous equations is applied. Generally, ensemble averaging (Enwald et al.,
7
1996) is applied to formulate averaged balance equations. The averaging of balance equations leads to
8
representation both phases as interpenetrating and continuous. The averaged equations then need extra
9
closure laws to close the set of equations. The closure laws for the continuous solids phase are derived using
10
the KTGF. Another set of closures are for the gas-solid interphase forces and heat exchanges. The closure
11
relations derived using the KTGF, interphase forces and heat exchanges are then incorporated into the
12
balance equations for the solids phase to achieve a close set of partial differential equations (PDEs).
13
3.1.1.
14
The governing equations for the E-E model consist of ensemble averaged Navier-Stokes equations (mass,
15
momentum and energy balances) for both the gas and catalyst phase. Several researchers have derived the
16
governing equations for two phase flows. Enwald et al. (1996) and van Wachem et al. (2001) have compared
17
the governing equations. Most of the gas-solid E-E simulation studies (as summarised in Table-2) have
18
referred their governing equations to Anderson and Jackson (1967). Therefore, in this study, we start from
19
the averaged balance equations of Anderson and Jackson (1967).
20
Mass balance equations:
21
The mass balance equation for the gas and solids phase can be written as:
23
cr
us
pt
ed
M
an
Governing equations
Ac ce
22
ip t
1
eq.(1)
eq.(2)
24
Where, εg and εs are the volume fraction of the gas and solids phases; ρg and ρs are the densities of the gas
25
and solids phases; ug and us are the velocities of the gas and solids phases; mgs and msg are the mass transfer
26
terms from the gas to solids phase and the solid to gas phase respectively; and Sg and Ss are the source terms.
27
In the E-E model, both phases are considered as interpenetrating continua; and both the phases share the
Page 18 of 18 Page 18 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
volume of flow domain. At any particular time therefore, a flow volume has εg of the gas phase and εs of the
2
solids phase with both volume fractions following:
3
For a simple gas-solid flow without mass transfer and mass sources, the continuity equations become:
ip t
4
eq.(3)
eq.(4)
us
cr
5
Momentum balance equations:
8
The conservation of momentum of the gas and solids phases is given by following equations:
ed
M
7
eq.(6)
Ac ce
pt
9
10
eq.(5)
an
6
eq.(7)
11
Where, P is the pressure that shared by both the phases;
and
are the stress tensor of the gas and solids
12
phases; g is the gravitational acceleration; FD, Flift and Fvm are the interphase forces such as drag force, lift
13
force and virtual mass forces respectively. Notably, in above equations, the force terms are positive for the
14
gas phase and negative for the solids phase.
15
Energy balance equations:
16
Energy conservation equation for the gaseous mixture phase:
Page 19 of 19 Page 19 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
2
eq.(8)
Energy conservation equation for the solids phase:
eq.(9)
ip t
3
Where, hg and hs are the specific enthalpies of the gas and solids phases respectively; qg and qs are the heat
5
fluxes in the gas and solids phases respectively; Sgh and Ssh are the energy source terms; and ΔQgs and ΔQsg
6
are the heat exchanges between the gas and solids phases respectively.
7
3.1.2.
8
3.1.2.1 Interphase forces
9
The momentum balance equations are coupled with interphase force terms, which represent the interactions
10
between the gas and dispersed solids phases. These forces are generally caused by a relative motion between
11
the gas and solids phases.
12
Drag force
13
The gas-solid drag force can be written as:
us
pt
ed
M
an
Constitutive equations
Ac ce
14
cr
4
eq.(10)
15
Where, β is the interphase momentum transfer coefficient and
is the difference between the
16
velocities of the gas and solids phases, usually termed as slip velocity. The value of the drag coefficient
17
depends on volume fraction, particle diameter and slip velocity; and in turn, particle Reynolds number.
18
Several drag coefficient functions (Ergun, 1952; Wen and Yu, 1966; Gibilaro et al., 1985; Syamlal and
19
O’Brien, 1987; Syamlal et al., 1993; Di Felice, 1994; Gidaspow, 1994; Li and Kawauk, 1994; Andrews IV et
20
al., 2005; Beetstra et al., 2007; Igci et al., 2008; Igci and Sundaresan, 2008; Tenneti et al., 2011) have been
Page 20 of 20 Page 20 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
proposed to cater different the gas-solid flow conditions. Comparison of the available drag models is
2
explained in section-3.2.1.
3
Virtual mass force
4
In the gas-solid flows in riser, the solids phase accelerates relative to the gas phase in the developing section.
5
The inertia of the gas-phase encountered by the accelerating solids phase exerts a virtual mass force on the
6
particles.
ip t
1
eq.(11)
The term
denotes the phase material time derivation of the form,
M
an
8
us
cr
7
eq.(12)
ed
9
The virtual mass effect is significant when the secondary phase density is much smaller than the primary
11
phase density (i.e. bubble column). But in the riser, the solids phase density is much higher than the gas
12
phase, and therefore, most of the gas-solid models do not account for the virtual mass force.
13
Lift force
14
The lift force acts on the suspended particle mainly due to velocity gradients in the gas-phase, relative
15
velocity between the gas and solids phase, and rotation of the particles.
Ac ce
16
pt
10
eq.(13)
17
Where,
18
particle spacing and the value of drag force is much larger, the lift force has been ignored in the previously
19
proposed flow models.
is the lift coefficient.Assuming that the FCC particle diameter is much smaller than the inter-
Page 21 of 21 Page 21 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
3.1.2.2. Gas phase stresses
2
The stresses in the gas phase are modelled by using either the laminar or turbulent model.
3
Considering laminar flow, the correlation for the gas phase stress tensor can be written as:
If the turbulent flow for the gas phase is considered then the stresses in the gas phase can be written as:
cr
5
eq.(14)
ip t
4
eq.(15)
an
us
6
Where, μg is the laminar viscosity; μt,g is the turbulent viscosity and k is the kinetic energy of the gas phase.
8
To close the equation of the turbulent stresses, a correlation for the turbulent viscosity and kinetic energy are
9
required.
M
7
Owing to their inherent complexity and scale of flow domain, the standard k- turbulent model has been
11
used (Dasgupta et al., 1998; Almuttahar and Taghipour, 2008; and Shah et al., 2011b) for the turbulent
12
properties. There are three different types of formulations that are possible for the standard k
13
model. The simplest k
14
for dissipation rate ( ). These equations are formulated for a single pseudo phase made of the mixture of the
15
gas and solids phases. This model then uses mixture properties to calculate the turbulence properties
16
(Dasgupta et al., 1998). The conservation equation for the kinetic energy and dissipation rate can be written
17
as:
turbulent
Ac ce
pt
ed
10
model has the two equations, one for the turbulent kinetic energy (k) and the other
Page 22 of 22 Page 22 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
and that for the dissipations rate can be written as:
ip t
2
eq.(16)
eq.(17)
cr
3
Where, m denotes the mixture properties; k is the turbulent kinetic energy; is the turbulent dissipation rate;
5
ρm is the mixture density (
us
4
M
an
); and μm is the mixture viscosity (
ed
6
9
eq.(18)
eq.(19)
Ac ce
8
pt
7
).
eq.(20)
10
In above equations, Cμ, C1,ε, and C2,ε are constants; whereas σk and σε are the turbulent Prandtl numbers
11
(Lauder and Spalding, 1972). Dasgupta et al. (1994) have used 0.09, 1.45, 2.0, 1.3 and 1.0 for Cμ, C1ε, C2ε, σk
12
and σε respectively.
13
More complex k- turbulent model for the gas-solid flow is derived by accounting for the effect of dispersed
14
solid particles on the gas phase turbulent kinetic energy and dissipation rate (Bolio et al., 1995; Cao and Page 23 of 23 Page 23 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Ahmadi, 1995). In this model, the equations for k and
are solved for the gas phase and the interphase
2
turbulent momentum transfer is accounted using a closure term. This model is known as the k- dispersed
3
model, and its equations can be given by
ip t
1
eq.(21)
us
cr
4
6
Where,
are source terms for the influence of the dispersed phases on the continuous phase. The
7
equations for these source terms are (Elghobashi and Abou-Arab, 1983; Simonin and Viollet, 1990):
ed
M
and
eq.(22)
an
5
eq.(23)
Ac ce
pt
8
9
eq.(24)
10
Where, kgs is the variance in the velocity of the gas phase and dispersed solids phase; and C3e is a constant
11
having a value of 1.2.
12
Additionally, if the presence of solids has a significant impact on gas phase turbulence, the above equations
13
of k and
14
phase model.
for the gas phase are also solved for the solids phase. This model is then known as the k- per
Page 24 of 24 Page 24 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
While the k- model is capable of useful predictions for engineering design and optimization, it requires
2
careful parametrization. There are two basic assumptions in the k- model, i.e. (i) isotropic turbulent
3
viscosity and (ii) simplifications due to Reynolds averaging. The assumption of isotropy is valid for many
4
flows, but not for those with strong separation or swirl. In flows where, the flow regime transitions occur, the
5
k- models do not provide the required accuracy. Therefore, other modelling approaches such as Reynolds
6
stress model (Lain et al., 2002; Lain and Garcia, 2006; Kuan et al., 2007 and Chu et al., 2011) and large eddy
7
simulation (LES) model (Mathiesen et al., 2000; Yamamoto et al., 2001; Huilin and Gidaspow, 2003; Ibsen
8
et al., 2004 and Vreman et al., 2009) have been used for the gas-solid flow. Such models are still
9
computationally expensive, particularly for full scale simulation of riser flow. Mostly, the RSM model is
10
used for gas-solid flow with strongly anisotropic flow of gas phase such as those in cyclone and pipe bends;
11
while the LES model has been used for dilute gas-solid flow in pipes.
12
In the RSM, the stress in the gas phase is given as sum of laminar and turbulent stresses. The laminar stress–
13
strain tensor can be given by eq. (14), while the turbulent stresses can be given by:
cr
us
an
M
ed pt
14
ip t
1
eq. (25)
Where, τt,g is the turbulent stresses in the gas phase; and
16
Similar to the k- two equation model, the RSM can also be implemented as either dispersed or mixture
17
model. In the RSM-dispersed model, the transport equation for the Reynolds-stress term is solved for only
18
the gas phase. The RSM transport equation can be given as:
Ac ce
15
is the Reynolds-stresses in the gas phase.
Page 25 of 25 Page 25 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1 2
eq.(26)
Where, u’ is the fluctuating velocity component; is the turbulent dissipation rate; and
4
stress terms due to turbulent interactions between phases.
cr
us M
is the Kronecker delta.
is given by modified model of Simonin and Viollet (1990):
ed
Where,
eq.(27)
eq.(28)
Ac ce
pt
7
can be given as (Cokljat et al. 2003):
an
5
6
is the Reynolds
ip t
3
8
Where, β is the drag coefficient; kgc is the gas-solid phase velocity covariant; k is the turbulent kinetic
9
energy; urel and udrift are relative and drift velocities respectively.
10
In the RSM-mixture model, the Reynolds-stress transport equation for the gas phase is solved, but without
11
term. Furthermore, properties of phase such as density and velocity are represented as mixture
12
properties. The RSM solves transport equations for six independent Reynolds stresses and one for the
13
turbulent dissipation. The transport equation of the turbulent dissipation is same as eq. (22). The RSM is
14
computationally more expensive than the k- model. The use of RSM avoids assumption of isotropy, but it
Page 26 of 26 Page 26 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
still involves Reynolds averaging which may cause error estimation of double and triple correlation of
2
turbulent velocities.
3
Because of the restrictive assumptions made in the k- and RSM models, some deviations (upto 25%) occur
4
in the prediction of the three components of the mean flows. However, the major impact of the assumption is
5
seen in the predictions of k and
6
of the rates of agglomeration, breakage, heat transfer and mass transfer need to use the local values of k and
7
and hence the need for accurate predictions. The conservation equations for k and consists of 5 terms : rates
8
of (i) convective transport (ii) turbulent transport (iii) viscous transport (iv) productions and (v) dissipation.
9
While solving these conservation equations, the values of individual terms take such values that the overall
10
balance is satisfied. However, normally the rates of dissipation and turbulent transport are under-predicted
11
upto 70%. On contrast the rates convective transport and production together get over-predicted so as to
12
match the balance. These incorrect predictions of individual terms results into incorrect predictions of k and
ip t
cr
us
an
M
ed
in the range of ± 20% to ± 200% (Zoheb Khan, 2016).
pt
13
profiles. It may be pointed out that most of the models for the predictions
In the LES model, the flow of the gas is divided in two parts, (i) large-scale flow resolved by the
15
computation and (ii) sub-grid scale (SGS) fluctuations. The effect of SGS fluctuations on the large-scale
16
flow is accounted for by applying SGS closure for turbulent stresses. The gas phase stress can be given as
17
(Yamamoto et al., 2001):
18
Ac ce
14
19
Where µ is the gas phase viscosity; Sij is the strain rate tensor,
20
of the gas phase.
eq.(29)
is the SGS stress and u is the velocity
Page 27 of 27 Page 27 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
The strain rate tensor can be given as:
2
eq.(30)
Using the Smagorinsky model (Smagorinsky, 1963), the SGS stress can be given as (Yamamoto et al.,
4
2001):
cr
eq.(31)
us
5
ip t
3
Where, Cs is the Smagorinsky constant;
8
be given as
is the filter size and fs is the damping function. The filter size can
M
7
pt
ed
and the damping coefficient can be given as:
eq.(33)
Ac ce
9
eq.(32)
an
6
10
Where, y+ is the distance from the wall in wall unit.
11
While we have provided an overview of LES formulation for the continuous phase, in multiphase flows,
12
even at laminar conditions, additional challenges arise due to small time and length scales involved in the
13
dispersed flow. These scales are further magnified in presence of turbulence. Thus it is important to capture
14
the microscopic phenomena that accounts for all the small scale interactions for development of accurate
15
multiphase model (Fox, 2012).
16
3.1.2.3. KTGF closure models
17
Using the continuum approach for the solids phase, the stresses in the solids phase are derived by making an
18
analogy between a random motion of particles and thermal motion of gas molecules. This approach is widely
Page 28 of 28 Page 28 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
known as the KTGF (Jenkins and Savage, 1983; Lun et al., 1984). Unlike the kinetic theory of the dense gas
2
(Chapman and Cowling, 1970), the KTGF accounts for inelasticity of particle-particle collisions. The
3
intensity of the particle velocity fluctuations determines the stresses, viscosity and pressure of the solids
4
phase. Taking the analogy with the thermodynamic temperature of the gas, the random motion of particles is
5
represented by “granular temperature”, which is proportional to the mean quadratic velocity of the particle:
eq.(34)
cr
ip t
6
eq.(35)
us
7
Where, Θ is the granular temperature; c is fluctuating velocity of particle and EΘ is the fluctuating energy due
9
to random motion of particles.
an
8
The solids phase stresses depends on the magnitude of this fluctuating particle velocity. Thus, conservation
11
of the pseudo thermal energy of the solids phase associated with the fluctuating velocity is required. This
12
conservation equation can be given as:
ed
M
10
eq.(36)
Ac ce
pt
13
14
Where,
15
dissipated energy due to particle collision;
16
energy; and
17
The collisional dissipation energy can be written as (Lun et al., 1984):
represents generation of the energy by the solid stress tensor;
represents
is diffusion coefficient (or conductivity) for the granular
represents exchange of fluctuating energy between the gas and solids phases.
Page 29 of 29 Page 29 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
2
eq.(37)
and the diffusion coefficient for the granular energy can be written as (Gidaspow, 1994):
eq.(38)
cr
ip t
3
4
Where,
5
The exchange of kinetic energy of random fluctuations in particle velocity to the gas phase is written as
6
(Gidaspow, 1994):
M
an
us
.
eq.(39)
ed
7
Syamlal et al. (1993) suggest an algebraic form of the granular temperature equations by neglecting
9
convection and diffusion terms. The algebraic form of the granular temperature equation can be written as
11
(Syamlal et al., 1993):
Ac ce
10
pt
8
eq.(40)
12
In above equations, ess is the coefficient of restitution for inelastic collisions between the particles; go,ss is the
13
radial distribution function, Ps is the solids pressure and
14
coefficient of restitution depends on the material properties, and its value for the FCC catalyst is not known.
15
As a result, this parameter works as an empirical tunning parameter. In previous simulations, the values of
16
the restitution coefficient between 0.7-0.99 have been used. The radial distribution function describes the
17
probability of finding two particles in close proximity. It is a correction factor that modifies the probability
18
of collisions between particles when the granular phase becomes dense. Its correlation given by Lun et al.
19
(1984) can be written as,
is the stress tensor for the solids phase. The
Page 30 of 30 Page 30 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
2
eq.(41)
The stress tensor for the solids phase can be written as:
eq.(42)
cr
ip t
3
The solids pressure (Ps) represents the normal solids phase forces due to particle-particle interactions. The
5
solids pressure term is composed of a kinetic term and a second term due to the particle collisions. The solids
6
pressure term prevents the solids volume fraction from exceeding the maxing packing limit. The correlation
7
for the solids pressure can be written as (Lun et al., 1984):
an ed
written as (Lun et al., 1984):
pt
11
The solids bulk viscosity (λs) is the resistance of particle suspension against the compression, and it can be
eq.(44)
Ac ce
10
eq.(43)
M
8
9
us
4
12
The solid shear viscosity (µs) is made up of the collisional, frictional and kinetic parts. All the models for
13
solid shear viscosity yield practically the same solid shear viscosity at solid volume fraction greater than
14
0.25. For lower volume factions of the solids, the models start deviating from one another (Ahuja and
15
Patwardhan, 2008). The Gidaspow model (Gidaspow, 1994) for the solid shear viscosity neglects the
16
inelastic nature of particle collisions in the kinetic contribution of the total stress and the Gidaspow’s model
17
has been used in several previous studies. The model equations for the granular shear viscosities and its
18
collisional and kinetic components can be written as,
19
eq.(45)
Page 31 of 31 Page 31 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(Gidaspow, 1994)
eq.(46)
eq.(47)
ip t
2
3.1.2.4. Frictional viscosity
4
At high solid volume fraction, frictional stress resulting from inter-particles contacts is significant, and
5
therefore, frictional part is considered to calculate the solid shear viscosity. The frictional viscosity is
6
generally calculated using Schaeffer’s (Schaeffer, 1987) model.
an
us
cr
3
(Schaeffer, 1987)
eq.(48)
M
7
Where, θ is the angle of internal friction; Pfri is the frictional pressure; and I2D is the second invariant of the
9
deviatoric stress tensor. The frictional pressure is given as (Johnson and Jackson, 1987):
ed
8
eq.(49)
Ac ce
pt
10
11
Where, εs,min is the minimum solid volume fraction from where frictional stresses become important. Fr, n
12
and p are empirical constants which depend on material properties. The values of Fr, n and p are often taken
13
to be 0.05, 2 and 3 respectively (Ocone et al., 1993). Syamlal et al. (1993) have suggested a model the
14
frictional pressure, which is:
15
eq.(50)
16
In the above equation, A and n are empirical constants and their values are taken to be 1025 and 10
17
respectively (Syamlal et al. 1993).
Page 32 of 32 Page 32 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Notably, several alternate models are available for the radial distribution function, solids pressure, bulk
2
viscosity, shear viscosity, frictional viscosity, granular conductivity and granular energy diffusion
3
coefficient. Van Wachem et al. (2001) and Ahuja and Patwardhan (2008) have reviewed the available
4
models by plotting their values against the solids volume fraction.
5
3.1.2.5. Heat exchange
6
In the energy conservation equations, ΔQgs and ΔQsg are heat exchanges between the gas and solids phases.
7
The heat exchange between the gas and solids phase can be given as,
cr
ip t
1
eq.(51)
us
8
Where, hsg (= hgs) is heat transfer coefficient; Asi is interfacial area between the gas and solids phases; Ts is
10
the temperature of the solids and Tg is the temperature of the gas phase. The heat transfer coefficient can be
11
correlated to the Nusselt number by,
M
an
9
eq.(52)
ed
12
Where, Kg is the thermal conductivity of the gas phase; Nus is the Nusselt number for the solids phase; and dp
14
is the diameter of the particles.
15
Generally, an empirical correlation of Ranz and Marshall (1952) is used to derive the heat transfer
16
coefficients.
18
Ac ce
17
pt
13
(Ranz and Marshall, 1952)
Where, Rep is the particle Reynolds number; and Pr is the Prandtl number.
19
20
eq.(53)
eq.(54)
Where, Cpg is the heat capacity of the gas.
Page 33 of 33 Page 33 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
The Ranz and Marshall model is based on the heat transfer between single particles moving in a stationary
2
gas. Table-3 summarises several other empirical models that have been derived for the fluidized bed
3
conditions. Notably, there is no model available for the heat transfer between the gas and solids under the
4
flow conditions of the riser. Table-3: Models for gas-solid heat transfer coefficient Ranz and Marshall (1952) Ranz and Marshall (1952) Kothari (1967)
ip t cr
Fixed bed
an
+
us
Gunn (1978)
Single particle
M
Kunni and Levenspiel (1991)
Fluidized bed (Rep = 1 -100) Fluidized bed (Rep = up to 105 Fluidized bed
5 3.2.
Effect of constitutive models on flow predictions
7
Table-2 shows that simulation studies have used different combination of the constitutive models, and the
8
different selections of the constitutive models have caused significant variations in flow predictions.
9
Therefore in this section, the impact of closure models on predictions of gas-solid riser flow is analysed.
Ac ce
pt
ed
6
10
3.2.1.
Effect of drag models
11
Available drag models for the riser flow can be divided in three categories, namely (i) conventional drag
12
models, (ii) multi-scale or cluster based drag models, and (iii) those derived from direct numerical
13
simulations. Table-4 summarises the available drag models with their correlations, methods of their
14
derivation and shortcomings. Conventionally, the Wen-Yu model (Wen and Yu, 1966), Gidaspow model
15
(Gidaspow, 1994), and Syamlal and O’Brian model (Syamlal and O’Brien, 1987 and Syamlal et al., 1993)
16
are used for riser simulations. The Gidaspow model is a combination of the Ergun (for εg<= 0.8) and the
17
Wen-Yu (for εg> 0.8) correlations. The Syamlal and O’Brian model is based on the velocity-volume fraction
18
correlations derived using sedimentation and fluidization experiments of Richardson and Zaki (1954) and
19
Garside and Al-Dibouni (1977). The conventional drag models are empirical correlations that derived from
20
experimental observations on settling velocities of a single particle in the liquid and pressure drops inside
21
packed bed of the solids. However, the flow in riser exhibits the formation of clusters (Yerushalmi et al., Page 34 of 34 Page 34 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1976; Yerushalmi and Cankurt, 1979; Hario and Kuroki, 1994), which coexist and interact with surrounding
2
dilute phase. Due to the presence of clusters, apparent diameter of the solids becomes higher, and terminal
3
velocity of the particles in the cluster then becomes higher. This leads to down flow of solids and higher
4
segregation of the solids near wall (Mueller and Reh, 1993; Nieuwland et al., 1996a, Noymer and
5
Glicksman, 2000; and Pandey et al., 2004). Furthermore, due to formation of the clusters, the gas tends to
6
flow around the clusters without penetrating inside them. This result in significant less resistance to flow (Li
7
and Kwauk, 1994) leading to reduction in gas-solid drag. In previous studies (Yang et al., 2003a,b; Andrew
8
IV et al., 2005; and Shah et al., 2011b,c), use of the conventional drag models in coarse grid riser simulations
9
has not captured axial and radial distribution phases accurately. These studies have attributed this drawback
10
to the use of coarse grid, where spatial resolution is insufficient to capture clusters. Agrawal et al. (2001)
11
have shown important of grid resolution to capture clusters. However, fine grid simulation of industrial-scale
12
riser is computationally intensive. Therefore, drag models specific to cater coarse grid riser simulations have
13
been derived by using different multiscale approaches. Wang (2009) has categorized methods for deriving
14
multiscale drag models for the Eulerian simulations of Geldart-A particles in six categories, i.e. empirical
15
correlation method, scaling factor method, structure-based method, modified Syamlal and O’Brien drag
16
correlation method, EMMS-model-based method, and correlative multi-scale method. Out of these methods,
17
the sub grid scale (SGS) approach (Sundaresan and co-workers) and energy minimization multi-scale
18
(EMMS) approach (Li and Kwauk and co-workers) have been often used in several previous studies.
19
In the SGS approach (Agrawal et al., 2001; Andrews IV et al., 2005; Igci et al., 2008; Igci and Sundaresan,
20
2008), the drag is derived by using numerical results from extremely fine grid simulations using a periodic
21
domain. Andrew IV et al. (2005) have conducted fine grid simulations of only a section of riser using
22
periodic boundary conditions, and the results have then been used to derive closures of the solids phase and
23
gas-solid drag. . Igci et al. (2008) have extended Andrew IV et al.’s work to show the effect of grid size,
24
known as filter size, on the closures derived from the fine grid simulations. These closures are then used in
25
coarse grid simulation of a full-scale riser (Benyahia, 2009). An alternative approach is the energy
26
minimization multi-scale (EMMS) model (Li and Kwauk, 1994; Xu and Li, 1998; Li et al., 1999; Yang et al.,
27
2003b). The EMMS approach divides the flow in three pseudo phases i.e. dilute phase, dense phase and
28
suspended clusters. The dilute phase represents individual particle in gas. The dense phase represents
29
particles and gas residing inside clusters; whereas the suspended clusters phase represents individual cluster
30
as a whole entity that moves in the gas phase. The EMMS calculates drag in each three pseudo phases by
31
using conventional Wen-Yu drag model. Drags in the dilute phase and dense phase are calculated by using
32
particle diameter, while drag in the suspended cluster phase is calculated by using cluster diameter. Addition
33
of three drags in three pseudo phases represents resultant EMMS drag in dilute gas-solid flow in riser.
34
However, calculation of the EMMS drag needs cluster diameter, fraction of solids forming the cluster and
35
slip velocities in three pseudo phases. Thus, EMMS drag largely relies on cluster parameters. Cluster
36
diameter in the EMMS is mostly calculated by empirical correlations and correlations based on different
37
hypothesis (Shah et al., 2011b); while fraction of solids and superficial velocity in cluster phase are
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 35 of 35 Page 35 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
calculated by solving mass and energy balance equations for each pseudo phase. For detailed equations and
2
calculations of the EMMS drag, readers are referred to various published articles of Li and co-workers (Li
3
and Kwauk, 1994; Xu and Li, 1998; Li et al., 1999; Yang et al., 2003a, b; Wang, 2008; and Shah et al.,
4
2011b, c). Recently, numerical observations from DNS have been used to derive gas-solid drag correlations
5
(Hill et al., 2001b; Benyahia et al., 2006; Beetstra et al., 2007; Tenneti et al., 2011). DNS studies are
6
conducted by using lattice Boltzmann method (LBM) or immersed boundary method (IBM). In these studies,
7
particle resolved simulations of a periodic domain with several randomly positioned particles are conducted
8
to capture interphase exchange forces at the boundary of each particle. The numerically calculated particle-
9
scale forces are then used to calculate an averaged drag force in the flow domain.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 36 of 36 Page 36 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Table-4: Gas-solid drag models Derived from Conventional drag models Gidaspow (1994)
Do not capture axial profile of the solid volume fraction (Yang et al., 2004; Wang et al., 2008; Benyahia, 2009; Shah et al., 2011; Baharanchi et al., 2015). Previous studies attributed this drawback to their inability to account for effect of clusters in coarse grid simulations. Actually, The Wen-Yu model is applicable to dilute flows (εg close to 1); while the Ergun is used for dense flows. Even the Syamlal-O’Brien is derived from a single particle drag. At intermediate range of solids A drag for multi-particle system is derived from single volume fraction (0.7 ≤ εs ≤ 0.95), these drag models particle drag, Dalla Valle (1948)’s equation for drag significantly over predict the drag force than the values coefficient for a single particle, Richardson and Zaki of multi-scale drag models. (1954)’s correlation between ratio (Vr) of terminal velocity of multi-particle to that of single particle with the gas volume fraction, and Garside and Al-Dibouni (1977)’s correlation between Vr and particle Reynolds number.
Uses Ergun model for εg less than 0.8 and Wen-Yu model for εg greater than 0.8. The Wen-Yu model has been derived from settling experiments of single particle in liquids, while the Ergun model has been derived pressure drop measurement in packed bed.
for εg> 0.8 (Wen and Yu, 1966) for εg<= 0.8 (Ergun, 1952b) If Res< 1000 Schiller &Nauman (1935)
ip t
= 0.44 If Res≥ 1000
us
cr
Syamlal-O’Brien (1987)
Multiscale drag models EMMS (Li and Kwauk, 1994)
for
> 0.748
<= 0.748
ed
for where,
M
an
if ε ≤ 0.85 if ε > 0.85 “a” and “b” are tunning constants that must be calibrated to match the minimum fluidization velocity for the solid particles in question (Syamalal and O’Brien 2003). Based on Umf calculated from the Ergun equation for 76 μm FCC particle and 1712 kg/m3 density, the values of “a” and “b” are 0.8 and 2.66.
is a correction factor and it can be calculated as,
Uses EMMS algorithm based on energy minimization of suspension and transportation in gas-solid riser flow. It requires cluster parameters such as fraction of solids in cluster, solid volume fraction of cluster, cluster diameter and slip velocity in clusters.
for 0.748 <εg<= 0.83
pt
for 0.83 <εg<= 0.86 for 0.86 <εg<= 0.94
Ac ce
for εg> 0.94 The equations for the corrections factors are for high solids flux (Gs = 489 kg/m2s, and ug = 5.2 m/s) as reported by Shah et al. (2014). These functions need to be recalculated using the EMMS algorithm (Li et al., 2011; Yang et al., 2003; Shah et al., 2011) for a given flow conditions and physical properties. SGS Anderew IV et al. (2005); Igci et al. (2008); Benyahia (2009)
Shortcomings
for
≤ 0.0.04083
Simulation results using very fine grid in periodic domain of a size comparable to control volume of the continuum simulations.
Capture both axial and radial profiles with some quantitative discrepancies (Yang et al., 2003; Benyahia, 2009; Igci et al., 2008 Shah et al., 2011). The EMMS requires cluster properties, of which detail measurements are not available. Furthermore, the EMMS correlation needs to be derived for riser with different flow conditions. Implementation of EMMS algorithm in each control volume using local flow conditions is computationally expensive. The SGS requires find grid simulations of a small flow domain and thus, the final drag function depend largely on grid size, as well as the simulated flow conditions (Igci et al., 2008; Benyahia, 2009). For riser simulation, this approach is computationally more expensive.
for 0.0.04083< ≤0.2589 for
>0.2589
Above drag equations are based on for high solids flux (Gs = 489 kg/m2s, and ug = 5.2 m/s) as reported by Benhyaia (2009) using dimensionless filter size gΔ/Vt2 = 2.056 where Δ = grid size = 0.014m, Vt = 18.46 m/s Models from DNS
Page 37 of 37 Page 37 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Beetstra et al. (2005)
LBM simulations of randomly arranged particles in periodic flow domain. Drag on individual particles is calculated at particle Reynolds number range of 21-105. The calculated drag values are used to derive a correlation between the drag, gas volume fraction and particle Reynolds number. IBM simulations of randomly arrange particles in periodic flow domain are conducted with flow conditions varies in a range of particle Reynolds number of 0-300. The calculated values of fluid-solid drag are used to derive drag correlation.
cr
ip t
Tenneti et al. (2011)
These drag models have not been extensively investigated by conducting riser simulations.
Ac ce
pt
ed
M
an
us
10
Drag values from Beetstra et al. (2005) are quite different from those of Tenneti et al. (2011) at high solids volume fraction and particle Reynolds number. Tenneti et al. (2011) have attributed the difference to insufficient grid resolution of Beetstra et al.’s simulations. Tenneti et al. (2011) have highlighted a need of high grid resolution to capture boundary layers between two particles, particularly at high solids volume fraction. Tenneti et al.’s drag model is based on particles arranged at random positions. Shah et al. (2013) have shown that drag values from particles arranged in cluster configurations are significantly lower than random configuration.
Page 38 of 38 Page 38 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
In Table-4, the correlations for the EMMS and SGS are available for the flow conditions of high flux riser
2
(Gs = 489 kg/m2s, ug = 5.2 m/s, ρs = 1712 kg/m3, Inlet solid volume fraction = 0.5, Dt = 0.2 m and H = 14.2
3
m) of Knowlton et al. (1995).Using the correlations of Table-4, values of the drag coefficients are calculated
4
for a range of volume fraction at a given slip velocity of 0.57 m/s. The comparison of the calculated drag
5
coefficients is shown in Figure-7. It can be seen that the drag coefficients from the conventional models and
6
DNS agree with each other over a range of the volume fraction. But at the lower gas volume fraction, the
7
drag coefficients from the DNS studies are higher than those from the Wen-Yu and Ergun models. The
8
values from the Syamlal-O’Brien model are lower than those from the Wen-Yu, Ergun, BVK (Beetstra, van
9
der Hoef and Kuiper) and PUReIBM (Tenneti et al., 2011) models. The drag coefficient calculated from the
10
EMMS model is significantly lower in an intermediate range of the gas volume fraction of 0.7-1. At
11
minimum fluidizing and maximum gas volume fraction, the mathematic formulation of the EMMS model
12
forces the drag values to approach the values calculated by the conventional drag models. The values from
13
the SGS model are significantly lower than other drag models over entire range of the gas volume fraction,
14
and the SGS values further drifting lower at the gas volume fraction close to 1. While explanation of the
15
behaviour of the SGS model is not clear, it is noteworthy that the SGS correlation is a weak function of the
16
particle Reynolds number (Benyahia, 2009) and it relies mainly on grid size (known as filter size) used in
17
fine grid periodic simulations. In previous studies, the lower drag values from both the EMMS and SGS
18
models have been attributed for their ability to account for the clusters or volume fraction heterogeneity.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Figure-7: Comparison of drag coefficients (Gs = 489 kg/m2s, ug = 5.2 m/s, Uslip = 0.57 m/s) 19
Yang et al. (2003a,b); Wang et al. (2008); Benyahia (2009), Shah et al., (2011b, c), Benyahia (2012), Wang
20
et al. (2012), Shah et al. (2015) and several others have compared predictions of axial and radial profiles
21
from the EMMS models with those from the Gidaspow drag model. Similarly, Andrews IV et al. (2005),
22
Benyahia (2009), Igci and Sundaresan (2008), and Igci and Sudaresan (2011) have compared predictions
23
from the SGS model with those from the Gidaspow and experimental data. Comparison by Shah et al. (2015) Page 39 of 39
Page 39 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
is shown in Figure-8. For a low solid flux riser, the predicted axial profile from the Gidaspow model shows
2
that the values do not vary along height. This axial profile from the Gidaspow drag model shows complete
3
disagreement with the experimental data both qualitatively and quantitatively. On the other hand, the EMMS
4
model gives a dense bottom with higher solid volume fraction and lower values at the top (Figure-8) forming
5
the S-shape axial profile with reasonable qualitative agreements with the experimental data. For high flux
6
solids riser (Figure-8b), the pressure drop profiles predicted by both the Gidaspow and EMMS model
7
qualitatively agrees with the experimental data. Quantitatively, the predictions from the EMMS model are
8
closer to the experimental data at the bottom; while at the top, the values from the Gidaspow model are in
9
close agreement with the experimental data.
Ac ce
pt
ed
M
an
us
cr
ip t
1
(b)
(a)
Page 40 of 40
Page 40 of 114
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
an
(c) (d) Figure-8: Effect of drag model: (a) and (b) axial profiles of solids volume fraction at low solids flux and pressure drop at high solids flux respectively; (c) and (d) radial profiles of solids volume fraction at low and high flux respectively (Shah et al., 2015). Shah et al. (2015) have also compared the predicted radial profiles from the EMMS and Gidaspow models
2
with the experimental data (Figure-8c and d). The radial profiles from the EMMS model also give reasonable
3
qualitative agreements with the experimental data for two (low and high solids flux) flow conditions.
4
Quantitatively, the predictions from the EMMS still show some discrepancies between the predicted values
5
and experimental data. Shah et al. (2015) attributed these discrepancies to assumptions and empirical
6
correlations involved in calculations of the cluster parameters. Benyahia (2009) have also found
7
discrepancies between their predictions from the SGS drag model and experimental data, and the author has
8
attributed the discrepancies to the dependence of the SGS drag on the selection of the filter size. The
9
numerical drag models proposed by Beetstra et al. (2007) and Tenneti et al. (2011) have still not been tested
ed
pt
Ac ce
10
M
1
by conducting riser simulations.
Page 41 of 41
Page 41 of 114
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(a)
an
(b) Figure-9: Effect of gas phase flow models, (a) radial profile of volume fraction of solids phase and (b) axial profile of volume fraction of solids phase. 3.2.2.
Effect of gas phase viscous model
2
Experimental data shows existence of both dense and dilute conditions inside the riser. The flow of the gas
3
phase in the dilute region is expected to be turbulent, while that in the dense region with narrow inter-particle
4
distances and in boundary layer around particles and clusters can be either laminar or turbulent. On the other
5
hand, clusters of the solid particles observed in riser can significantly dampen or enhance fluctuating velocity
6
of the gas phase. Hence, laminar model has also been used to model stresses in the gas phase. The gas phase
7
Reynolds number based on diameter of riser generally varies from 104 to 106; while the particle Reynolds
8
number that calculated from diameter of particle is less than 100. As a result, both the laminar and turbulent
9
models for the gas phase stresses are applied in CFD studies.
Ac ce
pt
ed
M
1
10
Almuttahar and Taghipour (2008) have compared predictions from the laminar and k- mixture turbulent
11
models. It can be seen in Figure-9a that radial distributions of the solids volume fraction predicted by the
12
laminar model are closer to the experimental values than those from the k- turbulent model; whereas, the
13
predicted axial profiles from both the laminar and turbulent models show wide discrepancies with the
14
experimental values (Figure-9b). Almuttahar and Taghipour (2008) have attributed these discrepancies to the
15
considered 2D geometry, and recommend 3D simulations for accurate predictions. Benyahia et al. (2005)
16
have investigated three different gas-solids turbulence models. Two of these models are those of Balzer et al. Page 42 of 42
Page 42 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(1996) and Cao and Ahmadi (1995) that account for the effect of gas turbulence through the use of a
2
modified k- model, whereas the third simulated model is of Agrawal et al. (2001) that does not model the
3
gas phase turbulence. Benyahia et al. (2005) have found only minor difference between the predicted radial
4
profiles from these models with and without gas phase turbulence.
5
Shah et al. (2015) have also compared predictions from the laminar, k-ε dispersed model and k-ε per phase
6
model for both high and low flux conditions. Shah et al. (2014) have found that predicted values from both
7
the laminar and k- dispersed models give reasonable agreement with the experimental data; while the per-
8
phase option slightly under predicted the values in the dense section and over predicted at the top section
9
(Figure-10a). For the high flux condition (Figure-10b), all three models predicted similar time-averaged axial
10
pressure drop profiles. Figure-10(c) and (d) shows the impact of three viscous stress models on instantaneous
11
contour plots of the solids volume fraction captured at 20s flow time. For both low and high flux conditions,
12
the contour plots from laminar and k- dispersed models give continuous flow structures with “cluster” and
13
“steamer” like solid structures appearing throughout the height. But the contour plots from k- per-phase
14
model give discontinuous solid structures, where zones of higher and lower solid volume fractions appear to
15
be stretched all the way from the wall to the centre of the riser. Furthermore for the low solids flux condition,
16
the contour plots from only laminar and k- dispersed model show a clear distinction between the dense
17
bottom and dilute top sections. Shah et al. (2015) have further compared instantaneous values of turbulent
18
kinetic energy of the gas phase and granular temperature of the solids phase at different heights. While Shah
19
et al. (2015) have shown different in predictions from the k- dispersed and per-phase models, they have not
20
given a detail analysis on reason for these differences.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 43 of 43
Page 43 of 114
(b)
Ac ce
pt
ed
M
(a)
an
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Laminar
k-ε dispersed
k-ε per-phase
Page 44 of 44
Page 44 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Laminar k-ε dispersed k-ε per-phase (d) High solids flux (Gs = 489 kg.m2s)
2
(c) Low solids flux (Gs = 14.3 kg.m s)
Figure-10: Effect of gas phase flow models on axial profiles, (a) volume fraction of solids phase at low solids flux condition; (b) pressure drop at high solids flux condition; and (c) and (d) Instantaneous contour plots of solid volume fraction (Shah et al., 2015). 3.2.3.
Effect of KTGF closure models
2
Table-2 shows that previous simulation studies have used the E-E model with the KTGF, where two types of
3
granular temperature conservation equations such as partial differential and algebraic have been used. In
4
either of these two selections, the model requires closure laws for the granular energy dissipation rate, radial
5
distribution function and solids bulk viscosity. If the partial differential equation for the granular temperature
6
is applied, then the closure law for the granular temperature conductivity is also required. Several models are
7
available for these KTGF closures. van Wachem et al. (2001) and Ahuja and Patwardhan (2008) have
8
explicitly calculated values of the solids phase KTGF properties from alternate closure models for a range of
9
the solids volume fraction of 0-0.65. For the solids volume fraction between 0-0.3, which is a typical range
10
of the solids volume fraction found in riser, these studies find only minor differences between the calculated
11
values from alternate closure models accept those from closure models proposed by Syamlal et al. (1993).
12
Consequently, the effects of alternate closure models for each of the KTGF parameters (such as the solid
13
pressure, frictional pressure, shear viscosity, granular viscosity, frictional viscosity, and radial distribution
14
function) on prediction of the riser hydrodynamics have not been investigated by conducting riser
15
simulations.
16
3.2.4.
17
3.2.4.1. Inlet boundary conditions
18
Boundary conditions have been used to specify numerical values at the boundary cells in CFD models. The
19
previous riser simulations (Table-2) employ three different types of flow domains such as periodic, 2D and
20
3D geometries. Using the periodic domain (Agrawal et al., 2001; Andrew IV et al., 2005; Benyahia et al.,
21
2009; Igci et al., 2008; Naren and Ranade, 2011), only a differential element of riser is considered and
22
computed values at the outlet boundary surface are passed to the inlet boundary at every time step. This type
23
of periodic flow domain is simulated to evaluate the effect of various closure models and flow conditions on
24
fully-developed radial profiles. In order to evaluate effect of model parameters on both axial and radial
25
profile, 2D and 3D flow domains with appropriate inlet and outlet boundary conditions are required. In these
26
simulations, inlets are configured with velocity-inlet boundary condition.
pt
ed
M
an
us
cr
ip t
1
Ac ce
Effect of boundary conditions
Page 45 of 45
Page 45 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
At gas inlet, velocity of the gas phase equal to operating superficial gas velocity is specified with volume
2
fraction of solids phase equal to zero, while at solid inlet; average solid velocity at inlet is calculated from a
3
circulating solid mass flux by,
eq.(55)
ip t
4
Where, Gs is the solid flux circulating in riser, Ariser is the cross sectional area of riser, and Ainlet is the cross
6
sectional area of solid inlet. This boundary condition also needs a value of solids volume fraction which is
7
generally taken as the value at minimum fluidizing condition (0.5) or maximum packing condition (0.63).
8
Pressure-outlet boundary condition is most commonly used for the outlet. By specifying gauge pressure
9
equal to zero, the pressure at outlet is set as the atmospheric pressure.
an
us
cr
5
Previous CFD studies have been conducted by using both 2D and 3D flow domains. Furthermore, it can also
11
be observed that the previous simulations more choose 2D flow domain, as they are computationally more
12
feasible. Representation of an actual 3D riser using 2D domain then requires appropriate selection of the inlet
13
boundary conditions. Benyahia et al. (2001) have simulated a 2D geometry with one solid inlet (as
14
schematically shown in Figure-11a), mimicking the inlet configurations of the 3D boundary conditions of
15
experimental setup. In Benyahia et al.’s simulation, most of the solid particles are found to be concentrated at
16
the inlet side of the riser causing unsymmetrical radial distribution of the solid density (Figure-11b), which
17
did not agree with the core-annulus profile. Recently, Li et al. (2014a,b) have also observed accumulation of
18
the solids at one side of riser with one side inlet in 2D simulations. Since a 2D domain with one side inlet
19
cannot directly capture the asymmetric entry of solids to a continuously flowing riser, the 2D riser
20
simulations have been conducted using different types of boundary conditions for the inlets and outlets (Pita
21
and Sundaresan, 1993; Neri and Gidaspow, 2000; Benyahia et al., 2001; Yang et al., 2003; Almuttahar and
22
Taghipour, 2008).Benyahia et al. (2001), Yang et al. (2003) and Shah et al. (2012) have used two-sided
23
solids inlets for 2D riser simulations (Figure-12a). Shah et al. (2012) have compared the predictions from the
24
two-sided solids inlets in 2D simulations with both predictions from 3D simulations and experimental data
25
(Figure-11c), and found reasonable qualitative agreements between predictions and experimental data.
Ac ce
pt
ed
M
10
Page 46 of 46
Page 46 of 114
Ac ce
pt
ed
M
(a)
an
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(c) (b) Figure-11: (a) Schematic of 3D and 2D boundary conditions; (b) and (c) radial profile using one side solids inlet and two-sided solids inlets respectively. 1
3.2.4.2. Wall boundary conditions
2
For wall boundary condition, three alternatives (no-slip, free-slip or partial-slip) boundary conditions can be
3
used for an individual phase. Generally, no-slip boundary condition is used for the gas phase; whereas, each
4
of the three alternates have been considered for solids phase by previous studies. The no-slip boundary
5
condition for solids phase is set by equating the tangential and normal velocities of the solids to zero. The
6
partial-slip boundary condition can be configured using correlations developed by Sinclair and Jackson
Page 47 of 47
Page 47 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(Sinclair and Jackson, 1989) for the wall shear and Johnson and Jackson (Johnson and Jackson, 1987) for the
2
granular energy dissipation at wall. The wall shear boundary condition for the solids phase is given by rate of
3
axial momentum transferred to wall by the particles in a thin layer adjacent to wall surface (Sinclair and
4
Jackson, 1989):
eq.(56)
cr
ip t
5
Where,
is the specularity coefficient. Its value equal to zero denotes free slip or specular wall, and unity
7
denotes no-slip wall.
8
The granular energy at wall can be obtained by using equation of Johnson and Jackson (1987). The granular
9
energy flux can be positive (wall as sink) or negative (wall as source) depending upon the relative magnitude
10
of granular energy dissipation due to non-elastic collision between particle and wall, and generation of
11
granular energy due to shear at wall. The granular energy dissipation due to inelastic collision with wall can
12
be written as,
ed
M
an
us
6
eq.(57)
Ac ce
pt
13
14
Thus, the partials-slip condition for the solids phase requires two additional parameters i.e. the specularity
15
and wall restitution coefficients.
16
The selection of the wall boundary conditions for the solids phase dictates velocity and volume fraction near
17
wall, and hence, appropriate wall boundary condition is very critical for capturing the core-annulus radial
18
profile. Table-10 shows that the partial-slip (Johnson and Jackson, 1987) wall boundary condition for the
19
solids phase is overwhelmingly used in previous studies. However, the use of the partial-slip configuration
20
needs a specularity coefficient and particle-wall coefficient of restitution to calculate friction between the
21
wall and particles and energy dissipation in the particle-wall collisions respectively. The specularity
22
coefficient represents a fraction of collisions that transfers momentum to wall; while the coefficient of the
23
restitution represents the fraction of the energy dissipated in the particle-wall collisions. Both these
24
parameters are dependent on material properties particles and wall. The restitution coefficient has been
Page 48 of 48
Page 48 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
measured from controlled experiments by Gorham and Kharaz (2000), Joseph et al. (2001), Kharaz et al.
2
(2001) and Mangwandi et al. (2007). However, values reported in these studies are for different particle sizes
3
and materials, not for FCC catalyst. Hui et al. (1984) and Li and Benyahia (2012) have mentioned that
4
experimental estimation of the specularity coefficient is not feasible. One way to determine the specularity
5
coefficient is to adjust this parameter to fit experimental data. Li and Benyahia (2012) have provided a model
6
for specularity coefficient as a function of the particle-wall restitution coefficient, the frictional coefficient
7
and the normalized slip velocity at the wall. Li and Benyahia’s model for specularity coefficient has not been
8
tested by conducting riser simulation.
ip t
1
cr
9
Benyahia et al. (2012); Almuttahar and Taghipour (2008); Wang et al. (2012) and Shah et al. (2015) have
11
investigated the effect of the specularity coefficient and coefficient of the restitution on flow predictions.
12
While Benyahia et al. (2012); Almuttahar and Taghipour (2008) and Wang et al. (2012) have found only
13
minor effect of the coefficient of the restitution on the predicted axial and radial profiles; the authors found
14
significant impact of the specularity coefficient. Shah et al. (2015) have investigated a wide range of
15
specularity coefficients from 0.1 to 0.0001 in riser with low and high solids flux conditions (Figure-12). As
16
shown in Figure-12(a) and (b), a lower specularity coefficient of 0.0001 gave axial profiles of the solids
17
volume fraction (in the low solids flux riser) and pressure drop (in the high solids flux riser) that reasonably
18
agree with the experimental data. Figure-12(b) and (c) shows the effect of the specularity coefficients on
19
radial profiles at low and high solids flux conditions (Shah et al., 2015). While the specularity coefficient of
20
0.0001 also predicted experimental values in dilute section of the low solids flux riser; in dense bottom
21
region of the low solids flux riser, a higher specularity coefficient of 0.1 gave good agreement between the
22
predictions and experimental data. Shah et al. (2014) have concluded that the impact of the specularity
23
coefficient on the flow predictions has varied with the flow conditions. Therefore, this parameter must be
24
adjusted with the flow conditions.
Ac ce
pt
ed
M
an
us
10
Page 49 of 49
Page 49 of 114
(a)
(b)
High solids flux
pt
ed
M
Low solids flux
an
us
cr
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Ac ce
(b) (c) Figure 12. Effect of wall boundary condition, (a) height of riser Vs. volume fraction of solid, and (b) height of riser Vs. pressure drop, (c) volume fraction of solid Vs. radial position (low solid flux riser, Gs = 14.3 kg/m2s), and (d) volume fraction of solid Vs. radial position (high solid flux riser, Gs = 489 kg/m2s) 1
3.3. Shortcomings
2
Following shortcomings of the E-E based CFD models of riser can be stated.
3
(1) The treatment of the solids phase as a fluid in the E-E approach requires several closure models that are
4
derived by using the KTGF. For each closure law, several models are proposed in literature. A selection
5
of appropriate closure models is critical to achieve close quantitative agreements between the predictions
6
and experimental data. Various combinations of the closure models (see Table-2) have been used in the
Page 50 of 50
Page 50 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
past to get closer agreements with the experimental data at various flow conditions. Clear guidelines for
2
selection of KTGF closures are not available.
3 (2) The dominant closure in the E-E model is the interphase drag. Both conventional models and multiscale
5
drag models have been used. For coarse grid simulations, conventional drag models could not predict
6
experimental observations quantitatively. The multiscale structure dependent drag models such as the
7
EMMS and SGS drag models have shown significant improvements in predictions of axial and radial
8
profiles at different flow conditions (Yang et al., 2003; Benyahia, 2009; Shah et al., 2011b,c; and Shah et
9
al., 2015). However, the SGS model is computationally expensive, and the application of the EMMS
10
drag model needs local values of cluster properties such as size and volume fraction at a given flow
11
conditions. So far, experimental observations on structure of clusters (shape, size and volume fraction) at
12
different flow conditions are not available.
an
us
cr
ip t
4
13
(3) Predicted radial profiles from different simulation studied show reasonable qualitative agreements with
15
experimental data with some quantitative discrepancies. To eliminate such discrepancies, most of
16
previous efforts have been devoted to investigate the gas-solid drag model; while the effect of the lift
17
force has not been investigated.
ed
pt
18
M
14
(4) Both the laminar and turbulent models for viscous stresses in the gas phase have been used in two-fluid
20
models of riser. Interestingly, both these models could give reasonable qualitative agreements with
21
experimental data in previous studies. Shah et al. (2015) showed that different types of the k- turbulent
22
models (dispersed and per phase models) for the gas phase gave similar time-averaged values, but
23
instantaneous flow structures are quite different. Further analysis suggested that turbulent kinetic energy
24
and dissipation rate predicted by the k- dispersed and per phase models are different by orders of
25
magnitude. Without detail time-averaged and instantaneous data on velocities and volume fraction, it is
26
difficult to conclude on suitability of use of the laminar or type of turbulent model.
Ac ce
19
Page 51 of 51
Page 51 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1 (5) As described in Table-2, almost all previous studies have used partial-slip wall boundary conditions,
3
which has been able to capture the core-annulus distributions of the solids phase. But, this wall boundary
4
conditions required empirical constant such as secularity coefficient, experimental value of which is also
5
not available. Therefore, most of the studies used this parameter for fine-tuning of the predictions. Li and
6
Benyahia (2012) have provided a model for specularity coefficient as a function of the particle-wall
7
restitution coefficient, the frictional coefficient and the normalized slip velocity at the wall. Such
8
correlation between the specularity coefficient and flow conditions needs further investigation by
9
conducting riser simulation.
cr
ip t
2
us
10
(6) Extensive simulation works have been conducted using 2D domain of riser due to computational
12
feasibility, while experiments are mostly conducted using 3D cylindrical riser. Furthermore, in the 2D
13
geometry, different types of assumptions such as two-sided inlets for the solids phase and bottom inlet
14
for the gas phase, single bottom inlet for gas and solids phases, and cycle boundary conditions to
15
maintain a constant solid holdup have been used. Shah et al. (2012) have showed that the inlet boundary
16
conditions in 2D domain significantly affect profiles in so-called fully developed region (at top section
17
of riser). Therefore, to avoid uncertainty arising from the assumed 2D boundary conditions, it is
18
necessary to evaluate effect of closure models in exact 3D geometry.
ed
M
an
11
3.4. Eulerian-Lagrangian model
20
In E-L model, the gas phase is represented as a continuous phase for which the averaged mass and
21
momentum balance equations are solved. These equations are similar to those discussed in the E-E model
22
(section-3.1.1.). The particle phase is considered as discrete particle. For individual particle, Newton’s force
23
balance equation is solved. The force balance equation around individual particle can be written as:
25
Ac ce
24
pt
19
(eq. 58)
(eq. 59)
Page 52 of 52
Page 52 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(eq. 60)
Where, xp is the position of particle p; Mp is the mass of particle, vp is the velocity of particle, Fp is the total
3
force on particle, Fd is the drag force on particle, Fc is the contact force on particle, Tp is the total torque on
4
particle, Ip is the moment of inertia of particle, Ω is the angular velocity of particle.
5
The drag force term on a particle residing in kth control volume can be written as (Garg and Dietiker, 2013):
cr
ip t
2
(eq. 61)
an
us
6
Where, Pg,k is the gas pressure in kth control volume, vp is the volume of parcel, β is the drag coefficient, Xp
8
is the location of parcel, ug(Xp) is the velocity of the gas phase at location Xp, and up is the velocity of parcel.
9
In the Eulerian equation for the gas phase, the drag force term can be calculated as:
ed
(eq. 62)
pt
10
M
7
Where, vk is the volume of kth control volume, and np is the total number of particles in kth control volume.
12
Table-5 lists CFD studies of riser using the E-L gas-solid flow model. It can be seen that there are five
13
different types of modelling approaches have been adopted within the E-L framework. These approaches are
14
(i) discrete particle model (DPM), (ii) discrete element model (DEM)-hard sphere, (iii) DEM- soft sphere,
15
(vi) multiphase particle in cell (MP-PIC) and (v) dense discrete particle model (DDPM). Helland et al.
16
(2000) and Mansoori et al. (2005) have adopted DPM approach by considering only the drag force and
17
gravity force in the force balance equation for the particle phase. Zhou et al. (2002), Zhou et al. (2007), He et
18
al. (2009) and Cheng and Jin (2010) have used DEM approach. The DEM approach considers the contact
19
forces due to particle-particle and particle-wall collisions. In this approach, the contact forces can be
20
calculated by either hard-sphere (Allen and Tildesley, 1989) or soft-sphere (Cundall and Stack, 1979) model.
21
In the hard-sphere model, collisions are considered to be binary and instantaneous. The time step is
22
determined by the minimum collision time between any one pair of particles. Consequently, the time step is
23
directly proportional to the mean free path or inversely proportional to the particle volume fraction. Due to
Ac ce
11
Page 53 of 53
Page 53 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
consideration of binary collisions, this approach is more suitable to dilute gas-solid flow (Garg et al., 2012).
2
For a dense gas-solid flow, the hard-sphere approach requires much smaller time-scale. Furthermore, a dense
3
flow is characterized by collisions and persistent contact between multiple particles. In the soft-sphere
4
model, the particles are allowed to overlap each other. The overlap between the two particles is represented
5
as a system of springs and dashpots in both normal and tangential directions. Here, the contact force on a
6
particle is calculated as a sum of conservative force represented by spring and dissipative forces represented
7
by dashpot. To calculate conservative and dissipative forces, stiffness of spring and damping coefficient of
8
dashpot are used. The spring stiffness, damping coefficient and overlap between particles are used to
9
calculate the time step. For detail formulation of the soft-sphere approach, readers are refer to Cundall and
cr
ip t
1
Stack (1979); Tsuji et al. (1993); Deen et al. (2007); Zhu et al. (2008) and Garg et al. (2012).
11
Vegendla et al (2011), Li et al. (2012), Chen et al. (2013), Wang et al. (2014) and Jiang et al. (2014) have
12
used multi-phase particle-in-cell (MP-PIC) method (Andrews andO’Rourke,1996; Snider,2001; Benyahia
13
and Galvin,2010) to simulate the gas-solid flow in riser. In MP-PIC method, particles are represented by
14
parcels where each parcel contains a certain number of real particles of the same diameter and density
15
(Snider, 2001). In implementation of MP-PIC, the effect of several particles forming a parcel is manifested
16
through solids volume fraction in the drag force term. The drag force on parcel residing in kth control volume
17
can be calculated similar to the drag force on each particle. However, the Eulerian drag force term can be
18
calculated as:
ed
M
an
us
10
eq.(63)
pt
19
Where, vk is the volume of kth control volume, NT is the total number of parcels in kth control volume and np is
21
the number of particles in a parcel.
22
The interactions between particles are considered through the solid phase normal stress but not directly
23
through particle contact forces (Li et al. 2012), while particle-wall interactions are modelled by using bounce
24
back rules where the particle-wall restitution factor is used for energy dissipation. In the MP-PIC method, the
25
collision force term is calculated as (Li et al., 2012 and Garg and Dietiker, 2013):
26
Ac ce
20
(eq. 64)
Page 54 of 54
Page 54 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
Where, Ps is the solids pressure which is calculated by,
2
(eq. 65)
Where, ps* and α are model parameters, whose recommended ranges are 1-100 Pa (Garg and Dietiker, 2013)
4
and 2-5 (Auzerais et al., 1988) respectively. Clote et al. (2010) have used slightly different approach where
5
solids pressure term appearing in above equation is calculated by used of the KTGF approach. This approach
6
is called as DDPM.
7
Different E-L modelling approaches can be evaluated in terms of (i) computational feasibility and (ii) their
8
ability to capture physical phenomena. Due to consideration of all real particles in riser, E-L approach
9
without contact forces (mentioned as DPM in Table-5) or DEM approach becomes computationally
10
impractical for real riser system. Consequently, most of DEM studies (see Table-5) have been conducted for
11
significantly less number of particles, scaled-down geometries or large sized particles. On the other hand,
12
MP-PIC approach considers flow of parcels number of which is significantly lower than number of real
13
particles. In addition, the MP-PIC approach does not take particle collision model explicitly and therefore, a
14
longer time step can be adopted for the particle phase. As a result, the MP-PIC approach becomes
15
computationally more feasible for riser simulations.
Ac ce
pt
ed
M
an
us
cr
ip t
3
Page 55 of 55
Page 55 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Table-5: CFD studies of risers using E-L models Modelling approaches (a) (b) (a) (b)
Drag model (a) (a) (c) (a)
Number of particles × 10-3 250 >7
Investigations
>7
(c)
(b)
40
Effect of inlet gas velocity Effect of coefficient of restitution and friction coefficient Effect of operational, particle properties and geometrical parameters on clusters
0.02
0.14 (2D)
100
1400
He et al. (2009) Cloete et al. (2010) Wu et al. (2010) Zhao et al. (2010) Vegendla et al (2011) Li et al. (2012)
0.05 × 0.015 0.076 0.05 0.1 0.1 0.05, 0.09
0.3 0.4 (2D) 2 (2D) 2-6 m 1 2.79, 10.5
335 67 500 500, 520 77 60, 54
2500 1500 1500 950, 2620 1550 1000, 930
(b) (d) (c) (c) (e) (e)
(b) (a) (c) (c) (b) (d)
40-200 -
Chen et al. (2013)
0.09, 0.2, 0.3
54, 76 , 140
930, 1712, 2600
(e)
(b)
-
Wang et al. (2014)
0.4
10.5, 14.2 (2D) 8.35 (3D) 3
(e)
(b)
400, 540, 890
Jiang et al. (2014)
0.42 × 0.92
5.8
200-1000
2620
(e)
(a)
-
Shi et al. (2015) Shi et al. (2015)
0.152 0.15
7.9 3
2550 2222
(e) (e)
(b) (b)
-
Modelling approaches:
cr
us
an
279, 291
150 160
ip t
0.37 (2D)
M
Particle ρs (kg/m3) 2400 2650 1020 2650
Dp (µm) 126 700, 1200 500 700
Solids residence time, effect of flow parameters Effect of number particles in parcel Effect of drag model and particle-wall restitution coefficient Axial and radial volume fraction profiles Distribution of particles in binary and polydispersed solids along height Distribution of solids concentration at inlet and in standpipe Analysis of solids RTDs and back mixing Effect of exit geometry
Drag model: Wen-Yu Gidaspow De Felice EMMS
Ac ce
DPM- Discrete particle model DEM (Discrete element method) - hard sphere DEM-soft sphere DDPM – Dense discrete phase model MP-PIC – Multi phase particle-in-cell
ed
Zhang et al. (2008)
Riser dimensions H (m) 1 0.378
pt
Helland et al. (2000) Zhou et al. (2002) Mansoori et al. (2005) Zhou et al. (2007)
Dt (m) 0.1 0.07 0.0305 0.07
Page 56 of 56
Page 56 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
3.5. Drift flux analysis
2
Various drag force formulations and their relative merits are described in section 3.2.1. While using any one
3
of the drag laws, it is necessary to have a priori knowledge of cluster size distribution (CSD). Obviously, the
4
CSD depends upon the type of particle, size distribution, age of the catalyst, solid loading, column diameter,
5
superficial gas velocity, physical properties of the raw material, etc. Thus the development of procedures for
6
the estimation of the CSD is a formidable task. However, this problem can be addressed (though not with
7
desired rigor) with the help of drift flux model and measurements by the gamma ray densitometry. In a
8
rigorous manner, the DNS has a capability of predicting cluster size distribution and the corresponding drag
9
forces and the subject is discussed in detail in section 4.
cr
ip t
1
In order to address the problem of relation between the CSD and the drag force, drift flux analysis can be
11
applied for the FCC riser for the representation of the combined effect of superficial gas velocity, solids
12
circulating flux, column dimensions and physical properties of gas and solids on the solid holdup and
13
velocity profiles. The drift flux model for two-phase flow, proposed by Zuber and Findlay (1969), correlates
14
the averaged volume fractions of solids with the superficial velocities.
an
us
10
eq. (66)
ed
M
15
pt
16
17
Where,
18
superficial velocity of gas, and c0 and c1 are the drift flux constants. The values of c0 and c1 indicate the
19
quality of radial profile of the solid phase and the slip velocity respectively. The original model of Zuber and
20
Findlay (1969) considers the continuous phase profile to be flat. However, it is known to have strong
21
circulatory flows (upflow in the central region and down flow near the column wall). Under these conditions
22
the drift flux constants take the following forms (Thakre and Joshi, 1999; Joshi, 2001; and Ekambara et al.,
23
2005).
Ac ce
is the time-averaged holdup of solids inside riser, Vs is the superficial velocity of solids, VG is the
24
25
;
eq. (67)
Page 57 of 57 Page 57 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Where Vslip is the relative velocity between two phases,
is the instantaneous volume fraction of the
2
dispersed phase, εg is the instantaneous volume fraction of the continuous phase and ug are the interstitial true
3
velocities of the solid and gas phases velocity of solids.
4
The values of the drift flux constants can be estimated by using experimental data. The cold flow
5
experiments have reported several sets of radial and axial profiles of volume fraction and velocities of the
6
solids phase at different superficial gas velocity and solid fluxes. These data can be used to plot
vs.
cr
ip t
1
to estimate c0 and c1 for a given combination of operating conditions and column dimensions. The
8
coefficient c0 represents the transverse volume fraction profile, while the coefficient c1 represents the slip
9
velocity. CFD predictions of FCC riser must give values of holdup and slip velocity that satisfy the values of
10
c0 and c1 for a given operating condition and column dimensions. As the slip can be adjusted by the drag
11
coefficient, the value of c1 can be used for the estimation of representative slip velocity of drag coefficient of
12
the solid particles (which includes the cluster size distributions) and the gas. This slip velocity can be used
13
for the estimation of drag force.
14
It is further recommended that the gamma ray densitometry be used for the measurement of radial εs profile
15
at 5 to 10 locations in the FCC riser. These be substituted in the CFD simulation for getting profiles of
16
and
17
drift flux constant c0 which enables the estimation of c1 at a given location and hence the axial variation of
18
the CSD and the drag force. It may be emphasized that, the abovementioned procedure is iterative and needs
19
to begin with some reasonable value of drag coefficient.
20
4.
21
In riser flow, clustering is an important phenomenon which includes continuous formation and breakup of
22
clusters and complex interactions between fluid turbulence and dispersed solids phase. The formation of
23
clusters dictates intrinsic hydrodynamics, and therefore, understanding the clustering phenomenon is critical.
24
Due to small spatial and temporal scale of the clusters, experimental investigation to resolve the clustering
25
phenomenon under flow conditions of riser is currently impractical. Consequently, several investigators have
Ac ce
pt
ed
M
an
us
7
at various axial locations. Thus, the knowledge of ug, us, εg and εs can be used for the estimation of
Direct numerical simulations (DNS)
Page 58 of 58 Page 58 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
adopted DNS approach. Two different numerical approaches, namely lattice Boltzmann method (LBM) and
2
immersed boundary method (IBM) have been previously used to conduct DNS of fluid-solid suspensions.
3
The LBM has been used by Hill et al. (2001a,b); van der Hoef et al. (2005); Beetstra et al. (2007); Yin and
4
Sundaresan (2009); and Shah et al. (2013). In the LBM, a flow domain is represented by number of lattices
5
and flow of fluid is calculated by updating velocity distribution at each lattice by using Boltzmann’s velocity
6
distribution function. Flow of particles is resolved by applying Newton’s force balance equation. Force
7
interactions between fluid and particle are then calculated from the velocity distributions at boundary nodes
8
and velocity of particles. The IBM has been used by Uhlmann (2005); Garg et al. (2010) and Tenneti et al.
9
(2011) to study the drag between the gas and solids phases. In the IBM, fluid is represented in an Eulerian
10
framework, whereas particles are represented in a Lagrangian framework. The Eulerian variables are defined
11
on a Cartesian mesh, and the Lagrangian variables are defined on a curvilinear mesh that moves freely
12
through the Cartesian mesh without being constrained to adapt to it in any way at all. The fluid-solid
13
interactions are accounted via a smoothed approximation to the Dirac delta function (Peskin, 2002). For
14
simplicity, previous DNS studies on the clustering can be divided into two categories, i.e. (i) those focus on
15
the impact of the clusters on the fluid-particle drag and (ii) those focus on the interactions between the
16
turbulent and clusters.
17
Hill et al. (2001a,b); van der Hoef et al. (2005); Beetstra et al. (2007); Yin and Sundaresan (2009); and
18
Tenneti et al. (2011) have adopted the DNS approach to investigate the drag force in solid-fluid suspensions
19
by simulating flow passing through a fixed assembly of randomly located mono- or bi-dispersed solid
20
particles inside a 3-D cubical flow domain with periodic boundary conditions. These studies calculated
21
effective drag force in gas-solid suspensions, and then, compared the calculated drag force with that obtained
22
from the conventional drag models such as the Wen-Yu, Syamlal-O’Brien, and Gidaspow models. Such
23
comparison has been made at various solid volume fraction and particle Reynlods numbers. The authors
24
observed wide discrepancies between the drag values from the DNS and the conventional models.
25
Consequently, Hill et al. (2001a,b), Beetstra et al. (2007) and Tenneti et al. (2011) have given a modified
26
drag correlations derived from the DNS results. Beetstra et al. (2006); Zhang et al. (2011); and Shah et al.
27
(2013) have conducted DNS of cluster configurations inside a cubical flow domain. Zhang et al. (2011) and
28
Shah et al. (2013) have observed that the calculated drag for cluster configurations is significantly lower than
29
that from random configurations. The reduction in drag due to the presence of cluster increases with an
30
increase in the solid volume fraction and particle Reynolds number. Main drawbacks of these studies are that
31
(i) cluster properties (sized, shape, number of particles) have been assumed, (ii) clusters have been treated as
32
stationary, and (iii) limited number of particles (32-1000) in a small flow domain (10-20 times particle
33
diameter) have been considered. To remove these limitations, simulations should be performed for a large
34
flow domain with several thousand free-flowing solid particles, which not only interact with surrounding
35
fluid but also experience particle-particle collisions.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 59 of 59 Page 59 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Xu and Subramanian (2010); Derksen (2012); Prevel et al. (2013) and Nicolai et al. (2014) have used DNS
2
approach to study influence of turbulence in fluid phase on the clustering phenomenon. Xu and Subramanian
3
(2010) have used IBM to simulate a flow with particle Reynolds number of 50 and the ratio of the particle
4
diameter to Kolmogorov scale of 5.5. In this study, numerical results for turbulent flow past fixed
5
assemblies, uniform and cluster configurations, of solid particles are analysed. It is found that clustered
6
configurations enhance the level of fluid-phase turbulent kinetic energy (TKE) more than uniform
7
configurations. This increase is attributed to a lower dissipation rate in the clustered particle configuration.
8
The simulations also reveal that the particle-fluid interactions result in anisotropic fluid-phase turbulence due
9
to the anisotropic nature of the interphase TKE transfer and dissipation tensors. Derksen (2012) has used the
10
LBM combined with the IBM to simulate fully coupled solid-liquid flow in a periodic, cubical 3-D domain
11
with homogeneous isotropic turbulence and 5000 spherical particles. In Derksen’s study, particles are made
12
sticky with a tendency to aggregate by a square-well potential, which is a function of a distance of interaction
13
and binding energy. In these simulations, the particles also interact with each other via interacting fluid and
14
through hard-sphere collisions. Derksen (2012) has investigated the impact of both the particle-particle
15
interactions and turbulent intensity on aggregate size distributions. It is found that a small increase in
16
velocity, volume fraction of the solids phase resulted in a single large aggregate with a size comparable to
17
size of the domain. In addition, friction coefficient related to hard-sphere collisions is found to have
18
influence on aggregate size distributions. Prevel et al. (2013) have performed DNS to investigate flow of
19
solids particles under the effect of hairpin vortices in a laminar boundary layer and the preferential
20
accumulation of particles close to wall. In this study, the Navier-Stokes equations are solved for the fluid
21
phase, along with Lagrangian tracking of 200,000 particles. In these simulations, hairpin vortices are
22
generated in a controlled way by a hemisphere protuberance mounted on the lower wall in an initial laminar
23
boundary layer. The authors observed preferential aggregation of outward moving particles in low speed
24
fluid regions between standing vortex and external vortex, and in high shear region between legs of the
25
hairpin vortex. They observe the ejection of the particles from high span wise vorticity region at the head of
26
hairpin vortices. Nicolai et al. (2014) have investigated flow of particles in homogenously sheared turbulent
27
to gain insight into structure and statistical properties of clusters. Nicolai et al. (2014) have performed both
28
experiments and DNS to record particle positions under turbulent flow conditions. In this experiment, a
29
homogeneous shear flow was generated in central part of recirculating water channel, and spherical glass
30
beads are injected in the channel. The instantaneous concentration is measured by imaging the positions of
31
particles using high resolution CCD camera. Both the experimental and numerical datasets are further
32
analysed by statistical methods such as the voronoi decomposition method and the Angular Distribution
33
Function (ADF). The particle snapshots indicate that particles tend to aggregate into thin clusters outside
34
vortex cores, and that the mean velocity gradient induces an evident preferential orientation on the particle
35
clusters. The voronoi analysis gave characteristic dimension of regions which may be identified as clusters
36
and voids. The use of the ADF approach gives further characterization of geometry of the clusters by
37
estimating anisotropy content of particle files at changing the scale of observations. Recently, Joshi and
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 60 of 60 Page 60 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Nandkumar (2015) have made several suggestions for future work which will enable additional
2
understanding in the mechanism of formation/break-up of clusters using DNS simulations.
3
The subject of transition from homogeneous to heterogeneous dispersions have been addressed by many
4
previous studies (Anderson and Jackson, 1967; Didwania and Homsy, 1982; Batchelor, 1988; Thorat et al.,
5
1998; Joshi et al., 2001; Sundaresan, 2003; Thorat and Joshi, 2004; Ghatage et al., 2014). All these studies
6
use the theory of linear stability which consists of the following steps (1) the starting point is one
7
dimensional governing equations of continuity and motion for multiphase system (2) use Reynolds averaging
8
procedure for turbulent dispersion and introduce dispersion coefficients for individual phases (3) eliminate
9
pressure term with the assumption that it is shared by the multiple phases (4) linearization (5) introduction of
10
purturbation variables (6) simplification of all the equations in the form of single equation having one
11
dependent purturbation variable (continuous phase hold-up) (7) application of the theory of linear stability.
12
Joshi et al (2001) have compared all the experimental data on transition in gas fluidized bed and have shown
13
very good agreement with the predictions of the theory of linear stability. Further, there has been a long time
14
need for understanding the origin of disturbance. In this context of transitions, Derksen and Sundaresan,
15
(2007) have performed DNS simulations and established the reasons for initiation of clustering. This
16
pioneering work needs to be taken forward for understanding transitions in FCC risers. The analysis needs to
17
be performed separately (1) for the entry region of complex geometry and (ii) the pneumatic transport region.
18
In addition to the modulation of turbulence, the effect of turbulence on particle aggregation has been
19
analysed (Xu and Subramaniam, 2010; Derksen, 2012). Accordingly, the core-annulus phenomenon in the
20
FCC has been explained. Wei et al. (2013), Shnip et al. (1992) and Derksen and Sundaresan (2007) have
21
made useful contributions on the prediction of onset of transition to heterogeneity from homogeneity.
22
Further, such a transition has been shown to follow enhancement in particle setting velocity due to possible
23
particle aggregation/generation of strong convection currents in the continuous phase. However, like
24
turbulence modulation, the phenomena of aggregation as well as the onset of transition have not been
25
quantitatively understood and substantial additional work is needed.
26
5.
27
CFD studies of a reactive flow in the FCC riser (listed in Table-9) can be classified in two categories, i.e. (i)
28
those who ignore droplet phase by assuming instantaneous vaporization of feedstock (Theologos and
29
Markatos, 1993; Theologos et al., 1997; Gao et al., 1999; Benyahia et al., 2003; Das Sharma et al., 2006;
30
Lopes et al., 2011; ad Li et al., 2013), and (ii) those who consider the flow of feed droplets and vaporization
31
(Gao et al., 2004; Mortignoni and Lasa, 2001; Nayak et al., 2005; Behjat et al., 2011; and Chang et al.,
32
2012). The studies in the first category use the E-E gas-solid flow models coupled with cracking kinetics
33
models; while those in the second category use the E-E gas-solid flow model coupled with droplet phase
34
model (either Eulerian or Lagrangian), droplet vaporization model, and cracking kinetics model. Simulation
Ac ce
pt
ed
M
an
us
cr
ip t
1
CFD models of a reactive flow in FCC riser
Page 61 of 61 Page 61 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
predictions of these models are often compared with the FCC riser performance data (conversion and yields)
2
available in literature.
3
Table-6
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 62 of 62 Page 62 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Table-6: CFD studies of a reactive flow in the FCC riser Riser Modelling approach geometry
4
Cracking kinetics
Momentum transfer
Mass transfer
Heat transfer
Molar expansion
Investigations
Capability of their CFD model Effect of number of feed nozzles Effect of droplet size Effect of water injection at 10m height Droplet concentration near feed inlet zone Effect of droplet size Capability of their CFD model Capability of their CFD model Radial distributions at various heights Effect of droplet size and CTO ratio Capability of their CFD model Performance of two stage riser system; effect of operating conditions Comparison between riser and downer reactors Comparison between phenomenological model with heterogeneous reactive 1-D flow model Profiles of catalyst volume fraction near entry section of riser Feed entry with a mixing section in two stage riser systems Effect of operating conditions and droplet size Effect of operating conditions Comparison between predictions from rotating fluidized bed and conventional riser Effect of catalyst size distributions in FCC riser Effect of reaction temperature and effect of CTO Effect of feedstock and catalyst temperature, effect of CTO
3-D 3-D 3-D 3-D
E-E E-E E-E E-E
Lump kinetics model 3 10 3 13
Deactivation (a) (b) (a) (b)
(a) (a) (a) (a)
NC NC NC NC
(a)+(d) (a)+(d)+(e) (a)+(d)+(e) (a)
NC NC NC Ideal gas
Gao et al. (2001)
3-D
E-E-E
13
(b)
(a)+(b)
Vaporization
(a)+(c)+(d)
Ideal gas
Benyahia et al. (2003) Chang and Zhou (2003) Das et al. (2003)
2-D 3-D 3-D
E-E E-E-E E-E
3 4 12
(a)
(a) (a) (a)
NC Vaporization NC
NC (a)+(b)+(c)+(d) (f)
NC Ideal gas NC
Nayak et al. (2005)
3-D
E-L-L
4, 10
(b)
(a)+(b)
Vaporization
(a)+(b)+(c)+(d)
Ideal gas
Souza et al. (2006) Lan et al. (2009)
2-D 2-D
Well mixed single phase E-E
6 14
(a)
NC
(a)+(d) (a)+(d)
NC Ideal gas
Wu et al. (2010)
2-D
E-L
4
(a)
(a)
NC
(a) + (d)
Ideal gas
Zhu et al. (2011)
1-D
E-E
4
(b)
(a)
NC
(a)+(d)
Ideal gas
Lopes et al. (2011)
3-D
E-E
4
(b)
(a)
NC
(a)+(d)
Ideal gas
Gan et al. (2011)
3-D
E-E
11
(b)
(a)
NC
(a)+(d)
Ideal gas
Behjat et al. (2011)
3-D
E-E-L
4
(a)
(a)+(b)+(c)
Vaporization
(a)+(b)+(c)+(d)
Ideal gas
Chang et al. (2012) Trujillo and De Wilde (2012)
3-D Rotating fluidized bed
E-E-E E-E
12 10
(a) (b)
(a)+(b)+(c) (a)
Vaporization NC
(a)+(b)+(c)+(d) (a)+(d)
Ideal gas Ideal gas
Li et al. (2013)
3-D
E-E
14
(b)
(a)
NC
(a)
Ideal gas
Chang et al. (2014)
3-D
E-E
9
(a)
NC
(a)
Ideal gas
Alvarez-Castro et al. (2015a, b)
3-D
E-E
12
(a)
NC
(a)
Ideal gas
cr
us an
(b) (b)
M
ed pt
Ac ce
(b)
ip t
Theologos and Markatos (1993) Theologos et al. (1997) Theologos et al. (1999) Gao et al. (1999a,b)
(b)
Modelling approach:
Deactivation:
Heat transfer:
E-E: Gas and solids phases E-E-E: Gas, solids and liquid phases E-E-L: Gas and solids as continua, droplets as a discrete phase E-L-L: Gas as a continuum, solids and droplets as discrete phases E-L: Gas as continuum and solids as discrete phase Well mixed single phase
Based on residence time of catalyst Based on coke concentration
Gas-solid heat transfer resistance Gas-liquid heat transfer resistance Solid-droplet heat transfer Heat of reactions Effect of droplet vaporization on enthalpy of hydrocarbon gases Imposed temperature profiles
Momentum transfer: Gas-solid drag Gas-droplet drag Solid-droplet collisions
Page 63 of 63 Page 63 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
5.1. FCC riser performance
2
Paraskos et al. (1976), Shah et al. (1977), Corella et al. (1986), and Bollas et al. (2002) have conducted pilot-
3
scale experiments of FCC plant; while Theologos et al. (1997), Ali et al. (1997), Derouin et al. (1997), Gao
4
et al. (2001), and Yang et al. (2009) have reported commercial FCC plant data. The available FCC plant data
5
can also be divided in two categories, one for which variation of conversion, yields and/or temperature are
6
reported along height of riser (Shah et al., 1977; Derouin et al., 1977) and the other for which the data only at
7
the riser exit are available (Paraskos et al., 1976; Corella et al., 1986; Ali et al., 1997; Bollas et al., 2002;
8
Yang et al., 2009).
M
an
us
cr
ip t
1
ed
(a) (b) Figure-13: (a) Conversion and yield Vs. dimensionless height; (b) Gas oil conversion Vs. gasoline yield (Shah et al., 1977: height = 30 m, diameter = 0.85 m, C/O = 7.06, pressure = 3.04 atm, catalyst dia = 68 μm, feed oil flow rate = 11.87 kg/s, catalyst flow rate = 83.09 kg/s; Derouin et al., 1997: height = 30 m, Diameter = 0.7-1 m, C/O = 5.5, pressure = 3.15atm, catalyst dia = 60-70μm, feed oil flow rate = 55-109 kg/s, catalyst flow rate = 300-600 kg/s) Figure-13(a) shows conversion of gas oil and yield of gasoline along height given in Shah et al. (1977) and
10
Derouin et al. (1997). Both the data sets show a steep rise in the conversion and yield in initial 20-40% of the
11
riser height, and then the profiles plateau. However, quantitative values in these two data sets are quite
12
different. Shah et al. (1977) have shown rise in conversion up to 45% and yield up to 67%; while Derouin et
13
al. (1997) have reported rise in conversion up to 72% and yield up to 48%. The discrepancies between the
14
two data sets, despite similar riser geometries, can be attributed to flow conditions. The feed and catalyst
15
flow rates in Derouin et al.’s riser system is almost 5-10 times higher than those used by Shah et al. (1977).
16
Furthermore, the Catalyst to oil ratio (CTO) in Shah et al.’s experiments is 7 compared to 5.5 in Derouin et
17
al.’s riser system. Figure-13(b) shows the gasoline yield against the gas oil conversion from the two data
18
sets. A significantly higher selectivity of the gasoline was reported by Shah et al. (1977). This suggests that
19
both higher CTO and lower flow rates of both feed and catalyst result in higher selectivity of the gasoline.
20
Paraskos et al. (1976) have conducted pilot-scale FCC riser experiments at isothermal condition, and
21
reported 14 data sets for varying CTO (6.82-8.6) and flowing space time (1.12-59.80). They report the
22
conversion at exit in range of 40% - 78% and yield in range of 31.3% - 53%. Ali et al. (1997) have reported
23
commercial FCC plant data, in which conversion of gas oil, and yields of gasoline, gas and coke components
Ac ce
pt
9
Page 64 of 64 Page 64 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
at exit are 73, 44, 23.45, 5.8 wt % respectively. Furthermore, Ali et al. (1997) also report an exit temperature
2
of 795 K. Gao et al. (2001) have reported yields of five-lumps at two exit temperatures in a commercial FCC
3
riser systems. In Gao et al.’s data, yields of heavy fuel oil, light fuel oil, gasoline, gas and coke are 42, 18.5,
4
27.9, 7.7 and 3.9 wt % at 504 ˚C and 404.4, 18.8, 29.7, 7.8 and 3.3 wt % at 496 ˚C.
5
5.2. Coupling of gas-solid flow model with droplet flow and vaporization models
6
In FCC riser, feedstock is injected through atomizing nozzles in the form of micron-sized droplets. After
7
entering, the feed droplets vaporize upon making contact with hot catalyst and steam. The vaporisation of
8
feedstock controls distribution of temperatures and concentration of gas oil in bottom section of riser.
9
Therefore, it is critical to model the droplet phase along with the gas-solid flow. In previous studies, the
10
droplet phase was modelled by using either Eulerian (Gao et al., 2001; Mortignoni and Lasa, 2001; and
11
Chang et al., 2012) or Langrangian approach (Nayak et al., 2005 and Behjat et al. 2011).
12
Eulerian approach for the droplet phase: If the Eulerian approach is used, the E-E gas-solid model is
13
modified with a third set of mass, momentum and energy balance equations for the droplet phase. The gas-
14
droplet and the solids-droplet momentum and heat transfer are modelled by using correlations for the
15
interphase drag (Table-4) and heat exchange (Table-3). Gao et al. (2001) and Chang et al. (2012) used
16
following mass balance equations for the three-phase flow.
19
cr
us
an
M
ed pt
18
Ac ce
17
ip t
1
eq.(68)
eq.(69)
eq.(70)
20
Where, d denotes the droplet phase; and mdg (
) is the mass transfer from the droplet phase to gas
21
phase due to droplet vaporization. Gao et al. (2001) and Chang et al. (2012) used following momentum
22
balance equations for the three-phase flow.
Page 65 of 65 Page 65 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
eq.(71)
2
eq.(72)
ip t
1
eq.(73)
Here, FD,sg (=
5
gas, droplet-gas, solid-droplet respectively. Chang et al. (2012) used the Syamlal-O’Brien drag model,
6
whereas Gao et al. (2001) used the Wen-Yu drag models to calculate the drag forces.
7
The energy balance equations for the three phases can be written as,
10
) are the drag force exchanges between the solid-
M
ed pt Ac ce
9
) and FD,sd (=
an
4
8
), FD,dg (=
us
cr
3
eq.(74)
eq.(75)
eq.(76)
11
Where, ΔQgs(= ̶ ΔQsg), ΔQdg (= ̶ ΔQgd) and ΔQds (= ̶ ΔQsd) are heat exchange terms between the gas and
12
catalyst phases and the gas and droplet phases respectively. Sgh is the energy source term in the gas phase
13
energy balance equation due to the endothermic heat of reaction. The equations of ΔQgs (= ̶ ΔQsg), are given
14
in the section 3.1.1.5. The heat exchange between the gas and droplet phase can be given as,
Page 66 of 66 Page 66 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
eq.(77)
Where, hgd (= hdg) is the heat transfer coefficient; Adi is the interfacial area between the gas and solids phases;
3
Tg is the temperature of the gas phase; Td is the temperature of the droplet phase and hfg is the heat of
4
vaporization of the droplet phase. The heat transfer coefficient can be correlated to the Nusselt number,
5
which is further correlated by the droplet Reynolds number and Prandtl number according Ranz and
6
Marshall (1952). Chang et al. (2012) have used correlation of Gunn (1978) to calculate the Nusselt number;
7
while Gao et al. (2001) and Theologos and Markatos (1993) have used a correlation of Kothari (1967) for the
8
particle Nusselt number.
9
Lagrangian approach for the droplet phase: In the Lagrangian approach (Nayak et al., 2005 and Behjat et al.
10
2011), force and energy balance equations for each droplet are solved separately. The force balance equation
11
on an individual droplet can be written as,
an
us
cr
ip t
2
eq.(78)
M
12
Where, md is the mass of the droplet, ud is the velocity of the droplet, t is the flow time, FD is the drag force
14
on the droplet, Fg is the gravity force, and Fothers are acceleration forces such as lift force, virtual mass force,
15
pressure force, etc. For each droplet, the position and velocity is updated at every time step.
16
Nayak et al. (2005) have considered forces on the droplet due to continuous phase pressure gradient.
17
Furthermore, Nayak et al. (2005) have used drag model of Morsi and Alexander (1972). Nayak et al.’s
18
equation for force balance is,
pt
Ac ce
19
ed
13
eq.(79)
20
Where, CD is the drag coefficient, vd is the volume of the droplet, ρcont is the density of continuous phase
21
(only the gas phase was considered as continuous phase by Nayak et al., 2005 and Behjat et al., 2011), ucont is
22
the velocity of the continuous phase, and ρd is the density of the droplet phase. As tracking individual droplet
23
is computationally infeasible, Nayak et al. (2005) and Behjat et al. (2011) have considered flow of packets of
24
several droplets. Then, an equivalent diameter of the packet is used to calculate interphase drag and heat
25
exchange terms.
Page 67 of 67 Page 67 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Regardless of approach used for the droplet phase, heat exchange between the continuous phase and the
2
droplet phase causes variation in droplet size and temperature, which is calculated by using a suitable droplet
3
vaporization model. Available droplet vaporization models can be classified in two categories, namely
4
homogeneous where droplets are assumed to receive heat only from the surrounding hot gases (Abramzon
5
and Sirignano, 1989; Renksizbulut and Bussmann, 1993; Buchanan, 1994; Miller et al., 1998; Abramzon and
6
Sazhin, 2006) and heterogeneous which involves the effect of direct collisions between droplet and hot solid
7
particles (Buchanan, 1994; Nayak et al., 2005). Recently, Nguyen et al. (2015) have presented a detailed
8
review of all droplet vaporization models. In this review, the models proposed by Ranz and Marshall (1952),
9
Buchanan (1994) and Nayak et al. (2005), which have been previously used in the FCC riser simulations,
ip t
1
have been compared.
11
The feed droplets are generally injected at a temperature lower than the vaporization temperature; as a result,
12
the vaporization of droplets occurs in three steps, namely inert heating, vaporization and boiling. For inert
13
heating step, the governing equation for heat transfer can be given as,
an
us
cr
10
eq.(80)
M
14
And the heat transfer coefficient for the inert heating step can be calculated by following equations.
16
According to Ranz and Marshall (1952) and Buchanan (1994),
18
19
20
21
pt Ac ce
17
ed
15
eq.(81)
And according to Nayak et al. (2005),
eq.(82)
For vaporization step, the equations for mass transfer can be written as,
eq.(83)
Page 68 of 68 Page 68 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
Where, K is the mass transfer coefficient; Ci,s and Ci, cont are the concentrations of vaporizing component at
2
surface of the droplet and in the bulk of the continuous phase respectively; and MWi is the molecular weight
3
of vaporizing component. According to Ranz and Marshall and Nayak et al.’s models, the mass transfer
4
coefficient can be calculated by,
cr
For the vaporization step, the equations for heat transfer can be written as,
7
And heat transfer for the vaporization step can be calculated by,
9
(Ranz and Marshall, 1952)
ed
M
an
8
us
6
eq.(84)
ip t
5
(Nayak et al. (2005)
Notably, Buchanan (1994) only considered the inert heating and boiling steps.
12
For the boiling step, equations can be written as,
Ac ce
11
13
(Ranz and Marshall, 1952; and Buchanan, 1994)
14
15
eq.(86)
eq.(87)
pt
10
eq.(85)
(Nayak et al., 2005)
eq.(88)
eq.(89)
In eq.(57), Nu can be calculated by,
Page 69 of 69 Page 69 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(Ranz and Marshall, 1952)
3
with
(Buchanan, 1994) – Heterogeneous
In eq.(58), Nueff can be calculated by,
an
4
(Buchanan, 1994) – Homogeneous
eq.(91)
ip t
with
cr
2
eq.(90)
eq.(92)
us
1
(Nayak et al. (2005)
M
5
eq.(93)
Buchanan (1994) has considered two hypothetical limiting cases for heterogeneous vaporization, i.e.
7
infinitely fast heat transfer during collision (limiting case-1) and hard-sphere collision (limiting case-2). In
8
the hard-sphere collision vaporization, Buchanan (1994) has considered elastic collision between droplets
9
and catalyst particles; while in the infinitely fast case, Buchanan (1994) has considered instantaneous transfer
10
of all possible heat from catalyst to droplets. The model equations for Buchanan’s heterogeneous case
11
mentioned above are for the hard-sphere collision model. For Buchanan’s infinitely fast case, the governing
12
equations can be given as,
14
pt
Ac ce
13
ed
6
for the inert heating step
eq.(94)
for the boiling step
eq.(95)
15
The model equations for inert heating and boiling steps are solved using finite element method with a time-
16
step size of 1 × 10-5 s to calculate vaporization times for the flow conditions given by Buchanan (1994). The Page 70 of 70 Page 70 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
calculations are conducted for the Buchanan-homogeneous, Buchanan-heterogeneous (case-1 and case-2),
2
Ranz and Marshall and Nayak et al. models. The comparison of the calculated vaporization times is shown in
3
Figure-14. To show the fidelity of the calculation method, the calculated values of this study for the
4
Buchanan models are also compared with those reported by Buchanan (1994). Good agreement between the
5
calculations of this study and those of Buchanan (1994) is achieved. Figure-14(a) shows that the inert heating
6
times for the Buchanan-homogeneous and Ranz and Marshall models are the same, but the vaporization time
7
of the Buchanan-homogeneous model is approximately three times longer than that of the Ranz and Marshall
8
model. This can be attributed to the significant reduction in the Nusselt number caused by a factor
9
of
ip t
1
cr
, whose value is 1.94 at the simulated flow conditions. This means that the Nusselt
number of Buchanan-homogeneous model is almost half of that in the Ranz and Marshall model. Figure-
11
14(b) shows the comparison of calculated vaporization times from the heterogeneous models. The
12
Buchanan-case-1 (infinitely fast vaporization) gives total vaporization time of 0.9 milliseconds; while the
13
Buchanan-case-2 (hard-collision model) yields vaporization time of 17 milliseconds. The Nayak et al. model
14
yields intermediate value of 5.7 milliseconds. The decrease in the droplet vaporization time using the
15
Buchanan case-2 can be attributed to the use of modified droplet Reynolds number equation (see eq.61),
16
which consists of particle density and particle volume fraction; and hence, the use of the modified droplet
17
Reynolds number results in significant higher heat transfer coefficient. In the calculation for the Nayak et al.
18
model, the Nusselt number equation is modified by using particle Reynolds number (see eq.56), and this also
19
resulted in the higher value for the heat transfer coefficient. As a result, the Nayak et al. model also gives
20
significantly shorter duration for the droplet vaporization. It is noteworthy that the Nayak et al. model has an
21
empirical constant, whose value, as recommended by Nayak et al. (2015), is taken to be 14 in our
22
calculations. But this empirical constant brings further uncertainty to Nayak et al vaporization model.
Ac ce
pt
ed
M
an
us
10
Figure-14: Droplet diameter Vs. time (a) homogeneous models (b) heterogeneous models (flow conditions of
Page 71 of 71 Page 71 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Buchanan et al., 1994: droplet dia. = 300 μm, temp of cont. phase = 811 K, temp of droplet = 700 K, gasdroplet slip velocity = 6.1 m/s, gas-solid slip velocity = 6.1 m/s)
The calculated vaporization times from the simulated models are in milliseconds, and this can be considered
2
as a valid reason for the instantaneous vaporization assumption by several CFD studies. However, the above
3
calculation is for single droplet under a particular flow conditions. In reality, multiple droplets and variation
4
in flow conditions would have significant impact on the droplet vaporization time. The comparison of the
5
models shows wide discrepancies in the computed vaporization times (Figure-14). In previous simulation
6
studies (Table-6), Behjat et al. (2011a, b) and Chang et al. (2012) have used homogeneous vaporization
7
model; while Nayak et al. (2005) have compared their own model with Buchanan’s heterogeneous
8
vaporization model. None of these studies have validated their predictions of droplet vaporization time or
9
height, as experimental data on droplet vaporization height or time in FCC riser are not available in
10
literature. Few studies have reported experimental data of homogeneous droplet vaporization under
11
controlled conditions. Experimental data are available for homogeneous vaporization of water (Ranz and
12
Marshall, 1952), decane (Wong and Lin, 1992) and heptane (Nomura et al., 1996). Nguyen et al. (2015) have
13
compared calculated values from the Ranz and Marshall model with the experimental data. Nguyen et al.
14
(2015) have found that the Ranz and Marshall model could predict the data for vaporization of water but it
15
significantly under predicts the vaporization time for hydrocarbon droplets. As experimental data on
16
homogeneous or heterogeneous vaporization of FCC feed droplets are not available, comparisons of
17
calculated values with those of Buchanan (1994) are considered in this section.
18
5.3. Coupling between gas-solid flow model and cracking kinetic model
19
Catalytic cracking involves conversion of heavier hydrocarbons in the vaporized feedstock to lighter
20
products in the presence of hot zeolite catalyst. The vaporized feed contains thousands of species and hence
21
the cracking results into a similar number of cracked products. The cracking reactions also produce coke and
22
heavy hydrocarbons that deposit on the surface of catalyst reducing its activity. Describing the kinetics of
23
such a complex process in entirety is rather difficult task. Therefore, traditionally, a lump kinetic approach
24
has been adopted to describe the kinetics of catalytic cracking.
25
Weekman Jr (1968) and Wojciechowski (1968) have proposed the first lumped kinetic model, a 3-lump
26
model, which divides reaction mixture in gas oil, gasoline and light gases-coke components. Weekman Jr
27
(1968) has also estimated rate constants and deactivation coefficient by using experimental data. Further,
28
Weekman Jr and Nace (1970) have estimated kinetic parameters of the 3-lump model using moving bed
29
experiments. They also compare the kinetic parameters from fixed bed and moving bed experiments. All
30
these estimates of the kinetic parameters are based on only one type of feedstock. Nace et al. (1971) have
31
reported range of 3-lump kinetic parameters based on experiments with different types of feedstock. Voltz et
32
al. (1971) have given correlations to calculate variations in rate constants with the concentration of aromatic
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 72 of 72 Page 72 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
and naphthenic components of feedstock. Paraskos et al. (1976) and Shah et al. (1977) and Corella et al.
2
(1986) have conducted pilot-scale riser experiments to estimate 3-lump kinetic parameters. In these three
3
pilot-scale experimental studies, residence times of the catalyst are significantly lower than those in the fixed
4
bed experiments of Nace et al. (1971). Furthermore, Paraskos et al. (1976), Shah et al. (1977) and Corella et
5
al. (1986) have used a governing (component balance) equation for the plug flow to calculate yields and
6
conversion. Due to different equations and operating conditions, the estimated 3-lump kinetic parameters
7
widely differ from each other. Lee et al. (1989) have proposed a 4-lump model and estimated kinetic
8
parameters by using micro-activity test (MAT) experiment data. Lee et al. (1989) have also investigated the
9
effect of temperature on the kinetic parameters, and consequently gave values of frequency factors and
10
activation energies for the 4-lump model. Gianetto et al. (1994), Pitault et al. (1994), Ancheyta-Juárez et al.
11
(1997), Abul-Hamayel (2003) have proposed different sets of the 4-lump kinetic parameters. Ancheyta-
12
Juárez et al. (1997) have further extended the work and provide rate constants for the 5-lump model. Bollas
13
et al. (2007) have also given kinetic parameters for the 5-lump model, but Bollas et al.’s model considers
14
additional cracking reaction of the LPG lump to the dry gas which is not considered by Ancheyta-Juárez et
15
al. (1997). Their model also considers tc–n type catalyst deactivation. To remove the dependence of the
16
kinetic parameters on charged feedstock, Jacobs et al. (1976) have proposed a 10-lump scheme; in which the
17
gas oil lump is further divided into heavy fuel oil consisting of paraffinic, naphthenic and aromatic lumps.
18
Jacobs et al. (1976) also use reaction rate equation different from previous 3-lump studies. Apart from the
19
lumped kinetic models, cracking kinetics models based on single-event kinetics (Froment, 1992; Feng et al.,
20
1993), molecular lump approach (Pitault et al., 1994), neural networks (McGreavy et al., 1994), and pseudo-
21
component approach (Liguras and Allen, 1989b, Gupta et al., 2007) are also available in the literature.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 73 of 73 Page 73 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Hari et al. (1995) 582.84
0.97-2.87 10.3-39.2 18.5-40.1
23.76 684 5688
15.12 845.28 2412
Ej (kJ/mol) 68.2495 89.2164 64.5750 52.7184 115.4580 72.2526 51.6726 117.7050
Abul-Hamayel (2003) Ej (kJ/mol) 121.-743 117.208 35-158 79.534 21-60 62.79 12-104 142.324 3-20 125.58
Ancheyta-Juarez et al. (1997) Ej (kJ/mol) 699 57.3482 125.28 52.7436 50.4 31.8136 33.48 65.7202 0.000072 66.5574
92-281
315
cr
(Pre-exp. Factor) 7.978×105 4.549×106 3.765×104 3.255×103 79.57 1.937×106 4.308×103 3.017×108
us
k12 (hr-1) (wt%-1) k13 (hr-1) (wt%-1) k14 (hr-1) (wt%-1) k23 (hr-1) k24 (hr-1) k1 (hr-1) (wt%-1) k2 (hr-1) α (hr-1)
Corella et al. (1985) 432
75.348
an
VGOGasoline VGOGas VGOCoke GasolineGas GasolineCoke Reacn 1+2+3 Reacn 3+4 Deactivation (ϕ=e-αtc)
Nace et al. (1970) 7.7-33.5
ip t
Table-7: 3-, 4-, and 5- lump kinetic parameters derived from controlled fixed bed experiments. 3-lump kinetic parameters Weekman (1968) VGOGasoline k12 (hr-1) (wt%-1) VGOGas + coke k13 (hr-1) (wt%-1) k23 (hr-1) GasolineGas + coke k1 (hr-1) (wt%-1) 143 and 194 Reacn 1+2 -αtc 18.8 and 21.8 Deactivation (ϕ=e ) α (hr-1) 4-lump kinetic parameters Lee et al. (1989)
5-lump kinetic parameters
ed
699 128.52 0.36 50.4 21.96 11.52 0.000072 7.2 315
pt
22
k12 (hr-1) (wt%-1) k13 (hr-1) (wt%-1) k14 (hr-1) (wt%-1) k15 (hr-1) (wt%-1) k23 (hr-1) k24 (hr-1) k25 (hr-1) K34 (hr-1) α (hr-1)
Ej (kJ/mol) 57.3208 52.3 49.3712 31.7984 73.22 45.1872 66.5256 39.748
Ac ce
VGOGasoline VGOLPG VGODry gas VGOCoke GasolineLPG GasolineDry gas GasolineCoke LPGDry gas Deactivation (ϕ=e-αtc)
M
Ancheyta-Juarez et al. (1997)
Page 74 of 74 Page 74 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Of all the kinetic models proposed, the lump kinetic models are most appropriate for CFD simulations; and
2
hence, they have been overwhelmingly used by previous CFD studies. However, selection of lump model
3
and its kinetic parameters is a tricky task as several sets of kinetic parameters for a given lump model are
4
reported in literature. These kinetic parameters are estimated from experiments using different feedstock,
5
catalyst, reactor type, and operating conditions. The kinetic parameters of 3, 4 and 5-lump models are
6
compared in Table-7. Notably, these listed parameters have been derived only from fixed-bed experiments
7
and similar operating conditions of temperature and CTO. Despite similar experimental conditions for their
8
estimation, large differences between the values of these kinetic parameters can be seen. The differences in
9
the kinetic parameters can be attributed to not only type of feedstock and catalyst properties but also
cr
governing component mass balance equations used in their estimation.
Ac ce
pt
ed
M
an
us
10
ip t
1
Figure-15: Model prediction (of gasoil conversion wt%) Vs. Experimental data of Nace et al. (1971) 11
Figure-15 shows a parity plot for calculated predictions using listed kinetic parameters of the 3-, 4- and 5-
12
lump models against the experimental data of Nace et al. (1971), who have reported gas oil conversion vs.
13
space time for 16 different feedstock at two different run times (1.25 and 5 min). The predictions (shown in
14
Figure-15) are calculated by using the equation for gas oil conversion given by Lee et al. (1989). As
15
expected, the predictions from the 3-lump parameters of Nace et al. (1971) are very close to the experimental
16
values. The 4-lump kinetic parameters of Lee et al. (1989) give predictions within ±20% errors, while the
17
predictions from all other sets of kinetic parameters show wide discrepancy. The error in predictions using
18
parameters of Weekman Jr (1968) and Lee et al. (1989) is due to the fact that these parameters are derived
19
only for one type of feedstock while the experimental data of Nace et al. (1971) is for different types of
20
feedstock. In addition, the parameters of Weekman Jr (1968) are derived from conversion and yield data vs.
Page 75 of 75
Page 75 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
liquid hourly space velocity, while experimental data of Nace et al. (1971) is available against weight hourly
2
space velocity. The discrepancy in using the kinetic parameters of Hari et al. (1995), Abul-Hamayel (2003)
3
and Ancheyta-Juárez et al. (1997) can be attributed to their use of governing equation that does not account
4
for the effect of space velocity sufficiently. In CFD simulation of riser, each cell has different flow
5
conditions; while the available kinetic parameters are only valid for range of operating conditions. The above
6
analysis shows that significant uncertainty in CFD predictions is associated by direct use of the kinetic
7
parameters of previous studies. Therefore, one should review of all available kinetic data its validity over a
8
range of operating conditions of interest.
9
In a CFD model, number of species in the gas phase depends on the selection of lump kinetic model. Change
10
in mass fraction of each species (lump) is calculated by using a separate component balance equation that
11
can be written as::
an
us
cr
ip t
1
eq.(96)
M
12
Where, yj is the mass fraction of the jth species; and rj is the net rate of the production of the jth species. The rj
14
can be calculated by:
ed
13
eq.(97)
pt
15
Where, Mw,i is the molecular weight of the jth species; k is the total number of cracking reactions; and ri,j is
17
the rate of production of the jth species in the ith reaction. In a generalised form, the rate of consumption of a
18
reactant jthlump in ith reaction can be given by
19
Ac ce
16
eq.(98)
20
Where, kij is the rate constant of consumption of jth lump in ith reaction; cj is the concentration of jth lump, n is
21
the order of the reaction; and
22
VGO cracking reactions and zero for all other cracking reactions.
is catalyst deactivation function. The value of n is assumed to be one for the
Page 76 of 76
Page 76 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
There are two other types of rate equations are available in the literature, which can be written as:
2
(for 3-lump model of Pachovsky and Wojciechowski, 1971)
eq.(99)
Where, n is one for the VGO cracking reactions; n is zero for all other cracking reactions, cj0 is initial
4
concentration of pure lump j. The value of n for the VGO cracking implies that the cracking of the VGO is
5
easier for the initially vaporized fraction and it becomes progressively difficult as the VGO travels upward in
6
the riser.
us
cr
ip t
3
(for 10-lump model of Jacob et al., 1976)
eq.(100)
an
7
Where, ρc/εg is the mass density of the catalyst relative to the gas volume; kh is the heavy aromatic ring
9
adsorption coefficient; and CAh is the concentration of aromatic ring in heavy fuel oil. In this rate equation,
10
influence of the adsorption of heavy aromatic ring on the catalyst active sites and consequently on the rate of
11
the reaction is considered.
12
In all three types of rate equations, the catalyst deactivation function ( ) can be given either as a function of
13
residence time of the catalyst (tc) or concentration of the coke on the catalyst. In addition, several different
14
types of deactivation functions are available in the literature (see Table-8).
Ac ce
pt
ed
M
8
Table-8: Models for catalyst deactivation Weekman Jr.(1969) α and n are the rate constants.
Jacob et al. (1976)
= partial pressure of oil at inlet; α, β, γ = constants Krambeck (1991) b = empirical constant Ccoke= coke concentration (kg coke/kg cat) Farag et al. (1993) Kc = deactivation constant (kg catalyst/kmol) Ccoke= coke concentration (kmol/m3)
Page 77 of 77
Page 77 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Pitault et al. (1995) Ccoke concentration in wt% A and B are empirical constants. As mentioned-above, rate constants used in the rate equations are estimated by lab-scale or pilot-scale
2
experimental data. Change in mass fraction of various species results in change in volume of the gas.
3
Furthermore, variation in the temperature and pressure conditions also brings in change in volume of the gas.
4
Volumetric expansion in gas phase is generally calculated by using the ideal gas law to calculate gas density
5
as:
cr
ip t
1
eq.(101)
an
us
6
Gas density calculated from mass fraction of various species is used in the mass and momentum balance
8
equations for the gas phase.
9
Coke deposition on catalyst surface causes change in catalyst activity and density. Change in the catalyst
10
activity due to the presence of coke is incorporated by using deactivation model as described in the kinetic
11
model section. In most CFD models (Table-9), the coke is assumed as a separate lump in the gas phase
12
instead of a separate solids phase.
13
The cracking reactions are endothermic, and their effect on the temperature is accounted for by a source term
14
in the enthalpy balance equation (Sgh in eq.43) for the gas phase. Heat of reaction is used to calculate the
15
source term. Chang et al. (2012) used the value of this source term equal to 9.127×103 KJ/kg of coke
16
produced; while Theologos and Markatos (1993) and Theologos et al. (1997) have used a constant value of
17
300 KJ/kg of hydrocarbon converted and 465 KJ/kg of feed respectively. Nayak et al. (2005) and Behjat et
18
al. (2011) have used reference enthalpy or heat of formation of each lump of the four-lump model to
19
calculate the heat of reaction in each control volume. In Nayak et al. (2005), the heat of formation values of
20
each lump are back calculated from heat of reactions estimated by Han and Chung (2000).
21
5.4. Predictions
22
The CFD models have been used to evaluate various designs and operational alternatives. Some of important
23
predictions can be summarised as follows.
Ac ce
pt
ed
M
7
Page 78 of 78
Page 78 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
(1) Effect of CTO: Nayak et al. (2005), Wu et al. (2010), Zhu et al. (2011) and Behjat et al. (2011b) have
2
investigated the effect of CTO on the performance of FCC riser. Nayak et al. (2005) have found that an
3
increase in CTO from 6 to 12 increases the conversion 70 to 79%. Zhu et al. (2011) have also found that
4
the increase in CTO from 5 to 9 increases the yield of gasoline from 45% to 60%. Zhu et al. (2011) also
5
report negligible the effect of CTO on the yield of coke. Wu et al. (2010) have investigated the effect of
6
CTO (in a range of 5-20) on the clustering phenomenon and residence time of catalyst in the FCC riser
7
by using Langrangian approach. They note that increase in CTO results in higher clustering, suggesting
8
higher back mixing in riser. Wu et al. (2010) have further observed that an increase in CTO from 5 to 15
9
almost doubles the residence time of the catalyst phase in riser.
cr
ip t
1
us
10
(2) Effect of feed nozzles: Theologos et al. (1997), Theologos et al. (1999) and Li et al. (2013) have
12
investigated the effect of feed nozzle configurations and droplet size on the performance of FCC riser.
13
Theologos et al. (1997) have investigated 3 and 12-nozzle configurations, and found that the 12-nozzle
14
configurations increases gasoline fraction at outlet from 45 wt % to 48 wt %. Theologos et al. (1999)
15
have investigated the effect of three different feed droplet sizes (30, 100 and 500 μm) on vaporization,
16
cracking reaction initiation, selectivity and overall performance. They have found that a higher degree of
17
atomization (lower droplet size) gives faster feed vaporization and initialization of the cracking
18
reactions. Nayak et al. (2005) have also found that an increase in the droplet size from less than 400 μm
19
to 2mm results in decrease in the conversion from 85% to less than 70%. As the initial droplet size is
20
reduced, more uniformity in temperature profiles is achieved with higher conversion rates and gasoline
21
selectivity. Li et al. (2013) have analysed the effect of jet velocity (41.7, 62.5, and 83.3 m/s), positions
22
(the distance between the nozzle opening and walls - 3.5, 4.5, and 5.5 m) and angle of feed injection
23
nozzles (15°, 30°, and 45°) on flow field and cracking reactions in the feedstock mixing zone of riser.
24
Their simulations suggest that a higher jet velocity significantly increases inhomogeneity in radial
25
distributions of velocity and volume fraction of the solids phase. A higher angle of nozzles gives more
26
homogeneous radial distribution of the gas phase temperature. However, the angle of nozzles does not
27
impact axial distribution of the solids volume fraction. Increasing the distance between nozzle opening
28
and walls also has marginal impact on flow field. Decreasing the injection angle from 60˚ to 45˚
29
increases the conversion of VGO from 94.8% to 95.7% and the yield of coke from 8.25% to 8.56%
30
(Chang et al., 2012).
Ac ce
pt
ed
M
an
11
31
Page 79 of 79
Page 79 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(3) Effect of temperature: Zhu et al. (2011) have studied the effects of inlet temperatures of the catalyst
2
phase (600, 850 and 1000K) on the yields and conversions. Low inlet temperature of 600K for the
3
catalyst phase gives the conversion of approximately 20% and yield of the gasoline of approximately
4
15%. Increase in the inlet temperature to 1000 K then resulted in the fractional conversion of
5
approximately 80% and the yield of the gasoline of approximately 60%.
ip t
6
(4) Effect of catalyst: Behjat et al. (2011) have studied the effect of catalyst deactivation on mass fraction of
8
various species. They compare predicted variation in the mass fraction along height using a CFD model
9
with and without deactivation of catalyst. The mass fraction of gas oil at outlet is almost 10% higher
10
when catalyst deactivation is considered, suggesting lower conversion of gas oil. Corresponding
11
reduction (approximately 10%) in mass fraction of gasoline at the outlet is also observed. Catalyst
12
deactivation affects light gases lump the most as its mass fraction decreases by approximately 50%.
an
us
cr
7
13
(5) Effect of riser geometry: Gao et al. (1999, 2001) have found that an optimum yield of desirable FCC
15
products is achieved at a riser height that is less than 10 m above feed inlet. Beyond this height,
16
excessive cracking reactions result in increased yield of by-products such as gas and coke at the expense
17
of desirable products. They have also suggested that injection of water as a reaction-terminating medium
18
above optimum riser height can be an effective option for optimizing the yield of desirable products. Lan
19
et al. (2009) have used CFD model to investigate performance of a two stage riser technology, where
20
two risers with different diameter and length replace conventional single riser. Gan et al. (2011) have
21
further investigated the impact of change in the diameter and nozzles in the two stage FCC risers. In this
22
study, the feed injection zone was fitted with a section with a diameter twice that of riser. Gan et al.
23
(2011) concluded that the changing diameter riser provides more mixing and is well suitable for the fast
24
reactions occurring at the feed injection section. Lopes et al. (2011) have analysed flow patterns near the
25
inlet section of the riser, and found a large-scale inhomogeneity in the distribution of solids volume
26
fraction near the catalyst inlet. Chang et al. (2012) have proposed use of an airlift loop mixer at the inlet
27
section of riser to improve hydrodynamics and mixing. Trujillo and Wilde (2012) have compared
28
performance of riser with a rotating fluidized bed in static geometry configuration, where solids are
29
injected in radially inward direction and VGO is injected in radially tangential direction. Due to intense
30
mixing in rotating fluidized bed configuration, significant process intensification is achieved in terms of
31
the bed height and corresponding rector volume required for a given conversion.
Ac ce
pt
ed
M
14
Page 80 of 80
Page 80 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
5.5. Shortcomings
2
(1) CFD models of FCC riser have not been thoroughly validated against experimental data. Bottleneck in
3
this aspect comes from unavailability of a comprehensive set of experimental data. Experiments which
4
give measurements on both hydrodynamics and conversions at multiple locations in a given riser systems
5
are not available in literature. Consequently, several models have been validated against the
6
experimental data of cold-flow riser and few models have been validated against the available plant data.
7
Without thorough validation, the CFD models have only been used for qualitative analysis of the effect
8
design alternatives or operating conditions on the FCC performance.
cr
ip t
1
us
9
(2) There is a lack of coherence between evaluation of closure models for gas-solid flow in riser and coupled
11
CFD models for reactive flow in FCC riser. Recommendations on the selection of closure models and
12
modifications of these models provided by the gas-solid flow modelling efforts have not been
13
implemented in the coupled CFD models of FCC riser. For example, it has been repeatedly shown that
14
the multiscale drag models such as the EMMS and SGS models significantly improve flow predictions.
15
Despite of this, most of the coupled CFD models (Table-6) have used the Gidaspow drag models.
M
an
10
ed
16
(3) Most of the coupled CFD models (Table-6) have ignore the droplet phase by assuming instantaneous
18
vaporization. While this assumption makes the model and computation less complex, ignoring the
19
droplet phase leads to inaccurate predictions in more critical bottom section of the riser.
Ac ce
20
pt
17
21
(4) The coupled CFD models of Theologos et al. (1997), Gao et al. (2000), Behjat et al. (2011) have
22
considered homogeneous fast vaporization of feed droplets, while Nayak et al. (2005) have implemented
23
a mechanistic heterogeneous model for droplet vaporization. The comparison (see Figure-13) of the
24
vaporization models for single droplet vaporization under FCC conditions shows wide differences in the
25
predicted vaporization times. The selection of the vaporization model therefore needs to be validated by
26
comparing predicted droplet vaporization height with real FCC riser. Without such validation, the
27
predictions of the impact of the droplet diameter or design configuration of feed nozzles on FCC
28
performance remain unconvincing.
29
Page 81 of 81
Page 81 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(5) The gas-solid heat transfer coefficient is mostly calculated using Ranz and Marshall (1952). Behjat et al.
2
(2011) have used a correlation of Gunn (1978), whereas Theologos and Markatos (1993) and Theologos
3
et al. (1997) have used a correlation of Kothari (1967). It is necessary to evaluate the effect of various
4
heat transfer coefficient models on predictions. Currently, systematic study on the effect of the heat
5
transfer models on FCC performance is not available.
ip t
6
(6) All the coupled CFD models (Table-6) have used lump kinetic models for cracking reactions. For each
8
lump kinetic model, several sets of estimated kinetic parameters are available in the literature. These
9
parameters have been derived by different catalyst, feedstock, operating conditions and flow conditions.
10
The comparison of the 3, 4, and 5-lump kinetic parameters (Figure-14) shows that different sets of
11
kinetic parameters do not predict experimental values of the gas oil conversion reported by Nace et al.
12
(1971). This clearly suggests that the available kinetic parameters are just tuned for a specific set(s) of
13
experimental data, and they do not represent intrinsic kinetic but they are highly influenced by other heat
14
and mass transfer processes.
M
an
us
cr
7
15
(7) Selection of appropriate set of kinetic parameters is not adequately explained in the previous CFD
17
models. Furthermore, none of the previous studies have investigated the sensitivity of predictions on
18
kinetic parameters. The selection of the kinetic parameters has been arbitrary in the previous efforts. For
19
example, Nayak et al. (2005) have used kinetic constants of Pitault et al. (1994) and activation energy of
20
Lee et al. (1989) for the four-lump kinetic model, though experimental data used to estimate the kinetic
21
constants and activation energy by Pitault et al. (1994) and Lee et al. (1989) have been quite different.
22
Nayak et al.’s modelling approach to couple cracking kinetics has been followed by Behjat et al. (2011)
23
and Lopez et al. (2011).
pt
Ac ce
24
ed
16
25
(8) The volume expansion due to vaporization of feedstock and cracking reactions can significantly change
26
the hydrodynamics in the riser. Many studies (Theologos and Markatos, 1993; Theologos et al., 1997;
27
Theologos et al., 1999; Chang and Zhou 2003; and Das et al., 2003) have ignored the volume expansion
28
in their model.
Page 82 of 82
Page 82 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
6. Recommendation for future work
2
This review identifies several shortcomings related to the selection of the KTGF closure models, modelling
3
of interphase exchange and turbulent forces, model dependency on several empirical constants, uncertainty
4
over the selection of the droplet vaporization model under riser flow conditions, uncertainty over the
5
selection of kinetic parameters for a given lump kinetic model and lack of experimental data for validation of
6
CFD models. All these shortcomings must be addressed in future research for development of accurate and
7
reliable computational model for FCC riser. Following suggestions are useful for future considerations:
ip t
1
(1)
KTGF: Riser simulations have been predominantly conducted by using the E-E model with the
us
9
cr
8
KTGF approach, which requires several closure models to represent rheology of the solids phase.
11
Several alternate models have been proposed for every closure term of solids phase properties.
12
Consequently, different models have used a different combination of closure laws. There are very
13
few comparative studies which have investigated the influence of available alternate models on flow
14
predictions of riser flow. In future, a detailed study that compares influence of all KTGF closure
15
models on flow predictions at various riser flow conditions should be carried out. Such study can
16
provide a guide line for selection of a set of closure models that are suitable for riser simulations.
18
(2)
Use of discrete particle model: On the basis of published CFD models, the KTGF approach appears
pt
17
ed
M
an
10
to be reasonable basis for estimating closure models for solids phase. Future work should include
20
discrete particle method (DPM) or discrete element method (DEM) to simulate the flow of the solids
21
phase. In DEM, the particle-particle and particle-walls collisions are captured by calculating contact
22
forces on each particle using the soft-sphere or hard-sphere approach. Both these approach use
23
spring stiffness and damping coefficient, which are derived from the material properties of the
24
solids. The DEM approach will ensure a realistic resolution of the solids phase without use of the
25
KTGF closure models. This approach is computationally more expensive; but due to recent advances
26
in high performance computing facilities and advances in DEM tools with ability to track millions of
27
particles in the flow domain, it is now possible to conduct DEM simulations of riser. In this
28
direction, DEM models for riser flow have been reported by Zhang et al. (2008), Wu et al. (2010)
29
and Chu and Yu (2012). However, these studies either use scaled-down geometry or consider
30
significantly lower number of particles then those found in actual FCC riser system.
Ac ce
19
31 Page 83 of 83
Page 83 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
(3)
Use of multiscale drag: Gas-solid drag model is the dominant closure model that significantly affects flow predictions. Previous studies have shown that use of multiscale drag models (such as EMMS
3
and SGS-filtered models) has considerably improved flow predictions (radial and axial profiles of
4
solids volume fraction) than those from the conventional drag models (Benayahia, 2009; Shah et al.,
5
2014). Despite their usefulness, the multiscale drag models have not been used by published coupled
6
CFD models for FCC riser. Therefore, future simulations of FCC riser should use multiscale drag
7
models instead of the conventional drag models.
ip t
2
9
(4)
cr
8
Resolution of clusters: Aggregation of particles, formation/break-up of clusters, and turbulentparticle interactions are key phenomena in upward gas-solid flow inside the riser. These phenomena
11
occur at meso-scale. Developing closure laws applicable to the riser that sufficiently account for
12
these meso-scale phenomena is a major challenge.
an
us
10
13
Recently, Shaffer et al. (2013) used high-speed PIV images to capture clusters away from the walls.
15
Shaffer et al. (2013) and Coco et al. (2013) have provided direct and indirect measurements for
16
cluster properties. As mentioned in Section-4, Nicolai et al. (2014) have performed experiments to
17
investigate particle aggregation under turbulent flow conditions. More measurements like these are
18
required to characterise clusters in the riser flow. Another useful approach is particle resolved DNS,
19
where fluid-particle and inter particle interactions can be investigated by applying first principles.
20
Current, numerical observations from DNS are limited to several thousand particles positioned at
21
fixed locations in a fictitious cubical flow domain with periodic boundary conditions. DNS with
22
large number of particles, bigger size of flow domain, and free movement of particles that
23
experience both drag and collisional forces should be performed. The suggested simulations would
24
provide numerical observations to characterize clusters with respect to flow conditions. Numerical
25
results from DNS can then become input to multiscale drag, as well as particle turbulent models used
26
in the FCC simulations..
Ac ce
pt
ed
M
14
27
28
29
(5)
LES for turbulent parameters: In the k- turbulent models, extreme assumption of isotropy is made,
where the turbulent viscosity is assumed to be same in all directions. In addition, the turbulence Page 84 of 84
Page 84 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
parameters are estimated on the basis of some ideal flow. In gas-solid flow, the effect of the solids
2
on the gas phase turbulent properties is significant. Capturing such effect in turbulent closure models
3
used in the two fluid model is critical, and hence it is proposed to undertake LES simulations. In
4
LES of gas-solid flow (Yomamoto et al., 2001; and Vreman et al., 2009), the gas phase is modelled
5
as a continuous phase whereas the solids are modelled by using Lagrangian approach. The gas phase
6
turbulent is modelled by using the LES model with the Smagorinski SGS model (section-3.1.2.2.) or
7
improved Smagorinski model (Boivin et al., 2000). Here, the average balance equation of the gas
8
phase and the force balance equations for the solids are coupled with the drag force term. This type
9
of simulations can be used to estimate the effect of the solids phase on turbulent properties of the gas
10
phase. Such estimates of turbulent properties can be made for different flow conditions at number
11
locations in riser (say at 10 points), which should be selected on the basis of various ranges of power
12
consumption per unit mass as well as solids volume fraction. Interpolated values of these turbulence
13
parameters can be used in the two fluid models for the riser simulation.
an
us
cr
ip t
1
14 (6)
Drift flux analysis: The cold-flow experiments and simulation studies have reported several sets of
M
15
radial and axial profiles of volume fraction and velocity of the solids phase. These concentration
17
profiles will enable the estimation of drift flux constants c0 and c1 (section-3.5.). The coefficient c0
18
represents the transverse volume fraction profile and the value of c1 represents the slip velocity.
19
These values will give additional basis for confirmation of drag coefficient.
21
(7)
For the detailed understanding of transport phenomenon in FCC risers, it is recommended to identify
Ac ce
20
pt
ed
16
22
the turbulent structures (together with the size and shape distribution) and their dynamics (together
23
with velocity and energy distributions). For this purpose, the procedures are given by Kulkarni et al.
24
(2001) and Joshi et al (2009). Mathpati and Joshi (2007) and Mathpati et al (2009) have described
25
the relationships between the structure dynamics and the transport phenomena.
26 27
(8)
Entry section: It is seen from experimental and simulation results that the concentration and
28
temperature profiles drop exponentially in the initial few meter height of the FCC riser. This
29
suggests that most of the variation in concentration and temperature occurs in the entry region,
30
which probably further indicates that majority of cracking reactions also occur in the entry section.
31
Therefore, the performance of riser strongly depends on the performance of entry section. Hence, a Page 85 of 85
Page 85 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
1
comprehensive research on the entry section is necessary. The future work should include (a)
2
transient 3D measurement of concentration and velocity profiles of gas and solids phase, (b)
3
dynamic pressure measurement, (c) time series analysis of all variables in order to extract flow
4
parameters, and (d) CFD simulation and validation with experimental measurements of the entry
5
section.
7
(9)
ip t
6
Pilot scale experiments: This review found a serious deficit of reliable plant experiment data in open literature. Future work should include an extensive experimental program for estimation of three
9
dimensional profiles in riser flow domain. For reactive flow, simultaneous measurements of
10
concentration, temperature, pressure, and catalyst volume fraction are required; while for
11
hydrodynamics, axial and radial profiles of both volume fraction and velocity in a single riser system
12
are desired.
an
us
cr
8
In conclusion, CFD modelling of FCC riser has come a long way, but it has still remained an art due to use
14
of various combinations of closure models and boundary conditions without rigorous experimental
15
validation. We feel that more fundamental work that consists of both numerical modelling and
16
experimentation is necessary to create greater understanding of underlying phenomena. Both simulations and
17
experiments can be conducted at various scales. At particle and cluster scale, DNS can provide more
18
observations on the interphase interactions and droplet vaporization. At equipment scale, pilot-scale
19
experimentations can provide much needed data for validation. Derivations of fundamental correlations for
20
interphase forces, turbulent interactions and droplet vaporization along with an extensive validation using
21
experiments are critical for development of more reliable computational models of FCC riser.
22
7. Nomenclature
23 24 25 26 27 28 29 30 31 32 33 34 35 36
Asi Ad Ariser Ainlet CAh CD0 Ci,s Ci,cont cμ, c1ε, c2ε cj cj0 c0 and c1 Cl Cp
Ac ce
pt
ed
M
13
Interfacial area between gas and solids, m2 Droplet surface area, m2 Cross sectional area of riser, m2 Cross sectional area of inlet of riser, m2 Concentration of heavy aromatic ring in heavy fuel oil, kg/m3 Drag coefficient of a particle Concentration of vaporizing component at the droplet surface, kmol/m3 Concentration of vaporizing component in the bulk of the continuous phase, kmol/m3 Constant in k-ϵ model Concentration of jth lump in gas phase, kg/m3 Initial concentration of jth lump in gas phase, kg/m3 Drift flux constants Lift coefficient Heat capacity, J/kg K Page 86 of 86
Page 86 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
ed
M
an
us
cr
ip t
Smagorinsky constant Diameter of riser, m Diameter of particle, m Diameter of droplet, m Diameter of solids, m Fluctuating energy due to random motion of particle, kg/m3 s Activation energy of ith reaction, J/mol Coefficient of restitution of particle collisions Coefficient of restitution of particle-wall collisions Contact force on particle, kg/m s2 Drag force, kg/m s2 Drag on a particle, kg/m s2 Gravitational force, kg/m s2 Force on particle, kg/m s2 Lift force, kg/m s2 Empirical constant in frictional pressure equation Vitual mass force, , kg/m s2 Damping coefficient Generation of turbulent kinetic energy due to shear in the gas phase, kg m2/s2 Solids flux, kg/m2s Gravitational acceleration, m/s2 Radial distribution function Height of riser, m Specific enthalpy, J/kg Heat transfer coefficient, J/kg m2 K Heat of vaporization, kJ/kg Thermal conductivity, J/m K Mass transfer coefficient, m/s Gas-solid phase velocity variant Rate constant of consumption of jth lump in ith reaction, m3reactant/m3cat s Heavy aromatic ring absorption coefficient Turbulent kinetic energy, m2/s2 Conductivity for granular temperature, kg/m s Moment of inertia of particle Unit tensor Second invariant of the deviatoric stress tensor Total number of parcels in kth control volume Nusselt number Effective Nusselt number (Nayak et al., 2005) Order of reaction Empirical constant in frictional pressure equation Number of particles in a parcel Molecular weight, kg/kmol Mass, kg Pressure, kg/m s2 Solids pressure due to friction, kg/m s2 Solids pressure, kg/m s2 Empirical constant in frictional pressure equation Model parameter in solids pressure equation Prandtl number Heat flux, kJ/m2 s Heat exchanges, kJ/kg s Universal gas constant, J/kmol K
pt
Cs Dt dp dd ds EΘ Ej ess ew Fc FD Fd,p Fg Fp Flift Fr Fv,m fs Gk,g Gs g go,ss H h h hfg K K Kgc Kij Kh k kΘs Ip I̿ I2D NT Nu Nueff n n np MW m P Pfri Ps p ps* Pr q ΔQ R
Ac ce
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
Page 87 of 87
Page 87 of 114
cr
Model parameter in solids pressure equation Volume fraction Maximum packing limit Minimum solid volume fraction from where friction between particle is dominant Time-averaged holdup of solids Density, kg/m3 Stress tensor, kg/m s2 Drag coefficient, kg/m2 s Filter size Kronecker delta Viscosity, kg/m s Turbulent viscosity, kg/m s Bulk viscosity of solids phase, kg/m s Component of bulk viscosity of solids phase due to collision between particles, kg/m s Component of bulk viscosity of solids phase due to friction, kg/m s Component of bulk viscosity of solids phase due to free motion of particles, kg/m s Turbulent dissipation rate, m2/s3 Interphase turbulent kinetic energy source term, m2/s2 Interphase turbulent energy dissipation source term, m2/s2 Reynolds stress due to turbulent interactions between phases Turbulent Prandtlnumber Granular temperature, m2/s2 Dissipation of granular energy due to inelastic collision between particles, kg/m3 s
pt
Greek letters α ε εs,max εs,min εs ρ τ β Δs δij μ μg,t μs μs,col μs,fri μs,kin ϵ Πk Πε Π(R,ij) σk, σε Θ γΘs
ed
M
an
us
Rij Reynold stresses, kg/m s2 Rep Particle Reynolds number Particle Reynolds number Res Droplet Reynolds number Red Rate of consumption of jth lump in ith reaction, kg/m3s rij S Source term in mass balances, kg/m3 Sh Source term in enthalpy balances, kg/m3 Strain rate Sij T Temperature, K Torque on particle, kg m3/s2 Tp t Time, s u Velocity, m/s Velocity of the gas phase at location Xp ug(Xp) up Velocity of parcel u’ Fluctuating velocity, m/s u Velocity, m/s Relative velocity, m/s urel Drift velocity, m/s udrift SGS stress, , kg/m s2 (u_i^' u_j^' ) ̅^SGS Us Slip velocity, m/s Volume of the droplet, m3 vd Volume of kth control volume, m3 vk Volume of parcel, m3 vp V Superficial velocity, m/s Slip velocity, m/s Vslip Location of parcel Xp y Mass fraction Distance from wall in wall unit y+
Ac ce
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
ip t
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Page 88 of 88
Page 88 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
ϕgs λs ϕ ϕ θ τs,w γ(Θs,w) Ω
Exchange of kinetic energy of random fluctuations in particle velocity to gas phase, kg m2/s2 Shear viscosity of solids phase, kg/m s Specularity coefficient Catalyst deactivation factor Angle of internal friction Wall shear, kg/m s2 Dissipation of granular energy due to inelastic collision between particles and wall, kg/m3 s Angular velocity of particle, rad/s
Subscripts g s SGS m p w d gs sg gd dg sd ds cont i
gas phase solids phase Sub grid scale mixture properties Parcel Wall Droplet From gas phase to solids phase From solids phase to gas phase From gas phase to droplet phase From droplet phase to gas phase From solids phase to droplet phase From droplet phase to solids phase Continuous phase Vaporizing phase
26
8. References
27
Abul-Hamayel, M.A., 2003. Kinetic modeling of high-severity fluidized catalytic cracking☆. Fuel 82, 1113-
28
1118.
29
Abramzon, B., Sirignano, W., 1989. Droplet vaporization model for spray combustion calculations.
30
International Journal of Heat and Mass Transfer 32, 1605-1618.
31
Abramzon, B., Sazhin, S., 2006. Convective vaporization of a fuel droplet with thermal radiation absorption.
32
Fuel 85, 32-46.
33
Agrawal, K., Loezos, P.N., Syamlal, M., Sundaresan, S., 2001. The role of meso-scale structures in rapid
34
gas–solid flows. Journal of Fluid Mechanics 445, 151-185.
35
Ahuja, G.N., Patwardhan, A.W., 2008. CFD and experimental studies of solids hold-up distribution and
36
circulation patterns in gas–solid fluidized beds. Chemical Engineering Journal 143, 147-160.
37
Ali, H., Rohani, S., 1997. Dynamic modeling and simulation of a riser-type fluid catalytic cracking unit.
38
Chemical Engineering & Technology 20, 118-130.
Ac ce
pt
ed
M
an
us
cr
ip t
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Page 89 of 89
Page 89 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Almuttahar, A., Taghipour, F., 2008. Computational fluid dynamics of high density circulating fluidized bed
2
riser: Study of modeling parameters. Powder Technology 185, 11-23.
3
Alvarez-Castro, H., Matos, E., Mori, M., Martignoni, W., Ocone, R., 2015. Evaluation of the Performance of
4
the Riser in the Fluid Catalytic Cracking Process. Petroleum Science and Technology 33, 579-587.
5
Alvarez-Castro, H., Matos, E., Mori, M., Martignoni, W., Ocone, R., 2015. Analysis of Process Variables via
6
CFD to Evaluate the Performance of a FCC Riser. International Journal of Chemical Engineering 2015.
7
Allen, M.P., Tildesley, D.J., 1989. Computer simulation of liquids. Oxford university press.
8
Ambler, P., Milne, B., Berruti, F., Scott, D., 1990. Residence time distribution of solids in a circulating
9
fluidized bed: experimental and modelling studies. Chemical Engineering Science 45, 2179-2186.
us
cr
ip t
1
Ancheyta-Juárez, J., López-Isunza, F., Aguilar-Rodríguez, E., Moreno-Mayorga, J.C., 1997. A strategy for
11
kinetic parameter estimation in the fluid catalytic cracking process. Industrial & Engineering Chemistry
12
Research 36, 5170-5174.
13
Anderson, T.B., Jackson, R., 1967. Fluid mechanical description of fluidized beds. Equations of motion.
14
Industrial & Engineering Chemistry Fundamentals 6, 527-539.
15
Andrews, M., O'rourke, P., 1996. The multiphase particle-in-cell (MP-PIC) method for dense particulate
16
flows. International Journal of Multiphase Flow 22, 379-402.
17
Andrews IV, A.T., Loezos, P.N., Sundaresan, S., 2005. Coarse-grid simulation of gas-particle flows in
18
vertical risers. Industrial & Engineering Chemistry Research 44, 6022-6037.
19
Arastoopour, H., Pakdel, P., Adewumi, M., 1990. Hydrodynamic analysis of dilute gas--solids flow in a
20
vertical pipe. Powder Technology 62, 163-170.
21
Auzerais, F., Jackson, R., Russel, W., 1988. The resolution of shocks and the effects of compressible
22
sediments in transient settling. Journal of Fluid Mechanics 195, 437-462.
23
Bader, R., Findlay, J., Knowlton, T., 1988. Gas/solids flow patterns in a 30.5-cm-diameter circulating
24
fluidized bed. Circulating fluidized bed technology II, 123-137.
25
Bai, D. R., Jin, Y., Yu, Z. Q., Zhu, J. X., 1992. The axial distribution of the cross-sectionally averaged
26
voidage in fast fluidized beds. Powder Technology 71, 51-58.
Ac ce
pt
ed
M
an
10
Page 90 of 90
Page 90 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Bai, D., Issangya, A.S., Grace, J.R., 1999. Characteristics of gas-fluidized beds in different flow regimes.
2
Industrial & Engineering Chemistry Research 38, 803-811.
3
Baharanchi, A.A., Gokaltun, S., Dulikravich, G., 2015. Performance improvement of existing drag models in
4
two-fluid modeling of gas–solid flows using a PR-DNS based drag model. Powder Technology 286, 257-
5
268.
6
Balzer, G., Simonin, O., Boelle, A., Lavieville, J., 1996. Unifying modelling approach for the numerical
7
prediction of dilute and dense gas-solid two phase flow. 5th International Conference on Circulating
8
Fluidized Beds, Beijing, China, 1996
9
Batchelor, G., 1988. A new theory of the instability of a uniform fluidized bed. Journal of Fluid Mechanics
us
cr
ip t
1
193, 75-110.
11
Beetstra, R., van der Hoef, M.A., Kuipers, J.A.M., 2006. A lattice-Boltzmann simulation study of the drag
12
coefficient of clusters of spheres. Computers & Fluids 35, 966-970.
13
Beetstra, R., van der Hoef, M.A., Kuipers, J.A.M., 2007. Drag force of intermediate Reynolds number flow
14
past mono-and bidisperse arrays of spheres. AIChE Journal 53, 489-501.
15
Behjat, Y., Shahhosseini, S., Marvast, M.A., 2011. CFD analysis of hydrodynamic, heat transfer and reaction
16
of three phase riser reactor. Chemical Engineering Research and Design 89, 978-989.
17
Benyahia, S., 2009. On the Effect of Subgrid Drag Closures. Industrial & Engineering Chemistry Research
18
49, 5122-5131.
19
Benyahia, S., 2012. Analysis of model parameters affecting the pressure profile in a circulating fluidized
20
bed. AIChE Journal 58, 427-439.
21
Benyahia, S., Arastoopour, H., Knowlton, T.M., Massah, H., 2001. Simulation of particles and gas flow
22
behavior in the riser section of a circulating fluidized bed using the kinetic theory approach for the
23
particulate phase. Powder Technology 112, 24-33.
24
Benyahia, S., Gonzalez Ortiz, A., Paredes, P., Ignacio, J., 2003. Numerical analysis of a reacting gas/solid
25
flow in the riser section of an industrial fluid catalytic cracking unit. International Journal of Chemical
26
reactor engineering 1.
Ac ce
pt
ed
M
an
10
Page 91 of 91
Page 91 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Benyahia, S., Syamlal, M., O'Brien, T.J., 2006. Extension of Hill-Koch-Ladd drag correlation over all ranges
2
of Reynolds number and solids volume fraction. Powder Technology 162, 166-174.
3
Benyahia, S., Syamlal, M., O'Brien, T.J., 2007. Study of the ability of multiphase continuum models to
4
predict core‐annulus flow. AIChE Journal 53, 2549-2568.
5
Benyahia, S., Galvin, J.E., 2010. Estimation of numerical errors related to some basic assumptions in discrete
6
particle methods. Industrial & Engineering Chemistry Research 49, 10588-10605.
7
Berruti, F., Pugsley, T., Godfroy, L., Chaouki, J., Patience, G., 1995. Hydrodynamics of circulating fluidized
8
bed risers: a review. The Canadian Journal of Chemical Engineering 73, 579-602.
9
Bhusarapu, S., Al-Dahhan, M.H., Dudukovic, M.P., Trujillo, S., O'Hern, T.J., 2005. Experimental Study of
10
the Solids Velocity Field in Gas- Solid Risers. Industrial & Engineering Chemistry Research 44, 9739-9749.
11
Bolio, E.J., Yasuna, J.A., Sinclair, J.L., 1995. Dilute turbulent gas‐solid flow in risers with particle‐particle
12
interactions. AIChE Journal 41, 1375-1388.
13
Bollas, G., Papadokonstadakis, S., Michalopoulos, J., Arampatzis, G., Lappas, A., Vasalos, I., Lygeros, A.,
14
2003. Using hybrid neural networks in scaling up an FCC model from a pilot plant to an industrial unit.
15
Chemical Engineering and Processing: Process Intensification 42, 697-713.
16
Bollas, G., Vasalos, I., Lappas, A., Iatridis, D., 2002. Modeling small-diameter FCC riser reactors. A
17
hydrodynamic and kinetic approach. Industrial & Engineering Chemistry Research 41, 5410-5419.
18
Bollas, G., Lappas, A., Iatridis, D., Vasalos, I., 2007. Five-lump kinetic model with selective catalyst
19
deactivation for the prediction of the product selectivity in the fluid catalytic cracking process. Catalysis
20
Today 127, 31-43.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 92 of 92
Page 92 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Brereton, C., Grace, J., 1993. Microstructural aspects of the behaviour of circulating fluidized beds.
2
Chemical Engineering Science 48, 2565-2572.
3
Buchanan, J.S., 1994. Analysis of heating and vaporization of feed droplets in fluidized catalytic cracking
4
risers. Industrial & Engineering Chemistry Research 33, 3104-3111.
5
Campbell, D., Wojciechowski, B., 1969. Theoretical patterns of selectivity in aging catalysts with special
6
reference to the catalytic cracking of petroleum. The Canadian Journal of Chemical Engineering 47, 413-
7
417.
8
Cao, J., Ahmadi, G., 1995. Gas-particle two-phase turbulent flow in a vertical duct. International Journal of
9
Multiphase Flow 21, 1203-1228.
us
cr
ip t
1
Chalermsinsuwan, B., Prajongkan, Y., Piumsomboon, P., 2013. Three-dimensional CFD simulation of the
11
system inlet and outlet boundary condition effects inside a high solid particle flux circulating fluidized bed
12
riser. Powder Technology 245, 80-93.
13
Chan, C.W., Seville, J.P., Parker, D.J., Baeyens, J., 2010. Particle velocities and their residence time
14
distribution in the riser of a CFB. Powder Technology 203, 187-197.
15
Chang, J., Meng, F., Wang, L., Zhang, K., Chen, H., Yang, Y., 2012. CFD investigation of hydrodynamics,
16
heat transfer and cracking reaction in a heavy oil riser with bottom airlift loop mixer. Chemical Engineering
17
Science 78, 128-143.
18
Chang, J., Cai, W., Zhang, K., Meng, F., Wang, L., Yang, Y., 2014. Computational investigation of the
19
hydrodynamics, heat transfer and kinetic reaction in an FCC gasoline riser. Chemical Engineering Science
20
111, 170-179.
21
Chang, S.L., Zhou, C., 2003. Simulation of FCC riser flow with multiphase heat transfer and cracking
22
reactions. Computational Mechanics 31, 519-532.
23
Chapman, S., Cowling, T.G., 1970. The mathematical theory of non-uniform gases 3rd ed., Cambridge
24
University Press.
25
Chen, C., Werther, J., Heinrich, S., Qi, H.-Y., Hartge, E.-U., 2013. CPFD simulation of circulating fluidized
26
bed risers. Powder Technology 235, 238-247.Chu, K.W., Yu, A.B., 2012. A Novel Circulating Fluidized Bed
27
to Improve Fluid-Solids Cotacting, Ninth International Conference on CFD in the Minerals and Process
28
Industries, CSIRO, Melbourne, Australia.
Ac ce
pt
ed
M
an
10
Page 93 of 93
Page 93 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Cheng, Y., Wei, F., Yang, G., Jin, Y., 1998. Inlet and outlet effects on flow patterns in gas-solid risers.
2
Powder Technology 98, 151-156.
3
Chu, K., Wang, B., Xu, D., Chen, Y., Yu, A., 2011. CFD–DEM simulation of the gas–solid flow in a cyclone
4
separator. Chemical Engineering Science 66, 834-847.
5
Cloete, S., Johansen, S., Braun, M., Popoff, B., Amini, S., 2010. Evaluation of a Lagrangian discrete phase
6
modeling approach for resolving cluster formation in CFB risers, 7th International Conference on Multiphase
7
Flow.
8
Cocco, R., Shaffer, F., Hays, R., Karri, S.R., Knowlton, T., 2010. Particle clusters in and above fluidized
9
beds. Powder Technology 203, 3-11.
us
cr
ip t
1
Cokljat, D., Slack, M., Vasquez, S., Bakker, A., Montante, G., 2006. Reynolds-stress model for Eulerian
11
multiphase. Progress in Computational Fluid Dynamics, An International Journal 6, 168-178.
12
Corella, J., Bilbao, R., Molina, J.A., Artigas, A., 1985. Variation with time of the mechanism, observable
13
order, and activation energy of catalyst deactivation by coke in the FCC process. Industrial & Engineering
14
Chemistry Process Design and Development 24, 625-636.
15
Corella, J., Fernandez, A., Vidal, J.M., 1986. Pilot plant for the fluid catalytic cracking process:
16
Determination of the kinetic parameters of deactivation of the catalyst. Industrial & engineering chemistry
17
product research and development 25, 554-562.
18
Cundall, P.A., Strack, O.D., 1979. A discrete numerical model for granular assemblies. Geotechnique 29, 47-
19
65.
20
Das, A.K., Baudrez, E., Marin, G.B., Heynderickx, G.J., 2003. Three-dimensional simulation of a fluid
21
catalytic cracking riser reactor. Industrial & engineering chemistry research 42, 2602-2617.
22
Das Sharma, S., Pugsley, T., Delatour, R., 2006. Three‐dimensional CFD model of the deaeration rate of
23
FCC particles. AIChE Journal 52, 2391-2400.
Ac ce
pt
ed
M
an
10
Page 94 of 94
Page 94 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Dasgupta, S., Jackson, R., Sundaresan, S., 1994. Turbulent gas‐particle flow in vertical risers. AIChE Journal
2
40, 215-228.
3
Dasgupta, S., Jackson, R., Sundaresan, S., 1998. Gas-particle flow in vertical pipes with high mass loading
4
of particles. Powder Technology 96, 6-23.
5
Deen, N.G., Van Sint Annaland, M., van der Hoef, M.A., Kuipers, J.A.M., 2007. Review of discrete particle
6
modeling of fluidized beds. Chemical Engineering Science 62, 28-44.
7
Derksen, J., 2012. Direct numerical simulations of aggregation of monosized spherical particles in
8
homogeneous isotropic turbulence. AIChE Journal 58, 2589-2600.
9
Derksen, J., Sundaresan, S., 2007. Direct numerical simulations of dense suspensions: wave instabilities in
an
us
cr
ip t
1
liquid-fluidized beds. Journal of Fluid Mechanics 587, 303-336.
11
Derouin, C., Nevicato, D., Forissier, M., Wild, G., Bernard, J.-R., 1997. Hydrodynamics of riser units and
12
their impact on FCC operation. Industrial & Engineering Chemistry Research 36, 4504-4515.
13
Di Felice, R., 1994. The voidage function for fluid-particle interaction systems. International Journal of
14
Multiphase Flow 20, 153-159.
15
Didwania, A., Homsy, G., 1982. Resonant sideband instabilities in wave propagation in fluidized beds.
16
Journal of Fluid Mechanics 122, 433-438.
17
Elghobashi, S., Abou‐Arab, T., 1983. A two‐equation turbulence model for two‐phase flows. Physics of
18
Fluids (1958-1988) 26, 931-938.
19
Enwald, H., Peirano, E., Almstedt, A.E., 1996. Eulerian two-phase flow theory applied to fluidization.
20
International Journal of Multiphase Flow 22, 21-66.
Ac ce
pt
ed
M
10
Page 95 of 95
Page 95 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Ekambara, K., Dhotre, M.T, and Joshi, J.B., 2005. CFD simulations of bubble column reactors: 1D, 2D and
2
3D approach. Chemical Engineering science 60, 6733-6746.
3
Ergun, S., 1952b. Fluid flow through packed columns. Chemical Engineering Progress 48, 89-94.
4
Farag, H., Ng, S., de Lasa, H., 1993. Kinetic modeling of catalytic cracking of gas oils using in situ traps
5
(FCCT) to prevent metal contaminant effects. Industrial & Engineering Chemistry Research 32, 1071-1080.
6
Feng, W., Vynckier, E., Froment, G.F., 1993. Single event kinetics of catalytic cracking. Industrial &
7
Engineering Chemistry Research 32, 2997-3005.
8
Fox, R.O., 2012. Large-eddy-simulation tools for multiphase flows. Annual Review of Fluid Mechanics 44,
9
47-76.
us
cr
ip t
1
Froment, G., 1992. Kinetics and reactor design in the thermal cracking for olefins production. Chemical
11
Engineering Science 47, 2163-2177.
12
Gan, J., Zhao, H., Berrouk, A.S., Yang, C., Shan, H., 2011. Numerical simulation of hydrodynamics and
13
cracking reactions in the feed mixing zone of a multiregime gas–solid riser reactor. Industrial & Engineering
14
Chemistry Research 50, 11511-11520.
15
Gao, J., Xu, C., Lin, S., Yang, G., Guo, Y., 1999. Advanced model for turbulent gas–solid flow and reaction
16
in FCC riser reactors. AIChE Journal 45, 1095-1113.
17
Gao, J., Xu, C., Lin, S., Yang, G., Guo, Y., 2001. Simulations of gas-liquid-solid 3-phase flow and reaction
18
in FCC riser reactors. AIChE Journal 47, 677-692.
19
Gao, J., Xu, C., Lin, S., Yang, G., Guo, Y., 2004. Simulations of gas‐liquid‐solid 3‐phase flow and reaction
20
in FCC riser reactors. AIChE Journal 47, 677-692.
21
Gao, X., Wu, C., Cheng, Y.W., Wang, L.J., Li, X., 2012. Experimental and numerical investigation of solid
22
behavior in a gas–solid turbulent fluidized bed. Powder Technology 228, 1-13.
Ac ce
pt
ed
M
an
10
Page 96 of 96
Page 96 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Garg, R., Tenneti, S., Mohd-Yusof, J., Subramaniam, S., 2010. Direct numerical simulation of gas–solids
2
flow based on the immersed boundary method. Computational Gas–Solids Flows and Reacting Systems:
3
Theory, Methods and Practice. IGI Global, 245-276.
4
Garg, R., Dietiker, J., 2013. Documentation of open-source MFIX–PIC software for gas-solids flows.
5
Garside, J., Al-Dibouni, M.R., 1977. Velocity-voidage relationships for fluidization and sedimentation in
6
solid-liquid systems. Industrial & Engineering Chemistry Process Design and Development 16, 206.
7
Ghatage, S.V., Peng, Z., Sathe, M.J., Doroodchi, E., Padhiyar, N., Moghtaderi, B., Joshi, J.B., Evans, G.M.,
8
2014. Stability analysis in solid–liquid fluidized beds: Experimental and computational. Chemical
9
Engineering Journal 256, 169-186.
us
cr
ip t
1
Gianetto, A., Farag, H.I., Blasetti, A.P., de Lasa, H.I., 1994. Fluid catalytic cracking catalyst for reformulated
11
gasolines. Kinetic modeling. Industrial & Engineering Chemistry Research 33, 3053-3062.
12
Gibilaro, L.G., Di Felice, R., Waldram, P.U., 1985. Generalized friction factor and drag coefficient
13
correlations for fluid-particle interactions. Chemical Engineering Science 40, 1817-1823.
14
Gidaspow, D., 1994. Multiphase flow and fluidization: continuum and kinetic theory descriptions. Academic
15
Press, San Diego, US.
16
Godfroy, L., Patience, G.S., Chaouki, J., 1999. Radial hydrodynamics in risers. Industrial & engineering
17
chemistry research 38, 81-89.
18
Gorham, D., Kharaz, A., 2000. The measurement of particle rebound characteristics. Powder Technology
19
112, 193-202.
20
Grace, J.R., Avidan, A.A., Knowlton, T.M., 1997. Circulating fluidized beds. Blackie academic &
21
professional, UK.
22
Gunn, D., 1978. Transfer of heat or mass to particles in fixed and fluidised beds. International Journal of
23
Heat and Mass Transfer 21, 467-476.
24
Gupta, R.K., Kumar, V., Srivastava, V., 2007. A new generic approach for the modeling of fluid catalytic
25
cracking (FCC) riser reactor. Chemical Engineering Science 62, 4510-4528.
Ac ce
pt
ed
M
an
10
Page 97 of 97
Page 97 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Gupta, R.K., Kumar, V., Srivastava, V., 2010. Modeling of fluid catalytic cracking riser reactor: a review.
2
International Journal of Chemical Reactor Engineering 8.
3
Han, I.S., Chung, C.-B., 2000. Dynamic modeling and simulation of fluidized catalytic cracking process,
4
Part: II: Property estimation and simulation. Chemical Engineering Science 56, 1973–1990.
5
Hari, C., Balaraman, K.S., Balakrishnan, A.R., 1995. Fluid catalytic cracking: Selectivity and product yield
6
patterns. Chemical Engineering & Technology 18, 364-369.
7
Harris, A., Davidson, J., Thorpe, R., 2003. The influence of the riser exit on the particle residence time
8
distribution in a circulating fluidised bed riser. Chemical Engineering Science 58, 3669-3680.
9
Harris, B., Davidson, J., 1994. Modelling options for circulating fluidized beds: a core/annulus deposition
us
cr
ip t
1
model. Circulating Fluidized Bed Technology IV, 32.
11
He, Y., Annaland, S., Deen, N.G., Kuipers, J.A.M., 2006. Gas-solid two-phase turbulent flow in a circulating
12
fluidized bed riser: an experimental and numerical study.
13
He, Y., Deen, N., Annaland, M.v.S., Kuipers, J., 2009. Gas− solid turbulent flow in a circulating fluidized
14
bed riser: numerical study of binary particle systems. Industrial & Engineering Chemistry Research 48,
15
8098-8108.
16
Helland, E., Occelli, R., Tadrist, L., 2000. Numerical study of cluster formation in a gas–particle circulating
17
fluidized bed. Powder Technology 110, 210-221.
18
Hill, R.J., Koch, D.L., Ladd, A.J., 2001a. The first effects of fluid inertia on flows in ordered and random
19
arrays of spheres. Journal of Fluid Mechanics 448, 213-241.
20
Hill, R.J., Koch, D.L., Ladd, A.J., 2001b. Moderate-Reynolds-number flows in ordered and random arrays of
21
spheres. Journal of Fluid Mechanics 448, 243-278.
22
Horio, M., Kuroki, H., 1994. Three-dimensional flow visualization of dilutely dispersed solids in bubbling
23
and circulating fluidized beds. Chemical Engineering Science 49, 2413-2421.
Ac ce
pt
ed
M
an
10
Page 98 of 98
Page 98 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Hrenya, C.M., Sinclair, J.L., 1997. Effects of particle‐phase turbulence in gas‐solid flows. AIChE Journal 43,
2
853-869.
3
Hui, K., Ungar, J., Haff, P., Jackson, R., 1984. Boundary conditions for high-shear grain flows. Journal of
4
Fluid Mechanics 145, 233.
5
Huilin, L., Gidaspow, D., Bouillard, J., Wentie, L., 2003. Hydrodynamic simulation of gas-solid flow in a
6
riser using kinetic theory of granular flow. Chemical Engineering Journal 95, 1-13.
7
Huilin, L., Qiaoqun, S., Yurong, H., Yongli, S., Ding, J., Xiang, L., 2005. Numerical study of particle cluster
8
flow in risers with cluster-based approach. Chemical Engineering Science 60, 6757-6767.
9
Heynderickx, G., De Wilde, J., Marin, G., 2011. Experimental and computational study of T-and L-outlet
an
us
cr
ip t
1
effects in dilute riser flow. Chemical Engineering Science 66, 5024-5044.
11
Ibsen, C.H., Solberg, T., Hjertager, B.H., Johnsson, F., 2002. Laser Doppler anemometry measurements in a
12
circulating fluidized bed of metal particles. Experimental Thermal and Fluid Science 26, 851-859.
13
Ibsen, C.H., Helland, E., Hjertager, B.H., Solberg, T., Tadrist, L., Occelli, R., 2004. Comparison of
14
multifluid and discrete particle modelling in numerical predictions of gas particle flow in circulating
15
fluidised beds. Powder Technology 149, 29-41.
16
Igci, Y., Andrews, I.V., Arthur, T., Sundaresan, S., Pannala, S., O'Brien, T., 2008. Filtered two fluid models
17
for fluidized gas particle suspensions. AIChE Journal 54, 1431-1448.
18
Igci, Y., Sundaresan, S., 2008. Coarse-grid simulation of fluidized gas-particle flows, AICHE annual
19
meeting, Philadelphia, PA.
20
Igci, Y., Sundaresan, S., 2011. Verification of filtered two‐fluid models for gas‐particle flows in risers.
21
AIChE Journal 57, 2691-2707.
Ac ce
pt
ed
M
10
Page 99 of 99
Page 99 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Ishii, H., Nakajima, T., Horio, M., 1989. The clustering annular flow model of circulating fluidized beds.
2
Journal of Chemical Engineering of Japan 22, 484-490.
3
Issangya, A.S., Grace, J.R., Bai, D., Zhu, J., 2000. Further measurements of flow dynamics in a high-density
4
circulating fluidized bed riser. Powder Technology 111, 104-113.
5
Jacob, S.M., Gross, B., Voltz, S.E., Weekman, V.W., 1976. A lumping and reaction scheme for catalytic
6
cracking. AIChE Journal 22, 701-713.
7
Jenkins, J.T., Savage, S.B., 1983. A theory for the rapid flow of identical, smooth, nearly elastic, spherical
8
particles. Journal of Fluid Mechanics 130, 187-202.
9
Jiang, Y., Qiu, G., Wang, H., 2014. Modelling and experimental investigation of the full-loop gas–solid flow
us
cr
ip t
1
in a circulating fluidized bed with six cyclone separators. Chemical Engineering Science 109, 85-97.
11
Jiradilok, V., Gidaspow, D., Damronglerd, S., Koves, W.J., Mostofi, R., 2006. Kinetic theory based CFD
12
simulation of turbulent fluidization of FCC particles in a riser. Chemical Engineering Science 61, 5544-5559.
13
Johnson, P.C., Jackson, R., 1987. Frictional–collisional constitutive relations for granular materials, with
14
application to plane shearing. Journal of Fluid Mechanics 176, 67-93.
15
Joseph, G., Zenit, R., Hunt, M., Rosenwinkel, A., 2001. Particle–wall collisions in a viscous fluid. Journal of
16
Fluid Mechanics 433, 329-346.
17
Joshi, J.B., 2001. Computational flow modelling and design of bubble column reactors. Chemical
18
Engineering Science 56, 5893-5933.
19
Joshi, J.B., Deshpande, N.S., Dinkar, M., Phanikumar, D.V., 2001. Hydrodynamic stability of multiphase
20
reactors. Advances in Chemical Engineering 26, 1-130.
21
Joshi, J.B., Tabib, M. V., Deshpande, S. S., and Mathpati, C.S. 2009. Dyamics of flow structures and
22
transport phenomena-1: Experimental and numerical techniquues for identification and energy content of
23
flow structures. Industrial & engineering chemistry research 48, 8244-8248.
24
Joshi, J.B., Nandakumar, K., 2015. Computational Modeling of Multiphase Reactors. Annual Review of
25
Chemical and Biomolecular Engineering 6347–378
Ac ce
pt
ed
M
an
10
Page 100 of 100
Page 100 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Kuan, B., Yang, W., Schwarz, M., 2007. Dilute gas–solid two-phase flows in a curved 90∘ duct bend: CFD
2
simulation with experimental validation. Chemical Engineering Science 62, 2068-2088.
3
Kharaz, A., Gorham, D., Salman, A., 2001. An experimental study of the elastic rebound of spheres. Powder
4
Technology 120, 281-291.
5
Kim, J.S., Tachino, R., Tsutsumi, A., 2008. Effects of solids feeder and riser exit configuration on
6
establishing high density circulating fluidized beds. Powder Technology 187, 37-45.
7
Knowlton, T., Geldart, D., Matsen, J., King, D., 1995. Comparison of CFB hydrodynamic models. , PSRI
8
challenge problem at the Eighth international fludization conference, France.
9
Krambeck, F. J., 1991. Continuous Mixtures in Fluid Catalytic Cracking and Extensions. Mobil Workshop
an
us
cr
ip t
1
on Chemical Reaction in Complex Mixtures; Van Nostrand Reinhold: New York, 42-59.
11
Kulkarni, A.A., Joshi, J.B., Ravikumar, V and Kulkarni, B.D.,2001. Application f multi-resolution analysis
12
for the simultaneous measurement of gas and liquid velocities and fractional gas hold-up in bubble column
13
using LDA. Chemical Engineering Science 56, 5037-5048
14
Kunii, D., Levenspiel, O., 1968. Bubbling bed model for kinetic processes in fluidized beds. Gas-solid mass
15
and heat transfer and catalytic reactions. Industrial & Engineering Chemistry Process Design and
16
Development 7, 481-492.
17
Laı́n, S., Sommerfeld, M., Kussin, J., 2002. Experimental studies and modelling of four-way coupling in
18
particle-laden horizontal channel flow. International journal of heat and fluid flow 23, 647-656.
19
Lain, S., Garcia, J., 2006. Study of four-way coupling on turbulent particle-laden jet flows. Chemical
20
engineering science 61, 6775-6785.
21
Lan, X., Xu, C., Wang, G., Wu, L., Gao, J., 2009. CFD modeling of gas–solid flow and cracking reaction in
22
two-stage riser FCC reactors. Chemical Engineering Science 64, 3847-3858.
23
Lauder, B., Spalding, D., 1972. Mathematical models of turbulence. Academic Press NY.
Ac ce
pt
ed
M
10
Page 101 of 101
Page 101 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Lee, L.-S., Chen, Y.-W., Huang, T.-N., Pan, W.-Y., 1989. Four-lump kinetic model for fluid catalytic
2
cracking process. The Canadian Journal of Chemical Engineering 67, 615-619.
3
Li, J., Cheng, C., Zhang, Z., Yuan, J., Nemet, A., Fett, F.N., 1999. The EMMS model--its application,
4
development and updated concepts. Chemical Engineering Science 54, 5409-5425.
5
Li, J., Fan, Y.-P., Lu, C.-X., Luo, Z.-H., 2013. Numerical Simulation of Influence of Feed Injection on
6
Hydrodynamic Behavior and Catalytic Cracking Reactions in a FCC Riser under Reactive Conditions.
7
Industrial & Engineering Chemistry Research 52, 11084-11098.
8
Li, J., Kwauk, M., 1994. Particle-fluid two-phase flow: The energy-minimization multi-scale method.
9
Metallurgical Industry Press, Beijing, China.
us
cr
ip t
1
Li, J., Tung, Y., Kwauk, M., 1988. Axial voidage profiles of fast fluidized beds in different operating
11
regions. Circulating fluidized bed technology II, 193-203.
12
Li, P., Lan, X., Xu, C., Wang, G., Lu, C., Gao, J., 2009. Drag models for simulating gas–solid flow in the
13
turbulent fluidization of FCC particles. Particuology 7, 269-277.
14
Li, T., Benyahia, S., 2012. Revisiting Johnson and Jackson boundary conditions for granular flows. AIChE
15
Journal 58, 2058-2068.
16
Li, F., Song, F., Benyahia, S., Wang, W., Li, J., 2012. MP-PIC simulation of CFB riser with EMMS-based
17
drag model. Chemical Engineering Science 82, 104-113.
18
Li, T., Benyahia, S., 2013. Evaluation of wall boundary condition parameters for gas–solids fluidized bed
19
simulations. AIChE Journal 59, 3624-3632.
20
Li, T., Gel, A., Pannala, S., Shahnam, M., Syamlal, M., 2014a. Reprint of “CFD simulations of circulating
21
fluidized bed risers, part I: Grid study”. Powder Technology.
22
Li, T., Pannala, S., Shahnam, M., 2014b. CFD simulations of circulating fluidized bed risers, part II,
23
evaluation of differences between 2D and 3D simulations. Powder Technology 254, 115-124.
24
Liguras, D.K., Allen, D.T., 1989b. Structural models for catalytic cracking. 2. Reactions of simulated oil
25
mixtures. Industrial & Engineering Chemistry Research 28, 674-683.
Ac ce
pt
ed
M
an
10
Page 102 of 102
Page 102 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Liu, S., Chen, Q., Wang, H., Jiang, F., Ismail, I., Yang, W., 2005. Electrical capacitance tomography for
2
gas–solids flow measurement for circulating fluidized beds. Flow Measurement and Instrumentation 16, 135-
3
144.
4
Lopes, G.C., Rosa, L., Mori, M., Nunhez, J.R., Martignoni, W.P., 2011. Three-dimensional modeling of fluid
5
catalytic cracking industrial riser flow and reactions. Computers & Chemical Engineering 35, 2159-2168.
6
Louge, M., Chang, H., 1990. Pressure and voidage gradients in vertical gas-solid risers. Powder Technology
7
60, 197-201.
8
Lu, B., Wang, W., Li, J., 2009. Searching for a mesh-independent sub-grid model for CFD simulation of
9
gas–solid riser flows. Chemical Engineering Science 64, 3437-3447.
us
cr
ip t
1
Lun, C.K.K., Savage, S.B., 1984. Kinetic theories for granular ow: inelastic particles in Couette ow and
11
slightly inelastic particles in a general ow eld. Journal of Fluid Mechanics 140, 223-256.
12
Lun, C.K.K., Sevage, S.B., Jeffrey, D.J., Chepurniy, N., 1984. Kinetic theories for granular flow: inelastic
13
particles in Couette flow and slightly inelastic particles in a general flowfield. Journal of Fluid Mechanics,
14
223.
15
Mabrouk, R., Chaouki, J., Guy, C., 2007. Effective drag coefficient investigation in the acceleration zone of
16
an upward gas–solid flow. Chemical Engineering Science 62, 318-327.
17
Mangwandi, C., Cheong, Y., Adams, M., Hounslow, M., Salman, A., 2007. The coefficient of restitution of
18
different representative types of granules. Chemical Engineering Science 62, 437-450.
19
Mansoori, Z., Saffar-Avval, M., Tabrizi, H.B., Dabir, B., Ahmadi, G., 2005. Inter-particle heat transfer in a
20
riser of gas–solid turbulent flows. Powder technology 159, 35-45.
21
Martignoni, W., De Lasa, H., 2001. Heterogeneous reaction model for FCC riser units. Chemical
22
Engineering Science 56, 605-612.
23
Mathiesen, V., Solberg, T., Hjertager, B.H., 2000. An experimental and computational study of multiphase
24
flow behavior in a circulating fluidized bed. International Journal of Multiphase Flow 26, 387-419.
25
Mathpati, C.S., and Joshi, J.B., 2007. Insight into theories of heat and mass transfer at the solid-fkuid
26
interface using direct numerical simulation and large eddy simulation. Industrial & engineering chemistry
27
research 46, 8343-8354.
Ac ce
pt
ed
M
an
10
Page 103 of 103
Page 103 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Mathpati, C.S., Tabib, M.V., Deshpande, S.S., and Joshi, J.B. 2009. Dynamics of flow structures and
2
transport phenomena-2: Reltionship with design objectives and design optimization. Industrial &
3
Engineering Chemistry Research 48, 8285-8311.
4
Manyele, S., Pärssinen, J., Zhu, J.-X., 2002. Characterizing particle aggregates in a high-density and high-
5
flux CFB riser. Chemical Engineering Journal 88, 151-161.
6
McGreavy, C., Lu, M.L., Wang, X.Z., Kam, E.K.T., 1994. Characterisation of the behaviour and product
7
distribution in fluid catalytic cracking using neural networks. Chemical Engineering Science 49, 4717-4727.
8
Miller, A., Gidaspow, D., 1992. Dense, vertical gas-solid flow in a pipe. AIChE Journal 38, 1801-1815.
9
Miller, R.S., Harstad, K., Bellan, J., 1998. Evaluation of equilibrium and non-equilibrium evaporation
10
models for many-droplet gas-liquid flow simulations. International Journal of Multiphase Flow 24, 1025-
11
1055.
12
Monazam, E.R., Breault, R.W., Shadle, L.J., 2015. Pressure and apparent voidage profiles for riser with an
13
abrupt exit (T-shape) in a CFB riser operating above fast fluidization regimes. Powder Technology.
14
Morsi, S., Alexander, A., 1972. An investigation of particle trajectories in two-phase flow systems. Journal
15
of Fluid Mechanics 55, 193-208.
16
Mueller, P., Reh, L., 1993. Particle drag and pressure drop in accelerated gas-solid flow. Preprint Volume for
17
Circulating Fluidized Beds IV, AIChE, Somerset, 193–198.
18
Nace, D.M., Voltz, S., Weekman Jr, V., 1971. Application of a kinetic model for catalytic cracking. Effects
19
of charge stocks. Industrial & Engineering Chemistry Process Design and Development 10, 530-538.
20
Naren, P., Ranade, V.V., 2011. Scaling laws for gas–solid riser flow through two-fluid model simulation.
21
Particuology 9, 121-129.
22
Nayak, S.V., Joshi, S.L., Ranade, V.V., 2005. Modeling of vaporization and cracking of liquid oil injected in
23
a gas–solid riser. Chemical Engineering Science 60, 6049-6066.
24
Neri, A., Gidaspow, D., 2000. Riser hydrodynamics: simulation using kinetic theory. AIChE Journal 46, 52-
25
67.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 104 of 104
Page 104 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Nguyen, T.T., Mitra, S., Pareek, V., Joshi, J.B., Evans, G., 2015. Comparison of vaporization models for
2
feed droplet in fluid catalytic cracking risers. Chemical Engineering Research and Design 101, 82-97.
3
Nicolai, C., Jacob, B., Gualtieri, P., Piva, R., 2014. Inertial particles in homogeneous shear turbulence:
4
experiments and direct numerical simulation. Flow, Turbulence and Combustion 92, 65-82.
5
Nieuwland, J.J., Annaland, S., Kuipers, J.A.M., Swaaij, W.P.M., 1996a. Hydrodynamic modelling of gas-
6
particle flows in riser reactors. AIChE Journal 42, 1569-1582.
7
Nieuwland, J.J., Meijer, R., Kuipers, J.A.M., Van Swaaij, W.P.M., 1996b. Measurements of solids
8
concentration and axial solids velocity in gas-solid two-phase flows. Powder Technology 87, 127-139.
9
Nomura, H., Ujiie, Y., Rath, H.J., Sato, J.i., Kono, M., 1996. Experimental study on high-pressure droplet
10
evaporation using microgravity conditions, Symposium (International) on Combustion. Elsevier, pp. 1267-
11
1273.
12
Noymer, P.D., Glicksman, L.R., 1998. Cluster motion and particle-convective heat transfer at the wall of a
13
circulating fluidized bed. International Journal of Heat and Mass Transfer 41, 147-158.
14
Ocone, R., Sundaresan, S., Jackson, R., 1993. Gas‐Particle flow in a duct of arbitrary inclination with
15
particle‐particle interactions. AIChE Journal 39, 1261-1271.
16
Pachovsky, R., Wojciechowski, B., 1971. Theoretical interpretation of gas oil conversion data on an x‐sieve
17
catalyst. The Canadian Journal of Chemical Engineering 49, 365-369.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 105 of 105
Page 105 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Pandey, P., Turton, R., Yue, P., Shadle, L., 2004. Nonintrusive particle motion studies in the near-wall
2
region of a pilot-scale circulating fluidized bed. Industrial & Engineering Chemistry Research 43, 5582-
3
5592.
4
Paraskos, J.A., Shah, Y.T., McKinney, J.D., Carr, N.L., 1976. A kinematic model for catalytic cracking in a
5
transfer line reactor. Industrial & Engineering Chemistry Process Design and Development 15, 165-169.
6
Pärssinen, J.H., Zhu, J.X., 2001. Particle velocity and flow development in a long and high-flux circulating
7
fluidized bed riser. Chemical Engineering Science 56, 5295-5303.
8
Patience, G.S., Chaouki, J., 1993. Gas phase hydrodynamics in the riser of a circulating fluidized bed.
9
Chemical Engineering Science 48, 3195-3205.
us
cr
ip t
1
Pantzali, M., Bayón, N.L., Heynderickx, G., Marin, G., 2013. Three-component solids velocity
11
measurements in the middle section of a riser. Chemical Engineering Science 101, 412-423.
12
Peskin, C.S., 2002. The immersed boundary method. Acta numerica 11, 479-517.
13
Pita, J.A., Sundaresan, S., 1993. Developing flow of a gas-particle mixture in a vertical riser. AIChE Journal
14
39, 541-552.
15
Pitault, I., Forissier, M., Bernard, J.R., 1995. Determination de constantes cinetiques du craquage catalytique
16
par la modelisation du test de microactivite (MAT). The Canadian Journal of Chemical Engineering 73, 498-
17
504.
18
Pitault, I., Nevicato, D., Forissier, M., Bernard, J.-R., 1994. Kinetic model based on a molecular description
19
for catalytic cracking of vacuum gas oil. Chemical Engineering Science 49, 4249-4262.
20
Prevel, M., Vinkovic, I., Doppler, D., Pera, C., Buffat, M., 2013. Direct numerical simulation of particle
21
transport by hairpin vortices in a laminar boundary layer. International Journal of Heat and Fluid Flow 43, 2-
22
14.
23
Pugsley, T.S., Berruti, F., 1996. A predictive hydrodynamic model for circulating fluidized bed risers.
24
Powder Technology 89, 57-69.
25
Ranz, W., Marshall, W., 1952. Evaporation from drops. Chemical Engineering Progress 48, 141-446.
Ac ce
pt
ed
M
an
10
Page 106 of 106
Page 106 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Renksizbulut, M., Bussmann, M., 1993. Multicomponent droplet evaporation at intermediate Reynolds
2
numbers. International Journal of Heat and Mass Transfer 36, 2827-2835.
3
Rhodes, M.J., Geldart, D., 1986. The hydrodynamics of re-circulating fluidized beds. Circulating Fluidized
4
Bed Technology, 193-200.
5
Richardson, J.F., Zaki, W.N., 1954. Sedimentation and fluidization. Trans. Inst. Chem. Eng 32, 1954.
6
Rodrigues, S.S., Forret, A., Montjovet, F., Lance, M., Gauthier, T., 2015. Riser hydrodynamic study with
7
different Group B powders. Powder Technology 272, 300-308.
8
Samuelsberg, A., Hjertager, B.H., 1996. Computational modeling of gas/particle flow in a riser. AIChE
9
Journal 42, 1536-1546.
us
cr
ip t
1
Sadeghbeigi, R., 2012. Fluid catalytic cracking handbook: An expert guide to the practical operation, design,
11
and optimization of FCC units. Elsevier.
12
Schaeffer, D.G., 1987. Instability in the evolution equations describing incompressible granular flow. Journal
13
of differential equations 66, 19-50.
14
Shaffer, F., Gopalan, B., Breault, B., Cocco, R., Hays, R., Karri, R., Knowlton, T., 2010. A New View of
15
Riser Flow Fields Using High Speed Particle Imaging Velocimetry (PIV), 2010 AIChE Annual Meeting.
16
AIChE, Salt Lake City, UT, US.
17
Shaffer, F., Gopalan, B., Breault, R.W., Cocco, R., Karri, S., Hays, R., Knowlton, T., 2013. High speed
18
imaging of particle flow fields in CFB risers. Powder Technology.
19
Shah, M.T., Utikar, R.P., Evans, G.M., Tade, M.O., Pareek, V.K., 2011a. Effect of Inlet Boundary
20
Conditions on Computational Fluid Dynamics (CFD) Simulations of Gas–Solid Flows in Risers. Industrial &
21
Engineering Chemistry Research 51, 1721-1728.
22
Shah, M.T., Utikar, R.P., Pareek, V.K., Tade, M.O., Evans, G.M., 2015. Effect of Closure Models on
23
Eulerian-Eulerian Gas-Solid Flow Predictions in Riser. Powder Technology 269, 247-258
24
Shah, M.T., Utikar, R.P., Tade, M.O., Evans, G., Pareek, V.K., 2011b. Simulation of gas-solid flows in riser
25
using energy minimization multiscale model: Effect of cluster diameter correlation. Chemical Engineering
26
Science 66, 3391-3300.
Ac ce
pt
ed
M
an
10
Page 107 of 107
Page 107 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Shah, M.T., Utikar, R.P., Tade, M.O., Evans, G.M., Pareek, V.K., 2013. Effect of a cluster on gas–solid drag
2
from lattice Boltzmann simulations. Chemical Engineering Science 102, 365-372.
3
Shah, M.T., Utikar, R.P., Tade, M.O., Pareek, V.K., 2011c. Hydrodynamics of an FCC riser using energy
4
minimization multiscale drag model. Chemical engineering journal 168, 812-821.
5
Shah, Y., Huling, G., Paraskos, J., McKinney, J., 1977. A kinematic model for an adiabatic transfer line
6
catalytic cracking reactor. Industrial & Engineering Chemistry Process Design and Development 16, 89-94.
7
Sharma, A.K., Tuzla, K., Matsen, J., Chen, J.C., 2000. Parametric effects of particle size and gas velocity on
8
cluster characteristics in fast fluidized beds. Powder Technology 111, 114-122.
9
Shi, X., Wu, Y., Lan, X., Liu, F., Gao, J., 2015. Effects of the riser exit geometries on the hydrodynamics
10
and solids back-mixing in CFB risers: 3D simulation using CPFD approach. Powder Technology 284, 130-
11
142.
12
Shi, X., Sun, R., Lan, X., Liu, F., Zhang, Y., Gao, J., 2015. CPFD simulation of solids residence time and
13
back-mixing in CFB risers. Powder Technology 271, 16-25.
14
Shnip, A.I., Kolhatkar, R.V., Swamy, D., Joshi, J.B., 1992. Criteria for the transition from the homogeneous
15
to the heterogeneous regime in two-dimensional bubble column reactors. International Journal of Multiphase
16
Flow 18, 705-726.
17
Shuai, W., Huilin, L., Guodong, L., Zhiheng, S., Pengfei, X., Gidaspow, D., 2011. Modeling of cluster
18
structure-dependent drag with Eulerian approach for circulating fluidized beds. Powder Technology 208, 98-
19
110.
20
Simonin, C., Viollet, P., 1990. Predictions of an oxygen droplet pulverization in a compressible subsonic
21
coflowing hydrogen flow. Numerical Methods for Multiphase Flows, FED91, 65-82.
22
Sinclair, J.L., Jackson, R., 1989. Gas particle flow in a vertical pipe with particle particle interactions. AIChE
23
Journal 35, 1473-1486.
24
Smagorinsky, J., 1963. General circulation experiments with the primitive equations: I. the basic
25
experiment*. Monthly Weather Review 91, 99-164.
26
Smolders, K., Baeyens, J., 2000. Overall solids movement and solids residence time distribution in a CFB-
27
riser. Chemical Engineering Science 55, 4101-4116.
Ac ce
pt
ed
M
an
us
cr
ip t
1
Page 108 of 108
Page 108 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Snider, D., 2001. An incompressible three-dimensional multiphase particle-in-cell model for dense particle
2
flows. Journal of Computational Physics 170, 523-549.
3
Soong, C., Tuzla, K., Chen, J., 1994. Identification of particle clusters in circulating fluidized bed.
4
Circulating Fluidized Bed Technology IV, 615-620.
5
Souza, J., Vargas, J., Von Meien, O., Martignoni, W., Amico, S., 2006. A two‐dimensional model for
6
simulation, control, and optimization of FCC risers. AIChE journal 52, 1895-1905.
7
Sugungun, M., Kolesnikov, I., Vinogradov, V., Kolesnikov, S., 1998. Kinetic modeling of FCC process.
8
Catalysis Today 43, 315-325.
9
Sundaresan, S., 2003. Instabilities in fluidized beds. Annual Review of Fluid Mechanics 35, 63-88.
M
an
us
cr
ip t
1
Syamlal, M., O’Brien, T.J., 1987. Derivation of a drag coefficient from velocity-voidage correlation. US
11
Dept. of Energy, Office of Fossil Energy, National Energy Technology Laboratory, Morgantown, West
12
Virginia April.
13
Syamlal, M., Rogers, W., O'Brien, T.J., 1993. MFIX documentation: Theory guide. Technical Note,
14
DOE/METC-95/1013.
15
Syamlal, M., O'Brien, T.J., 2003. Fluid dynamic simulation of O3 decomposition in a bubbling fluidized bed.
16
AIChE Journal 49, 2793-2801.
17
Tenneti, S., Garg, R., Subramaniam, S., 2011. Drag law for monodisperse gas–solid systems using particle-
18
resolved direct numerical simulation of flow past fixed assemblies of spheres. International Journal of
19
Multiphase Flow 37, 1072-1092.
20
Thakare, S.S., and Joshi, J.B., 1999. CFD simulation of flow in bubble column reactors importance of drag
21
force formulation. Chemical Engineering Science 54, 5055-5060.
22
Theologos, K.N., Markatos, N.C., 1993. Advanced modeling of fluid catalytic cracking riser-type reactors.
23
AIChE Journal 39, 1007-1017.
Ac ce
pt
ed
10
Page 109 of 109
Page 109 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Theologos, K.N., Nikou, I.D., Lygeros, A.I., Markatos, N.C., 1997. Simulation and design of fluid catalytic-
2
cracking riser-type reactors. AIChE Journal 43, 486-494.
3
Theologos, K.N., Lygeros, A.I., Markatos, N.C., 1999. Feedstock atomization effects on FCC riser reactors
4
selectivity. Chemical Engineering Science 54, 5617-5625.
5
Thorat, B.N., Joshi, J.B., 2004. Regime transition in bubble columns: experimental and predictions.
6
Experimental and Thermal and Fluid Science 28, 423-430.
7
Thorat, B.N., Shevade, A.V., Bhilegaonkar, K.N., Aglawe, R., Parasu Veera, U., Thakre, S.S., Pandit, A.B.,
8
Sawant, S.B., Joshi, J.B, 1998. Effect of sparger design and height to diameter ratio on fractional gas hold-up
9
in bubble columns. Chemical Engineering Research and Design 76, 823-834.
us
cr
ip t
1
Tsuji, Y., Kawaguchi, T., Tanaka, T., 1993. Discrete particle simulation of two-dimensional fluidized bed.
11
Powder Technology 77, 79-87.
12
Trujillo, W.R., De Wilde, J., 2012. Fluid catalytic cracking in a rotating fluidized bed in a static geometry: a
13
CFD analysis accounting for the distribution of the catalyst coke content. Powder Technology 221, 36-46.
14
Uhlmann, M., 2005. An immersed boundary method with direct forcing for the simulation of particulate
15
flows. Journal of Computational Physics 209, 448-476.
16
Vegendla, S.P., Heynderickx, G.J., Marin, G.B., 2011. Comparison of Eulerian–Lagrangian and Eulerian–
17
Eulerian method for dilute gas–solid flow with side inlet. Computers & Chemical Engineering 35, 1192-
18
1199.
19
Van de Velden, M., Baeyens, J., Smolders, K., 2007. Solids mixing in the riser of a circulating fluidized bed.
20
Chemical Engineering Science 62, 2139-2153.
21
van der Hoef, M., Beetstra, R., Kuipers, J., 2005. Lattice-Boltzmann simulations of low-Reynolds-number
22
flow past mono-and bidisperse arrays of spheres: results for the permeability and drag force. Journal of Fluid
23
Mechanics 528, 233-254.
24
van Wachem, B.G.M., Schouten, J.C., Van den Bleek, C.M., Krishna, R., Sinclair, J.L., 2001. Comparative
25
analysis of CFD models of dense gas–solid systems. AIChE Journal 47, 1035-1051.
Ac ce
pt
ed
M
an
10
Page 110 of 110
Page 110 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Voltz, S.E., Nace, D.M., Weekman Jr, V.W., 1971. Application of a kinetic model for catalytic cracking.
2
Some correlations of rate constants. Industrial & Engineering Chemistry Process Design and Development
3
10, 538-541.
4
Vreman, B., Geurts, B.J., Deen, N., Kuipers, J., Kuerten, J., 2009. Two-and four-way coupled Euler–
5
Lagrangian large-eddy simulation of turbulent particle-laden channel flow. Flow, Turbulence and
6
Combustion 82, 47-71.
7
Wang, J., Ge, W., Li, J., 2008. Eulerian simulation of heterogeneous gas–solid flows in CFB risers: EMMS-
8
based sub-grid scale model with a revised cluster description. Chemical Engineering Science 63, 1553-1571.
9
Wang, J., 2008. High-resolution Eulerian simulation of RMS of solid volume fraction fluctuation and
us
cr
ip t
1
particle clustering characteristics in a CFB riser. Chemical Engineering Science 63, 3341-3347.
11
Wang, J., 2009. A review of Eulerian simulation of Geldart A particles in gas-fluidized beds. Industrial &
12
Engineering Chemistry Research 48, 5567-5577.
13
Wang, Q., Niemi, T., Peltola, J., Kallio, S., Yang, H., Lu, J., Wei, L., 2014. Particle size distribution in
14
CPFD modeling of gas–solid flows in a CFB riser. Particuology.
15
Wang, W., Lu, B., Zhang, N., Shi, Z., Li, J., 2010. A review of multiscale CFD for gas–solid CFB modeling.
16
International Journal of Multiphase Flow 36, 109-118.
17
Wang, X., Jiang, F., Xu, X., Fan, B., Lei, J., Xiao, Y., 2010. Experiment and CFD simulation of gas–solid
18
flow in the riser of dense fluidized bed at high gas velocity. Powder Technology 199, 203-212.
19
Weekman Jr, V., 1968. Model of catalytic cracking conversion in fixed, moving, and fluid-bed reactors.
20
Industrial & Engineering Chemistry Process Design and Development 7, 90-95.
21
Weekman Jr, V.W., 1969. Kinetics and dynamics of catalytic cracking selectivity in fixed-bed reactors.
22
Industrial & Engineering Chemistry Process Design and Development 8, 385-391.
23
Weekman, V.W., Nace, D.M., 1970. Kinetics of catalytic cracking selectivity in fixed, moving, and fluid bed
24
reactors. AIChE Journal 16, 397-404.
25
Wei, F., Lin, H., Cheng, Y., Wang, Z., Jin, Y., 1998. Profiles of particle velocity and solids fraction in a
26
high-density riser. Powder Technology 100, 183-189.
Ac ce
pt
ed
M
an
10
Page 111 of 111
Page 111 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Wei, F., Yang, G.Q., Jin, Y., Yu, Z.Q., 1995. The characteristics of cluster in a high density circulating
2
fluidized bed. The Canadian Journal of Chemical Engineering 73, 650-655.
3
Wei, M., Wang, L., Li, J., 2013. Unified stability condition for particulate and aggregative fluidization—
4
Exploring energy dissipation with direct numerical simulation. Particuology 11, 232-241.
5
Wen, C.Y., Yu, Y.H., 1966. Mechanics of fluidization, Chemical Engineering Progress Symposium Series
6
62, 100–111..
7
Wojciechowski, B., 1968. A theoretical treatment of catalyst decay. The Canadian Journal of Chemical
8
Engineering 46, 48-52.
9
Wong, S.-C., Lin, A.-C., 1992. Internal temperature distributions of droplets vaporizing in high-temperature
us
cr
ip t
1
convective flows. Journal of fluid Mechanics 237, 671-687.
11
Wu, C., Cheng, Y., Ding, Y., Jin, Y., 2010. CFD–DEM simulation of gas–solid reacting flows in fluid
12
catalytic cracking (FCC) process. Chemical Engineering Science 65, 542-549.
13
Xu, G., Li, J., 1998. Analytical solution of the energy-minimization multi-scale model for gas-solid two-
14
phase flow. Chemical Engineering Science 53, 1349-1366.
15
Xu, J., Zhu, J.X., 2011. Visualization of particle aggregation and effects of particle properties on cluster
16
characteristics in a CFB riser. Chemical Engineering Journal 168, 376-389.
17
Xu, Y., Subramaniam, S., 2010. Effect of particle clusters on carrier flow turbulence: A direct numerical
18
simulation study. Flow, Turbulence and Combustion 85, 735-761.
19
Yamamoto, Y., Potthoff, M., Tanaka, T., Kajishima, T., Tsuji, Y., 2001. Large-eddy simulation of turbulent
20
gas–particle flow in a vertical channel: effect of considering inter-particle collisions. Journal of fluid
21
Mechanics 442, 303-334.
22
Yan, C., Fan, Y., Lu, C., Zhang, Y., Liu, Y., Cao, R., Gao, J., Xu, C., 2009. Solids mixing in a fluidized bed
23
riser. Powder Technology 193, 110-119.
24
Yang, Y.-L., Jin, Y., Yu, Z.-Q., Wang, Z.-W., 1992. Investigation on slip velocity distributions in the riser of
25
dilute circulating fluidized bed. Powder technology 73, 67-73.
Ac ce
pt
ed
M
an
10
Page 112 of 112
Page 112 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Yang, N., Wang, W., Ge, W., Li, J., 2003a. Analysis of flow structure and calculation of drag coefficient for
2
concurrent-up gas-solid flow. Chinese Journal of Chemical Engineering 11, 79-84.
3
Yang, N., Wang, W., Ge, W., Li, J., 2003b. CFD simulation of concurrent-up gas-solid flow in circulating
4
fluidized beds with structure-dependent drag coefficient. Chemical Engineering Journal 96, 71-80.
5
Yang, B.l., Zhou, X.W., Yang, X.H., Chen, C., Wang, L.Y., 2009. Multi‐scale study on the secondary
6
reactions of fluid catalytic cracking gasoline. AIChE journal 55, 2138-2149.
7
Yen, L., Wrench, R., Ong, A., 1988. Reaction kinetic correlation equation predicts fluid catalytic cracking
8
coke yields. Oil & Gas Journal;(United States) 86.
9
Yerushalmi, J., Cankurt, N., Geldart, D., Liss, B., 1976 Flow regimes in vertical gas-solid contact systems.
an
us
cr
ip t
1
AIChE Symposium Series. 69, 1-13..
11
Yerushalmi, J., Cankurt, N., 1979. Further studies of the regimes of fluidization. Powder Technology 24,
12
187-205.
13
Yin, X., Sundaresan, S., 2009. Fluid‐particle drag in low‐Reynolds‐number polydisperse gas–solid
14
suspensions. AIChE Journal 55, 1352-1368.
15
Yamamoto, Y., Potthoff, M., Tanaka, T., Kajishima, T., Tsuji, Y., 2001. Large-eddy simulation of turbulent
16
gas–particle flow in a vertical channel: effect of considering inter-particle collisions. Journal of Fluid
17
Mechanics 442, 303-334.
18
Zuber, N., Findlay, J., 1965. Average volumetric concentration in two-phase flow systems. Journal of Heat
19
Transfer 87, 453-468.
20
Zhang, M., Qian, Z., Yu, H., Wei, F., 2003. The solid flow structure in a circulating fluidized bed
21
riser/downer of 0.42-m diameter. Powder Technology 129, 46-52.
Ac ce
pt
ed
M
10
Page 113 of 113
Page 113 of 114
Draft manuscript: - Computational Flow Modelling of FCC Riser: A Review
Zhang, M., Chu, K., Wei, F., Yu, A., 2008. A CFD–DEM study of the cluster behavior in riser and downer
2
reactors. Powder Technology 184, 151-165.
3
Zhang, Y., Ge, W., Wang, X., Yang, C., 2011. Validation of EMMS-based drag model using lattice
4
Boltzmann simulations on GPUs. Particuology 9, 365-373.
5
Zhang, R., Yang, H., Wu, Y., Zhang, H., Lu, J., 2013. Experimental study of exit effect on gas–solid flow
6
and heat transfer inside CFB risers. Experimental Thermal and Fluid Science 51, 291-296.
7
Zhou, H., Flamant, G., Gauthier, D., Lu, J., 2002. Lagrangian approach for simulating the gas-particle flow
8
structure in a circulating fluidized bed riser. International Journal of Multiphase Flow 28, 1801-1821.
9
Zhou, H., Flamant, G., Gauthier, D., 2007. Modelling of the turbulent gas–particle flow structure in a two-
us
cr
ip t
1
dimensional circulating fluidized bed riser. Chemical Engineering Science 62, 269-280.
11
Zhao, Y., Cheng, Y., Wu, C., Ding, Y., Jin, Y., 2010. Eulerian–Lagrangian simulation of distinct clustering
12
phenomena and RTDs in riser and downer. Particuology 8, 44-50.
13
Zhou, Q., Wang, J., 2014. CFD study of mixing and segregation in CFB risers: Extension of EMMS drag
14
model to binary gas-solid flow. Chemical Engineering Science.
15
Zhu, H., Zhou, Z., Yang, R., Yu, A., 2008. Discrete particle simulation of particulate systems: a review of
16
major applications and findings. Chemical Engineering Science 63, 5728-5770.
17
Zhu, C., Jun, Y., Patel, R., Wang, D., Ho, T.C., 2011. Interactions of flow and reaction in fluid catalytic
18
cracking risers. AIChE Journal 57, 3122-3131.
19
Khan, Z., 2016. Flow visualization and CFD modelling in multiphase reactors. Ph.D. Thesis, Institute of
20
Chemical Technology, Mumbai, India.
21
Zuber, N., Findlay, J., 1965. Average volumetric concentration in two-phase flow systems. Journal of Heat
22
Transfer 87, 453-468.
Ac ce
pt
ed
M
an
10
Page 114 of 114
Page 114 of 114