The effect of confinement on the motion of a single clean bubble

2008 ◽  
Vol 616 ◽  
pp. 419-443 ◽  
Author(s):  
B. FIGUEROA-ESPINOZA ◽  
R. ZENIT ◽  
D. LEGENDRE
Keyword(s):  
Reynolds Number ◽  
Lift Force ◽  
Reasonable Agreement ◽  
Bubble Radius ◽  
Reynolds Numbers ◽  
Bubble Surface ◽  
Total Drag ◽  
High Reynolds ◽  
Main Effects ◽  
Gap Width

The effect of confining a gas bubble between two parallel walls was investigated for the inertia-dominated regime characterized by high Reynolds and low Weber numbers. Single bubble experiments were performed with non-polar liquids such that the bubble surface could be considered clean; hence, shear free. The drag coefficient was found to be the result of two main effects: the Reynolds number and the confinement. The total drag could be written as the product of the corresponding unconfined drag, which depended mainly on the Reynolds number, and a function F(s)=1 + κs3. The confinement parameter s was defined as the ratio of the bubble radius to the gap width. The value of the constant κ depended on the way in which the bubbles moved within the gap, which was found to be either in a rectilinear (κ≈8) or oscillatory trajectory (κ≈80). For Re < 70, and a range of values of the confinement parameter, the bubbles followed a rectilinear path. For this regime, numerical simulations were performed to obtain the drag force on the bubble directly; a reasonable agreement was found with experiments. Moreover, a comparison of these results with a potential-flow-based model indicated that the vorticity produced at the walls induced a significant part of the drag. For Re > 70, oscillations were observed in the bubble trajectory. In all cases, the oscillation occurred in a zigzag manner. Near the transition the bubbles oscillated but did not reach the walls; for larger Reynolds numbers, the bubbles collided repeatedly with the walls as they ascended. The instability, which is different from the well-known unconfined path instability, resulted from the reversal of sign of the wall-induced lift force: for low Reynolds number, the walls have a stabilizing effect because of the repulsive nature of the lift force between the walls and the bubble, while for high Reynolds number the lift is attractive and trajectories become unstable. Considering a model for the lift force of a bubble moving near a wall, the conditions for the transition were identified. A reasonable agreement between the model and experiments was found.

The lift force on a spherical bubble in a viscous linear shear flow

1998 ◽  
Vol 368 ◽  
pp. 81-126 ◽  
Author(s):  
DOMINIQUE LEGENDRE ◽  
JACQUES MAGNAUDET
Keyword(s):  
Reynolds Number ◽  
Shear Rate ◽  
Shear Flow ◽  
Lift Force ◽  
Lift Coefficient ◽  
Reynolds Numbers ◽  
Bubble Surface ◽  
Vorticity Field ◽  
Linear Shear ◽  
High Reynolds

The three-dimensional flow around a spherical bubble moving steadily in a viscous linear shear flow is studied numerically by solving the full Navier–Stokes equations. The bubble surface is assumed to be clean so that the outer flow obeys a zero-shear-stress condition and does not induce any rotation of the bubble. The main goal of the present study is to provide a complete description of the lift force experienced by the bubble and of the mechanisms responsible for this force over a wide range of Reynolds number (0.1[les ]Re[les ]500, Re being based on the bubble diameter) and shear rate (0[les ]Sr[les ]1, Sr being the ratio between the velocity difference across the bubble and the relative velocity). For that purpose the structure of the flow field, the influence of the Reynolds number on the streamwise vorticity field and the distribution of the tangential velocities at the surface of the bubble are first studied in detail. It is shown that the latter distribution which plays a central role in the production of the lift force is dramatically dependent on viscous effects. The numerical results concerning the lift coefficient reveal very different behaviours at low and high Reynolds numbers. These two asymptotic regimes shed light on the respective roles played by the vorticity produced at the bubble surface and by that contained in the undisturbed flow. At low Reynolds number it is found that the lift coefficient depends strongly on both the Reynolds number and the shear rate. In contrast, for moderate to high Reynolds numbers these dependences are found to be very weak. The numerical values obtained for the lift coefficient agree very well with available asymptotic results in the low- and high-Reynolds-number limits. The range of validity of these asymptotic solutions is specified by varying the characteristic parameters of the problem and examining the corresponding evolution of the lift coefficient. The numerical results are also used for obtaining empirical correlations useful for practical calculations at finite Reynolds number. The transient behaviour of the lift force is then examined. It is found that, starting from the undisturbed flow, the value of the lift force at short time differs from its steady value, even when the Reynolds number is high, because the vorticity field needs a finite time to reach its steady distribution. This finding is confirmed by an analytical derivation of the initial value of the lift coefficient in an inviscid shear flow. Finally, a specific investigation of the evolution of the lift and drag coefficients with the shear rate at high Reynolds number is carried out. It is found that when the shear rate becomes large, i.e. Sr=O(1), a small but consistent decrease of the lift coefficient occurs while a very significant increase of the drag coefficient, essentially produced by the modifications of the pressure distribution, is observed. Some of the foregoing results are used to show that the well-known equality between the added mass coefficient and the lift coefficient holds only in the limit of weak shears and nearly steady flows.

