An experimental study of fluid flow and heat transfer in decaying swirl through a heated annulus

The mean and turbulent structures of turbulent swirling flow in a heated annulus have been measured. Both forced and free vortex swirling flows were generated, and the outer wall of the test section was heated uniformly. The maximum swirl number was 1.39, Reynolds numbers were up to 200000, and heat input was 10.5 kW. Mean and turbulent velocity components, air and wall temperatures, and wall static pressures were all measured. Hot-film techniques were developed to measure turbulence. From these parameters, the flow and temperature fields, pressure distribution, and heat transfer coefficients were determined. The mechanisms of heat transfer were identified.     

双柱单锥型液-液旋流管内流场的激光诊断

应用激光测速仪，对一种双柱单锥型液 液旋流管内的流动结构，进行了全场范围内的多工况流动诊断研究．揭示出其切向速度由内旋流区和外旋流区构成，其中内旋流区中的速度分布符合准强制涡关系，外旋流区中的速度分布符合准自由涡关系；轴向速度由上行流动区和下行流动区构成，两者之间在直管段以零速点作分界，在锥体段则以零速区作过渡并伴随有一定的回流出现，且该过渡区或回流区的大小随锥体截面的收缩而减小，直到进入直管段后消失；各湍流量的分布以管芯处最大向外逐渐减小，沿轴向是直管段中的湍流度大于锥体段中的湍流度，而且湍流度在旋流管内具有各向异性的特性．     

Simultaneous Measurement of Velocity and Fluctuating Pressure in a Turbulent Wing-tip Vortex Using Triple Hot-film Sensor and Miniature Total Pressure Probe

We have developed a probe-system for simultaneous measurement of three velocity components and pressure in turbulent flows. A miniature total pressure probe is placed adjacent to the sensors of a triple hot-film probe in order to achieve the spatial resolution which is equivalent to that of the triple hot-film probe itself. The instantaneous static pressure is calculated from measured velocity and total pressure by means of a newly developed processing method based on the Bernoulli equation for unsteady flows. The measurements were undertaken in a turbulent wing-tip vortex flow. The look-up table method is employed for the calibration of the hot-film probe so accurate velocity data could be obtained over a wide range of the flow-attack angles. It is also demonstrated that the present probe-system is capable of measuring fluctuations in both velocity and pressure in the 20?C650 Hz frequency range. The distribution of the fluctuating pressure obtained by this indirect method is in good agreement with the results from direct measurements of static pressure, demonstrating the promising performance of the present method. Furthermore, an improvement in the ability to make measurements of the velocity?Cpressure correlation across the wing-tip vortex is achieved. This improvement is possible because the effects of lateral velocity components are properly taken into account in the present formulation. The investigation regarding the transport equation budget for turbulent kinetic energy shows an anomalous structure of turbulence in this flow, mainly due to the meandering of the vortex, and the measurement of pressure diffusion is found to play an important role in the characterization of this kind of flow.     

Velocity–pressure correlation measurements in complex free shear flows

Simultaneous measurements of fluctuating velocity and pressure were performed in various turbulent free shear flows including a turbulent mixing layer and the wing-tip vortex trailing from a NACA0012 half-wing. Two different methods for fluctuating static pressure measurement were considered: a direct method using a miniature Pitot tube and an indirect method where static pressure was calculated from total pressure. The pressure obtained by either of these methods was correlated with the velocity measured by an X-type hot-wire probe. The results from these two techniques agreed with each other in the turbulent mixing layer. In the wing-tip vortex case, however, some discrepancies were found, although overall characteristics of the pressure-related statistics were adequately captured by both methods.     

Swirling-strength based large eddy simulation of turbulent flow around single square cylinder at low Reynolds numbers

In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation(LES) is considered as a better simulation model. However, the sub-grid scale(SGS) models reported so far have not ascertained under what flow conditions the LES can lapse into the direct numerical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteristics of zero-incident incompressible turbulent flows around a single square cylinder(SC)at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re = 2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square(RMS) values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls.The instantaneous factor of swirling strength intermittency(FSI) exhibits some laminated structure involved with vortex shedding.     

Three-dimensional flow structures and associated turbulence in the tip region of a waterjet pump rotor blade

