Lattice Boltzmann simulation of flow past a circular cylinder near a moving wall

The two‐dimensional flows past a circular cylinder near a moving wall are simulated by lattice Boltzmann method. The wall moves at the inlet velocity and the Reynolds number ranges from 300 to 500. The influence of the moving wall on the flow patterns is demonstrated and the corresponding mechanism is illustrated by using instability theory. The correlations among flow features based on gap ratio are interpreted. Force coefficients, velocity profile and vortex structure are analyzed to determine the critical gap ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

二维水平通道内流动沸腾换热的格子Boltzmann模拟

利用格子Boltzmann方法模拟二维水平通道内水的流动沸腾过程,获得不同壁面过热度下流型特点和不同因素对换热过程的影响规律。结果表明,随着壁面过热度升高,流道内流型依次经历从泡状流、弹状流到反环流的转变,平均热流密度和平均换热系数先增大后减小。入口流速降低会使流道内出现受限气泡流,核态沸腾受到抑制。提高入口流速能够有效促进气泡脱离,壁面平均换热系数随入口流速增大而增大,但增长速率有所减小。减小通道宽度有利于汽化现象发生,核态沸腾得到强化,壁面平均换热系数有所提高。     

Lattice Boltzmann simulation of heat transfer and fluid flow in a microchannel with nanofluids

