A numerical study of global frequency selection in the time-mean wake of a circular cylinder

2010 ◽  
Vol 645 ◽  
pp. 435-446 ◽  
Author(s):  
J. S. LEONTINI ◽  
M. C. THOMPSON ◽  
K. HOURIGAN

A series of direct numerical simulations, both in two- and three-dimensions, of the flow past a circular cylinder for Reynolds numbers Re ≤ 600 has been conducted. From these simulations, the time-mean (and, for the three-dimensional simulations, the spanwise spatial-mean) flow has been calculated. A global linear stability analysis has been conducted on these mean flows, showing that the mean cylinder wake for Re ≤ 600 is marginally stable and the eigenfrequency of the leading global mode closely predicts the eventual saturated vortex shedding frequency. A local stability analysis has also been conducted. For this, a series of streamwise velocity profiles has been extracted from the mean wake and the stability of these profiles has been analysed using the Rayleigh stability equation. The real and imaginary instability frequencies gained from these profiles have then been used to find the global frequency selected by the flow using a saddle-point criterion. The results confirm the success of the saddle-point criterion when the mean flow is quasi-parallel in the vicinity of the saddle point; however, the limitations of the method when the mean flow exhibits higher curvature are also elucidated.

1979 ◽  
Vol 92 (2) ◽  
pp. 269-301 ◽  
Author(s):  
R. E. Britter ◽  
J. C. R. Hunt ◽  
J. C. Mumford

The flow of grid-generated turbulence past a circular cylinder is investigated using hot-wire anemometry over a Reynolds number range from 4·25 × 103 to 2·74 × 104 and a range of intensities from 0·025 to 0·062. Measurements of the mean velocity distribution, and r.m.s. intensities and spectral energy densities of the turbulent velocity fluctuations are presented for various radial and circumferential positions relative to the cylinder, and for ratios of the cylinder radius a to the scale of the incident turbulence Lx ranging from 0·05 to 1·42. The influence of upstream conditions on the flow in the cylinder wake and its associated induced velocity fluctuations is discussed.For all measurements, detailed comparison is made with the theoretical predictions of Hunt (1973). We conclude the following. The amplification and reduction of the three components of turbulence (which occur in different senses for the different components) can be explained qualitatively in terms of the distortion by the mean flow of the turbulent vorticity and the ‘blocking’ or ‘source’ effect caused by turbulence impinging on the cylinder surface. The relative importance of the first effect over the second increases as a/Lx increases or the distance from the cylinder surface increases.Over certain ranges of the variables involved, the measurements are in quantitative agreement with the predictions of the asymptotic theory when a/Lx [Lt ] 1, a/Lx [Gt ] 1 or |k| a [Gt ] 1 (where k is the wavenumber).The incident turbulence affects the gross properties of the flow in the cylinder wake, but the associated velocity fluctuations are probably statistically independent of those in the incident flow.The dissipation of turbulent energy is greater in the straining flow near the cylinder than in the approach flow. Some estimates for this effect are proposed.


Author(s):  
Alexander Vakhrushev ◽  
Abdellah Kharicha ◽  
Ebrahim Karimi-Sibaki ◽  
Menghuai Wu ◽  
Andreas Ludwig ◽  
...  

AbstractA numerical study is presented that deals with the flow in the mold of a continuous slab caster under the influence of a DC magnetic field (electromagnetic brakes (EMBrs)). The arrangement and geometry investigated here is based on a series of previous experimental studies carried out at the mini-LIMMCAST facility at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR). The magnetic field models a ruler-type EMBr and is installed in the region of the ports of the submerged entry nozzle (SEN). The current article considers magnet field strengths up to 441 mT, corresponding to a Hartmann number of about 600, and takes the electrical conductivity of the solidified shell into account. The numerical model of the turbulent flow under the applied magnetic field is implemented using the open-source CFD package OpenFOAM®. Our numerical results reveal that a growing magnitude of the applied magnetic field may cause a reversal of the flow direction at the meniscus surface, which is related the formation of a “multiroll” flow pattern in the mold. This phenomenon can be explained as a classical magnetohydrodynamics (MHD) effect: (1) the closure of the induced electric current results not primarily in a braking Lorentz force inside the jet but in an acceleration in regions of previously weak velocities, which initiates the formation of an opposite vortex (OV) close to the mean jet; (2) this vortex develops in size at the expense of the main vortex until it reaches the meniscus surface, where it becomes clearly visible. We also show that an acceleration of the meniscus flow must be expected when the applied magnetic field is smaller than a critical value. This acceleration is due to the transfer of kinetic energy from smaller turbulent structures into the mean flow. A further increase in the EMBr intensity leads to the expected damping of the mean flow and, consequently, to a reduction in the size of the upper roll. These investigations show that the Lorentz force cannot be reduced to a simple damping effect; depending on the field strength, its action is found to be topologically complex.


