Bistable flows in precessing spheroids

Precession driven flows are found in any rotating container filled with liquid, when the rotation axis itself rotates about a secondary axis that is fixed in an inertial frame of reference. Because of its relevance for planetary fluid layers, many works consider spheroidal containers, where the uniform vorticity component of the bulk flow is reliably given by the well-known equations obtained by Busse in 1968. So far however, no analytical result on the solutions is available. Moreover, the cases where multiple flows can coexist have not been investigated in details since their discovery by Noir et al. (2003). In this work, we aim at deriving analytical results on the solutions, aiming in particular at, first estimating the ranges of parameters where multiple solutions exist, and second studying quantitatively their stability. Using the models recently proposed by Noir \& C{\'e}bron (2013), which are more generic in the inviscid limit than the equations of Busse, we analytically describe these solutions, their conditions of existence, and their stability in a systematic manner. We then successfully compare these analytical results with the theory of Busse (1968). Dynamical model equations are finally proposed to investigate the stability of the solutions, which allows to describe the bifurcation of the unstable flow solution. We also report for the first time the possibility that time-dependent multiple flows can coexist in precessing triaxial ellipsoids. Numerical integrations of the algebraic and differential equations have been efficiently performed with the dedicated script FLIPPER (supplementary material)

Precession-driven flows in non-axisymmetric ellipsoids

We study the flow forced by precession in rigid non-axisymmetric ellipsoidal containers. To do so, we revisit the inviscid and viscous analytical models that have been previously developed for the spheroidal geometry by, respectively, Poincaré (Bull. Astronomique, vol. XXVIII, 1910, pp. 1-36) and Busse (J. Fluid Mech., vol. 33, 1968, pp. 739-751), and we report the first numerical simulations of flows in such a geometry. In strong contrast with axisymmetric spheroids, where the forced flow is systematically stationary in the precessing frame, we show that the forced flow is unsteady and periodic. Comparisons of the numerical simulations with the proposed theoretical model show excellent agreement for both axisymmetric and non-axisymmetric containers. Finally, since the studied configuration corresponds to a tidally locked celestial body such as the Earth's Moon, we use our model to investigate the challenging but planetary-relevant limit of very small Ekman numbers and the particular case of our Moo

Libration driven multipolar instabilities

We consider rotating flows in non-axisymmetric enclosures that are driven by libration, i.e. by a small periodic modulation of the rotation rate. Thanks to its simplicity, this model is relevant to various contexts, from industrial containers (with small oscillations of the rotation rate) to fluid layers of terrestial planets (with length-of-day variations). Assuming a multipolar $n$-fold boundary deformation, we first obtain the two-dimensional basic flow. We then perform a short-wavelength local stability analysis of the basic flow, showing that an instability may occur in three dimensions. We christen it the Libration Driven Multipolar Instability (LDMI). The growth rates of the LDMI are computed by a Floquet analysis in a systematic way, and compared to analytical expressions obtained by perturbation methods. We then focus on the simplest geometry allowing the LDMI, a librating deformed cylinder. To take into account viscous and confinement effects, we perform a global stability analysis, which shows that the LDMI results from a parametric resonance of inertial modes. Performing numerical simulations of this librating cylinder, we confirm that the basic flow is indeed established and report the first numerical evidence of the LDMI. Numerical results, in excellent agreement with the stability results, are used to explore the non-linear regime of the instability (amplitude and viscous dissipation of the driven flow). We finally provide an example of LDMI in a deformed spherical container to show that the instability mechanism is generic. Our results show that the previously studied libration driven elliptical instability simply corresponds to the particular case $n=2$ of a wider class of instabilities. Summarizing, this work shows that any oscillating non-axisymmetric container in rotation may excite intermittent, space-filling LDMI flows, and this instability should thus be easy to observe experimentally

Libration-driven multipolar instabilities

We consider rotating flows in non-axisymmetric enclosures that are driven by libration, i.e. by a small periodic modulation of the rotation rate. Thanks to its simplicity, this model is relevant to various contexts, from industrial containers (with small oscillations of the rotation rate) to fluid layers of terrestrial planets (with length-of-day variations). Assuming a multipolar $n$ -fold boundary deformation, we first obtain the two-dimensional basic flow. We then perform a short-wavelength local stability analysis of the basic flow, showing that an instability may occur in three dimensions. We christen it the libration-driven multipolar instability (LDMI). The growth rates of the LDMI are computed by a Floquet analysis in a systematic way, and compared to analytical expressions obtained by perturbation methods. We then focus on the simplest geometry allowing the LDMI, a librating deformed cylinder. To take into account viscous and confinement effects, we perform a global stability analysis, which shows that the LDMI results from a parametric resonance of inertial modes. Performing numerical simulations of this librating cylinder, we confirm that the basic flow is indeed established and report the first numerical evidence of the LDMI. Numerical results, in excellent agreement with the stability results, are used to explore the nonlinear regime of the instability (amplitude and viscous dissipation of the driven flow). We finally provide an example of LDMI in a deformed spherical container to show that the instability mechanism is generic. Our results show that the previously studied libration-driven elliptical instability simply corresponds to the particular case $n= 2$ of a wider class of instabilities. Summarizing, this work shows that any oscillating non-axisymmetric container in rotation may excite intermittent, space-filling LDMI flows, and this instability should thus be easy to observe experimentall

Shallow water waves generated by subaerial solid landslides

Subaerial landslides are common events, which may generate very large water waves. The numerical modelling and simulation of these events are thus of primary interest for forecasting and mitigation of tsunami disasters. In this paper, we aim at describing these extreme events using a simplified shallow water model to derive relevant scaling laws. To cope with the problem, two different numerical codes are employed: one, SPHysics, is based on a Lagrangian meshless approach to accurately describe the impact stage whereas the other, Gerris, based on a two-phase finite-volume method is used to study the propagation of the wave. To validate Gerris for this very particular problem, two numerical cases of the literature are reproduced: a vertical sinking box and a 2-D wedge sliding down a slope. Then, to get insights into the problem of subaerial landslide-generated tsunamis and to further validate the codes for this case of landslides, a series of experiments is conducted in a water wave tank and successfully compared with the results of both codes. Based on a simplified approach, we derive different scaling laws in excellent agreement with the experiments and numerical simulation

Inviscid instabilities in rotating ellipsoids on eccentric Kepler orbits

We consider the hydrodynamic stability of homogeneous, incompressible and rotating ellipsoidal fluid masses. The latter are the simplest models of fluid celestial bodies with internal rotation and subjected to tidal forces. The classical problem is the stability of Roche–Riemann ellipsoids moving on circular Kepler orbits. However, previous stability studies have to be reassessed. Indeed, they only consider global perturbations of large wavelength or local perturbations of short wavelength. Moreover many planets and stars undergo orbital motions on eccentric Kepler orbits, implying time-dependent ellipsoidal semi-axes. This time dependence has never been taken into account in hydrodynamic stability studies. In this work we overcome these stringent assumptions. We extend the hydrodynamic stability analysis of rotating ellipsoids to the case of eccentric orbits. We have developed two open-source and versatile numerical codes to perform global and local inviscid stability analyses. They give sufficient conditions for instability. The global method, based on an exact and closed Galerkin basis, handles rigorously global ellipsoidal perturbations of unprecedented complexity. Tidally driven and libration-driven elliptical instabilities are first recovered and unified within a single framework. Then we show that new global fluid instabilities can be triggered in ellipsoids by tidal effects due to eccentric Kepler orbits. Their existence is confirmed by a local analysis and direct numerical simulations of the fully nonlinear and viscous problem. Thus a non-zero orbital eccentricity may have a strong destabilising effect in celestial fluid bodies, which may lead to space-filling turbulence in most of the parameters range

The BinaMIcS project: understanding the origin of magnetic fields in massive stars through close binary systems

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.Comment: To appear in New Windows on Massive Stars, proceedings of the IAU Symposium 30

Rotating double-diffusive convection in stably stratified planetary cores

In planetary fluid cores, the density depends on temperature and chemical composition, which diffuse at very different rates. This leads to various instabilities, bearing the name of double-diffusive convection (DDC). We investigate rotating DDC (RDDC) in fluid spheres. We use the Boussinesq approximation with homogeneous internal thermal and compositional source terms. We focus on the finger regime, in which the thermal gradient is stabilizing whereas the compositional one is destabilizing. First, we perform a global linear stability analysis in spheres. The critical Rayleigh numbers drastically drop for stably stratified fluids, yielding large-scale convective motions where local analyses predict stability. We evidence the inviscid nature of this large-scale double-diffusive instability, enabling the determination of the marginal stability curve at realistic planetary regimes. In particular, we show that in stably stratified spheres, the Rayleigh numbers Ra at the onset evolve like Ra ∼ Ek−1, where Ek is the Ekman number. This differs from rotating convection in unstably stratified spheres, for which Ra ∼ Ek−4/3. The domain of existence of inviscid convection thus increases as Ek−1/3. Secondly, we perform non-linear simulations. We find a transition between two regimes of RDDC, controlled by the strength of the stratification. Furthermore, far from the RDDC onset, we find a dominating equatorially antisymmetric, large-scale zonal flow slightly above the associated linear onset. Unexpectedly, a purely linear mechanism can explain this phenomenon, even far from the instability onset, yielding a symmetry breaking of the non-linear flow at saturation. For even stronger stable stratification, the flow becomes mainly equatorially symmetric and intense zonal jets develop. Finally, we apply our results to the early Earth core. Double diffusion can reduce the critical Rayleigh number by four decades for realistic core conditions. We suggest that the early Earth core was prone to turbulent RDDC, with large-scale zonal flows

Fossil field decay due to nonlinear tides in massive binaries

Context. Surface magnetic fields have been detected in 5–10% of isolated massive stars, hosting outer radiative envelopes. They are often thought to have a fossil origin, resulting from the stellar formation phase. Yet, magnetic massive stars are scarcer in (close) short-period binaries, as reported by the BinaMIcS (Binarity and Magnetic Interaction in various classes of Stars) Collaboration. Aims. Different physical conditions in the molecular clouds giving birth to isolated stars and binaries are commonly invoked. In addition, we propose that the observed lower magnetic incidence in close binaries may be due to nonlinear tides. Indeed, close binaries are probably prone to tidal instability, a fluid instability growing upon the equilibrium tidal flow via nonlinear effects. Yet, stratified effects have hitherto been largely overlooked. Methods. We theoretically and numerically investigate tidal instability in rapidly rotating, stably stratified fluids permeated by magnetic fields. We use the short-wavelength stability method to propose a comprehensive (local) theory of tidal instability at the linear onset, discussing damping effects. Then, we propose a mixing-length theory for the mixing generated by tidal instability in the nonlinear regime. We successfully assess our theoretical predictions against proof-of-concept, direct numerical simulations. Finally, we compare our predictions with the observations of short-period, double-lined spectroscopic binary systems. Results. Using new analytical results, cross-validated by a direct integration of the stability equations, we show that tidal instability can be generated by nonlinear couplings of inertia-gravity waves with the equilibrium tidal flow in short-period massive binaries, even against the Joule diffusion. In the nonlinear regime, a fossil magnetic field can be dissipated by the turbulent magnetic diffusion induced by the saturated tidal flows. Conclusions. We predict that the turbulent Joule diffusion of fossil fields would occur in a few million years for several short-period massive binaries. Therefore, turbulent tidal flows could explain the observed dearth of some short-period magnetic binaries