Search

In this work we study the heat transport in inhomogeneous bubbly flow. The experiments were performed in a rectangular bubble column heated from one side wall and cooled from the other, with millimetric bubbles introduced through one half of the injection section (close to the hot wall or close to the cold wall). We characterise the global heat transport while varying two parameters: the gas volume fraction $\alpha = 0.4\% - 5.1 \%$, and the Rayleigh number $Ra_H=4\times10^9-2.2\times10^{10}$. As captured by imaging and characterised using Laser Doppler Anemometry (LDA), different flow regimes occur with increasing gas flow rates. In the generated inhomogeneous bubbly flow there are three main contributions to the mixing: (\textit{i}) transport by the buoyancy driven recirculation, (\textit{ii}) bubble induced turbulence (BIT) and (\textit{iii}) shear-induced turbulence (SIT). The strength of these contributions and their interplay depends on the gas volume fraction which is reflected in the measured heat transport enhancement. We compare our results with the findings for heat transport in homogeneous bubbly flow from Gvozdi{\'c} \emph{et al.} (2018). We find that for the lower gas volume fractions ($\alpha4\%$, when the contribution of SIT becomes stronger, but so does the competition between all three contributions, the homogeneous injection is more efficient.

We report on turbulent thermal convection experiments in a rotating cylinder with a diameter () to height () aspect ratio of. Using nitrogen and pressurised sulphur hexafluoride we cover Rayleigh numbers (Ra) from to at Prandtl numbers. For these Ra we measure the global vertical heat flux (i.e. the Nusselt number - Nu), as well as statistical quantities of local temperature measurements, as a function of the rotation rate, i.e. the inverse Rossby number - 1/Ro. In contrast to measurements in fluids with a higher Pr we do not find a heat transport enhancement, but only a decrease of Nu with increasing 1/Ro. When normalised with Nu(0) for the non-rotating system, data for all different Ra collapse and, for sufficiently large 1/Ro, follow a power law. Furthermore, we find that both the heat transport and the temperature field qualitatively change at rotation rates and. We interpret the first transition at as change from a large-scale circulation roll to the recently discovered boundary zonal flow (BZF). The second transition at rotation rate is not associated with a change of the flow morphology, but is rather the rotation rate for which the BZF is at its maximum. For faster rotation the vertical transport of warm and cold fluid near the sidewall within the BZF decreases and hence so does Nu.

Thermal convection is of major importance in various astro- and geophysical systems, exemplary are buoyancy driven flows in the atmosphere or in the stellar interior. It has been studied for decades in an idealized model system - the Rayleigh-Bénard convection (RBC) - which consists of a horizontal fluid layer heated at the bottom and cooled at the top. Within the Oberbeck-Boussinesq approximation this system is controlled by two parameters only. These are the Rayleigh number (Ra), which represents the thermal driving and the Prandtl number (Pr) that relates the momentum and thermal diffusivities of the fluid. Convection flows in geo- and astrophysics are often influenced by Coriolis forces due to the rotation of the planet or the star. In RBC-system, Coriolis forces are introduced by rotating the convection cell around its vertical axis. The rotation is expressed by an additional dimensionless control parameter, i.e., the inverse Rossby number 1/Ro. We study experimentally the influence of rotation on the heat transport and the temperature field at very large Ra in the High Pressure Convection Facility (HPCF) in Göttingen. The facility consists of a cylindrical cell of 1.10m diameter and 2.20m height that is filled with pressurized sulfur hexafluoride (SF6) at up to 19bar. The height of the cell and the large density of SF6 enable us to reach very large Ra (up to 8×1014) at 0.74We find a monotonic decrease of the heat transport with increasing rotation rate. Furthermore, we measure quantities of the flow close to the lateral side walls of the convection cylinder. For large rotation rates we analyze this as part of the recently proposed “Boundary Zonal Flow” (BZF), where the vertical heat transport is enhanced and warm (cold) up (down) flow self-organizes in a periodic manner. In the experiment we observe the BZF most notably in the probability density function of the temperature, which develops a bimodal Gaussian distribution. We also find that the periodic warm-cold structure drifts in anti-cyclonic direction and thus form traveling waves of the temperature field.

For rapidly rotating turbulent Rayleigh--B\'enard convection in a slender cylindrical cell, experiments and direct numerical simulations reveal a boundary zonal flow (BZF) that replaces the classical large-scale circulation. The BZF is located near the vertical side wall and enables enhanced heat transport there. Although the azimuthal velocity of the BZF is cyclonic (in the rotating frame), the temperature is an anticyclonic traveling wave of mode one whose signature is a bimodal temperature distribution near the radial boundary. The BZF width is found to scale like $Ra^{1/4}Ek^{2/3}$ where the Ekman number $Ek$ decreases with increasing rotation rate.

A new vertical water tunnel with global temperature control and the possibility for bubble and local heat & mass injection has been designed and constructed. The new facility offers the possibility to accurately study heat and mass transfer in turbulent multiphase flow (gas volume fraction up to $8\%$) with a Reynolds-number range from $1.5 \times 10^4$ to $3 \times 10^5$ in the case of water at room temperature. The tunnel is made of high-grade stainless steel permitting the use of salt solutions in excess of 15$\%$ mass fraction. The tunnel has a volume of 300 liters. The tunnel has three interchangeable measurement sections of $1$ m height but with different cross sections ($0.3 \times 0.04 m^2$, $0.3 \times 0.06 m^2$, $0.3 \times 0.08 m^2$). The glass vertical measurement sections allow for optical access to the flow, enabling techniques such as laser Doppler anemometry, particle image velocimetry, particle tracking velocimetry, and laser-induced fluorescent imaging. Local sensors can be introduced from the top and can be traversed using a built-in traverse system, allowing for e.g. local temperature, hot-wire, or local phase measurements. Combined with simultaneous velocity measurements, the local heat flux in single phase and two phase turbulent flows can thus be studied quantitatvely and precisely.

We present results on the global and local characterisation of heat transport in homogeneous bubbly flow. Experimental measurements were performed with and without the injection of $\sim 2.5$ mm diameter bubbles (corresponding to $Re_b \approx 600$) in a rectangular water column heated from one side and cooled from the other. The gas volume fraction $\alpha$ was varied in the range $0\% - 5\%$, and the Rayleigh number $Ra_H$ in the range $4.0 \times 10^9 - 1.2 \times 10^{11}$. We find that the global heat transfer is enhanced up to 20 times due to bubble injection. Interestingly, for bubbly flow, for our lowest concentration $\alpha = 0.5\% $ onwards, the Nusselt number $\overline{Nu}$ is nearly independent of $Ra_H$, and depends solely on the gas volume fraction~$\alpha$. We observe the scaling $\overline{Nu} \propto \alpha^{0.45}$, which is suggestive of a diffusive transport mechanism. Through local temperature measurements, we show that the bubbles induce a huge increase in the strength of liquid temperature fluctuations, e.g. by a factor of 200 for $\alpha = 0.9\%$. Further, we compare the power spectra of the temperature fluctuations for the single- and two-phase cases. In the single-phase cases, most of the spectral power of the temperature fluctuations is concentrated in the large-scale rolls. However, with the injection of bubbles, we observe intense fluctuations over a wide range of scales, extending up to very high frequencies. Thus, while in the single-phase flow the thermal boundary layers control the heat transport, once the bubbles are injected, the bubble-induced liquid agitation governs the process from a very small bubble concentration onwards., Comment: Journal of Fluid Mechanics (accepted)

Taylor-Couette flow with independently rotating inner (i) and outer (o) cylinders is explored numerically and experimentally to determine the effects of the radius ratio {\eta} on the system response. Numerical simulations reach Reynolds numbers of up to Re_i=9.5 x 10^3 and Re_o=5x10^3, corresponding to Taylor numbers of up to Ta=10^8 for four different radius ratios {\eta}=r_i/r_o between 0.5 and 0.909. The experiments, performed in the Twente Turbulent Taylor-Couette (T^3C) setup, reach Reynolds numbers of up to Re_i=2x10^6$ and Re_o=1.5x10^6, corresponding to Ta=5x10^{12} for {\eta}=0.714-0.909. Effective scaling laws for the torque J^{\omega}(Ta) are found, which for sufficiently large driving Ta are independent of the radius ratio {\eta}. As previously reported for {\eta}=0.714, optimum transport at a non-zero Rossby number Ro=r_i|{\omega}_i-{\omega}_o|/[2(r_o-r_i){\omega}_o] is found in both experiments and numerics. Ro_opt is found to depend on the radius ratio and the driving of the system. At a driving in the range between {Ta\sim3\cdot10^8} and {Ta\sim10^{10}}, Ro_opt saturates to an asymptotic {\eta}-dependent value. Theoretical predictions for the asymptotic value of Ro_{opt} are compared to the experimental results, and found to differ notably. Furthermore, the local angular velocity profiles from experiments and numerics are compared, and a link between a flat bulk profile and optimum transport for all radius ratios is reported., Comment: Submitted to JFM, 28 pages, 17 figures

In the present work it will be shown how the curvature of the outer cylinder affects laser Doppler anemometry measurements inside a Taylor?Couette apparatus. The measurement position and the measured velocity are altered by curved surfaces. Conventional methods for curvature correction are not applicable to our setup, and it will be shown how a ray-tracer can be used to solve this complication. By using a ray-tracer the focal position can be calculated, and the velocity can be corrected. The results of the ray-tracer are verified by measuring an a priori known velocity field, and after applying refractive corrections good agreement with theoretical predictions are found. The methods described in this paper are applied to measure the azimuthal velocity profiles in high Reynolds number Taylor?Couette flow for the case of outer cylinder rotation.

Strongly turbulent Taylor-Couette flow with independently rotating inner and outer cylinders with a radius ratio of \eta = 0.716 is experimentally studied. From global torque measurements, we analyse the dimensionless angular velocity flux Nu_\omega(Ta, a) as a function of the Taylor number Ta and the angular velocity ratio a = -\omega_o/\omega_i in the large-Taylor-number regime 10^{11} \lesssim Ta \lesssim 10^{13}. We analyse the data with the common power-law ansatz for the dimensionless angular velocity transport flux Nu_\omega(Ta, a) = f(a)Ta^\gamma, with an amplitude f(a) and an exponent \gamma. The data are consistent with one effective exponent \gamma = 0.39\pm0.03 for all a. The amplitude of the angular velocity flux f(a) = Nu_\omega(Ta, a)/Ta^0.39 is measured to be maximal at slight counter-rotation, namely at an angular velocity ratio of a_opt = 0.33\pm0.04. This value is theoretically interpreted as the result of a competition between the destabilizing inner cylinder rotation and the stabilizing but shear-enhancing outer cylinder counter-rotation. With the help of laser Doppler anemometry, we provide angular velocity profiles and identify the radial position r_n of the neutral line. While for moderate counter-rotation -0.40 \omega_i \lesssim \omega_o < 0, the neutral line still remains close to the outer cylinder and the probability distribution function (p.d.f.) of the bulk angular velocity is observed to be monomodal. For stronger counter-rotation the neutral line is pushed inwards towards the inner cylinder; in this regime the p.d.f. of the bulk angular velocity becomes bimodal, reflecting intermittent bursts of turbulent structures beyond the neutral line into the outer flow domain, which otherwise is stabilized by the counter-rotating outer cylinder. Finally, a hypothesis is offered allowing a unifying view for all these various results., Comment: 30 pages, 22 figures

Three-dimensional particle tracking velocimetry (PTV) and phase-sensitive constant temperature anemometry in pseudo-turbulence – i.e. flow solely driven by rising bubbles – were performed to investigate bubble clustering and to obtain the mean bubble rise velocity, distributions of bubble velocities and energy spectra at dilute gas concentrations (α ≤ 2.2 %). To characterize the clustering the pair correlation function G(r, θ) was calculated. The deformable bubbles with equivalent bubble diameter db = 4–5 mm were found to cluster within a radial distance of a few bubble radii with a preferred vertical orientation. This vertical alignment was present at both small and large scales. For small distances also some horizontal clustering was found. The large number of data points and the non-intrusiveness of PTV allowed well-converged probability density functions (PDFs) of the bubble velocity to be obtained. The PDFs had a non-Gaussian form for all velocity components and intermittency effects could be observed. The energy spectrum of the liquid velocity fluctuations decayed with a power law of −3.2, different from the ≈ −5/3 found for homogeneous isotropic turbulence, but close to the prediction −3 by Lance & Bataille (J. Fluid Mech., vol. 222, 1991, p. 95) for pseudo-turbulence.

We report measurements of the temperature variance σ(2)(z,r) and frequency power spectrum P(f,z,r) (z is the distance from the sample bottom and r the radial coordinate) in turbulent Rayleigh-Bénard convection (RBC) for Rayleigh numbers Ra = 1.6 × 10(13) and 1.1 × 10(15) and for a Prandtl number Pr ≃ 0.8 for a sample with a height L = 224 cm and aspect ratio D/L=0.50 (D is the diameter). For z/L ≲ 0.1 σ(2)(z,r) was consistent with a logarithmic dependence on z, and there was a universal (independent of Ra, r, and z) normalized spectrum which, for 0.02 ≲ fτ(0) ≲ 0.2, had the form P(fτ(0)) = P(0)(fτ(0))(-1) with P(0) = 0.208 ± 0.008 a universal constant. Here τ(0) = sqrt[2R] where R is the radius of curvature of the temperature autocorrelation function C(τ) at τ = 0. For z/L ≃ 0.5 the measurements yielded P(fτ(0))∼(fτ(0))(-α) with α in the range from 3/2 to 5/3. All the results are similar to those for velocity fluctuations in shear flows at sufficiently large Reynolds numbers, suggesting the possibility of an analogy between the flows that is yet to be determined in detail.

The research issues addressed in this mostly experimental thesis concern highly turbulent Taylor-Couette (TC) flow (Re>105, implying Ta>1011). We study it on a fundamental level to aid our understanding of (TC) turbulence and to make predictions towards astrophysical disks, and at a practical level as applications can be found in bubble-induced skin-friction drag reduction on ships. In PART I we introduce the new TC facility of our Physics of Fluids group, called the Twente turbulent Taylor-Couette (T3C) facility. It features two independently rotating cylinders of variable radius ratio with accurate rotation rate and temperature control, torque sensing, bubble injection and it is equipped with several local sensors. It is able to reach Reynolds numbers up to 3:4�106. In PART II we focus on highly turbulent single-phase TC flow. We measure the global torque as a function of the driving parameters and we provide local angular velocity measurements. The results are interpreted as the transport of angular velocity, based on the model proposed by Eckhardt, Grossmann & Lohse (2007). Furthermore, we study the turbulence transport in quasi-Keplerian profiles, mimicking astrophysical disks. In PART III we study the effect of bubbles on highly turbulent TC flow, focusing not only on global drag reduction, but also on the local bubble distribution and angular velocity profiles. We find that drag reduction is a boundary layer effect and that the deformability of bubbles is crucial for strong drag reduction in bubbly turbulent TC flow.

The flow structure of strongly turbulent Taylor-Couette flow with Reynolds numbers up to Re(i)=2×10(6) of the inner cylinder is experimentally examined with high-speed particle image velocimetry (PIV). The wind Reynolds numbers Re(w) of the turbulent Taylor-vortex flow is found to scale as Re(w)∝Ta(1/2), exactly as predicted by Grossmann and Lohse [Phys. Fluids 23, 045108 (2011).] for the ultimate turbulence regime, in which the boundary layers are turbulent. The dimensionless angular velocity flux has an effective scaling of Nu(ω)∝Ta(0.38), also in correspondence with turbulence in the ultimate regime. The scaling of Nu(ω) is confirmed by local angular velocity flux measurements extracted from high-speed PIV measurements: though the flux shows huge fluctuations, its spatial and temporal average nicely agrees with the result from the global torque measurements.

A new turbulent Taylor-Couette system consisting of two independently rotating cylinders has been constructed. The gap between the cylinders has a height of 0.927 m, an inner radius of 0.200 m, and a variable outer radius (from 0.279 to 0.220 m). The maximum angular rotation rates of the inner and outer cylinder are 20 and 10 Hz, respectively, resulting in Reynolds numbers up to 3.4 x 10^6 with water as working fluid. With this Taylor-Couette system, the parameter space (Re_i, Re_o, {\eta}) extends to (2.0 x 10^6, {\pm}1.4 x 10^6, 0.716-0.909). The system is equipped with bubble injectors, temperature control, skin-friction drag sensors, and several local sensors for studying turbulent single-phase and two-phase flows. Inner cylinder load cells detect skin-friction drag via torque measurements. The clear acrylic outer cylinder allows the dynamics of the liquid flow and the dispersed phase (bubbles, particles, fibers, etc.) inside the gap to be investigated with specialized local sensors and nonintrusive optical imaging techniques. The system allows study of both Taylor-Couette flow in a high-Reynolds-number regime, and the mechanisms behind skin-friction drag alterations due to bubble injection, polymer injection, and surface hydrophobicity and roughness., Comment: 13 pages, 14 figures

We analyze the global transport properties of turbulent Taylor-Couette flow in the strongly turbulent regime for independently rotating outer and inner cylinders, reaching Reynolds numbers of the inner and outer cylinders of Re_i = 2 x 10^6 and Re_o = 1.4 x 10^6, respectively. For all Re_i, Re_o, the dimensionless torque G scales as a function of the Taylor number Ta (which is proportional to the square of the difference between the angular velocities of the inner and outer cylinders) with a universal effective scaling law G \propto Ta^{0.88}, corresponding to Nu_omega \propto Ta^{0.38} for the Nusselt number characterizing the angular velocity transport between the inner and outer cylinders. The exponent 0.38 corresponds to the ultimate regime scaling for the analogous Rayleigh-Benard system. The transport is most efficient for the counterrotating case along the diagonal in phase space with omega_o \approx -0.4 omega_i., 3 pages, 4 figures

In turbulent Taylor-Couette flow, the injection of bubbles reduces the overall drag. On the other hand, rough walls enhance the overall drag. In this work, we inject bubbles into turbulent Taylor-Couette flow with rough walls (with a Reynolds number up to 4 x 10(5), finding an enhancement of the dimensionless drag as compared to the case without bubbles. The dimensional drag is unchanged. As in the rough-wall case no smooth boundary layers can develop, the results demonstrate that bubbly drag reduction is a pure boundary layer effect.