Mathematical modeling is performed to simulate forced convection flow of 47 nm- Al₂O₃/water nanofluids in a microchannel using the lattice Boltzmann method (LBM). Single channel flow and conjugate heat transfer problem are taken into consideration and the heat transfer rate using a nanofluid is examined. Simulations are conducted at low Reynolds numbers (2 ≤ Re ≤ 16). The computed average Nusselt number, which is associated with the thermal conductivity of nanofluid, is in the range of 0.6 ￡ [`(Nu)] ￡ 13 0.6 \le \overline{Nu} \le 13 . Results indicate that the average Nusselt number increases with the increase of Reynolds number and particle volume concentration. The fluid temperature distribution is more uniform with the use of nanofluid than that of pure water. Furthermore, great deviations of computed Nusselt numbers using different models associated with the physical properties of a nanofluid are revealed. The results of LBM agree well with the classical CFD method for predictions of flow and heat transfer in a single channel and a microchannel heat sink concerning the conjugate heat transfer problem, and consequently LBM is robust and promising for practical applications.     

Lattice Boltzmann simulation of phase change and heat transfer characteristics in the multi-layer deposition

The metal droplets deposition method(MDDM) is a rapid prototyping technology, implemented via metallurgy bonding within droplets. The anisotropy of heat transfer and re-melting is caused by an asymmetric deposition process. A lattice Boltzmann method(LBM) model is established to predict the heat transfer and phase change in the multi-layer deposition. The prediction model is verified by the experimental temperature profiles in existing literature. The monitoring points are set to compare the tem...     

Analytical study of heat transfer from circular cylinder in liquid metals

In this study the influence of a thin hydrodynamic boundary layer on the heat transfer from a single circular cylinder in liquid metals having low Prandtl number (0.004–0.03) is investigated under isothermal and isoflux boundary conditions. Two separate analytical heat transfer models, viscous and inviscid, are developed to clarify the discrepancy between previous results. For both models, integral approach of the boundary layer analysis is employed to derive closed form expressions for the calculation of the average heat transfer coefficients. For an inviscid model, the energy equation is solved using potential flow velocity only whereas for a viscous model, a fourth-order velocity profile is used in the hydrodynamic boundary layer and potential flow velocity is used outside the boundary layer. The third-order temperature profile is used inside the thermal boundary layer for both models. It is shown that the inviscid model gives higher heat transfer coefficients whereas viscous flow model gives heat transfer results in a fairly good agreement with the previous experimental/numerical results.     

Numerical simulation of natural convective heat transfer and fluid flow around a heated cylinder inside an enclosure

 A collocated, non-orthogonal grid based finite volume technique has been applied for investigating the two dimensional natural convective flow and heat transfer around a heated cylinder kept in a square enclosure. The effects of different enclosure wall thermal boundary conditions, fluid Prandtl number and the ratio between enclosure and cylinder dimensions (aspect ratio) upon the flow and thermal features, have been systematically studied. It is observed that the patterns of recirculatory flow and thermal stratification in the fluid are significantly modified, if any of these parameters is varied. The overall heat transfer rates are also affected due to the changes in the flow and temperature patterns. The study presents useful observations regarding the variation of local Nusselt number along each wall, for the different cases considered. Received on 2 August 2000 / Published online: 29 November 2001     

Momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel

Unsteady momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel is numerically investigated using FLUENT for the ranges of Reynolds numbers as 10≤Re≤500, of the blockage ratio as 0.1≤β≤0.4, and of the gap ratio as 0.125≤γ≤1 for a constant value of the Prandtl number of 0.744. The transition of the flow from steady to unsteady (characterized by critical Re) is determined as a function of γ and β. The effect of γ on the mean drag and lift coefficients, Strouhal number (St), and Nusselt number (Nu _w ) is studied. Critical Re was found to increase with decreasing γ for all values of β. and St were found to increase with decreasing values of γ for fixed β and Re. The effect of decrease in γ on was found to be negligible for all blockage ratios investigated.     

Lattice Boltzmann simulation of the effects of cavity structures and heater thermal conductivity on nucleate boiling heat transfer

The boiling heat transfer technology with cavity surfaces can provide higher heat flux under lower wall superheat, which is of great significance for the cooling of electronic chips and microelectromechanical devices. In this paper, the boiling characteristics of the cavity surfaces are investigated based on the lattice Boltzmann(LB) method, focusing on the effects of cavity shapes, sizes, and heater thermal conductivity on the heat transfer performance. The results show that the triangular cavi...     

Forced convective mass transfer from circular disks

Overall forced convective mass transfer rates were measured from naphthalene circular disks of two inches diameter suspended in an air stream. The data consisting of fifty two experiments were correlated with following expression with a standard deviation of 2.0%, Sh=0.266 (Re)^0.60 (Sc)^1/3 for the range of 270< Re< 34900.-Data obtained for mass transfer from frontal areas of disks were correlated with a standard deviation of 2.63% by following relationship. Sh=0.67 (Re)_0.54 (Re)^1/3 for the range of 250< Re< 34800.

---

Zusammenfassung Die Gesamtbeträge der Stoffübertragung bei Zwangskonvektion wurden an einer Naphthalinscheibe von 50 mm Durchmesser gemessen, die in einem Luftstrom aufgehängt war. Die in 52 Versuchen gewonnenen Meßdaten wurden durch den folgenden Ausdruck mit einer Standardabweichung von 2% korreliert. Sh=0,266 (Re)^0,60(Sc)^1/3, gültig für 270< Re< 34 900. — Die Versuchsergebnisse für den Stoffübergang von den Stirnseiten der Scheibe konnten durch folgenden Ausdruck Sh=0,67(Re)^0,54(Re)^1/3, gültig für 250 < Re < 34 800 mit einer Standardabweichung von 2,63%, korreliert werden.

---

Nomenclature A Constant in Eq. (3) - B Constant in Eq. (3) - k_c Continuous phase mass transfer coefficient (cm/hr) - N Overall mass transfer rates (gm · mole/hr) - P^* Equilibrium vapour pressure of naphthalene (mm · Hg) - p₀ Vapour pressure of naphthalene in air stream of wind tunnel (mm · Hg) - q Number of observations - R Gas constant, 62.368 x 10³ (cm³) (mm · Hg)/(°k) (gm · moles) - Re =ud/ Reynolds number - S Surface area (cm₂) - S² (Y) Estimate of variance, bases on the deviations of data points from their predicted values - Sc =/D_v Schmidt number - Sh=k_cd/D_v Sherwood number - T Absolute temperature (°K) - U Velocity of air in the wind tunnel. (cm/sec.) - ² True variance calculated with infinite degree of freedom     

PIV measurements of the wake behind a rotationally oscillating circular cylinder

The near-wake behind a circular cylinder undergoing rotational oscillatory motion with a relatively high forcing frequency has been investigated experimentally. Experiments were carried out varying the ratio of the forcing frequency f_f to the natural vortex shedding frequency f_n in the range of 0.0 (stationary) to 1.6 at an oscillation amplitude of θ_A=30° and Reynolds number of Re=4.14×10³. Depending on the frequency ratio (F_R=f_f /f_n), the near-wake flow could be divided into three regimes—non-lock-on (F_R=0.4), transition (F_R=0.8, 1.6) and lock-on (F_R=1.0) regimes—with markedly different flow structures. When the frequency ratio was less than 1.0 (F_R⩽1.0), the rotational oscillatory motion of the cylinder decreased the length of the vortex formation region and enhanced the mutual interaction between large-scale vortices across the wake centerline. The entrainment of ambient fluid seemed to play an important role in controlling the near-wake flow and shear-layer instability. In addition, strong vortex motion was observed throughout the near-wake region. The flow characteristics changed markedly beyond the lock-on flow regime (F_R=1.0) due to the high frequency forcing. At F_R=1.6, the high frequency forcing decreased the size of the large-scale vortices by suppressing the lateral extent of the wake. In addition, the interactions between the vortices shed from both sides of the cylinder were not so strong at this forcing frequency. As a consequence, the flow entrainment and momentum transfer into the wake center region were reduced. The turbulent kinetic energy was large in the region near the edge of the recirculation region, where the vortices shed from both sides of the cylinder cross the wake centerline for all frequency ratios except for the case of F_R=1.6. The temporally resolved quantitative flow information extracted in the present work is useful for understanding the effects of open-loop active flow control on the near-wake flow structure.     

Visual observations of the flow around a half-submerged oscillating circular cylinder

Visual observations are made on the flow around a horizontal circular cylinder which is half-submerged in still water and forced to oscillate vertically. The ends of the cylinder have great influence on the wave pattern and flow field. Progressive plane waves are generated at small forcing amplitudes, but cross-waves are superimposed on the progressive plane waves at large forcing amplitudes. The wavelength of the cross-waves in the direction parallel to the cylinder axis increases with the forcing amplitude. The crests of the cross-waves are in parallel lines which are oblique to the cylinder axis. The angle at which the parallel lines meet the cylinder axis decreases as the forcing amplitude is increased. Three kinds of steady flows are induced in the water: surface flow, undersurface flow, and vertical jet.     

Exact solution of the problem of convective heat exchange for a circular cylinder at small Prandtl numbers

The problem of the heat exchange of a circular cylinder in an incompressible flow at small Prandtl numbers Pr ≪ 1 is solved. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 43–48, January–February, 1994.     

Flow structure of wake behind a rotationally oscillating circular cylinder

Flow around a circular cylinder oscillating rotationally with a relatively high forcing frequency has been investigated experimentally. The dominant parameters affecting this experiment are the Reynolds number (Re), oscillation amplitude (θ_A), and frequency ratio F_R=f_f/f_n, where f_f is the forcing frequency and f_n is the natural frequency of vortex shedding. Experiments were carried out under conditions of Re=4.14×10³, 0°θ_A60° and 0.0F_R2.0. Rotational oscillation of the cylinder significantly modified the flow structure in the near-wake. Depending on the frequency ratio F_R, the cylinder wake showed five different flow regimes, each with a distinct wake structure. The vortex formation length and the vortex shedding frequency were greatly changed before and after the lock-on regime where vortices shed at the same frequency as the forcing frequency. The lock-on phenomenon always occurred at F_R=1.0 and the frequency range of the lock-on regime expanded with increasing oscillation amplitude θ_A. In addition, the drag coefficient was reduced when the frequency ratio F_R was less than 1.0 (F_R<1.0) while fixing the oscillation amplitude at θ_A=30°. When the oscillation amplitude θ_A was used as a control parameter at a fixed frequency ratio F_R=1.0 (lock-on regime), the drag reduction effect was observed at all oscillation amplitudes except for the case of θ_A=30°. This type of active flow control method can be used effectively in aerodynamic applications while optimizing the forcing parameters.     

A turbulence model for the heat transfer near stagnation point of a circular cylinder

A one-equation low-Reynolds number turbulence model has been applied successfully to the flow and heat transfer over a circular cylinder in turbulent cross flow. The turbulence length-scale was found to be equal 3.7y up to a distance 0.05 and then constant equal to 0.185 up to the edge of the boundary layer (wherey is the distance from the surface and is the boundary layer thickness).The model predictions for heat transfer coefficient, skin friction factor, velocity and kinetic energy profiles were in good agreement with the data. The model was applied for Re 250,000 and Tu0.07.Nomenclature µ,C _D Constants in the turbulence kinetic energy equation - C ₁,C ₂ Constants in the turbulence length-scale equation - Skin friction coefficient atx - D Cylinder diameter - F Dimensionless flow streamwise velocityu/u _e - k Turbulence kinetic energy =1/2 the sum of the squared three fluctuating velocities - K Dimensionless turbulence kinetic energyk/u _e ^/2 - I Dimensionless temperature (T–T _w )/(T –T _w ) - l Turbulence length-scale - l _e Turbulence length-scale at outer region - Nu _D Nusselt number - p Pressure - Pr Prandtl number - Pr _t Turbulent Prandtl number - Pr _k Constant in the turbulence kinetic energy equation - R Cylinder radius - Re _D Reynolds number u D/µ - Re _x Reynolds number u x/µ - R _K Reynolds number of turbulence - T Mean temperature - T Mean temperature at ambient - T _s Mean temperature at surface - Tu Cross flow turbulence intensity, - u Mean flow streamwise velocity - u Fluctuating streamwise velocity - u _e Mean flow velocity at far field distance - u Mean flow velocity at ambient - u* Friction velocity - v Mean velocity normal to surface - V Dimensionless mean velocity normal to surface - x,x ₁ Distance along the surface - y Distance normal to surface - Dimensionless pressure gradient parameter - Boundary layer thickness atu=0.9995u _e - Transformed coordinate iny direction - Fluid molecular viscosity - _t Turbulent viscosity - _eff + _t - µ Fluid molecular viscosity at ambient - Kinematic viscosity/ - Density - Density at ambient - _w Wall shear stress - _w,0 Wall shear stress at zero free stream turbulence     

Unsteady natural convective flow past a moving vertical cylinder with heat and mass transfer

Transient natural convection boundary layer flow of an incompressible viscous fluid past an impulsively moving semi- infinite vertical cylinder is considered. The temperature and concentration of the cylinder surface are taken to be uniform. The unsteady, nonlinear and coupled governing equations of the flow are solved using an implicit finite difference scheme. The finite difference scheme is unconditionally stable and accurate. Numerical results are presented with various sets of parameters for both air and water. Transient effects of velocity, temperature and concentration profiles are analyzed. Local and average skin friction, rates of heat and mass transfer are shown graphically. Received on 1 November 1999     

A numerical study of the steady forced convection heat transfer from an unconfined circular cylinder

Forced convection heat transfer from an unconfined circular cylinder in the steady cross-flow regime has been studied using a finite volume method (FVM) implemented on a Cartesian grid system in the range as 10 ≤ Re ≤ 45 and 0.7 ≤ Pr ≤ 400. The numerical results are used to develop simple correlations for Nusselt number as a function of the pertinent dimensionless variables. In addition to average Nusselt number, the effects of Re, Pr and thermal boundary conditions on the temperature field near the cylinder and on the local Nusselt number distributions have also been presented to provide further physical insights into the nature of the flow. The rate of heat transfer increases with an increase in the Reynolds and/or Prandtl numbers. The uniform heat flux condition always shows higher value of heat transfer coefficient than the constant wall temperature at the surface of the cylinder for the same Reynolds and Prandtl numbers. The maximum difference between the two values is around 15–20%.     

Formation of vortex street and vortex pair from a circular cylinder oscillating in water

The flow around a circular cylinder undergoing sinusoidal oscillating movement in still water is investigated by phase-locked PIV measurements. The pattern and development of large-scale vortex structures in the flow are studied from the velocity vectors and vorticity contours obtained at eight successive phases of an oscillating cycle. Experiments are performed at three Keulegan–Carpenter numbers; KC=12, 6.28 and 4.25. Results at KC=12 reveal the mechanism of vortex formation and the development of the shed vortices into a vortex street at a lateral direction to the line of cylinder movement. The role of a biased flow stream and the length of the cylinder stroke in the formation of the vortex street are discussed. At the lower KC numbers, a symmetric pair of vortices is found attached to the leeward face of the cylinder. The vortex pair exhibits an increasing degree of asymmetry when KC increases from 4.25 to 6.28. An explanation in terms of the length of the cylinder strokes and the degree of flow asymmetry is offered for the transition of flow regimes from a vortex pair to a vortex street. The present results are compared with the observations made in previous experimental and numerical studies in the literature.     

Three-dimensional instability of an oscillating viscous flow past a circular cylinder

IntroductionTheunsteadyflowpastacircularcylinderhasreceivedagreatdealofattentionowingmainlytoitstheoreticalandpracticalsignificance .Theflowgeneratedbytheoscillationofthecylinder,oroscillatingflowsaroundthecylinder,canbecharacterizedbytwoparameters.OneistheKeulegan_Carpenternumber,definedasKC =UmT/D ,andtheotheristheReynoldsnumberRe=UmD/ν,orafrequencyparameter,definedasβ=D2 / (νT) =Re/KC) ,whichisoftenusedtoreplacetheReynoldsnumberasthesecondparameter.Here,Umisthemaximumvelocityofth…