Developing and fully developed turbulent flow in ribbed channels

Wall-mounted roughness features, such as ribs, are often placed along the walls of a channel to increase the convective surface area and to augment heat transfer and mixing by increasing turbulence. Depending on the relative roughness size and orientation, the ribs also have varying degrees of increased pressure losses. Designs that use ribs to promote heat transfer encompass the full range of having only a few streamwise ribs, which do not allow fully developed flow conditions, to multiple streamwise ribs, which do allow the flow to become fully developed. The majority of previous studies have focused on perturbing the geometry of the rib with little attention to the spatially and temporally varying flow characteristics and their dependence on the Reynolds number. A staggered rib-roughened channel study was performed using time-resolved digital particle image velocimetry (TRDPIV). Both the developing (entry region) and a fully developed region were interrogated for three Reynolds numbers of 2,500, 10,000, and 20,000. The results indicate that the flow was more sensitive to Reynolds number at the inlet than within the fully developed region. Despite having a similar mean-averaged flowfield structure over the full Reynolds number range investigated, the population and distribution of coherent structures and turbulent dissipation within the fully developed region were also found to be Reynolds number dependent. Exploring the time-accurate flow characteristics revealed that in addition to vortices shed from the rib shear layer, the region of the rib wake was governed by a periodic process of bursting of the wake vortices resulting in the intermittent ejection of the inter-rib recirculation region into the core flow. This periodic process was the driving mechanism resulting in mixing and heat transfer augmentation. A quadrant-splitting burst analysis was also performed to determine the characteristic frequency and duration of inter-rib bursting as well as the wake shedding frequency, both of which were determined to be Reynolds number dependent.     

Flooded countercurrent two-phase flow in horizontal tubes and channels

A criterion for flooding in the countercurrent flow of two fluids in horizontal tubes and channels is developed. It exhibits a dependence upon the density ratio of the fluids beyond that present in the criterion of Wallis (1969). Experiments were carried out with air and water in a horizontal tube at atmospheric pressure and these, together with others reported in the literature, are shown to be in fair agreement with the prediction of the criterion, though it is emphasised that more experimental work is desirable.Work reported in the literature with miscible fluids with a density ratio close to unity confirms the extra dependence upon the density ratio.     

Influence of particles on turbulent flow in channels

In this work, the equations of balance of the second single-point moments of pulsations of the carrier phase are used to analyze the influence of particles on the intensity of pulsation motion. Besides dissipation due to pulsation phase slip, these equations take account the effect of dissipation of small-scale vortices on the particles and also of the additional transfer of pulsations by particles due to the particles being drawn into the pulsating motion and the migration of particles across the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 40–48, January–February, 1985.     

Turbulent friction in the flow of liquid in tubes and channels

A calculating relationship is presented for turbulent flow; it takes a unique form over the whole cross section of the flow. A relationship is also derived between turbulent friction and the mean velocity profile on the basis of the equation for the maximum turbulent friction, which follows directly from the equation of motion. The proportionality factor in this relationship is obtained with due allowance for twelve boundary conditions relating to the turbulent flow, the mean velocity, and their derivatives. The resultant turbulent-friction profiles agree with the experimental data of Laufer. The profile parameters may be related to the Reynolds number.Leningrad. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 140–145, March–April, 1972.     

Development of secondary free-convection currents in forced turbulent flow in horizontal tubes

The problem of the development of secondary free-convection currents in forced turbulent flow in horizontal tubes for relatively weak thermal gravitation influence is analytically solved. The results of the solution are compared to experimental data.     

Simulation of turbulent flow and heat transfer in channels between rod bundles

The turbulent flow and heat transfer in triangular rod bundles are investigated theoretically with CFD code FLUENT. The unsteady Reynolds Stress Model is adopted as turbulence modeling. The wall function is used for near wall boundary layer. The calculation results were in agreement with experimental data. The effects of the Reynolds number and pitch to diameter ratio on the flow and heat transfer in the lattice are significant. The traditional theoretical models could not predict the flow and heat transfer in the lattice. The P/D = 1.03 is a critical point. In this case, the flow and heat transfer in the lattice is the most desirable and most efficient, and the nuclear power could also reach its maximum. The variation of large scale coherent structure with pitch to diameter ratio is consistent with the variation of the Nusselt number with pitch to diameter ratio.     

