No Arabic abstract
Coherent magnetic fields in disc galaxies are thought to be generated by a large-scale (or mean-field) dynamo operating in their interstellar medium. A key driver of mean magnetic field growth is the turbulent electromotive force (EMF), which represents the influence of correlated small-scale (or fluctuating) velocity and magnetic fields on the mean field. The EMF is usually expressed as a linear expansion in the mean magnetic field and its derivatives, with the dynamo tensors as expansion coefficients. Here, we adopt the singular value decomposition (SVD) method to directly measure these turbulent transport coefficients in a simulation of the turbulent interstellar medium that realizes a large-scale dynamo. Specifically, the SVD is used to least-square fit the time series data of the EMF with that of the mean field and its derivatives, to determine these coefficients. We demonstrate that the spatial profiles of the EMF reconstructed from the SVD coefficients match well with that taken directly from the simulation. Also, as a direct test, we use the coefficients to simulate a 1-D mean-field dynamo model and find an overall similarity in the evolution of the mean magnetic field between the dynamo model and the direct simulation. We also compare the results with those which arise using simple regression and the ones obtained previously using the test-field (TF) method, to find reasonable qualitative agreement. Overall, the SVD method provides an effective post-processing tool to determine turbulent transport coefficients from simulations.
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 ($phi$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $alpha$ 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 $phiphi$-component of the $alpha$ 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.
The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dipolar and quadrupolar modes. We compute the turbulent transport coefficients, describing the effects of the turbulent velocity field on the mean magnetic field, using the test-field method. The test-field analysis provides a plausible explanation of the missing cycle in terms of the reduction of $alpha_{phiphi}$ in advance of the reduced surface activity, and enhanced downward turbulent pumping during the event to confine some of the magnetic field at the bottom of the convection zone, where local maximum of magnetic energy is observed during the event. At the same time, however, a quenching of the turbulent magnetic diffusivities is observed, albeit differently distributed in depth compared to the other transport coefficients. Therefore, dedicated mean-field modelling is required for verification.
The ordered magnetic field observed via polarized synchrotron emission in nearby disc galaxies can be explained by a mean-field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical-flux initial conditions are potentially able to influence this dynamo via the occurrence of the magneto-rotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean-field dynamo. This is motivated by the observation that different saturation behavior was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three-dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing-periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean-field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched $alphaOmega$~dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxys star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub-grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields.
This paper introduces the functional tensor singular value decomposition (FTSVD), a novel dimension reduction framework for tensors with one functional mode and several tabular modes. The problem is motivated by high-order longitudinal data analysis. Our model assumes the observed data to be a random realization of an approximate CP low-rank functional tensor measured on a discrete time grid. Incorporating tensor algebra and the theory of Reproducing Kernel Hilbert Space (RKHS), we propose a novel RKHS-based constrained power iteration with spectral initialization. Our method can successfully estimate both singular vectors and functions of the low-rank structure in the observed data. With mild assumptions, we establish the non-asymptotic contractive error bounds for the proposed algorithm. The superiority of the proposed framework is demonstrated via extensive experiments on both simulated and real data.