Find in Library
Search millions of books, articles, and more
Indexed Open Access Databases
Computation of Kinematic and Magnetic <i>α</i>-Effect and Eddy Diffusivity Tensors by Padé Approximation
oleh: Sílvio M.A. Gama, Roman Chertovskih, Vladislav Zheligovsky
| Format: | Article |
|---|---|
| Diterbitkan: | MDPI AG 2019-06-01 |
Deskripsi
We present examples of Padé approximations of the <inline-formula> <math display="inline"> <semantics> <mi>α</mi> </semantics> </math> </inline-formula>-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 <inline-formula> <math display="inline"> <semantics> <mi>α</mi> </semantics> </math> </inline-formula>-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 <inline-formula> <math display="inline"> <semantics> <mrow> <mn>1</mn> <mo>/</mo> <mi>η</mi> </mrow> </semantics> </math> </inline-formula> around <inline-formula> <math display="inline"> <semantics> <mrow> <mn>1</mn> <mo>/</mo> <mi>η</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics> </math> </inline-formula>. This yields the values of the dominant growth rate to satisfactory accuracy for <inline-formula> <math display="inline"> <semantics> <mi>η</mi> </semantics> </math> </inline-formula>, 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.