Experimental and numerical study of turbulent flow and heat transfer inside hexagonal duct

Flow and heat transfer characteristics in transition and turbulent regions are studied experimentally and numerically in a horizontal smooth regular hexagonal duct under constant wall temperature boundary condition covering a range of Reynolds number from 2.3 × 10³ to 52 × 10³. Two types of k-omega (standard and shear stress transport (SST)) and three types of k-ε (standard, renormalization (RNG), and realizable) turbulence model are employed for transition and turbulent regions, respectively. Both average and fully developed Darcy friction factor and Nusselt number are presented as a function of Reynolds number. It is seen that k-omega SST and k-ε realizable turbulence models gave the best agreement with the experimental data in transition and turbulent regions, respectively. All the experimental results are correlated within an accuracy of ±13 % and ±7 % for Nusselt number and Darcy friction factor, respectively. Results obtained in this study are compared with circular duct results using hydraulic diameter.     

Large-eddy simulation of turbulent flow and heat transfer in plane and rib-roughened channels

Large-eddy simulation results are presented and discussed for turbulent flow and heat transfer in a plane channel with and without transverse square ribs on one of the walls. They were obtained with the finite-difference code Harwell-FLOW3D, Release 2, by using the PISOC pressure-velocity coupling algorithm, central differencing in space, and Crank-Nicolson time stepping. A simple Smagorinsky model, with van Driest damping near the walls, was implemented to model subgrid scale effects. Periodic boundary conditions were imposed in the streamwise and spanwise directions. The Reynolds number based on hydraulic diameter (twice the channel height) ranged from 10 000 to 40 000. Results are compared with experimental data, k-? predictions, and previous large-eddy simulations.     

Variation of fiber orientation in turbulent flow inside a planar contraction with different shapes

In this study, we have investigated the influence of shape of planar contractions on the orientation distribution of stiff fibers suspended in turbulent flow. To do this, we have employed a model for the orientational diffusion coefficient based on the data obtained by high-speed imaging of suspension flow at the centerline of a contraction with flat walls. This orientational diffusion coefficient depends only on the contraction ratio and turbulence intensity. Our measurements show that the turbulence intensity decays exponentially independent of the contraction angle. This implies that the turbulence variation in the contraction is independent of the shape, consistent with the results by the rapid distortion theory and the experimental results of axisymmetric contractions. In order to determine the orientation anisotropy, we have solved a Fokker–Planck type equation governing the orientation distribution of fibers in turbulent flow. Although the turbulence variation and the orientational diffusion are independent of the contraction shape, the results show that the variation of the orientation anisotropy is dependent on shape. This can be explained by the variation of the rotational Péclet number, Pe_r, inside the contractions. This quantity is a measure of the importance of the mean rate of the strain relative to the orientational diffusion. We have shown that when Pe_r < 10 turbulence can significantly influence the evolution of the orientation anisotropy. Since in contractions with identical inlet conditions the streamwise position where Pe_r = 10 depends on the shape, the orientation anisotropy is dependent on the variation of rate of strain in a given contraction. We demonstrate the shape effect by considering contraction with flat walls as well as three contractions with different mean rate of strain variation.     

Two-phase flow behavior inside a header connected to multiple parallel channels

The main objective of this work is to examine the flow distribution of two-phase mixture to parallel channels and to investigate the flow behavior at header-channel junctions simulating the corresponding parts of compact heat exchangers. The cross-section of the header and the channels were fixed to 14 mm × 14 mm and 12 mm × 1.6 mm, respectively. The mass flux and the mass quality ranges were 70–165 kg/m² s and 0.3–0.7, respectively. Air and water were used as the test fluids. The flow distribution at the fore part of the header (region A) is affected only by the upstream flow configuration and the rate of liquid flow separation decreased a flowing downwards. On the other hand, in the rear part, the downstream effect predominates over the upstream effect due to strong flow recirculation near the end plate. In this part, the liquid separation increased (region B) and then decreased (region C) as the mixture proceeds downwards. The validity of the existing models for branching flows at parallel T-junction was tested, and turned out to be appropriate for region A. However, the models were not applicable to the rear part due to a strong flow recirculation. Moreover, the effect of the membranes in channels was investigated, but that was minor.     