2006 ◽  
Vol 24 (12) ◽  
pp. 3241-3256 ◽  
Author(s):  
C. M. Huang ◽  
S. D. Zhang ◽  
F. Yi

Abstract. To quantitatively study the effects of nonlinear interactions on tide structure, a nonlinear numerical tidal model is developed, and the reliability and convergence of the adopted algorithm and coding are checked by numerical experiments. Under the same conditions as those employed by the GSWM-00 (Global Scale Wave Model 2000), our model provides the nonlinear quasi-steady solution of the migrating semidiurnal tide, which differs from the GSWM-00 result (the linear steady solution) in the MLT region, especially above 100 km. Additionally, their amplitude difference displays a remarkable month-to-month variation, and its significant magnitudes occur during the month with strong semidiurnal tide. A quantitative analysis suggests that the main cause for the amplitude difference is that the initial migrating 12-h tide will interact with the mean flow as well as the nonlinearity-excited 6-h tide, and subsequently yield a new 12-h tidal part. Furthermore, our simulations also show that the mean flow/tidal interaction will significantly alter the background wind and temperature fields. The large magnitudes of the tidal amplitude difference and the background alteration indicate that the nonlinear processes involved in tidal propagations should be comprehensively considered in the description of global atmospheric dynamics in the MLT region. The comparisons among our simulations, the GSWMs and some observations of tides suggest that the nonlinearity-induced tidal structure variation could be a possible mechanism to account for some discrepancies between the GSWMs and the observations.


2019 ◽  
Vol 104 (2-3) ◽  
pp. 533-552
Author(s):  
J. Casacuberta ◽  
K. J. Groot ◽  
Q. Ye ◽  
S. Hickel

AbstractMicro-ramps are popular passive flow control devices which can delay flow separation by re-energising the lower portion of the boundary layer. We compute the laminar base flow, the instantaneous transitional flow, and the mean flow around a micro-ramp immersed in a quasi-incompressible boundary layer at supercritical roughness Reynolds number. Results of our Direct Numerical Simulations (DNS) are compared with results of BiLocal stability analysis on the DNS base flow and independent tomographic Particle Image Velocimetry (tomo-PIV) experiments. We analyse relevant flow structures developing in the micro-ramp wake and assess their role in the micro-ramp functionality, i.e., in increasing the near-wall momentum. The main flow feature of the base flow is a pair of streamwise counter-rotating vortices induced by the micro-ramp, the so-called primary vortex pair. In the instantaneous transitional flow, the primary vortex pair breaks up into large-scale hairpin vortices, which arise due to linear varicose instability of the base flow, and unsteady secondary vortices develop. Instantaneous vortical structures obtained by DNS and experiments are in good agreement. Matching linear disturbance growth rates from DNS and linear stability analysis are obtained until eight micro-ramp heights downstream of the micro-ramp. For the setup considered in this article, we show that the working principle of the micro-ramp is different from that of classical vortex generators; we find that transitional perturbations are more efficient in increasing the near-wall momentum in the mean flow than the laminar primary vortices in the base flow.


2015 ◽  
Vol 783 ◽  
pp. 379-411 ◽  
Author(s):  
I. Marusic ◽  
K. A. Chauhan ◽  
V. Kulandaivelu ◽  
N. Hutchins

