-
Numerical evidence for a small-scale dynamo approaching solar magnetic Prandtl numbers
Authors:
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Frederick A. Gent,
Matthias Rheinhardt
Abstract:
Magnetic fields on small scales are ubiquitous in the universe. Though they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers ($Pr_M$) such as are present in the Sun and othe…
▽ More
Magnetic fields on small scales are ubiquitous in the universe. Though they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers ($Pr_M$) such as are present in the Sun and other cool stars. We have performed high-resolution simulations of isothermal forced turbulence employing the lowest $Pr_M$ values so far achieved. Contrary to earlier findings, the SSD turns out to be not only possible for $Pr_M$ down to 0.0031, but even becomes increasingly easier to excite for $Pr_M$ below $\simeq\,$0.05. We relate this behaviour to the known hydrodynamic phenomenon referred to as the bottleneck effect. Extrapolating our results to solar values of $Pr_M$ indicates that an SSD would be possible under such conditions.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
Compressible test-field method and its application to shear dynamos
Authors:
Maarit J. Käpylä,
Matthias Rheinhardt,
Axel Brandenburg
Abstract:
In this study we present a compressible test-field method (CTFM) for computing $α$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in…
▽ More
In this study we present a compressible test-field method (CTFM) for computing $α$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in simplified MHD (SMHD). We present several test cases using velocity and magnetic fields of the Roberts geometry, and also compare with the imposed-field method. We show that, for moderate imposed field strengths, the nonlinear CTFM (nCTFM) gives results in agreement with the imposed-field method. Comparison of different flavors of the nCTFM in the shear dynamo case also agree up to equipartition field strengths. Some deviations between the CTFM and SMHD variants exist. As a relevant physical application, we study non-helically forced shear flows, which exhibit large-scale dynamo action, and present a re-analysis of low Reynolds number, moderate shear systems, where we previously neglected the pressure gradient in the momentum equation, and found no coherent shear-current effect. Another key difference is that in the earlier study we used magnetic forcing to mimic small-scale dynamo action, while here it is self-consistently driven by purely kinetic forcing. The kinematic CTFM with general validity forms the core of our analysis. We still find no coherent shear-current effect, but do recover strong large-scale dynamo action that, according to our analysis, is driven through the incoherent effects.
△ Less
Submitted 3 February, 2022; v1 submitted 2 June, 2021;
originally announced June 2021.
-
Investigating global convective dynamos with mean-field models: full spectrum of turbulent effects required
Authors:
Jörn Warnecke,
Matthias Rheinhardt,
Mariangela Viviani,
Frederick Gent,
Simo Tuomisto,
Maarit J. Käpylä
Abstract:
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport…
▽ More
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport coefficients, based on the MF framework, from an exemplary GCD simulation using the test-field method. These coefficients are then used as an input into an MF model. We find a good agreement between the MF and GCD solutions, which validates our theoretical approach. This agreement requires all turbulent effects to be included, even those which have been regarded as unimportant so far. Our results suggest that simple dynamo models, as are commonly used in the solar and stellar community, relying on very few, precisely fine-tuned turbulent effects, may not be representative of the full dynamics of dynamos in global convective simulations and astronomical objects.
△ Less
Submitted 5 October, 2021; v1 submitted 17 May, 2021;
originally announced May 2021.
-
Scalable communication for high-order stencil computations using CUDA-aware MPI
Authors:
Johannes Pekkilä,
Miikka S. Väisälä,
Maarit J. Käpylä,
Matthias Rheinhardt,
Oskar Lappi
Abstract:
Modern compute nodes in high-performance computing provide a tremendous level of parallelism and processing power. However, as arithmetic performance has been observed to increase at a faster rate relative to memory and network bandwidths, optimizing data movement has become critical for achieving strong scaling in many communication-heavy applications. This performance gap has been further accent…
▽ More
Modern compute nodes in high-performance computing provide a tremendous level of parallelism and processing power. However, as arithmetic performance has been observed to increase at a faster rate relative to memory and network bandwidths, optimizing data movement has become critical for achieving strong scaling in many communication-heavy applications. This performance gap has been further accentuated with the introduction of graphics processing units, which can provide by multiple factors higher throughput in data-parallel tasks than central processing units. In this work, we explore the computational aspects of iterative stencil loops and implement a generic communication scheme using CUDA-aware MPI, which we use to accelerate magnetohydrodynamics simulations based on high-order finite differences and third-order Runge-Kutta integration. We put particular focus on improving intra-node locality of workloads. Our GPU implementation scales strongly from one to $64$ devices at $50\%$--$87\%$ of the expected efficiency based on a theoretical performance model. Compared with a multi-core CPU solver, our implementation exhibits $20$--$60\times$ speedup and $9$--$12\times$ improved energy efficiency in compute-bound benchmarks on $16$ nodes.
△ Less
Submitted 10 May, 2022; v1 submitted 2 March, 2021;
originally announced March 2021.
-
Hunting down the cause of solar magnetism
Authors:
M. Viviani,
A. Prabhu,
J. Warnecke,
L. Duarte,
J. Pekkilä,
M. Rheinhardt,
M. J. Käpylä
Abstract:
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions…
▽ More
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Relic gravitational waves from the chiral magnetic effect
Authors:
Axel Brandenburg,
Yutong He,
Tina Kahniashvili,
Matthias Rheinhardt,
Jennifer Schober
Abstract:
Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity…
▽ More
Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity and characteristic wavenumber) and the speed characterizing the maximum magnetic field strength expected from the CME. When the latter exceeds the former (regime I), the regime applicable to the early universe, we obtain an inverse cascade with moderate GW energy that scales with the third power of the magnetic energy. When the generation speed exceeds the CME limit (regime II), the GW energy continues to increase without a corresponding increase of magnetic energy. In the early kinematic phase, the GW energy spectrum (per linear wavenumber interval) has opposite slopes in both regimes and is characterized by an inertial range spectrum in regime I and a white noise spectrum in regime II. The occurrence of these two slopes is shown to be a generic consequence of a nearly monochromatic exponential growth of the magnetic field. The resulting GW energy is found to be proportional to the fifth power of the limiting CME speed and the first power of the generation speed.
△ Less
Submitted 22 April, 2021; v1 submitted 20 January, 2021;
originally announced January 2021.
-
Interaction of large- and small-scale dynamos in isotropic turbulent flows from GPU-accelerated simulations
Authors:
Miikka S. Väisälä,
Johannes Pekkilä,
Maarit J. Käpylä,
Matthias Rheinhardt,
Hsien Shang,
Ruben Krasnopolsky
Abstract:
Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both…
▽ More
Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both helical and non-helical forcing was employed to separate the effects of LSD and SSD in a periodic domain. Reynolds numbers (Rm) up to $\approx 250$ were examined and multiple resolutions used for convergence checks. We ran our simulations with the Astaroth code, designed to accelerate 3D stencil computations on graphics processing units (GPUs) and to employ multiple GPUs with peer-to-peer communication. We observed a speedup of $\approx 35$ in single-node performance compared to the widely used multi-CPU MHD solver Pencil Code. We estimated the growth rates both from the averaged magnetic fields and their power spectra. At low Rm, LSD growth dominates, but at high Rm SSD appears to dominate in both helically and non-helically forced cases. Pure SSD growth rates follow a logarithmic scaling as a function of Rm. Probability density functions of the magnetic field from the growth stage exhibit SSD behaviour in helically forced cases even at intermediate Rm. We estimated mean-field turbulence transport coefficients using closures like the second-order correlation approximation (SOCA). They yield growth rates similar to the directly measured ones and provide evidence of $α$ quenching. Our results are consistent with the SSD inhibiting the growth of the LSD at moderate Rm, while the dynamo growth is enhanced at higher Rm.
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
Authors:
A. Brandenburg,
A. Johansen,
P. A. Bourdin,
W. Dobler,
W. Lyra,
M. Rheinhardt,
S. Bingert,
N. E. L. Haugen,
A. Mee,
F. Gent,
N. Babkovskaia,
C. -C. Yang,
T. Heinemann,
B. Dintrans,
D. Mitra,
S. Candelaresi,
J. Warnecke,
P. J. Käpylä,
A. Schreiber,
P. Chatterjee,
M. J. Käpylä,
X. -Y. Li,
J. Krüger,
J. R. Aarnes,
G. R. Sarson
, et al. (12 additional authors not shown)
Abstract:
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the…
▽ More
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
On the existence of shear-current effects in magnetized burgulence
Authors:
Maarit J. Käpylä,
Javier Álvarez Vizoso,
Matthias Rheinhardt,
Axel Brandenburg,
Nishant K. Singh
Abstract:
The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbolη$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficie…
▽ More
The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbolη$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficients from systems with magnetic background turbulence, we present an extension of the test-field method (TFM), applicable to our setup where the pressure gradient is dropped from the momentum equation: the nonlinear TFM (NLTFM). Our momentum equation is related to Burgers' equation and the resulting flows are referred to as magnetized burgulence. We use both stochastic kinetic and magnetic forcings to mimic cases without and with simultaneous small-scale dynamo action (SSD). When we force only kinetically, negative $η_{yx}$ are obtained with exponential growth in both the radial and azimuthal mean magnetic field components. Using isotropic magnetic forcing, the field growth is no longer exponential, while NLTFM yields positive $η_{yx}$. By employing an alternative forcing from which wavevectors having small components are removed, the exponential growth is recovered, but the NLTFM results do not change significantly. Analyzing the dynamo excitation conditions for the coherent SC and incoherent $α$ and SC effects shows that the incoherent effects are the main drivers of the dynamo in the majority of cases. We find no evidence for MSC-effect-driven dynamos in our simulations.
△ Less
Submitted 12 January, 2021; v1 submitted 10 June, 2020;
originally announced June 2020.
-
Stellar dynamos in the transition regime: multiple dynamo modes and anti-solar differential rotation
Authors:
M. Viviani,
M. J. Käpylä,
J. Warnecke,
P. J. Käpylä,
M. Rheinhardt
Abstract:
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions…
▽ More
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $αΩ$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $α$ effect is comparable or even larger than the $Ω$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $α^2Ω$ or $α^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping.
△ Less
Submitted 22 November, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence
Authors:
Petri J. Käpylä,
Matthias Rheinhardt,
Axel Brandenburg,
Maarit J. Käpylä
Abstract:
(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbule…
▽ More
(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions.
△ Less
Submitted 17 February, 2020; v1 submitted 3 January, 2019;
originally announced January 2019.
-
$f$-mode strengthening from a localized bipolar subsurface magnetic field
Authors:
Nishant K. Singh,
Harsha Raichur,
Maarit J. Käpylä,
Matthias Rheinhardt,
Axel Brandenburg,
Petri J. Käpylä
Abstract:
Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. So…
▽ More
Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the $f$-mode is significantly larger at high horizontal wavenumbers $k$, but the fanning is weaker for the localized subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the $f$-mode in the $kω$ diagram increases proportionally in such a way that the normalized $f$-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localized magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the $f$-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of $f$-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots.
△ Less
Submitted 12 April, 2019; v1 submitted 27 August, 2018;
originally announced August 2018.
-
Extended subadiabatic layer in simulations of overshooting convection
Authors:
Petri J. Käpylä,
Matthias Rheinhardt,
Axel Brandenburg,
Rainer Arlt,
Maarit J. Käpylä,
Andreas Lagg,
Nigul Olspert,
Jörn Warnecke
Abstract:
We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kr…
▽ More
We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.
△ Less
Submitted 23 August, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Turbulent transport coefficients in spherical wedge dynamo simulations of solar-like stars
Authors:
Jörn Warnecke,
Matthias Rheinhardt,
S. Tuomisto,
Petri J. Käpylä,
Maarit J. Käpylä,
Axel Brandenburg
Abstract:
We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine t…
▽ More
We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $α$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $φφ$-component of the $α$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case.
△ Less
Submitted 16 January, 2018; v1 submitted 14 January, 2016;
originally announced January 2016.
-
Fanning out of the $f$-mode in presence of nonuniform magnetic fields
Authors:
Nishant K. Singh,
Axel Brandenburg,
Matthias Rheinhardt
Abstract:
We show that in the presence of a harmonically varying magnetic field the fundamental or $f$-mode in a stratified layer is altered in such a way that it fans out in the diagnostic $kω$ diagram, but with mode power also within the fan. In our simulations, the surface is defined by a temperature and density jump in a piecewise isothermal layer. Unlike our previous work (Singh et al. 2014) where a un…
▽ More
We show that in the presence of a harmonically varying magnetic field the fundamental or $f$-mode in a stratified layer is altered in such a way that it fans out in the diagnostic $kω$ diagram, but with mode power also within the fan. In our simulations, the surface is defined by a temperature and density jump in a piecewise isothermal layer. Unlike our previous work (Singh et al. 2014) where a uniform magnetic field was considered, we employ here a nonuniform magnetic field together with hydromagnetic turbulence at length scales much smaller than those of the magnetic fields. The expansion of the $f$-mode is stronger for fields confined to the layer below the surface. In some of those cases, the $kω$ diagram also reveals a new class of low frequency vertical stripes at multiples of twice the horizontal wavenumber of the background magnetic field. We argue that the study of the $f$-mode expansion might be a new and sensitive tool to determining subsurface magnetic fields with longitudinal periodicity.
△ Less
Submitted 1 July, 2014;
originally announced July 2014.
-
Quenching and anisotropy of hydromagnetic turbulent transport
Authors:
Bidya Binay Karak,
Matthias Rheinhardt,
Axel Brandenburg,
Petri J. Kapyla,
Maarit J. Kapyla
Abstract:
Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $α$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Ro…
▽ More
Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $α$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the {\it quenched} flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows these quenching exponents become about 4, 1.5, and 2.3, respectively. For the $α$ effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the $\bf Ω\times \bf J$ effect nearly the same quenching exponent is obtained as for $α$. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents.
△ Less
Submitted 20 August, 2014; v1 submitted 17 June, 2014;
originally announced June 2014.
-
Properties of $p$- and $f$-modes in hydromagnetic turbulence
Authors:
Nishant K. Singh,
Axel Brandenburg,
S. M. Chitre,
Matthias Rheinhardt
Abstract:
With the ultimate aim of using the fundamental or $f$-mode to study helioseismic aspects of turbulence-generated magnetic flux concentrations, we use randomly forced hydromagnetic simulations of a piecewise isothermal layer in two dimensions with reflecting boundaries at top and bottom. We compute numerically diagnostic wavenumber-frequency diagrams of the vertical velocity at the interface betwee…
▽ More
With the ultimate aim of using the fundamental or $f$-mode to study helioseismic aspects of turbulence-generated magnetic flux concentrations, we use randomly forced hydromagnetic simulations of a piecewise isothermal layer in two dimensions with reflecting boundaries at top and bottom. We compute numerically diagnostic wavenumber-frequency diagrams of the vertical velocity at the interface between the denser gas below and the less dense gas above. For an Alfvén-to-sound speed ratio of about 0.1, a 5% frequency increase of the $f$-mode can be measured when $k_xH_{\rm p}=3$-$4$, where $k_x$ is the horizontal wavenumber and $H_{\rm p}$ is the pressure scale height at the surface. Since the solar radius is about 2000 times larger than $H_{\rm p}$, the corresponding spherical harmonic degree would be 6000-8000. For weaker fields, a $k_x$-dependent frequency decrease by the turbulent motions becomes dominant. For vertical magnetic fields, the frequency is enhanced for $k_xH_{\rm p}\approx4$, but decreased relative to its nonmagnetic value for $k_xH_{\rm p}\approx9$.
△ Less
Submitted 21 March, 2015; v1 submitted 11 April, 2014;
originally announced April 2014.
-
Mean-field dynamo action from delayed transport
Authors:
Matthias Rheinhardt,
Ebru Devlen,
Karl-Heinz Rädler,
Axel Brandenburg
Abstract:
We analyze the nature of dynamo action that produces horizontally averaged magnetic fields in two particular flows that were studied by Roberts (1972, Phil. Trans. R. Soc. A 271, 411), namely his flows II and III. They have zero kinetic helicity either pointwise (flow II), or on average (flow III). Using direct numerical simulations, we determine the onset conditions for dynamo action at moderate…
▽ More
We analyze the nature of dynamo action that produces horizontally averaged magnetic fields in two particular flows that were studied by Roberts (1972, Phil. Trans. R. Soc. A 271, 411), namely his flows II and III. They have zero kinetic helicity either pointwise (flow II), or on average (flow III). Using direct numerical simulations, we determine the onset conditions for dynamo action at moderate values of the magnetic Reynolds number. Using the test-field method, we show that the turbulent magnetic diffusivity is then positive for both flows. However, we demonstrate that for both flows large-scale dynamo action occurs through delayed transport. Mathematically speaking, the magnetic field at earlier times contributes to the electromotive force through the off-diagonal components of the alpha tensor such that a zero mean magnetic field becomes unstable to dynamo action. This represents a qualitatively new mean-field dynamo mechanism not previously described.
△ Less
Submitted 20 January, 2014;
originally announced January 2014.
-
Testing turbulent closure models with convection simulations
Authors:
J. E. Snellman,
P. J. Käpylä,
M. J. Käpylä,
M. Rheinhardt,
B. Dintrans
Abstract:
We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors,…
▽ More
We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors, corresponding to locating the computational domain at different latitudes ranging from the pole to the equator of the star. We perform three-dimensional numerical simulations in the same parameter regimes for comparison. The free parameters appearing in the closure models are calibrated by two fitting methods using simulation data. Unique determination of the closure parameters is possible only in the non-rotating case or when the system is placed at the pole. In the other cases the fit procedures yield somewhat differing results. The quality of the closure is tested by substituting the resulting coefficients back into the closure model and comparing with the simulation results. To eliminate the possibilities that the results obtained depend on the aspect ratio of the simulation domain or suffer from too small Rayleigh numbers we performed runs varying these parameters. The simulation data for the Reynolds stress and heat fluxes broadly agree with previous compressible simulations. The closure works fairly well with slow and fast rotation but its quality degrades for intermediate rotation rates. We find that the closure parameters depend not only on rotation rate but also on latitude. The weak dependence on Rayleigh number and the aspect ratio of the domain indicates that our results are generally valid
△ Less
Submitted 24 September, 2014; v1 submitted 21 September, 2012;
originally announced September 2012.
-
Mean-field closure parameters for passive scalar turbulence
Authors:
J. E. Snellman,
M. Rheinhardt,
P. J. Käpylä,
M. J. Mantere,
A. Brandenburg
Abstract:
Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the qua…
▽ More
Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution.
We vary the Reynolds (Re) and Péclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.
△ Less
Submitted 20 December, 2011;
originally announced December 2011.
-
Magnetic feature tracking, what determines the speed?
Authors:
G. Guerrero,
M. Rheinhardt,
M. Dikpati
Abstract:
Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_θ$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_θ|$. As the gradient's sign (inferred by local helioseismology) is sti…
▽ More
Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_θ$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_θ|$. As the gradient's sign (inferred by local helioseismology) is still unclear, cases with a negative gradient are studied in this paper. Under this condition, the velocity of the magnetic tracers turns out to be larger than the plasma velocity, in disagreement with the observations. Alternative mechanisms for explaining them independently are proposed. For the turbulent magnetic pumping it is shown that it has to be unrealistically strong to reconcile the model with the observations.
△ Less
Submitted 15 November, 2011;
originally announced November 2011.
-
Modeling spatio-temporal nonlocality in mean-field dynamos
Authors:
M. Rheinhardt,
A. Brandenburg
Abstract:
When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving…
▽ More
When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean-field equations are solved for oscillatory alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.
△ Less
Submitted 13 October, 2011;
originally announced October 2011.
-
Plasma flow vs. magnetic feature-tracking speeds in the Sun
Authors:
G. Guerrero,
M. Rheinhardt,
A. Brandenburg,
M. Dikpati
Abstract:
We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensiona…
▽ More
We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.
△ Less
Submitted 24 July, 2011;
originally announced July 2011.
-
Mean-field diffusivities in passive scalar and magnetic transport in irrotational flows
Authors:
Karl-Heinz Rädler,
Axel Brandenburg,
Fabio Del Sordo,
Matthias Rheinhardt
Abstract:
Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding…
▽ More
Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding the magnetic mean-field diffusivity in an electrically conducting turbulently moving compressible fluid. These phenomena occur if the irrotational part of the motion dominates the vortical part, the Pèclet or magnetic Reynolds number is not too large, and, in addition, the variation of the flow pattern is slow. For both the passive scalar and the magnetic cases several further analytical results on mean-field diffusivities and related quantities found within the second-order correlation approximation are presented, as well as numerical results obtained by the test-field method, which applies independently of this approximation. Particular attention is paid to non-local and non-instantaneous connections between the turbulence-caused terms and the mean fields. Two examples of irrotational flows, in which interesting phenomena in the above sense occur, are investigated in detail. In particular, it is demonstrated that the decay of a mean scalar in a compressible fluid under the influence of these flows can be much slower than without any flow, and can be strongly influenced by the so-called memory effect, that is, the fact that the relevant mean-field coefficients depend on the decay rates themselves.
△ Less
Submitted 8 November, 2011; v1 submitted 8 April, 2011;
originally announced April 2011.
-
The fratricide of alpha-Omega dynamos by their alpha-squared siblings
Authors:
Alexander Hubbard,
Matthias Rheinhardt,
Axel Brandenburg
Abstract:
Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the α - Ω mode is generally expected to dominate because it is the fastest growing. In an α - Ω dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be t…
▽ More
Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the α - Ω mode is generally expected to dominate because it is the fastest growing. In an α - Ω dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory α - Ω type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both α - Ω and α2 dynamo action is possible. Results. After initial dominance and saturation, the α - Ω mode is found to be destroyed by an orthogonal α2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of α - Ω dynamos by its α2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of α - Ω dynamos from an α2 dynamo has not yet been found.
△ Less
Submitted 13 February, 2011;
originally announced February 2011.
-
Spontaneous chiral parity breaking by hydromagnetic buoyancy
Authors:
Piyali Chatterjee,
Dhrubaditya Mitra,
Axel Brandenburg,
Matthias Rheinhardt
Abstract:
Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and…
▽ More
Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.
△ Less
Submitted 24 August, 2011; v1 submitted 4 November, 2010;
originally announced November 2010.
-
Alpha effect due to buoyancy instability of a magnetic layer
Authors:
Piyali Chatterjee,
Dhrubaditya Mitra,
Matthias Rheinhardt,
Axel Brandenburg
Abstract:
A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with th…
▽ More
A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.
△ Less
Submitted 4 November, 2010;
originally announced November 2010.
-
Turbulent transport in hydromagnetic flows
Authors:
A. Brandenburg,
P. Chatterjee,
F. Del Sordo,
A. Hubbard,
P. J. Käpylä,
M. Rheinhardt
Abstract:
The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha ef…
▽ More
The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.
△ Less
Submitted 29 April, 2010;
originally announced April 2010.
-
Test-field method for mean-field coefficients with MHD background
Authors:
M. Rheinhardt,
A. Brandenburg
Abstract:
Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using both the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describi…
▽ More
Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using both the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both forces are chosen to have Roberts flow-like geometry. Results: Examples with an MHD background are studied where the previously used quasi-kinematic test-field method breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the same results as the imposed-field method, where the field-aligned component of the actual electromotive force from the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, $α$ is found to be quenched proportional to the fourth power of the field strength, regardless of the type of background studied.
△ Less
Submitted 5 April, 2010;
originally announced April 2010.
-
Nonlinear magnetic diffusivity and alpha tensors in helical turbulence
Authors:
A. Brandenburg,
K. -H. Rädler,
M. Rheinhardt,
K. Subramanian
Abstract:
The effect of a dynamo-generated mean magnetic field of Beltrami type on the mean electromotive force is studied. In the absence of the mean magnetic field the turbulence is assumed to be homogeneous and isotropic, but it becomes inhomogeneous and anisotropic with this field. Using the testfield method the dependence of the alpha and turbulent diffusivity tensors on the magnetic Reynolds number…
▽ More
The effect of a dynamo-generated mean magnetic field of Beltrami type on the mean electromotive force is studied. In the absence of the mean magnetic field the turbulence is assumed to be homogeneous and isotropic, but it becomes inhomogeneous and anisotropic with this field. Using the testfield method the dependence of the alpha and turbulent diffusivity tensors on the magnetic Reynolds number Rm is determined for magnetic fields that have reached approximate equipartition with the velocity field. The tensor components are characterized by a pseudoscalar alpha and a scalar turbulent magnetic diffusivity etat. Increasing Rm from 2 to 600 reduces etat by a factor ~5, suggesting that the quenching of etat is, in contrast to the 2-dimensional case, only weakly dependent on Rm. Over the same range of Rm, however, alpha is reduced by a factor ~14, which can qualitatively be explained by a corresponding increase of a magnetic contribution to the alpha effect with opposite sign. The level of fluctuations of alpha and etat is only 10% and 20% of the respective kinematic reference values.
△ Less
Submitted 14 October, 2008; v1 submitted 9 May, 2008;
originally announced May 2008.
-
Magnetic diffusivity tensor and dynamo effects in rotating and shearing turbulence
Authors:
A. Brandenburg,
K. -H. Rädler,
M. Rheinhardt,
P. J. Käpylä
Abstract:
The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from…
▽ More
The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and opposite to the sign of the shear. This implies that dynamo action from the shear--current effect is impossible, except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the average, the components of the alpha tensor display Gaussian fluctuations around zero. These fluctuations are strong enough to drive an incoherent alpha--shear dynamo. The incoherent shear--current effect, on the other hand, is found to be subdominant.
△ Less
Submitted 28 February, 2008; v1 submitted 22 October, 2007;
originally announced October 2007.
-
Mean-field concept and direct numerical simulations of rotating magnetoconvection and the geodynamo
Authors:
M. Schrinner,
K. -H. Raedler,
D. Schmitt,
M. Rheinhardt,
U. R. Christensen
Abstract:
A comparison is made between mean-field models and direct numerical simulations of rotating magnetoconvection and the geodynamo. The mean-field coefficients are calculated with the fluid velocity taken from the direct numerical simulations. The magnetic fields resulting from mean-field models are then compared with the mean magnetic field from the direct numerical simulations.
A comparison is made between mean-field models and direct numerical simulations of rotating magnetoconvection and the geodynamo. The mean-field coefficients are calculated with the fluid velocity taken from the direct numerical simulations. The magnetic fields resulting from mean-field models are then compared with the mean magnetic field from the direct numerical simulations.
△ Less
Submitted 27 September, 2006;
originally announced September 2006.
-
Mean--field electrodynamics: Critical analysis of various analytical approaches to the mean electromotive force
Authors:
Karl-Heinz Raedler,
Matthias Rheinhardt
Abstract:
There are various analytical approaches to the mean electromotive force $\cal E =< \vec{u}\times \vec{b}>$ crucial in mean--field electrodynamics, with $\vec{u}$ and
$\vec{b}$ being velocity and magnetic field fluctuations. In most cases the traditional approach, restricted to the second--order correlation approximation, has been used. Its validity is only guaranteed for a range of conditions,…
▽ More
There are various analytical approaches to the mean electromotive force $\cal E =< \vec{u}\times \vec{b}>$ crucial in mean--field electrodynamics, with $\vec{u}$ and
$\vec{b}$ being velocity and magnetic field fluctuations. In most cases the traditional approach, restricted to the second--order correlation approximation, has been used. Its validity is only guaranteed for a range of conditions, which is narrow in view of many applications, e.g., in astrophysics. With the intention to have a wider range of applicability other approaches have been proposed which make use of the so--called $τ$--approximation, reducing correlations of third order in $\vec{u}$ and $\vec{b}$ to such of second order. After explaining some basic features of the traditional approach a critical analysis of the approaches of that kind is given. It is shown that they lead in some cases to results which are in clear conflict with those of the traditional approach. It is argued that this indicates shortcomings of the $τ$--approaches and poses serious restrictions to their applicability. These shortcomings do not result from the basic assumption of the $τ$--approximation. Instead, they seem to originate in some simplifications made in order to derive $\cal E$ without really solving the equations governing $\vec{u}$ and $\vec{b}$. A starting point for a new approach is described which avoids the conflict.
△ Less
Submitted 22 November, 2006; v1 submitted 12 June, 2006;
originally announced June 2006.
-
Magnetars versus Radio Pulsars: MHD Stability in Newborn Highly Magnetized Neutron Stars
Authors:
U. Geppert,
M. Rheinhardt
Abstract:
We study the stability/establishment of dipolar magnetostatic equilibrium configurations in new--born neutron stars (NSs) in dependence on the rotational velocity $Ω$ and on the initial angle $α$ between rotation and magnetic axis. The NS is modeled as a sphere of a highly magnetized ($B \sim 10^{15}$G) incompressible fluid of uniform density which rotates rigidly. For the initial dipolar backgr…
▽ More
We study the stability/establishment of dipolar magnetostatic equilibrium configurations in new--born neutron stars (NSs) in dependence on the rotational velocity $Ω$ and on the initial angle $α$ between rotation and magnetic axis. The NS is modeled as a sphere of a highly magnetized ($B \sim 10^{15}$G) incompressible fluid of uniform density which rotates rigidly. For the initial dipolar background magnetic field, which defines the magnetic axis, two different configurations are assumed. We solve the 3D non--linear MHD equations by use of a spectral code. The problem in dimensionless form is completely defined by the initial field strength (for a fixed field geometry), the magnetic Prandtl number $\Pm$, and the normalized rotation rate. The evolution of the magnetic and velocity fields is considered for initial magnetic field strengths characterized by the ratio of ohmic diffusion and initial \Alf{} travel times $\ttOhm/\ttAO \approx 1000$, for $\Pm = 0.1, 1, 10$, and the ratio of rotation period and initial \Alf{} travel time, $P/\ttAO = 0.012, 0.12, 1.2, 12$. We find hints for the existence of a unique stable dipolar magnetostatic configuration for any specific $α$, independent of the initial field geometry. Comparing NSs possessing the same field structure at the end of their proto--NS phase, it turns out that sufficiently fast rotating NSs ($P\la6 $ms) with $α\la 45^0$ retain their magnetar field, while the others lose almost all of their initial magnetic energy by transferring it into magnetic and kinetic energy of relatively small--scaled fields and continue their life as radio pulsars with a dipolar surface field of $10^{12...13}$G.
△ Less
Submitted 6 June, 2006;
originally announced June 2006.
-
The proto--neutron--star dynamo -- viability and impediments
Authors:
M. Rheinhardt,
U. Geppert
Abstract:
We study convective motions taken from hydrodynamic simulations of rotating proto--neutron stars (PNSs) with respect to their ability to excite a dynamo instability which may be responsible for the giant neutron star magnetic fields. Since it is impossible to simulate the magnetic field evolution employing the actual magnetic Reynolds numbers ($\Rm$) resulting from the hydrodynamic simulations,…
▽ More
We study convective motions taken from hydrodynamic simulations of rotating proto--neutron stars (PNSs) with respect to their ability to excite a dynamo instability which may be responsible for the giant neutron star magnetic fields. Since it is impossible to simulate the magnetic field evolution employing the actual magnetic Reynolds numbers ($\Rm$) resulting from the hydrodynamic simulations, (smallest) critical $\Rm$s and the corresponding field geometries are derived on the kinematic level by rescaling the velocity amplitudes. It turns out that the actual values of $\Rm$ are by many orders of magnitude larger than the critical values found. A dynamo might therefore start to act vigorously very soon after the onset of convection. But as in general dynamo growth rates are non--monotonous functions of $\Rm$ the later fate of the magnetic field is uncertain. Hence, no reliable statements on the existence and efficiency of PNS dynamos can be drawn without considering the interplay of magnetic field and convection from the beginning. Likewise, in so far as convection inside the PNS is regarded to be essential in re--launching the supernova explosion, a revision of its role in this respect could turn out to be necessary.
△ Less
Submitted 14 February, 2005;
originally announced February 2005.
-
Comment on ``The linear instability of magnetic Taylor-Couette flow with Hall effect''
Authors:
M. Rheinhardt,
U. Geppert
Abstract:
In the paper we comment on (Rüdiger & Shalybkov, Phys. Rev. E. 69, 016303 (2004) (RS)), the instability of the Taylor--Couette flow interacting with a homogeneous background field subject to Hall effect is studied. We correct a falsely generalizing interpretation of results presented there which could be taken to disprove the existence of the Hall--drift induced magnetic instability described in…
▽ More
In the paper we comment on (Rüdiger & Shalybkov, Phys. Rev. E. 69, 016303 (2004) (RS)), the instability of the Taylor--Couette flow interacting with a homogeneous background field subject to Hall effect is studied. We correct a falsely generalizing interpretation of results presented there which could be taken to disprove the existence of the Hall--drift induced magnetic instability described in Rheinhardt and Geppert, Phys. Rev. Lett. 88, 101103. It is shown that in contrast to what is suggested by RS, no additional shear flow is necessary to enable such an instability with a non--potential magnetic background field, whereas for a curl--free one it is. In the latter case, the instabilities found in RS in situations where neither a hydrodynamic nor a magneto--rotational instability exists are demonstrated to be most likely magnetic instead of magnetohydrodynamic. Further, some minor inaccuracies are clarified.
△ Less
Submitted 17 January, 2005;
originally announced January 2005.
-
Spot-like Structures of Neutron Star Surface Magnetic Fields
Authors:
U. Geppert,
M. Rheinhardt,
J. Gil
Abstract:
There is growing evidence, based on both X-ray and radio observations of isolated neutron stars, that besides the large--scale (dipolar) magnetic field, which determines the pulsar spin--down behaviour, small--scale poloidal field components are present, which have surface strengths one to two orders of magnitude larger than the dipolar component. We argue in this paper that the Hall--effect can…
▽ More
There is growing evidence, based on both X-ray and radio observations of isolated neutron stars, that besides the large--scale (dipolar) magnetic field, which determines the pulsar spin--down behaviour, small--scale poloidal field components are present, which have surface strengths one to two orders of magnitude larger than the dipolar component. We argue in this paper that the Hall--effect can be an efficient process in producing such small--scale field structures just above the neutron star surface. It is shown that due to a Hall--drift induced instability, poloidal magnetic field structures can be generated from strong subsurface toroidal fields, which are the result of either a dynamo or a thermoelectric instability acting at early times of a neutron star's life. The geometrical structure of these small--scale surface anomalies of the magnetic field resembles that of some types of ``star--spots''. The magnetic field strength and the length--scales are comparable with values that can be derived from various observations.
△ Less
Submitted 17 November, 2003; v1 submitted 5 November, 2003;
originally announced November 2003.
-
The Occurrence of the Hall--Instability in Crusts of Isolated Neutron Stars
Authors:
M. Rheinhardt,
D. Konenkov,
U. Geppert
Abstract:
In former papers we showed that during the decay of a neutron star's magnetic field under the influence of the Hall--drift, an unstable rise of small--scale field structures at the expense of the large--scale background field may happen. This linear stability analysis was based on the assumption of a uniform density throughout the neutron star crust, whereas in reality the density and all transp…
▽ More
In former papers we showed that during the decay of a neutron star's magnetic field under the influence of the Hall--drift, an unstable rise of small--scale field structures at the expense of the large--scale background field may happen. This linear stability analysis was based on the assumption of a uniform density throughout the neutron star crust, whereas in reality the density and all transport coefficients vary by many orders of magnitude. Here, we extend the investigation of the Hall--drift induced instability by considering realistic profiles of density and chemical composition, as well as background fields with more justified radial profiles. Two neutron star models are considered differing primarily in the assumption on the core matter equation of state. For their cooling history and radial profiles of density and composition we use known results to infer the conductivity profiles. These were fed into linear calculations of a dipolar field decay starting from various initial configurations. At different stages of the decay, snapshots of the magnetic fields at the equator were taken to yield background field profiles for the stability analysis. The main result is that the Hall instability may really occur in neutron star crusts. Characteristic growth times are in the order of \lesssim 10^4 ... 10^6 yrs depending on cooling age and background field strength. The influence of the equation of state and of the initial field configuration is discussed.
△ Less
Submitted 2 January, 2004; v1 submitted 31 October, 2003;
originally announced October 2003.
-
Non--linear magnetic field decay in neutron stars -- Theory and observations
Authors:
U. Geppert,
M. Rheinhardt
Abstract:
There exists both theoretical and observational evidence that the magnetic field decay in neutron stars may proceed in a pronounced non--linear way during a certain episode of the neutron star's life. In the presence of a strong magnetic field the Hall--drift dominates the field evolution in the crust and/or the superfluid core of neutron stars. Analysing observations of $P$ and $\dot{P}$ for su…
▽ More
There exists both theoretical and observational evidence that the magnetic field decay in neutron stars may proceed in a pronounced non--linear way during a certain episode of the neutron star's life. In the presence of a strong magnetic field the Hall--drift dominates the field evolution in the crust and/or the superfluid core of neutron stars. Analysing observations of $P$ and $\dot{P}$ for sufficiently old isolated pulsars we gain strong hints for a significantly non--linear magnetic field decay. Under certain conditions with respect to the geometry and strength of a large--scale magnetic background field an instability is shown to occur which rapidly raises small--scale magnetic field modes. Their growth rates increase with the background field strength and may reach $\sim 10^4$ times the ohmic decay rate. Consequences for the rotational and thermal evolution as well as for the cracking of the crust of neutron stars are discussed.
△ Less
Submitted 2 July, 2002;
originally announced July 2002.
-
A Hall--Drift Induced Magnetic Field Instability
Authors:
M. Rheinhardt,
U. Geppert
Abstract:
In the presence of a strong magnetic field and under conditions as realized in the crust and the superfluid core of neutron stars the Hall--drift dominates the field evolution. We show by a linear analysis that for a sufficiently strong large--scale background field depending at least quadratically on position in a plane conducting slab an instability occurs which rapidly generates small--scale…
▽ More
In the presence of a strong magnetic field and under conditions as realized in the crust and the superfluid core of neutron stars the Hall--drift dominates the field evolution. We show by a linear analysis that for a sufficiently strong large--scale background field depending at least quadratically on position in a plane conducting slab an instability occurs which rapidly generates small--scale fields. Their growth rates depend on the choice of the boundary conditions, increase with the background field strength and may reach $10^3$ times the ohmic decay rate. The effect of that instability on the rotational and thermal evolution of neutron stars is discussed.
△ Less
Submitted 12 February, 2002;
originally announced February 2002.