Numerical simulation of the modulated flow pattern for vertical upflows by the phase separation concept

The multiphase heat transfer could be enhanced by creating thin liquid film on the wall. The phase separation concept is called due to the separated flow paths of liquid and gas over the tube cross section to yield thin liquid film. Our proposed heat transfer tube consists of an annular region close to the wall and a core region, interfaced by a suspending mesh cylinder in the tube. The heat transfer tube is a multiscale system with micron scale of mesh pores, miniature scale of annular region and macroscale of tube diameter and length. Great effort has been made to link from micron scale to macroscale. The Volume of Fluid (VOF) method simulates air/water two-phase flow for vertical upflow. The three-dimensional system was successfully converted to a two-dimensional one by using three equivalent criteria for mesh pores. The non-uniform base grid generation and dynamic grid adaption method capture the bubble interface. The numerical results successfully reproduce our experimental results. The numerical findings identify the following mechanisms for the enhanced heat transfer: (a) counter-current flow exists with upward flow in the annular region and downward flow in the core region; (b) void fractions are exact zero in the core region and higher in the annular region; (c) the liquid film thicknesses are decreased to 1/6–1/3 of those in the bare tube section; (d) the gas–liquid mixture travels much faster in the annular region than in the bare tube; (e) three-levels of liquid circulation exists: meter-scale bulk liquid circulation, moderate-scale liquid circulation around a single-elongated-ring-slug-bubble, and microliquid circulation following the ring-slug-bubble tails. These liquid circulations promote the fluid mixing over the whole tube length and within the radial direction. The modulated parameters of void fractions, velocities and liquid film thicknesses in the annular region and three-levels of liquid circulation are greatly beneficial for the multiphase heat transfer enhancement.     

Numerical modeling of cryogenic chilldown process in terrestrial gravity and microgravity

In this paper, a numerical model is developed to predict the cryogenic chilldown process of a vertical tube for both terrestrial and microgravity conditions. The flow field is covered by four distinct regions, which are single-phase vapor region, dispersed flow region, inverted annular flow region, and single-phase liquid region. Heat transfer mechanisms are dictated by the flow patters. A two-fluid model is employed to analyze the dispersed flow region and the inverted annular film boiling region. Gravity effect on the chilldown process is also investigated. Model results indicate that film boiling heat transfer decreases with decreasing gravity level for the bottom flooding condition. The model results show a good agreement with the experimental data.     

Experimental study of R134A condensation heat transfer inside the horizontal micro-fin tubes

The Investigation of the two-phase flow patterns and their transitions during the condensation has gained increasing interest and importance from the well-known phenomenon that the heat transfer characteristics are strongly dependent on the flow patterns. Therefore, it is very important to study on which heat transfer enhancement approach is suitable for an individual flow pattern inside a condenser, so that an accurate heat transfer mechanism can be understood that is consistent with the flow patterns. The condensation heat transfer for R134a in the two kinds of in-tube three-dimensional (3-D) micro-fin tubes with different geometries is experimentally investigated. Based on the flow pattern observations, the flow patterns in the Soliman flow regime map are divided into two-flow regimes; one with the vapor-shear-dominant annular regime and the other with the gravitational-force-dominant stratified-wavy regime. The flow regime transition criterion between the annular regime and the stratified-wavy regime is at Fr equal to 2. In the annular regime, the heat transfer coefficients h of the two kinds of in-tube 3-D micro-fin tubes decreases as the vapor quality x decreases. The regressed condensation heat transfer correlation from the experimental data of the annular flow region is obtained. The dispersibility of the experimental data is inside the limits of ±25%. In the stratified-wavy regime, the average heat transfer coefficient h of the two kinds of in-tube 3-D micro-fin tubes increases as the mass flux increases and the number of micro fins in the 3-D micro-fin tube is not the controlling factor for the performance of a condensation heat transfer. The regressed condensation heat transfer correlation of the stratified-wavy flow regime is experimentally obtained. The dispersibility of the experimental data is inside the limits of ±22%. Combined with the criteria of flow pattern transitions, the correlations can be used for the design of a condenser with 3-D micro-fin tubes.     

Enhancement of turbulent heat transfer during interaction of an impinging axisymmetric mist jet with a target

The flow structure and heat transfer of a mist jet with a low mass concentration of droplets (within 1%) impinging onto a flat surface aligned normal to the jet are studied numerically. The mathematical model is based on solving a system of Reynolds-averaged Navier-Stokes equations for a two-phase flow with the kinetic equation of the probability density function for coordinates, velocity, and temperature of particles. Addition of droplets is demonstrated to enhance heat transfer substantially, as compared with an impinging single-phase air jet in the region directly adjacent to the stagnation point of the jet.     

Research on heat transfer of two-phase dispersed flow in narrow annular channel

