arXiv:1712.03787v1 [physics.flu-dyn] 11 Dec 2017 yia flqi easadslrsse bet.A a As the reaches objects. flow system the solar when and consequence, metals liquid of typical yaoDSaeteeoersrce omdrtl low moderately Pm typically to State-of-the-art Pm, restricted of values therefore flow. are turbulent DNS fully dynamo gi- the the of require resolve scales would accurately small values to resources Pm the computational realistic gantic of nu- at (DNS) the effect simulations On numerical dynamo direct [1–3]. side, fields merical magnetic succeeded dynamo have producing experiments such in three only and needed, al:bcueo h ihpwrnee ossanatur- a sustain Rm to needed above power flow high bulent the numeri- of because and experimentally cally: both challenge a constitutes h ag 10 range the au Rm value ai icst:vle fP fteodro 10 of order the of Pm of values viscosity: matic 4 n Pm and [4] nti us o cln eain,poigtednm ef- dynamo the probing 8–10]. relations, [5, scaling observations for quest with this Pm DNS In of relating values for motivated essential physically is the to results merical oeOmcdffso.I iesols om hshap- this Rm value form, critical a dimensionless above In pens over- diffusion. Ohmic and the come induction per- when electromagnetic field through magnetic in amplify turbations sets mag- to vigorous instability into sufficiently energy dynamo is flow kinetic The fluid energy. the instability netic of an part effect, converts dynamo the that from originate to lieved ubrR = Rm number h udand fluid the with diffusivity, magnetic nimdaedffiut rssfo h o au fthe = of Pm value low number the Prandtl from magnetic arises difficulty immediate An oiyadlnt clso h o and flow the of scales length and locity eiigadtsigsaiglw htetaoaenu- extrapolate that laws scaling testing and Deriving h antcfil fms srpyia bet sbe- is objects astrophysical most of field magnetic The c nvri´ ai ieo,UiestePer tMreCur Marie et Universit´e Pierre Diderot, Universit´e Paris h iei enlsnme e= Re number Reynolds kinetic the , ≥ 5 rniint h ubln auainrgm cusa lo a at occurs lat regime the t saturation Pm regime, energy turbulent scaling magnetic the turbulent the to low-Pm that transition show a We to DNS. regime 3D scaling fully to reach of out h eutn e frdcdeutosfrP ont 2 to down direction Pm third for equations the c reduced in spatial of two mode t set Fourier on of resulting single vicinity depends the a the then of on field consists bypas focus velocity field To we the Pm, Pm. flow: low number discrepan rotating very Prandtl this at magnetic of saturation the origin dynamo of The values enough magnitude. low of orders t leads several values by parameter realistic numeri to direct scaling-laws viscous 3D state-of-the-art in viscosity-dependent − 0 µ hl h auae antceeg sidpneto visco of independent is energy magnetic saturated the While . Uℓ/η 0 5i peia emty[5–8]. geometry spherical in 05 10 ≃ h antcpreblt fvacuum. of permeability magnetic the 1 6 10 aoaor ePyiu Statistique, Physique de Laboratoire n h o sflytruet This turbulent. fully is flow the and where , − 3 anbrnSeshasayanan Kannabiran c hc xliswyi a enoeloe ynmrclstud numerical by overlooked been has it why explains which , ag xeietlfclte are facilities experimental large , 2 σ nvri´ ai-aly E aly 19 Gif-sur-Yvet 91191 Saclay, CEA Universit´e Paris-Saclay, ≥ evc ePyiu el’ de Physique de Service h lcrclcnutvt of conductivity electrical the U 0 rniint ubln yaosaturation dynamo turbulent to Transition c . 1i rpeproi boxes periodic triple in 01 ftemgei Reynolds magnetic the of and ν/η ℓ r h yia ve- typical the are where , η O 1 = 1 threshold (1) /µ Dtd eebr1,2017) 12, December (Dated: Uℓ/ν ν 0 1 σ skine- is , − ∗ sthe is 5 ttCnese E,CR M 3680, UMR CNRS CEA, Condens´e, Etat ´ sin is aieGallet Basile are cl oml u´rer,CR M 8550, UMR Sup´erieure, CNRS Normale Ecole ´ h icu n.Temgei nryaoethreshold above [11]: scaling-law energy “viscous” magnetic the The follows then one. and viscous force Lorentz the the between dy- balance flows a dynamo high-Pm through laminar saturate Indeed, of examples theoretical number. and value Prandtl namos the magnetic on the depends of crucially that prefactor mensional n xet h antceeg osaelnal with linearly (Rm scale threshold to from onset, energy departure to magnetic the close the planets: expects some laboratory one dynamo possibly all bifurca- to and relevant dynamo experiments is which supercritical regime a a of tion, vicinity for the Consider instance threshold. above the produced affects energy strongly magnetic flow because background desirable, the highly of is turbulence values the Pm lower much at fect 2 hntemgei rnt ubri o enough, low is number Prandtl magnetic Pm scaling-law the the for when evidence computa- provide (2) reasonable they cap- at [13]: to cost properties rather tional but statistical tur- flows, their meant MHD the ture not of describe are models quantitatively models shell to phenomenological on Such based one bulence. the is dictions hntevsossaiglw()b atrPm factor energy a by magnetic (1) higher law scaling much viscous a the than to corresponds which airSoe qain[2.Ti ed ote“turbu- [11]: the energy the magnetic to the of leads for term This scaling-law lent” the advective [12]. nonlinear between equation the balance Navier-Stokes and a this force through Lorentz that achieved indicate is experiments saturation laboratory contrast, By nitrsigapoc ots hs hoeia pre- theoretical these test to approach interesting An e 4reLood 50 ai,Fac,and France, Paris, 75005 Lhomond, rue 24 ie, nudrsiaino h antcenergy magnetic the of underestimation an o . a iuain DS.Etaoaigsuch Extrapolating (DNS). simulations cal · 2 10 au ftemgei rnt number, Prandtl magnetic the of value w n lxnrsAlexakis Alexandros and , [14]. 1 e en needn fvsoiy The viscosity. of independent being ter − epromnmrclsmltosof simulations numerical perform We . iyi yaoeprmns tremains it experiments, dynamo in sity yi htfly3 N antreach cannot DNS 3D fully that is cy 5 hsprmtrrgm scurrently is regime parameter This . odntsol,wietemagnetic the while only, oordinates ednm hehl narapidly a in threshold dynamo he astosfo ihP viscous high-Pm a from ransitions hslmtto n investigate and limitation this s | ρµ B | ρµ | B 0 2 η ℓ | 0 2 2 2 e France. te, η e ofar. so ies ℓ 2 ∝ 2 ∝ m(Rm Pm (Rm − − 1 Rm Rm − c ) Rm c , ) . c ,wt di- a with ), − 1 . (2) (1) 2

The natural next step towards quantitative numeri- where c.c. denotes the complex conjugate and k the verti- cal dynamo models would be to reproduce the turbu- cal wavenumber. Inserting such a Fourier mode into the lent scaling regime directly from numerical solutions of induction equation leads to the evolution equation for b: the Navier-Stokes and induction equations, a task which b ∇ e u b ∇2 2 b remains beyond reach of state-of-the-art fully 3D DNS. ∂t = ( ⊥ + ik z) × ( × )+ η( ⊥ − k ) , (4) Indeed, it has been recently shown that all spherical dy- where ∇⊥ = (∂x, ∂y, 0) and ez is the unit vector along z namo simulations obey the viscous scaling law [8, 10]: be- [17–19]. cause of the moderately low Pm, they are restricted to the We focus on the weakly nonlinear regime in the vicinity viscosity-dominated regime, which makes their extrapo- of the dynamo threshold, for which we can keep only lation to Earth like parameters questionable. A central the first unstable magnetic mode of the form (3). The question is therefore: how much smaller need the mag- instability saturates through the action of the Lorentz netic be to start observing the turbulent force. The latter being quadratic in the magnetic field scaling regime (2)? (3), it contains two harmonics in z: We address this question by focusing on rapidly rotat- ing flows in the vicinity of the dynamo threshold. Global • A z-independent component (harmonic 0), through rotation is a relevant ingredient of many planetary dy- which the magnetic field directly feeds back onto namos, which strongly affects the dynamo characteris- the z-invariant base flow. tics [15]. For rapid global rotation, we are able to reduce • A second harmonic component, with vertical de- the full MHD system to a set of quasi-2D equations gov- ±2 pendence e ikz. It forces a second harmonic in erning the interaction between a two-dimensional three- 2ikz the velocity field, u2 = uˆ2(x,y,t)e + c.c., the component (2D3C) flow and a vertically-dependent dy- amplitude of which results from a balance between namo magnetic field. This approach allows us to by- the Coriolis and Lorentz forces: pass the current limitations of 3D DNS: we focus on the 2 high-Reynolds-number regime where nonlinear advection B u2 ∼ . (5) strongly dominates over viscous effects, in a quasi-2D sys- ρµ0ℓΩ tem of manageable computational cost. The vertical average of the Coriolis force being absorbed by the pressure gradient, the z-independent flow follows 0.5 the equation: 2 0.4 ∂tu + (u · ∇⊥)u = −∇⊥p − γu⊥ + ν∇⊥u + f(x, y) (6) 1 ∇ e b b∗ 0.3 + [[( ⊥ + ik z) × ] × + c.c.] ,

c ρµ0

Rm u 0.2 where ⊥ = (u, v, 0), γ is a linear Ekman friction coef- ficient [20], and the last term is the vertical average of 0.1 the Lorentz force. The latter induces a correction to the 2 turbulent 2D3C base flow of order B /(ρµ0U). By con- 0.0 u 10−5 10−4 10−3 10−2 10−1 100 101 102 103 trast, equation (5) indicates that 2 is smaller than this P m correction by a factor equal to the U/ℓΩ, which is asymptotically small in the present study: we FIG. 1. Threshold magnetic Rmc for dy- 2 −3 therefore discard u2 in the following. The weakly nonlin- namo action as a function of Pm, for γℓ /η = 5.1 10 . ear regime in the vicinity of the dynamo threshold then corresponds to the interaction between the first unstable Reduced equations. We consider a flow driven by a ver- vertical Fourier mode of the magnetic field and the z- f tically invariant body force (x, y) in a frame rotating at invariant 2D3C flow. Their coupled evolution obeys the a rate Ω around the vertical z axis. It was recently proven reduced system of equations (4-6). that the corresponding 2D turbulent flow is stable to 3D Numerical experiments. We solve equations (4-6) in- perturbations provided the Rossby number is sufficiently side a domain (x, y) ∈ [0, 2πL]2 using standard pseudo- low [16]. We focus on this parameter range, the turbulent spectral methods [18, 21]. The body-force has the u 2D3C velocity field = [u(x,y,t), v(x,y,t), w(x,y,t)] be- Roberts-flow geometry: ing the base flow of the present dynamo study. For the kinematic dynamo problem, the invariance of the flow f(x, y)= F [cos(y/ℓ), sin(x/ℓ), cos(x/ℓ)+sin(y/ℓ)] , (7) along the vertical direction allows us to decompose the magnetic field into vertical Fourier modes and to study where the scale of the forcing is ℓ = L/4. This value their evolution independently: was shown to be close to optimal for reducing the dy- namo threshold [22], which is desirable to reach the low- B(x,y,z,t)= b(x,y,t)eikz + c.c. , (3) magnetic-Prandtl-number regime. 3

on F,ν,γ. This flow is then used as the starting point of the MHD simulations. We solve the MHD problem for increasing values of F , i.e., for increasing values of the Rm while keeping Pm fixed (ν,γ and η are fixed). Above the threshold value Rmc for dynamo action, the dynamo instability sets in and magnetic perturbations grow. As shown in figure 1, Rmc varies little with Pm, with Rmc ∈ [0.30;0.42] over seven decades of Pm. Such a constant value of Rmc at low Pm is rather generic and has been reported for several other forcing geometries [25, 26]. By contrast, the be- havior of Rmc in the transitional region of intermediate Pm strongly depends on the structure of the forcing: for some fully 3D flows with weak scale separation, Rmc dis- plays a two-fold increase at intermediate Pm [24]. The weak variation of Rmc with Pm in figure 1 is therefore attributed to both the scale separation and the 2D3C nature of the flow [22, 27]. After the initial growth phase, the magnetic field sat- urates through the feedback of the Lorentz force onto the 2D3C flow. In figure 2, we show snapshots of the velocity and magnetic fields in the saturated state, for Pm = 4.25 · 10−5. In agreement with standard α2- [23], the vertical velocity and vorticity are concentrated at the forcing scale ℓ while the dynamo magnetic field is at large scale, with horizontal compo- nents more energetic than the vertical one. We compute the space and time averaged magnetic energy in the sat- urated state, to produce bifurcation curves such as the ones shown in figure 3. We repeat this procedure for various values of the magnetic Prandtl number Pm, and from each bifurcation curve we extract the slope S(Pm) relating the magnetic energy to the departure from onset:

|B|2 ℓ2 2 = S(Pm) × (Rm − Rmc) , (8) ρµ0η where the angular brackets denote 3D space and time − FIG. 2. Snapshots of the saturated state for Pm = 4.25·10 5 , average. S(Pm) is the central quantity of the present 2 −3 γℓ /η = 5.1 · 10 and Rm = 0.44. Top: vertical vorticity in study. units of U/ℓ. Middle: vertical velocity in units of U. Bot- From viscous to turbulent saturation. In figure 4, we 2 tom: dimensionless magnetic field B ℓ/pρµ0η in the plane show S(Pm) over seven decades of Pm. For large Pm, the z = 0. The arrows indicate the horizontal components while flow is laminar and has a low Reynolds number near the color codes for the vertical one. These arrows correspond to − a typical magnitude 4 · 10 2 of the dimensionless horizontal dynamo threshold. Accordingly, the dynamo saturation magnetic field. Their average direction rotates with z. obeys the viscous scaling-law (1), i.e., S(Pm) ∼ Pm for large Pm. The precise expression of S(Pm) can be deter- mined analytically in the limit of scale separation kℓ ≪ 1 and corresponds to the usual viscous alpha quenching. The vertical wavenumber is set to k = 0.2/L, which 2 corresponds to the lowest wavenumber inside a domain of For small friction γℓ ≪ ν, following the standard weakly height 10πL. The influence of this parameter on the dy- nonlinear approach [28, 29] we obtain: namo threshold and magnetic field geometry was studied 2 in detail in a previous publication [18], the phenomenol- S(Pm) = Pm . (9) 2 r ogy being that of standard α dynamos [23, 24]. kℓ We first solve the purely hydrodynamic problem until This analytical prediction is displayed in figure 4 with a statistically steady state is reached. We denote as U a solid line and captures perfectly the high-Pm behavior the root-mean-square velocity of this state that depends of S(Pm). 4

−3 5 ×10 independent of both ν and γ (see inset of figure 4) and the dominant balance in the Navier-Stokes equation (6) is between the Lorentz force and the nonlinear term. To 4 our knowledge, this study constitutes the first numerical observation of the turbulent scaling regime of dynamo

2 3 saturation in a solution of the coupled Navier-Stokes and ℓ 2 i η

2 induction equations. This is because extremely low val- 0 | B ρµ | ues of Pm are needed to observe such turbulent satu- h 2 ration: the plateau in figure 4 arises for Pm . 10−3, an order of magnitude below the smallest values of Pm 1 achieved in state-of-the-art fully 3D DNS. Discussion. We have reported the numerical observa- 0 tion of the turbulent scaling regime for dynamo satura- 0.00 0.01 0.02 tion in a solution of the MHD equations. Our work there- (Rm − Rm ) c fore quantifies for the first time previous theoretical esti- mates based on dimensional analysis [11] and shell models FIG. 3. Magnetic energy as a function of the departure [14]: in the present setup, the turbulent scaling regime from onset, for several values of the magnetic Prandtl num- 2 −3 sets in for values of the magnetic Prandtl number an or- ber. The friction is γℓ /η = 5.1 10 and symbols are: •, − − − Pm = 1.4 10 3; ▽, Pm = 7.0 10 3; ⋆, Pm = 4.0 10 2; ⋄, der of magnitude lower than currently achieved by state- Pm = 8.9 10−2. The dashed linear fits allow to extract the of-the-art fully 3D DNS. This explains the mismatch be- slope S(Pm) of each bifurcation curve. tween spherical dynamo simulations, which obey the vis- cous scaling law [10], and experimental dynamos, which 4 10 follow the turbulent one [12]. Because the turbulent scal- −1 10 ing regime is likely to be the generic situation for many Pm = 1.70 10−4 3 natural and experimental dynamos, this study paves the 10 way for quantitative numerical estimates of the magnetic field in astrophysical objects and laboratory experiments. 2 10 This work was granted access to the HPC resources of S ∼ Pm1: viscous GENCI-CINES (Project No.x2014056421, x2015056421) 1 −2 10 10 saturation −3 2 −2 and MesoPSL financed by the Region Ile de France and 10 γ ℓ /η 10 S the project EquipMeso ANR-10-EQPX-29-01. KS ac- 0 10 knowledges support from LabEx ENS-ICFP: ANR-10- LABX-0010/ANR-10-IDEX-0001-02 PSL and from Ecole´ −1 10 Doctorale ˆIle de France (EDPIF). BG acknowledges sup- port from LabEx PALM ANR-10-LABX-0039. −2 10 S ∼ Pm0: turbulent saturation

−3 10 −4 −2 0 2 10 10 10 10 ∗ Pm [email protected] [1] R. Stieglitz and U. M¨uller, Physics of Fluids (1994- FIG. 4. Slope S of the magnetic energy above onset. High- present) 13, 561 (2001). Pm solutions are time-independent (blue circles) and obey [2] A. Gailitis, O. Lielausis, E. Platacis, S. Dement’ev, the quantitative prediction (9) from viscous alpha quenching A. Cifersons, G. Gerbeth, T. Gundrum, F. Stefani, (thick solid line). For low Pm, the flow is time-dependent (red M. Christen, and G. Will, Physical Review Letters 86, − triangles). For Pm . 10 3, the dynamo saturation follows 3024 (2001). the turbulent scaling-law (2), represented as a dashed eye- [3] R. Monchaux, M. Berhanu, M. Bourgoin, M. Moulin, guide. The main figure corresponds to a friction coefficient P. Odier, J.-F. Pinton, R. Volk, S. Fauve, N. Mordant, γℓ2/η = 5.1 10−3. The inset highlights the independence of S F. P´etr´elis, et al., Physical review letters 98, 044502 on friction in the turbulent saturation regime. (2007). [4] P. Mininni, Physical Review E 76, 026316 (2007). [5] U. Christensen and J. Aubert, Geophysical Journal In- In contrast with such viscous dynamos, S(Pm) reaches ternational 166, 97 (2006). [6] U. Christensen, Space science reviews 152, 565 (2010). a plateau at low Pm, with values orders of magnitude [7] C. A. Jones, Annual Review of Fluid Mechanics 43, 583 larger than predicted by the laminar theory. This corre- (2011). sponds to the turbulent scaling-law (2), for which S(Pm) [8] E. King and B. Buffett, Earth and Planetary Science Let- is independent of Pm. More precisely, in this regime S is ters 371, 156 (2013). 5

[9] P. Davidson, Geophysical Journal International 195, 67 [20] J. Pedlosky, Geophysical fluid dynamics (Springer, 1979). (2013). [21] D. O. Gomez, P. D. Mininni, and P. Dmitruk, Phys. Scr. [10] L. Oruba and E. Dormy, Geophysical Journal Interna- 116, 123 (2005). tional 198, 828 (2014). [22] M. Sadek, A. Alexakis, and S. Fauve, [11] F. P´etr´elis and S. Fauve, The European Physical Jour- Phys. Rev. Lett. 116, 074501 (2016). nal B-Condensed Matter and Complex Systems 22, 273 [23] H. Moffatt, Magnetic field generation in electrically con- (2001). ducting fluids (Cambridge University Press, 1978). [12] F. P´etr´elis, N. Mordant, and S. Fauve, Geophysical and [24] Y. Ponty and F. Plunian, Physical review letters 106 Astrophysical 101, 289 (2007). (2011). [13] F. Plunian, R. Stepanov, and P. Frick, Physics Reports [25] Y. Ponty, P. Mininni, D. Montgomery, J.-F. Pinton, 523 (2013). H. Politano, and A. Pouquet, Physical review letters [14] R. Stepanov and F. Plunian, Journal of Turbulence 7 94 (2005). (2006). [26] A. Iskakov, A. Schekochihin, S. Cowley, J. McWilliams, [15] J.-F. Donati, in Stellar Rotation, Vol. 215 (2004) p. 258. and P. M.R.E., Physical review letters 98 (2007). [16] B. Gallet, Journal of Fluid Mechanics 783, 412 (2015). [27] K. Seshasayanan, V. Dallas, and A. Alexakis, arXiv [17] S. G. L. Smith and S. Tobias, Journal of Fluid Mechanics preprint arXiv:1701.08996 (2017). 498, 1 (2004). [28] A. D. Gilbert and P.-L. Sulem, Geophysical & Astrophys- [18] K. Seshasayanan and A. Alexakis, ical Fluid Dynamics 51, 243 (1990). J. Fluid Mech. 799, 246 (2016). [29] S. Fauve and F. P´etr´elis, Peyresq lectures on nonlinear [19] K. Seshasayanan and A. Alexakis, Journal of Fluid Me- phenomena (World Scientific, 2003). chanics 806, 627648 (2016).