In this paper we study the spatial evolution of zero-pressure-gradient (ZPG) turbulent boundary layers from their origin to a canonical high-Reynolds-number state. A prime motivation is to better understand under what conditions reliable scaling behaviour comparisons can be made between different experimental studies at matched local Reynolds numbers. This is achieved here through detailed streamwise velocity measurements using hot wires in the large University of Melbourne wind tunnel. By keeping the unit Reynolds number constant, the flow conditioning, contraction and trip can be considered unaltered for a given boundary layer’s development and hence its evolution can be studied in isolation from the influence of inflow conditions by moving to different streamwise locations. Careful attention was given to the experimental design in order to make comparisons between flows with three different trips while keeping all other parameters nominally constant, including keeping the measurement sensor size nominally fixed in viscous wall units. The three trips consist of a standard trip and two deliberately ‘over-tripped’ cases, where the initial boundary layers are over-stimulated with additional large-scale energy. Comparisons of the mean flow, normal Reynolds stress, spectra and higher-order turbulence statistics reveal that the effects of the trip are seen to be significant, with the remnants of the ‘over-tripped’ conditions persisting at least until streamwise stations corresponding to $Re_{x}=1.7\times 10^{7}$ and $x=O(2000)$ trip heights are reached (which is specific to the trips used here), at which position the non-canonical boundary layers exhibit a weak memory of their initial conditions at the largest scales $O(10{\it\delta})$, where ${\it\delta}$ is the boundary layer thickness. At closer streamwise stations, no one-to-one correspondence is observed between the local Reynolds numbers ($Re_{{\it\tau}}$, $Re_{{\it\theta}}$ or $Re_{x}$ etc.), and these differences are likely to be the cause of disparities between previous studies where a given Reynolds number is matched but without account of the trip conditions and the actual evolution of the boundary layer. In previous literature such variations have commonly been referred to as low-Reynolds-number effects, while here we show that it is more likely that these differences are due to an evolution effect resulting from the initial conditions set up by the trip and/or the initial inflow conditions. Generally, the mean velocity profiles were found to approach a constant wake parameter ${\it\Pi}$ as the three boundary layers developed along the test section, and agreement of the mean flow parameters was found to coincide with the location where other statistics also converged, including higher-order moments up to tenth order. This result therefore implies that it may be sufficient to document the mean flow parameters alone in order to ascertain whether the ZPG flow, as described by the streamwise velocity statistics, has reached a canonical state, and a computational approach is outlined to do this. The computational scheme is shown to agree well with available experimental data.


2016 ◽  
Vol 138 (11) ◽  
Author(s):  
Anna E. Craig ◽  
John O. Dabiri ◽  
Jeffrey R. Koseff

Experimental data are presented for large arrays of rotating, variable-height cylinders in order to study the dependence of the three-dimensional mean flows on the height heterogeneity of the array. Elements in the examined arrays were spatially arranged in the same staggered paired configuration, and the heights of each element pair varied up to ±37.5% from the mean height (kept constant across all arrays), such that the arrays were vertically structured. Four vertical structuring configurations were examined at a nominal Reynolds number (based on freestream velocity and cylinder diameter) of 600 and nominal tip-speed ratios of 0, 2, and 4. It was found that the vertical structuring of the array could significantly alter the mean flow patterns. Most notably, a net vertical flow into the array from above was observed, which was augmented by the arrays' vertical structuring, showing a 75% increase from the lowest to highest vertical flows (as evaluated at the maximum element height, at a single rotation rate). This vertical flow into the arrays is of particular interest as it represents an additional mechanism by which high streamwise momentum can be transported from above the array down into the array. An evaluation of the streamwise momentum resource within the array indicates up to a 56% increase in the incoming streamwise velocity to the elements (from the lowest to highest ranking arrays, at a single rotation rate). These arrays of rotating cylinders may provide insight into the flow kinematics of arrays of vertical axis wind turbines (VAWTs). In a physical VAWT array, an increase in incoming streamwise flow velocity to a turbine corresponds to a (cubic) increase in the power output of the turbine. Thus, these results suggest a promising approach to increasing the power output of a VAWT array.


2011 ◽  
Vol 679 ◽  
pp. 383-414 ◽  
Author(s):  
K. OBERLEITHNER ◽  
M. SIEBER ◽  
C. N. NAYERI ◽  
C. O. PASCHEREIT ◽  
C. PETZ ◽  
...  

The spatio-temporal evolution of a turbulent swirling jet undergoing vortex breakdown has been investigated. Experiments suggest the existence of a self-excited global mode having a single dominant frequency. This oscillatory mode is shown to be absolutely unstable and leads to a rotating counter-winding helical structure that is located at the periphery of the recirculation zone. The resulting time-periodic 3D velocity field is predicted theoretically as being the most unstable mode determined by parabolized stability analysis employing the mean flow data from experiments. The 3D oscillatory flow is constructed from uncorrelated 2D snapshots of particle image velocimetry data, using proper orthogonal decomposition, a phase-averaging technique and an azimuthal symmetry associated with helical structures. Stability-derived modes and empirically derived modes correspond remarkably well, yielding prototypical coherent structures that dominate the investigated flow region. The proposed method of constructing 3D time-periodic velocity fields from uncorrelated 2D data is applicable to a large class of turbulent shear flows.


