scholarly journals Baroclinic instability of axially symmetric flow over sloping bathymetry

2016 ◽  
Vol 799 ◽  
pp. 265-296 ◽  
Author(s):  
Aviv Solodoch ◽  
Andrew L. Stewart ◽  
James C. McWilliams

Observations and models of deep ocean boundary currents show that they exhibit complex variability, instabilities and eddy shedding, particularly over continental slopes that curve horizontally, for example around coastal peninsulas. In this article the authors investigate the source of this variability by characterizing the properties of baroclinic instability in mean flows over horizontally curved bottom slopes. The classical two-layer quasi-geostrophic solution for linear baroclinic instability over sloping bottom topography is extended to the case of azimuthal mean flow in an annular channel. To facilitate comparison with the classical straight channel instability problem of uniform mean flow, the authors focus on comparatively simple flows in an annulus, namely uniform azimuthal velocity and solid-body rotation. Baroclinic instability in solid-body rotation flow is analytically analogous to the instability in uniform straight channel flow due to several identical properties of the mean flow, including vanishing strain rate and vorticity gradient. The instability of uniform azimuthal flow is numerically similar to straight channel flow instability as long as the mean barotropic azimuthal velocity is zero. Non-zero barotropic flow generally suppresses the instability via horizontal curvature-induced strain and Reynolds stress work. An exception occurs when the ratio of the bathymetric to isopycnal slopes is close to (positive) one, as is often observed in the ocean, in which case the instability is enhanced. A non-vanishing mean barotropic flow component also results in a larger number of growing eigenmodes and in increased non-normal growth. The implications of these findings for variability in deep western boundary currents are discussed.

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.


2018 ◽  
Vol 40 ◽  
pp. 05039
Author(s):  
Priscilla Williams ◽  
Vesselina Roussinova ◽  
Ram Balachandar

This paper focuses on the turbulence structure in a non-uniform, gradually varied, sub-critical open channel flow (OCF) on a rough bed. The flow field is analysed under accelerating, near-uniform and decelerating conditions. Information for the flow and turbulence parameters was obtained at multiple sections and planes using two different techniques: two-component laser Doppler velocimetry (LDV) and particle image velocimetry (PIV). Different outer region velocity scaling methods were explored for evaluation of the local friction velocity. Analysis of the mean velocity profiles showed that the overlap layer exists for all flow cases. The outer layer of the decelerated velocity profile was strongly affected by the pressure gradient, where a large wake was noted. Due to the prevailing nature of the experimental setup it was found that the time-averaged flow quantities do not attained equilibrium conditions and the flow is spatially heterogeneous. The roughness generally increases the friction velocity and its effect was stronger than the effect of the pressure gradient. It was found that for the decelerated flow section over a rough bed, the mean flow and turbulence intensities were affected throughout the flow depth. The flow features presented in this study can be used to develop a model for simulating flow over a block ramp. The effect of the non-uniformity and roughness on turbulence intensities and Reynolds shear stresses was further investigated.


2009 ◽  
Vol 638 ◽  
pp. 423-452 ◽  
Author(s):  
EMMANUEL MIGNOT ◽  
D. HURTHER ◽  
E. BARTHELEMY

