The Kelvin-Hemholtz (K-H) instability in the corona EUV jet is studied via 2.5 MHD numerical simulations. The jet results from magnetic reconnection due to the interation of the new emerging magnetic field and the pre-existing magnetic field in the corona. Our results show that the Alfven Mach number along the jet is about 5-14 just before the instability occurs, and it is even higher than 14 at some local areas. During the K-H instability process, several vortex-like plasma blobs of high temperature and high density appear along the jet, and magnetic fields have also been rolled up and the magnetic configuration including anti-parallel magnetic fields forms, which leads to magnetic reconnection at many X-points and current sheet fragments inside the vortex-like blob. After magnetic islands appear inside the main current sheet, the total kinetic energy of the reconnection outflows decreases, and cannot support the formation of the vortex-like blob along the jet any longer, then the K-H instability eventually disappears. We also present the results about how the guide field and the flux emerging speed affect the K-H instability. We find that the strong guide field inhibits the shock formation in the reconnection upward outflow regions but helps secondary magnetic islands appear earlier in the main current sheet, and then apparently suppresses the K-H instability. As the speed of the emerging magnetic field decreases, the K-H instability appears later, the highest temperature inside the vortex blob gets lower and the vortex structure gets smaller.

Context: Long-term variability in solar cycles represents a challenging constraint for theoretical models. Mean-field Babcock-Leighton dynamos that consider non-instantaneous rising flux tubes have been shown to exhibit long-term variability in their magnetic cycle. However a relation that parameterizes the rise-time of non-axisymmetric magnetic flux tubes in terms of stellar parameters is still missing. Aims: We aim to find a general parameterization of the rise-time of magnetic flux tubes for solar-like stars. Methods: By considering the influence of magnetic tension on the rise of non-axisymmetric flux tubes, we predict the existence of a control parameter referred as $\Gamma_{\alpha_1}^{\alpha_2}$. This parameter is a measure of the balance between rotational effects and magnetic effects (buoyancy and tension) acting on the magnetic flux tube. We carry out two series of numerical experiments (one for axisymmetric rise and one for non-axisymmetric rise) and demonstrate that $\Gamma_{\alpha_1}^{\alpha_2}$ indeed controls the rise-time of magnetic flux tubes. Results: We find that the rise-time follows a power law of $\Gamma_{\alpha_1}^{\alpha_2}$ with an exponent that depends on the azimuthal wavenumber of the magnetic flux loop. Conclusions: Compressibility does not impact the rise of magnetic flux tubes, while non-axisymmetry does. In the case of non-axisymmetric rise, the tension force modifies the force balance acting on the magnetic flux tube. We identified the three independent parameters required to predict the rise-time of magnetic flux tubes, that is, the stellar rotation rate, the magnetic flux density of the flux tube, and its azimuthal wavenumber. We combined these into one single relation that is valid for any solar-like star. We suggest using this generalized relation to constrain the rise-time of magnetic flux tubes in Babcock-Leighton dynamo models.

(Abridged) Based on the rapidly increasing all-sky data of Faraday rotation measures and polarised synchrotron radiation, the Milky Way's magnetic field is now modelled with an unprecedented level of detail and complexity. We aim to complement this heuristic approach with a physically motivated, quantitative Galactic dynamo model -- a model that moreover allows for the evolution of the system as a whole, instead of just solving the induction equation for a fixed static disc. Building on the framework of mean-field magnetohydrodynamics and extending it to the realm of a hybrid evolution, we perform three-dimensional global simulations of the Galactic disc. Closure coefficients embodying the mean-field dynamo are calibrated against resolved box simulations of supernova-driven interstellar turbulence. The emerging dynamo solutions comprise a mixture of the dominant axisymmetric S0 mode, with even parity, and a subdominant A0 mode, with odd parity. Notably, such a superposition of modes creates a strong localised vertical field on one side of the Galactic disc. We moreover find significant radial pitch angles, which decay with radius -- explained by flaring of the disc. In accordance with previous work, magnetic instabilities appear to be restricted to the less-stirred outer Galactic disc. Their main effect is to create strong fields at large radii such that the radial scale length of the magnetic field increases from 4 kpc (for the case of a mean-field dynamo alone) to about 10 kpc in the hybrid models. There remain aspects (e.g., spiral arms, X-shaped halo fields, fluctuating fields) that are not captured by the current model and that will require further development towards a fully dynamical evolution. Nevertheless, the work presented demonstrates that a hybrid modelling of the Galactic dynamo is feasible and can serve as a foundation for future efforts.

