- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Chemical Engineering
Volume 2014 (2014), Article ID 706149, 7 pages
CFD Simulation of the Discharge Flow from Standard Rushton Impeller
1Institute of Hydrodynamics AS CR, v. v. i., Pod Patankou 30/5, 166 12 Prague, Czech Republic
2Department of Process Engineering, Faculty of Mechanical Engineering, Czech Technical University in Prague, Technicka 4, 166 07 Prague, Czech Republic
Received 5 August 2013; Revised 22 November 2013; Accepted 23 November 2013; Published 29 January 2014
Academic Editor: Alírio Rodrigues
Copyright © 2014 Bohuš Kysela 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.
The radial discharge jet from the standard Rushton turbine was investigated by the CFD calculations and compared with results from the Laser Doppler Anemometry (LDA) measurements. The Large Eddy Simulation (LES) approach was employed with Sliding Mesh (SM) model of the impeller motion. The obtained velocity profiles of the mean ensemble-averaged velocity and r.m.s. values of the fluctuating velocity were compared in several distances from the impeller blades. The calculated values of mean ensemble-averaged velocities are rather in good agreement with the measured ones as well as the derived power number from calculations. However, the values of fluctuating velocities are obviously lower from LES calculations than from LDA measurements.
The flow inside the agitated vessel has a key role in the mixing processes. Only the CFD modeling gives us the complex information about the whole flow field in contrast with the results of the experimental measurements. The enormous progress of the computational equipment has allowed using more exacting turbulence models for solution of the flow in the agitated vessel. Nowadays not only Reynolds Averaged Navier-Stokes (RANS) models (for example, [1–7], , and Reynolds Stress Model [4, 5]) are commonly used, but also other more sophisticated methods become topical (for example, Detached Eddy Simulation , Large Eddy Simulation [8–15] and even Direct Numerical Simulation ). However, all calculations also need some validation by the experimental results or by the analytical models. The radial impellers are most often used in experiments and calculations [16–21], namely, Rushton turbine, and there are also analytical models, where the impeller discharge stream is modeled as a turbulent jet [1, 2, 16, 19, 20].
The aim of this study is the description of the turbulent velocity field in the discharge stream from the standard Rushton turbine impeller in the pilot plant mixing vessel with baffles at the wall. Investigation will be carried out experimentally (LDA technique) as well as by means of CFD simulation, where the LES approach was used with Sliding Mesh model for the impeller movement because Sliding Mesh approach has a potential to be used as a design tool to screen different configurations, but it has not been sufficiently validated for turbulent regime .
2. CFD Calculations
A commercial ANSYS FLUENT v.13.0 solver of the finite volume method was employed. The turbulence was modelled by Large Eddy Simulation (LES) with Sliding Mesh (SM) simulation for the impeller movement. The solver was pressure based and for pressure-velocity coupling the PISO method was used. The Smagorinsky-Lilly model was used as the subgrid-scale model with second-order implicit scheme. The schemes for spatial discretization were Gradient—Least Squares Cell based, Pressure—Second Order, and Momentum—Second Order Upwind. The boundary conditions were set: water level to the symmetry and others to the no slip wall, where the part of the impeller shaft outside of the sliding region is defined as wall with impeller speed velocity. The walls of the vessel and baffles are provided by the boundary layer mesh; see Figure 1. The sliding region has cylindrical shape with distance /10 from the impeller cylindrical envelope; see Figure 2. The solved hexahedral meshes consists of 2 465 228 cells (LES 1) and 7 435 557 cells (LES 2), respectively. The larger mesh was refined, namely, in bulk vessel region to attained the maximal cell size under 2 mm, which corresponds with maximal size of the measurement LDA volume; see Section 3. The time step must not exceed 1/60 of one revolution  that corresponds with 0.0032 s for 300 rpm. Hence, the time step 0.001 s was used. The calculated time was 60 s while the flow development required min 20 revolutions [8, 13] which represent the first 2 s of the simulation. Calculated instantaneous flow field over 60 s in the measured plane is depicted in Figure 3.
Measurements of the velocity profiles were carried out in a pilot plant flat bottomed mixing vessel with four baffles at its wall (see Figure 4), with water as the working liquid (density kg m−3, dynamic viscosity mPa s) under the constant impeller speed 300 rpm (impeller Reynolds number ). A standard Rushton impeller [22, 23] was used for the investigation (see Figure 5).
Laser Doppler Anemometry (LDA) one-component measurements of the radial velocity were performed in the impeller discharge stream (see Figure 4) in the vertical plane between two adjacent baffles in axial profiles with several distances from the impeller blade. The dimensionless radial coordinates were 1.2; 1.4; 1.6; 1.8; 2.0; 2.2. One-component LDA system setup consists of Coherent INNOVA 305 Ion-Argon laser supply with power 5 W and separated beam only for one component measurement on wavelength 514.5 nm; DANTEC FiberFlow transmitting optics and P80 DANTEC BSA processor. The optic parameters were focal length 310 mm, diameter of the beam 1.2 mm, fringe spacing 2.67 μm, number of fringes 63, and the measured volume was ellipsoid with sizes of the axes mm. The used frequency shift was 40 MHz and velocity span 7.51 m/s. The setup was supervised by BSA FLOW SOFTWARE v3.0 installed on standard PC where the data was processed. S-HGS (silver coated-hollow glass Spheres) with mean diameter 10 μm and density 1.1 gcm−3 were used as trace particles. The measurement was performed through the glass flat bottom of the vessel to eliminate optical effects of the cylindrical walls.
4. Results and Discussion
Profiles of the radial mean ensemble-averaged velocity component in the dimensionless form where the radial velocity component is normalized by the impeller tip speed are depicted in Figures 6, 7, 8, 9, 10, and 11. The dimensionless coordinate is the distance from the impeller disk axis normalized by the half-height of the impeller blade. The profiles are depicted for six values of the dimensionless radius .
The depicted results from CFD calculations and results in measured points were compared also by calculations of variance. Variance between the measured and the calculated data of dimensionless radial component of the mean ensemble-averaged velocity was calculated using the formula where is the number of compared points on the profile with index LDA means the value obtained from LDA measurements and with index CFD is the value interpolated from calculated profile from CFD. Standard deviation is expressed as the square root of the variance. The results of variance (standard deviation) in Table 1 signify that the higher discrepancy of the profiles is in the region where the zone of establishment is changing to the zone of the established flow [1, 2]. The increase of the standard deviation with increasing dimensionless radius is probably caused by the different shape of the discharge stream and it seems to be different also at the vertical position of the stream which depends on the impeller off-bottom clearance.
The r.m.s. values of the fluctuation velocity were treated as well as the mean velocities to the dimensionless form. The two values of the dimensionless radius and 1.6 are shown in Figure 12, where the trailing vortices have an impact on the fluctuation velocity. The zone is titled the zone of establishment (ZFE) . The calculated values of the fluctuations are rather lower than the measured ones. The results in the next zone titled the zone of established flow (ZEF) are depicted in Figure 13. There are compared values of the dimensionless radius and 2.2. It seems that agreement between computed and experimentally determined values of the fluctuating velocity is better in ZEF than in ZFE.
The power number was calculated from impeller torque which was obtained from the force balance on the impeller surface provided by the CFD calculations :
For Rushton turbine, the results of the power number could be compared with the empirical correlation [22, 23], relating the power number to the relative blade thickness and the relative vessel diameter : where mm is the thickness of the separating disc of a standard Rushton impeller and quantity m. The power number derived from calculations (2) was and it is in a good agreement with power number calculated by (3), where .
The flow in the discharge stream from the standard Rushton turbine was calculated by the Large Eddy Simulation approach. The comparison of the mean radial ensemble-averaged velocity profiles obtained from LDA measurements are in good agreement with the calculated results from both LES cases. The r.m.s. values of fluctuating velocity show the similar shape of profiles, but the calculations mostly underestimated the values obtained by the LDA measurements. The ensemble-averaging shows the dependency on the spatial resolution of the calculations or the measurement method; namely, the r.m.s. values of fluctuations are strongly affected by a spatial averaging. Calculated values of standard deviation from the mean velocities are increasing from lower dimensionless radius to the value where it probably indicates the boundary between the zone of the flow establishment and the zone of the established flow. The same trend is shown in the comparisons of the r.m.s. values of the fluctuating velocities. The power number derived from the impeller torque calculations of the presented LES numerical modeling is very close to the value estimated from empirical correlation based on experimental measurements.
|Width of radial baffle (m)|
|:||Impeller off-bottom clearance (m)|
|:||Impeller diameter (m)|
|:||Height of agitated charge above vessel bottom (m)|
|:||Impeller blade width (m)|
|:||Impeller torque (Nm)|
|:||Number of experimental points|
|:||Impeller speed (s−1)|
|:||Impeller power number|
|:||Dimensional longitudinal (radial) coordinate (m)|
|:||Dimensionless longitudinal (radial) coordinate|
|:||Impeller Reynolds number|
|:||Thickness of the impeller separating disc (m)|
|:||Mixing vessel diameter (m)|
|:||Impeller tip speed (ms−1)|
|:||Dimensionless radial ensemble-averaged mean velocity|
|:||Dimensionless root mean square value of the radial component of fluctuating velocity.|
|:||Dynamic viscosity (Pa s)|
|:||Density (kg m−3)|
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
This research has been subsidized by the Research Project no. GA CR P101/12/2274 and RVO: 67985874.
- R. Ben-Nun and M. Sheintuch, “Characterizing turbulent jet properties of radial discharge impeler: the zone of flow establishment (ZFE) and ZEF,” AIChE Journal, 2013.
- R. Ben-Nun and M. Sheintuch, “Charcterizing turbulent jet properties of radial discharge impeller: potential core, spreading rate and averaged flow field parameters,” in Proceedings of the 9th European Congress of Chemical Engineering, The Hague, The Netherlands, April 2013.
- S. Yeoh, G. Papadakis, and M. Yianneskis, “Numerical simulation of turbulent flow characteristics in a stirred vessel using the LES and RANS approaches with the sliding/deforming mesh methodology,” Chemical Engineering Research and Design, vol. 82, no. 7, pp. 834–848, 2004.
- J. B. Joshi, N. K. Nere, C. V. Rane et al., “CFD simulation of stirred tanks: comparison of turbulence models. Part I: radial flow impellers,” Canadian Journal of Chemical Engineering, vol. 89, no. 1, pp. 23–82, 2011.
- J. B. Joshi, N. K. Nere, C. V. Rane et al., “CFD simulation of stirred tanks: comparison of turbulence models (Part II: axial flow impellers, multiple impellers and multiphase dispersions),” Canadian Journal of Chemical Engineering, vol. 89, no. 4, pp. 754–816, 2011.
- M. Coroneo, G. Montante, A. Paglianti, and F. Magelli, “CFD prediction of fluid flow and mixing in stirred tanks: numerical issues about the RANS simulations,” Computers and Chemical Engineering, vol. 35, no. 10, pp. 1959–1968, 2011.
- A. Barker, R. D. Laroche, M. H. Wang, and R. V. Calabrese, “Sliding mesh simulation of laminar flow in stirred reactors,” Chemical Engineering Research and Design, vol. 75, no. 1, pp. 42–44, 1997.
- J. Gimbun, C. D. Rielly, Z. K. Nagy, and J. J. Derksen, “Detached eddy simulation on the turbulent flow in a stirred tank,” AIChE Journal, vol. 58, no. 10, pp. 3224–3241, 2012.
- J. Derksen and H. E. A. Van Den Akker, “Large eddy simulations on the flow driven by a Rushton turbine,” AIChE Journal, vol. 45, no. 2, pp. 209–221, 1999.
- J. Derksen, “Long-time solids suspension simulations by means of a large-eddy approach,” Chemical Engineering Research and Design, vol. 84, no. 1, pp. 38–46, 2006.
- A. Bakker and L. M. Oshinowo, “Modelling of turbulence in stirred vessels using large eddy simulation,” Chemical Engineering Research and Design, vol. 82, no. 9, pp. 1169–1178, 2004.
- M. Jahoda, M. Moštěk, A. Kukuková, and V. Machoň, “CFD modelling of liquid homogenization in stirred tanks with one and two impellers using large eddy simulation,” Chemical Engineering Research and Design, vol. 85, no. 5, pp. 616–625, 2007.
- Z. Li, M. Hu, Y. Bao, and Z. Gao, “Particle image velocimetry experiments and large eddy simulations of merging flow characteristics in dual Rushton turbine stirred tanks,” Industrial and Engineering Chemistry Research, vol. 51, no. 5, pp. 2438–2450, 2012.
- Z. Li, Y. Bao, and Z. Gao, “PIV experiments and large eddy simulations of single-loop flow fields in Rushton turbine stirred tanks,” Chemical Engineering Science, vol. 66, no. 6, pp. 1219–1231, 2011.
- J. J. Gillissen and H. E. Van den Akker, “Direct numerical simulation of the turbulent flow in a baffled tank driven by a Rushton turbine,” AIChE Journal, vol. 58, no. 12, pp. 3878–3890, 2012.
- V. Kolář, P. Filip, and A. Curev, “Hydrodynamics of radially discharging impeller stream in agitated vessels,” Chemical Engineering Communications, vol. 27, no. 5-6, pp. 313–326, 1984.
- J. Drbohlav, I. Fořt, K. Máca, and J. Ptáček, “Turbulent characteristics of discharge flow from turbine impeller,” Collection of Czechoslovak Chemical Communications, vol. 43, no. 12, pp. 3148–3161, 1978.
- I. Fořt, H. O. Möckel, J. Drbohlav, and M. Hrach, “The flow of a liquid in a stream from the standard turbine impeller,” Collection of Czechoslovak Chemical Communications, vol. 44, no. 3, pp. 700–710, 1979.
- A. Obeid, I. Fořt, and J. Bertrand, “Hydrodynamic characteristics of flow in systems with turbine impeller,” Collection of Czechoslovak Chemical Communications, vol. 48, no. 2, pp. 568–577, 1983.
- J. Talaga and I. Fořt, “The velocity field in the discharge stream from a rushton turbine impeller,” in Proceedings of the 14th European Conference on Mixing, Warszava, Poland, September 2012.
- B. C. Venneker, J. J. Derksen, and H. E. A. Van den Akker, “Turbulent flow of shear-thinning liquids in stirred tanks-The effects of Reynolds number and flow index,” Chemical Engineering Research and Design, vol. 88, no. 7, pp. 827–843, 2010.
- W. Bujalski, A. W. Nienow, S. Chatwin, and M. Cooke, “The dependency on scale of power numbers of Rushton disc turbines,” Chemical Engineering Science, vol. 42, no. 2, pp. 317–326, 1987.
- K. R. Beshay, J. Kratěna, I. Fořt, and O. Brůha, “Power input of high-speed rotary impellers,” Acta Polytechnica, vol. 41, no. 6, pp. 18–23, 2001.