This study examines the structure of shear stress and turbulent kinetic energy (TKE) flux across the roughness layer of a uniform, fully rough gravel-bed channel flow (ks+ ≫ 100, δ/k = 20) using high-resolution acoustic Doppler velocity profiler measurements. The studied gravel-bed roughness layer exhibits a complex random multi-scale roughness structure in strong contrast with conceptualized k- or d-type roughness in standard rough-wall flows. Within the roughness layer, strong spatial variability of all time-averaged flow quantities are observed affecting up to 40% of the boundary layer height. This variability is attributed to the presence of bed zones with emanating bed protuberances (or gravel clusters) acting as local flow obstacles and bed zones of more homogenous roughness of densely packed gravel elements. Considering the strong spatial mean flow variability across the roughness layer, a spatio-temporal averaging procedure, called double averaging (DA), has been applied to the analysed flow quantities. Three aspects have been addressed: (a) the DA shear stress and DA TKE flux in specific bed zones associated with three classes of velocity profiles as previously proposed in Mignot, Barthélemy & Hurther (J. Fluid Mech., vol. 618, 2009, p. 279), (b) the global and per class DA conditional statistics of shear stress and associated TKE flux and (c) the contribution of large-scale coherent shear stress structures (LC3S) to the TKE flux across the roughness layer. The mean Reynolds and dispersive shear structure show good agreement between the protuberance bed zones associated with the S-shape/accelerated classes and recent results obtained in standard k-type rough-wall flows (Djenidi et al., Exp. Fluids, vol. 44, 2008, p. 37; Pokrajac, McEwan & Nikora, Exp. Fluids, vol. 45, 2008, p. 73). These gravel-bed protuberances act as local flow obstacles inducing a strong turbulent activity in their wake regions. The conditional statistics show that the Reynolds stress contribution is fairly well distributed between sweep and ejection events, with threshold values ranging from H = 0 to H = 8. However, the TKE flux across the roughness layer primarily results from the residual shear stress between ejection and sweep of very high magnitude (H = 10–20) and of small turbulent scale. Although LC3S are seen to penetrated the interfacial roughness layer, their TKE flux contribution is found to be negligible compared to the very energetic small-scale sweep events. These sweeps are dominantly produced in the bed zones of local gravel protuberances where the velocity profiles are inflexional of S-shape type and the mean flow properties are of mixing-layer flow type as previously shown in Mignot et al. (2009).


2012 ◽  
Vol 712 ◽  
pp. 169-202 ◽  
Author(s):  
A. Busse ◽  
N. D. Sandham

AbstractThe effects of rough surfaces on turbulent channel flow are modelled by an extra force term in the Navier–Stokes equations. This force term contains two parameters, related to the density and the height of the roughness elements, and a shape function, which regulates the influence of the force term with respect to the distance from the channel wall. This permits a more flexible specification of a rough surface than a single parameter such as the equivalent sand grain roughness. The effects of the roughness force term on turbulent channel flow have been investigated for a large number of parameter combinations and several shape functions by direct numerical simulations. It is possible to cover the full spectrum of rough flows ranging from hydraulically smooth through transitionally rough to fully rough cases. By using different parameter combinations and shape functions, it is possible to match the effects of different types of rough surfaces. Mean flow and standard turbulence statistics have been used to compare the results to recent experimental and numerical studies and a good qualitative agreement has been found. Outer scaling is preserved for the streamwise velocity for both the mean profile as well as its mean square fluctuations in all but extremely rough cases. The structure of the turbulent flow shows a trend towards more isotropic turbulent states within the roughness layer. In extremely rough cases, spanwise structures emerge near the wall and the turbulent state resembles a mixing layer. A direct comparison with the study of Ashrafian, Andersson & Manhart (Intl J. Heat Fluid Flow, vol. 25, 2004, pp. 373–383) shows a good quantitative agreement of the mean flow and Reynolds stresses everywhere except in the immediate vicinity of the rough wall. The proposed roughness force term may be of benefit as a wall model for direct and large-eddy numerical simulations in cases where the exact details of the flow over a rough wall can be neglected.


2007 ◽  
Vol 37 (9) ◽  
pp. 2267-2289 ◽  
Author(s):  
Richard G. Williams ◽  
Chris Wilson ◽  
Chris W. Hughes

