DYNAMICS OF AN IN-LINE TUBE ARRAY SUBJECTED TO STEAM–WATER CROSS-FLOW. PART II: UNSTEADY FLUID FORCES

In this paper, experimental results of unsteady fluid-force measurements are reported. Important deviations of the measured fluid forces from their single-phase flow counterparts were uncovered. Most importantly, the resulting force coefficients are not simple functions of the reduced flow velocity U/fD, as is the case for single-phase flow. Test results at 0·5 MPa challenge the basic assumption of the existence of a time-invariant linear transfer function between tube displacement and the resulting fluid forces. Time–frequency analysis using Wigner–Ville transforms shows that the phase difference between tube displacement and the fluid force (an indicator of stabilizing or destabilizing fluid effects) undergoes significant variation under what may be considered steady flow conditions. This variation may explain the previously reported phenomenon of intermittent fluidelastic instability in two-phase flows.     

DYNAMICS OF AN IN-LINE TUBE ARRAY SUBJECTED TO STEAM–WATER CROSS-FLOW. PART III: FLUIDELASTIC INSTABILITY TESTS AND COMPARISON WITH THEORY

This paper presents the results of comprehensive flow-induced vibration tests conducted on an in-line array in steam–water two-phase flow. The responses of three essentially isolated flexible cylinders at different depths within the array were simultaneously measured. The main test parameters were, ambient pressure (and saturation temperature), in the range 0·5–5·8 MPa, void fraction, 0·70–0·96, and phase flow velocity. Tests reported here were conducted simultaneously with the damping tests reported in Part I of this study. At the highest pressures (3·0 and 5·8 MPa), strong instabilities, in homogeneous flow akin to single-phase flow occurred. The test tube located in the central region of the array was the most susceptible to instability. This was attributed partly to reduced two-phase damping deep in the array, while differences in local fluid forces at different locations in the array are not ruled out. The flow at 0·5 MPa was a nonhomogeneous intermittent slug-type flow. Strong turbulence excitation obscured clear fluidelastic instability; intermittent instability was, however, ascertained. Stability boundary calculations were done using unsteady fluid forces presented in Part II of this series of papers. Results for the case of P=5·8 MPa show good agreement with the measured instability boundary.     

Modelling the time-dependent flow over riblets in the viscous wall region

The flow over drag reducing riblets is examined computationally using a time-dependent model of the viscous wall region. The flow at the upper bound of the computational domain (y ⁺?40) is described using a streamwise eddy model consisting of two scales. A control-volume finite-element method utilizing triangular meshes is used to exactly fit the riblet cross-sectional geometry. Observations of the transient flow conditions suggest that the riblets limit the lateral spread of fluid inrushes towards the wall and retain low momentum fluid in the riblet valleys effectively isolating much of the wall from such inrushes. The generation of intermittent secondary vortices within the riblet valleys also occurs; however, these appear to be quite weak and fairly short-lived.     

Fence probe measurement of flow reversal in separating turbulent boundary layers

A small fence probe was evaluated for measurements in the time-dependent flow reversal region of the transition from boundary layer to separated flow. For moderate and high Reynolds numbers, the fence probe is demonstrated to be a usable tool for the measurement of the reverse flow associated with separation. Although the present probe pressure transducer system was limited to approximately 200?Hz, pulses of positive and negative shear stress were readily detected. At or near the location of zero surface shear stress, the measurements were limited by the signal-to-noise ratio. For the separated flow investigated, a marked reduction in the pressure gradient occurred when the fence probe indicated approximately 20?% reversal for the higher Reynolds numbers. The reversal increased to 24?% for the lower Reynolds numbers. The measurements indicate that flow reversal alone may not be adequate to identify the degree of separation. Upstream of turbulent boundary layer (intermittent) separation, the duration of the reversed shear stress was found to be very short (0.002?C0.007?s), suggesting a local, small-scale, impulse-type separation. At and beyond the location of intermittent separation, the shear stress reversal duration was an order of magnitude longer. Estimates of the maximum and minimum surface shear stress in the separation region were also obtained with the fence probe.     

Investigation of oil-water flow regimes and pressure drops in mini-channels