We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPD's surface density $\sim30{\rm\,g\,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 \times 10^{-3} {\rm M_\oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

In this paper we investigate, by means of two-dimensional magnetohydrodynamic simulations, the impact of temperature-dependent resistivity and thermal conduction on the development of plasmoid instabilities in reconnecting current sheets in the solar corona. We find that the plasma temperature in the current sheet region increases with time and it becomes greater than that in the inflow region. As secondary magnetic islands appear, the highest temperature is not always found at the reconnection $X$-points, but also inside the secondary islands. One of the effects of anisotropic thermal conduction is to decrease the temperature of the reconnecting $X-$points and transfer the heat into the $O-$points, the plasmoids, where it gets trapped. In the cases with temperature-dependent magnetic diffusivity, $\eta \sim T^{-3/2}$, the decrease in plasma temperature at the $X-$points leads to: (i) increase in the magnetic diffusivity until the characteristic time for magnetic diffusion becomes comparable to that of thermal conduction; (ii) increase in the reconnection rate; and, (iii) more efficient conversion of magnetic energy into thermal energy and kinetic energy of bulk motions. These results provide further explanation of the rapid release of magnetic energy into heat and kinetic energy seen during flares and coronal mass ejections. In this work, we demonstrate that the consideration of anisotropic thermal conduction and Spitzer-type, temperature-dependent magnetic diffusivity, as in the real solar corona, are crucially important for explaining the occurrence of fast reconnection during solar eruptions.

We study the evolution of isolated and merging interstellar bubbles of three stars (25, 32 and 60 M\odot) in a homogeneous background medium with a density of 10 mp /ccm via 3D-hydrodynamic simulations with standard ISM thermodynamics (optically thin radiative cooling and photo-electric heating) and time dependent energy and mass input according to stellar evolutionary tracks. We vary the position of the three stars relative to each other to compare the energy response for cases of isolated, merging and initially cospatial bubbles. Due to mainly the Vishniac instability, our simulated bubbles develop thick shells and filamentary internal structures in column density. The shell widths reach tens of per cent of the outer bubble radius, which compares favourably to observations. More energy is retained in the ISM for more closely packed groups, by up to a factor of three and typically a factor of two for intermediate times after the first supernova. Once the superbubble is established, different positions of the contained stars make only a minor difference to the energy tracks. For our case of three massive stars, the energy deposition varies only very little for distances up to about 30 pc between the stars. Energy injected by supernovae is entirely dissipated in a superbubble on a timescale of about 1 Myr, which increases slightly with the superbubble size at the time of the explosion. The Vishniac instability may be responsible for the broadening of the shells of interstellar bubbles. Massive star winds are significant energetically due to their - in the long run - more efficient, steady energy injection and because they evacuate the space around the massive stars. For larger scale simulations, the feedback effect of close groups of stars or clusters may be subsumed into one effective energy input with insignificant loss of energy accuracy.

Supernovae are the dominant energy source for driving turbulence within the interstellar plasma. Until recently, their effects on magnetic field amplification in disk galaxies remained a matter of speculation. By means of self-consistent simulations of supernova-driven turbulence, we find an exponential amplification of the mean magnetic field on timescales of a few hundred million years. The robustness of the observed fast dynamo is checked at different magnetic Reynolds numbers, and we find sustained dynamo action at moderate Rm. This indicates that the mechanism might indeed be of relevance for the real ISM. Sensing the flow via passive tracer fields, we infer that SNe produce a turbulent alpha effect which is consistent with the predictions of quasilinear theory. To lay a foundation for global mean-field models, we aim to explore the scaling of the dynamo tensors with respect to the key parameters of our simulations. Here we give a first account on the variation with the supernova rate.

