The instability of rhombic cell flows

Instability of two-dimensional periodic flows with rhombic cell structure represented by the stream function Ψ=cos kx+cosy is investigated. Stability characteristics are obtained for the Reynolds number R=1, 2, 3 and 4 and the ratio of the diagonals of the cell . Variation of the critical Reynolds number R_c with k is obtained, and the square cell flow (k=1) is found to be most stable (R_c=√2). It is found that R_c → 1 as k → 0, which leads to a finite gap between this limiting R_c and R_c=√2 for K=0 (Ψ=cos y).     

Experimental investigation of heat transfer and pressure drop characteristics of flow through spirally fluted tubes

Spirally fluted tubes are used extensively in the design of tubular heat exchangers. In previous investigations, results for tubes with flute depths e/D_vi < 0.2 were reported, with most correlations applicable for Re ≥ 5000. This paper presents the results of an experimental investigation of the heat transfer and pressure drop characteristics of spirally fluted tubes with the following tube and flow parameter ranges: flute depth e/D_vi = 0.1−0.4, flute pitch p/D_vi = 0.4−7.3, helix angle θ/90° = 0.3−0.65, Re = 500−80,000, and Pr = 2−7. The heat transfer coefficients inside the fluted tube were obtained from measured values of the overall heat transfer coefficient using a nonlinear regression scheme. The friction factor data obtained consisted of 507 data points. The proposed correlation for the friction factor predicts 96% of the database within ±20%. The heat transfer correlation for the range 500 ≤ Re ≤ 5000 predicts 76% of the database (178 data points) within ±20%, and the correlation for the higher Re range predicts 97% of the 342 data points within ±20%. Comparison of heat transfer and friction data show that these tubes are most effective in the laminar and transition flow regimes. The present results show that the increase of flute depth in the range considered does not improve heat transfer.     

Three-dimensional numerical study and field synergy principle analysis of wavy fin heat exchangers with elliptic tubes

Three dimensional numerical studies were performed for laminar heat transfer and fluid flow characteristics of wavy fin heat exchangers with elliptic/circular tubes by body-fitted coordinates system. The simulation results of circular tube were compared with the experiment data, then circular and elliptic (e = b/a = 0.6) arrangements with the same minimum flow cross-sectional area were compared. A max relative heat transfer gain of up to 30% is observed in the elliptic arrangement, and corresponding friction factor only increased by about 10%. The effects of five factors on wavy fin and elliptic tube heat exchangers were examined: Reynolds number (based on the smaller ellipse axis, 500  4000), eccentricity (b/a, 0.6  1.0), fin pitch (F_p/2b, 0.05  0.4), fin thickness (F_t/2b, 0.006  0.04) and tube spanwise pitch (S₁/2b, 1.0  2.0). The results show that with the increasing of Reynolds number and fin thickness, decreasing of the eccentricity and spanwise tube pitch, the heat transfer of the finned tube bank are enhanced with some penalty in pressure drop. There is an optimum fin pitch (F_p/2b = 0.1) for heat transfer, but friction factor always decreases with increase of fin pitch. And when F_p/2b is larger than 0.25, it has little effects on heat transfer and pressure drop. The results were also analyzed from the view point of field synergy principle. It was found that the effects of the five factors on the heat transfer performance can be well described by the field synergy principle.     

A study of properties of vortex stretching and enstrophy generation in numerical and laboratory turbulence

Dynamically relevant alignments are used in order to show that regions with weak vorticity are not structureless, non-Gaussian and dynamically not passive. for example, the structure of vorticity in quasi-homogeneous/isotropic turbulent flows is associated with strong alignment between vorticity ω and the eigenvectors of the rate of strain tensor λ_i (especially — but not only — between ω and λ₂) rather than with intense vorticity only. Consequently, much larger regions of turbulent flow than just those with intense vorticity are spatially structured. The whole flow field — even with the weakest measurable enstrophy — is strongly non-Gaussian, which among other things is manifested in strong alignment between vorticity and the vortex stretching vector W_i ≡ ω_jS_ij. It is shown that the quasi-two-dimensional regions corresponding to large cos(ω, λ₂) are qualitatively different from purely two-dimensional ones, e.g. in that they possess essentially nonvanishing enstrophy generation, which is larger than its mean for the whole field.     