Oil-water flow regimes were studied in 2.1 mm and 3.7 mm borosilicate glass tubes; both tubes exhibit Eötvös numbers less than one and therefore surface tension forces may be more important in these mini-channels compared to larger diameter tubes. A closed-loop, adiabatic experimental apparatus was constructed and validated using water. This study focused on tap water and two mineral oils (i.e., Parol 70 and 100) with a density of 840 kg/m³ but a factor of two difference in viscosity. Experiments included a wide range of oil superficial velocities (e.g., 0.84–6.84 m/s for D = 2.1 mm and 0.27–3.30 m/s for D = 3.7 mm) and water superficial velocities (e.g., 0.21–7.69 m/s for D = 2.1 mm and 0.07–4.96 m/s for D = 3.7 mm). Stratified, annular, intermittent, and dispersed flow regimes were observed in both tubes, although the annular flow regime was more prevalent in the smaller tube. Pressure drops increased with decreasing tube diameter and were flow regime dependent. Flow maps were created for these mini-channels and equations adapted from Brauner and Maron (1999) were used to predict the flow regime transitions. The effects of viscosity were modest, although increased oil viscosity enhanced stability of oil-water flows.     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

Visualization and structure of confined, milliscale, unsteady impinging slot jets and associated vortices

Flow characteristics of confined, laminar milliscale slot jets are investigated from visualizations, as they impinge upon a flat target plate, with a fully developed velocity profile at the nozzle exit. The effects of Reynolds number Re and normalized nozzle-to-plate distance H/B are considered for a nozzle width B of 1.0 mm. Transition from a stable symmetric jet to an unsteady oscillating jet is observed as the Reynolds number increases (with H/B constant), where the Reynolds number associated with this transition decreases as the normalized nozzle-to-plate distance H/B increases. Instantaneous visualizations show unsteady lateral distortions of jet columns at experimental conditions corresponding to the presence of continuous sinusoidal oscillations, intermittent oscillating motion of the jet column, and jet flow fluctuation/flapping motion. Also apparent in flow visualization sequences are smoke signatures associated with instantaneous vortex structures, which form as secondary flows develop in fluid which, initially, is just adjacent to and within the jet column. Associated jet and vortex structural changes are described as different modes of unsteadiness are present, including characterization of jet column unsteadiness using jet column oscillation frequency, and lateral and streamwise extents of jet distortion.     

Oil mist transport process in a long pipeline on turbulent flow transition region

Internal gas velocity fluctuations and their effects on the mist diffusion process were examined in a long horizontal pipe to understand oil mist transportation, particularly in the laminar-to-turbulent flow transition region. Three hot-wire anemometers and aerosol concentration monitors were used to deduce these effects as the two-phase mist flow gradually developed in the stream-wise direction. We found significant axial mist diffusion at Reynolds numbers (Re) < 1000 because of passive scalar transport by Poiseuille flow. However, this diffusion was restricted by the non-zero inertia of the mist at a Stokes number, O(10⁻⁵), relying on the Brownian motion of the mist. At Re > 2400, a sharp mist waveform was maintained by a turbulent flow with active radial mixing. New data were obtained within the range of 1000 < Re < 2400, which cannot be explained by interpolation between the above-mentioned two states. The mist concentration displays multiple temporal peaks at Re < 2000 owing to perturbations of localized turbulence as well as radial anisotropy as being conveyed more than 2000-diameters in distance. This behavior is caused by intermittent disturbances induced by the pipe wall roughness, which sharply distorts the wall-aligned laminar mist layer left by parabolic axial stretching of local laminar flow.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

The turbulent/non-turbulent interface in an adverse pressure gradient turbulent boundary layer

The turbulent/non-turbulent interface (TNTI) in an adverse pressure gradient (APG, β = 1.45) turbulent boundary layer (TBL) is explored here by using direct numerical simulation (DNS) data; β is the Clauser pressure gradient parameter. For comparison, the DNS data for a zero pressure gradient (ZPG) TBL is included. The interface is extracted with an approach based on enstrophy criteria. Depending on the enstrophy, the outer boundary layer flow can be classified into the free stream, boundary layer wake, and intermittent flow regimes. The fractal dimension of the interface is obtained by using the box-counting algorithm, and was found to be constant over a long range of box sizes. The TNTI shows a monofractal behavior. The geometric complexity of a TNTI can be determined in terms of the genus, which is defined as the number of handles in a geometric object. We examine the volume and projection area of the genus of the TNTI to analyze the entrainment process. The geometric complexity of the APG TBL interface and the local entrainment are greater than those of the ZPG TBL, as is evident in the increases in the genus near the interface. The local entrainment velocity is dominantly affected by the viscous diffusion at the interface.     

