No Arabic abstract
The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells -- a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed.
Double-diffusive convection driven by both thermal and compositional buoyancy in a rotating spherical shell can exhibit a rather large number of behaviours often distinct from that of the single diffusive system. In order to understand how the differences in thermal and compositional molecular diffusivities determine the dynamics of thermo-compositional convection we investigate numerically the linear onset of convective instability in a double-diffusive setup. We construct an alternative equivalent formulation of the non-dimensional equations where the linearised double-diffusive problem is described by an effective Rayleigh number, $text{Ra}$, measuring the amplitude of the combined buoyancy driving, and a second parameter, $alpha$, measuring the mixing of the thermal and compositional contributions. This formulation is useful in that it allows for the analysis of several limiting cases and reveals dynamical similarities in the parameters space which are not obvious otherwise. We analyse the structure of the critical curves in this $text{Ra}-alpha$ space, explaining asymptotic behaviours in $alpha$, transitions between inertial and diffusive regimes, and transitions between large scale (fast drift) and small scale (slow drift) convection. We perform this analysis for a variety of diffusivities, rotation rates and shell aspect ratios showing where and when new modes of convection take place.
Convection and magnetic field generation in the Earth and planetary interiors are driven by both thermal and compositional gradients. In this work numerical simulations of finite-amplitude double-diffusive convection and dynamo action in rapidly rotating spherical shells full of incompressible two-component electrically-conducting fluid are reported. Four distinct regimes of rotating double-diffusive convection identified in a recent linear analysis (Silva et al., 2019, Geophys. Astrophys. Fluid Dyn., doi:10.1080/03091929.2019.1640875) are found to persist significantly beyond the onset of instability while their regime transitions remain abrupt. In the semi-convecting and the fingering regimes characteristic flow velocities are small compared to those in the thermally- and compositionally-dominated overturning regimes, while zonal flows remain weak in all regimes apart from the thermally-dominated one. Compositionally-dominated overturning convection exhibits significantly narrower azimuthal structures compared to all other regimes while differential rotation becomes the dominant flow component in the thermally-dominated case as driving is increased. Dynamo action occurs in all regimes apart from the regime of fingering convection. While dynamos persist in the semi-convective regime they are very much impaired by small flow intensities and very weak differential rotation in this regime which makes poloidal to toroidal field conversion problematic. The dynamos in the thermally-dominated regime include oscillating dipolar, quadrupolar and multipolar cases similar to the ones known from earlier parameter studies. Dynamos in the compositionally-dominated regime exhibit subdued temporal variation and remain predominantly dipolar due to weak zonal flow in this regime. These results significantly enhance our understanding of the primary drivers of planetary core flows and magnetic fields.
We consider the effect of stratification on systematic, large-scale flows generated in anelastic convection. We present results from three-dimensional numerical simulations of convection in a rotating plane layer in which the angle between the axis of rotation and gravity is allowed to vary. This model is representative of different latitudes of a spherical body. We consider two distinct parameter regimes: (i) weakly rotating and (ii) rapidly rotating. In each case, we examine the effect of stratification on the flow structure and heat transport properties focussing on the difference between Boussinesq and anelastic convection. Furthermore, we show that regimes (i) and (ii) generate very different large-scale flows and we investigate the role stratification has in modifying these flows. The stratified flows possess a net helicity not present in the Boussinesq cases which we suggest, when combined with the self-generated shear flows, could be important for dynamo action.
We investigate the linear properties of the steady and axisymmetric stress-driven spin-down flow of a viscous fluid inside a spherical shell, both within the incompressible and anelastic approximations, and in the asymptotic limit of small viscosities. From boundary layer analysis, we derive an analytical geostrophic solution for the 3D incompressible steady flow, inside and outside the cylinder $mathcal{C}$ that is tangent to the inner shell. The Stewartson layer that lies on $mathcal{C}$ is composed of two nested shear layers of thickness $O(E^{2/7})$ and $O(E^{1/3})$. We derive the lowest order solution for the $E^{2/7}$-layer. A simple analysis of the $E^{1/3}$-layer laying along the tangent cylinder, reveals it to be the site of an upwelling flow of amplitude $O(E^{1/3})$. Despite its narrowness, this shear layer concentrates most of the global meridional kinetic energy of the spin-down flow. Furthermore, a stable stratification does not perturb the spin-down flow provided the Prandtl number is small enough. If this is not the case, the Stewartson layer disappears and meridional circulation is confined within the thermal layers. The scalings for the amplitude of the anelastic secondary flow have been found to be the same as for the incompressible flow in all three regions, at the lowest order. However, because the velocity no longer conforms the Taylor-Proudman theorem, its shape differs outside the tangent cylinder $mathcal{C}$, that is, where differential rotation takes place. Finally, we find the settling of the steady-state to be reached on a viscous time for the weakly, strongly and thermally unstratified incompressible flows. Large density variations relevant to astro- and geophysical systems, tend to slightly shorten the transient.
We theoretically investigate the effect of random fluctuations on the motion of elongated microswimmers near hydrodynamic transport barriers in externally-driven fluid flows. Focusing on the two-dimensional hyperbolic flow, we consider the effects of translational and rotational diffusion as well as tumbling, i.e. sudden jumps in the swimmer orientation. Regardless of whether diffusion or tumbling are the primary source of fluctuations, we find that noise significantly increases the probability that a swimmer crosses one-way barriers in the flow, which block the swimmer from returning to its initial position. We employ an asymptotic method for calculating the probability density of noisy swimmer trajectories in a given fluid flow, which produces solutions to the time-dependent Fokker-Planck equation in the weak-noise limit. This procedure mirrors the semiclassical approximation in quantum mechanics and similarly involves calculating the least-action paths of a Hamiltonian system derived from the swimmers Fokker-Planck equation. Using the semiclassical technique, we compute (i) the steady-state orientation distribution of swimmers with rotational diffusion and tumbling and (ii) the probability that a diffusive swimmer crosses a one-way barrier. The semiclassical results compare favorably with Monte Carlo calculations.