Abstract Signatures of eddy variability and vorticity forcing are diagnosed in the atmosphere and ocean from weather center reanalysis and altimetric data broadly covering the same period, 1992–2002. In the atmosphere, there are localized regions of eddy variability referred to as storm tracks. At the entrance of the storm track the eddies grow, providing a downgradient heat flux and accelerating the mean flow eastward. At the exit and downstream of the storm track, the eddies decay and instead provide a westward acceleration. In the ocean, there are similar regions of enhanced eddy variability along the extension of midlatitude boundary currents and the Antarctic Circumpolar Current. Within these regions of high eddy kinetic energy, there are more localized signals of high Eady growth rate and downgradient eddy heat fluxes. As in the atmosphere, there are localized regions in the Southern Ocean where ocean eddies provide statistically significant vorticity forcing, which acts to accelerate the mean flow eastward, provide torques to shift the jet, or decelerate the mean flow. These regions of significant eddy vorticity forcing are often associated with gaps in the topography, suggesting that the ocean jets are being locally steered by topography. The eddy forcing may also act to assist in the separation of boundary currents, although the diagnostics of this study suggest that this contribution is relatively small when compared with the advection of planetary vorticity by the time-mean flow.


2009 ◽  
Vol 66 (4) ◽  
pp. 837-863 ◽  
Author(s):  
Pablo Zurita-Gotor ◽  
Geoffrey K. Vallis

Abstract This paper investigates the equilibration of baroclinic turbulence in an idealized, primitive equation, two-level model, focusing on the relation with the phenomenology of quasigeostrophic turbulence theory. Simulations with a comparable two-layer quasigeostrophic model are presented for comparison, with the deformation radius in the quasigeostrophic model being set using the stratification from the primitive equation model. Over a fairly broad parameter range, the primitive equation and quasigeostrophic results are in qualitative and, to some degree, quantitative agreement and are consistent with the phenomenology of geostrophic turbulence. The scale, amplitude, and baroclinicity of the eddies and the degree of baroclinic instability of the mean flow all vary fairly smoothly with the imposed parameters; both models are able, in some parameter ranges, to produce supercritical flows. The criticality in the primitive equation model, which does not have any convective parameterization scheme, is fairly sensitive to the external parameters, most notably the planet size (i.e., the f /β ratio), the forcing time scale, and the factors influencing the stratification. In some parameter settings of the models, although not those that are most realistic for the earth’s atmosphere, it is possible to produce eddies that are considerably larger than the deformation scales and an inverse cascade in the barotropic flow with a −5/3 spectrum. The vertical flux of heat is found to be related to the isentropic slope.


Author(s):  
Sebastian Schimek ◽  
Bernhard Ćosić ◽  
Jonas P. Moeck ◽  
Steffen Terhaar ◽  
Christian Oliver Paschereit

The current paper investigates the nonlinear interaction of the flow field and the unsteady heat release rate and the role of swirl fluctuations. The test rig consists of a generic swirl-stabilized combustor fed with natural gas and equipped with a high-amplitude forcing device. The influence of the phase between axial and azimuthal velocity oscillations is assessed on the basis of the amplitude and phase relations between the velocity fluctuations at the inlet and the outlet of the burner. These relations are determined in the experiment with the multimicrophone-method and a two component laser Doppler velocimeter (LDV). Particle image velocimetry (PIV) and OH*-chemiluminescence measurements are conducted to study the interaction between the flow field and the flame. For several frequency regimes, characteristic properties of the forced flow field and flame are identified, and a strong amplitude dependence is observed. It is found that the convective time delay between the swirl generator and the flame has an important influence on swirl-number oscillations and the flame dynamics in the low-frequency regime. For mid and high frequencies, significant changes in the mean flow field and the mean flame position are identified for high forcing amplitudes. These affect the interaction between coherent structures and the flame and are suggested to be responsible for the saturation in the flame response at high forcing amplitudes.


2007 ◽  
Vol 572 ◽  
pp. 179-188 ◽  
Author(s):  
THOMAS BOECK ◽  
DMITRY KRASNOV ◽  
EGBERT ZIENICKE