Multiphase drag reduction: Effect of eliminating slugs

Two-phase pressure drop measurements were taken for air/water mixtures in a 0.052-m diameter horizontal pipe with special focus on the superficial liquid velocity range of 0.03–1.2 m/s at superficial gas velocities of 3.8, 5.2, and 6.6 m/s. It was found that the addition of 400 ppm of sodium dodecyl sulfate (SDS) to the water reduced the pressure drop by 25–40% when compared to equal flow rates without SDS. The pressure drop reduction occurred where the SDS eliminated the occurrence of the intermittent flow present with water. It was also found that the same concentration of SDS had virtually no effect on single phase liquid pressure drop. The pressure drop reduction appears to be due solely to the suppression of intermittent flow patterns.     

Unsteady separated and reattaching turbulent flow over a two-dimensional square rib

The spatio-temporal characteristics of the separated and reattaching turbulent flow over a two-dimensional square rib were studied experimentally. Synchronized measurements of wall-pressure fluctuations and velocity fluctuations were made using a microphone array and a split-fiber film, respectively. Profiles of time-averaged streamwise velocity and wall-pressure fluctuations showed that the shear layer separated from the leading edge of the rib sweeps past the rib and directly reattaches on the bottom wall (x/H=9.75) downstream of the rib. A thin region of reverse flow was formed above the rib. The shedding large-scale vortical structures (fH/U₀=0.03) and the flapping separation bubble (fH/U₀=0.0075) could be discerned in the wall-pressure spectra. A multi-resolution analysis based on the maximum overlap discrete wavelet transform (MODWT) was performed to extract the intermittent events associated with the shedding large-scale vortical structures and the flapping separation bubble. The convective dynamics of the large-scale vortical structures were analyzed in terms of the autocorrelation of the continuous wavelet-transformed wall pressure, cross-correlation of the wall-pressure fluctuations, and the cross-correlation between the wall pressure at the time-averaged reattachment point and the streamwise velocity field. The convection speeds of the large-scale vortical structures before and after the reattachment point were U_c=0.35U₀ and 0.45U₀, respectively. The flapping motion of the separation bubble was analyzed in terms of the conditionally averaged reverse-flow intermittency near the wall region. The instantaneous reattachment point in response to the flapping motion was obtained; these findings established that the reattachment zone was a 1.2H-long region centered at x/H=9.75. The reverse-flow intermittency in one period of the flapping motion demonstrated that the thin reverse flow above the rib is influenced by the flapping motion of the separation bubble behind the rib.     

Two-phase flow characteristics in gas–liquid microreactors

Multiphase chemical microreactors require a detailed knowledge of the flow conditions inside the reaction system. This paper reports flow visualization measurements of the two-phase gas–liquid flow pattern and the liquid velocity distribution inside liquid plugs of an intermittent flow. Rectangular cross-section silicon microchannels with hydraulic diameters between 187.5 and 218 μm are fabricated. Laser Induced Fluorescence (LIF) is used to determine the flow pattern. To analyze the influence of the liquid properties and the channel diameter on the two-phase flow pattern, we present flow regime maps using different channel geometries and fluids. A universal flow pattern map based on dimensional analysis is presented. In contrast to microchannel flows, a great number of correlations for flow characteristics for multiphase flow in (round) pipes with diameters >1 mm exist. We compare our experimental results from optical flow visualizations in microreactors with common flow correlations and regime maps for macro- and microchannels. The recirculation motion in the liquid segments of an intermittent gas–liquid flow is analyzed using micron-resolution particle image velocimetry (μPIV). The velocity distribution influences the mixing and the mass transport towards the reactive phase interface dealing with two-phase chemical reactions. For straight microchannels hardly any mass transport over the center line is quantified. For enhanced mixing geometrical adaptations are suggested.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Intermittent patterns in mixed convection

We have performed a Laser Doppler anemometry study of the mixed convective flow of water in a horizontal, rectangular duct of 4.15 mm height, 115 mm length and 15.05 mm width. This study has revealed a particular flow regime that lies approximately between 1920相似文献   

Influence of wall proximity on characteristics of wake behind a square cylinder: PIV measurements and POD analysis

