5 results
From three-dimensional to quasi-two-dimensional: transient growth in magnetohydrodynamic duct flows
- Oliver G. W. Cassells, Tony Vo, Alban Pothérat, Gregory J. Sheard
-
- Journal:
- Journal of Fluid Mechanics / Volume 861 / 25 February 2019
- Published online by Cambridge University Press:
- 19 December 2018, pp. 382-406
-
- Article
- Export citation
-
This study seeks to elucidate the linear transient growth mechanisms in a uniform duct with square cross-section applicable to flows of electrically conducting fluids under the influence of an external magnetic field. A particular focus is given to the question of whether at high magnetic fields purely two-dimensional mechanisms exist, and whether these can be described by a computationally inexpensive quasi-two-dimensional model. Two Reynolds numbers of $5000$ and $15\,000$ and an extensive range of Hartmann numbers $0\leqslant \mathit{Ha}\leqslant 800$ were investigated. Three broad regimes are identified in which optimal mode topology and non-modal growth mechanisms are distinct. These regimes, corresponding to low, moderate and high magnetic field strengths, are found to be governed by the independent parameters; Hartmann number, Reynolds number based on the Hartmann layer thickness $R_{H}$ and Reynolds number built upon the Shercliff layer thickness $R_{S}$, respectively. Transition between regimes respectively occurs at $\mathit{Ha}\approx 2$ and no lower than $R_{H}\approx 33.\dot{3}$. Notably for the high Hartmann number regime, quasi-two-dimensional magnetohydrodynamic models are shown to be excellent predictors of not only transient growth magnitudes, but also the fundamental growth mechanisms of linear disturbances. This paves the way for a precise analysis of transition to quasi-two-dimensional turbulence at much higher Hartmann numbers than is currently achievable.
Stability of the wakes of cylinders with triangular cross-sections
- Zhi Y. Ng, Tony Vo, Gregory J. Sheard
-
- Journal:
- Journal of Fluid Mechanics / Volume 844 / 10 June 2018
- Published online by Cambridge University Press:
- 12 April 2018, pp. 721-745
-
- Article
- Export citation
-
The stability of the wakes of cylinders with triangular cross-sections at incidence is investigated using Floquet stability analysis to elucidate the effects of cylinder inclination on the dominant flow instability. The upper limit of the Reynolds numbers (scaled by the height projected by the cylinder in this study) at which the wake of the two-dimensional base flow is time periodic is $Re\approx 140$ for most cylinder inclinations, exceeding which the flow becomes aperiodic, restricting the range of Reynolds numbers permitted for the stability analysis. Two different instability modes are predicted to manifest as the first-occurring mode at various cylinder inclinations – a regular mode possessing perturbation structures consistent with mode A dominates the wakes of cylinders at inclinations $\unicode[STIX]{x1D6FC}\lesssim 34.6^{\circ }$ and $\unicode[STIX]{x1D6FC}\gtrsim 55.4^{\circ }$, with a subharmonic mode consistent with mode C emerging as the primary mode in the wakes of the cylinder at the intermediate range of inclinations. For all inclinations, the mode B branch is not detected within the range of Reynolds numbers examined. The peak instability growth rates corresponding to mode A for all cylinder inclinations describe a linear variation with $(Re-Re_{A})/Re_{A}$, where $Re_{A}$ is the mode A transition Reynolds number, while those corresponding to mode C vary only approximately linearly. The generalized trend most pertinently shows mode C to develop more rapidly than mode A at inclinations which permit it. Examination of the near wake of the two-dimensional time-periodic base flow demonstrates the dependence of the development and intensity of mode C on imbalances in the flow solution over each shedding period, directly implying that the two-dimensional base flow solutions deviate from the half-period-flip map as the cylinder inclination is increased. The degree of asymmetry of the two-dimensional base flow relative to the ideal half-period-flip map is quantified using several measures. The results show distinctly different trends in these asymmetry measures between inclinations where modes A or C are dominant, agreeing with results from the stability analysis. The nature of the predicted instability modes at transition are also investigated by applying the Stuart–Landau equation, showing the transitions to be supercritical for all cylinder inclinations, with mode C being consistently more strongly supercritical than mode A.
Non-axisymmetric flows in a differential-disk rotating system
- Tony Vo, Luca Montabone, Peter L. Read, Gregory J. Sheard
-
- Journal:
- Journal of Fluid Mechanics / Volume 775 / 25 July 2015
- Published online by Cambridge University Press:
- 25 June 2015, pp. 349-386
-
- Article
- Export citation
-
The non-axisymmetric structure of an unstable Stewartson shear layer generated via a differential rotation between flush disks and a cylindrical enclosure is investigated numerically using both three-dimensional direct numerical simulation and a quasi-two-dimensional model. Previous literature has only considered the depth-independent quasi-two-dimensional model due to its low computational cost. The three-dimensional model implemented here highlights the supercritical instability responsible for the polygonal deformation of the shear layer in the linear and nonlinear growth regimes and reveals that linear stability analysis is capable of accurately determining the preferred azimuthal wavenumber for flow conditions near the onset of instability. This agreement is lost for sufficiently forced flows where nonlinear effects encourage the coalescence of vortices towards lower-wavenumber structures. Time-dependent flows are found for large Reynolds numbers defined based on the Stewartson layer thickness and azimuthal velocity differential. However, this temporal behaviour is not solely characterized by Reynolds number but is rather a function of both the Rossby and Ekman numbers. At high Ekman and Rossby numbers, unsteady flow emerges through a small-scale azimuthal destabilization of the axial jets within the Stewartson layers; at low Ekman numbers, unsteady flow emerges through a modulation in the strength of one of the axial vortices rolled up by non-axisymmetric instability of the Stewartson layer.
Effect of enclosure height on the structure and stability of shear layers induced by differential rotation
- Tony Vo, Luca Montabone, Gregory J. Sheard
-
- Journal:
- Journal of Fluid Mechanics / Volume 765 / 25 February 2015
- Published online by Cambridge University Press:
- 15 January 2015, pp. 45-81
-
- Article
- Export citation
-
The 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.
Linear stability analysis of a shear layer induced by differential coaxial rotation within a cylindrical enclosure
- Tony Vo, Luca Montabone, Gregory J. Sheard
-
- Journal:
- Journal of Fluid Mechanics / Volume 738 / 10 January 2014
- Published online by Cambridge University Press:
- 05 December 2013, pp. 299-334
-
- Article
- Export citation
-
The 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.