Stereo particle image velocimetry measurements focus on the flow structure and turbulence within the tip leakage vortex (TLV) of an axial waterjet pump rotor. Unobstructed optical access to the sample area is achieved by matching the optical refractive index of the transparent pump with that of the fluid. Data obtained in closely spaced planes enable us to reconstruct the 3D TLV structure, including all components of the mean vorticity and strain-rate tensor along with the Reynolds stresses and associated turbulence production rates. The flow in the tip region is highly three-dimensional, and the characteristics of the TLV and leakage flow vary significantly along the blade tip chordwise direction. The TLV starts to roll up along the suction side tip corner of the blade, and it propagates within the passage toward the pressure side of the neighboring blade. A shear layer with increasing length connects the TLV to the blade tip and initially feeds vorticity into it. During initial rollup, the TLV involves entrainment of a few vortex filaments with predominantly circumferential vorticity from the blade tip. Being shed from the blade, these filaments also have high circumferential velocity and appear as swirling jets. The circumferential velocity in the TLV core is also substantially higher than that in the surrounding passage flow, but the velocity peak does not coincide with the point of maximum vorticity. When entrainment of filaments stops in the aft part of the passage, newly forming filaments wrap around the core in helical trajectories. In ensemble-averaged data, these filaments generate a vortical region that surrounds the TLV with vorticity that is perpendicular to that in the vortex core. Turbulence within the TLV is highly anisotropic and spatially non-uniform. Trends can be traced to high turbulent kinetic energy and turbulent shear stresses, e.g., in the shear layer containing the vortex filaments and the contraction region situated along the line where the leakage backflow meets the throughflow, causing separation of the boundary layer at the pump casing. Upon exposure to adverse pressure gradients in the aft part of the passage, at 0.65–0.7 chord fraction in the present conditions, the TLV bursts into a broad turbulent array of widely distributed vortex filaments.     

Phase-locked analysis of velocity fluctuations in a turbulent free swirling jet after vortex breakdown

Large-scale organized vortical structures were studied experimentally in a free swirling jet of air experiencing vortex precession (PVC) at ambient conditions. Detailed measurements were performed in the region near the nozzle exit using phase-locked LDV and PIV, at a Reynolds number of Re ?? 24,400 and a swirl parameter S ?? 1.0. The investigation allowed reconstruction of the time-averaged flowfield, with the associated distribution of turbulent fluctuations, the phase-locked structure of the jet and the associated precessing vortex structure. An original joint analysis of power spectra and probability density functions of velocity data led to quantification of the PVC effect on turbulent fluctuations. This analysis showed that the PVC contribution can be properly separated from the background random turbulence, reproducing the results of phase-locked measurements. It is found that the background turbulence in the near field is substantially weaker if compared to the coherent fluctuations induced by vortex precession.     

Direct and large eddy simulation of stratified turbulence

Simulations of geophysical turbulent flows require a robust and accurate subgrid-scale turbulence modeling. To evaluate turbulence models for stably stratified flows, we performed direct numerical simulations (DNSs) of the transition of the three-dimensional Taylor–Green vortex and of homogeneous stratified turbulence with large-scale horizontal forcing. In these simulations we found that energy dissipation is concentrated within thin layers of horizontal tagliatelle-like vortex sheets between large pancake-like structures. We propose a new implicit subgrid-scale model for stratified fluids, based on the Adaptive Local Deconvolution Method (ALDM). Our analysis proves that the implicit turbulence model ALDM correctly predicts the turbulence energy budget and the energy spectra of stratified turbulence, even though dissipative structures are not resolved on the computational grid.     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

Buoyancy effects on turbulent swirling flows

A three-parameter model of turbulence applicable to free boundary layers has been developed and applied for the prediction of axisymmetric turbulent swirling flows in uniform and stagnant surroundings under the action of buoyancy forces. The turbulent momentum and heat fluxes appearing in the time-averaged equations for the mean motion have been determined from algebraic expressions, derived by neglecting the convection and diffusion terms in the differential transport equations for these quantities, which relate the turbulent fluxes to the kinetic energy of turbulence, k, the dissipation length scale of turbulence, L, and the temperature covariance, T ′². Differential transport equations have been used to determine these latter quantities. The governing equations have been solved using fully implicit finite difference schemes. The turbulence model is capable of reproducing the gross features of pure jet flows, buoyant flows and swirling flows for weak and moderate swirl. The behaviour of a turbulent buoyant swirling jet has been found to depend solely on exit swirl and Froude numbers. The predicted results indicate that the incorporation of buoyancy can cause significant changes in the behaviour of a swirling jet, particularly when the buoyancy strength is high. The jet exhibits similarity behaviour in the initial region for weak swirl and weak buoyancy strengths only, and the asymptotic case of a swirling jet under the action of buoyancy forces is a pure plume in the far field. The predicted results have been found to be in satisfactory agreement with the available experimental data and in good qualitative agreement with other predicted results.     

