PIV measurement of particle motion in spiral gas–solid two-phase flow

With a concise review on some basic and novel algorithms and methods for the techniques of particle-imaging velocimetry (PIV), the paper reports an application of the PIV techniques to the investigation of particle motion in a gas–solid two-phase spiral flow in a horizontal tube. Axial velocities of the transported particles are obtained. Some important features of particle motion governing high transportation efficiency of the spiral flow are revealed by investigating probability density distribution of particle locations in a pipe cross-section.     

On the use of Csanady’s formulae in a turbulent gas–solid channel flow

The paper examines the use of expressions proposed by Csanady to predict the influence of the crossing trajectory and continuity effects on the decorrelation time scales of the fluid along solid particle trajectories in horizontal and downward vertical channel flows. The model is evaluated using data provided by a direct numerical simulation (DNS) of the carrier phase combined with a Lagrangian simulation of discrete particle (LS). Two particle relaxation times and two values of the gravity acceleration are considered. The results show the possibility of using Csanady’s expressions in a turbulent channel flow provided that the spatial and temporal correlations anisotropy is included in the model. As in isotropic homogeneous turbulence, a decrease of the decorrelation time scales is found to be more important in the directions perpendicular to the mean relative velocity.     

Simulation of gas–solid particle flows over a wide range of concentration

A two-fluid model of gas–solid particle flows that is valid for a wide range of the solid-phase volume concentration (dilute to dense) is presented. The governing equations of the fluid phase are obtained by volume averaging the Navier–Stokes equations for an incompressible fluid. The solid-phase macroscopic equations are derived using an approach that is based on the kinetic theory of dense gases. This approach accounts for particle–particle collisions. The model is implemented in a control-volume finite element method for simulations of the flows of interest in two-dimensional, planar or axisymmetric, domains. The chosen mathematical model and the proposed numerical method are applied to three test problems and one demonstration problem. © 1998 John Wiley & Sons, Ltd.     

Parallelization of pseudo-particle modeling and its application in simulating gas–solid fluidization

Pseudo-Particle Modeling (PPM) is a particle method proposed by Ge and Li in 1996 [Ge, W., & Li, J. (1996). Pseudo-particle approach to hydrodynamics of particle–fluid systems. In M. Kwauk & J. Li (Eds.), Proceedings of the 5th international conference on circulating fluidized bed (pp. 260–265). Beijing: Science Press] and has been used to explore the microscopic mechanism in complex particle–fluid systems. But as a particle method, high computational cost remains a main obstacle for its large-scale application; therefore, parallel implementation of this method is highly desirable. Parallelization of two-dimensional PPM was carried out by spatial decomposition in this paper. The time costs of the major functions in the program were analyzed and the program was then optimized for higher efficiency by dynamic load balancing and resetting of particle arrays. Finally, simulation on a gas–solid fluidized bed with 102,400 solid particles and 1.8 × 10⁷ pseudo-particles was performed successfully with this code, indicating its scalability in future applications.     

Convective heat transfer in the laminar–turbulent transition region with molten salt in a circular tube

In order to understand the heat transfer characteristics of molten salt and testify the validity of the well-known empirical convective heat transfer correlations, experimental study on transition convective heat transfer with molten salt in a circular tube was conducted. Molten salt circulations were realized and operated in a specially designed system over 1000 h. The average forced convective heat transfer coefficients of molten salt were determined by least-squares method based on the measured data of flow rates and temperatures. Finally, a heat transfer correlation of transition flow with molten salt in a circular tube was obtained and good agreement was observed between the experimental data of molten salt and the well-known correlations presented by Hausen and Gnielinski, respectively.     

High-order LES of turbulent heat transfer in a rotor–stator cavity

