BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Emmanuel Dormy

Emmanuel Dormy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.geo-phastro-ph.EPphysics.comp-phastro-ph.SRmath.APmath-phmath.MPmath.NAnlin.CDnlin.PSNumerical Analysis

Catalog footprint

What is connected

19works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2021arXiv

Weak branch and multimodal convection in rapidly rotating spheres at low Prandtl number

The focus of this study is to investigate primary and secondary bifurcations to weakly nonlinear flows (weak branch) in convective rotating spheres in a regime where only strongly nonlinear oscillatory sub- and super-critical flows (strong branch) were previously found in [E. J. Kaplan, N. Schaeffer, J. Vidal, and P. Cardin, Phys. Rev. Lett. 119, 094501 (2017)]. The relevant regime corresponds to low Prandtl and Ekman numbers, indicating a predominance of Coriolis forces and thermal diffusion in the system. We provide the bifurcation diagrams for rotating waves (RWs) computed by means of continuation methods and the corresponding stability analysis of these periodic flows to detect secondary bifurcations giving rise to quasiperiodic modulated rotating waves (MRWs). Additional direct numerical simulations (DNS) are performed for the analysis of these quasiperiodic flows for which Poincaré sections and kinetic energy spectra are presented. The diffusion time scales are investigated as well. Our study reveals very large initial transients (more than 30 diffusion time units) for the nonlinear saturation of solutions on the weak branch, either RWs or MRWs, when DNS are employed. In addition, we demonstrate that MRWs have multimodal nature involving resonant triads. The modes can be located in the bulk of the fluid or attached to the outer sphere and exhibit multicellular structures. The different resonant modes forming the nonlinear quasiperiodic flows can be predicted with the stability analysis of RWs, close to the Hopf bifurcation point, by analyzing the leading unstable Floquet eigenmode.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The vortex method for 2D ideal flows in exterior domains

The vortex method is a common numerical and theoretical approach used to implement the motion of an ideal flow, in which the vorticity is approximated by a sum of point vortices, so that the Euler equations read as a system of ordinary differential equations. Such a method is well justified in the full plane, thanks to the explicit representation formulas of Biot and Savart. In an exterior domain, we also replace the impermeable boundary by a collection of point vortices generating the circulation around the obstacle. The density of these point vortices is chosen in order that the flow remains tangent at midpoints between adjacent vortices and that the total vorticity around the obstacle is conserved. In this work, we provide a rigorous justification of this method for any smooth exterior domain, one of the main mathematical difficulties being that the Biot--Savart kernel defines a singular integral operator when restricted to a curve (here, the boundary of the domain). We also introduce an alternative method---the fluid charge method---which, as we argue, is better conditioned and therefore leads to significant numerical improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A systematic method to enforce conservativity on semi-Lagrangian schemes

Semi-Lagrangian schemes have proven to be very efficient to model advection problems. However most semi-Lagrangian schemes are not conservative. Here, a systematic method is introduced in order to enforce the conservative property on a semi-Lagrangian advection scheme. This method is shown to generate conservative schemes with the same linear stability range and the same order of accuracy as the initial advection scheme from which they are derived. We used a criterion based on the column-balance property of the schemes to assess their conservativity property. We show that this approach can be used with large CFL numbers and third order schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multi-stage high order semi-Lagrangian schemes for incompressible flows in Cartesian geometries

Efficient transport algorithms are essential to the numerical resolution of incompressible fluid flow problems. Semi-Lagrangian methods are widely used in grid based methods to achieve this aim. The accuracy of the interpolation strategy then determines the properties of the scheme. We introduce a simple multi-stage procedure which can easily be used to increase the order of accuracy of a code based on multi-linear interpolations. This approach is an extension of a corrective algorithm introduced by Dupont \& Liu (2003, 2007). This multi-stage procedure can be easily implemented in existing parallel codes using a domain decomposition strategy, as the communications pattern is identical to that of the multi-linear scheme. We show how a combination of a forward and backward error correction can provide a third-order accurate scheme, thus significantly reducing diffusive effects while retaining a non-dispersive leading error term.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strong-Field Spherical Dynamos

