the kinematic dynamo effect is generally believed to be facilitated by a separation of spatial scales between the flow and the magnetic field .
this idea has its roots in the discovery of the @xmath5-effect ( steenbeck , 1966 ; krause and rdler , 1980 ) for non - parity - invariant flows ( an instance of which are helical flows ) .
indeed , recent results indicate that helicity and scale separation can lead to experimental dynamos ( gailitis , 2000 ; busse , 2000 , stieglitz and mller , 2001 ) .
lanotte  ( 1999 ) showed that scale separation can also be useful for flows possessing parity invariance and thus lacking @xmath5-effect .
they found a first instance of a three - dimensional parity - invariant flow , the `` modified taylor  green flow '' , which can amplify a large - scale magnetic field by a negative eddy diffusivity mechanism , as conjectured by kraichnan ( 1976 ) .
examples of negative eddy diffusivities when the magnetic field depends on all three coordinates , while the flow depends only on two ( two and a half dimensional dynamos ) were actually known since the work of roberts ( 1972 ) .
new instances of large - scale dynamos with the roberts flow were found more recently by tilgner ( 1997 ) and plunian and rdler ( 2000 , 2001 ) .
emergence of the @xmath5-effect and the negative eddy diffusivity effect can be demonstrated by asymptotic theory involving , in principle , very large scale separation .
both the aforementioned experiments and numerical simulations with helical ( galanti , 1992 ) flows indicate that even a moderate scale separation of a factor two can help for obtaining a dynamo ( a physical interpretation for the case of the abc flows has been proposed by archontis ( 2000 ) ) .    the following argument due to fauve ( 1999 ) shows , how crucial for experimental dynamos it is to lower the critical magnetic reynolds number @xmath6 . here
, @xmath7 is the overall scale of the apparatus , @xmath8 is the typical velocity and @xmath9 is the ( molecular ) magnetic diffusivity . given the typically very small value of the magnetic prandtl number , the flow is expected to be strongly turbulent .
hence , the kinetic energy dissipated per unit time in the apparatus is @xmath10 , where @xmath11 is the fluid density and @xmath12 is the typical spatial scale of the velocity .
thus , @xmath13 , where @xmath14 is the scale - separation factor . therefore , for a given fluid ( @xmath9 ) , geometry ( @xmath15 ) and size ( @xmath7 ) , the required pump power behaves as the cube of the attainable @xmath16 .
it is generally believed that the cost of an experiment scales roughly as its volume , i.e. @xmath17 ; thus , if there are practical limitations on the maximum pump power available , the cost will scale as the ninth power of the attainable @xmath16 .
all this makes it desirable to design apparatus with the lowest possible critical magnetic reynolds number .
our intention is to begin a systematic investigation of the kinematic dynamo effect for parity - invariant flows with large or moderate scale separation with emphasis on ( i ) the genericity of the phenomenon of negative eddy diffusivity and ( ii ) the lowering of the critical magnetic reynolds number .
note , that parity - invariant flows are at least as `` natural '' as helical ones .
indeed , the navier - stokes equations are parity - invariant .
hence , even if this particular symmetry , like most other symmetries , is broken by the mechanisms generating the flow , we can expect the symmetry to be dynamically restored in a statistical sense at high reynolds numbers ( frisch , 1995 ) .
note also that the existence of three - dimesnional flow with negative eddy diffusivity is a recent discovery ; contrary to the alpha effect which can be derived in purely analytic fashion by working , e.g. , at low magnetic reynolds numbers or with a delta - correlated ( in time ) velocity , only ( very cpu - demanding ) numerical solutions of the equations arising from multi - scale techniques can demonstrate the existence of negative eddy - diffusivities and can show that it is not a pathologically unusual effect .
we shall work here with time - independent flows , which we may think of as being averages of turbulent flows .
we are of course aware that time dependence of the flow and , generally more realistic choices of the basic flow could somewhat modify our conclusions ; such issues are left for future work .
the paper is organized as follows .    in sections
[ s : expansion ] and [ s : abundance ] we consider the case where the spatial scale of magnetic field is much larger than that of the flow .
then @xmath15 is a natural small parameter , enabling one to apply asymptotic methods .
multiscale expansions in hydrodynamics may be viewed as adaptations of general homogenization methods used to determine effective large - scale transport coefficients of periodic or random materials ( bensoussan , 1978 ; kozlov , 1978 ; zhikov , 1979 ) .
usually the effective coefficients can not be obtained without solving suitable auxiliary small - scale problems . in three space dimensions
this may require significant computational resources . in lanotte  ( 1999 ) equations for the large - scale component of magnetic fields
are derived along the same lines as in previous works on negative eddy viscosity in ordinary incompressible hydrodynamics ( dubrulle and frisch , 1991 ; gama , 1994 ; wirth , 1995 ) .
here we use a similar but a distinct approach based on the work of vishik ( 1986 , 1987 ) , which yields a complete decomposition of magnetic modes and their eigenvalues ( in arbitrarily high - order terms of such expansions are derived ) .    to investigate ,
how commonly the phenomenon of negative magnetic eddy diffusivity is exhibited by fully three - dimensional flows , we follow in a strategy similar to the one , which gama  ( 1994 ) applied to estimate the abundance of two - dimensional flows with negative eddy viscosity . by solving numerically the relevant auxiliary problems we determine the minimum over directions of the ( anisotropic ) magnetic eddy diffusivities for two ensembles of flows with random fourier components
this way we determine the fraction of flows such that , when decreasing molecular magnetic diffusivity , a large - scale negative eddy diffusivity instability appears before the flow becomes a small - scale dynamo .
an interpretation of `` outliers '' with strongly negative eddy diffusivity , found in computations , is given in .
since it is important for our purpose to analyze as many instances as possible , we have used a new efficient algorithm for solving the elliptic auxiliary problems arising from the two - scale expansion .
this method is an outgrowth of the chebyshev iterative method developed by podvigina and zheligovsky ( 1997 ) for self - adjoint negative definite linear operators .
the operators arising in hydrodynamical problems are not self - adjoint , but they naturally split into a dominant self - adjoint negative definite part ( the molecular diffusion term ) and the remaining transport term .
we can then combine chebyshev methods with operator splitting to reduce the harm stemming from complex eigenvalues .
this method and its theoretical foundations are presented in detail in zheligovsky ( 2001 ) , available only in russian . for the reader
s convenience , the essentials are given in .
in we turn to cases with moderate scale separation ( i.e. when @xmath15 is finite ) .
first , we consider instances of flows for which we already know that there is a negative magnetic eddy diffusivity , and we determine windows in @xmath15 over which dynamo action is present .
second , we examine the case @xmath18 , that is , we assume that the magnetic field has a spatial period twice that of the flow , and study how the critical magnetic reynolds number @xmath16 is affected , compared to the case with no scale separation ( @xmath19 ) . in both types of computations we exploit the fact that , for space - periodic boundary conditions , eigenmodes have the form @xmath20 where @xmath21 is a prescribed wavevector .
this enables one to solve a problem involving only small scales , but parameterized by @xmath15 and @xmath21 ( cf .
roberts , 1972 ; frisch , 1986 ) .
solution of the eigenvalue problem is more cpu - intensive than the solution of auxiliary problems involved in the asymptotic case ( because efficient preconditioning techniques are available for the latter ) .
numerical methods and codes used for the eigenvalue computations are adapted from zheligovsky ( 1993 ) .
in we summarize our main results .
here we briefly review the asymptotic expansion arising in the kinematic dynamo problem , when the scale ratio @xmath15 is small .
the conducting fluid flow @xmath22 is assumed to be time - independent , and the problem therefore reduces to the eigenvalue problem for the magnetic induction operator @xmath23 in the space of solenoidal vector fields : @xmath24 here @xmath25 is magnetic molecular diffusivity , @xmath26 is a magnetic mode and @xmath27 is the associated eigenvalue ( the real part of which is here referred to as the growth rate of the mode ; negative growth rates correspond to decaying modes ) .
the flow is assumed incompressible : @xmath28 @xmath29-periodic in space and parity - invariant : @xmath30    the magnetic mode @xmath31 is supposed to depend on the fast variable @xmath32 and on the slow variable @xmath33 , and is @xmath29-periodic in both fast and slow variables ( the velocity @xmath34 is independent of slow variables ) .
introduction of the two spatial scales results in a `` decomposition of the gradient '' @xmath35 in and .
here and in what follows the subscripts @xmath32 and @xmath36 refer to differential operators in fast and slow variables , respectively .
let @xmath37 and @xmath38 denote the mean ( over the periodicity cube ) and the fluctuating part of a vector field , respectively : @xmath39 parity invariance implies @xmath40 .
solutions to the eigenproblem - are sought in the form of asymptotic series @xmath41 @xmath42 as shown in , substitution of these expansions in leads to a hierarchy of equations . by application of standard techniques for singular perturbation expansions , terms in the expansion eigenvalues
are successively determined from solvability conditions , obtained here by just averaging the equations .
in particular , it can be shown that @xmath43 , and the leading term @xmath44 in the series can be determined from the following eigenvalue problem @xmath45 subject to the solenoidality condition @xmath46    thus , the leading terms of the expansions in powers of @xmath15 of the mean magnetic field and of the eigenvalue turn out to be , respectively , an eigenfunction and the associated eigenvalue of a partial differential operator comprised only of second order derivatives , which may be regarded as representing an anisotropic diffusion . however , this operator is not necessarily dissipative : it may have eigenvalues with positive real parts . in this context one
then may speak of negative eddy diffusivity .    as shown in eqn .
, the coefficients @xmath47 of the homogenised equation are given by @xmath48 where the tensor field @xmath49 is obtained by solving a set of two auxiliary problems : @xmath50 and @xmath51 which will be referred to as the first and the second auxiliary problem , respectively . here
, @xmath52 is the unit vector in the @xmath53-th cartesian coordinate direction and superscripts denote cartesian components of vectors . for solutions to and to exist , we require that zero should not be an eigenvalue of the magnetic induction operator @xmath54 restricted to small - scale solenoidal vector fields with a vanishing mean ; by `` small - scale '' we mean independent of the slow variable . generically , this condition is satisfied .
problems and are formulated with a @xmath29-periodicity boundary condition for the vector fields @xmath55 and @xmath56 , which should also have a vanishing mean .
it follows that the @xmath57 s are parity - antiinvariant and solenoidal ( @xmath58 ) , and that the @xmath49 s are parity - invariant ( @xmath59  ) .
since the homogenized operator defined by - has constant coefficients , its eigenvectors are fourier harmonics : @xmath60 , where @xmath61 and @xmath21 are constant vectors .
the wavevector @xmath21 is taken real because we are only interested in modes uniformly bounded in the entire space . substitution in - yields @xmath62 @xmath63
following gama  ( 1994 ) , we now ask how common is the phenomenon of negative diffusivity .
time - independent periodic flows , satisfying - , are represented here by their fourier series @xmath64 with random amplitude harmonics . for the velocity to be real , @xmath65 , and
parity invariance requires re@xmath66=0 .
we study two ensembles of flows with different types of fall - offs of the energy spectrum @xmath67 .
the latter is defined as the total energy contained in harmonics , whose wavevectors @xmath68 belong to the @xmath69-th spherical shell , i.e. @xmath70 , where @xmath71 is integer . when @xmath72 , the flow is said to have exponential spectrum ; when @xmath73 , it is said to have hyperbolic spectrum .
these have respectively rapid and slow fall - offs .
the spectra are band - limited and set to zero beyond wavenumber @xmath74 for the exponential case and beyond @xmath75 for the hyperbolic case .
the value @xmath76 is used throughout in the exponential case .
the flows are produced by the following procedure ( applied only to one half of the set of wavevectors , the second one being related by complex conjugation ) . for @xmath68 in the ball @xmath77 ,
random imaginary vectors @xmath78 are generated with the imaginary part of each cartesian component uniformly distributed in the interval @xmath79 $ ] ( the ran2 random number generator by press , 1992 , is used ) .
then the vectors @xmath78 are projected onto the plane perpendicular to @xmath68 to gain incompressibility . finally , the obtained vectors are normalized , so that the flow has the desired energy spectrum and the root - mean - square velocity is one , when the normalized vectors are employed as the fourier coefficients @xmath66 in .
each ensemble consists of 100 flows . for each flow
, computations are carried out for three values of molecular diffusivity , @xmath80  , 0.2 and 0.1  .
note that when @xmath9 is too large , no dynamo action can take place at large or small scales because the eddy diffusivity will basically be equal to the molecular diffusitivity .
( in contrast , an alpha - effect dynamo is still possible if the system is large enough ) .
if @xmath9 is too small , small - scale dynamo instabilities will have much larger growth rates than large - scale instabilities . in both extreme cases it is meaningless to study negative eddy diffusivity large - scale instabilities .
the range of values of @xmath9 chosen , between 0.1 and 0.3 , is an empirically determined tradeoff between these extremes .
all computations are performed with the resolution of @xmath81 fourier harmonics , employing standard spectral methods based on fast fourier transforms and dealiasing .
( each case requires from 10 to 30  minutes of cpu on one dec alpha  ev6 processor , depending on the value of @xmath9 . ) the power spectra of the solutions to and turn out to decrease monotonically very quickly for @xmath82 . for @xmath83 in the hyperbolic case ( resp . exponential case ) , the decrease is typically by 7 - 8 orders of magnitude ( 14 - 15 orders of magnitude ) from maximum to minimum .
note that our code reproduces the negative eddy diffusivity values obtained by an alternative numerical method by lanotte  ( 1999 ) for the modified taylor  green flow .
100 10 eddy diffusivity histograms    table 1 .
statistics of the eddy diffusivities for the considered ensembles of flows .    [ cols="^,^,^,^,^,^,^ " , ]     after the homogenized coefficients are determined , the minimum magnetic eddy diffusivity @xmath84 is calculated as the minimum eigenvalue of the problem - .
( note , that for a unit wavevector the eigenvalue is just the opposite of the eddy diffusivity . ) for a given @xmath21 this eigenproblem yields two eigenvalues satisfying a quadratic equation ; the eigenvalues can thus be complex conjugate ( a similar phenomenon was observed for three - dimensional eddy viscosities ; see wirth , 1995 ) .
then the minimum of the real part of the complex eddy diffusivity , which determines the maximum rate of instability of a large - scale magnetic perturbation , is sought . however , in our computations complex eddy diffusivities appear quite rarely .
histograms of the obtained magnetic eddy diffusivities are shown in .
the statistics of eddy diffusivities shown in table  1 indicate that the lower the molecular diffusivity , the higher the chance that a flow will have a negative `` eddy correction '' ( correction to the molecular value coming from the flow ) and also the higher the chance that the resulting eddy viscosity will be negative .
note however that the eddy correction ( given , at least implicitly , by - ) has a complicated , generally non - monotonic , nonlinear dependence on the molecular diffusivity .
table  1 also shows that flows with the exponential energy spectrum are more efficient in lowering the minimum of the eddy diffusivity than those with the hyperbolic spectrum .
this suggests that low wavenumber harmonics of the underlying flow are more important for lowering than high wavenumber harmonics .
examples of magnetic eigenmodes for @xmath83 are shown in figures  [ f : sample1 ] and [ f : sample2 ] .
they correspond to two samples having exponential velocity spectra .
the upper of each figure shows regions where @xmath85 is less than 0.2% of its maximum , hence revealing the immediate neighborhood of stagnation points .
note that in consequence of parity invariance about the origin and of @xmath29-periodicity , within a basic @xmath29-periodicity cube there are in total eight points about which the flow is parity - invariant .
these stagnation points are located at the vertices of a cube with side @xmath86 and they form a uniform @xmath86-periodic grid in space .
the magnetic structures seen in the lower left and right of each figure are clearly located near some of those stagnation points .
these structures are magnetic flux ropes of the kind considered by galloway and zheligovsky ( 1994 ) .
note that only the fluctuating part of the magnetic eigenmodes are shown ( both the periodicity - box averaged magnetic field @xmath87 and the overall slow spatial factor @xmath88 are omitted ) . in @xmath87
is very small ( about two orders of magnitude smaller than the fluctuating field ) , whereas for the case of it is of the same order of magnitude as the fluctuating field .
100 10 eigenmode sample 1 ( exponential spectrum )    100 10 eigenmode sample 2 ( exponential spectrum )    100 10 eddy diffusivity vs growth rate ( exponential spectrum )    100 10 eddy diffusivity vs growth rate ( hyperbolic spectrum )    a physically important question is : which instability arises first , when magnetic molecular diffusivity is decreased , a large - scale or a small - scale one ?
similarly to what was done by gama(1994 ) , for each considered pair of a flow and a molecular diffusivity , it is checked that small - scale magnetic fields with a vanishing mean are stable , i.e. have negative growth rate ( see figures  [ f : dotsntrianglese ] and  [ f : dotsntrianglesh ] ) .
100 10 infinitely low eddy diffusivity    several `` outliers '' with relatively large negative eddy diffusivities can be seen on the histograms shown on for @xmath83  .
this can be explained as follows .
we see on figures  [ f : dotsntrianglese ] and  [ f : dotsntrianglesh ] that in the case of outliers the growth rate of _ small - scale _ magnetic modes ( plotted horizontally ) is negative , but close to zero .
furthermore , we have checked that the corresponding eigenvalue is real . in other words , we are close to a non - oscillatory small - scale instability threshold .
now , we observe that the expression of the eddy diffusivity , derived in and , involves the inverse of the small - scale magnetic induction operator .
it is thus not surprising that large eddy diffusivities can be obtained .
depending on whether a certain inequality is satisfied at the threshold , this eddy diffusivity can be either large negative or large positive ( see details in ) .
we have checked numerically that the outlier eddy diffusivities are approximately inversely proportional to the small dominant eigenvalue ( see ) .
note , that if we are exactly at the threshold , our two - scale expansion breaks down .
a modified expansion can be carried out which takes into account the enlargement of the kernel of the magnetic induction operator @xmath89 . in the presence of a large - scale perturbation
this can lead to an eigenvalue @xmath90 rather than the @xmath91-dependence obtained in the general parity - invariant case .
in other words , the large - scale dynamics can be governed by a first - order differential operator .
for intermediate values of @xmath15 , corresponding to moderate scale separation , direct evaluation of the most unstable magnetic modes is desirable .
computations are simplified by the fact that eigenmodes have the form , all the dependence on the slow variable @xmath33 being absorbed into an exponential factor @xmath92 .
since the flow is assumed to depend only on the fast variable @xmath32 , the eigenvalue problem is reduced to a small - scale one for a _ modified operator _ by changing @xmath93 in the magnetic induction operator and in the solenoidality condition .    to find the most unstable magnetic mode one must , in principle , allow for arbitrary orientations of the vector @xmath21 .
throughout this section we assume only `` binary '' @xmath21 s such that each component takes only the values 0 or 1 , so that @xmath94 may be regarded as the maximum scale - separation factor over cartesian directions .
since our simulations are computationally demanding , we choose for each flow only one @xmath21-direction as follows . for a given flow and for the molecular diffusivity @xmath83
, we obtain the most unstable @xmath21-directions _ in the asymptotic regime _
@xmath95 from the results of .
it turns out that for the cases reported in this section , these directions are approximately parallel to binary vectors with one or two non - zero components .
these binary vectors are then employed when @xmath15 is finite for all values of @xmath9 studied .
100 10 growth rate vs scale ratio ( a )    100 10 growth rate vs scale ratio ( b )    100 10 growth rate vs scale ratio iii    100 10 growth rate vs scale ratio iv    figures  [ f : growthvsepsilona]-[f : growthvsepsilond ] show the dominant growth rates for four different flows selected out of the ensemble of flows with hyperbolic spectrum , which sustain negative eddy diffusivity .
actually , when @xmath21 has integer components , it is only necessary to make computations for @xmath96 .
indeed , first note , that the modified magnetic induction operators for two values of @xmath15 differing by 1 have the same spectrum and their eigenfunctions differ by a factor @xmath97 .
second , by parity invariance of the velocity , for two opposite values of @xmath15 the spectra are also the same and eigenfunctions are related by parity @xmath98 .
( from this one easily deduces that the spectrum of the modified operator is complex conjugate , eigenvalues @xmath27 and @xmath99 being associated to eigenfunctions @xmath26 and @xmath100 , respectively . )
hence , outside the interval @xmath96 the plots can be continued by symmetry about the vertical line @xmath18 and by 1-periodicity .
the results reported here are obtained with the resolution of @xmath81 fourier harmonics .
for the `` worst '' case , @xmath101 , they differ only in the third significant digit from those obtained with the resolution @xmath102 .
note , that energy spectra of the computed magnetic modes decrease by at least 4 orders of magnitude for @xmath103 and by at least 5 for @xmath80  , indicating that the employed resolution is adequate .
( for @xmath101 the decrease is only by 3 orders of magnitude , which may seem marginal .
however , it is well - known that eigenvalues are less sensitive to insufficient resolution , than the fine structure of eigenmodes . )
it is clear from figures  [ f : growthvsepsilona]-[f : growthvsepsilond ] that for all instances where the eddy diffusivity is negative ( i.e.  the graph go to zero by negative values as @xmath104 ) there is also an instability for any finite @xmath15 .
the converse is , however , not always true ; for example , for the last flow ( figure  [ f : growthvsepsilond ] ) and @xmath105 there is a dynamo for @xmath106 , which disappears below @xmath107 ( the eddy diffusivity is positive ) .
furthermore , inspection of figures  [ f : growthvsepsilonb]-[f : growthvsepsilond ] reveals that there can be one or more windows in scale separation , where the dynamo ( positive growth rate ) is present , whereas in figures  [ f : growthvsepsilona ] the dynamo appears to be present for all scale separations , provided the eddy diffusivity is negative . in the two windows
belong to two different analytical branches of eigenvalues and eigenmodes , with the transition occurring near @xmath108  .
note , that all the growth rates vanish ( quadratically ) when @xmath95 .
the reason is that the small-@xmath15 eigenvalues , obtained by the two - scale expansion of , are analytic continuations of the zero eigenvalue in the small - scale problem .
the presence of the eigenvalue zero is due to the structure of the magnetic induction operator , which is the curl of a first - order operator .
( hence , its adjoint , acting also in the space of divergenceless functions @xmath109 has constant vectors in its kernel . )
conversely , it can be demonstrated by spatial averaging of @xmath110 , that for a mode with a non - vanishing mean the associated eigenvalue equals 0 .
( note , that in we restricted the operator to functions of vanishing spatial mean and thus the growth rate zero was absent from figures  [ f : dotsntrianglese ] and [ f : dotsntrianglesh ] . )    100 10 critical rm s : small - scale vs double period    we now turn to the case , where there is a finite scale separation of a factor two ( @xmath18 ) .
the question we address is : how does a mere doubling of the spatial period of magnetic modes ( relative to that of the basic flow ) affect generation ?
more precisely : to lower the critical magnetic reynolds number , does it help to introduce such a scale separation ?
we systematically compare critical numbers when @xmath19 ( no scale separation ) and @xmath18 for an ensemble of flows with hyperbolic energy spectrum .
we define our magnetic reynolds numbers @xmath6 in terms of the _ largest _ scale available @xmath111 ( the inverse of the minimum wavenumber of the magnetic field ) and of the r.m.s . velocity @xmath8 . in computations
@xmath112 and hence @xmath113 .
when @xmath18 we face again the problem , mentioned earlier , of choosing the vectors @xmath21 . for each flow
we now search for the critical @xmath16 , allowing all 7 possible binary @xmath114 .
more precisely , we find an approximate ( minimum ) critical @xmath16 and the corresponding @xmath21 by running our code with @xmath102 fourier harmonics , and then refine the critical @xmath16 with @xmath81 fourier harmonics for the obtained @xmath21 .
the secant method is used to find stationary modes .
iterations arere terminated when the real part of the dominant eigenvalue was below @xmath115 in magnitude .
this usually requires 4 - 5 iterations .
as far as the accuracy of the calculation is concerned , note that in all @xmath18 runs with @xmath116 the energy spectrum decreases by at least 4 orders of magnitude , an indication that @xmath81 harmonics provide sufficient resolution . in the most demanding runs , with @xmath19 and @xmath117 , the energy spectra decrease only by 23 orders of magnitude but ,
as pointed out before , eigenvalues are computed nevertheless quite accurately ( critical @xmath16 s , computed with the @xmath102 and @xmath81 resolutions , differ in the third significant digit )
.    shows all the critical @xmath16 s calculated for 28 cases .
two cases not shown have critical @xmath118 and @xmath119 when @xmath18 , but do not display dynamo action for any @xmath120 when @xmath19 .
the results demonstrate that scale separation lowers the critical @xmath16 in all but 9 cases out of the 30 studied .
note that nore  ( 1997 ) have studied a flow ( a variant of taylor  green ) for which a scale separation @xmath18 seems to lower significantly the critical magnetic reynolds number , as is the case for the helical abc flows studied by galanti  ( 1992 ) , archontis ( 2000 ) and dorch ( 2000 ) .
however the nore  ( 1997 ) definition of the magnetic reynolds number uses the scale of the flow and not that of the magnetic field .
if instead our definition is used , their flow appears to be just a borderline case .
three main results are presented in this work .
first , in , we demonstrated that parity - invariant flows possessing a negative magnetic eddy diffusivity are quite common .
specifically , we generated two random ensembles of flows ; both had time - independent randomly chosen fourier components , one with an exponentially decreasing energy spectrum , the other  with a power - law decrease ( referred to as `` hyperbolic '' because the exponent is -1 ) .
one hundred independent selections were made in each case and three different molecular magnetic diffusivities were used , all corresponding to small - scale magnetic reynolds numbers not exceeding ten .
none of the 600 cases investigated was a small - scale dynamo .
for the lowest diffusivity considered , 86% of the `` exponential '' spectrum flows and 53% of the `` hyperbolic '' spectrum flows produced a negative eddy diffusivity and , hence , can act as a large - scale dynamo ( for more details , see table  1 ) .
second , several instances with particularly negative eddy diffusivities were encountered .
the theory of these `` outliers '' ( given in details in ) relates them to a near - critical state , where an eigenvalue of the small - scale dynamo problem is about to pass through zero ( there is also an additional inequality ensuring that a large _ negative _ value is obtained ) .
finally , we presented numerical evidence that , for a parity - invariant flow , even a moderate scale separation between the flow and the magnetic field is generally beneficial for magnetic field generation . specifically , we showed that the critical magnetic reynolds number based on the largest scale present ( that of the magnetic field ) is often lowered when the scale of twice the spatial period of the velocity is introduced ( 70% of the instances studied had this property ; see ) .
for a few flows we also studied the dependence on the scale - separation factor , when it is varied from a 1:2 ratio to a 1:20 ratio .
one instance ( shown in ) has a window lacking dynamo action between low and large values , where it is present .    for helical flow , such as abc
, the lowering of the critical magnetic reynolds number , when a moderate separation of scale is assumed , is understood in terms of symmetry breaking , stagnation points and other geometrical features ( archontis 2000 ) .
it would be of interest to find which geometrical features are important for the lowering mechanism when the flow is parity - invariant .    * acknowledgments*.    we are grateful to the authors of lanotte  ( 1999 ) for making available to us the source of the code used to obtain the results reported in their paper , to s.  fauve , s.  gama , a.  lanotte , a.  noullez , y.  ponty and m.  vergassola for useful discussions and to an anonymous referee for several suggestions .
computational facilities were provided by the program `` simulations interactives et visualisation en astronomie et mcanique ( sivam ) '' at the observatoire de la cte dazur .
visits of o.m .
podvigina and v.a .
zheligovsky were supported by the french ministry of education .
this work was also supported by the european union under contract hprn - ct-2000 - 00162 .
[ a : expansion ]    we construct here complete asymptotic expansions of magnetic eigenmodes and of their eigenvalues for small values of the scale ratio @xmath15 .
assumptions concerning the flow are listed in .
the kinematic dynamo problem reduces to the eigenvalue problem @xmath121 here @xmath26 is a magnetic mode , @xmath27 is its eigenvalue , and @xmath122 acting in the space of solenoidal vector fields , is the magnetic induction operator ( it acts like @xmath89 , but in a larger domain ) .
the magnetic mode @xmath31 is supposed to be solenoidal and @xmath29-periodic in both fast and slow variables .
solutions to the eigenproblem are sought in the form of asymptotic series @xmath123 and @xmath124 .
denote @xmath125 where @xmath37 and @xmath38 refer to the mean and fluctuating part ( see  ) , respectively .
the expansion is constructed under the assumption that zero is not an eigenvalue of the operator @xmath89 acting in the space of small - scale periodic solenoidal vector fields with vanishing mean .
( by `` small - scale '' we understand `` depending only on the fast variables '' . )
the decomposition of the gradient @xmath126 , applied to the solenoidality condition , expanded in powers of @xmath15 and separated into mean and fluctuating part , gives @xmath127 @xmath128 for all @xmath129 .
( @xmath130 is assumed for all @xmath131 . )    substituting the expansions and the decomposition of the gradient into , we obtain : @xmath132\epsilon^n=0 .
\label{serieseqn}\end{aligned}\ ] ] we proceed by successively equating to zero the mean and the fluctuating part of each term of this series . in our context , averaging is equivalent to obtaining the solvability condition for a partial differential equation in fast variables .
transform the leading @xmath134 term of using vector algebra identities and : @xmath135 averaging of this equation yields @xmath136 .
then the solution of can be expressed as @xmath137 where the vector fields @xmath57 are the solutions of the first auxiliary problem ; the solutions exist since it is assumed that the kernel of @xmath89 is empty .
( the first and second auxiliary equations were stated in . )    by taking the divergence of the first auxiliary equation , one obtains @xmath138 . since the velocity is parity - invariant , the domain of @xmath89 is a direct sum of two invariant subspaces , one of which is comprised of parity - invariant vector fields ( i.e.  @xmath139  ) , and the second  of the parity - antiinvariant ones ( @xmath140  ) .
thus , @xmath141 is a parity - antiinvariant matrix : @xmath142    @xmath143 .
the @xmath144 term of , after the use of , takes the form @xmath145 & & -({\bf v}\cdot\nabla_{\bf y}){\bf h}_0=\lambda_1({\bf h}_0+{\bf g}_0 ) .
\label{eigeneq1}\end{aligned}\ ] ] substituting in and averaging , we obtain @xmath146    so far , developments strictly followed those of vishik ( 1987 ) , but here the point of divergence is reached .
since @xmath34 is parity - invariant and @xmath141 is parity - antiinvariant , we have @xmath147 , and therefore implies @xmath148 . upon inserting into the fluctuating part of , we obtain @xmath149 from which one finds @xmath150 here , the nine vector fields @xmath49 are solutions to the second auxiliary equation . note @xmath59 by the symmetry properties of @xmath34 , @xmath57 and of the operator @xmath89 . by taking divergence of the second auxiliary equation
one finds @xmath151 .
hence , for @xmath152 is automatically satisfied .
the @xmath154 term of takes the form @xmath155 & & + \,\,\nabla_{\bf y}\times({\bf v}\times{\bf g}_1 ) -({\bf v}\cdot\nabla_{\bf y}){\bf
h}_0+{\bf g}_0 ) . \label{eigeneq2}\end{aligned}\ ] ] upon substitution of , averaging produces the second - order equation @xmath156 which defines the eddy diffusivity operator .
so , the leading terms of the expansions of the mean magnetic field , @xmath157 , and of the associated eigenvalue , @xmath44 , satisfy the eigenvalue equation for this operator .
its eigenvectors are fourier harmonics : @xmath158 , where @xmath61 and @xmath21 are constant vectors satisfying and .
note that the fluctuating part of the leading term of the eigenmode , @xmath159 , is now completely determined by .    at this point
our construction goes beyond that of lanotte  ( 1999 ) and , in the spirit of vishik ( 1987 ) , we show how to extend the analysis to higher orders in @xmath15 .    the fluctuating part of reads @xmath160 @xmath161 @xmath162 here , only the terms involving @xmath163 and @xmath164 are not yet determined .
however , as one can immediately see , @xmath32-dependent prefactors in front of these terms are identical to those , which are in in front of @xmath165 and @xmath166 , respectively .
this implies a representation @xmath167 where the vector field @xmath168 can be uniquely determined from @xmath169 @xmath170 @xmath171 .
as already stated , @xmath44 is one of the two eigenvalues of - . in what follows
we assume that the other eigenvalue , @xmath172 , is distinct from @xmath44 and we denote by @xmath173 the associated eigenvector .
equations arising at higher orders in @xmath15 are now solved recursively .
assume that all the equations up to order @xmath174 have been solved , so that the following information has been obtained :    @xmath175 vector fields @xmath176 and @xmath177 for all @xmath178 ;    @xmath175 representations of @xmath177 of the form @xmath179 for @xmath180 and @xmath181 with known vector fields @xmath182 , @xmath183 ;    @xmath175 quantities @xmath184 for @xmath185 .
upon substitution of for @xmath180 the average of the equation at order @xmath186 becomes @xmath187 where the right - hand side at this stage is a known vector field .
projecting this equation out in the direction of @xmath157 one can uniquely determine @xmath188 .
in the complementary invariant subspace the operator @xmath189 is invertible , and hence @xmath190 can be determined up to an arbitrary multiple of @xmath157 .
this additive term can be neglected : one can demand that @xmath190 belong to the complementary subspace .
this is a normalization condition related to the fact that eigenvectors can be multiplied by arbitrary ( analytic ) functions of @xmath15 . as a result ,
@xmath191 is now determined by with @xmath181 .    upon use of for @xmath180
the fluctuating part of the equation at order @xmath186 becomes @xmath192 @xmath193 @xmath194 @xmath195 only the terms involving @xmath196 and derivatives of @xmath197 are not yet determined .
like in , @xmath32-dependent prefactors in front of these terms are identical to those in front of @xmath165 and derivatives of @xmath198 in , respectively .
hence , yields a representation of the form for @xmath199 .
an equation for @xmath200 can be obtained from by replacing @xmath199 by @xmath200 and omitting all terms involving @xmath196 and derivatives of @xmath201 .
thus , a complete formal decomposition of magnetic modes and their eigenvalues has been constructed .
a step - by - step analysis of the solution reveals that @xmath202 ( with suitable scalars @xmath203 ) and @xmath204 thus the expanded eigenmode admits a representation in the form of a plane wave in the slow variable .
( for this reason it was sufficient to demand , for constructions presented under heading @xmath171 , that the two eigenvalues @xmath44 and @xmath172 of @xmath205 be distinct in the subspace of solenoidal vector fields @xmath206 with @xmath207const , instead of requiring that @xmath189 be invertible in the whole domain . )
the plane - wave representation does not come as a surprise , since the domain of the magnetic induction operator @xmath208 splits into invariant subspaces , each comprised of vector fields of the form and categorized by wavevectors @xmath21 .    if the homogenized operator @xmath205 turns out to be elliptic , one can prove , following vishik ( 1986 , 1987 ) , that any point @xmath44 of the spectrum of @xmath205 is associated with an analytic branch of eigenvalues of the original magnetic induction operator @xmath208 .
the proof relies on the general perturbation theory for linear operators ( kato , 1980 ) ; it will not be presented here .
[ a : loweddydiff ]    in this appendix we discuss a mechanism for appearance of negative eddy diffusivities with arbitrarily large magnitudes .
let @xmath209 denote the basis of small - scale @xmath29-periodic magnetic modes with vanishing mean and @xmath210  the associated eigenvalues : @xmath211 ( assuming in this section , for the sake of simplicity , that @xmath89 is diagonalizable ) .
suppose that when molecular diffusivity @xmath9 is decreased , small - scale magnetic fields with vanishing mean lose stability in such a way that the dominant eigenvalue of the magnetic induction operator @xmath89 passes through zero .
let @xmath212 be the critical diffusivity : @xmath213 generically , the loss of stability occurs either in the proper subspace comprised of parity - invariant fields , or in that of parity - antiinvariant fields , implying two possible variants of the mechanism .
suppose small - scale parity - invariant magnetic fields are the first to lose stability .
solutions @xmath214 to the first auxiliary problem continuously depend on @xmath9 for @xmath215 .
the second auxiliary problem can be reexpressed as @xmath216 where @xmath217 and @xmath218 are decompositions of the r.h.s . of and of @xmath219 in the basis of magnetic modes ; summation over all parity - invariant eigenvectors
is assumed in the two sums . from , @xmath220 .
since generically the coefficients @xmath221 do not vanish , we have in the limit @xmath222 @xmath223 then the equation for the eigenvalue @xmath44 asymptotically reduces to @xmath224 where @xmath225 . equations and yield two eigenvalues .
it is easily checked that one of them , denoted @xmath226 , vanishes and thus is of no interest for us ; the associated eigenvector @xmath173 is orthogonal to the wavevector @xmath21 and satisfies @xmath227 the second eigenvalue is @xmath228 the associated eigenvector being @xmath229 ( we consider a generic case , where @xmath230 does not vanish ) .
therefore , if @xmath231 , then , in the limit @xmath222 , the minimum magnetic eddy diffusivity becomes arbitrarily large negative : @xmath232 alternatively , suppose now that the loss of stability occurs in the subspace of small - scale parity - antiinvariant magnetic fields .
then the mechanism is similar . in this case already
the solutions @xmath214 to the first auxiliary problem tend to infinity .
indeed , can be represented as @xmath233 where @xmath234 and @xmath235 are decompositions of the r.h.s . of the first auxiliary problem and of @xmath236 in the basis of magnetic modes ; summation over parity - antiinvariant eigenvectors is assumed in both sums .
thus , @xmath237 .
since generically @xmath238 , we have in the limit @xmath222 @xmath239 let @xmath240 denote the solution to @xmath241 then @xmath242 and asymptotically reduces to @xmath243 where @xmath244 equations and yield two eigenvalues .
again , one can easily check that one of them , @xmath226 , vanishes ; the associated eigenvector @xmath173 is orthogonal to @xmath21 and satisfies @xmath245 the second is @xmath246 with the associated vector @xmath247 ( in the generic case @xmath248 ) .
therefore , if @xmath249 , when @xmath222 the minimum magnetic eddy diffusivity becomes again arbitrarily large negative : @xmath250
[ a : chebyshev ]    we discuss here three chebyshev iterative methods suitable for large systems of linear or non - linear equations .
they employ extremal properties of chebyshev polynomials and their roots . the first method ( ) , discussed in ,
is applicable to problems with negative - definite self - adjoint jacobians .
an extension ( ) , discussed in , is an upgrade of , which eliminates these restrictions . the second method ( ) , discussed in ,
is designed for discretizations of elliptic differential equations , such as the auxiliary problems and .  and
were introduced by podvigina and zheligovsky ( 1997 ) and  is presented in detail in zheligovsky ( 2001 ) . for general background on chebyshev iterative methods ,
see axelsson ( 1996 ) .
let @xmath254 be a solution to .
we first consider the case , where the jacobian matrix @xmath255 is diagonalizable and its spectrum is ( real and ) negative .
then the solution can be computed by iterations @xmath256 @xmath257 provided @xmath258 is chosen in the basin of attraction of @xmath259 . here
@xmath260 where @xmath261 for any @xmath262 ( the optimal value @xmath263 is usually unknown ) , and where the @xmath264 s are chebyshev polynomials of the first kind .
the @xmath265 s can be evaluated recursively : @xmath266    consider how the discrepancy vectors @xmath267 evolve in the course of iterations - .
introduce the eigenvalues @xmath268 of @xmath255 ( here,@xmath269 ) and the associated eigenvectors @xmath270 , and decompose the initial discrepancy in this basis : @xmath271 it is then easily shown that @xmath272 ( the @xmath273 correction does not appear if the mapping @xmath274 is linear ) .    at this point ,
the standard strategy is to use the least deviation principle for chebyshev polynomials : the iterative process - with @xmath275 ensures maximum uniform decay of the discrepancy vector component in the proper subspace of @xmath255 associated with the part of the spectrum in the interval @xmath276 .
the alternative procedure used by podvigina and zheligovsky ( 1997 ) relies , in addition , on the extremal property of the largest root of a chebyshev polynomial : for @xmath277 , the iterative process gives optimal reduction of the discrepancy vector component in the proper subspace of @xmath255 associated with eigenvalues @xmath268 closest to zero .
therefore , a suitable algorithm , called here , consists of a succession of iterative sequences - with @xmath275 and @xmath277 .    as usually , lowering of the condition number of the system allows to improve efficiency of the algorithm .
if , in addition , @xmath255 is self - adjoint , such lowering can be achieved by replacing with an equivalent system @xmath278 where @xmath279 @xmath280 being a self - adjoint and positive - definite linear preconditioning operator .
it may be checked that for such a @xmath280 the mapping @xmath281 has the required properties stated at the beginning of this subsection .      in the case of an arbitrary smooth mapping @xmath274 , to which we turn now , a modified system @xmath282 can be considered , where @xmath283 @xmath284 being arbitrary positive definite self - adjoint operators and@xmath285 being the operator adjoint to @xmath286 .
the method discussed above can be applied for the solution of the modified system , since @xmath287 has all the required properties .
\2 . in any of the following iterative sequences
@xmath277 , unless @xmath27 was zero in the previous sequence and that sequence was short ( i.e. it consisted of less than @xmath289 $ ] iterations , where @xmath290 $ ] denotes the integer part of a number ) , in which case in the next sequence @xmath288 .
\4 . in an iterative sequence with @xmath277 the errors @xmath291 can initially grow monotonically , if the respective errors @xmath292 decrease monotonically .
( an increase of @xmath293 at the first iteration indicates that @xmath294 is underestimated . )
\5 . after the initial phase allowed by rules 3 or 4
is completed , an iterative sequence is terminated if both discrepancy reduction rate values , @xmath295 and @xmath296 , are smaller than at the previous iteration .
we present here an alternative to , better suited to problems arising from discretization of elliptic partial differential equations , when they involve not just leading - order elliptic operators ( such as the laplacian ) , but also lower - order non - self - adjoint perturbations .
we do not attempt to present this method in full generality ( see zheligovsky , 2001 ) .
in what follows the resolution of the discretization , assumed to be large , is the control parameter , which we denote by @xmath303 .
construction of the modified mapping can be avoided under the conditions , that the basic mapping @xmath304 can be split : @xmath305 and if a linear invertible preconditioning operator @xmath306 can be selected , such that :    \1 .
all eigenvalues @xmath307 of the linear operator @xmath308 are negative .
( if @xmath306 is a self - adjoint positive - definite operator and @xmath309 is a self - adjoint negative - definite operator , this condition is automatically verified . )
the ratio @xmath310 is small for large @xmath303 . here
@xmath311 , and @xmath294 is a `` tight '' upper bound for the spectral radius of @xmath312 , i.e.  a bound , which stays of the order of the spectral radius for large resolutions .
the operators @xmath313 and @xmath314 nearly commute : @xmath315 .
( note , that this condition is not that restrictive : it holds automatically in the absence of preconditioning , if @xmath313 is for example the laplacian . )
a decomposition arises naturally if the problem is a discretization of a dissipative partial differential equation : then the mapping @xmath316 can be identified with the higher - order elliptic part of the equation , and @xmath230  with the remaining terms .
auxiliary problems and arising from the two - scale expansion fall into this class .    under the conditions stated above a modification of the method discussed in , named ,
can be applied .
solution proceeds by a succession of chebyshev iterative sequences - , but the first iteration in all sequences is made differently .
suppose a termination condition ( for example , the one considered in the previous subsection ) for a chebyshev sequence number @xmath317 is verified for the first time in this sequence for the iterate @xmath318 ( we then call @xmath69 the termination signal number ) . then the sequence is continued by @xmath319 further iterations , and we start the next sequence , by selecting @xmath320 subsequent iterations are defined as in : @xmath321 ( the quantities @xmath322 and @xmath265 were defined by - and @xmath281  by . )    in this iterative process , errors behave in the following way .
discrepancy vectors @xmath323 can be represented in terms of chebyshev polynomials of the first ( @xmath264 ) and of the second ( @xmath324 ) kind : @xmath325 where @xmath326 , @xmath327 is the identity operator , and vectors @xmath328 and @xmath329 , which we call the basic error vectors , can be determined from @xmath330 and @xmath331 .
( the terms @xmath332 do not emerge if the mapping @xmath274 is linear ; in what follows they are neglected . )    now , suppose that for the ( @xmath333)-nd sequence the termination signal number @xmath69 was the same as for the ( @xmath317)-st iteration ( in computations , successive chebyshev sequences often do have equal lengths )
. then it can be shown that      @xmath175 within the succession of iterative sequences , the basic error vectors evolve as follows : @xmath335 where @xmath336 @xmath337 , @xmath307 being eigenvalues of @xmath312 , and the coefficients @xmath338 can be found from the decomposition of @xmath339 in the basis of eigenvectors @xmath340 of @xmath312 : @xmath341 thus , reduction of the basic error vectors @xmath328 is predominantly controlled by @xmath312 ( note that @xmath342 ) . the standard first chebyshev iteration is recovered when one sets @xmath343 and hence @xmath344 . it may be shown that our procedure leads to factors @xmath345 , and thus to enhanced convergence , for eigenvalues sufficiently small in absolute value ( provided @xmath277 or @xmath27 is small ) . for the other eigenvalues good convergence
is ensured as usually .
( actually , because of the presence of the non - self - adjoint perturbation @xmath346 , fast convergence is guaranteed only for eigenvalues large in absolute value ; for those of intermediate magnitude slowdown might occur , but has not been observed in any of the cases we studied . ) of course , these convergence arguments hold only in so far as the resolution @xmath303 is high enough and thus @xmath347 is small enough to justify keeping only the leading terms in .    in the application of the generalized chebyshev method to the auxiliary problems and , we take @xmath348 ( more precisely , @xmath349 is the fourier
galerkin discretization of this operator ) .
we select the preconditioning operator @xmath350 for some @xmath351 by analogy with podvigina and zheligovsky ( 1997 ) and podvigina ( 1999a , b ) .
then @xmath352 is a pseudodifferential operator of order @xmath353 , and @xmath346  of order @xmath354 .
hence , condition 2 is inherited from the non - discretized versions of the operators for @xmath355 .
condition 3 is satisfied for @xmath356 .
indeed , in the commutator @xmath357 the leading symbols cancel out and therefore its order is by 1 less than the order of each of the operators @xmath358 and @xmath359 ; thus @xmath360 , @xmath361 , where @xmath303 is the maximum wavenumber of the fourier harmonics , retained in the abridged series , this completing the argument .
since condition 3 is also trivially satisfied for @xmath362 , it is natural to consider the larger interval @xmath363 , and we find the convergence to be the best for @xmath364 .
when applied to and ,  proved superior to .
gailitis a. , lielausis o.a . , dementev s. , platacis e. , cifersons a. , gerbeth g. , gundrum t. , stefani f. , christen m. , hanel h. and will g. `` detection of a flow induced magnetic field eigenmode in the riga dynamo facility '' , _ phys
lett . _ * 84 * , 43654368 ( 2000 ) .
podvigina o.m . and zheligovsky v.a .
`` an optimized iterative method for numerical solution of large systems of equations based on the extremal property of zeroes of chebyshev polynomials '' , _ j. sci . computing _ * 12 * , 433464 ( 1997 ) .
steenbeck m. , krause f. , rdler k .- h . `` a calculation of the mean electromotive force in an electrically conducting fluid in turbulent motion , under the influence of coriolis forces '' , _
z. naturforsch _ * 21a * , 369376 ( 1966 ) .
stieglitz r. and mller u. `` an experimental demonstration of a homogeneous two scale dynamo '' , _ proceed .
nato arw `` dynamo and dynamics , a mathematical challenge '' _ , august 21 - 26 , 2000 , institut dtudes scientifiques de cargse , to appear ( 2001 ) .
zheligovsky v.a .
`` a chebyshev iterative method with operator splitting for evaluation of roots of large systems of equations '' ( in russian ) , accepted _ proc .
international workshop on nonlinear problems of hydrodynamic stability , moscow , 2000_. institute of mechanics , lomonosov moscow state university , moscow univ . press ( 2001 ) . also available at + http://www.obs-nice.fr/etc7/vlad.ps.gz