2012 ◽  
Vol 1 (33) ◽  
pp. 41
Author(s):  
Kian Yew Lim ◽  
Ole Secher Madsen ◽  
Hin Fatt Cheong

An experimental study involving near-orthogonal wave-current interaction in a wave basin is reported in this paper. Due to previous shortcomings associated with 2D bottom configurations, i.e. occurrence of ripple-induced turning of flows close to the bed, the present experiments were conducted with the bottom covered by closely packed ceramic marbles (mean diameter of 1.25cm). Three types of flows were generated over this bottom: current-alone, wave-alone and combined wave-current flow. For current-alone and wave-current cases, the log-profile analysis was used to resolve the equivalent Nikuradse sand grain roughness, kn, while the energy dissipation method was used to estimate kn for wave-alone case. The results show that kn obtained for current- and wave-alone tests is roughly 2.2 times the diameter of the marbles. For orthogonal wave-current flows, the kn value, when used in combination with the Grant-Madsen (GM) model to reproduce the experimental apparent roughness, is found to be smaller than the measured current-alone and wave-alone kn. Similar under-prediction of bottom roughness is also observed when the GM model is compared with a numerical study, thus supporting the conjecture that when the current is weak compared to the waves, simple theoretical models like GM are not sufficiently sensitive to the angle of wave-current interaction. Experiments with currents at angles of 60° and 120° to the wave direction yield apparent roughness smaller than the 90° case, which is counter-intuitive since one would expect the mean flow to experience a stronger wave-induced turbulence when it is more aligned with the wave direction. This result indicates a possible contamination from wave-induced mass transport to the mean flow profile for non-orthogonal combined flow cases, and therefore highlights the need for other alternatives to the log-profile analysis when attempting to resolve kn from current velocity profiles from combined wave-current flows.


1999 ◽  
Vol 122 (1) ◽  
pp. 179-183 ◽  
Author(s):  
Robert E. Spall ◽  
Blake M. Ashby

Solutions to the incompressible Reynolds-averaged Navier–Stokes equations have been obtained for turbulent vortex breakdown within a slightly diverging tube. Inlet boundary conditions were derived from available experimental data for the mean flow and turbulence kinetic energy. The performance of both two-equation and full differential Reynolds stress models was evaluated. Axisymmetric results revealed that the initiation of vortex breakdown was reasonably well predicted by the differential Reynolds stress model. However, the standard K-ε model failed to predict the occurrence of breakdown. The differential Reynolds stress model also predicted satisfactorily the mean azimuthal and axial velocity profiles downstream of the breakdown, whereas results using the K-ε model were unsatisfactory. [S0098-2202(00)01601-1]


2016 ◽  
Vol 800 ◽  
pp. 327-357 ◽  
Author(s):  
P. Meliga ◽  
E. Boujo ◽  
F. Gallaire

We use the adjoint method to compute sensitivity maps for the limit-cycle frequency and amplitude of the Bénard–von Kármán vortex street in the wake of a circular cylinder. The sensitivity analysis is performed in the frame of the semi-linear self-consistent model recently introduced by Mantič et al. (Phys. Rev. Lett., vol. 113, 2014, 084501), which allows us to describe accurately the effect of the control on the mean flow, but also on the finite-amplitude fluctuation that couples back nonlinearly onto the mean flow via the formation of Reynolds stress. The sensitivity is computed with respect to arbitrary steady and synchronous time-harmonic body forces. For a small amplitude of the control, the theoretical variations of the limit-cycle frequency predict well those of the controlled flow, as obtained from either self-consistent modelling or direct numerical simulation of the Navier–Stokes equations. This is not the case if the variations are computed in the simpler mean flow approach overlooking the coupling between the mean and fluctuating components of the flow perturbation induced by the control. The variations of the limit-cycle amplitude (that falls out the scope of the mean flow approach) are also correctly predicted, meaning that the approach can serve as a relevant and systematic guideline to control strongly unstable flows exhibiting non-small, finite amplitudes of oscillation. As an illustration, we apply the method to control by means of a small secondary control cylinder and discuss the obtained results in the light of the seminal experiments of Strykowski & Sreenivasan (J. Fluid Mech., vol. 218, 1990, pp. 71–107).


Sign in / Sign up

Export Citation Format

Share Document