Application of Wray-Agarwal turbulence model for numerical simulation of gas-solid flows in CFB risers

2021 ◽  
pp. 1-25
Author(s):  
Yali Shao ◽  
Ramesh K. Agarwal ◽  
Xudong Wang ◽  
Baosheng Jin
Keyword(s):  
Turbulence Model ◽  
Circulating Fluidized Bed ◽  
Volume Fraction ◽  
Solid Phase ◽  
Solid Volume Fraction ◽  
Turbulence Models ◽  
Two Phase ◽  
Efficient Simulation ◽  
Pressure Profiles ◽  
First Time

Abstract In recent decades, increasing attention has been focused on accurate modeling of circulating fluidized bed (CFB) risers to provide valuable guidance to design, optimization and operation of reactors. Turbulence model plays an important role in accurate prediction of complex gas-solid flows. Recently developed Wray-Agarwal (WA) model is a one-equation turbulence model with the advantages of high computational efficiency and competitive accuracy with two-equation models. In this paper for the first time, Eulerian-Eulerian approach coupled with different turbulence models including WA model, standard κ-ε model and shear stress transport (SST) κ-ω model is employed to simulate two-phase flows of gas phase and solid phase in two CFB risers, in order to assess accuracy and efficiency of WA model compared to other well-known two-equation models. Predicted gas-solid flow dynamic characteristics including the gas-solid volume fraction distributions in radial and axial directions, pressure profiles and solid mass flux distributions are compared with data obtained from experiment in detail. The results demonstrate WA model is very promising for accurate and efficient simulation of gas-solid multiphase flows.

Simulation Studies of Gas-Solid in the Riser of a Circulating Fluidized Bed

2005 ◽  
Author(s):  
Ahmad Hussain ◽  
Farid Nasir Ani ◽  
Amer Nordin Darus ◽  
Azeman Mustafa ◽  
Arshad A. Salema
Keyword(s):  
Gas Phase ◽  
Fluidized Bed ◽  
Turbulence Model ◽  
Circulating Fluidized Bed ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Design Parameters ◽  
Cfd Model ◽  
End Effects ◽  
Wide Range

A numerical parametric study was performed on the influence of various riser exit geometries on the hydrodynamics of gas-solid two-phase flow in the riser of a Circulating Fluidized Bed (CFB). A Eulerian continuum formulation was applied to both phases. A two fluid framework has been used to simulate fully developed gas-solid flows in vertical riser. A two dimensional Computational Fluid Dynamics (CFD) model of gas-particle flow in the CFB has been investigated using the code FLUENT. The turbulence was modeled by a k-ε turbulence model in the gas phase. The simulations were done using the geometrical configuration of a CFB test rig at the Universiti Teknologi Malaysia (UTM). The CFB riser column has 265 mm (width), 72 mm (depth) and 2.7 m height. The riser is made up of interchangeable Plexiglas columns. The computational model was used to simulate the riser over a wide range of operating and design parameters. In addition, several numerical experiments were carried out to understand the influence of riser end effects, particle size, gas solid velocity and solid volume fraction on the simulated flow characteristics. The CFD model with a k-ε turbulence model for the gas phase and a fixed particle viscosity in the solids phase showed good mixing behaviour. These results were found to be useful in further development of modeling of gas solid flow in the riser.

Comparative Assessment of Turbulence Models for Prediction of Flow-Induced Corrosion Damages

2011 ◽  
Author(s):  
Kaushik Das ◽  
Debashis Basu ◽  
Todd Mintz
Keyword(s):  
Corrosion Rate ◽  
Cross Sections ◽  
Volume Fraction ◽  
Solid Phase ◽  
Turbulence Models ◽  
Reynolds Stress Model ◽  
Continuous Phase ◽  
Comparative Assessment ◽  
Solid Particles ◽  
Protective Film

