A bathtub vortex under the influence of a protruding cylinder in a rotating tank

2013 ◽  
Vol 733 ◽  
pp. 134-157 ◽  
Author(s):  
Yin-Chung Chen ◽  
Shih-Lin Huang ◽  
Zi-Ya Li ◽  
Chien-Cheng Chang ◽  
Chin-Chou Chu

AbstractNumerical simulations and laboratory experiments were jointly conducted to investigate a bathtub vortex under the influence of a protruding cylinder in a rotating tank. In the set-up, a central drain hole is placed at the bottom of the tank and a top-down cylinder is suspended from the rigid top lid, with fluid supplied from the sidewall for mass conservation. The cylinder is protruded to produce the Taylor column effect. The flow pattern depends on the Rossby number ($\mathit{Ro}= U/ fR$), the Ekman number ($\mathit{Ek}= \nu / f{R}^{2} )$ and the height ratio, $h/ H$, where $R$ is the radius of the cylinder, $f$ is the Coriolis parameter, $\nu $ is the kinematic viscosity of the fluid, $h$ is the vertical length of the cylinder and $H$ is the height of the tank. It is found appropriate to choose $U$ to be the average inflow velocity of fluid entering the column beneath the cylinder. Steady-state solutions obtained by numerically solving the Navier–Stokes equations in the rotating frame are shown to have a good agreement with flow visualizations and particle tracking velocimetry (PTV) measurements. It is known that at $\mathit{Ro}\sim 1{0}^{- 2} $, the central downward flow surrounded by the neighbouring Ekman pumping forms a classic one-celled bathtub vortex structure when there is no protruding cylinder ($h/ H= 0$). The influence of a suspended cylinder ($h/ H\not = 0$) leads to several findings. The bathtub vortex exhibits an interesting two-celled structure with an inner Ekman pumping (EP) and an outer up-drafting motion, termed Taylor upwelling (TU). The two regions of up-drafting motion are separated by a notable finite-thickness structure, identified as a (thin-walled) Taylor column. The thickness ${ \delta }_{T}^{\ast } $ of the Taylor column is found to be well correlated to the height ratio and the Ekman number by ${\delta }_{T} = { \delta }_{T}^{\ast } / R= {(1- h/ H)}^{- 0. 32} {\mathit{Ek}}^{0. 095} $. The Taylor column presents a barrier to the fluid flow such that the fluid from the inlet may only flow into the inner region through the narrow gaps, one above the Taylor column and one beneath it (conveniently called Ekman gaps). As a result, five types of routes along which the fluid may flow to and exit at the drain hole could be identified for the multi-celled vortex structure. Moreover, the flow rates associated with the five routes were calculated and compared to help understand the relative importance of the component flow structures. The weaker influence of the Taylor column effect on the bathtub vortex at $\mathit{Ro}\sim 1$ or even higher $\mathit{Ro}\sim 1{0}^{2} $ is also discussed.

Author(s):  
Jens A. Melheim ◽  
Stefan Horender ◽  
Martin Sommerfeld

Numerical calculations of a particle-laden turbulent horizontal mixing-layer based on the Eulerian-Lagrangian approach are presented. Emphasis is given to the determination of the stochastic fluctuating fluid velocity seen by the particles in anisotropic turbulence. The stochastic process for the fluctuating velocity is a “Particle Langevin equation Model”, based on the Simplified Langevin Model. The Reynolds averaged Navier-Stokes equations are closed by the standard k-epsilon turbulence model. The calculated concentration profile and the mean, the root-mean-square (rms) and the cross-correlation terms of the particle velocities are compared with particle image velocimetry (PIV) measurements. The numerical results agree reasonably well with the PIV data for all of the mentioned quantities. The importance of the modeled vortex structure “seen” by the particles is discussed.


Author(s):  
Tong-Miin Liou ◽  
Shih-Hui Chen

Computations and measurements of time mean velocities, total fluctuation intensities, and Reynolds stresses are presented for spatially periodic flows past an array of bluff bodies aligned along the channel axis. The Reynolds number based on the channel hydraulic diameter and cross-sectional bulk mean velocity, the pitch to rib-height ratio, and the rib-height to channel-height ratio were 2 × 104, 10, and 0.133, respectively. The unsteady phase-averaged Navier-Stokes equations were solved using a Reynolds stress model with wall function and wall-related pressure strain treatment to reveal the feature of examined unsteady vortex shedding flow. Laser Doppler velocimetry measurements were performed to measure the velocity filed. Code verifications were performed through comparisons with others’ measured developing single-rib flow and our measured fully developed rib-array flow. The computed results and measured data are found in reasonable agreement, which justifies the turbulence model adopted. The calculated phase-averaged flow field clearly displays the vortex shedding behind the rib and is characterized in terms of shedding Strouhal number, vortex trajectory, vortex celerity, and vortex travelling distance in a phase cycle. Furthermore, the difference between the computed developing single-rib flow and fully developed rib-array flow is addressed.


