A unified approach to a class of slow dynamos

Abstract

Dynamo action in a highly conducting fluid with small magnetic diffusivity η is particularly sensitive to the topology of the flow. The sites of rapid magnetic field regeneration, when they occur, appear to be located at the stagnation points or in regions where the particle paths are chaotic. Elsewhere only slow dynamo action is to be expected. Two such examples are the nearly axially symmetric dynamo of Braginsky and the generalisation to smooth velocity fields of the Ponomarenko dynamo. Here a method of solution is developed, which applies to both these examples and is applicable to other situations, where magnetic field lines are close to either closed or spatially periodic contours. Particular attention is given to field generation in the neighbourhood of resonant surfaces where growth rates may be intermediate between the slow diffusive and fast convective time scales. The method is applied to the case of the two-dimensional ABC-flows, where it is shown that such intermediate dynamo action can occur on resonant surfaces.     

The motion of magnetic fronts in spiral galaxies

Abstract

We study the nonlinear asymptotic thin disc approximation to the mean field dynamo equations, as applicable to spiral galaxies. The circumstances in which sharp magnetic field structures (fronts) can propagate radially are investigated, and an expression for the speed of propagation derived. We find that the speed of an interior front is proportional to η//R _? (where η is the diffusivity and R_t the galactic radius), whereas an exterior front moves with speed of order , where γ is the local growth rate of the dynamo. Numerical simulations are presented, that agree well with our asymptotic results. Further, we perform numerical experiments using the 'no-z' approximation for thin disc dynamos, and show that the propagation of magnetic fronts in this approximation can also be understood in terms of our asymptotic results.     

Kinematic stationary geodynamo models with separated toroidal and poloidal motions

Abstract

This paper is concerned with a three-dimensional spherical model of a stationary dynamo that consists of a convective layer with a simple poloidal flow of the S^2c ₂ kind between a rotating inner body core and solid outer shell. The rotation of the inner core and the outer shell means that there are regions of concentrated shear or differential rotation at the convective layer boundaries. The induction equation for the inside of the convective layer was solved numerically by the Bullard-Gellman method, the eigenvalue of the problem being the magnetic Reynolds number of the poloidal flow (R _M2) and it was assumed that the magnetic Reynolds number of the core (R _M1) and of the shell (R _M3) were prescribed parameters. Hence R _M2 was studied as a function of R _M1 and R _M3, along with the orientation of the rotation axis, the radial dependence of the poloidal velocity and the relative thickness of the layers for the three different situations, (i) the core alone rotating, (ii) the shell alone rotating and (iii) the core and the shell rotating together. In all three cases it was found that, at definite orientations of the rotation axis, there is a good convergence of both the eigenvalues and the eigenfunctions of the problem as the number of spherical harmonics used to represent the problem increases. For R _M1 =R _M3= 10³, corresponding to the westward drift velocity and the parameters of the Earth's core, the critical values of R _M2 are found to be three orders of magnitude lower than R _M1, R _M3 so that the poloidal flow velocity sufficient for maintaining the dynamo process is 10-20 m/yr. With only the core or the shell rotating, the velocity field generally differs little from the axially symmetric case. However, for R _M2 (or R _M3) lying in the range 10² to 10⁵, the self-excitation condition is found to be of the form R _M2˙R ^½ _M1=constant (or R _M2˙R^½ _M3=constant) and the solution does not possess the properties of the Braginsky near-axisymmetric dynamo. We should expect this, in particular, in the Braginsky limit R _M2˙R^?½; _M1=constant.

An analysis of known three-dimensional dynamo models indicates the importance of the absence of mirror symmetry planes for the efficient generation of magnetic fields.     

A nonlinear dynamo driven by rapidly rotating convection

In Kim et al. (Kim, E., Hughes, D.W. and Soward, A.M., “An investigation into high conductivity dynamo action driven by rotating convection”, Geophys. Astrophys. Fluid Dynam. 91, 303–332 ().) we investigated kinematic dynamo action driven by rapidly rotating convection in a cylindrical annulus. Here we extend this work to consider self-consistent nonlinear dynamo action in which the back-reaction of the Lorentz force on the flow is taken into account. In particular, we investigate, as a function of magnetic Prandtl number, the evolution of an initially weak magnetic field in two different types of convective flow – one chaotic and the other integrable. On saturation, the latter shows a systematic dependence on the magnetic Prandtl number whereas the former appears not to. In addition, we show how, in keeping with the findings of Cattaneo et al. (Cattaneo, F., Hughes, D.W. and Kim, E., “Suppression of chaos in a simplified nonlinear dynamo model”, Phys. Rev. Lett. 76, 2057–2060 ().), saturation of the growth of the magnetic field is brought about, for the originally chaotic flow, by a strong suppression of chaos.     