A correlation of the isentropic exponents of real gases

In previous publications, three isentropic exponents, k_pv, k_Tv, k_pT, have been introduced, which when used in place of the classical isentropic exponent k = c_p/c_v in the ideal gas isentropic change equations, the latter can describe very accurately the isentropic change of real gases. The present work provides a general method for determining the values of k_pv, k_Tv, k_pT within the ranges of reduced pressure p_r = 0 to 10 and of reduced temperature T_r = 1 to 3.5, thus allowing the calculation of the isentropic flow of those real gases for which no detailed thermodynamic data are available. The accuracy obtained is limited only by the accuracy of the generalized Lee-Kesler equation of state, which is employed in the method developed.     

多层平行裂隙型多孔介质通道内流体流动传热特性

基于Brinkman-extended Darcy模型和局部热平衡模型,对多层平行裂隙型多孔介质通道内的流动传热特性进行研究.获得了多层平行裂隙型多孔介质通道内各区域的速度场、温度场、摩擦系数及努塞尔数解析解,并分析了裂隙层数、达西数、空心率、有效热导率之比等对通道内流动传热特性的影响.结果表明:达西数较小时,通道多孔介质层内会出现不随高度变化的达西速度,此达西速度会随裂隙层数的增加而增大,但却不受各裂隙层下多孔介质层位置变化的影响.增加裂隙层数会减弱空心率对压降的影响,会使通道内流体压降升高,但升高程度会逐渐降低.增大热导率之比或减小空心率会使多裂隙通道内出现阶梯式温度分布,而在较小热导率之比或较大空心率时多裂隙情况下的温度分布曲线会趋于一致.此外,当热导率之比较小时,多层裂隙通道内的传热效果在任何空心率下都要优于单裂隙情况,当热导率之比较大时,存在临界空心率使各裂隙层数通道内的传热效果相同,且多裂隙通道内继续增加裂隙层数对传热强度影响不大.     

Hole diameter effect on flow characteristics of wake behind porous fences having the same porosity

The hole diameter effect on the flow characteristics of wake behind porous fences has been investigated experimentally in a circulating water channel having a test section of 300^w×200^h×1200^l (mm). Three porous fences having different hole diameters of d=1.4,2.1,2.8 mm were tested in this study, but they have the same =38.5% geometric porosity. One thousand instantaneous velocity fields for each fence were measured consecutively by the hybrid PTV system employing a high-speed CCD camera. Free stream velocity was fixed at 10 cm/sec and the corresponding Reynolds number based on the fence height was Re=2,985. Consequently, the fence with the smallest hole diameter d=1.4 mm (d_1.4) decreases the streamwise velocity component and increases the vertical velocity component. Among the three hole diameters tested in this study, the d_1.4 fence has the largest turbulence intensity in the shear layer developed from the fence top. Regardless of the hole diameter, however, all three fences having the same porosity reduce the reduction of turbulent intensity in the lower region below the fence height (y/H<1).     

Impingement heat transfer at a circular cylinder due to a submerged slot jet of water

An experimental investigation was carried out on the heat transfer due to a submerged slot jet of water impinging on a circular cylinder in crossflow. The cylinder diameter and the slot width are of the same order of magnitude, specifically D_s = 2.0 and 3.0 mm and D_c = 2.5 and 3.0 mm. The experimental apparatus allowed variation of the slot width, the cylinder diameter, and the distance from nozxle exit to heater. Conditions of impingement from the bottom (ascending flow) were taken into consideration as well as impingement from above (descending flow). The Nusselt number was determined as a function of Reynolds and Prandtl numbers in the range 1.5 × 10³ < Re < 2.0 × 10⁴, 2.7 < Pr < 7.0, and 1.5 ≤ z/D_s ≤ 10. The experimental data were correlated with a simple equation that fits 90% of the data with a precision of 20%.     