Narrow channel heat transfer technique is a new developing heat transfer technique in recent years. As the temperature of droplet, steam and wall are decided by forced convection heat transfer between the steam and the wall, between the droplet and the wall, between the steam and the droplet and radiation heat transfer, which makes heat transfer mechanism of dispersed flow be difficultly interpretative. Dispersed flow in narrow annular channel is analyzed in the paper, investigating the influence of all kinds of heat transfer processes on dispersed flow, building annular channel dispersed flow model using thermodynamic non-equilibrium model. Calculation results show heat transfer is mainly controlled by heat transfer process between steam and wall. When temperature is low, radiation can be ignored on heat transfer coefficient calculation. The calculation of model can provide a reference for engineering application of steam generator, refrigeration system and so on.     

Experiments and homogeneous turbulence model of boiling flowin narrow channels

The boiling heat transfer experiments have been carried out in vertical narrow annular channels with pure water. A two-dimensional homogeneous turbulence model of boiling flow has been developed and solved numerically to yield pressure gradient, and velocity, thermal and turbulence fields, together with local heat transfer coefficient along the length of the tube. Predictions are compared with the data of experiments and agreed well with it. The model results show that the heat transfer coefficient increases as the gap size decreases in annular channels. This model can be used to predict heat transfer of boiling flow in narrow channels.     

Peristaltic flow and heat transfer in a vertical porous annulus, with long wave approximation

In this paper, we study the interaction of peristalsis with heat transfer for the flow of a viscous fluid in a vertical porous annular region between two concentric tubes. Long wavelength approximation (that is, the wavelength of the peristaltic wave is large in comparison with the radius of the tube) is used to linearise the governing equations. Using the perturbation method, the solutions are obtained for the velocity and the temperature fields. Also, the closed form expressions are derived for the pressure-flow relationship and the heat transfer at the wall. The effect of pressure drop on flux is observed to be almost negligible for peristaltic waves of large amplitude; however, the mean flux is found to increase by 10-12% as the free convection parameter increases from 1 to 2. Also, the heat transfer at the wall is affected significantly by the amplitude of the peristaltic wave. This warrants further study on the effects of peristalsis on the flow and heat transfer characteristics.     

Algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow, focusing in particular on momentum and heat transfer (so-called ‘convective boiling’) through the annular liquid film. In contrast with single-phase wall-bounded flow theory, shear-driven annular liquid films are assumed here to behave as fluid-bounded flows, mostly interacting with the shearing gas-entrained droplets core flow. Besides providing velocity and temperature profiles through the liquid film, the turbulence model proposed here predicts key parameters such as the average liquid film thickness, the void fraction and the convective boiling heat transfer coefficient with accuracies comparable or better than those of leading design correlations. This turbulence model is part of a unified annular flow modeling suite that includes methods to predict the entrained liquid fraction and the axial frictional pressure gradient. The underlying heat transfer database covers nine fluids (water, two hydrocarbons and six refrigerants) for vertical and horizontal tubes of 1.03-14.4 mm diameter and pressures of 0.1-7.2 MPa. Importantly, this study shows that there appears to be no macro-to-microscale transition when it comes to annular flow. Simply better physical modeling is required to span this range.     

Pressure drop and heat transfer characteristics of turbulent flow in annular tubes with internal wave-like longitudinal fins

Measured were pressure drop and heat transfer characteristics with uniform axial heat input using air as the working fluid in both the entrance and fully developed regions of annular tubes with wave-like longitudinal fins. Five series of experiments were performed for turbulent flow and heat transfer in the annular tubes with number of waves equal to 4, 8, 12, 16 and 20, respectively. The test tube has a double-pipe structure with the inner blocked tubes as an insertion. The wave-like fins are in the annulus and span its full width. The friction factor and Nusselt number in the fully developed region were obtained. The friction factor and Nusselt number can be well corrected by a power-law correction in the Reynolds number range tested. In order to evaluate the thermal performance of the longitudinal finned tubes over a plain circular tube, comparisons were made under three conditions: (1) identical pumping power; (2) identical pressure drop and (3) identical mass flow. It was found that under the three constraints all the wave-like finned tubes can enhance heat transfer with the tube with wave number 20 being superior. Finally, discussion on the enhancement mechanism is conducted and a general correlation for the fully developed heat transfer is provided, which can cover all the fifty data of the five tubes with a mean deviation of 9.3%.     

Analytical considerations of flow boiling heat transfer in metal-foam filled tubes