On a Physically Realistic Fast Dynamo

As a step towards a physically realistic model of a fast dynamo, we study numerically a kinematic dynamo driven by convection in a rapidly rotating cylindrical annulus. Convection maintains the quasi-geostrophic balance whilst developing more complicated time-dependence as the Rayleigh number is increased. We incorporate the effects of Ekman suction and investigate dynamo action resulting from a chaotic flow obtained in this manner. We examine the growth rate as a function of magnetic Prandtl number Pm, which is proportional to the magnetic Reynolds number. Even for the largest value of Pm considered, a clearly identifiable asymptotic behaviour is not established. Nevertheless the available evidence strongly suggests a fast dynamo process.     

Numerical simulations of transport of non-ergodic solute plumes in heterogeneous aquifers

Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S _{i i} (t)>, and the plume centroid covariance, R _{i i} (t) (i=1,2), were simulated for the variance of Y=log T, σ _Y ²=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ _Y ²=0.1, all simulated <S _{i i} (t)>−S _{i i} (0) and R _{i i} (t) agree well with the first-order theoretical values, where S _{i i} (0) are the initial values of S _{i i} (t). For σ _Y ²=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S ₁₁(t)>−S ₁₁(0) and R ₁₁(t), agree well with the first-order theoretical results but the simulated transverse moments <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) are significantly larger than the first-order values. For the same two larger values of σ _Y ² but the line sources parallel to μ, the simulated <S ₁₁(t)>−S ₁₁(0) are larger than but the simulated R ₁₁ are smaller than the first-order values, and both simulated <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) stay larger than the first-order values. For a fixed value of σ _Y ², the summations of <S _{i i} (t)>−S _{i i} (0) and R _{i i} , i.e., X _{i i} (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X ₁₁ are in good agreement with the first-order theory but the simulated X ₂₂ are significantly larger than the first-order values. The simulated X ₂₂, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S ₂₂(t)>−S ₂₂(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ _Y ² is small.     

Numerical simulations of transport of non-ergodic solute plumes in heterogeneous aquifers

Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S _{i i} (t)>, and the plume centroid covariance, R _{i i} (t) (i=1,2), were simulated for the variance of Y=log T, σ _Y ²=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ _Y ²=0.1, all simulated <S _{i i} (t)>−S _{i i} (0) and R _{i i} (t) agree well with the first-order theoretical values, where S _{i i} (0) are the initial values of S _{i i} (t). For σ _Y ²=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S ₁₁(t)>−S ₁₁(0) and R ₁₁(t), agree well with the first-order theoretical results but the simulated transverse moments <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) are significantly larger than the first-order values. For the same two larger values of σ _Y ² but the line sources parallel to μ, the simulated <S ₁₁(t)>−S ₁₁(0) are larger than but the simulated R ₁₁ are smaller than the first-order values, and both simulated <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) stay larger than the first-order values. For a fixed value of σ _Y ², the summations of <S _{i i} (t)>−S _{i i} (0) and R _{i i} , i.e., X _{i i} (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X ₁₁ are in good agreement with the first-order theory but the simulated X ₂₂ are significantly larger than the first-order values. The simulated X ₂₂, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S ₂₂(t)>−S ₂₂(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ _Y ² is small.     

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.     

Stellar Dynamo Waves: Asymptotic Configurations

We investigate the parameter space of a Parker dynamo with a simple alpha quenching nonlinearity, taking as governing parameters the dynamo number D (D<0) and the ratio of diffusion times in the radial and latitudinal directions in the convective zone. The latter parameter, μ, is connected with the aspect ratio (dimensionless thickness) of the convective zone. We isolate two asymptotic configuration of the dynamo waves excited by the Parker dynamo in the limiting case of strong generation. Apart from the standard case with the solar type dynamo wave travelling from mid-latitudes to the equator, we describe a form of dynamo activity which is basically an anharmonic standing wave. The first situation occurs when μ increases with |D|. With μ fixed and |D| increasing, the second asymptotic configuration occurs. We discuss possibilities of identifying these asymptotic configurations with various types of stellar activity as traced by stellar CaII data.     

Analysis of the nonlinear storage–discharge relation for hillslopes through 2D numerical modelling

Storage–discharge curves are widely used in several hydrological applications concerning flow and solute transport in small catchments. This article analyzes the relation Q(S) (where Q is the discharge and S is the saturated storage in the hillslope), as a function of some simple structural parameters. The relation Q(S) is evaluated through two‐dimensional numerical simulations and makes use of dimensionless quantities. The method lies in between simple analytical approaches, like those based on the Boussinesq formulation, and more complex distributed models. After the numerical solution of the dimensionless Richards equation, simple analytical relations for Q(S) are determined in dimensionless form, as a function of a few relevant physical parameters. It was found that the storage–discharge curve can be well approximated by a power law function Q/(LK_s) = a(S/(L²(? ? θ_r)))^b, where L is the length of the hillslope, K_s the saturated conductivity, ? ? θ_r the effective porosity, and a, b two coefficients which mainly depend on the slope. The results confirm the validity of the widely used power law assumption for Q(S). Similar relations can be obtained by performing a standard recession curve analysis. Although simplified, the results obtained in the present work may serve as a preliminary tool for assessing the storage–discharge relation in hillslopes. Copyright © 2012 John Wiley & Sons, Ltd.     

Fast dynamo action in the Ponomarenko dynamo

Abstract

An analysis of small-scale magnetic fields shows that the Ponomarenko dynamo is a fast dynamo; the maximum growth rate remains of order unity in the limit of large magnetic Reynolds number. Magnetic fields are regenerated by a “stretch-diffuse” mechanism. General smooth axisymmetric velocity fields are also analysed; these give slow dynamo action by the same mechanism.     

Nonlinear dynamics of a stellar dynamo in a spherical shell

Abstract

A simple mean-field model of a nonlinear stellar dynamo is considered, in which dynamo action is supposed to occur in a spherical shell, and where the only nonlinearity retained is the influence of the Lorentz forces on the zonal flow field. The equations are simplified by truncating in the radial direction, while full latitudinal dependence is retained. The resulting nonlinear p.d.e.'s in latitude and time are solved numerically, and it is found that while regular dynamo wave type solutions are stable when the dynamo number D is sufficiently close to its critical value, there is a wide variety of stable solutions at larger values of D. Furthermore, two different types of dynamo can coexist at the same parameter values. Implications for fields in late-type stars are discussed.     

An antidynamo theorem for partly symmetric flows

Abstract

It is shown that magnetic fields generated by flows v _r,(r,t)e_r+v_T where v_T is an arbitrary toroidal component (e_r˙v_T≡V≡v_T≡0), cannot be maintained indefinitely against ohmic dissipation. The poloidal field variable max |r ² B _r| is shown to decay strictly monotonically with an undetermined decay rate. A bound on the growth of the toroidal field norm ∥T∥₁ is established solely dependent on the rate of conversion of poloidal to toroidal field, so that when the poloidal field is negligible then ∥T∥₁ decays strictly monotonically. The main application of these results is to models of stellar evolution based on axisymmetric differential rotation and spherically symmetric contraction. This symmetric velocity theorem overlaps with two already known theorems, namely the toroidal velocity theorem where v _r≡0 and the radial velocity theorem where v_T≡0. The new theorem does not entirely include the already established ones, principal differences being in the rates of decay and the field variables for which the decay is proven (see Table 1).     

The Displacement and Strain Field of Three-dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains ε_xx (r, t), ε_yy (r, t) and ε_zz (r, t) and the bulk strain θ (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (α, β, γ) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.     

Space-Time scales of internal waves

Abstract

We have contrived a model E(αω) α μ^?1ω^?p+1(ω ²?ω _i ²)^?+ for the distribution of internal wave energy in horizontal wavenumber, frequency-space, with wavenumber α extending to some upper limit μ(ω) α ω ^r-1 (ω ²?ω _i ²)^½, and frequency ω extending from the inertial frequency ω _i to the local Väisälä frequency n(y). The spectrum is portrayed as an equivalent continuum to which the modal structure (if it exists) is not vital. We assume horizontal isotropy, E(α, ω) = 2παE(α₁, α₂, ω), with α₁, α₂ designating components of α. Certain moments of E(α₁, α₂, ω) can be derived from observations. (i) Moored (or freely floating) devices measuring horizontal current u(t), vertical displacement η(t),…, yield the frequency spectra F _(u,η,…)(ω) = ∫∫ (U ², Z ²,…)E(α₁, ∞₂, ω) dα₁ dα₂, where U, Z,… are the appropriate wave functions. (ii) Similarly towed measurements give the wavenumber spectrum F _(…)(α1) = ∫∫… dα₂ dω. (iii) Moored measurements horizontally separated by X yield the coherence spectrum R(X, ω) which is related to the horizontal cosine transform ∫∫ E(α1, α2 ω) cos α₁ Xdα₁ dα₁. (iv) Moored measurements vertically separated by Y yield R(Y, ω) and (v) towed measurements vertically separated yield R(Y, α₁), and these are related to similar vertical Fourier transforms. Away from inertial frequencies, our model E(α, ω) α ω ^?p-r for α ≦ μ ω ω ^r, yields F(ω) ∞ ω ^?p, F(α₁) ∞ α₁ ^?q, with q = (p + r ? 1)/r. The observed moored and towed spectra suggest p and q between 5/3 and 2, yielding r between 2/3 and 3/2, inconsistent with a value of r = 2 derived from Webster's measurements of moored vertical coherence. We ascribe Webster's result to the oceanic fine-structure. Our choice (p, q, r) = (2, 2, 1) is then not inconsistent with existing evidence. The spectrum is E(∞, ω) ∞ ω ^?1(ω ²?ω _i ^{2 ?1}, and the α-bandwith μ ∞ (ω ²?ω _i ²)⁺ is equivalent to about 20 modes. Finally, we consider the frequency-of-encounter spectra F([sgrave]) at any towing speed S, approaching F(ω) as S ≦ S _o, and F(α₁) for α₁ = [sgrave]/S as S ≧ S _o, where S _o = 0(1 km/h) is the relevant Doppler velocity scale.     

Magnetic field generation by the motion of a highly conducting fluid

Abstract

Using an asymptotic expansion of Green's function for the problem of magnetic field generation by 3D steady flow of highly conducting fluid a general antidynamo theorem is proved in the case of no exponential stretching of liquid particles. Explicit formulae connecting the spectrum of the magnetic modes with the geometry of the Lagrangian trajectories are obtained. The existence of the fast dynamo action for special flows with exponential stretching is proved under the condition of smoothness of the fields of stretching and non-stretching directions.     

Oscillatory large-scale dynamos from Cartesian convection simulations

We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo is expected to be of α² Ω type, which is generally expected to be relevant for the Sun, whereas the latter case refers to α² dynamos that are more likely to occur in more rapidly rotating stars whose differential rotation is small. We perform a parameter study where the shear flow and the rotational influence are varied to probe the relative importance of both types of dynamos. Oscillatory solutions are preferred both in the kinematic and saturated regimes when the negative ratio of shear to rotation rates, q?≡??S/Ω, is between 1.5 and 2, i.e. when shear and rotation are of comparable strengths. Other regions of oscillatory solutions are found with small values of q, i.e. when shear is weak in comparison to rotation, and in the regime of large negative qs, when shear is very strong in comparison to rotation. However, exceptions to these rules also appear so that for a given ratio of shear to rotation, solutions are non-oscillatory for small and large shear, but oscillatory in the intermediate range. Changing the boundary conditions from vertical field to perfect conductor ones changes the dynamo mode from oscillatory to quasi-steady. Furthermore, in many cases an oscillatory solution exists only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. However, the cases with rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic boundary conditions. The strengths of total and large-scale components of the magnetic field in the saturated state, however, are sensitive to the chosen boundary conditions.     

The α-effect in 2D fast dynamos

We look at the large-scale dynamo properties of spatially periodic, time dependent, helical 2D flows of the form u(x, t)?=?(? _y ?ψ?(x, y, t), ?? _x ?ψ?(x, y, t), ?ψ (x, y, t). These flows act as kinematic fast dynamos and are able to generate a mean magnetic field uniform and constant in the xy-plane but whose direction varies periodically along z with wavenumber k. Using Mean Field Electrodynamics, the generation mechanism can be understood in terms of a k-dependent α-effect, which depends on the magnetic Reynolds number, R _m . We calculate this effect for different motions and investigate how its limit as k?→?0 depends on R _m and on the properties of the flows such as their spatial structure or correlation time. This work generalises earlier studies based on 2D steady flows to motions with time dependence.