Proper orthogonal decomposition of the flow in geometries containing a narrow gap

Geometries containing a narrow gap are characterized by strong quasi-periodical flow oscillations in the narrow gap region. The above mentioned phenomena are of inherently unstable nature and, even if no conclusive theoretical study on the subject has been published, the evidence shown to this point suggests that the oscillations are connected to interactions between eddy structures of turbulent flows on opposite sides of the gap. These coherent structures travel in the direction of homogeneous turbulence, in a fashion that strongly recalls a vortex street. Analogous behaviours have been observed for arrays of arbitrarily shaped channels, within certain range of the geometric parameters. A modelling for these phenomena is at least problematic to achieve since they are turbulence driven. This work aims to address the use of Proper Orthogonal Decomposition (POD) to reduce the Navier–Stokes equations to a set of ordinary differential equations and better understand the dynamics underlying these oscillations. Both experimental and numerical data are used to carry out the POD.     

Numerical aspects of calculation of confined swirling flows with internal recirculation

Predictions were performed for two different confined swirling flows with internal recirculation zones. The convection terms in the elliptic governing equations were discretized using three different finite differencing schemes: hybrid, quadratic upwind interpolation and skew upwind differencing. For each flow case, calculations were carried out with these schemes and successively refined grids were employed. For the turbulent flow case the k-ε turbulence model was used. The predicted cases were a laminar swirling flow investigated by Bornstein and Escudier, and a turbulent low-swirl case studied by Roback and Johnson. In both cases an internal recirculation zone was present. The laminar case is well predicted when account is taken of the estimated radial velocity component at the chosen inlet plane. The quadratic upwind interpolation and skew upwind schemes predict the main features of the internal recirculation zone also with a coarse grid. The turbulent case is well predicted with the coarse as well as the finer grids, the skew upwind and quadratic upwind interpolation schemes yielding results very close to the measurements. It is concluded that the skew upwind scheme reaches grid independence slightly before the quadratic upwind scheme, both considerably earlier than the hybrid scheme.     

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.     

Particle image velocimetry investigation of the flow-field of a 3D turbulent annular swirling decaying flow induced by means of a tangential inlet

Particle image velocimetry (PIV) has been used in order to measure the three mean components and turbulence intensities of the velocity vector in a swirling decaying flow induced by a tangential inlet in an annulus. This kind of flow motion is found to be very complex, exhibiting three-dimensional and non-axisymmetric characteristics coupled with a free decay of the swirling intensity along the flow path, thereby making it difficult to study. A method allowing the measurement of the three components of the velocity flow-field with a standard PIV system with two-dimensional acquisitions, is presented. The evolution of each velocity component between the inlet and the outlet of the annulus is obtained. Furthermore, the PIV technique is extended to the measurement of turbulent characteristics such as turbulent intensities and dimensionless turbulent energy. The main characteristics of the swirling flow are discussed and the swirl number is estimated as a function of the axial distance from the tangential inlet. Received: 6 July 1998/Accepted: 20 March 1999     

Vortex dynamics of turbulence–coherent structure interaction

We study the interaction between a coherent structure (CS) and imposed external turbulence by employing direct numerical simulations (DNS) designed for unbounded flows with compact vorticity distribution. Flow evolution comprises (i) the reorganization of turbulence into finer-scale spiral filaments, (ii) the growth of wave-like perturbations within the vortex core, and (iii) the eventual arrest of production, leading to the decay of ambient turbulence. The filaments, preferentially aligned in the azimuthal direction, undergo two types of interactions: parallel filaments pair to form higher-circulation “threads”, and anti-parallel threads form dipoles that self-advect radially outwards. The consequent radial transport of angular momentum manifests as an overshoot of the mean circulation profile—a theoretically known consequence of faster-than-viscous vortex decay. It is found that while the resulting centrifugal instability can enhance turbulence production, vortex decay is arrested by the dampening of the instability due to the “turbulent mixing” caused by instability-generated threads. Ensemble-averaged turbulence statistics show strong fluctuations within the core; these are triggered by the external turbulence, and grow even as the turbulence decays. This surprising growth on a normal-mode-stable vortex results from algebraic amplification through “linear transient growth”. Transient growth is examined by initializing DNS with the “optimal” modes obtained from linear analysis. The simulations show that the growth of transient modes reproduces the prominent dynamics of CS-turbulence interaction: formation of thread-dipoles, growth of core fluctuations, and appearance of bending waves on the column’s core. At the larger Reynolds numbers prevailing in practical flows, transient growth may enable accelerated vortex decay through vortex column breakdown.     