Numerical models of the geodynamo are usually classified in two categories: those denominated dipolar modes, observed when the inertial term is small enough, and multipolar fluctuating dynamos, for stronger forcing. We show that a third dynamo branch corresponding to a dominant force balance between the Coriolis force and the Lorentz force can be produced numerically. This force balance is usually referred to as the strong-field limit. This solution co-exists with the often described viscous branch. Direct numerical simulations exhibit a transition from a weak-field dynamo branch, in which viscous effects set the dominant length scale, and the strong-field branch in which viscous and inertial effects are largely negligible. These results indicate that a distinguished limit needs to be sought to produce numerical models relevant to the geodynamo and that the usual approach of minimizing the magnetic Prandtl number (ratio of the fluid kinematic viscosity to its magnetic diffusivity) at a given Ekman number is misleading.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Equatorial Ekman Layer

The steady incompressible viscous flow in the wide gap between spheres rotating about a common axis at slightly different rates (small Ekman number E) has a long and celebrated history. The problem is relevant to the dynamics of geophysical and planetary core flows, for which, in the case of electrically conducting fluids, the possible operation of a dynamo is of considerable interest. A comprehensive asymptotic study, in the limit E<<1, was undertaken by Stewartson (J. Fluid Mech. 1966, vol. 26, pp. 131-144). The mainstream flow, exterior to the E^{1/2} Ekman layers on the inner/outer boundaries and the shear layer on the inner sphere tangent cylinder C, is geostrophic. Stewartson identified a complicated nested layer structure on C, which comprises relatively thick quasi-geostrophic E^{2/7} (inside C) and E^{1/4} (outside C) layers. They embed a thinner E^{1/3} ageostrophic shear layer (on C), which merges with the inner sphere Ekman layer to form the E^{2/5} Equatorial Ekman layer of axial length E^{1/5}. Under appropriate scaling, this $E^{2/5}$--layer problem may be formulated, correct to leading order, independent of E. Accordingly, the Ekman boundary layer and ageostrophic shear layer become features of the far-field (as identified by the large value of the scaled axial co-ordinate z) solution. We present a numerical solution, which uses a non-local integral boundary condition at finite $z$ to account for the far-field behaviour. Adopting z^{-1} as a small parameter we extend Stewartson's similarity solution for the ageostrophic shear layer to higher orders. This far-field solution agrees well with that obtained from our numerical model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward an asymptotic behaviour of the ABC dynamo

The ABC flow was originally introduced by Arnol'd to investigate Lagrangian chaos. It soon became the prototype example to illustrate magnetic-field amplification via fast dynamo action, i.e. dynamo action exhibiting magnetic-field amplification on a typical timescale independent of the electrical resistivity of the medium. Even though this flow is the most classical example for this important class of dynamos (with application to large-scale astrophysical objects), it was recently pointed out (Bouya Ismaël and Dormy Emmanuel, Phys. Fluids, 25 (2013) 037103) that the fast dynamo nature of this flow was unclear, as the growth rate still depended on the magnetic Reynolds number at the largest values available so far $(\text{Rm} = 25000)$ . Using state-of-the-art high-performance computing, we present high-resolution simulations (up to 40963) and extend the value of $\text{Rm}$ up to $ 5\cdot10^5$ . Interestingly, even at these huge values, the growth rate of the leading eigenmode still depends on the controlling parameter and an asymptotic regime is not reached yet. We show that the maximum growth rate is a decreasing function of $\text{Rm}$ for the largest values of $\text{Rm}$ we could achieve (as anticipated in the above-mentioned paper). Slowly damped oscillations might indicate either a new mode crossing or that the system is approaching the limit of an essential spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Influence of the mass distribution on the magnetic field topology

