Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 1. Isothermal flows

This paper presents a review of authors' collective works in the field of two-phase flow modeling done in the past few decades. The paper is aimed at the construction of mathematical models for simulation of particle-laden turbulent flows. A kinetic equation was obtained for the probability density function (PDF) of the particle velocity distribution in turbulent flows. The proposed kinetic equation describes both the interaction of particles with turbulent eddies of the carrier phase and particle-particle collisions. This PDF equation is used for the derivation of different schemes describing turbulent momentum transfer in the dispersed particle phase. The turbulent characteristics of the gaseous phase are calculated on the basis of the k - turbulence model with a modulation effect of particles on the turbulence.

The constructed models have been applied to the calculation of various two-phase gas-particle turbulent flows in jets and channels as well as particle deposition in tubes and separators. For validating the theoretical and numerical results, a wide range of comparisons with experimental data from Russian and foreign sources has been done.     

Heat and mass transfer in a system of plane turbulent jets with diffusion combustion

We consider turbulent motion of premixed chemically active gases in an infinite system of plane turbulent jets in the presence of diffusive combustion. The proposed calculation method permits determining the distribution of all the parameters in the mixing zone, including the longitudinal pressure. Numerical examples of the calculation of hydrogen combustion in air are presented.The study of heat and mass transfer in jet flows presents major difficulties at the present time. Therefore all the existing methods for calculating jet flows with heat and mass transfer and chemical processes [1–5] are based on an extension of the known semiempirical theories of free turbulence to the more complex cases of flow with chemical reactions. The present study is no exception in this sense; it covers an investigation of the motion in an infinite system of plane turbulent jets with diffusive combustion.     

Periodic flow and heat transfer using unstructured meshes

A numerical scheme has been developed for computing fluid flow and heat transfer in periodically repeating geometries. Unstructured solution-adaptive meshes are used in a cell-centred finite volume formulation. The SIMPLE algorithm is used for pressure‒velocity coupling. For periodic flows the static pressure is decomposed into a periodic component and one that varies linearly in the streamwise direction. The latter is computed from the imposition of overall mass balance at the periodic boundary. A subiteration between the periodic pressure correction equation and the correction to the linear component is used. For heat transfer a formulation using the physical rather than the scaled temperature is employed. The scheme is applied to both laminar and turbulent computations of periodic flow and heat transfer in a variety of heat exchanger geometries; comparison with published computations and experimental data is found to be satisfactory. © 1997 John Wiley & Sons, Ltd.     

A separated flow model for predicting two-phase pressure drop and evaporative heat transfer for vertical annular flow

A separated flow model has been developed that is applicable to vertical annular two-phase flow in the purely convective heat transfer regime. Conservation of mass, momentum, and energy are used to solve for the liquid film thickness, pressure drop, and heat transfer coefficient. Closure relationships are specified for the interfacial friction factor, liquid film eddy-viscosity, turbulent Prandtl number, and entrainment rate. Although separated flow models have been reported previously, their use has been limited, because they were tested over a limited range of flow and thermal conditions. The unique feature of this model is that it has been tested and calibrated against a vast array of two-phase pressure drop and heat transfer data, which include upflow, downflow, and microgravity flow conditions. The agreements between the measured and predicted pressure drops and heat transfer coefficients are, on average, better or comparable to the most reliable empirical correlations. This separated flow model is demonstrated to be a reliable and practical predictive tool for computing two-phase pressure drop and heat transfer rates. All of the datasets have been obtained from the open literature.     

Calculation of the momentum and heat transfer in turbulent gas-particle jet flows

The equations for the second moments of the dispersed-phase velocity and temperature fluctuations are used for calculating gas-suspension jet flows within the framework of the Euler approach. The advantages of introducing the equations for the second moments of the particle velocity fluctuations has previously been quite convincingly demonstrated with reference to the calculation of two-phase channel boundary flows [9–11]. The flows considered below have a low solid particle volume concentration, so that interparticle collisions can be neglected and, consequently, the stochastic motion of the particles is determined exclusively by their involvement in the fluctuating motion of the carrier flow. In addition to the equations for the turbulent energy of the gas and its dissipation, the calculation scheme includes the equations for the turbulent energy and turbulent heat transfer of the solid phase; however, the model constructed does not contain additional empirical constants associated with the presence of the particles in the flow.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.3, pp. 69–80, May–June, 1992.     

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.     

An analytical wall-function for recirculating and impinging turbulent heat transfer

