Flow instabilities in the wide-gap spherical Couette system

2013 ◽  
Vol 738 ◽  
pp. 184-221 ◽  
Author(s):  
Johannes Wicht

AbstractThe spherical Couette system is a spherical shell filled with a viscous fluid. Flows are driven by the differential rotation between the inner and the outer boundary that rotate with $\Omega $ and $\Omega + \mathrm{\Delta} \Omega $ about a common axis. This setup has been proposed for second-generation dynamo experiments. We numerically explore the different instabilities emerging for rotation rates up to $\Omega = (1/ 3)\times 1{0}^{7} $, venturing also into the nonlinear regime where oscillatory and chaotic solutions are found. The results provide a comprehensive overview of the possible flow regimes. For low values of $\Omega $ viscosity dominates and an equatorial jet in meridional circulation and zonal flow develops that becomes unstable as the differential rotation is increased beyond a critical value. For intermediate $\Omega $ and an inner boundary rotating slower than the outer one, new double-roll and helical instabilities are found. For large $\Omega $ values Coriolis effects enforce a nearly two-dimensional fundamental flow where a Stewartson shear layer develops at the tangent cylinder. This shear layer is the source of nearly geostrophic non-axisymmetric instabilities that resemble columnar Rossby modes. At first, the instabilities differ significantly depending on whether the inner boundary rotates faster $( \mathrm{\Delta} \Omega \gt 0)$ or slower $( \mathrm{\Delta} \Omega \lt 0)$ than the outer one. For very large outer boundary rotation rates, however, both instabilities once more become comparable. Fast inertial waves similar to those observed in recent spherical Couette experiments prevail for larger $\Omega $ values and $ \mathrm{\Delta} \Omega \lt 0$ in when $ \mathrm{\Delta} \Omega $ and $\Omega $ are of comparable magnitude. For larger differential rotations $ \mathrm{\Delta} \Omega \gg \Omega $, however, the equatorial jet instability always takes over.

2018 ◽  
Vol 843 ◽  
pp. 211-243 ◽  
Author(s):  
A. Barik ◽  
S. A. Triana ◽  
M. Hoff ◽  
J. Wicht

The spherical Couette system, consisting of a viscous fluid between two differentially rotating concentric spheres, is studied using numerical simulations and compared with experiments performed at BTU Cottbus-Senftenberg, Germany. We concentrate on the case where the outer boundary rotates fast enough for the Coriolis force to play an important role in the force balance, and the inner boundary rotates slower or in the opposite direction as compared to the outer boundary. As the magnitude of differential rotation is increased, the system is found to transition through three distinct hydrodynamic regimes. The first regime consists of the emergence of the first non-axisymmetric instability. Thereafter one finds the onset of ‘fast’ equatorially antisymmetric inertial modes, with pairs of inertial modes forming triadic resonances with the first instability. A further increase in the magnitude of differential rotation leads to the flow transitioning to turbulence. Using an artificial excitation, we study how the background flow modifies the inertial mode frequency and structure, thereby causing departures from the eigenmodes of a full sphere and a spherical shell. We investigate triadic resonances of pairs of inertial modes with the fundamental instability. We explore possible onset mechanisms through numerical experiments.


2002 ◽  
Vol 452 ◽  
pp. 263-291 ◽  
Author(s):  
E. DORMY ◽  
D. JAULT ◽  
A. M. SOWARD

We consider axisymmetric magnetohydrodynamic motion in a spherical shell driven by rotating the inner boundary relative to the stationary outer boundary – spherical Couette flow. The inner solid sphere is rigid with the same electrical conductivity as the surrounding fluid; the outer rigid boundary is an insulator. A force-free dipole magnetic field is maintained by a dipole source at the centre. For strong imposed fields (as measured by the Hartmann number M), the numerical simulations of Dormy et al. (1998) showed that a super-rotating shear layer (with angular velocity about 50% above the angular velocity of the inner core) is attached to the magnetic field line [Cscr ] tangent to the outer boundary at the equatorial plane of symmetry. At large M, we obtain analytically the mainstream solution valid outside all boundary layers by application of Hartmann jump conditions across the inner- and outer-sphere boundary layers. We formulate the large-M boundary layer problem for the free shear layer of width M−1/2 containing [Cscr ] and solve it numerically. The super-rotation can be understood in terms of the nature of the meridional electric current flow in the shear layer, which is fed by the outer-sphere Hartmann layer. Importantly, a large fraction of the current entering the shear layer is tightly focused and effectively released from a point source at the equator triggered by the tangency of the [Cscr ]-line. The current injected by the source follows the [Cscr ]-line closely but spreads laterally due to diffusion. In consequence, a strong azimuthal Lorentz force is produced, which takes opposite signs either side of the [Cscr ]-line; order-unity super-rotation results on the equatorial side. In fact, the point source is the small equatorial Hartmann layer of radial width M−2/3 ([Lt ]M−1/2) and latitudinal extent M−1/3. We construct its analytic solution and so determine an inward displacement width O(M−2/3) of the free shear layer. We compare our numerical solution of the free shear layer problem with our numerical solution of the full governing equations for M in excess of 104. We obtain excellent agreement. Some of our more testing comparisons are significantly improved by incorporating the shear layer displacement caused by the equatorial Hartmann layer.


