Tilting at wave beams: a new perspective on the St. Andrew’s Cross

2017 ◽  
Vol 830 ◽  
pp. 660-680 ◽  
Author(s):  
T. Kataoka ◽  
S. J. Ghaemsaidi ◽  
N. Holzenberger ◽  
T. Peacock ◽  
T. R. Akylas
Keyword(s):  
Wave Field ◽  
Finite Length ◽  
Line Source ◽  
Cutoff Frequency ◽  
Three Dimensional ◽  
Resonance Phenomenon ◽  
Buoyancy Frequency ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Viscous Effects

The generation of internal gravity waves by a vertically oscillating cylinder that is tilted to the horizontal in a stratified Boussinesq fluid of constant buoyancy frequency, $N$, is investigated. This variant of the widely studied horizontal configuration – where a cylinder aligned with a plane of constant gravitational potential induces four wave beams that emanate from the cylinder, forming a cross pattern known as the ‘St. Andrew’s Cross’ – brings out certain unique features of radiated internal waves from a line source tilted to the horizontal. Specifically, simple kinematic considerations reveal that for a cylinder inclined by a given angle $\unicode[STIX]{x1D719}$ to the horizontal, there is a cutoff frequency, $N\sin \unicode[STIX]{x1D719}$, below which there is no longer a radiated wave field. Furthermore, three-dimensional effects due to the finite length of the cylinder, which are minor in the horizontal configuration, become a significant factor and eventually dominate the wave field as the cutoff frequency is approached; these results are confirmed by supporting laboratory experiments. The kinematic analysis, moreover, suggests a resonance phenomenon near the cutoff frequency as the group-velocity component perpendicular to the cylinder direction vanishes at cutoff; as a result, energy cannot be easily radiated away from the source, and nonlinear and viscous effects are likely to come into play. This scenario is examined by adapting the model for three-dimensional wave beams developed in Kataoka & Akylas (J. Fluid Mech., vol. 769, 2015, pp. 621–634) to the near-resonant wave field due to a tilted line source of large but finite length. According to this model, the combination of three-dimensional, nonlinear and viscous effects near cutoff triggers transfer of energy, through the action of Reynolds stresses, to a circulating horizontal mean flow. Experimental evidence of such an induced mean flow near cutoff is also presented.

scholarly journals Three-Dimensional Turbulent Vortex Shedding From a Surface-Mounted Square Cylinder: Predictions With Large-Eddy Simulations and URANS

2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi
Keyword(s):  
Experimental Data ◽  
Turbulence Model ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
Large Eddy

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.

scholarly journals Wave field and zonal flow of a librating disk

2015 ◽  
Vol 782 ◽  
pp. 178-208 ◽  
Author(s):  
Stéphane Le Dizès
Keyword(s):  
Wave Field ◽  
Rotation Rate ◽  
Zonal Flow ◽  
Shear Layers ◽  
Buoyancy Frequency ◽  
Frequency Interval ◽  
Mean Flow ◽  
Disk Radius ◽  
Harmonic Response ◽  
The Mean

In this work, we provide a viscous solution of the wave field generated by librating a disk (harmonic oscillation of the rotation rate) in a stably stratified rotating fluid. The zonal flow (mean flow correction) generated by the nonlinear interaction of the wave field is also calculated in the weakly nonlinear framework. We focus on the low dissipative limit relevant for geophysical applications and for which the wave field and the zonal flow exhibit generic features (Ekman scaling, universal structures, etc.). General expressions are obtained which depend on the disk radius $a^{\ast }$, the libration frequency ${\it\omega}^{\ast }$, the rotation rate ${\it\Omega}^{\ast }$ of the frame, the buoyancy frequency $N^{\ast }$ of the fluid, its kinematic diffusion ${\it\nu}^{\ast }$ and its thermal diffusivity ${\it\kappa}^{\ast }$. When the libration frequency is in the inertia-gravity frequency interval ($\min ({\it\Omega}^{\ast },N^{\ast })<{\it\omega}^{\ast }<\max ({\it\Omega}^{\ast },N^{\ast })$), the presence of conical internal shear layers is observed in which the spatial structures of the harmonic response and of the mean flow correction are provided. At the point of focus of these internal shear layers on the rotation axis, the largest amplitudes are obtained: the angular velocity of the harmonic response and the mean flow correction are found to be $O({\it\varepsilon}E^{-1/3})$ and $({\it\varepsilon}^{2}E^{-2/3})$ respectively, where ${\it\varepsilon}$ is the libration amplitude and $E={\it\nu}^{\ast }/({\it\Omega}^{\ast }a^{\ast 2})$ is the Ekman number. We show that the solution in the internal shear layers and in the focus region is at leading order the same as that generated by an oscillating source of axial flow localized at the edge of the disk (oscillating Dirac ring source).