Context. Three-dimensional spherical dynamo simulations carried out within the framework of the anelastic approximation have revealed that the established distinction between dipolar and multipolar dynamos tends to be less clear than it was in Boussinesq studies. This result was first interpreted as a direct consequence of the existence of a larger number of models with a high equatorial dipole contribution, together with an intermediate dipole field strength. However, this finding has not been clearly related to specific changes that would have been introduced by the use of the anelastic approximation. Aims. In this paper, we primarily focus on the effects of choosing of a different mass distribution. Indeed, it is likely to have as large consequences as would taking a stratified reference state into account, especially when comparing our results to previous Boussinesq studies. Methods. Our investigation is based on a systematic parameter study of weakly stratified anelastic dynamo models. Results. We show that the tendencies highlighted in previous anelastic dynamo simulations are already present in the Boussinesq limit. Thus they cannot be systematically related to anelastic effects. Actually, a central mass distribution can result in changes in the magnetic field topology that are mainly due to the concentration of convective cells close to the inner sphere.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predictive Scaling Laws for Spherical Rotating Dynamos

State of the art numerical models of the Geodynamo are still performed in a parameter regime extremely remote from the values relevant to the physics of the Earth's core. In order to establish a connection between dynamo modeling and the geophysical motivation, {it is necessary to use} scaling laws. Such scaling laws establish the dependence of essential quantities (such as the magnetic field strength) on measured or controlled quantities. They allow for a direct confrontation of advanced models with geophysical {constraints}. (...) We show that previous empirical scaling laws for the magnetic field strength essentially reflect the statistical balance between energy production and dissipation for saturated dynamos. Such power based scaling laws are thus necessarily valid for any dynamo in statistical equilibrium and applicable to any numerical model, irrespectively of the dynamo mechanism. We show that direct numerical fits can provide contradictory results owing to biases in the parameters space covered in the numerics and to the role of a priori hypothesis on the fraction of ohmic dissipation. We introduce predictive scaling laws, i.e. relations involving input parameters of the governing equations only. We guide our reasoning on physical considerations. We show that our predictive scaling laws can properly describe the numerical database and reflect the dominant forces balance at work in these numerical simulations. We highlight the dependence of the magnetic field strength on the rotation rate. Finally, our results stress that available numerical models operate in a viscous dynamical regime, which is not relevant to the Earth's core.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transition between viscous dipolar and inertial multipolar dynamos

We investigate the transition from steady dipolar to reversing multipolar dynamos. The Earth has been argued to lie close to this transition, which could offer a scenario for geomagnetic reversals. We show that the transition between dipolar and multipolar dynamos is characterized by a three terms balance (as opposed to the usually assumed two terms balance), which involves the non-gradient parts of inertial, viscous and Coriolis forces. We introduce from this equilibrium the sole parameter ${\rm Ro}\,{\rm E}^{-1/3} \equiv {\rm Re}\,{\rm E}^{2/3}$, which accurately describes the transition for a wide database of 132 fully three dimensional direct numerical simulations of spherical rotating dynamos (courtesy of U. Christensen). This resolves earlier contradictions in the literature on the relevant two,terms balance at the transition. Considering only a two terms balance between the non-gradient part of the Coriolis force and of inertial forces, provides the classical ${\rm Ro}/{\ell_u}$ (Christensen and Aubert, 2006). This transition can be equivalently described by ${\rm Re} \, {\ell^{2}_u}$, which corresponds to the two terms balance between the non-gradient part of inertial forces and viscous forces (Soderlund {\it et al.}, 2012).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy Fluxes during Dynamo Reversals

Using direct numerical simulations of the equations of magnetohydrodynamics, we study reversals of the magnetic field generated by the flow of an electrically conducting fluid in a sphere. We show that at low magnetic Prandtl numbers, Pm=0.5, the decrease of magnetic energy, ohmic dissipation and power of the Lorentz force during a reversal is followed by an increase of the power injected by the force driving the flow and an increase of viscous dissipation. Cross correlations show that the Lorentz energy flux is in advance with respect to the other energy fluxes. We also observe that during a reversal, the maximum of the magnetic energy density migrates from one hemisphere to the other and comes back to its initial position, in agreement with recent experimental observations. For larger magnetic Prandtl numbers (Pm= 1, 2), the magnetic field reversals do not display these trends and strongly differ one from another.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Intermittency in spherical Couette dynamos