The present study makes a comparative assessment of different turbulence models in simulating the flow-assisted corrosion (FAC) process for pipes with noncircular cross sections and bends, features regularly encountered in heat exchangers and other pipeline networks. The case study investigates material damage due to corrosion caused by dissolved oxygen (O2) in a stainless steel pipe carrying an aqueous solution. A discrete solid phase is also present in the solution, but the transport of the solid particles is not explicitly modeled. It is assumed that the volume fraction of the solid phase is low, so it does not affect the continuous phase. Traditional two-equation models are compared, such as isotropic eddy viscosity, standard k-ε and k-ω models, shear stress transport (SST) k-ω models, and the anisotropic Reynolds Stress Model (RSM). Computed axial and radial velocities, and turbulent kinetic energy profiles predicted by the turbulence models are compared with available experimental data. Results show that all the turbulence models provide comparable results, though the RSM model provided better predictions in certain locations. The convective and diffusive motion of dissolved O2 is calculated by solving the species transport equations. The study assumes that solid particle impingement on the pipe wall will completely remove the protective film formed by corrosion products. It is also assumed that the rate of corrosion is controlled by diffusion of O2 through the mass transfer boundary layer. Based on these assumptions, corrosion rate is calculated at the internal pipe walls. Results indicate that the predicted O2 corrosion rate along the walls varies for different turbulence models but show the same general trend and pattern.

Numerical Simulation of Modified Trailing Edge of Francis Turbine Stay Vane Based on Solid-Liquid Two-Phase Flow

2014 ◽  
Vol 700 ◽  
pp. 643-646
Author(s):  
Dong Wang ◽  
Si Qing Zhang ◽  
Yun Long Zhang
Keyword(s):  
Numerical Simulation ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Francis Turbine ◽  
Phase Flow ◽  
Oblique Angle ◽  
Trailing Edge ◽  
Two Phase ◽  
High Solid ◽  
Solid Liquid

In order to investigate the silt abrasion of modified trailing edge of stay vane in Francis turbine, the numerical simulation of trailing edge with different geometries were carried out based on the solid-liquid two-phase flow by means of Computation Fluid Dynamics. The results show that low solid volume fraction distributes on the chamfered surface of trailing edge, and high solid volume fraction distributes on the end of oblique surface. The smaller the modified angle is, the larger the distribution area of high solid volume fraction is, which show the trailing edge with smaller oblique angle may suffer from silt abrasion. Therefore, in order to solve the vibration caused by Karman vortex the trailing edge has to be sharpened, the oblique angle of trailing edge should not be too small. At end of trailing edge needs to ensure a certain thickness, especially the trailing edge near the lower ring can be thicker, which can meet the anti-abrasion requirements.

scholarly journals Computational investigation of liquid-solid slurry flow through an expansion in a rectangular duct

2014 ◽  
Vol 62 (3) ◽  
pp. 234-240 ◽  
Author(s):  
Gianandrea Vittorio Messa ◽  
Stefano Malavasi
Keyword(s):  
Volume Fraction ◽  
Particle Density ◽  
Fluid Model ◽  
Solid Volume Fraction ◽  
Solid Particles ◽  
Rectangular Duct ◽  
Two Phase ◽  
Horizontal Pipe ◽  
Two Fluid Model ◽  
Two Fluid

Abstract The flow of a mixture of liquid and solid particles at medium and high volume fraction through an expansion in a rectangular duct is considered. In order to improve the modelling of the phenomenon with respect to a previous investigation (Messa and Malavasi, 2013), use is made of a two-fluid model specifically derived for dense flows that we developed and implemented in the PHOENICS code via user-defined subroutines. Due to the lack of experimental data, the two-fluid model was validated in the horizontal pipe case, reporting good agreement with measurements from different authors for fully-suspended flows. A 3D system is simulated in order to account for the effect of side walls. A wider range of the parameters characterizing the mixture (particle size, particle density, and delivered solid volume fraction) is considered. A parametric analysis is performed to investigate the role played by the key physical mechanisms on the development of the two-phase flow for different compositions of the mixture. The main focuses are the distribution of the particles in the system and the pressure recovery

Computational Study of Slurry Flow in Pipes to Determine Profiles Sensed in Near Infrared Slurry Measurements