Stability of internal gravity wave beams to three-dimensional modulations

2013 ◽  
Vol 736 ◽  
pp. 67-90 ◽  
Author(s):  
T. Kataoka ◽  
T. R. Akylas
Keyword(s):  
Modulational Instability ◽  
Three Dimensional ◽  
Harmonic Component ◽  
Internal Gravity Wave ◽  
Beam Width ◽  
Threshold Value ◽  
Resonant Interaction ◽  
Buoyancy Frequency ◽  
Mean Flow ◽  
Time Harmonic

AbstractThe linear stability of uniform, plane internal wave beams with locally confined spatial profile, in a stratified fluid of constant buoyancy frequency, is discussed. The associated eigenvalue problem is solved asymptotically, assuming perturbations of long wavelength relative to the beam width. In this limit, instability is found only for oblique disturbances which vary in the along-beam and the horizontal transverse directions. The mechanism of instability is a first-harmonic–mean resonant interaction between the underlying wave beam and three-dimensional perturbations that comprise a time-harmonic component, with the beam frequency, and a mean flow. Progressive beams which transport energy in one direction, in particular, are unstable if the beam steepness exceeds a certain threshold value, whereas purely standing beams are unstable even at infinitesimal steepness. A distinguishing feature of this three-dimensional modulational instability is the generation of circulating horizontal mean flows at large distances from the vicinity of the beam.

A Transport Model for the Deterministic Stresses Associated With Turbomachinery Blade Row Interactions

2000 ◽  
Vol 122 (4) ◽  
pp. 593-603 ◽  
Author(s):  
Allan G. van de Wall ◽  
Jaikrishnan R. Kadambi ◽  
John J. Adamczyk
Keyword(s):  
Turbine Blade ◽  
Transport Model ◽  
Three Dimensional ◽  
Tip Clearance ◽  
Two Dimensional ◽  
Mean Flow ◽  
Viscous Effects ◽  
Vortical Structures ◽  
Blade Row ◽  
The Mean

The unsteady process resulting from the interaction of upstream vortical structures with a downstream blade row in turbomachines can have a significant impact on the machine efficiency. The upstream vortical structures or disturbances are transported by the mean flow of the downstream blade row, redistributing the time-average unsteady kinetic energy (K) associated with the incoming disturbance. A transport model was developed to take this process into account in the computation of time-averaged multistage turbomachinery flows. The model was applied to compressor and turbine geometry. For compressors, the K associated with upstream two-dimensional wakes and three-dimensional tip clearance flows is reduced as a result of their interaction with a downstream blade row. This reduction results from inviscid effects as well as viscous effects and reduces the loss associated with the upstream disturbance. Any disturbance passing through a compressor blade row results in a smaller loss than if the disturbance was mixed-out prior to entering the blade row. For turbines, the K associated with upstream two-dimensional wakes and three-dimensional tip clearance flows are significantly amplified by inviscid effects as a result of the interaction with a downstream turbine blade row. Viscous effects act to reduce the amplification of the K by inviscid effects but result in a substantial loss. Two-dimensional wakes and three-dimensional tip clearance flows passing through a turbine blade row result in a larger loss than if these disturbances were mixed-out prior to entering the blade row. [S0889-504X(00)01804-3]

Measurement of the Reynolds stresses and the mean-flow field in a three-dimensional pressure-driven boundary layer

1982 ◽  
Vol 119 ◽  
pp. 121-153 ◽  
Author(s):  
Udo R. Müller
Keyword(s):  
Boundary Layer ◽  
Shear Stress ◽  
Flow Field ◽  
Three Dimensional ◽  
Turbulent Shear ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Mean Velocity ◽  
Flow Acceleration ◽  
The Mean