We investigate dynamo action in three-dimensional numerical simulations of turbulent spherical Couette flows. Close to the onset of dynamo action, the magnetic field exhibits an intermittent behavior, characterized by a series of short bursts of the magnetic energy separated by low-energy phases. We show that this behavior corresponds to the so-called on-off intermittency. This behavior is here reported for dynamo action with realistic boundary conditions. We investigate the role of magnetic boundary conditions in this phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bistability between equatorial and axial dipoles during magnetic field reversals

Numerical simulations of the geodynamo in presence of an heterogeneous heating are presented. We study the dynamics and the structure of the magnetic field when the equatorial symmetry of the flow is broken. If the symmetry breaking is sufficiently strong, the m = 0 axial dipolar field is replaced by an hemispherical magnetic field, dominated by an oscillating m = 1 magnetic field. Moreover, for moderate symmetry breaking, a bistability between the axial and the equatorial dipole is observed. In this bistable regime, the axial magnetic field exhibits chaotic switches of its polarity, involving the equatorial dipole during the transition period. This new scenario for magnetic field reversals is discussed within the framework of the Earth's dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dipole Collapse and Dynamo Waves in Global Direct Numerical Simulations

Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipole to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models obtained by direct numerical simulations. The control parameters, the aspect ratio and the mechanical boundary conditions have been varied to build up this sample of models. Both, strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes can in general be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case, the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos co-exist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated 'single mode'. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Revisiting the ABC flow dynamo

The ABC flow is a prototype for fast dynamo action, essential to the origin of magnetic field in large astrophysical objects. Probably the most studied configuration is the classical 1:1:1 flow. We investigate its dynamo properties varying the magnetic Reynolds number Rm. We identify two kinks in the growth rate, which correspond respectively to an eigenvalue crossing and to an eigenvalue coalescence. The dominant eigenvalue becomes purely real for a finite value of the control parameter. Finally we show that even for Rm = 25000, the dominant eigenvalue has not yet reached an asymptotic behaviour. Its still varies very significantly with the controlling parameter. Even at these very large values of Rm the fast dynamo property of this flow cannot yet be established.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The dynamo bifurcation in rotating spherical shells

We investigate the nature of the dynamo bifurcation in a configuration applicable to the Earth's liquid outer core, i.e. in a rotating spherical shell with thermally driven motions. We show that the nature of the bifurcation, which can be either supercritical or subcritical or even take the form of isola (or detached lobes) strongly depends on the parameters. This dependence is described in a range of parameters numerically accessible (which unfortunately remains remote from geophysical application), and we show how the magnetic Prandtl number and the Ekman number control these transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Chaotic magnetic field reversals in turbulent dynamos

We present direct numerical simulations of reversals of the magnetic field generated by swirling flows in a spherical domain. In agreement with a recent model, we observe that coupling dipolar and quadrupolar magnetic modes by an asymmetric forcing of the flow generates field reversals. In addition, we show that this mechanism strongly depends on the value of the magnetic Prandtl number.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetostrophic MRI in the Earth's Outer Core

We show that a simple, modified version of the Magnetorotational Instability (MRI) can develop in the outer liquid core of the Earth, in the presence of a background shear. It requires either thermal wind, or a primary instability, such as convection, to drive a weak differential rotation within the core. The force balance in the Earth's core is very unlike classical astrophysical applications of the MRI (such as gaseous disks around stars). Here, the weak differential rotation in the Earth core yields an instability by its constructive interaction with the planet's much larger rotation rate. The resulting destabilising mechanism is just strong enough to counteract stabilizing resistive effects, and produce growth on geophysically interesting timescales. We give a simple physical explanation of the instability, and show that it relies on a force balance appropriate to the Earth's core, known as magnetostrophic balance.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the ill-posedness of the Prandtl equation

The concern of this paper is the Cauchy problem for the Prandtl equation. This problem is known to be well-posed for analytic data, or for data with monotonicity properties. We prove here that it is linearly ill-posed in Sobolev type spaces. The key of the analysis is the construction, at high tangential frequencies, of unstable quasimodes for the linearization around solutions with non-degenerate critical points. Interestingly, the strong instability is due to vicosity, which is coherent with well-posedness results obtained for the inviscid version of the equation. A numerical study of this instability is also provided.

0 citations0 reviewsNo code linkedOpen access