The performance of the analytical wall-function (AWF) of Craft et al. [Craft, T.J., Gerasimov, A.V., Iacovides, H., Launder, B.E., 2002, Progress in the generalisation of wall-function treatments. Int. J. Heat Fluid Flow 23, 148–160.] is improved for predicting turbulent heat transfer in recirculating and impinging flows. Since constant parameters of the eddy viscosity formula were used to derive the AWF, the prediction accuracy of the original AWF tends to deteriorate in complex flows where those parameters need changing according to the local turbulence. To overcome such shortcomings, the present study introduces a functional behaviour on the strain parameter into the coefficient of the eddy viscosity of the AWF. The presently modified version of the AWF is validated in turbulent heat transfer of pipe flows, channel flows, back-step flows, pipe flows with abrupt expansion and plane impinging slot jets. The results confirm that the present modification successfully improves the performance of the original AWF for all the flows and heat transfer tested.     

Thermographic heat transfer measurements in separated flows

A measurement technique to determine the surface heat transfer distribution in complex turbulent flows is described. For this purpose, a constant wall heat flux test surface has been designed. To measure the surface temperature of the test plate, an infrared camera was used. The instrumentation allows the determination of the heat transfer with high accuracy and detailed spatial resolution. In examining combustor-type separated flow, the capabilities of the technique are demonstrated and its accuracy is verified by appropriate conventional techniques.     

Turbulent mixed convection heat transfer to liquid sodium

The influences of buoyancy on turbulent heat transfer to a liquid metal flowing in a vertical pipe are considered. A theoretical model is presented which provides a criterion for the conditions under which such influences become significant and which predicts the impairment of heat transfer for upward flow and enhancement for downward flow. The variation with Peclet number of the maximum impairment of heat transfer and conditions under which it occurs are established. A generalization of the model leading to an equation for the entire mixed convection region is proposed. From this an equation for turbulent free convection to liquid metals is obtained.     

Heat transfer characteristics of water and APG surfactant solution in a micro-channel heat sink

The steady increase in internal heat production of cost and high performance electronic components has lead researchers to seek improved ways to remove the heat generated. Single-phase liquid flow has been considered as a potential solution for solving this cooling problem. However, when considering that any solution needs to be of low cost and low mass fluxes and yet retain low temperature gradients across the electronic components, it seems that two-phase boiling flow is preferred. Surfactant solutions have been introduced in connection with enhancement of the boiling processes. We investigated the effects of surfactant solution flows through a micro-channel heat sink. The experimental setup included a high-speed IR radiometer and a CCD camera that were used to characterize the test module. The module consisted of inlet and outlet manifolds that distributed surfactant solutions through an array of 26 parallel micro-channels. The experimental results have shown that there exists an optimal solution concentration and mass flux for enhancing heat removal. Surfactant solution boiling flows were also found to stabilize the maximum and average surface temperatures for a wide range of applied heat fluxes. In addition, the use of surfactant solutions at low mass fluxes has led to CHF enhancement when compared to regular water flows. In the last part of this work, possible explanations for the observed non-ionic surfactant effects are presented.     

Effects of free stream turbulence on the distribution of heat transfer around turbine blade sections

Local convective heat transfer coefficients to a number of modern gas turbine blade sections have been measured under a wide range of mainstream conditions, from notionally steady flows to highly perturbed turbulent flows. The paper discusses the results and, through a detailed analysis of the pertinent boundary layer flow parameters and their relation to the observed experimental results, tests criteria for the occurrence of transition from laminar to turbulent boundary layers, a factor which all the data from this work confirm as critical in predicting the quantitative effects of mainstream turbulence on heat transfer rates. Artificially induced mainstream turbulence, which is endemic in the flows in a real turbine, enhances significantly the heat transfer rates, especially to the leading edge regions and on the pressure surface, particularly when the acceleration is tending to suppress transition. The results presented here confirm existing criteria for laminarisation and the applicability of some of those available for predicting laminar-turbulent transition. The observations also demonstrate how surface geometry can influence the stability of the flows, and the uncertainties which remain in assessing the effect of Goertler vortices and their role in the convective heat transfer process.     

An analytical skin friction and heat transfer model for compressible, turbulent, internal flows