Influence of wall proximity on characteristics of the wake behind a two-dimensional square cylinder was experimentally studied in the present work. A low-speed recirculation water channel was established for the experiment; the Reynolds number based on the free-stream velocity and cylinder width (D) was kept at Re_D = 2250. Four cases with different gap width, e.g., G/D = 0.1, 0.2, 0.4 and 0.8, were chosen for comparison. Two experimental techniques, e.g., the standard PIV with high image-density CCD camera and TR-PIV with a high-speed camera were employed in measuring the wake field, enabling a comprehensive view of the time-averaged wake pattern at high spatial resolution and the instantaneous flow field at high temporal resolution, respectively. For the four cases, the difference in spatial characteristics of the wake in the vicinity of the plane wall was analyzed in terms of the time-averaged quantities measured by the standard PIV, e.g., the streamline pattern, the vector field, the streamwise velocity fluctuation intensity and the reverse-flow intermittency. The proper orthogonal decomposition (POD) method was extensively used to decompose the TR-PIV measurements, giving a close-up view of the energetic POD modes buried in the wake. The low-order flow model of the wake at G/D = 0.8 and 0.4 was constructed by using the linear combination of the first two POD modes and the time-mean flow field, which reflected well the vortex shedding process in the sense of the phase-dependent patterns. The intermittent appearance of the weakly separated region near the wall was found at G/D = 0.4. On going from G/D = 0.8 to 0.4, the remarkable variation of the instantaneous wake in the longitudinal direction confirmed that the wall constraint stretches the vortices in the plane of the wall and transfers the energy to the longitudinal component at the expense of the lateral one.     

Intermittency and stochastic modeling of hydrodynamic pressure fluctuations in the near field of compressible jets

The intermittent statistics of the pressure fluctuations measured in the near field of a compressible jet are investigated under several flow conditions. An experiment is carried out in a semi-anechoic chamber on a single stream compressible jet at Mach numbers varying from 0.5 to 0.9 and measuring the fluctuating pressure in several positions in the near field. The main quantities analyzed are the intermittent time and the energy amplitude of events that are extracted from the experimental data through a wavelet-based tracking algorithm. As an extension of a previous paper (Camussi et al., 2017), low order statistical moments (mean and variance) and Probability Distribution Functions are parametrized in terms of three relevant quantities characterizing the jet flow physics: the Mach number, the radial distance from the jet axis and the axial position. It is observed that the non-dimensional statistical quantities are weakly dependent upon the flow conditions allowing for simple stochastic models to be introduced on the basis of suitable fittings of averaged statistical properties and of the Probability Distribution Functions.     

Analysis of dryout in horizontal and inclined tubes

The infrared thermography technique was used to study the thermal and hydrodynamic phenomena in intermittent two-phase air–water flow in horizontal and inclined tubes at atmospheric pressure. The study was aimed at elucidating the relationship between the hydrodynamic parameters and dryout phenomena. It focuses on the empirical evaluation of the wall temperature distribution in a uniformly heated pipe. The results reveal the existence of dryout phenomena in horizontal pipe flow only. The flow parameter based on the frequency, length and velocity of elongated bubble is presented for the prediction of dryout.     

Inertial Particles in Homogeneous Shear Turbulence: Experiments and Direct Numerical Simulation

The properties of the transport of heavy inertial particles in a uniformly sheared turbulent flow have been investigated by combining experimental and numerical data at particle Stokes number St?≈?0.3?÷?0.5 respectively. As in isotropic turbulence, particles are observed to avoid zones of intense enstrophy and to cluster in strain-dominated regions, resulting in highly intermittent spatial distributions. Moreover, the anisotropy of the mean flow is found to imprint a clear preferential orientation of the particle clusters in the direction of the maximum mean strain. These features are observed both in the numerics and in the experiments, and have been consistently quantified by a number of complementary statistical tools, such as the Voronoï tessellations and the pair correlation function. The latter quantity has been generalized in the form of the Angular Distribution Function and has allowed to evaluate the anisotropy content of the particle field at each scale. The behavior of this observable exhibits the same trend in the two datasets and suggests that, owing to increased inertia, the particle distribution starts to recover isotropy at scales smaller than the carrier velocity field. A proper rescaling of the two datasets in terms of their respective values of the shear scale allows to account for differences in the Reynolds number of experiments and numerics in the range of scales dominated by the mean shear.