Research Article  Open Access
Shitanshu Gohel, Shitalkumar Joshi, Mohammed Azhar, Marc Horner, Gustavo Padron, "CFD Modeling of Solid Suspension in a Stirred Tank: Effect of Drag Models and Turbulent Dispersion on Cloud Height", International Journal of Chemical Engineering, vol. 2012, Article ID 956975, 9 pages, 2012. https://doi.org/10.1155/2012/956975
CFD Modeling of Solid Suspension in a Stirred Tank: Effect of Drag Models and Turbulent Dispersion on Cloud Height
Abstract
Many chemical engineering processes involve the suspension of solid particles in a liquid. In dense systems, agitation leads to the formation of a clear liquid layer above a solid cloud. Cloud height, defined as the location of the clear liquid interface, is a critical measure of process performance. In this study, solidliquid mixing experiments were conducted and cloud height was measured as a function operating conditions and stirred tank configuration. Computational fluid dynamics simulations were then performed using an EulerianGranular multiphase model. The effects of hindered and unhindered drag models and turbulent dispersion force on cloud height were investigated. A comparison of the experimental and computational data showed excellent agreement over the full range of conditions tested.
1. Introduction
The suspension of solid particles in a liquid is a key requirement of many industrial processes. Examples relevant to this work include crystallisation, dissolution, and adsorption processes. Each of these can be characterized as complex multiphase processes that are often facilitated by mechanically stirred tanks. Maximizing contact between the solid and liquid phases facilitates mass transfer and reaction, therefore, assessing the ability of the process to suspend particles is a key objective. The quality of suspension can be quantified by three parameters: just suspension velocity, solids distribution, and cloud height. The latter is defined as the location of the interface between solidsrich liquid and clear liquid regions. Cloud height is critical because of limited mixing between the suspended solids and the upper clear liquid layer. The current study aims to develop and validate a computational model of a solids suspension process using experimental measurements of cloud height.
The hydrodynamics in stirred vessels are complex, threedimensional, and turbulent. The interplay between solid and liquid materials, vessel design, impeller type and location, and level of agitation all determine the efficiency of the mixing process. Numerous empirical models have been developed to relate tank performance to operating conditions and geometry [1, 2]. However, these models lack wider application, owing to the complexity of the flow field. This has helped to motivate the need for more detailed analysis.
Computational fluid dynamics (CFD) methods are powerful tools for improving our understanding of mixing in stirred tanks. In the last few years, encouraging results have been obtained in the simulation of solid suspension systems [3]. The solids distribution has been traditionally simulated using either Lagrangian [4–6] or Eulerian multiphase modeling approaches [7–14]. The Lagrangian model treats solids as individual particles subject to hydrodynamic forces that are approximated using singleparticle empirical models. The Eulerian model, on the other hand, treats the solids as an additional continuous phase. Interphase transport models are added to account for the interaction between the solids and the fluid phase.
The Granular model is an extension of the Eulerian model wherein the solids are modeled as a pseudofluid. The pseudofluid physical properties, such as viscosity, solids pressure and stress, are derived from the kinetic theory of granular flow. According to Montante et al. [15], the Granular model has been found to be more appropriate for modeling processes exhibiting high local solid concentrations.
Turbulence modeling is also critical to the success of CFD analysis in solidliquid systems. Turbulence plays an important role in developing singlephase flow structures, and it is generally accepted that it will also have some influence on the solids distribution in stirred tanks. Previous CFD investigations have highlighted the importance of turbulenceassisted solids dispersion in stirred tanks [16–20]. However, a quantitative evaluation of the effect of turbulent dispersion on cloud height has not been considered in detail.
To address the above, this study aims to develop a CFD modeling strategy for predicting cloud height in a solid suspension process. An EulerianGranular multiphase model is used to model the particle suspension process. Experimental data is used to justify the selection of model inputs, specifically the drag model and turbulent dispersion force. Each of these were evaluated independently to illustrate their relative importance. Finally, a comparison of the model results with corresponding experimental data establishes the validity of the modeling approach across a wide range of experimental conditions.
2. Experimental Setup and Procedure
The experimental work was carried out at BHR Group as part of the fluid mixing processes (FMPs) industrial research program. The tests were performed in a transparent Perspex tank with a diameter of 0.61 m, equipped with four standard baffles and a torispherical base, as shown in Figure 1. The impellers were centrally mounted pitched blade turbines (PBTs) with a diameter of 0.2 m and a blade width of 0.055 m. One set of tests was carried out with a singleimpeller located at a distance from the tank bottom of 0.15 m. For the dual impeller tests, a second impeller was installed on the same shaft at an offbottom distance of 0.39 m. The liquid height was kept constant at 1.5 in all experiments.
The solid phase was comprised of sand particles with a nominal average diameter of 180 m and a density of 2,630 kg/m^{3}. The liquid phase was water with a density of 1,000 kg/m^{3} and a viscosity of 0.001 Pas. Solids concentrations of 10% (w/w) and 15% (w/w) were tested. These values correspond to 4% and 6.2% by volume, respectively.
The solids cloud height was obtained by measuring the distance from the bottom of the tank to the cloud surface. Measurements were determined visually from outside the tank. Given the turbulent nature of the flow, the cloud surface was not flat. Therefore, the minimum and maximum positions of the surface at a fixed radial position on the wall of the tank were recorded. The average cloud height could then be calculated for each experimental condition. Measurements were carried out at impeller speeds ranging from 150 to 450 RPM in 50 RPM increments.
For both impeller systems, the minimum impeller speed for just suspension was also determined at each solids loading. This is the minimum speed at which the impeller(s) operate to ensure that the solids are suspended, that is, not settled at the vessel bottom. was determined by placing a mirror underneath the tank and observing the bed of solids at the bottom while the impeller speed was increased in small increments. The criterion used to determine whether the solids were suspended followed that of Zwietering [1], who defined as “the impeller speed at which no material remains stationary on the vessel base for more than 12 seconds.”
Figure 2 shows the results for the fractional average cloud height versus impeller speed normalized by . At impeller speeds below , the cloud height is considerably higher with the dual impeller system than with the single impeller system (up to about 70% higher). The difference becomes less significant at higher impeller speeds.
With the singleimpeller system, there is no significant difference between the two different solid concentrations. There is a small difference at high impeller speeds with the dualimpeller system. However, this could be due to experimental error (typically 5 to 10%).
3. Modeling and Numerical Simulations
3.1. EulerGranular Theory
The simulations were performed using the EulerGranular multiphase model available in the commercial CFD software, ANSYS Fluent 13.0 (ANSYS Inc., Canonsburg, PA, USA). This model assumes that each phase coexists at every point in space in the form of interpenetrating continua. The continuity and momentum equations are solved for all phases and the coupling between phases is obtained through pressure and interphase exchange coefficients.
The conservation equation for each phase is [21]: where terms with the subscript “” account for exchange between the fluid and solid phases. One fluid phase and one solid phase were modeled in this work. They are denoted individually with subscripts and , respectively, in the rest of this section.
In the EulerianGranular model, the solid phase momentum equation includes an additional solid pressure term. The solids pressure and solids stress tensor terms are determined from the kinetic theory of granular flow as described by Ding and Gidaspow [22]. The EulerGranular model also includes an additional transport equation for the granular temperature.
In (2), accounts for the effects of lift force, virtual mass force and any other external forces acting on the solid phase. The term represents fluid solid interphase momentum exchange which is expressed as: where is the interphase exchange coefficient, which is defined in terms of the drag model proposed by Ding and Gidaspow [22]: For turbulent flows, the drag force is a function of average velocity. A correction term is included in the interphase momentum exchange equation to account for the effect of turbulent dispersion due to instantaneous fluctuations, as proposed by Deutsch and Simonin [23]. The resulting interphase momentum exchange term is then given by: Here represents phasic average velocities and is the drift velocity: In this expression, are the fluidsolid turbulent dispersion coefficients, and is the dispersion Prandtl number (which is set equal to 0.75 in this study).
3.2. Simulation Setup
3.2.1. Geometry and Mesh
The geometric models and meshes were constructed in the commercial software, GAMBIT 2.4 (ANSYS, Inc., Canonsburg, PA, USA). Full hexahedral meshes were constructed for each of the two tank configurations (Figure 3). The computational mesh for the single impeller tank had 47, 72, and 122 cells in the radial, angular, and axial directions, respectively; while the mesh for the dual impeller configuration had 47, 72, and 177 cells in the respective directions. These settings resulted in 402k hexahedral elements for the singleimpeller tank and 576k hexahedral elements for the dualimpeller tank. (Note that multiplying the mesh sizings in each direction will result in a mesh count that is higher than the final mesh count quoted in the text. This is because there is no mesh at impeller and shaft locations and because there are small regions of paved mesh which are not regularly spaced.) The mesh spacing was more or less uniform in each direction, with mesh refinement on wall surfaces and in the impeller zones.
A griddependency study was carried out to evaluate mesh suitability for the single impeller configuration. The procedure involved comparing results from the baseline mesh to a refined grid with 3.2 million hexahedral cells, representing an almost 8fold increase over the baseline mesh. The effects on one flow parameter—impeller power number—and one multiphase model output—cloud height—were extracted and compared (see Table 1). The Gidaspow drag model and turbulent dispersion force were included in all cases. The maximum difference in the impeller power numbers was less than 4%, and the maximum difference in the predicted cloud height was less than 2%. Therefore, the resolution of the baseline mesh was deemed satisfactory.

