Source author record

Miikka S. Väisälä

Miikka S. Väisälä appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR physics.comp-ph physics.flu-dyn astro-ph.CO Distributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

4works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Field

Misalignment between rotation and magnetic field has been suggested to be one type of physical mechanisms which can easen the effects of magnetic braking during collapse of cloud cores leading to formation of protostellar disks. However, its essential factors are poorly understood. Therefore, we perform a more detailed analysis of the physics involved. We analyze existing simulation data to measure the system torques, mass accretion rates and Toomre Q parameters. We also examine the presence of shocks in the system. While advective torques are generally the strongest, we find that magnetic and gravitational torques can play substantial roles in how angular momentum is transferred during the disk formation process. Magnetic torques can shape the accretion flows, creating two-armed magnetized inflow spirals aligned with the magnetic field. We find evidence of an accretion shock that is aligned according to the spiral structure of the system. Inclusion of ambipolar diffusion as explored in this work has shown a slight influence in the small scale structures but not in the main morphology. We discuss potential candidate systems where some of these phenomena could be present.

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.

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.

preprint2009arXiv

Influence of Ohmic diffusion on the excitation and dynamics of MRI

In this paper we make an effort to understand the interaction of turbulence generated by the magnetorotational instability (MRI) with turbulence from other sources, such as supernova explosions (SNe) in galactic disks. First we perform a linear stability analysis (LSA) of non-ideal MRI to derive the limiting value of Ohmic diffusion that is needed to inhibit the growth of the instability for different types of rotation laws. With the help of a simple analytical expression derived under first-order smoothing approximation (FOSA), an estimate of the limiting turbulence level and hence the turbulent diffusion needed to damp the MRI is derived. Secondly, we perform numerical simulations in local cubes of isothermal nonstratified gas with external forcing of varying strength to see whether the linear result holds for more complex systems. Purely hydrodynamic calculations with forcing, rotation and shear are made for reference purposes, and as expected, non-zero Reynolds stresses are found. In the magnetohydrodynamic calculations, therefore, the total stresses generated are a sum of the forcing and MRI contributions. To separate these contributions, we perform reference runs with MRI-stable shear profiles (angular velocity increasing outwards), which suggest that the MRI-generated stresses indeed become strongly suppressed as function of the forcing. The Maxwell to Reynolds stress ratio is observed to decrease by an order of magnitude as the turbulence level due to external forcing exceeds the predicted limiting value, which we interpret as a sign of MRI suppression. Finally, we apply these results to estimate the limiting radius inside of which the SN activity can suppress the MRI, arriving at a value of 14 kpc.