Conjugate natural convection from an array of discrete heat sources: part 2 — a numerical parametric study

Coupled conduction and natural convection transport within a discretely heated cavity have been investigated numerically. One vertical wall of the cavity is composed of discrete, isoflux heat sources mounted in a substrate of finite thermal conductivity. The opposite vertical wall and the horizontal walls are assumed to be isothermal and adiabatic, respectively. The governing steady-state partial differential equations for the fluid and solid region are solved simultaneously using a control volume formulation, coupled with an additive correction multigrid procedure that increases the convergence rate of the solution. The fluid Prandtl number and heater/fluid thermal conductivity ratio are fixed at 25 and 2350, respectively, corresponding to a dielectric fluid (FC-77) and heaters manufactured from silicon. With increasing modified Rayleigh number (10⁴ < Ra_Lz^* < 109), the cavity flow becomes more boundary layer-like along the vertical walls, and multiple fluid cells develop in the central region. Thermal spreading in the substrate increases with decreasing modified Rayleigh number and with increasing values of the substrate/fluid thermal conductivity ratio (10⁻¹ <- R_s ≤ 10³). For large R_s, the discrete heat sources lose their thermal identity, and the streamlines and isotherms resemble those associated with a differentially heated cavity. Thermal spreading in the substrate also has a significant effect on circulation in the cavity and on maximum surface temperatures.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

Numerical study of steady/unsteady flow and heat transfer in porous media using a characteristics-based matrix-free implicit FV method on unstructured grids

In this study, a high-resolution characteristic-based finite-volume (FV) method on unstructured grids [Int. J. Numer. Method Eng. 50 (2001) 11; Int. J. Heat Fluid Flow 21 (2000) 432] is extended by a matrix-free implicit dual-time stepping scheme for the numerical simulation of steady and unsteady flow and heat transfer with porous media. The method has been used to study the characteristics of a complex problem: flow and heat transfer in a channel with multiple discrete porous blocks, which was originally proposed by Huang and Vafai [J. Thermophys. Heat Transfer 8 (3) (1994) 563]. In addition, flow and heat transfer in a channel partially or fully filled with porous layers and containing solid protruding blocks with constant heat flux on its lower surface are also investigated in details. Hydrodynamic and heat transfer results are reported for both steady and transient flow cases. In particular, the effects of Darcy and Reynolds numbers on heat transfer augmentation and pressure loss are studied. An in-depth discussion of the formation and variation of recirculation is presented and the existence of optimum porous insert is demonstrated. At high Reynolds numbers the flow in the porous channel exhibits a cyclic characteristics although unlike the non-porous channel flow, the cyclic vortex development is only restricted to a small area behind the last solid block, while temperature changes more slowly and does not exhibit cyclic variations over a long period of time. It is shown that for all the cases studied altering some parametric values can have significant and interesting effects on both flow pattern as well as heat transfer characteristics.     

Fluid flow in a rotating cylindrical container with a rotating disk at the fluid surface

Fluid flow in a rotating cylindrical container of radius R_w and height H with a co-axially rotating disk of radius R_d at the fluid surface is numerically investigated. The container and the disk rotate with angular velocities Ω_w and Ω_d, respectively. We solve the axisymmetric Navier-Stokes equations using a finite-volume method. The effects of the relative directions and magnitudes of the disk and container rotations are studied. The calculations are carried out with various ratios of Ω_w and Ω_d for H/R_w = 2 and R_d/R_w = 0.7. Streamlines and velocity vectors in the meridional plane and azimuthal velocities are obtained. The flow fields in the meridional plane are discussed with relation to azimuthal velocities in the interior of the container. The numerical results are also compared with experimental data.     