3.2.2. Physical Models and Boundary Conditions
The realizable kε turbulence model was used to resolve the turbulent solidliquid flow field. For comparison purposes, interphase momentum exchange terms were calculated using the Gidaspow (hindered) and SchillerNaumann (unhindered) drag models. The effect of including turbulenceassisted solids dispersion in the interphase exchange term was also evaluated [23, 24].
The multiple reference frame (MRF) approach was used to model impeller rotation. Although an approximation, this approach has been shown to produce satisfactory results, especially when the interaction between baffles and impellers is weak. The free surface was assumed to be fixed, thus a symmetry boundary condition was prescribed. A noslip boundary condition was applied for all other wall surfaces.
Highorder discretization schemes were chosen to minimize numerical diffusion. The GreenGauss Node Based Gradient option was used for gradient calculations at cell interfaces and the QUICK discretization scheme was used for momentum, volume fraction and turbulence equations. All cases were initialized with a uniformly distributed solid phase. The steady state solver was utilized to solve for all flow variables.
4. Results
4.1. Flow and Solid Distributions
Liquid velocities provide insight into the mixing processes occurring inside stirred tanks. Figures 4 and 5 show model results for liquid velocities and solid phase volume fraction on a tank midplane. Recirculation zones are evident near the tank bottom. These recirculation zones are predicted to be larger in the dual impeller configuration and increase in size with RPM. A distinct interface between clear liquid at the top of the tank and the cloudy region is observed at low RPM. The model also predicts that solids are more uniformly distributed throughout the entire reactor at higher RPM.
(a)
(b)
(a)
(b)
4.2. Effect of Drag
Drag forces induce significant interphase momentum transfer and tend to dominate other interphase processes. Therefore, drag force model selection is critical to the accuracy of the solidliquid suspension model. The Gidaspow and SchillerNaumann drag models were investigated. The Gidaspow drag model has been specifically tailored to take into consideration particle induced hindrance. For this reason, it is expected to give good predictions for high solid volume fractions (as was observed under certain operating conditions studied here). Additional simulations with the SchillerNaumann model were performed to help elucidate the importance of particle hindrance in this system.
Figure 6 shows volume rendered (volume rendered postprocessing permits visualization of the solidphase distribution throughout the entire domain by increasing the opacity of the foreground results) contour plots of solids concentration and cloud height for the single and the dual impeller configurations. The doubleheaded arrows marked “EXP” represent the location of the experimentally observed cloud height with the arrow thickness representing the observed ±8% variation in the experimental data. When postprocessing the CFD model, the cloud height was defined as the maximum coordinate of an isosurface of solid phase volume fraction (note that the shaft axis is aligned with axis). Following Kasat et al. [12], the location of this isosurface is defined as the average solid phase volume fraction. It can be seen from Figure 6 that the Gidaspow model is a better predictor of cloud height versus the SchillerNaumann model. Experimental results further substantiate the Gidaspow drag model as a better predictor of cloud height (see Figure 7). These results suggest that particle hindrance effects are an important aspect of the model, even at the relatively low solids loading studied here.
4.3. Effect of the Turbulent Dispersion Force (TDF)
The contribution of the turbulent dispersion force is significant when the size of the turbulent eddies is larger than the particle size [12]. Figure 8 shows how the inclusion of TDF substantially increases the cloud height. Moreover, the predicted cloud heights give better comparison with the experimental data when TDF is included. In an extreme example, it was observed that omitting TDF leads to a 250% difference in the cloud height prediction (see Figure 8(b)). Similar studies performed with the SchillerNaumann drag model showed an identical trend, although the predicted cloud heights were observed to be slightly lower than that predicted by the Gidaspow. Therefore, TDF must be included in the model for the reactor configuration and conditions studied here.
4.4. Experimental Validation
Results from the previous sections establish the model requirements based on a limited set of experiments. Applying the modeling strategy over the full range of experimental operating conditions establishes model validation. The experimental study included single and dualimpeller tanks operating at six to seven speeds and for two different solids concentrations (10% (w/w) and 15% (w/w)), resulting in a total of 26 experimental measurements. Numerical simulations were performed for each of these experimental points.
The profiles of the normalized cloud heights are plotted in Figures 9(a) and 9(b). It can be seen that the predicted cloud height for all operating conditions are in good agreement with the experimental data. For the singleimpeller system, the simulations overpredict cloud height at lower RPM for both 10% (w/w) and 15% (w/w) solids loading. This may be due to the fact that the single impeller may be not able to fully suspend solids at low RPM. Modifications to the TDF model to include effect of local solid concentration may overcome this difficulty. This analysis is left for future study.
(a)
(b)
5. Conclusion
A CFD modeling strategy was developed to predict cloud height in stirred tanks with single and dualaxial PBT impellers. The CFD model was based on EulerianGranular multiphase theory. The effects of two different drag models, Gidaspow and SchillerNaumann, and a drift velocitybased turbulent dispersion force were analyzed. The modeling strategy was developed for a limited set of experimental conditions and then applied across the full range solids loading, agitation rates, and reactor configurations. It was determined that the inclusion of turbulent dispersion force was critical to accurate prediction of cloud height for the conditions studied. The drag model was observed to have a much lower influence. The resulting modeling strategy was then compared to the full experimental range of cloud height measurements. The model was found to be predictive of cloud height across a broad range of experimental conditions, the exception being the singleimpeller configuration at low agitation rates. It is the authors’ opinion that solids concentrationbased inclusion of the TDF can be studied further to improve the model prediction at lower agitation rates.
Nomenclature
:  Distance of lower impeller from the tank bottom, m 
:  Distance of upper impeller from the tank bottom, m 
:  Drag coefficient 
:  Impeller diameter, m 
:  Particle diameter, m 
:  fluidsolid turbulent dispersion coefficient, 
:  Lift, virtual mass, and any other external forces acting on the solid phase, 
:  Liquid fill level, m 
:  Cloud height, m 
:  Interphase exchange coefficient, 
:  Fluidsolid interphase mass exchange, 
:  Impeller rotational speed, rpm 
:  Just suspension speed, rpm 
:  Pressure shared by all phases, 
Re:  Relative Reynolds number 
:  Fluidsolid interphase momentum exchange, 
:  Tank diameter, m 
:  Velocity of fluid phase, m/s 
:  Interphase velocity, m/s 
:  Velocity of th phase, m/s 
:  Velocity of solid phase, m/s 
:  Phasic average velocity of fluid phase, m/s 
:  Phasic average velocity of solid phase, m/s 
:  Drift velocity, m/s. 
:  Volume fraction of phase 
:  Stressstrain tensor for phase , 
:  Kinematic viscosity of liquid phase, kg/ms 
:  Dispersion Prandtl number 
:  Density of phase , . 
:  Fluid phase 
:  Exchange between fluid and solid phases 
:  Phase number 
:  Solid phase 
:  Turbulent. 
Acknowledgments
The authors gratefully acknowledge the invaluable help provided by Prem Andrade and Mohan Srinivasa (ANSYS, Pune) for their guidance during the CFD simulations.
References
 T. N. Zwietering, “Suspending of solid particles in liquid by agitators,” Chemical Engineering Science, vol. 8, no. 34, pp. 244–253, 1958. View at: Google Scholar
 K. J. Bittorf and S. M. Kresta, “Prediction of cloud height for solid suspensions in stirred tanks,” in Proceedings of the CHISA Conference, Prague, Czech Republic, August 2002. View at: Google Scholar
 A. Ochieng and M. S. Onyango, “CFD simulation of solids suspension in stirred tanks: review,” Hemijska Industrija, vol. 64, no. 5, pp. 365–374, 2010. View at: Publisher Site  Google Scholar
 S. Decker and M. Sommerfeld, “Calculation of particle suspension in agitated vessels with the EulerLagrange approach,” in Proceedings of the IChemE Symposium Series, no. 140, pp. 71–82, July 1996. View at: Google Scholar
 J. J. Derksen, “Numerical simulation of solids suspension in a stirred tank,” AIChE Journal, vol. 49, no. 11, pp. 2700–2714, 2003. View at: Publisher Site  Google Scholar
 H. Barrue, C. Xureb, and J. Bertrand, “A computational study on solid suspension in stirred reactor,” in Proceedings of the ECCE 1, p. 1843, Florence, Italy, 1997. View at: Google Scholar
 G. Montante and F. Magelli, “Modelling of solids distribution in stirred tanks: analysis of simulation strategies and comparison with experimental data,” International Journal of Computational Fluid Dynamics, vol. 19, no. 3, pp. 253–262, 2005. View at: Publisher Site  Google Scholar
 G. Micale, F. Grisafi, L. Rizzuti, and A. Brucato, “CFD simulation of particle suspension height in stirred vessels,” Chemical Engineering Research and Design, vol. 82, no. 9, pp. 1204–1213, 2004. View at: Publisher Site  Google Scholar
 A. R. Khopkar, G. R. Kasat, A. B. Pandit, and V. V. Ranade, “Computational fluid dynamics simulation of the solid suspension in a stirred slurry reactor,” Industrial and Engineering Chemistry Research, vol. 45, no. 12, pp. 4416–4428, 2006. View at: Publisher Site  Google Scholar
 M. Spidla, M. Mostek, V. Sinevic, M. Jahoda, and V. Machon, “Experimental assessment and CFD simulations of local solid concentration profiles in a pilotscale stirred tank,” Chemical Papers, vol. 59, no. 6, pp. 386–393, 2005. View at: Google Scholar
 G. Micale, G. Montante, F. Grisafi, A. Brucato, and J. Godfrey, “CFD simulation of particle distribution in stirred vessels,” Transactions of IChemE, Part A, Chemical Engineering Research and Design, vol. 78, no. 3, pp. 435–444, 2000. View at: Publisher Site  Google Scholar
 G. R. Kasat, A. R. Khopkar, A. B. Pandit, and V. V. Ranade, “CFD simulation of liquidphase mixing in solidliquid stirred reactor,” Chemical Engineering Science, vol. 63, no. 15, pp. 3877–3885, 2008. View at: Publisher Site  Google Scholar
 M. V. Sardeshpande, V. A. Juvekar, and V. V. Ranade, “Hysteresis in cloud heights during solid suspension in stirred tank reactor: experiments and CFD simulations,” AIChE Journal, vol. 56, no. 11, pp. 2795–2804, 2010. View at: Publisher Site  Google Scholar
 M. Ljungqvist and A. Rasmuson, “Numerical simulation of the twophase flow in an axially stirred vessel,” Chemical Engineering Research and Design, vol. 79, no. 5, pp. 533–546, 2001. View at: Publisher Site  Google Scholar
 G. Montante, D. Rondini, A. Bakker, and F. Magelli, “CFD predictions of solid concentration distributions in a baffled stirred vessel agitated with multiple PBT impellers,” in Proceedings of the CHISA Conference, Prague, Czech Republic, August 2002. View at: Google Scholar
 A. Ochieng and A. E. Lewis, “CFD simulation of solids offbottom suspension and cloud height,” Hydrometallurgy, vol. 82, no. 12, pp. 1–12, 2006. View at: Publisher Site  Google Scholar
 R. Angst, E. Harnack, M. Singh, and M. Kraume, “Grid and model dependency of the solid/liquid twophase flow CFD simulation of stirred reactors,” in Proceedings of the 11th European Conference of Mixing, p. 347, Bambarg, Germany, 2003. View at: Google Scholar
 A. Ochieng and A. E. Lewis, “Nickel solids concentration distribution in a stirred tank,” Minerals Engineering, vol. 19, no. 2, pp. 180–189, 2006. View at: Publisher Site  Google Scholar
 D. F. Fletcher and G. J. Brown, “Numerical simulation of solid suspension via mechanical agitation: effect of the modelling approach, turbulence model and hindered settling drag law,” International Journal of Computational Fluid Dynamics, vol. 23, no. 2, pp. 173–187, 2009. View at: Publisher Site  Google Scholar
 H. Barrue, J. Bertrand, B. Cristol, and C. Xuereb, “Eulerian simulation of dense solidliquid suspension in multistage stirred vessel,” Journal of Chemical Engineering of Japan, vol. 34, no. 5, pp. 585–594, 2001. View at: Publisher Site  Google Scholar
 ANSYS FLUENT 13.0 Documentation, ANSYS, Inc., Canonsburg, Pa, USA.
 J. Ding and D. Gidaspow, “Bubbling fluidization model using kinetic theory of granular flow,” AIChE Journal, vol. 36, no. 4, pp. 523–538, 1990. View at: Google Scholar
 E. Deutsch and O. Simonin, “Large eddy simulation applied to the motion of particles in stationary homogeneous fluid turbulence,” in Turbulence Modifications in Multiphase Flows, vol. 110, pp. 35–42, ASME Fluids Engineering Division, 1991. View at: Google Scholar
 O. Simonin, E. Deutsch, and J. P. Minier, “Eulerian prediction of the fluid/particle correlated motion in turbulent twophase flows,” Advances in turbulence 4, Applied Scientific Research, vol. 51, no. 12, pp. 275–283, 1993. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 Shitanshu Gohel et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.