scholarly journals Reversal of the lift force on an oblate bubble in a weakly viscous linear shear flow

2009 ◽  
Vol 628 ◽  
pp. 23-41 ◽  
Author(s):  
RICHARD ADOUA ◽  
DOMINIQUE LEGENDRE ◽  
JACQUES MAGNAUDET
Keyword(s):  
Shear Flow ◽  
Lift Force ◽  
Reynolds Numbers ◽  
Bubble Surface ◽  
Relative Shear ◽  
Linear Shear ◽  
Wake Instability ◽  
Reversal Mechanism ◽  
Bubble Rising ◽  
High Reynolds

We compute the flow about an oblate spheroidal bubble of prescribed shape set fixed in a viscous linear shear flow in the range of moderate to high Reynolds numbers. In contrast to predictions based on inviscid theory, the numerical results reveal that for weak enough shear rates, the lift force and torque change sign in an intermediate range of Reynolds numbers when the bubble oblateness exceeds a critical value that depends on the relative shear rate. This effect is found to be due to the vorticity generated at the bubble surface which, combined with the velocity gradient associated with the upstream shear, results in a system of two counter-rotating streamwise vortices whose sign is opposite to that induced by the classical inviscid tilting of the upstream vorticity around the bubble. We show that this lift reversal mechanism is closely related to the wake instability mechanism experienced by a spheroidal bubble rising in a stagnant liquid.

scholarly journals Aerodynamics of smooth and rough square-section prisms at incidence in very high Reynolds-number cross-flows

2021 ◽  
Vol 62 (3) ◽  
Author(s):  
Nils Paul van Hinsberg
Keyword(s):  
Reynolds Number ◽  
Cross Sections ◽  
Reynolds Numbers ◽  
Circular Cylinders ◽  
Surface Boundary ◽  
Experimental Proof ◽  
Surface Boundary Layer ◽  
Pitch Moment ◽  
The Mean ◽  
High Reynolds

Abstract The aerodynamics of smooth and slightly rough prisms with square cross-sections and sharp edges is investigated through wind tunnel experiments. Mean and fluctuating forces, the mean pitch moment, Strouhal numbers, the mean surface pressures and the mean wake profiles in the mid-span cross-section of the prism are recorded simultaneously for Reynolds numbers between 1$$\times$$ × 10$$^{5}$$ 5 $$\le$$ ≤ Re$$_{D}$$ D $$\le$$ ≤ 1$$\times$$ × 10$$^{7}$$ 7 . For the smooth prism with $$k_s$$ k s /D = 4$$\times$$ × 10$$^{-5}$$ - 5 , tests were performed at three angles of incidence, i.e. $$\alpha$$ α = 0$$^{\circ }$$ ∘ , −22.5$$^{\circ }$$ ∘ and −45$$^{\circ }$$ ∘ , whereas only both “symmetric” angles were studied for its slightly rough counterpart with $$k_s$$ k s /D = 1$$\times$$ × 10$$^{-3}$$ - 3 . First-time experimental proof is given that, within the accuracy of the data, no significant variation with Reynolds number occurs for all mean and fluctuating aerodynamic coefficients of smooth square prisms up to Reynolds numbers as high as $$\mathcal {O}$$ O (10$$^{7}$$ 7 ). This Reynolds-number independent behaviour applies to the Strouhal number and the wake profile as well. In contrast to what is known from square prisms with rounded edges and circular cylinders, an increase in surface roughness height by a factor 25 on the current sharp-edged square prism does not lead to any notable effects on the surface boundary layer and thus on the prism’s aerodynamics. For both prisms, distinct changes in the aerostatics between the various angles of incidence are seen to take place though. Graphic abstract

Temperature Control of Blow-Down, Supersonic Tunnels

