BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthias Rheinhardt

Matthias Rheinhardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRphysics.flu-dynphysics.plasm-phastro-ph.GAphysics.comp-phastro-phastro-ph.EPastro-ph.HEastro-ph.IMcond-mat.stat-mechDistributed, Parallel, and Cluster Computingnlin.AOphysics.ao-ph

Catalog footprint

What is connected

14works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Compressible test-field method and its application to shear dynamos

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scalable communication for high-order stencil computations using CUDA-aware MPI

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the existence of shear-current effects in magnetized burgulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interaction of large- and small-scale dynamos in isotropic turbulent flows from GPU-accelerated simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence

(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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Properties of $p$- and $f$-modes in hydromagnetic turbulence

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fanning out of the $f$-mode in presence of nonuniform magnetic fields

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mean-field dynamo action from delayed transport

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quenching and anisotropy of hydromagnetic turbulent transport

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mean-field diffusivities in passive scalar and magnetic transport in irrotational flows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spontaneous chiral parity breaking by hydromagnetic buoyancy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The fratricide of alpha-Omega dynamos by their alpha-squared siblings

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Alpha effect due to buoyancy instability of a magnetic layer

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Mean--field electrodynamics: Critical analysis of various analytical approaches to the mean electromotive force

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.

0 citations0 reviewsNo code linkedOpen access