Experimental study of flow and heat transfer in rib-roughened rectangular channels

Secondary flow patterns, pressure drop and heat transfer in rib-roughened rectangular channels have been investigated. The aspect ratio of the channels is 1–8, and ribs are attached to the wide channel walls in order to set up swirling motions. The geometries tested consist of channels having cross ribs, parallel ribs, cross V-ribs, parallel V-ribs, and multiple V-ribs (Swirl Flow Tube). The flow patterns were investigated using smoke wire visualization and LDV measurements. The smoke wire experiments have been performed at Re=1100 and the LDV measurements at Re=3000 at periodic fully developed conditions. The heat transfer and pressure drop are described by j and f factors for Reynolds numbers from 500 to 15 000. The distributions of axial mean velocity and turbulent fluctuations are strongly influenced by the secondary flows. Large mean velocities and small fluctuations are found in regions where the secondary flow is directed towards a surface, while small mean velocities and large fluctuations are found in regions where the secondary flow is directed away from a surface. The Swirl Flow Tube provides a significant increase in the j factor at Reynolds numbers from 1000 to 2000, but unfortunately also an increase in the f factor. At higher Reynolds numbers, the j and f factors of the Swirl Flow Tube are of the same order of magnitude as for the other rib-roughened channels. It is found that the flow direction in a channel with parallel V-ribs has important influence on the j/f ratio. At Reynolds numbers above 4000, this channel provides the highest j/f ratio if the V-ribs are pointing upstream; while it provides the lowest j/f ratio of all rib configurations, if the V-ribs are pointing downstream.     

Steady laminar flow through thin curved pipes

Let R, τ denote, respectively, the radius of curvature and radius of torsion of the pipe (centre-line) and let a be a typical cross-sectional diameter.

The major part of the present paper addresses the case of flows through pipes of constant cross-section; (Re)²(a/R), Re(a/τ), (a/R) and (a/τ) all being small. Re is the Reynolds number for the flow. It is found that, even without further specifications of the details of the pipe, many important results can be obtained about the secondary flow which occurs and the pressure losses resulting from it. For example, it is shown that an important feature of such flows is valid for any corss-sectional shape; this was not obvious from previous works which treated only special cases having significant symmetries. Also, a new method for calculating the modified axial pressure gradient is presented which reduces dramatically the amount of work required therefor.

The remainder of the paper presents some results for similar flows through pipes of varying cross-section.     

A numerical method for forced convection in porous and homogeneous fluid domains coupled at interface by stress jump

A numerical method was developed for flows involving an interface between a homogeneous fluid and a porous medium. It is based on the finite volume method with body‐fitted and multi‐block grids. The Brinkman–Forcheimmer extended model was used to govern the flow in the porous medium region. At its interface, the flow boundary condition imposed is a shear stress jump, which includes the inertial effect, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The forced convection through a porous insert over a backward‐facing step is investigated. The results are presented with flow configurations for different Darcy numbers, 10^?2 to 10^?5, porosity from 0.2 to 0.8, Reynolds number from 10 to 800, and the ratio of insert length to channel height from 0.1 to 0.3. The heat transfer is improved by using porous insert. To enhance the heat transfer with minimal frictional losses, it is preferable to have a medium length of insert with medium Darcy number, and larger Reynolds number. The interfacial stress jump coefficients β and β₁ were varied from ?1 to 1, and within this range the average and local lower‐wall Nusselt numbers are not sensitive to the parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

Vortex shedding behind rings and discs

The wake structure of discs and bluff rings has been investigated experimentally in a wind tunnel. The rings have an inner diameter d_i, and an outer diameter d_o and are classified according to the parameter (d_o + d_i)/(d_o − d_i) = d/w. the ratio of mean diameter to ring width. As d/w → ∞ the flow approaches that around a two dimensional bluff body whereas as d/w tends to unity the body approaches a solid disc. A distinct change in the vortex shedding pattern is found around d/w = 5. Below this critical value velocity fluctuations in the wake have a weak periodic component which is 180° out of phase across a diameter of the body. Above d/w = 5. regular and coherent axisymmetric vortex ring shedding is observed with shedding occurring alternately from the inner and outer circumferences of the bluff body. Flow visualization and conditional averaging of hot-wire data are used to investigate the vortex structure.     