1956 ◽  
Vol 60 (541) ◽  
pp. 67-70
Author(s):  
T. A. Thomson
Keyword(s):  
Heat Transfer ◽  
Boundary Layer ◽  
Reynolds Number ◽  
Mach Number ◽  
Reynolds Numbers ◽  
Stagnation Temperature ◽  
Blow Down ◽  
High Reynolds Numbers ◽  
Experimental Errors ◽  
High Reynolds

The blow-down type of intermittent, supersonic tunnel is attractive because of its simplicity and because relatively high Reynolds numbers can be obtained for a given size of test section. An adverse characteristic, however, is the fall of stagnation temperature during runs, which can affect experiments in several ways. The Reynolds number varies and the absolute velocity is not constant, even if the Mach number and pressure are; heat-transfer cannot be studied under controlled conditions and the experimental errors arising from the effect of heat-transfer on the boundary layer vary in time. These effects can become significant in quantitative experiments if the tunnel is large and the variation of temperature very rapid; the expense required to eliminate them might then be justified.

Reynolds Number Effects on the Freewheeling Behavior of a High Solidity Vertical River Kinetic Turbine

2010 ◽  
Author(s):  
Amir Hossein Birjandi ◽  
Eric Bibeau
Keyword(s):  
Reynolds Number ◽  
Reynolds Numbers ◽  
Speed Ratio ◽  
Stream Velocity ◽  
Blade Profile ◽  
Low Reynolds Numbers ◽  
Tip Speed Ratio ◽  
Reynolds Number Effects ◽  
High Reynolds ◽  
The University

A four-bladed, squirrel-cage, and scaled vertical kinetic turbine was designed, instrumented and tested in the water tunnel facilities at the University of Manitoba. With a solidity of 1.3 and NACA0021 blade profile, the turbine is classified as a high solidity model. Results were obtained for conditions during freewheeling at various Reynolds numbers. In this study, the freewheeling tip speed ratio, which relates the ratio of maximum blade speed to the free stream velocity at no load, was divided into three regions based on the Reynolds number. At low Reynolds numbers, the tip speed ratio was lower than unity and blades were in a stall condition. At the end of the first region, there was a sharp increase of the tip speed ratio so the second region has a tip speed ratio significantly higher than unity. In this region, the tip speed ratio increases almost linearly with Reynolds number. At high Reynolds numbers, the tip speed ratio is almost independent of Reynolds number in the third region. It should be noted that the transition between these three regions is a function of the blade profile and solidity. However, the three-region behavior is applicable to turbines with different profiles and solidities.

Stability of non-parabolic flow in a flexible tube

1999 ◽  
Vol 395 ◽  
pp. 211-236 ◽  
Author(s):  
V. SHANKAR ◽  
V. KUMARAN
Keyword(s):  
Reynolds Number ◽  
Velocity Profile ◽  
Asymptotic Analysis ◽  
High Reynolds Number ◽  
Reynolds Numbers ◽  
Flexible Tube ◽  
Developing Flow ◽  
Parabolic Flow ◽  
High Reynolds ◽  
The Stability

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such ‘non-parabolic’ flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

Turbulent Annular Pipe Flow in Subcritical Transition Regime: Direct Numerical Simulation of a Large Domain

2015 ◽  
Author(s):  
Takahiro Ishida ◽  
Takahiro Tsukahara
Keyword(s):  
Reynolds Number ◽  
Poiseuille Flow ◽  
Reynolds Numbers ◽  
Radius Ratio ◽  
Transition Regime ◽  
Transitional Regime ◽  
Large Domain ◽  
High Reynolds ◽  
Annular Pipe ◽  
Transition Scenario

We performed direct numerical simulations of annular Poiseuille flow (APF) with a radius ratio of η (= rin/rout) = 0.8, in order to investigate the subcritical transition scenario from the developed turbulent state to the laminar state. In previous studies on annular Poiseuille flow, the flows at high Reynolds numbers were well examined and various turbulence statistics were obtained for several η, because of their dependence on η. Since the transitional APF is still unclear, we investigate annular Poiseuille flows in the transitional regime through the large-domain simulations in a range of the friction Reynolds number from Reτ = 150 down to 56. At a transitional Reynolds number, weak-fluctuation regions occur intermittently and regularly in the flow field, and the localized turbulence appears in the form of banded patterns same as in plane Poiseuille flow (PPF). The flow system of APF with a high radius ratio η ≈ 1 can be regarded as PPF and, hence, the transition regime in high radius-ratio of APF and in PPF should be analogous. However, in APF, the banded structure takes on helical shape around the inner cylinder, since APF is a closed system in the spanwise (azimuthal) direction. In this paper, the (dis-)similarity between APF and PPF is discussed.