Mean flow properties of turbulent magnetohydrodynamic channel flow with electrically insulating channel walls are studied using high-resolution direct numerical simulations. The Lorentz force due to the homogeneous wall-normal magnetic field is computed in the quasi-static approximation. For strong magnetic fields, the mean velocity profile shows a clear three-layer structure consisting of a viscous region near each wall and a plateau in the middle connected by logarithmic layers. This structure reflects the significance of viscous, turbulent, and electromagnetic stresses in the streamwise momentum balance dominating the viscous, logarithmic, and plateau regions, respectively. The width of the logarithmic layers changes with the ratio of Reynolds- and Hartmann numbers. Turbulent stresses typically decay more rapidly away from the walls than predicted by mixing-length models.


2020 ◽  
Vol 50 (1) ◽  
pp. 3-33
Author(s):  
Shih-Nan Chen ◽  
Chiou-Jiu Chen ◽  
James A. Lerczak

AbstractThis study examines the utility of Eady-type theories as applied to understanding baroclinic instability in coastal flows where depth variations and bottom drag are important. The focus is on the effects of nongeostrophy, boundary dissipation, and bottom slope. The approach compares theoretically derived instability properties against numerical model calculations, for experiments designed to isolate the individual effects and justified to have Eady-like basic states. For the nongeostrophic effect, the theory of Stone (1966) is shown to give reasonable predictions for the most unstable growth rate and wavelength. It is also shown that the growing instability in a fully nonlinear model can be interpreted as boundary-trapped Rossby wave interactions—that is, wave phase locking and westward phase tilt allow waves to be mutually amplified. The analyses demonstrate that both the boundary dissipative and bottom slope effects can be represented by vertical velocities at the lower boundary of the unstable interior, via inducing Ekman pumping and slope-parallel flow, respectively, as proposed by the theories of Williams and Robinson (1974; referred to as the Eady–Ekman problem) and Blumsack and Gierasch (1972). The vertical velocities, characterized by a friction parameter and a slope ratio, modify the bottom wave and thus the scale selection. However, the theories have inherent quantitative limitations. Eady–Ekman neglects boundary layer responses that limit the increase of bottom stress, thereby overestimating the Ekman pumping and growth rate reduction at large drag. Blumsack and Gierasch’s (1972) model ignores slope-induced horizontal shear in the mean flow that tilts the eddies to favor converting energy back to the mean, thus having limited utility over steep slopes.


1970 ◽  
Vol 40 (2) ◽  
pp. 273-306 ◽  
Author(s):  
Michael E. McIntyre

Perturbation series are developed and mathematically justified, using a straightforward perturbation formalism (that is more widely applicable than those given in standard textbooks), for the case of the two-dimensional inviscid Orr-Sommerfeld-like eigenvalue problem describing quasi-geostrophic wave instabilities of parallel flows in rotating stratified fluids.The results are first used to examine the instability properties of the perturbed Eady problem, in which the zonal velocity profile has the formu=z+ μu1(y,z) where, formally, μ [Lt ] 1. The connexion between baroclinic instability theories with and without short wave cutoffs is clarified. In particular, it is established rigorously that there is instability at short wavelengths in all cases for which such instability would be expected from the ‘critical layer’ argument of Bretherton. (Therefore the apparently conflicting results obtained earlier by Pedlosky are in error.)For the class of profiles of formu=z+ μu1(y) it is then shown from an examination of theO(μ) eigenfunction correction why, under certain conditions, growing baroclinic waves will always produce a counter-gradient horizontal eddy flux of zonal momentum tending to reinforce the horizontal shear of such profiles. Finally, by computing a sufficient number of the higher corrections, this first-order result is shown to remain true, and its relationship to the actual rate of change of the mean flow is also displayed, for a particular jet-like form of profile withfinitehorizontal shear. The latter detailed results may help to explain at least one interesting feature of the mean flow found in a recent numerical solution for the wave régime in a heated rotating annulus.


Sign in / Sign up

Export Citation Format

Share Document