Within the interstellar medium, supernovae are thought to be the prevailing agents in driving turbulence. Until recently, their effects on magnetic field amplification in disk galaxies remained uncertain. Analytical models based on the uncorrelated-ensemble approach predicted that any created field would be expelled from the disk before it could be amplified significantly. By means of direct simulations of supernova-driven turbulence, we demonstrate that this is not the case. Accounting for galactic differential rotation and vertical stratification, we find an exponential amplification of the mean field on timescales of several hundred million years. We especially highlight the importance of rotation in the generation of helicity by showing that a similar mechanism based on Cartesian shear does not lead to a sustained amplification of the mean magnetic field.

Supernovae are known to be the dominant energy source for driving turbulence in the interstellar medium. Yet, their effect on magnetic field amplification in spiral galaxies is still poorly understood. Previous analytical models, based on the evolution of isolated, non-interacting supernova remnants, predicted a dominant vertical pumping that would render dynamo action improbable. In the present work, we address the issue of vertical transport, which is thought to be the key process that inhibits dynamo action in the galactic context. We aim to demonstrate that supernova driving is a powerful mechanism to amplify galactic magnetic fields. We conduct direct numerical simulations in the framework of resistive magnetohydrodynamics. Our local box model of the interstellar medium comprises optically-thin radiative cooling, an external gravitational potential, and background shear. Dynamo coefficients for mean-field models are measured by means of passive test fields. Our simulations show that supernova-driven turbulence in conjunction with shear leads to an exponential amplification of the mean magnetic field. We found turbulent pumping to be directed inward and approximately balanced by a galactic wind.

Observations in polarized emission reveal the existence of large-scale coherent magnetic fields in a wide range of spiral galaxies. Radio-polarization data show that these fields are strongly inclined towards the radial direction, with pitch angles up to $35\degr$ and thus cannot be explained by differential rotation alone. Global dynamo models describe the generation of the radial magnetic field from the underlying turbulence via the so called $\alpha$-effect. However, these global models still rely on crude assumptions about the small-scale turbulence. To overcome these restrictions we perform fully dynamical MHD simulations of interstellar turbulence driven by supernova explosions. From our simulations we extract profiles of the contributing diagonal elements of the dynamo $\alpha$-tensor as functions of galactic height. We also measure the coefficients describing vertical pumping and find that the ratio $\hat{\gamma}$ between these two effects has been overestimated in earlier analytical work, where dynamo action seemed impossible. In contradiction to these models based on isolated remnants we always find the pumping to be directed inward. In addition we observe that $\hat{\gamma}$ depends on whether clustering in terms of super-bubbles is taken into account. Finally, we apply a test field method to derive a quantitative measure of the turbulent magnetic diffusivity which we determine to be ~ 2 kpc kms.

We describe the implementation of the shearingbox approach into the Godunov-type central-upwind/constraint-transport magnetohydrodynamics code NIRVANA. This will allow for applications which require sheared-periodic boundary conditions as typically used in local Cartesian simulations of differentially rotating systems. We present the algorithm in detail and discuss necessary modifications in the numerical fluxes in order to preserve conserved quantities and to fulfill other analytical constraints as good as seem feasible within the numerical scheme. We check the source terms which come with the shearingbox formulation by investigating the conservation of the epicyclic mode energy. We also perform more realistic simulations of the magneto-rotational instability with initial zero-net-flux vertical magnetic field and compare the obtained stresses and energetics with previous non-conservative results exploring the same parameter regime.

The equations for fully compressible rotating magnetoconvection are numerically solved in a Cartesian box assuming conditions roughly suitable for the geodynamo. The mean electromotive force describing the generation of mean magnetic flux by convective turbulence in the rotating fluid is directly calculated from the simulations, and the corresponding alpha-coefficients are derived. Due to the very weak density stratification the alpha-effect changes its sign in the middle of the box. It is positive at the top and negative at the bottom of the convection zone. For strong magnetic fields we also find a clear downward advection of the mean magnetic field. Both of the simulated effects have been predicted by quasi-linear computations (Soward, 1979; Kitchatinov and Ruediger, 1992). Finally, the possible connection of the obtained profiles of the EMF with mean-field models of oscillating alpha^2-dynamos is discussed.