2011 ◽  
Author(s):  
V. Pasangulapati ◽  
N. R. Kesana ◽  
G. Sharma ◽  
F. W. Chambers ◽  
M. E. McNally ◽  
...  
Keyword(s):  
Near Infrared ◽  
Volume Fraction ◽  
Solid Phase ◽  
Computational Study ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Optical Measurements ◽  
Concentration Profiles ◽  
Horizontal Pipe ◽  
Volume Fractions

It is desired to perform accurate Near Infrared sensor measurements of slurries flowing in pipes leaving large batch reactors. A concern with these measurements is the degree to which the slurry sensed is representative of the material in the reactor and flowing through the pipe. Computational Fluid Dynamics (CFD) has been applied to the flow in the pipe to determine the flow fields and the concentration profiles seen by the sensors. The slurry was comprised of a xylene liquid phase and an ADP (2-amino-4, 6-dimethylpyrimidine) solid phase with a density ratio of 1.7. Computations were performed for a horizontal pipe with diameter 50.8 mm, length 2.032 m, and 1.76 m/s and 3.26 m/s mixture velocities. The corresponding pipe Reynolds numbers were 1.19E+05 and 2.21E+05. The flow through a slotted cylindrical probe inserted radially in the pipe also was considered. Spherical slurry particles with diameters from 10 μm to 1000 μm were considered with solid volume fractions of 12%, 24%, and 35%. Computations were performed with ANSYS FLUENT 12 software using the Realizable k-ε turbulence model and the enhanced wall treatment function. Comparisons of computed vertical profiles of solid volume fraction to results in the literature showed good agreement. Symmetric, nearly flat solid volume fraction profiles were observed for 38 μm particles for all three initial solid volume fractions. Asymmetric solid volume fraction profiles with greater values toward the bottom were observed for the larger particles. Changes in the profiles of turbulent kinetic energy also were observed. These changes are important for optical measurements which depend upon the mean concentration profiles as well as the turbulent motion of the slurry particles.

Comparison of Two Two-Equation Turbulence Model Used for the Numerical Simulation of Underwater Ammunition Fuze Turbine Flow Field

2012 ◽  
Vol 591-593 ◽  
pp. 1968-1972
Author(s):  
De Zhang Shen ◽  
He Zhang ◽  
Hao Jie Li
Keyword(s):  
Numerical Simulation ◽  
Turbulence Model ◽  
Turbulence Models ◽  
Axial Pressure ◽  
Two Phase ◽  
Speed Increase ◽  
Pressure Coefficients ◽  
Surface Pressure Distribution ◽  
Calculation Results ◽  
Distribution Trends

To figure out the problem of turbulence simulation of underwater ammunition fuze turbine numerical simulation, respectively, realizable k-ε turbulence model and SST k-ω turbulence model are used for two-phase flow numerical simulation of the turbine rotation. The analysis compared the calculation results of the two turbulence models. The results showed that: the cavitation scale obtained from realizable k-ε turbulence model is shorter than that of SST k-ω turbulence model; turbine surface pressure distribution trends are similar of this two model, the results of realizable k-ε turbulence model are bigger than SST k-ω turbulence model; the turbine axial pressure coefficients using realizable k-ε turbulence model are also bigger than that of SST k-ω turbulence model, and the deviation increases with the speed increase.

Computational study of granular shear flows of dry flexible fibres using the discrete element method

2015 ◽  
Vol 775 ◽  
pp. 24-52 ◽  
Author(s):  
Y. Guo ◽  
C. Wassgren ◽  
B. Hancock ◽  
W. Ketterhagen ◽  
J. Curtis
Keyword(s):  
Discrete Element Method ◽  
Shear Flows ◽  
Volume Fraction ◽  
Solid Phase ◽  
Computational Study ◽  
Solid Volume Fraction ◽  
Discrete Element ◽  
Effective Coefficients ◽  
Phase Stresses ◽  
Element Method

