Effect of enclosure height on the structure and stability of shear layers induced by differential rotation

2015 ◽  
Vol 765 ◽  
pp. 45-81 ◽  
Author(s):  
Tony Vo ◽  
Luca Montabone ◽  
Gregory J. Sheard

AbstractThe structure and stability of Stewartson shear layers with different heights are investigated numerically via axisymmetric simulation and linear stability analysis, and a validation of the quasi-two-dimensional model is performed. The shear layers are generated in a rotating cylindrical tank with circular disks located at the lid and base imposing a differential rotation. The axisymmetric model captures both the thick and thin nested Stewartson layers, which are scaled by the Ekman number ($\mathit{E}\,$) as $\mathit{E}\,^{1/4}$ and $\mathit{E}\,^{1/3}$ respectively. In contrast, the quasi-two-dimensional model only captures the $\mathit{E}\,^{1/4}$ layer as the axial velocity required to invoke the $\mathit{E}\,^{1/3}$ layer is excluded. A direct comparison between the axisymmetric base flows and their linear stability in these two models is examined here for the first time. The base flows of the two models exhibit similar flow features at low Rossby numbers ($\mathit{Ro}$), with differences evident at larger $\mathit{Ro}$ where depth-dependent features are revealed by the axisymmetric model. Despite this, the quasi-two-dimensional model demonstrates excellent agreement with the axisymmetric model in terms of the shear-layer thickness and predicted stability. A study of various aspect ratios reveals that a Reynolds number based on the theoretical Ekman layer thickness is able to describe the transition of a base flow that is reflectively symmetric about the mid-plane to a symmetry-broken state. Additionally, the shear-layer thicknesses scale closely to the expected ${\it\delta}_{vel}\propto A\mathit{E}\,^{1/4}$ and ${\it\delta}_{vort}\propto A\mathit{E}\,^{1/3}$ for shear layers that are not affected by the confinement ($A\mathit{E}\,^{1/4}\lesssim 0.34$ in this system, the ratio of tank height to shear-layer radius). The linear stability analysis reveals that the ratio of Stewartson layer radius to thickness should be greater than $45$ for the stability of the flow to be independent of aspect ratio. Thus, for sufficiently small $A\mathit{E}\,^{1/4}$ and $A\mathit{E}\,^{1/3}$, the flow characteristics remain similar and the linear stability of the flow can be described universally when the azimuthal wavelength is scaled against $A$. The analysis also recovers an asymptotic scaling for the normalized azimuthal wavelength which suggests that ${\it\lambda}_{{\it\theta},c}^{\ast }\propto (|\mathit{Ro}|/\mathit{E}\,^{2})^{-1/5}$ for geometry-independent shear layers at marginal stability.

1985 ◽  
Vol 157 ◽  
pp. 449-470 ◽  
Author(s):  
P. Blondeaux ◽  
G. Seminara

A two-dimensional model of flow and bed topography in sinuous channels with erodible boundaries is developed and applied in order to investigate the mechanism of meander initiation. By reexamining the problem recently tackled by Ikeda, Parker & Sawai (1981), a previously undiscovered ‘resonance’ phenomenon is detected which occurs when the values of the relevant parameters fall within a neighbourhood of certain critical values. It is suggested that the above resonance controls the bend growth, and it is shown that it is connected in some sense with bar instability. In fact, by performing a linear stability analysis of flow in straight erodible channels, resonant flow in sinuous channels is shown to occur when curvature ‘forces’ a ‘natural’ solution represented by approximately steady perturbations of the alternate bar type. A comparison with experimental observations appears to support the idea that resonance is associated with meander formation.


2013 ◽  
Vol 738 ◽  
pp. 299-334 ◽  
Author(s):  
Tony Vo ◽  
Luca Montabone ◽  
Gregory J. Sheard

AbstractThe generation of distinct polygonal configurations via the instability of a Stewartson shear layer is numerically investigated. The shear layer is induced using a rotating cylindrical tank with differentially forced disks located at the top and bottom boundaries. The incompressible Navier–Stokes equations are solved on a two-dimensional semi-meridional plane. Axisymmetric base flows are consistently found to reach a steady state for a wide range of flow conditions, and details of the vertical structure are revealed. An axially invariant two-dimensional flow is ascertained for small $\vert \mathit{Ro}\vert $, which substantiates the Taylor–Proudman theorem. Sufficient increases in $\vert \mathit{Ro}\vert $ forcing develops flow features that break this quasi-two-dimensionality. The onset of this breaking occurs earlier with increasing $\vert \mathit{Ro}\vert $ for $\mathit{Ro}\gt 0$ compared with $\mathit{Ro}\lt 0$. The thickness scaling of the vertical Stewartson layers are in agreement with previous analytical results. Growth rates of the most unstable azimuthal wavenumber from a global linear stability analysis are obtained. The threshold between axisymmetric and non-axisymmetric flow follows a power law, and both positive- and negative-$\mathit{Ro}$ regimes are found to adopt the same threshold for instability, namely $\vert \mathit{Ro}\vert \geq 18. 1{E}^{0. 767} $. This relationship corresponds to a constant critical internal Reynolds number of ${\mathit{Re}}_{i, c} \simeq 22. 5$. A review of reported critical internal Reynolds number and their characteristic length scales yields a consistent instability onset given by $\vert \mathit{Ro}\vert / {E}^{3/ 4} = 15. 4{\unicode{x2013}} 16. 6$; here we find $\vert \mathit{Ro}\vert / {E}^{3/ 4} = 15. 8$. At the onset of linear instability, the initially circular shear layer deforms, resulting in a polygonal structure consistent with barotropic instability. Dominant azimuthal wavenumbers range from $3$ to $7$ at the onset of instability for the parameter space explored. Empirical relationships for the preferential wavenumber have been obtained. Additional instability modes have been discovered that favour higher wavenumbers, and these exhibit structures localized to the disk–tank interfaces.