Tip Vortex Cavitation Inception Scaling for High Reynolds Number Applications

2009 ◽  
Vol 131 (7) ◽  
Author(s):  
Young T. Shen ◽  
Scott Gowing ◽  
Stuart Jessup
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Boundary Layer Thickness ◽  
Reynolds Numbers ◽  
Present Theory ◽  
Full Scale ◽  
Tip Vortex ◽  
Tip Vortex Cavitation ◽  
Cavitation Inception ◽  
High Reynolds

Tip vortices generated by marine lifting surfaces such as propeller blades, ship rudders, hydrofoil wings, and antiroll fins can lead to cavitation. Prediction of the onset of this cavitation depends on model tests at Reynolds numbers much lower than those for the corresponding full-scale flows. The effect of Reynolds number variations on the scaling of tip vortex cavitation inception is investigated using a theoretical flow similarity approach. The ratio of the circulations in the full-scale and model-scale trailing vortices is obtained by assuming that the spanwise distributions of the section lift coefficients are the same between the model-scale and the full-scale. The vortex pressure distributions and core sizes are derived using the Rankine vortex model and McCormick’s assumption about the dependence of the vortex core size on the boundary layer thickness at the tip region. Using a logarithmic law to describe the velocity profile in the boundary layer over a large range of Reynolds number, the boundary layer thickness becomes dependent on the Reynolds number to a varying power. In deriving the scaling of the cavitation inception index as the ratio of Reynolds numbers to an exponent m, the values of m are not constant and are dependent on the values of the model- and full-scale Reynolds numbers themselves. This contrasts traditional scaling for which m is treated as a fixed value that is independent of Reynolds numbers. At very high Reynolds numbers, the present theory predicts the value of m to approach zero, consistent with the trend of these flows to become inviscid. Comparison of the present theory with available experimental data shows promising results, especially with recent results from high Reynolds number tests. Numerical examples of the values of m are given for different model- to full-scale sizes and Reynolds numbers.

Spectral energy transfer in high Reynolds number turbulence

1977 ◽  
Vol 79 (2) ◽  
pp. 337-359 ◽  
Author(s):  
K. N. Helland ◽  
C. W. Van Atta ◽  
G. R. Stegen
Keyword(s):  
Reynolds Number ◽  
Energy Transfer ◽  
Energy Spectrum ◽  
High Reynolds Number ◽  
Reynolds Numbers ◽  
Second Order ◽  
Spectral Energy ◽  
Local Isotropy ◽  
High Reynolds ◽  
Good Agreement

The spectral energy transfer of turbulent velocity fields has been examined over a wide range of Reynolds numbers by experimental and empirical methods. Measurements in a high Reynolds number grid flow were used to calculate the energy transfer by the direct Fourier-transform method of Yeh & Van Atta. Measurements in a free jet were used to calculate energy transfer for a still higher Reynolds number. An empirical energy spectrum was used in conjunction with a local self-preservation approximation to estimate the energy transfer at Reynolds numbers beyond presently achievable experimental conditions.Second-order spectra of the grid measurements are in excellent agreement with local isotropy down to low wavenumbers. For the first time, one-dimensional third-order spectra were used to test for local isotropy, and modest agreement with the theoretical conditions was observed over the range of wavenumbers which appear isotropic according to second-order criteria. Three-dimensional forms of the measured spectra were calculated, and the directly measured energy transfer was compared with the indirectly measured transfer using a local self-preservation model for energy decay. The good agreement between the direct and indirect measurements of energy transfer provides additional support for both the assumption of local isotropy and the assumption of self-preservation in high Reynolds number grid turbulence.An empirical spectrum was constructed from analytical spectral forms of von Kármán and Pao and used to extrapolate energy transfer measurements at lower Reynolds number to Rλ = 105 with the assumption of local self preservation. The transfer spectrum at this Reynolds number has no wavenumber region of zero net spectral transfer despite three decades of $k^{-\frac{5}{3}}$. behaviour in the empirical energy spectrum. A criterion for the inertial subrange suggested by Lumley applied to the empirical transfer spectrum is in good agreement with the $k^{-\frac{5}{3}}$ range of the empirical energy spectrum.

Sign in / Sign up

Export Citation Format

Share Document