2010 ◽  
Vol 645 ◽  
pp. 145-185 ◽  
Author(s):  
A. M. SOWARD ◽  
E. DORMY

We consider the steady axisymmetric motion of an electrically conducting fluid contained within a spherical shell and permeated by a centred axial dipole magnetic field, which is strong as measured by the Hartmann number M. Slow axisymmetric motion is driven by rotating the inner boundary relative to the stationary outer boundary. For M ≫ 1, viscous effects are only important in Hartmann boundary layers adjacent to the inner and outer boundaries and a free shear-layer on the magnetic field line that is tangent to the outer boundary on the equatorial plane of symmetry. We measure the ability to leak electric current into the solid boundaries by the size of their relative conductance ɛ. Since the Hartmann layers are sustained by the electric current flow along them, the current inflow from the fluid mainstream needed to feed them increases in concert with the relative conductance, because of the increasing fraction ℒ of the current inflow leaked directly into the solids. Therefore the nature of the flow is sensitive to the relative sizes of ɛ−1 and M.The current work extends an earlier study of the case of a conducting inner boundary and an insulating outer boundary with conductance ɛo = 0 (Dormy, Jault & Soward, J. Fluid Mech., vol. 452, 2002, pp. 263–291) to other values of the outer boundary conductance. Firstly, analytic results are presented for the case of perfectly conducting inner and outer boundaries, which predict super-rotation rates Ωmax of order M1/2 in the free shear-layer. Successful comparisons are made with numerical results for both perfectly and finitely conducting boundaries. Secondly, in the case of a finitely conducting outer boundary our analytic results show that Ωmax is O(M1/2) for ɛo−1 ≪ 1 ≪ M3/4, O(ɛo2/3M1/2) for 1 ≪ ɛo−1 ≪ M3/4 and O(1) for 1 ≪ M3/4 ≪ ɛo−1. On increasing ɛo−1 from zero, substantial electric current leakage into the outer boundary, ℒo ≈ 1, occurs for ɛo−1 ≪ M3/4 with the shear-layer possessing the character appropriate to a perfectly conducting outer boundary. When ɛo−1 = O(M3/4) the current leakage is blocked near the equator, and the nature of the shear-layer changes. So, when M3/4 ≪ ɛo−1, the shear-layer has the character appropriate to an insulating outer boundary. More precisely, over the range M3/4 ≪ ɛo−1 ≪ M the blockage spreads outwards, reaching the pole when ɛo−1 = O(M). For M ≪ ɛo−1 current flow into the outer boundary is completely blocked, ℒo ≪ 1.


2020 ◽  
Vol 498 (3) ◽  
pp. 3758-3781 ◽  
Author(s):  
Adam S Jermyn ◽  
Shashikumar M Chitre ◽  
Pierre Lesaffre ◽  
Christopher A Tout

ABSTRACT We derive the scaling of differential rotation in both slowly and rapidly rotating convection zones using order of magnitude methods. Our calculations apply across stars and fluid planets and all rotation rates, as well as to both magnetized and purely hydrodynamic systems. We find shear |R∇Ω| of order the angular frequency Ω for slowly rotating systems with Ω ≪ |N|, where N is the Brünt–Väisälä frequency, and find that it declines as a power law in Ω for rapidly rotating systems with Ω ≫ |N|. We further calculate the meridional circulation rate and baroclinicity and examine the magnetic field strength in the rapidly rotating limit. Our results are in general agreement with simulations and observations and we perform a detailed comparison with those in a companion paper.


2016 ◽  
Vol 20 (4) ◽  
pp. 944-968
Author(s):  
Kai Zhang ◽  
Ming Li ◽  
Jingzhi Li

AbstractRemoving geometric details from the computational domain can significantly reduce the complexity of downstream task of meshing and simulation computation, and increase their stability. Proper estimation of the sensitivity analysis error induced by removing such domain details, called defeaturing errors, can ensure that the sensitivity analysis fidelity can still be met after simplification. In this paper, estimation of impacts of removing arbitrarily constrained domain details to the analysis of incompressible fluid flows is studied with applications to fast analysis of incompressible fluid flows in complex environments. The derived error estimator is applicable to geometric details constrained by either Dirichlet or Neumann boundary conditions, and has no special requirements on the outer boundary conditions. Extensive numerical examples were presented to demonstrate the effectiveness and efficiency of the proposed error estimator.


1987 ◽  
Vol 178 ◽  
pp. 521-534 ◽  
Author(s):  
D. R. Fearn ◽  
M. R. E. Proctor

