Source author record

Sander G. Huisman

Sander G. Huisman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2026arXiv

Settling of chiral particles in a turbulent flow

Chiral particles are experimentally investigated while settling inwater with various turbulence intensity levels. The locations and orientations of the particles are tracked over time, allowing the close investigation of the particles' settling dynamics. The generated turbulent flow is measured using laser Doppler anemometry (LDA), and the turbulence strength varies between experiments in the range $0 \leq Re_λ\leq 250$. Starting with quiescent particle settling, the chiral particle's orientation dynamics are studied, revealing a preferred alignment and a strong translation-rotation coupling. The particle chirality determines the preferred rotation direction, though the alignment and translation-rotation coupling gradually vanish with increasing turbulence. We identify multiple settling modes for the chiral particles, which are characterised by the evolution of the rotation angles. Finally, a theoretical model assuming a simplified chiral particle in Stokes flow clarifies the emergence of each settling mode.

preprint2026arXiv

Transition from classical to ultimate melting

Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulence to ultimate turbulence, with its enhanced transport properties. Our results thus provide a quantitative understanding of the flow physics of the melting process and thereby enable a better extrapolation and prediction of melt rates on large scales such as relevant in geophysics, oceanography, and climate science.

preprint2025arXiv

Collective effects of neighbouring melting ice objects

We present a study on the melting dynamics of neighbouring ice bodies by means of idealised simulations, focusing on collective effects, with the goal of obtaining fundamental insight into how collective interactions influence the melting of ice. Two neighbouring (vertically or horizontally aligned), square-shaped, and equally sized ice objects (size on the order of centimetres) are immersed in quiescent fresh water at a temperature of 20°C. By performing two-dimensional direct numerical simulations, and using the phase-field method to model the phase change, the collective melting of these objects is studied. When the objects are horizontally aligned, no significant influence of the neighbouring object on the melting time is observed. On the other hand, when vertically aligned, though the melting of the upper object is mostly unaffected, the melting time and the morphology of the lower ice body strongly depends on the initial inter-object distance. We report that the melting of the bottom object can be enhanced by more than 10%, or delayed more than 20%, displaying a non-monotonic dependence on the initial object size. We show that this behaviour results from a non-trivial competition between layering of cold fluid, which lowers the heat transfer, and convective flows, which favour mixing and heat transfer. For this melting in mixed convection, we were able to collapse our data onto a single curve.

preprint2022arXiv

Micro-droplet nucleation through solvent exchange in a turbulent buoyant jet

Solvent exchange is a process involving mixing between a good solvent with dissolved solute and a poor solvent. The process creates local oversaturation which causes the nucleation of minute solute droplets. Such ternary systems on a macro-scale have remained unexplored in the turbulent regime. We experimentally study the solvent exchange process by injecting mixtures of ethanol and trans-anethole into water, forming a turbulent buoyant jet in upward direction. Locally, turbulent mixing causes oversaturation of the trans-anethole following turbulent entrainment. We optically measure the concentration of the nucleated droplets using a light attenuation technique and find that the radial concentration profile is sub-Gaussian. In contrast to the entrainment-based models, the spatial evolution of the oversaturation reveals continuous droplet nucleation downstream and radially across the jet, which we attribute to the limited mixing capacity of the jet. Though we are far from a full quantitative understanding, this work extends the knowledge on solvent exchange into the turbulent regime, and brings in a novel type of flow, broadening the scope of multicomponent, multiphase turbulent jets with phase transition.

preprint2022arXiv

Recent developments of turbulent emulsions in Taylor-Couette flow