Flow boiling in metal-foam filled tube was analytically investigated based on a modified microstructure model, an original boiling heat transfer model and fin analysis for metal foams. Microstructure model of metal foams was established, by which fiber diameter and surface area density were precisely predicted. The heat transfer model for flow boiling in metal foams was based on annular pattern, in which two phase fluid was composed by vapor region in the center of the tube and liquid region near the wall. However, it was assumed that nucleate boiling performed only in the liquid region. Fin analysis and heat transfer network for metal foams were integrated to obtain the convective heat transfer coefficient at interface. The analytical solution was verified by its good agreement with experimental data. The parametric study on heat transfer coefficient and boiling mechanism was also carried out.     

Visualisation of flow boiling heat transfer in a microtube

This paper presents the results of the flow boiling patterns and heat transfer coefficients of FC-72 in a small tube. The internal diameter of the tube is 0.48 mm, with a heated length of 73 mm. The mass flow rate varies from 50 to 3,000 kg/m²-s. The microtube is made of Pyrex in order to obtain the visualisation of the flow pattern along the heated channel. Different types of flow pattern have been observed: bubbly flow, deformed bubbly flow, bubbly/slug flow, slug flow, slug/annular flow, and annular flow. The experiments show the presence of flow instabilities in a large portion of the tests at low mass flow rates and low subcooling. Flow patterns in presence of flow instabilities are mainly characterized by bubbly/slug flow and slug/annular flow. Heat transfer rates have been studied in all flow pattern conditions. The two groups of results, with flow instabilities and without flow instabilities, show similar heat transfer behaviour. The heat transfer characteristics of the pipes have been studied in comparison with mass flux and vapour quality.     

Radial variation of adiabatic and diabatic spiral vortex flow in a wide annular gap

Transitions occurring after the onset of spiral vortex flow in a wide concentric annular gap of radius ratio 0.8, formed by a stationary outer cylinder and a rotatable inner cylinder, have been studied experimentally. By isothermal heating of the annular surface, it was possible to consider diabatic as well as adiabatic conditions. At an axial Reynolds number of 500 and for a range of Taylor numbers up to 10⁷, power spectra and auto-correlograms were taken at two radial positions near to the inner and outer annular surfaces; these are compared with previous results taken at mid-gap under adiabatic conditions. Measurements of turbulence intensity across the gap were made also. Probability histograms and signal traces for diabatic flow near to the outer annular surface are presented. It has been shown that the vortex transitions affect the thermal boundary layer and, consequently, the heat transfer rates at the outer surface. A positive radial thermal gradient was seen to have little effect on the flow. The imposed axial flow was shown to be stabilising under adiabatic conditions but destabilising under diabatic conditions.     

Entropy analysis for third grade fluid flow with Vogel model viscosity in annular pipe

The steady-state flow of a third grade fluid between concentric circular cylinders is considered and entropy generation due to fluid friction and heat transfer in the annular pipe is examined. Depending upon the fluid viscosity, entropy generation in the flow system varies. The third grade fluid is employed to account for the non-Newtonian effect while Vogel model is accommodated for temperature-dependent viscosity. The analysis is based on perturbation technique. The closed form solutions for velocity, temperature and entropy fields are presented. Entropy generation due to fluid friction and heat transfer in the flow system is formulated. The influence of viscosity parameters A and B on the entropy generation number is investigated. It is found that entropy generation number reduces with increasing viscosity parameter A, which is more pronounced in the region close to the annular pipe inner wall and opposite is true for increasing viscosity parameter B.     

MHD liquid metal flow and heat transfer between vertical coaxial cylinders under horizontal magnetic field

Direct numerical simulations are presented of MHD liquid metal flow and heat transfer in vertical annuli. Three annular gaps and four ratios of annular height to annular gap are considered. The walls of the external and internal cylinders are isothermal with the temperature of the outer cylinder being higher and, thus, buoyancy is the driving force. The results show that the fluid motion increases as the aspect ratio and the annular gap become larger. The presence of the magnetic field results to fluid deceleration and, thus, to flow stabilization. Additionally, non symmetric flow patterns develop, due to the magnetic field, resulting in differently sized normal and parallel wall layers, namely the Hartmann and the Roberts layers, respectively. For all annular gaps considered, the highest spatially averaged heat transfer rates are obtained for aspect ratios equal to 1.     

Numerical simulation of flow boiling for organic fluid with high saturation temperature in vertical porous coated tube

A semi-analytical model is developed for the prediction of flow boiling heat transfer inside vertical porous coated tubes. The model assumes that the forced convection and nucleate boiling coexist together in the annular flow regime. Conservations of mass, momentum, and energy are used to solve for the liquid film thickness and temperature. The heat flux due to nucleate boiling consists of those inside and outside micro-tunnels. To close the equations, a detailed analysis of various forces acting on the bubble is presented to predict its mean departure diameter. The active nucleation site density of porous layer is determined from the pool boiling correlation by introducing suppression factor. The flow boiling heat transfer coefficients of organic fluid (cumene) with high saturation temperature in a vertical flame-spraying porous coated tube are studied numerically. It is shown that the present model can predict most of the experimental values within ±20%. The numerical results also indicate that the nucleate boiling contribution to the overall heat transfer coefficient decreases from 50% to 15% with vapor quality increasing from 0.1 to 0.5.     