An experimental study of a steady, incompressible, three-dimensional turbulent boundary layer approaching separation is reported. The flow field external to the boundary layer was deflected laterally by turning vanes so that streamwise flow deceleration occurred simultaneous with cross-flow acceleration. At 21 stations profiles of the mean-velocity components and of the six Reynolds stresses were measured with single- and X-hot-wire probes, which were rotatable around their longitudinal axes. The calibration of the hot wires with respect to magnitude and direction of the velocity vector as well as the method of evaluating the Reynolds stresses from the measured data are described in a separate paper (Müller 1982, hereinafter referred to as II). At each measuring station the wall shear stress was inferred from a Preston-tube measurement as well as from a Clauser chart. With the measured profiles of the mean velocities and of the Reynolds stresses several assumptions used for turbulence modelling were checked for their validity in this flow. For example, eddy viscosities for both tangential directions and the corresponding mixing lengths as well as the ratio of resultant turbulent shear stress to turbulent kinetic energy were derived from the data.

The nonlinear stability of a free shear layer in the viscous critical layer régime

1980 ◽  
Vol 293 (1408) ◽  
pp. 643-672 ◽  
Keyword(s):  
Shear Layer ◽  
Critical Layer ◽  
Reynolds Stresses ◽  
Free Shear Layer ◽  
Mean Flow ◽  
Viscous Effects ◽  
Flow Distortion ◽  
Threshold Amplitude ◽  
The Mean ◽  
Valid Solution

The nonlinear evolution of weakly amplified waves in a hyperbolic tangent free shear layer is described for spatially and temporally growing waves when the shear layer Reynolds number is large and the critical layer viscous. An artificial body force is introduced in order to keep the mean flow parallel. Jump conditions on the perturbation velocity and mean vorticity are derived across the critical layer by applying the method of matched asymptotic expansions and it is shown that viscous effects outside the critical layer have to be taken into account in order to obtain a uniformly valid solution. Consequently the true neutral wavenumber and frequency are lower than their inviscid counterparts. When only the harmonic fluctuations are considered, it is known that the Landau constant is negative so that linearly amplified disturbances reach an equilibrium amplitude. It is shown that when the mean flow distortion generated by Reynolds stresses is also included, the Landau constant becomes positive. Thus, in both the spatial and temporal case, linearly amplified waves are further destabilized and damped waves are unstable above a threshold amplitude.

Direct Numerical Simulations of Planar and Cylindrical Density Currents

2005 ◽  
Vol 73 (6) ◽  
pp. 923-930 ◽  
Author(s):  
Mariano I. Cantero ◽  
S. Balachandar ◽  
Marcelo H. García ◽  
James P. Ferry
Keyword(s):  
Stokes Equations ◽  
Three Dimensional ◽  
Boussinesq Approximation ◽  
The Body ◽  
Density Current ◽  
Density Currents ◽  
Mean Flow ◽  
Viscous Effects ◽  
Spreading Velocity ◽  
Slip Boundary

The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof (Gr) numbers are investigated (105 and 1.5×106) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial three-dimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.

Multiple Jets in a Crossflow: Detailed Measurements and Numerical Simulations

1997 ◽  
Vol 119 (2) ◽  
pp. 330-342 ◽  
Author(s):  
P. Ajersch ◽  
J.-M. Zhou ◽  
S. Ketler ◽  
M. Salcudean ◽  
I. S. Gartshore
Keyword(s):  
Numerical Simulations ◽  
Film Cooling ◽  
Three Dimensional ◽  
Light Sheet ◽  
Algebraic Model ◽  
Laser Doppler Velocimeter ◽  
Width Ratio ◽  
Video Tape ◽  
Reynolds Stresses ◽  
Mean Flow