The present work examines the turbulent flow in an enclosed rotor–stator system subjected to heat transfer effects. Besides their fundamental importance as three-dimensional prototype flows, such flows arise in many industrial applications but also in many geophysical and astrophysical settings. Large eddy simulations (LES) are here performed using a spectral vanishing viscosity technique. The LES results have already been favorably compared to velocity measurements in the isothermal case (Séverac, E., Poncet, S., Serre, E., Chauve, M.P., 2007. Large eddy simulation and measurements of turbulent enclosed rotor–stator flows. Phys. Fluids, 19, 085113) for a large range of Reynolds numbers 10⁵Re=Ωb²/ν10⁶, in an annular cavity of large aspect ratio G=(b-a)/H=5 and weak curvature parameter R_m=(b-a)/(b+a)=1.8 (a,b the inner and outer radii of the rotor and H the interdisk spacing). The purpose of this paper is to extend these previous results in the non-isothermal case using the Boussinesq approximation to take into account the buoyancy effects. Thus, the effects of thermal convection have been examined for a turbulent flow Re=10⁶ of air in the same rotor–stator system for Rayleigh numbers up to Ra=10⁸. These LES results provide accurate, instantaneous quantities which are of interest in understanding the physics of turbulent flows and heat transfers in an interdisk cavity. Even at high Rayleigh numbers, the structure of the iso-values of the instantaneous normal temperature gradient at the disk surfaces resembles the one of the iso-values of the tangential velocity with large spiral arms along the rotor and more thin structures along the stator. The averaged results show small effects of density variation on the mean and turbulent fields. The turbulent Prandtl number is a decreasing function of the distance to the wall with 1.4 close to the disks and about 0.3 in the outer layers. The local Nusselt number is found to be proportional to the local Reynolds number to the power 0.7. The evolution of the averaged Bolgiano length scale L_B with the Rayleigh number indicates that temperature fluctuations may have a large influence on the dynamics only at the largest scales of the system for Ra10⁷, since L_B remains lower than the thermal boundary layer thicknesses.     

Numerical simulation and experimental validation of heat transfer within rotating flows for three‐dimensional non‐axisymmetric,turbulent conditions

A control volume type numerical methodology for the analysis of steady three‐dimensional rotating flows with heat transfer, in both laminar and turbulent conditions, is implemented and experimentally tested. Non‐axisymmetric momentum and heat transfer phenomena are allowed for. Turbulent transport is alternatively represented through three existing versions of the k–ε model that were adjusted to take into account the turbulence anisotropy promoted by rotation, streamline curvature and thermal buoyancy. Their relative performance is evaluated by comparison of calculated local and global heat balances with those obtained through measurements in a laboratory device. A modified version of the Lam and Bremhorst, low Reynolds number model is seen to give the best results. A preliminary analysis focused on the flow structure and the transfer of heat is reported. Copyright © 2002 John Wiley & Sons, Ltd.     

Two-phase flow distribution of air–water annular flow in a parallel flow heat exchanger

The air and water flow distribution are experimentally studied for a round header – flat tube geometry simulating a parallel flow heat exchanger. The number of branch flat tube is 30. The effects of tube outlet direction, tube protrusion depth as well as mass flux, and quality are investigated. The flow at the header inlet is identified as annular. For the downward flow configuration, the water flow distribution is significantly affected by the tube protrusion depth. For flush-mounted configuration, most of the water flows through frontal part of the header. As the protrusion depth increases, more water is forced to the rear part of the header. The effect of mass flux or quality is qualitatively the same as that of the protrusion depth. Increase of the mass flux or quality forces the water to rear part of the header. For the upward flow configuration, however, most of the water flows through rear part of the header. The protrusion depth, mass flux, or quality does not significantly alter the flow pattern. Possible explanations are provided based on the flow visualization results. Negligible difference on the water flow distribution was observed between the parallel and the reverse flow configuration.     

Numerical simulation and experimental validation of gas–solid flow in the riser of a dense fluidized bed reactor

Gas–solid flow in the riser of a dense fluidized bed using Geldart B particles (sand), at high gas velocity (7.6–15.5 m/s) and with comparatively high solid flux (140–333.8 kg/m² s), was investigated experimentally and simulated by computational fluid dynamics (CFD), both two- and three-dimensional and using the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd and EMMS drag models. The results predicted by EMMS drag model showed the best agreement with experimental results. Calculated axial solids hold-up profiles, in particular, are well consistent with experimental data. The flow structure in the riser was well represented by the CFD results, which also indicated the cause of cluster formation. Complex hydrodynamical behaviors of particle cluster were observed. The relative motion between gas and solid phases and axial heterogeneity in the three subzones of the riser were also investigated, and were found to be consistent with predicted flow structure. The model could well depict the difference between the two exit configurations used, viz., semi-bend smooth exit and T-shaped abrupt exit. The numerical results indicate that the proposed EMMS method gives better agreement with the experimental results as compared with the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas–solid two-phase flow.     

A Lagrangian–Eulerian model of particle dispersion in a turbulent plane mixing layer