1979 ◽  
Vol 94 (4) ◽  
pp. 711-727 ◽  
Author(s):  
M. R. Foster

Even very small Coriolis forces are shown to alter significantly the nature of the upstream wake of an object in slow (small Froude number) translation through a non-diffusive stratified fluid. If the Ekman number is of order one, the far upstream extent of the wake is reduced. If the fluid rotation is sufficient to make the Ekman number small, the contraction of the wake is much greater. We study a particular case in detail; the Ekman number is small enough to make horizontal boundary layers Ekman layers. In this case, the wake is confined to the vicinity of the object, the upstream flow arising from a combination of Ekman pumping and baroclinic vorticity generation. The upstream flow is described by an eigenfunction whose amplitude is dependent on object geometry. If the object is a semi-infinite rectangular parallelepiped, that amplitude is determined by detailed examination of the shear layer at the face of the parallelepiped and its interaction with the Ekman layer on the top surface of the object


2001 ◽  
Vol 8 (4/5) ◽  
pp. 223-239 ◽  
Author(s):  
V. I. Vlasenko ◽  
K. Hutter

Abstract. Results of an experimental and theoretical study of the interaction of a first mode internal solitary wave with a localised bottom topography (sill) are presented. Laboratory experiments have been performed in a 10m long and 0.33m wide channel filled with a stratified fluid. The interface between the two layers (fresh and salt water) is diffuse and has a finite thickness. Soliton-type disturbances of the interface having characteristics of the first baroclinic mode are generated at one channel end. They move along the channel and encounter an underwater obstacle (sill) in the middle of the channel, where they break into reflected and transmitted waves. Two types of internal waves are produced by the interaction: a fast first mode internal soliton and a slower (by a factor of approximately 3) second mode soliton-like wave. A numerical model, based on the two-dimensional Navier-Stokes equations in the Boussinesq approximation, is used tore produce the laboratory experiment. The detailed analysis of the horizontal and vertical structures of transmitted and reflected waves showed that the fast reflected and transmitted waves observed in the experiment can be interpreted as a first mode internal solitary wave whose characteristics are very close to those of the K-dV solitons. It is also demonstrated that the slow speed waves, generated during the interaction behind the first fast wave have vertical and horizontal structures very close to the second mode internal K-dV solitons.


2020 ◽  
Vol 494 (3) ◽  
pp. 3141-3155 ◽  
Author(s):  
Umin Lee

ABSTRACT We calculate small amplitude gravitational and thermal tides of uniformly rotating hot Jupiters composed of a nearly isentropic convective core and a geometrically thin radiative envelope. We treat the fluid in the convective core as a viscous fluid and solve linearized Navier–Stokes equations to obtain tidal responses of the core, assuming that the Ekman number, Ek, is a constant parameter. In the radiative envelope, we take account of the effects of radiative dissipations on the responses. The properties of tidal responses depend on thermal time-scales τ* in the envelope and Ekman number, Ek, in the core and on whether the forcing frequency ω is in the inertial range or not, where the inertial range is defined by |ω| ≤ 2Ω for the rotation frequency Ω. If Ek ≳ 10−7, the viscous dissipation in the core is dominating the thermal contributions in the envelope for τ* ≳ 1 d. If Ek ≲ 10−7, however, the viscous dissipation is comparable to or smaller than the thermal contributions and the envelope plays an important role to determine the tidal torques. If the forcing is in the inertial range, frequency resonance of the tidal forcing with core inertial modes significantly affects the tidal torques, producing numerous resonance peaks of the torque. Depending on the sign of the torque in the peaks, we suggest that there exist cases in which the resonance with core inertial modes hampers the process of synchronization between the spin and orbital motion of the planets.


1994 ◽  
Vol 116 (4) ◽  
pp. 808-814 ◽  
Author(s):  
Kwan Yeop Kim ◽  
Jae Min Hyun

A numerical and analytical study is made of spin-up from rest of a two-layer liquid in a rapidly rotating cylinder. The overall system Ekman number is small. The density of the top layer is smaller than that of the bottom layer (ρ1/ρ2<1.0), but the ratio of the individual layer kinematic viscosities is arbitrary (v1/v2<1.0 or v1/v2>1.0). The highlights of the analytical model, which is based on amended formulations of the Wedemeyer-Gerber-Homicz flow configurations, are briefly recapitulated. Comprehensive numerical solutions are secured to the time-dependent Navier–Stokes equations. The numerical solutions are validated by comparing the maximum interface displacements with the available experimental data as well as the analytical model predictions. Descriptions are made of the prominent characteristics of the interface shape for the two regimes of v1/v2<1.0 and v1/v2 > 1.0. Details of the azimuthal and meridional flow structures are illustrated by exploiting the numerical solutions. The computed meridional flows are compatible with the basic assumptions embedded in the development of the analytical model. Sequential plots of the radial profiles of azimuthal velocities are presented. These show that the global spin-up process is substantially accomplished over (En−1/2Ω−1), where En denotes the value of the smaller Ekman number of the two layers. The numerical study gives credence to the reliability and accuracy of the simplified analytical model.