An analytical skin friction model for compressible, turbulent, internal, fully developed flow involving adiabatic and non-adiabatic, smooth and rough flows has been developed by extending the incompressible law-of-the-wall relation to compressible cases. The formula recovers Prandtl's incompressible law of friction for pipes (within 2%) for incompressible flow. The model also shows good correlation with available data for compressible, adiabatic flows and flows involving cold wall heat transfer (within 15%). Comparison with hot wall data is only moderate (15–30%). Finally, using Reynold's analogy, the Stanton number and Nusselt numbers may be estimated.     

Pulsating flow in a heat exchanger

The problem of heat transfer in industrial processes, heat exchangers, and combustion chambers is formulated for a case where flow inside the chamber consists of a periodic motion imposed on a fully developed turbulent flow. It is shown that the velocity pulsations induce harmonic oscillations in temperature, thus breaking the temperature field into a steady mean part and a harmonic part. The interaction between the velocity and temperature oscillations introduces an extra term into the energy equation which reflects the effect of pulsations in producing higher heat transfer rates. The analysis shows that when the mean temperature is fully developed with constant heat flux at the wall, there is no effect of the velocity pulsations on the total heat transfer rate along the chamber. For the case where the mean temperature profile is not fully developed, analytical solutions are obtained for asymptotic values of the pulsations frequency. The results show the temperature gradient and its dependence on the frequency. These results are used to evaluate the feasibility of pulsating the flow in a heat exchanger for obtaining higher rates of heat transfer.     

Stationary modes of combustion in a fine-scale turbulent stream

The majority of models of the turbulent combustion of gases are based mainly on intuitive concepts concerning the processes occurring in the flame. The characteristics of a turbulent flame are estimated from considerations of dimensionality and similarity. A detailed review of works on turbulent combustion is given in [1]. Problems on the calculation of the combustion rate in a turbulent stream as a proper value of the equations of heat and mass transfer and of the corresponding boundary conditions have recently been raised. Here too one must rest on assumptions of a semiempirical nature, which in large measure is connected with the inadequate level of development of turbulence theory. In the present work the equation of propagation of the zone of chemical reactions in the stream is averaged statistically by analogy with studies of turbulent flows. Correct averaging is possible at scales of hydrodynamic disturbances smaller than the flame thickness (fine-scale turbulence). The temperature pulsations are related with the size of the heat flux using the theory of mixing lengths. The main influence is specific to effects arising during averaging of the heat release function. Two stationary modes, distinguished by the normal propagation velocity ₁, are isolated within the framework of the Cauchy problem with a given initial mixture temperature and zero heat flux in the burned gas. A heat conduction mode occurs with a stream velocity > ₁ and an induction mode with < ₁. An expression is found for ₁ which reflects the principal effects in the flame and which in the limit coincides with the equation of Zel'dovich and Frank-Kamenetskii for a laminar flame. In those cases when the distorting effect of the heat release function is small, the turbulence affects the combustion rate through mechanisms of intensification of transport processes.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 118–124, September–October, 1973.     

Effects of velocity fluctuations on heat transfer enhancement

The three-dimensional velocity fluctuation effects on heat transfer enhancement were experimentally investigated using a wind tunnel system and cylinders placed upstream of the test section in the wind tunnel. The cylinders with different diameters were used as turbulators to generate vortical flow motions with three-dimensional velocity fluctuations. A heated plate, part of the tunnel wall, was placed far downstream of the cylinders such that it was subjected mainly to flows with velocity fluctuations but with negligible steady vortical motions. These studies included three-component velocity measurements to characterize the near-wall and cross-section velocity fields and to obtain the turbulent kinetic energy. The temperatures were measured by thermocouples on the heated plate to obtain the convection heat transfer coefficients and the Nusselt numbers. Results indicate that the heat transfer was enhanced by the velocity fluctuations, which is attributed to the modification of boundary layer velocity profiles without the modification of boundary layer thickness. The resulting normalized Nusselt number was approximately a parabolic function of a dimensionless parameter, the product of Reynolds number and normalized turbulent kinetic energy.     

Analysis of Mixing Models for use in Simulations of Turbulent Spray Combustion

