scholarly journals Optimal Taylor–Couette flow: direct numerical simulations

2013 ◽  
Vol 719 ◽  
pp. 14-46 ◽  
Author(s):  
Rodolfo Ostilla ◽  
Richard J. A. M. Stevens ◽  
Siegfried Grossmann ◽  
Roberto Verzicco ◽  
Detlef Lohse

AbstractWe numerically simulate turbulent Taylor–Couette flow for independently rotating inner and outer cylinders, focusing on the analogy with turbulent Rayleigh–Bénard flow. Reynolds numbers of $R{e}_{i} = 8\times 1{0}^{3} $ and $R{e}_{o} = \pm 4\times 1{0}^{3} $ of the inner and outer cylinders, respectively, are reached, corresponding to Taylor numbers $Ta$ up to $1{0}^{8} $. Effective scaling laws for the torque and other system responses are found. Recent experiments with the Twente Turbulent Taylor–Couette (${T}^{3} C$) setup and with a similar facility in Maryland at very high Reynolds numbers have revealed an optimum transport at a certain non-zero rotation rate ratio $a= - {\omega }_{o} / {\omega }_{i} $ of about ${a}_{\mathit{opt}} = 0. 33$. For large enough $Ta$ in the numerically accessible range we also find such an optimum transport at non-zero counter-rotation. The position of this maximum is found to shift with the driving, reaching a maximum of ${a}_{\mathit{opt}} = 0. 15$ for $Ta= 2. 5\times 1{0}^{7} $. An explanation for this shift is elucidated, consistent with the experimental result that ${a}_{\mathit{opt}} $ becomes approximately independent of the driving strength for large enough Reynolds numbers. We furthermore numerically calculate the angular velocity profiles and visualize the different flow structures for the various regimes. By writing the equations in a frame co-rotating with the outer cylinder a link is found between the local angular velocity profiles and the global transport quantities.

2015 ◽  
Vol 774 ◽  
pp. 342-362 ◽  
Author(s):  
Freja Nordsiek ◽  
Sander G. Huisman ◽  
Roeland C. A. van der Veen ◽  
Chao Sun ◽  
Detlef Lohse ◽  
...  

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 ${\it\eta}=0.716$, an aspect ratio of ${\it\Gamma}=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. This regime is linearly stable and is characterized by radially increasing specific angular momentum. We present several Rayleigh-stable profiles for shear Reynolds numbers $\mathit{Re}_{S}\sim O(10^{5})$, for both ${\it\Omega}_{i}>{\it\Omega}_{o}>0$ (quasi-Keplerian regime) and ${\it\Omega}_{o}>{\it\Omega}_{i}>0$ (sub-rotating regime), where ${\it\Omega}_{i,o}$ is the inner/outer cylinder rotation rate. None of the velocity profiles match the non-vortical laminar Taylor–Couette profile. The deviation from that profile increases as solid-body rotation is approached at fixed $\mathit{Re}_{S}$. Flow super-rotation, an angular velocity greater than those 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 are 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.


2014 ◽  
Vol 747 ◽  
pp. 1-29 ◽  
Author(s):  
Rodolfo Ostilla-Mónico ◽  
Sander G. Huisman ◽  
Tim J. G. Jannink ◽  
Dennis P. M. Van Gils ◽  
Roberto Verzicco ◽  
...  