Analysis of Thermal Flow in a Rotating Porous U-Turn Duct Using Lattice Boltzmann Method

The lattice Boltzmann method is carried out to investigate the heat transfer enhancement in a U-turn duct which is partially filled with a porous media. The porous layer is inserted at the core of the duct and is modeled using the Brinkman–Forchheimer assumptions. In order to validate the results, first a channel flow problem without any porous layer is compared with available data. Second, the porous Couette flow and partially porous channel flow are successfully compared with the studies of other researchers. Then, fluid flow in a clear U-turn duct is studied looking carefully at the velocity, curvature and rotation effects. Finally, the effects of porous layer thickness on the rate of heat transfer and pressure drop are investigated. Parametric studies are conducted to evaluate the effects of various parameters (i.e., Reynolds number, Darcy number, rotation number), highlighting their influences on the thermo-hydrodynamics behavior of the flow. The optimum values of porous layer thickness are presented for specific flow parameters.     

Inertial instability of the Stewartson -layer

Inertial stability of a vertical shear layer (Stewartson E^1/4-layer) on the sidewall of a cylindrical tank with respect to stationary axisymmetric perturbations is inverstigated by means of a linear theory. The stability is determined by two non-dimensional parameters, the Rossby number Ro = U/2ΩL and Ekman number E = v/ΩH², where U and L = (E/4)1/4H are the characteristic velocity and width of the shear layer, respectively, Ω the angular velocity of the basic rotation, v the kinematic viscosity and H the depth of the tank.

For a given Ekman number, the flow is more unstable for larger values of the Rossby number. For E = 10⁻⁴, which is a typical value of the Ekman number realized in rotating tank experiments, the critical Rossby number Ro_c for instability and the critical axial wavenumber m_c non-dimensionalized by L⁻¹ are found to be 1.3670 and 8.97, respectively. The value of Ro_c increases and that of m_c decreases with increasing E.     

DNS of heat transfer in turbulent and transitional channel flow obstructed by rectangular prisms

Direct numerical simulation (DNS) of heat transfer in a channel flow obstructed by rectangular prisms has been performed for Re_τ = 80–20, where Re_τ is based on the friction velocity, the channel half width and the kinematic viscosity. The molecular Prandtl number is set to be 0.71. The flow remains unsteady down to Re_τ = 40 owing to the disturbance induced by the prism. For Re_τ = 30 and 20, the flow results in a steady laminar flow. In the vicinity of the prism, the three-dimensional complex vortices are generated and heat transfer is enhanced. The Reynolds number effect on the time-averaged vortex structure and the local Nusselt number are investigated. The mechanism of the heat transfer enhancement is discussed. In addition, the mean flow parameters such as the friction factor and the Nusselt number are examined in comparison with existing DNS and experimental data.     

Effects of Brownian Motion and Thermophoresis on the Mixed Convection of Nanofluid in a Porous Channel Including Flow Reversal

This article is devoted to combined convection heat transfer of nanofluids through a vertical channel filled with a homogeneous and isotropic porous medium. The flow is assumed to be fully developed and the “Brinkman extended Darcy” model is used for the flow in the porous media and “clear compatible” viscous dissipation model is considered. Also the model utilized for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The governing momentum, energy, and nanopartices volume fraction equations are solved both analytically and numerically. The effects of the influential dimensionless parameters such as Brownian and thermophoresis parameters, mixed convection parameter (Gr/Re), Brinkman, Darcy and Lewis numbers on dimensionless velocity and temperature distributions and pressure drop are studied. Also, the results of the Nusselt number for the both left and right walls are presented and discussed.