1970 ◽  
Vol 68 (2) ◽  
pp. 485-491 ◽  
Author(s):  
M. J. Lighthill

AbstractWhen Rossby number is small but Ekman number is very much smaller, study of the flow field far from a body moving at right angles to the axis of rotation of a large body of fluid indicates that the region of influence should not be a Taylor column parallel to the axis, but a trailing Taylor column, bent backwards on both sides of the body at a small angle (proportional to Rossby number) to the axis. The paper reviews the physical significance of, and experimental evidence for, this conclusion.


2021 ◽  
Vol 6 (4) ◽  
pp. 1061-1077
Author(s):  
Edgar Martinez-Ojeda ◽  
Francisco Javier Solorio Ordaz ◽  
Mihir Sen

Abstract. The actuator-cylinder model was implemented in OpenFOAM by virtue of source terms in the Navier–Stokes equations. Since the stand-alone actuator cylinder is not able to properly model the wake of a vertical-axis wind turbine, the steady incompressible flow solver simpleFoam provided by OpenFOAM was used to resolve the entire flow and wakes of the turbines. The source terms are only applied inside a certain region of the computational domain, namely a finite-thickness cylinder which represents the flight path of the blades. One of the major advantages of this approach is its implicitness – that is, the velocities inside the hollow cylinder region feed the stand-alone actuator-cylinder model (AC); this in turn computes the volumetric forces and passes them to the OpenFOAM solver in order to be applied inside the hollow cylinder region. The process is repeated in each iteration of the solver until convergence is achieved. The model was compared against experimental works; wake deficits and power coefficients are used in order to assess the validity of the model. Overall, there is a good agreement of the pattern of the power coefficients according to the positions of the turbines in the array. The actual accuracy of the power coefficient depends strongly on the solidity of the turbine (actuator cylinder related) and both the inlet boundary turbulence intensity and turbulence length scale (RANS simulation related).


2021 ◽  
Vol 24 (6) ◽  
pp. 27-41
Author(s):  
V. V. Vyshinsky ◽  
K. T. Zoan

Wind boundary layer flow over the mountain landscape and large structures located around runways (RWs) creates coherent vortex structures (CVSs) that can cross a glideslope and airspace in the vicinity of an airport. The aircraft, encountering a vortex structure, experiences significant changes of the aerodynamic forces and moments, what is especially hazardous due to proximity to terrain. From a mathematical point of view, the solution of this problem presents a challenge due to extremely large space – time scale of the phenomenon, the lack of relevant atmospheric models, as well as comprehensive initial – boundary conditions in numerical modeling. In this paper, a composite solution is constructed: the CVSs area generation is computed in sufficient details within the framework of the grid method. Based on the data obtained in the approximation of analytical functions, an initial vortex structure is formed, the evolution and stochastics of which are modeled within the potential approximation by means of Rankine vortices. The evaluation of the forces and moments increment from the impact of vortex structures on the aircraft was carried out by the panel method using the engineering approach. As an example, the CVSs, resulting from wind flow around the mountainous area of the Son Tra Peninsula, that is located short of RWs 35R-17L and 35L-17R of Da Nang airport, are investigated. To improve the computational grids quality and verify the method of solving the boundary value problem for the Reynolds-averaged Navier-Stokes equations, we used the criteria based on the principle of maximum pressure, requiring Q-parameter positivity property in the vortices cores and flow separation regions. A CVS related aviation event, involving a passenger aircraft MC-21, is studied. The aircraft, after takeoff from RW 35R-17L setting the course close to the direction of the vortex wind structure axis from the Son Tra Peninsula, encountered the mountainous area CVS.


Author(s):  
Yu-Tai Lee ◽  
Chunill Hah ◽  
James Loellbach

This paper summarizes a numerical investigation of the fundamental structure of the rotor tip-clearance vortex and its interaction with a passage trailing-edge vortex in a single-stage stator-rotor pump. The flow field of a highly-loaded rotor measured in a high Reynolds number pump facility (HIREP) is used for comparison. The numerical solution was obtained by solving the three-dimensional Reynolds averaged Navier-Stokes equations. The calculated results are visualized in order to understand the details of the tip-vortex structure. The study shows that the tip geometry should be accurately represented to predict the tip-vortex structure correctly.


Sign in / Sign up

Export Citation Format

Share Document