We investigate the dynamical consequences of an axisymmetric velocity field with a poloidal magnetic field driven by a prescribed e.m.f. E. The problem is motivated by previous investigations of dynamically driven dynamos in the magnetostrophic range. A geostrophic zonal flow field is added to a previously described velocity, and determined by the requirement that Taylor's constraint (Taylor 1963) (guaranteeing dynamical self-consistency of the fields) be satisfied. Several solutions are exhibited, and it is suggested that self-consistent solutions can always be found to this ‘forced’ problem, whereas the usual α-effect dynamo formalism in which E is a linear function of the magnetic field leads to a difficult transcendentally nonlinear characteristic value problem that may not always possess solutions.


1972 ◽  
Vol 53 (4) ◽  
pp. 647-655 ◽  
Author(s):  
M. R. Foster

The flow induced by the differential rotation of a cylindrical depression of radius a in one of two parallel rigid planes rapidly rotating about their common normal at speed Q is studied. A Taylor column bounded by the usual Stewartson layers arises, but the shear-layer structure is rather different from any previously studied. The Ei-layers (E = v/ωa2) smooth the discontinuity in the geostrophic flow, but the way in which this is accomplished is related to the possible singu-larities of the E1/3-layer solutions. The fact that the 1/4-layer is partially free and partially attached to a vertical boundary accounts for the new joining conditions for the 1/4-layer. The drag on a right circular cylindrical bump in uniform flow is given in addition to some general comments on the applicability of these joining conditions to the motion of an axisymmetric object of quite general shape.


Fractals ◽  
2011 ◽  
Vol 19 (03) ◽  
pp. 311-316 ◽  
Author(s):  
ROCCO PORTARO ◽  
MOHAMED FAYED ◽  
AMY-LEE GUNTER ◽  
HAMID AIT ABDERRAHMANE ◽  
HOI DICK NG

In this study, we illustrate the fractal nature of the wake shed by a periodically flapping filament. Such wake structure is a combination of primary vortex shedding resulting in the von Kármán vortex street, a series of concentrated vortex dipoles formed when the trailing edges of filaments reach their maximum amplitudes and small eddies form along the shear layer connected with the concentrated vortices due to the shear layer instability. The vortex dynamics of the flapping filament are visualized and imaged experimentally using a soap-film flow tunnel with a high-speed camera and a low pressure sodium lamp as a light source. The wake fractal geometry is measured using the standard box-counting method and it is shown that the fractal dimension of the soap pattern boundaries in the wake is D = 1.38 ± 0.05, which agrees well with those measured for fully developed turbulences and other shear flow phenomena. The invariant of the fractality in the wake induced by the flapping filament thus provides another illustration of the geometrical self-similarity and nonlinear dynamics of chaotic fluid flows.


2009 ◽  
Vol 641 ◽  
pp. 335-358 ◽  
Author(s):  
CHRISTOPHER J. DAVIES ◽  
DAVID GUBBINS ◽  
PETER K. JIMACK

We investigate thermally driven convection in a rotating spherical shell subject to inhomogeneous heating on the outer boundary, extending previous results to more rapid rotation rates and larger amplitudes of the boundary heating. The analysis explores the conditions under which steady flows can be obtained, and the stability of these solutions, for two boundary heating modes: first, when the scale of the boundary heating corresponds to the most unstable mode of the homogeneous problem; second, when the scale is larger. In the former case stable steady solutions exhibit a two-layer flow pattern at moderate rotation rates, but at very rapid rotation rates no steady solutions exist. In the latter case, stable steady solutions are always possible, and unstable solutions show convection rolls that cluster into nests that are out of phase with the boundary anomalies and remain trapped for many thermal diffusion times.


2008 ◽  
Vol 604 ◽  
pp. 175-197 ◽  
Author(s):  
DENYS SCHMITT ◽  
T. ALBOUSSIÈRE ◽  
D. BRITO ◽  
P. CARDIN ◽  
N. GAGNIÈRE ◽  
...  

The magnetostrophic regime, in which Lorentz and Coriolis forces are in balance, has been investigated in a rapidly rotating spherical Couette flow experiment. The spherical shell is filled with liquid sodium and permeated by a strong imposed dipolar magnetic field. Azimuthally travelling hydromagnetic waves have been put in evidence through a detailed analysis of electric potential differences measured on the outer sphere, and their properties have been determined. Several types of wave have been identified depending on the relative rotation rates of the inner and outer spheres: they differ by their dispersion relation and by their selection of azimuthal wavenumbers. In addition, these waves constitute the largest contribution to the observed fluctuations, and all of them travel in the retrograde direction in the frame of reference bound to the fluid. We identify these waves as magneto-inertial waves by virtue of the close proximity of the magnetic and inertial characteristic time scales of relevance in our experiment.


Sign in / Sign up

Export Citation Format

Share Document