Author(s):  
Mark Frederick ◽  
Kiran Manoharan ◽  
Joshua Dudash ◽  
Brian Brubaker ◽  
Santosh Hemchandra ◽  
...  

Combustion instability, the coupling between flame heat release rate oscillations and combustor acoustics, is a significant issue in the operation of gas turbine combustors. This coupling is often driven by oscillations in the flow field. Shear layer roll-up, in particular, has been shown to drive longitudinal combustion instability in a number of systems, including both laboratory and industrial combustors. One method for suppressing combustion instability would be to suppress the receptivity of the shear layer to acoustic oscillations, severing the coupling mechanism between the acoustics and the flame. Previous work suggested that the existence of a precessing vortex core (PVC) may suppress the receptivity of the shear layer, and the goal of this study is to first, confirm that this suppression is occurring, and second, understand the mechanism by which the PVC suppresses the shear layer receptivity. In this paper, we couple experiment with linear stability analysis to determine whether a PVC can suppress shear layer receptivity to longitudinal acoustic modes in a nonreacting swirling flow at a range of swirl numbers. The shear layer response to the longitudinal acoustic forcing manifests as an m = 0 mode since the acoustic field is axisymmetric. The PVC has been shown both in experiment and linear stability analysis to have m = 1 and m = −1 modal content. By comparing the relative magnitude of the m = 0 and m = −1,1 modes, we quantify the impact that the PVC has on the shear layer response. The mechanism for shear layer response is determined using companion forced response analysis, where the shear layer disturbance growth rates mirror the experimental results. Differences in shear layer thickness and azimuthal velocity profiles drive the suppression of the shear layer receptivity to acoustic forcing.


2019 ◽  
Vol 9 (22) ◽  
pp. 4840
Author(s):  
Yue Chen

This paper starts with a generalized Burton, Cabrera and Frank (BCF) model by considering the energetic contribution of the adjacent terraces to the step chemical potential. We use the linear stability analysis of the quasistatic free-boundary problem for a two-dimensional step separated by broad terraces to study the step-meandering instabilities. The results show that the equilibrium adatom coverage has influence on the morphological instabilities.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


2009 ◽  
Vol 622 ◽  
pp. 291-320 ◽  
Author(s):  
M. D. GRIFFITH ◽  
T. LEWEKE ◽  
M. C. THOMPSON ◽  
K. HOURIGAN

Pulsatile inlet flow through a circular tube with an axisymmetric blockage of varying size is studied both numerically and experimentally. The geometry consists of a long, straight tube and a blockage, semicircular in cross-section, serving as a simplified model of an arterial stenosis. The stenosis is characterized by a single parameter, the aim being to highlight fundamental behaviours of constricted pulsatile flows. The Reynolds number is varied between 50 and 700 and the stenosis degree by area between 0.20 and 0.90. Numerically, a spectral element code is used to obtain the axisymmetric base flow fields, while experimentally, results are obtained for a similar set of geometries, using water as the working fluid. For low Reynolds numbers, the flow is characterized by a vortex ring which forms directly downstream of the stenosis, for which the strength and downstream propagation velocity vary with the stenosis degree. Linear stability analysis is performed on the simulated axisymmetric base flows, revealing a range of absolute instability modes. Comparisons are drawn between the numerical linear stability analysis and the observed instability in the experimental flows. The observed flows are less stable than the numerical analysis predicts, with convective shear layer instability present in the experimental flows. Evidence is found of Kelvin–Helmholtz-type shear layer roll-ups; nonetheless, the possibility of the numerically predicted absolute instability modes acting in the experimental flow is left open.


2010 ◽  
Vol 643 ◽  
pp. 363-394 ◽  
Author(s):  
M. RIEUTORD ◽  
L. VALDETTARO

We investigate the properties of forced inertial modes of a rotating fluid inside a spherical shell. Our forcing is tidal like, but its main property is that it is on the large scales. By numerically solving the linear equations of this problem, including viscosity, we first confirm some analytical results obtained on a two-dimensional model by Ogilvie (J. Fluid Mech., vol. 543, 2005, p. 19); some additional properties of this model are uncovered like the existence of narrow resonances associated with periodic orbits of characteristics. We also note that as the frequency of the forcing varies, the dissipation varies drastically if the Ekman number E is low (as is usually the case). We then investigate the three-dimensional case and compare the results to the foregoing model. The three-dimensional solutions show, like their two-dimensional counterpart, a spiky dissipation curve when the frequency of the forcing is varied; they also display small frequency intervals where the viscous dissipation is independent of viscosity. However, we show that the response of the fluid in these frequency intervals is crucially dominated by the shear layer that is emitted at the critical latitude on the inner sphere. The asymptotic regime, where the dissipation is independent of the viscosity, is reached when an attractor has been excited by this shear layer. This property is not shared by the two-dimensional model where shear layers around attractors are independent of those emitted at the critical latitude. Finally, resonances of the three-dimensional model correspond to some selected least damped eigenmodes. Unlike their two-dimensional counter parts these modes are not associated with simple attractors; instead, they show up in frequency intervals with weakly contracting webs of characteristics. Besides, we show that the inner core is negligible when its relative radius is less than the critical value 0.4E1/5. For these spherical shells, the full sphere solutions give a good approximation of the flows.


Sign in / Sign up

Export Citation Format

Share Document