The fluid mechanics and heat transfer characteristics of film cooling are three-dimensional and highly complex. To understand this problem better, an experimental study was conducted in a low-speed wind tunnel on a row of six rectangular jets injected at 90 deg to the crossflow (mainstream flow). The jet-to-crossflow velocity ratios (blowing ratios) examined were 0.5, 1.0, and 1.5, and the jet spacing-to-jet width ratio was 3.0. No significant temperature difference between jet and crossflow air was introduced. Mean velocities and six flow stresses were measured using a three-component laser-Doppler velocimeter operating in coincidence mode. Seeding of both jet and cross-stream air was achieved with a commercially available smoke generator. Flow statistics are reported in the form of vector plots, contours, and x-y graphs, showing velocity, turbulence intensity, and Reynolds stresses. To complement the detailed measurements, flow visualization was accomplished by transmitting the laser beam through a cylindrical lens, thereby generating a narrow, intense sheet of light. Jet air only was seeded with smoke, which was illuminated in the plane of the light sheet. Therefore, it was possible to record on video tape the trajectory and penetration of the jets in the crossflow. Selected still images from the recordings are presented. Numerical simulations of the observed flow field were made by using a multigrid, segmented, k–ε CFD code. Special near-wall treatment included a nonisotropic formulation for the effective viscosity, a low-Re model for k, and an algebraic model for the length scale. Comparisons between the measured and computed velocities show good agreement for the nonuniform mean flow at the jet exit plane. Velocities and stresses on the jet centerline downstream of the orifice are less well predicted, probably because of inadequate turbulence modeling, while values off the centerline match those of the experiments much more closely.

Strato-rotational instability without resonance

2018 ◽  
Vol 846 ◽  
pp. 815-833 ◽  
Author(s):  
Chen Wang ◽  
Neil J. Balmforth
Keyword(s):  
Three Dimensional ◽  
Normal Modes ◽  
Phase Speed ◽  
Flow Speed ◽  
Rotational Instability ◽  
Buoyancy Frequency ◽  
Critical Levels ◽  
Mean Flow ◽  
Linear Eigenvalue Problem ◽  
Resonant Interactions

Strato-rotational instability (SRI) is normally interpreted as the resonant interactions between normal modes of the internal or Kelvin variety in three-dimensional settings in which the stratification and rotation are orthogonal to both the background flow and its shear. Using a combination of asymptotic analysis and numerical solution of the linear eigenvalue problem for plane Couette flow, it is shown that such resonant interactions can be destroyed by certain singular critical levels. These levels are not classical critical levels, where the phase speed $c$ of a normal mode matches the mean flow speed $U$, but are a different type of singularity where $(c-U)$ matches a characteristic gravity-wave speed $\pm N/k$, based on the buoyancy frequency $N$ and streamwise horizontal wavenumber $k$. Instead, it is shown that a variant of SRI can occur due to the coupling of a Kelvin or internal wave to such ‘baroclinic’ critical levels. Two characteristic situations are identified and explored, and the conservation law for pseudo-momentum is used to rationalize the physical mechanism of instability. The critical level coupling removes the requirement for resonance near specific wavenumbers $k$, resulting in an extensive continuous band of unstable modes.

Calculation of turbulence-driven secondary motion in non-circular ducts

1984 ◽  
Vol 140 ◽  
pp. 189-222 ◽  
Author(s):  
A. O. Demuren ◽  
W. Rodi
Keyword(s):  
Eddy Viscosity ◽  
Three Dimensional ◽  
Numerical Procedure ◽  
Momentum Equation ◽  
Shear Stresses ◽  
Reynolds Stresses ◽  
Calculation Methods ◽  
Mean Flow ◽  
Square Duct ◽  
Secondary Motion

Experiments on and calculation methods for flow in straight non-circular ducts involving turbulence-driven secondary motion are reviewed. The origin of the secondary motion and the shortcomings of existing calculation methods are discussed. A more refined model is introduced, in which algebraic expressions are derived for the Reynolds stresses in the momentum equations for the secondary motion by simplifying the modelled Reynolds-stress equations of Launder, Reece & Rodi (1975), while a simple eddy-viscosity model is used for the shear stresses in the axial momentum equation. The kinetic energy k and the dissipation rate ε of the turbulent motion which appear in the algebraic and the eddy-viscosity expressions are determined from transport equations. The resulting set of equations is solved with a forward-marching numerical procedure for three-dimensional shear layers. The model, as well as a version proposed by Naot & Rodi (1982), is tested by application to developing flow in a square duct and to developed flow in a partially roughened rectangular duct investigated experimentally by Hinze (1973). In both cases, the main features of the mean-flow and the turbulence quantities are simulated realistically by both models, but the present model underpredicts the secondary velocity while the Naot-Rodi model tends to overpredict it.

Sign in / Sign up

Export Citation Format

Share Document