In this study, shear flows of dry flexible fibres are numerically modelled using the discrete element method (DEM), and the effects of fibre properties on the flow behaviour and solid-phase stresses are explored. In the DEM simulations, a fibre is formed by connecting a number of spheres in a straight line using deformable and elastic bonds. The forces and moments induced by the bond deformation resist the relative normal, tangential, bending and torsional movements between two bonded spheres. The bond or deforming stiffness determines the flexibility of the fibres and the bond damping accounts for the energy dissipation in the fibre vibration. The simulation results show that elastically bonded fibres have smaller effective coefficients of restitution than rigidly connected fibres. Thus, smaller solid-phase stresses are obtained for flexible fibres, particularly with bond damping, compared with rigid fibres. Frictionless fibres tend to align with a small angle from the flow direction as the solid volume fraction increases, and fibre deformation is minimized due to the alignment. However, jamming, with a corresponding sharp stress increase, large fibre deformation and dense contact force network, occurs for fibres with friction at high solid volume fractions. It is also found that jamming is more prevalent in dense flows with larger fibre friction coefficient, rougher surface, larger stiffness and larger aspect ratio.

scholarly journals Buoyancy-induced convection from a pair of heated and cooled horizontal circular cylinders inside an adiabatic tilted cavity filled with alumina/water nanofluids

2019 ◽  
Vol 30 (6) ◽  
pp. 3163-3181
Author(s):  
Massimo Corcione ◽  
Emanuele Habib ◽  
Alessandro Quintino ◽  
Elisa Ricci ◽  
Vincenzo Andrea Spena
Keyword(s):  
Heat Transfer ◽  
Heat Transfer Performance ◽  
Volume Fraction ◽  
Solid Phase ◽  
Suspended Solid ◽  
Circular Cylinders ◽  
Transfer Performance ◽  
Two Phase ◽  
Content Type ◽  
Tilting Angle

Purpose This paper aims to investigate numerically buoyancy-induced convection from a pair of differentially heated horizontal circular cylinders set side by side in a nanofluid-filled adiabatic square enclosure, inclined with respect to gravity so that the heated cylinder is located below the cooled one, using a two-phase model based on the double-diffusive approach assuming that the Brownian diffusion and thermophoresis are the only slip mechanisms by which the solid phase can develop a significant relative velocity with respect to the liquid phase. Design/methodology/approach The system of the governing equations of continuity, momentum and energy for the nanofluid, and continuity for the nanoparticles, is solved by a computational code based on the SIMPLE-C algorithm. Numerical simulations are performed for Al2O3 + H2O nanofluids using the average volume fraction of the suspended solid phase, the tilting angle of the enclosure, the nanoparticle size, the average nanofluid temperature and the inter-cylinder spacing, as independent variables. Findings The main results obtained may be summarized as follows: at high temperatures, the nanofluid heat transfer performance relative to that of the pure base liquid increases with increasing the average volume fraction of the suspended solid phase, whereas at low temperatures it has a peak at an optimal particle loading; the relative heat transfer performance of the nanofluid has a peak at an optimal tilting angle of the enclosure; the relative heat transfer performance of the nanofluid increases notably as the average temperature is increased, and just moderately as inter-cylinder spacing is increased and the nanoparticle size is decreased. Originality/value The two-phase computational code used in the present study incorporates three empirical correlations for the evaluation of the effective thermal conductivity, the effective dynamic viscosity and the coefficient of thermophoretic diffusion, all based on a high number of literature experimental data.

Three-Dimensional Numerical Simulation of Convective Melting of Solid Particles in a Fluid

1999 ◽  
Author(s):  
Y. L. Hao ◽  
Y.-X. Tao
Keyword(s):  
Volume Fraction ◽  
Solid Phase ◽  
Rectangular Cross Section ◽  
Three Dimensional ◽  
Phase Changes ◽  
Solid Particles ◽  
Preliminary Calculation ◽  
Two Phase ◽  
Closed Set ◽  
Gravity Effects

Abstract A physical model of two-phase flow and heat-mass transfer with the phase changes based on the theory of interacting continua is proposed. All terms in the conservation equations are analyzed and the constitutive equations are presented. A closed set of governing equations describing the convective melting of solid particles in a fluid is obtained. The numerical method is developed for the solution of velocity, temperature, and volume fraction of solid phase for the three-dimensional melting in a rectangular cross-section channel. Preliminary calculation, including gravity effects, shows that the result is reasonable. This study provides a basis for the theoretical and experimental investigation of convective melting of solid particles in a fluid.

Sign in / Sign up

Export Citation Format

Share Document