Emulsions are common in many natural and industrial settings. Recently, much attention has been put on understanding the dynamics of turbulent emulsions. This paper reviews some recent studies of emulsions in turbulent Taylor-Couette flow, mainly focusing on the statistics of the dispersed phase and the global momentum transport of the system. We first study the size distribution and the breakup mechanism of the dispersed droplets for turbulent emulsions with a low volume-fraction (dilute) of the dispersed phase. For systems with a high volume-fraction of the dispersed phase (dense), we address the detailed response of the global transport (effective viscosity) of the turbulent emulsion and its connection to the droplet statistics. Finally, we will discuss catastrophic phase inversions, which can happen when the volume fraction of the dispersed phase exceeds a critical value during dynamic emulsification. We end the manuscript with a summary and an outlook including some open questions for future research. This article is part of the theme issue `Taylor-Couette and Related Flows on the Centennial of Taylor's Seminal Philosophical Transactions Paper'.

preprint2022arXiv

Supergravitational turbulent thermal convection

High-Rayleigh number convective turbulence is ubiquitous in many natural phenomena and in industries, such as atmospheric circulations, oceanic flows, flows in the fluid core of planets, and energy generations. In this work, we present a novel approach to boost the Rayleigh number in thermal convection by exploiting centrifugal acceleration and rapidly rotating a cylindrical annulus to reach an effective gravity of 60 times Earth's gravity. We show that in the regime where the Coriolis effect is strong, the scaling exponent of Nusselt number versus Rayleigh number exceeds one-third once the Rayleigh number is large enough. The convective rolls revolve in prograde direction, signifying the emergence of zonal flow. The present findings open a new avenue on the exploration of high-Rayleigh number turbulent thermal convection and will improve the understanding of the flow dynamics and heat transfer processes in geophysical and astrophysical flows and other strongly rotating systems.

preprint2021arXiv

Broken mirror symmetry of tracer's trajectories in turbulence

Topological properties of physical systems play a crucial role in our understanding of nature, yet their experimental determination remains elusive. We show that the mean helicity, a dynamical invariant in ideal flows, quantitatively affects trajectories of fluid elements: the linking number of Lagrangian trajectories depends on the mean helicity. Thus, a global topological invariant and a topological number of fluid trajectories become related, and we provide an empirical expression linking them. The relation shows the existence of long-term memory in the trajectories: the links can be made of the trajectory up to a given time, with particles positions in the past. This property also allows experimental measurements of mean helicity.

preprint2020arXiv

Bubbly drag reduction using a hydrophobic inner cylinder in Taylor-Couette turbulence

In this study we experimentally investigate bubbly drag reduction in a highly turbulent flow of water with dispersed air at $5.0 \times 10^{5} \leq \text{Re} \leq 1.7 \times 10^{6}$ over a non-wetting surface containing micro-scale roughness. To do so, the Taylor-Couette geometry is used, allowing for both accurate global drag and local flow measurements. The inner cylinder - coated with a rough, hydrophobic material - is rotating, whereas the smooth outer cylinder is kept stationary. The crucial control parameter is the air volume fraction $α$ present in the working fluid. For small volume fractions ($α< {4}\,\%$), we observe that the surface roughness from the coating increases the drag. For large volume fractions of air ($α\geq 4\,\%$), the drag decreases compared to the case with both the inner and outer cylinders uncoated, i.e. smooth and hydrophilic, using the same volume fraction of air. This suggests that two competing mechanisms are at place: on the one hand the roughness invokes an extension of the log-layer - resulting in an increase in drag - and on the other hand there is a drag-reducing mechanism of the hydrophobic surface interacting with the bubbly liquid. The balance between these two effects determines whether there is overall drag reduction or drag enhancement. For further increased bubble concentration $α= {6}\,\%$ we find a saturation of the drag reduction effect. Our study gives guidelines for industrial applications of bubbly drag reduction in hydrophobic wall-bounded turbulent flows.

preprint2020arXiv

Catastrophic phase inversion in high-Reynolds number turbulent Taylor--Couette flow

Emulsions are omnipresent in the food industry, health care, and chemical synthesis. In this Letter the dynamics of meta-stable oil-water emulsions in highly turbulent ($10^{11}\leq\text{Ta}\leq 3\times 10^{13}$) Taylor--Couette flow, far from equilibrium, is investigated. By varying the oil-in-water void fraction, catastrophic phase inversion between oil-in-water and water-in-oil emulsions can be triggered, changing the morphology, including droplet sizes, and rheological properties of the mixture, dramatically. The manifestation of these different states is exemplified by combining global torque measurements and local in-situ laser induced fluorescence (LIF) microscopy imaging. Despite the turbulent state of the flow and the dynamic equilibrium of the oil-water mixture, the global torque response of the system is found to be as if the fluid were Newtonian, and the effective viscosity of the mixture was found to be several times bigger or smaller than either of its constituents.

preprint2020arXiv

Double maxima of angular momentum transport in $η=0.91$ TC turbulence

We use experiments and direct numerical simulations to probe the phase-space of low-curvature Taylor--Couette (TC) flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $η=r_i/r_o=0.91$. Non-dimensional shear drivings (Taylor numbers $\text{Ta}$) in the range $10^7\leq\text{Ta}\leq10^{11}$ are explored for both co- and counter-rotating configurations. In the $\text{Ta}$ range $10^8\leq\text{Ta}\leq10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either as "co-" and "counter-rotating" due to their location or as "broad" and "narrow" due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann \emph{et al.}~(2016), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with $\text{Ta}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around $\text{Ta}\approx10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

preprint2020arXiv

Effect of axially varying sandpaper roughness on bubbly drag reduction in Taylor-Couette turbulence

We experimentally investigate the influence of alternating rough and smooth walls on bubbly drag reduction (DR). We apply rough sandpaper bands of width $s$ between $48.4\,mm$ and $148.5\,mm$, and roughness height $k = 695\,μm$, around the smooth inner cylinder (IC) of the Twente Turbulent Taylor-Couette facility. Between sandpaper bands, the IC is left uncovered over similar width $s$, resulting in alternating rough and smooth bands, a constant pattern in axial direction. We measure the DR in water that originates from introducing air bubbles to the fluid at (shear) Reynolds numbers $\textit{Re}_s$ ranging from $0.5 \times 10^6$ to $1.8 \times 10^6$. Results are compared to bubbly DR measurements with a completely smooth IC and an IC that is completely covered with sandpaper of the same roughness $k$. The outer cylinder is left smooth for all variations. Results are also compared to bubbly DR measurements where a smooth outer cylinder is rotating in opposite direction to the smooth IC. This counter rotation induces secondary flow structures that are very similar to those observed when the IC is composed of alternating rough and smooth bands. For the measurements with roughness, the bubbly DR is found to initially increase more strongly with $\textit{Re}_s$, before levelling off to reach a value that no longer depends on $\textit{Re}_s$. This is attributed to a more even axial distribution of the air bubbles, resulting from the increased turbulence intensity of the flow compared to flow over a completely smooth wall at the same $\textit{Re}_s$. The air bubbles are seen to accumulate at the rough wall sections in the flow. Here, locally, the drag is largest and so the drag reducing effect of the bubbles is felt strongest. Therefore, a larger maximum value of bubbly DR is found for the alternating rough and smooth walls compared to the completely rough wall.

preprint2020arXiv

Heat transfer enhancement in Rayleigh-Bénard convection using a single passive barrier

In this numerical study on Rayleigh-Bénard convection we seek to improve the heat transfer by passive means. To this end we introduce a single tilted conductive barrier centered in an aspect ratio one cell, breaking the symmetry of the geometry and to channel the ascending hot and descending cold plumes. We study the global and local heat transfer and the flow organization for Rayleigh numbers $10^5 \leq Ra \leq 10^9$ for a fixed Prandtl number of $Pr=4.3$. We find that the global heat transfer can be enhanced up to $18\%$, and locally around $800\%$. The averaged Reynolds number is always decreased when a barrier is introduced, even for those cases where the global heat transfer is increased. We map the entire parameter space spanned by the orientation and the size of a single barrier for $Ra=10^8$.

preprint2019arXiv

Controlling secondary flow in Taylor-Couette turbulence through spanwise-varying roughness

Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $η=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices.

preprint2019arXiv

Drag reduction in boiling Taylor-Couette turbulence

We create a highly controlled lab environment-accessible to both global and local monitoring-to analyse turbulent boiling flows and in particular their shear stress in a statistically stationary state. Namely, by precisely monitoring the drag of strongly turbulent Taylor-Couette flow (the flow in between two co-axially rotating cylinders, Reynolds number $\textrm{Re}\approx 10^6$) during its transition from non-boiling to boiling, we show that the intuitive expectation, namely that a few volume percent of vapor bubbles would correspondingly change the global drag by a few percent, is wrong. Rather, we find that for these conditions a dramatic global drag reduction of up to 45% occurs. We connect this global result to our local observations, showing that for major drag reduction the vapor bubble deformability is crucial, corresponding to Weber numbers larger than one. We compare our findings with those for turbulent flows with gas bubbles, which obey very different physics than vapor bubbles. Nonetheless, we find remarkable similarities and explain these.

preprint2016arXiv

Columnar structure formation of a dilute suspension of settling spherical particles in a quiescent fluid

The settling of heavy spherical particles in a column of quiescent fluid is investigated. The performed experiments cover a range of Galileo numbers ($110 \leq \text{Ga} \leq 310$) for a fixed density ratio of $Γ= ρ_p/ρ_f = 2.5$. In this regime the particles are known (M. Jenny, J. Dušek and G. Bouchet, Journal of Fluid Mechanics 508, 201 (2004).) to show a variety of motions. It is known that the wake undergoes several transitions for increasing $\text{Ga}$ resulting in particle motions that are successively: vertical, oblique, oblique oscillating, and finally chaotic. Not only does this change the trajectory of single, isolated, settling particles, but it also changes the dynamics of a swarm of particles as collective effects become important even for dilute suspensions, with volume fraction up to $Φ_V = \mathcal{O}\left(10^{-3}\right)$, which are investigated in this work. Multi-camera recordings of settling particles are recorded and tracked over time in 3 dimensions. A variety of analysis are performed and show a strong clustering behavior. The distribution of the cell areas of the Voronoï tessellation in the horizontal plane are compared to that of a random distribution of particles and shows clear clustering. Moreover, a negative correlation was found between the Voronoï area and the particle velocity; clustered particles fall faster. In addition, the angle between two adjacent particles and the vertical is calculated and compared to a homogeneous distribution of particles, clear evidence of vertical alignment of particles is found. The experimental findings are compared to simulations.

preprint2016arXiv

Exploring the phase space of multiple states in highly turbulent Taylor-Couette flow

We investigate the existence of multiple turbulent states in highly turbulent Taylor-Couette flow in the range of $\mathrm{Ta}=10^{11}$ to $9\cdot10^{12}$, by measuring the global torques and the local velocities while probing the phase space spanned by the rotation rates of the inner and outer cylinder. The multiple states are found to be very robust and are expected to persist beyond $\mathrm{Ta}=10^{13}$. The rotation ratio is the parameter that most strongly controls the transitions between the flow states; the transitional values only weakly depend on the Taylor number. However, complex paths in the phase space are necessary to unlock the full region of multiple states. Lastly, by mapping the flow structures for various rotation ratios in a Taylor-Couette setup with an equal radius ratio but a larger aspect ratio than before, multiple states were again observed. Here, they are characterized by even richer roll structure phenomena, including, for the first time observed in highly turbulent TC flow, an antisymmetrical roll state.

preprint2016arXiv

Multiple states in highly turbulent Taylor-Couette flow

The ubiquity of turbulent flows in nature and technology makes it of utmost importance to fundamentally understand turbulence. Kolmogorov's 1941 paradigm suggests that for strongly turbulent flows with many degrees of freedom and its large fluctuations, there would only be \emph{one} turbulent state as the large fluctuations would explore the entire higher-dimensional phase space. Here we report the first conclusive evidence of multiple turbulent states for large Reynolds number $\text{Re}=\mathcal{O}(10^6)$ (Taylor number $\text{Ta}=\mathcal{O}(10^{12})$) Taylor-Couette flow in the regime of ultimate turbulence, by probing the phase space spanned by the rotation rates of the inner and outer cylinder. The manifestation of multiple turbulent states is exemplified by providing combined global torque and local velocity measurements. This result verifies the notion that bifurcations can occur in high-dimensional flows i.e. very large $\text{Re}$) and questions Kolmogorov's paradigm.

preprint2016arXiv

The boiling Twente Taylor-Couette (BTTC) facility: Temperature controlled turbulent flow between independently rotating, coaxial cylinders

A new Taylor-Couette system has been designed and constructed with precise temperature control. Two concentric independently rotating cylinders are able to rotate at maximum rates of $f_i = \pm20$ Hz for the inner cylinder and $f_o = \pm10$ Hz for the outer cylinder. The inner cylinder has an outside radius of $r_i = 75$ mm, and the outer cylinder has an inside radius of $r_o = 105$ mm, resulting in a gap of $d=30$ mm. The height of the gap $L =549$ mm, giving a volume of $V=9.3$l. The geometric parameters are $η= r_i/r_o = 0.714$ and $Γ= L/d = 18.3$. With water as working fluid at room temperature the Reynolds numbers that can be achieved are $\text{Re}_i = ω_i r_i (r_o-r_i)/ν= 2.8 \times 10^5$ and $\text{Re}_o = ω_o r_o (r_o-r_i)/ν= 2 \times 10^5$, or a combined Reynolds number of up to $\text{Re} = (ω_i r_i -ω_o r_o)(r_o-r_i)/ν= 4.8 \times 10^5$. If the working fluid is changed to the fluorinated liquid FC-3284 with kinematic viscosity 0.42 cSt, the combined Reynolds number can reach $\text{Re} = 1.1 \times 10^6$. The apparatus features precise temperature control of the outer and inner cylinder separately, and is fully optically accessible from the side and top. The new facility offers the possibility to accurately study the process of boiling inside a turbulent flow and its effect on the flow.

preprint2015arXiv

Azimuthal velocity profiles in Rayleigh-stable Taylor-Couette flow and implied axial angular momentum transport

We present azimuthal velocity profiles measured in a Taylor-Couette apparatus, which has been used as a model of stellar and planetary accretion disks. The apparatus has a cylinder radius ratio of $η= 0.716$, an aspect-ratio of $Γ= 11.74$, and the plates closing the cylinders in the axial direction are attached to the outer cylinder. We investigate angular momentum transport and Ekman pumping in the Rayleigh-stable regime. The regime is linearly stable and is characterized by radially increasing specific angular momentum. We present several Rayleigh-stable profiles for shear Reynolds numbers $Re_S \sim O(10^5) \,$, both for $Ω_i > Ω_o > 0$ (quasi-Keplerian regime) and $Ω_o > Ω_i > 0$ (sub-rotating regime) where $Ω_{i,o}$ is the inner/outer cylinder rotation rate. None of the velocity profiles matches the non-vortical laminar Taylor-Couette profile. The deviation from that profile increased as solid-body rotation is approached at fixed $Re_S$. Flow super-rotation, an angular velocity greater than that of both cylinders, is observed in the sub-rotating regime. The velocity profiles give lower bounds for the torques required to rotate the inner cylinder that were larger than the torques for the case of laminar Taylor-Couette flow. The quasi-Keplerian profiles are composed of a well mixed inner region, having approximately constant angular momentum, connected to an outer region in solid-body rotation with the outer cylinder and attached axial boundaries. These regions suggest that the angular momentum is transported axially to the axial boundaries. Therefore, Taylor-Couette flow with closing plates attached to the outer cylinder is an imperfect model for accretion disk flows, especially with regard to their stability.

preprint2015arXiv

Taylor-Couette turbulence at radius ratio $η=0.5$: scaling, flow structures and plumes

Using high-resolution particle image velocimetry we measure velocity profiles, the wind Reynolds number and characteristics of turbulent plumes in Taylor-Couette flow for a radius ratio of 0.5 and Taylor number of up to $6.2\cdot10^9$. The extracted angular velocity profiles follow a log-law more closely than the azimuthal velocity profiles due to the strong curvature of this $η=0.5$ setup. The scaling of the wind Reynolds number with the Taylor number agrees with the theoretically predicted 3/7-scaling for the classical turbulent regime, which is much more pronounced than for the well-explored $η=0.71$ case, for which the ultimate regime sets in at much lower Ta. By measuring at varying axial positions, roll structures are found for counter-rotation while no clear coherent structures are seen for pure inner cylinder rotation. In addition, turbulent plumes coming from the inner and outer cylinder are investigated. For pure inner cylinder rotation, the plumes in the radial velocity move away from the inner cylinder, while the plumes in the azimuthal velocity mainly move away from the outer cylinder. For counter-rotation, the mean radial flow in the roll structures strongly affects the direction and intensity of the turbulent plumes. Furthermore, it is experimentally confirmed that in regions where plumes are emitted, boundary layer profiles with a logarithmic signature are created.

preprint2013arXiv

Logarithmic boundary layers in highly turbulent Taylor-Couette flow

We provide direct measurements of the boundary layer properties in highly turbulent Taylor-Couette flow up to $\text{Ta}=6.2 \times 10^{12}$ using high-resolution particle image velocimetry (PIV). We find that the mean azimuthal velocity profile at the inner and outer cylinder can be fitted by the von Kármán log law $u^+ = \frac 1κ\ln y^+ +B$. The von Kármán constant $κ$ is found to depend on the driving strength $\text{Ta}$ and for large $\text{Ta}$ asymptotically approaches $κ\approx 0.40$. The variance profiles of the local azimuthal velocity have a universal peak around $y^+ \approx 12$ and collapse when rescaled with the driving velocity (and not with the friction velocity), displaying a log-dependence of $y^+$ as also found for channel and pipe flows [1,2].

preprint2013arXiv

Statistics of turbulent fluctuations in counter-rotating Taylor-Couette flows

The statistics of velocity fluctuations of turbulent Taylor-Couette flow are examined. The rotation rate of the inner and outer cylinder are varied while keeping the Taylor number fixed to $1.49 \times 10^{12}$ ($\mathcal{O}(\text{Re})=10^6$). The azimuthal velocity component of the flow is measured using laser Doppler anemometry (LDA). For each experiment $5\times10^6$ datapoints are acquired and carefully analysed. Using extended self-similarity (ESS) \cite{ben93b} the longitudinal structure function exponents are extracted, and are found to weakly depend on the ratio of the rotation rates. For the case where only the inner cylinder rotates the results are in good agreement with results measured by Lewis and Swinney \cite{lew99} using hot-film anemometry. The power spectra shows clear -5/3 scaling for the intermediate angular velocity ratios $-ω_o/ω_i \in \{0.6, 0.8, 1.0\}$, roughly -5/3 scaling for $-ω_o/ω_i \in \{0.2, 0.3, 0.4, 2.0\}$, and no clear scaling law can be found for $-ω_0/ω_i = 0$ (inner cylinder rotation only); the local scaling exponent of the spectra has a strong frequency dependence. We relate these observations to the shape of the probability density function of the azimuthal velocity and the presence of a neutral line.

preprint2012arXiv

Applying Laser Doppler Anemometry inside a Taylor-Couette geometry - Using a ray-tracer to correct for curvature effects

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.

preprint2012arXiv

Optimal Taylor-Couette turbulence

Strongly turbulent Taylor-Couette flow with independently rotating inner and outer cylinders with a radius ratio of η= 0.716 is experimentally studied. From global torque measurements, we analyse the dimensionless angular velocity flux Nu_ω(Ta, a) as a function of the Taylor number Ta and the angular velocity ratio a = -ω_o/ω_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_ω(Ta, a) = f(a)Ta^γ, with an amplitude f(a) and an exponent γ. The data are consistent with one effective exponent γ= 0.39\pm0.03 for all a. The amplitude of the angular velocity flux f(a) = Nu_ω(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 ω_i \lesssim ω_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.

preprint2011arXiv

Torque Scaling in Turbulent Taylor-Couette Flow with Co- and Counterrotating Cylinders

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.

preprint2011arXiv

Ultimate Turbulent Taylor-Couette 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 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.