Three-dimensional temporally resolved measurements of turbulence–flame interactions using orthogonal-plane cinema-stereoscopic PIV

A new orthogonal-plane cinema-stereoscopic particle image velocimetry (OPCS-PIV) diagnostic has been used to measure the dynamics of three-dimensional turbulence–flame interactions. The diagnostic employed two orthogonal PIV planes, with one aligned perpendicular and one aligned parallel to the streamwise flow direction. In the plane normal to the flow, temporally resolved slices of the nine-component velocity gradient tensor were determined using Taylor’s hypothesis. Volumetric reconstruction of the 3D turbulence was performed using these slices. The PIV plane parallel to the streamwise flow direction was then used to measure the evolution of the turbulence; the path and strength of 3D turbulent structures as they interacted with the flame were determined from their image in this second plane. Structures of both vorticity and strain-rate magnitude were extracted from the flow. The geometry of these structures agreed well with predictions from direct numerical simulations. The interaction of turbulent structures with the flame also was observed. In three dimensions, these interactions had complex geometries that could not be reflected in either planar measurements or simple flame–vortex configurations.     

Effect of volume energy supply on swirled flows in a subsonic cocurrent stream

The results of a numerical investigation of the effect of thermal energy supply on a swirling viscous heat-conducting gas flow in a subsonic cocurrent stream are presented. The initial stage of development of the swirling flow in the neighborhood of the vortex axis with constant circulation in the outer flow region is considered for two different distributions of the streamwise velocity vector component which simulate a swirling jet-type flow and a wake flow with a streamwise velocity deficit. The effect of local volume energy supply in the neighborhood of the vortex axis, the circulation of the azimuthal velocity component, and the longitudinal pressure gradient in the inviscid stream on the development of the swirling flow and the process of breakdown of cocurrent vortex flows is investigated. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 47–53, November–December, 1998. The work was carried out with financial support from the Russian Foundation for Basic Research (project No. 96-01-00586).     

Investigation of turbulence development in incompressible jets with zonal detached eddy simulation (ZDES) and synthetic turbulent inflow

Hybrid RANS/LES simulations of two incompressible jets are performed with the Zonal Detached Eddy Simulation (ZDES). Two functioning modes of the ZDES for the selection of RANS and DES areas are evaluated, namely the user-defined mode (mode 1) and the global- or automatic- mode (mode 2). The RANS-to-LES transition occurs quickly downstream of the nozzle exit and is found to involve the same physics as a laminar to turbulent transition with vortex pairing near the nozzle exit. The effect of the delay in the RANS-to-LES transition on the jet flow development is analyzed. In particular, the delay in the formation of small-scale turbulent structures results in too high turbulence levels in the mixing layers. Furthermore, it is shown, for two cases, that the injection of synthetic turbulence at the nozzle inlet, originally targeted at reproducing the experimental turbulence level in the jet core, has a significant impact on the mixing layer as it accelerates the RANS-to-LES transition, reduces the spatial wavelength of the vortex pairing and promotes the production of fine-scale turbulence which leads to a better agreement with experiments.     

充分发展圆管湍流的实验研究

采用粒子数字图像测速（digital particle image velocimetry,DPIV)和定量流动显示技术（quantitative flow visualization,QFA)对充分发展的圆管湍流进行了研究。测量结果和直接数值模拟（direct numerical simulation,DNS)结果进行了比较，结果表明作者开发的DPIV技术取得了满意的精度。在此基础上对圆管湍流的动力学机理进行了研究，分析了上抛和下扫在湍流生成中的贡献以及流动显示结构内的脉动速度分布，测量结果显示在圆管湍流的近壁区存在横向强脉冲现象和流动显示所能观察到的结构为上抛占主导地位的结构。     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.