Mathematical modeling of near-wall flows of two-phase mixture with evaporating droplets

ＩｎｔｒｏｄｕｃｔｉｏｎＯｐｔｉｍｉｚａｔｉｏｎｏｆｖａｒｉｏｕｓｔｈｅｒｍａｌｐｏｗｅｒｄｅｖｉｃｅｓ (ｓｔｅａｍｇｅｎｅｒａｔｏｒｓ,ｈｅａｔｅｘｃｈａｎｇｅｒｓ ,ｅｔｃ .)ａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｔｈｉｎｆｉｌｍｔｅｃｈｎｏｌｏｇｉｅｓ(ｆｏｒｅｘａｍｐｌｅ ,ｗｉｔｈｔｈｅｕｓｅｏｆｔｗｏ＿ｐｈａｓｅｊｅｔｓ)ｐｒｏｍｐｔｓｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｉｎｇｏｆｎｅａｒ＿ｗａｌｌｆｌｏｗｓｏｆｇａｓ＿ｄｒｏｐｌｅｔｍｉｘｔｕｒｅｓ .Ｄｕｅｔｏｔｈｅｇｒｅａｔｄｉｖｅｒｓｉｔｙｉｎ…     

A model for droplet entrainment in heated annular flow

The ability to accurately predict droplet entrainment in annular two-phase flow is required to effectively calculate the interfacial mass, momentum, and energy transfer, which characterizes nuclear reactor safety, system design, analysis, and performance. Most annular flow entrainment models in the open literature are formulated in terms of dimensionless groups, which do not directly account for interfacial instabilities. However, many researchers agree that there is a clear presence of interfacial instability phenomena having a direct impact on droplet entrainment. The present study proposes a model for droplet entrainment, based on the underlying physics of droplet entrainment from upward co-current annular film flow that is characteristic to light water reactor safety analysis. The model is developed based on a force balance and stability analysis that can be implemented into a transient three-field (continuous liquid, droplet, and vapor) two-phase heat transfer and fluid flow systems analysis computer code.     

Wärmeübergang bei turbulenter Ringspaltströmung und konstanter thermischer Randbedingung zweiter Art

An analysis has been carried out to determine the heat transfer characteristics for fully developed turbulent flow in concentric annular ducts with constant thermal boundary conditions of second kind. The energy equation has been solved analytically by separation of variables. Results for heat transfer in the thermal entrance region, in the thermal fully developed flow and thermal entrance-length results are presented over a wide range of Reynolds number (10⁴?Re?10⁶) and Prandtl-number (0?Pr?100) and compared to available empirical data.     

Heat transfer in laminar and turbulent flows in the thermal entrance region of concentric annuli: Axial heat conduction effects in the fluid

The present analytical study investigates the influence of axial heat conduction within the flow on the heat transfer in the thermal entrance region of a concentric annular duct with laminar and turbulent internal flow. The solution is based on a decomposition of the elliptic energy equation into a pair of first order partial differential equations. By using a new defined vector norm it is possible to obtain a selfadjoint eigenvalue problem for the extended Graetz problem even though the original convective diffusion operator is non-selfadjoint. The obtained exact analytical solutions for the Graetz problem with axial heat conduction are as simple to compute as the related solutions of the parabolic problem. Received on 28 October 1996     

Frictional pressure drop during vapour–liquid flow in minichannels: Modelling and experimental evaluation

Condensation in minichannels is widely used in air-cooled condensers for the automotive and air-conditioning industry, in heat pipes and other applications for system thermal control. The knowledge of pressure drops in such small channels is important in order to optimize heat transfer surfaces. This paper presents a model for calculation of the frictional pressure gradient during condensation or adiabatic liquid–gas flow inside minichannels with different surface roughness. In order to account for the effects of surface roughness, new experimental frictional pressure gradient data associated to single-phase flow and adiabatic two-phase flow of R134a inside a single horizontal mini tube with rough wall has been used in the modelling. It is a Friedel (1979) [Friedel, L., 1979. Improved friction pressure drop correlations for horizontal and vertical two-phase pipe flow. In: Proceedings of the European Two-Phase Flow Group Meeting, Ispra, Paper E2] based model and it takes into account mass velocity, vapor quality, fluid properties, reduced pressure, tube diameter, entrainment ratio and surface roughness. With respect to the flow pattern prediction capability, it has been built for shear dominated flow regimes inside pipes, thus, annular, annular-mist and mist flow are here predicted. However, the suggested procedure is extended to the intermittent flow in minichannels and it is also applied with success to horizontal macro tubes.