Model of particle deposition from turbulent gas-solid flow in channels with absorbent walls

A nonlocal model of particle deposition is developed without resorting to empirical information on the fluctuating motion of the particles. The effects of particle inertia are described by a system of differential equations for the moments of the dispersed phase velocity. The model is tested on examples of flows in channels with smooth walls and with grassy roughness.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 58–65, January–February, 1992.     

Low Reynolds number flow inside straight micro channels with irregular cross sections

This paper presents an analysis of the compound effect of finite temperature differences and fluid friction on the existence of an optimum laminar flow regime in singly connected micro channels with complex free flow area cross sections. A widespread conviction has been established that the two competing irreversibility sources in a channel flow with heat transfer lead to the existence of an optimum flow regime. The results presented in this paper clearly shows the opposite. When an objective function is represented by the entropy generation rate per unit heat capacity rate of the fluid stream, the thermodynamic optimum flow regime represents a rather rare occurrence in the laminar region of irregularly shaped ducts. The presence of an extremum is more probable for very small diameters, the ones of an order of magnitude of O(≤10⁻³ m). The analysis is performed for selected ranges of relevant geometric, flow, and thermal parameters of a set of straight micro channels with irregular free flow area cross-sections. The following geometries of the free flow area cross section were investigated: (i) sine duct, (ii) circular duct, (iii) elliptical duct, (iv) moon-shaped ducts, and (v) four-cuspped duct. The range of Reynolds numbers has been established between O(10²) and O(10⁴). The existence of the objective function minimum is confirmed for ducts with an irregular cross section only for very small hydraulic diameters. These minima are relatively weak, and as a general rule, the sets of optimum parameters are close to the onset of turbulence or possibly even in the transitional or turbulent regions. Received on 10 November 1998     

Experiments on two-phase flow distribution inside parallel channels of compact heat exchangers

Uneven distribution in heat exchangers is a cause of reduction in both thermal and fluid-dynamic performances. Many papers have dealt with single-phase flow and both flow distribution data and analytical or numerical models are available for header design. With regard to two-phase flow, phase separation in manifolds with several outlets is so complicated that, to date, there is no general way to predict the distribution of two-phase mixtures at header-channel junctions. The design of headers for new generation compact heat exchangers and multi-microchannel evaporators is still based on an empirical approach, as a number of variables act together: geometrical parameters and orientation of the manifolds and of the channels, operating conditions, fluid physical properties.     

The influence of rotation on turbulent flow and heat transfer in an annulus between independently rotating tubes

Experimental and numerical investigations of turbulent flow and heat transfer have been performed in a concentric annulus between independently rotating tubes. Numerical predictions, applying a Reynolds stress turbulence model, are compared with experimental fluid flow and heat transfer results for the case of a heated outer tube and an adiabatic inner tube. Compared to the above mentioned boundary conditions for the conservation equation of energy, differences in heat transfer in case of a heated inner tube and an adiabatic outer one, are examined by analysis, applying a mixing length turbulence model. Numerical investigations with both kinds of models about the influence of annulus radius ratio make evident that due to different superimpositions of centrifugal force and additional shear stress there is a wide variation of effects on fluid flow and heat transfer caused by the rotation of the inner and the outer tube.     

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%.     

Forced convective film condensation inside vertical tubes

A theoretical study of forced convective film condensation inside vertical tubes is presented. We propose a unified procedure for predicting the pressure gradient and condensation heat transfer coefficient of a vapor flowing turbulently in the core and associated with laminar or turbulent film on the tube wall. The analysis for the vapor flows is performed under the condition that the velocity profiles are locally self-similar. The laminar and turbulent film models equate the gravity, pressure and viscous forces, and consider the effect of interfacial shear. The transition from laminar to turbulent film depends not only on the liquid Reynolds number but also on the interfacial shear stress. In this work we also proposed a new eddy viscosity model which is divided into three regions: the inner region in liquid condensate near the wall; the interface region including both liquid and vapor; and the outer region for the vapor core. Comparisons of the theory with some published experimental data showed good agreement.