The present study concerns the investigation of different mixing models for use in the transported probability density function (PDF) modeling of turbulent (reacting) spray flows. The modeling of the turbulent mixing and other characteristic scalar variables such as gas enthalpy using transported (joint) PDFs has become an important method to describe turbulent (reacting) spray flows since the evaporation process causes the PDF of the mixture fraction to deviate from the widely used β function, which is typically used in models for turbulent gas flows. In the PDF transport equation, the molecular mixing does not appear in closed form so that modeling strategies are required. For gas combustion, the interaction-by-exchange-with-the-mean (IEM) model, the modified Curl (MC) model, and the Euclidean minimum spanning tree (EMST) models are used. More recently, a new mixing model, the PSP model, which is based on parameterized scalar profiles has been developed. The present study focuses on the use and analysis of the IEM, MC and PSP models for turbulent spray flames. For this purpose, the models are reconsidered with respect to the evaporation process that must be included and evaluated if spray combustion is considered. For model evaluation, turbulent ethanol/air spray flames are simulated, and the results are compared to experimental data by A. Masri, University of Sydney, Australia.     

Condensation heat transfer characteristics of R-22, R-134a and R-410A in small diameter tubes

The condensation heat transfer of pure refrigerants, R-22, R-134a and a binary refrigerant R-410A flowing in small diameter tubes was investigated experimentally. The condenser is a countflow heat exchanger which refrigerant flows in the inner tube and cooling water flows in the annulus. The heat exchanger is smooth, horizontal copper tube of 1.77, 3.36 and 5.35 mm inner diameter, respectively. The length of heat exchanger is 1220, 2660 and 3620 mm, respectively. The experiments were conducted at mass flux of 200–400 kg/m² s and saturation temperature of 40°C. The main results were summarized as follows: in case of single-phase flow, the single-phase Nusselt Number measured by experimental data was higher than that calculated by Gnielinski and Wu and Little correlation. The new single-phase correlation based on the experimental data was proposed in this study. In case of two-phase flow, the condensation heat transfer coefficient of R-410A for three tubes was slightly higher than that of R-22 and R-134a at the given mass flux. The condensation heat transfer coefficient of R-22 showed almost a similar value to that of R-134a. The condensation heat transfer coefficient for R-22, R-134a and R-410A increased with increasing mass flux and decreasing tube diameter. Most of the existing correlations which were proposed in the large diameter tube failed to predict condensation heat transfer. Therefore, the new condensation heat transfer correlation based on the experimental data was proposed in the present study.     

Effect of inertial particles with different specific heat capacities on heat transfer in particle-laden turbulent flow

The effect of inertial particles with different specific heat on heat transfer in particle-laden turbulent channel flows is studied using the direct numerical simulation(DNS) and the Lagrangian particle tracking method. The simulation uses a two-way coupling model to consider the momentum and thermal interactions between the particles and turbulence. The study shows that the temperature fields display differences between the particle-laden flow with different specific heat particles and the particle-free flow,indicating that the particle specific heat is an important factor that affects the heat transfer process in a particle-laden flow. It is found that the heat transfer capacity of the particle-laden flow gradually increases with the increase of the particle specific heat. This is due to the positive contribution of the particle increase to the heat transfer. In addition,the Nusselt number of a particle-laden flow is compared with that of a particle-free flow.It is found that particles with a large specific heat strengthen heat transfer of turbulent flow, while those with small specific heat weaken heat transfer of turbulent flow.     

Numerical analysis of swirling turbulent droplet-laden flow and heat transfer in a sudden pipe expansion

The effect of swirling intensity on the structure and heat transfer of a turbulent gas–droplet flow after a sudden pipe expansion has been numerically simulated. Air is used as the carrier phase, and water, ethanol, and acetone are used as the dispersed phase. The Eulerian approach is applied to simulate the dynamics and heat transfer in the dispersed phase. The gas phase is described by a system of Reynolds-averaged Navier-Stokes (RANS) equations, taking into account the effect of droplets on mean transport and turbulent characteristics in the carrier phase. Gas phase turbulence is predicted using the second-moment closure. A swirling droplet-laden flow is characterized by an increase in the number of small particles on the pipe axis due to their accumulation in the zone of flow recirculation and the action of the turbulent migration (turbophoresis) force. A rapid dispersion of fine droplets over the pipe cross-section is observed without swirling. With an increase in swirling intensity, a significant reduction in the length of the separation region occurs. The swirling of a two-phase flow with liquid droplets leads to an increase in the level of turbulence for all three types of liquid droplets investigated in this work due to their intensive evaporation. It is shown that the addition of droplets leads to a significant increase in heat transfer in comparison with a single-phase swirling flow. The greatest effect of flow swirling on heat transfer intensification in a two-phase gas-droplet flow is obtained for the droplets of ethanol and water and smallest effect is for the acetone droplets.