Computation of Kinematic and Magnetic <i>α</i>-Effect and Eddy Diffusivity Tensors by Padé Approximation

Abstract

We present examples of Padé approximations of the α -effect and eddy viscosity/diffusivity tensors in various flows. Expressions for the tensors derived in the framework of the standard multiscale formalism are employed. Algebraically, the simplest case is that of a two-dimensional parity-invariant six-fold rotation-symmetric flow where eddy viscosity is negative, indicating intervals of large-scale instability of the flow. Turning to the kinematic dynamo problem for three-dimensional flows of an incompressible fluid, we explore the application of Padé approximants for the computation of tensors of magnetic α -effect and, for parity-invariant flows, of magnetic eddy diffusivity. We construct Padé approximants of the tensors expanded in power series in the inverse molecular diffusivity 1 / η around 1 / η = 0 . This yields the values of the dominant growth rate to satisfactory accuracy for η , several dozen times smaller than the threshold, above which the power series is convergent. We do computations in Fortran in the standard “double” (real*8) and extended “quadruple” (real*16) precision, and perform symbolic calculations in Mathematica.

Keywords

• incompressible fluid
• magnetic mode
• alpha-effect
• eddy diffusivity
• eddy viscosity
• Padé approximant