A Lagrangian–Eulerian model for the dispersion of solid particles in a two‐dimensional, incompressible, turbulent flow is reported and validated. Prediction of the continuous phase is done by solving an Eulerian model using a control‐volume finite element method (CVFEM). A Lagrangian model is also applied, using a Runge–Kutta method to obtain the particle trajectories. The effect of fluid turbulence upon particle dispersion is taken into consideration through a simple stochastic approach. Validation tests are performed by comparing predictions for both phases in a particle‐laden, plane mixing layer airflow with corresponding measurements formerly reported by other authors. Even though some limitations are detected in the calculation of particle dispersion, on the whole the validation results are rather successful. Copyright © 2002 John Wiley & Sons, Ltd.     

Motion and heat and mass transfer in a disperse admixture in turbulent nonisothermal jets of a gas and a low-temperature plasma

The motion and heat and mass transfer of particles of a disperse admixture in nonisothermal jets of a gas and a low-temperature plasma are simulated with allowance for the migration mechanism of particle motion actuated by the turbophoresis force and the influence of turbulent fluctuations of the jet flow velocity on heat and mass transfer of the particle. The temperature distribution inside the particle at each time step is found by solving the equation of unsteady heat conduction. The laws of scattering of the admixture and the laws of melting and evaporation of an individual particle are studied, depending on the injection velocity and on the method of particle insertion into the jet flow. The calculated results are compared with data obtained with ignored influence of turbulent fluctuations on the motion and heat and mass transfer of the disperse phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 95–108, May–June, 2008.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

Solid‐particle erosion in the tube end of the tube sheet of a shell‐and‐tube heat exchanger

Erosion is one of the major problems in many industrial processes, and in particular, in heat exchangers. The effects of flow velocity and sand particle size on the rate of erosion in a typical shell‐and‐tube heat exchanger were investigated numerically using the Lagrangian particle‐tracking method. Erosion and penetration rates were obtained for sand particles of diameters ranging from 10 to 500 µm and for inlet flow velocities ranging from 0.197 to 2.95 m/s. A flow visualization experiment was conducted with the objective of verifying the accuracy of the continuous phase calculation procedure. Comparison with available experimental data of penetration rates was also conducted. These comparisons resulted in a good agreement. The results show that the location and number of eroded tubes depend mainly on the particle size and velocity magnitude at the header inlet. The rate of erosion depends exponentially on the velocity. The particle size shows negligible effect on the erosion rate at high velocity values and the large‐size particles show less erosion rates compared to the small‐size particles at low values of inlet flow velocities. The results indicated that the erosion and penetration rates are insignificant at the lower end of the velocity range. However, these rates were found to increase continuously with the increase of the inlet flow velocity for all particle sizes. The particle size creating the highest erosion rate was found to depend on the flow velocity range. Copyright © 2005 John Wiley & Sons, Ltd.     

Forced convective flow and heat transfer of upward cocurrent air–water slug flow in vertical plain and swirl tubes

This experimental study comparatively examined the two-phase flow structures, pressured drops and heat transfer performances for the cocurrent air–water slug flows in the vertical tubes with and without the spiky twisted tape insert. The two-phase flow structures in the plain and swirl tubes were imaged using the computerized high frame-rate videography with the Taylor bubble velocity measured. Superficial liquid Reynolds number (Re_L) and air-to-water mass flow ratio (AW), which were respectively in the ranges of 4000–10000 and 0.003–0.02 were selected as the controlling parameters to specify the flow condition and derive the heat transfer correlations. Tube-wise averaged void fraction and Taylor bubble velocity were well correlated by the modified drift flux models for both plain and swirl tubes at the slug flow condition. A set of selected data obtained from the plain and swirl tubes was comparatively examined to highlight the impacts of the spiky twisted tape on the air–water interfacial structure and the pressure drop and heat transfer performances. Empirical heat transfer correlations that permitted the evaluation of individual and interdependent Re_L and AW impacts on heat transfer in the developed flow regions of the plain and swirl tubes at the slug flow condition were derived.     

Effect of heat transfer on the flow of a second‐grade fluid in divergent/convergent channel

This paper studies the effects of a second‐grade fluid on the flow and heat transfer characteristics in a divergent/convergent channel. The momentum and energy equations are first given in a nondimensional form and then solved analytically using the method of homotopy analysis method. Convergence of derived series solutions is shown. Graphical results for the velocity and the temperatures are presented and discussed for various emerging parameters. Copyright © 2009 John Wiley & Sons, Ltd.