AbstractTaylor–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 $\mathit{Re}_i=9.5\times 10^3$ and $\mathit{Re}_o=5\times 10^3$, corresponding to Taylor numbers of up to $\mathit{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 ($\mathrm{T^3C}$) set-up, reach Reynolds numbers of up to $\mathit{Re}_i=2\times 10^6$ and $\mathit{Re}_o=1.5\times 10^6$, corresponding to $\mathit{Ta}=5\times 10^{12}$ for $\eta =0.714\mbox{--}0.909$. Effective scaling laws for the torque $J^{\omega }(\mathit{Ta})$ are found, which for sufficiently large driving $\mathit{Ta}$ are independent of the radius ratio $\eta $. As previously reported for $\eta =0.714$, optimum transport at a non-zero Rossby number $\mathit{Ro}=r_i |\omega _i-\omega _o |/[2(r_o-r_i)\omega _o]$ is found in both experiments and numerics. Here $\mathit{Ro}_{opt}$ is found to depend on the radius ratio and the driving of the system. At a driving in the range between $\mathit{Ta}\sim 3\times 10^{8}$ and $\mathit{Ta}\sim 10^{10}$, $\mathit{Ro}_{opt}$ saturates to an asymptotic $\eta $-dependent value. Theoretical predictions for the asymptotic value of $\mathit{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.


2014 ◽  
Vol 750 ◽  
pp. 555-577 ◽  
Author(s):  
Stephanie Welsh ◽  
Evy Kersalé ◽  
Chris A. Jones

AbstractTaylor–Couette flow in a compressible perfect gas is studied. The onset of instability is examined as a function of the Reynolds numbers of the inner and outer cylinder, the Mach number of the flow and the Prandtl number of the gas. We focus on the case of a wide gap, with radius ratio 0.5. We find new modes of instability at high Prandtl number, which can allow oscillatory axisymmetric modes to onset first. We also find that onset can occur even when the angular momentum increases outwards, so that the classical Rayleigh criterion can be violated in the compressible case. We have also considered the case of counter-rotating cylinders, where the $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}m=0$ and $m=1$ modes can onset simultaneously to give a codimension 2 bifurcation, leading to the formation of complex flow patterns. In compressible flow we also find codimension 3 points. The Mach number and the critical inner and outer Reynolds numbers can be adjusted so that the two neutral curves for the $m=0$ and $m=1$ modes touch rather than cross. Complex codimension 3 points occur more readily in the compressible case than in the Boussinesq case, and they are expected to lead to a rich nonlinear behaviour.


2012 ◽  
Vol 706 ◽  
pp. 118-149 ◽  
Author(s):  
Dennis P. M. van Gils ◽  
Sander G. Huisman ◽  
Siegfried Grossmann ◽  
Chao Sun ◽  
Detlef Lohse

AbstractStrongly 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 ${\mathit{Nu}}_{\omega } (\mathit{Ta}, a)$ as a function of the Taylor number $\mathit{Ta}$ and the angular velocity ratio $a= \ensuremath{-} {\omega }_{o} / {\omega }_{i} $ in the large-Taylor-number regime $1{0}^{11} \lesssim \mathit{Ta}\lesssim 1{0}^{13} $ and well off the inviscid stability borders (Rayleigh lines) $a= \ensuremath{-} {\eta }^{2} $ for co-rotation and $a= \infty $ for counter-rotation. We analyse the data with the common power-law ansatz for the dimensionless angular velocity transport flux ${\mathit{Nu}}_{\omega } (\mathit{Ta}, a)= f(a)\hspace{0.167em} {\mathit{Ta}}^{\gamma } $, with an amplitude $f(a)$ and an exponent $\gamma $. The data are consistent with one effective exponent $\gamma = 0. 39\pm 0. 03$ for all $a$, but we discuss a possible $a$ dependence in the co- and weakly counter-rotating regimes. The amplitude of the angular velocity flux $f(a)\equiv {\mathit{Nu}}_{\omega } (\mathit{Ta}, a)/ {\mathit{Ta}}^{0. 39} $ is measured to be maximal at slight counter-rotation, namely at an angular velocity ratio of ${a}_{\mathit{opt}} = 0. 33\pm 0. 04$, i.e. along the line ${\omega }_{o} = \ensuremath{-} 0. 33{\omega }_{i} $. 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 in particular identify the radial position ${r}_{n} $ of the neutral line, defined by $ \mathop{ \langle \omega ({r}_{n} )\rangle } \nolimits _{t} = 0$ for fixed height $z$. For these large $\mathit{Ta}$ values, the ratio $a\approx 0. 40$, which is close to ${a}_{\mathit{opt}} = 0. 33$, is distinguished by a zero angular velocity gradient $\partial \omega / \partial r= 0$ in the bulk. While for moderate counter-rotation $\ensuremath{-} 0. 40{\omega }_{i} \lesssim {\omega }_{o} \lt 0$, the neutral line still remains close to the outer cylinder and the probability distribution function 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 probability distribution function 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 and consistent interpretation for all these various results.


2017 ◽  
Vol 831 ◽  
pp. 330-357 ◽  
Author(s):  
A. Froitzheim ◽  
S. Merbold ◽  
C. Egbers

Fully turbulent Taylor–Couette flow between independently rotating cylinders is investigated experimentally in a wide-gap configuration ($\unicode[STIX]{x1D702}=0.5$) around the maximum transport of angular momentum. In that regime turbulent Taylor vortices are present inside the gap, leading to a pronounced axial dependence of the flow. To account for this dependence, we measure the radial and azimuthal velocity components in horizontal planes at different cylinder heights using particle image velocimetry. The ratio of angular velocities of the cylinder walls $\unicode[STIX]{x1D707}$, where the torque maximum appears, is located in the low counter-rotating regime ($\unicode[STIX]{x1D707}_{max}(\unicode[STIX]{x1D702}=0.5)=-0.2$). This point coincides with the smallest radial gradient of angular velocity in the bulk and the detachment of the neutral surface from the outer cylinder wall, where the azimuthal velocity component vanishes. The structure of the flow is further revealed by decomposing the flow field into its large-scale and turbulent contributions. Applying this decomposition to the kinetic energy, we can analyse the formation process of the turbulent Taylor vortices in more detail. Starting at pure inner cylinder rotation, the vortices are formed and strengthened until $\unicode[STIX]{x1D707}=-0.2$ quite continuously, while they break down rapidly for higher counter-rotation. The same picture is shown by the decomposed Nusselt number, and the range of rotation ratios, where turbulent Taylor vortices can exist, shrinks strongly in comparison to investigations at much lower shear Reynolds numbers. Moreover, we analyse the scaling of the Nusselt number and the wind Reynolds number with the shear Reynolds number, finding a communal transition at approximately $Re_{S}\approx 10^{5}$ from classical to ultimate turbulence with a transitional regime lasting at least up to $Re_{S}\geqslant 2\times 10^{5}$. Including the axial dispersion of the flow into the calculation of the wind amplitude, we can also investigate the wind Reynolds number as a function of the rotation ratio $\unicode[STIX]{x1D707}$, finding a maximum in the low counter-rotating regime slightly larger than $\unicode[STIX]{x1D707}_{max}$. Based on our study it becomes clear that the investigation of counter-rotating Taylor–Couette flows strongly requires an axial exploration of the flow.


Author(s):  
S. Gilchrist ◽  
C. Y. Ching ◽  
D. Ewing

An experimental investigation was performed to determine the effect that surface roughness has on the heat transfer in an axial Taylor-Couette flow. The experiments were performed using an inner rotating cylinder in a stationary water jacket for Taylor numbers of 106 to 5×107 and axial Reynolds numbers of 900 to 2100. Experiments were performed for a smooth inner cylinder, a cylinder with two-dimensional rib roughness and a cylinder with three-dimensional cubic protrusions. The heat transfer results for the smooth cylinder were in good agreement with existing experimental data. The change in the Nusselt number was relatively independent of the axial Reynolds number for the cylinder with rib roughness. This result was similar to the smooth wall case but the heat transfer was enhanced by 5% to 40% over the Taylor number range. The Nusselt number for the cylinder with cubic protrusions exhibited an axial Reynolds number dependence. For a low axial Reynolds number of 980, the Nusselt number increased with the Taylor number in a similar way to the other test cylinders. At higher axial Reynolds numbers, the heat transfer was initially independent of the Taylor number before increasing with Taylor number similar to the lower Reynolds number case. In this higher axial Reynolds number case the heat transfer was enhanced by up to 100% at the lowest Taylor number of 1×106 and by approximately 35% at the highest Taylor number of 5×107.


2014 ◽  
Vol 748 ◽  
pp. 756-767 ◽  
Author(s):  
B. Martínez-Arias ◽  
J. Peixinho ◽  
O. Crumeyrolle ◽  
I. Mutabazi

AbstractTorque measurements in Taylor–Couette flow, with large radius ratio and large aspect ratio, over a range of velocities up to a Reynolds number of 24 000 are presented. Following a specific procedure, nine states with distinct numbers of vortices along the axis were found and the aspect ratios of the vortices were measured. The relationship between the speed and the torque for a given number of vortices is reported. In the turbulent Taylor vortex flow regime, at relatively high Reynolds number, a change in behaviour is observed corresponding to intersections of the torque–speed curves for different states. Before each intersection, the torque for a state with a larger number of vortices is higher. After each intersection, the torque for a state with a larger number of vortices is lower. The exponent, from the scaling laws of the torque, always depends on the aspect ratio of the vortices. When the Reynolds number is rescaled using the mean aspect ratio of the vortices, only a partial collapse of the exponent data is found.


2004 ◽  
Vol 67 (11) ◽  
pp. 2410-2415 ◽  
Author(s):  
L. J. FORNEY ◽  
J. A. PIERSON ◽  
Z. YE

A novel reactor is described with flow characteristics that approach that of ideal plug flow but with a residence time that is uncoupled from the hydrodynamics or boundary layer characteristics. The design described consists of an inner cylinder that rotates within a stationary but larger outer cylinder. At low rotation rates, a laminar, hydrodynamic configuration called Taylor-Couette flow is established, which consists of a system of circumferential vortices within the annular fluid gap. The latter constitutes a spatially periodic flow that is the hydrodynamic equivalent to cross flow over a tube bank or lamp array. These vortices provide radial mixing, reduce the boundary layer thickness, and are independent of the axial flow rate and thus the fluid residence time. An additional feature of the rotating design is the repetitive exposure of the fluid parcels to a minimum number of lamps, which substantially reduces the maintenance requirements. Inactivation data for Escherichia coli (ATCC 15597) were recorded in commercial apple and grape juice that are relatively opaque to UV radiation. With initial E. coli concentrations of approximately 106 CFU/ml, Taylor-Couette flow was found to provide a 3- to 5-log improvement in the inactivation efficiency compared with simple channel flow between concentric cylinders.


Sign in / Sign up

Export Citation Format

Share Document