Search for non-helical disc dynamos in simulations

The possibility of non-helical large scale dynamo action is investigated using three-dimensional simulations of global accretion discs as well as idealized local simulations without rotation and only shear. Particular emphasis is placed on a certain correlation between vorticity and azimuthal velocity gradient which has been predicted to drive large scale dynamo action, independent of the presence or absence of kinetic helicity. In the global disc simulations two types of behaviours are found: those which do show this type of velocity correlation and those which do not. The former ones are typically also the cases where the resistivity is larger. The latter ones show signs typical of dynamo action based on the usual helicity effect. In the idealized simulations without rotation and just shear the above correlation is found to be particularly strong. In both cases there is, as expected, a systematic flux of magnetic helicity through the midplane. However, very little helicity flux leaves the domain through the top and bottom boundaries. The idealized simulations reveal that much of this systematic flux comes from the rotational component of the helicity flux and does not contribute to its divergence.Comment: 14 pages, 17 figures, 5 tables, submitted to Astron. & Astrophy

Domain structure of epitaxial Co films with perpendicular anisotropy

Epitaxial hcp Cobalt films with pronounced c-axis texture have been prepared by pulsed lased deposition (PLD) either directly onto Al2O3 (0001) single crystal substrates or with an intermediate Ruthenium buffer layer. The crystal structure and epitaxial growth relation was studied by XRD, pole figure measurements and reciprocal space mapping. Detailed VSM analysis shows that the perpendicular anisotropy of these highly textured Co films reaches the magnetocrystalline anisotropy of hcp-Co single crystal material. Films were prepared with thickness t of 20 nm < t < 100 nm to study the crossover from in-plane magnetization to out-of-plane magnetization in detail. The analysis of the periodic domain pattern observed by magnetic force microscopy allows to determine the critical minimum thickness below which the domains adopt a pure in-plane orientation. Above the critical thickness the width of the stripe domains is evaluated as a function of the film thickness and compared with domain theory. Especially the discrepancies at smallest film thicknesses show that the system is in an intermediate state between in-plane and out-of-plane domains, which is not described by existing analytical domain models

The inverse cascade and nonlinear alpha-effect in simulations of isotropic helical hydromagnetic turbulence

A numerical model of isotropic homogeneous turbulence with helical forcing is investigated. The resulting flow, which is essentially the prototype of the alpha^2 dynamo of mean-field dynamo theory, produces strong dynamo action with an additional large scale field on the scale of the box (at wavenumber k=1; forcing is at k=5). This large scale field is nearly force-free and exceeds the equipartition value. As the magnetic Reynolds number R_m increases, the saturation field strength and the growth rate of the dynamo increase. However, the time it takes to built up the large scale field from equipartition to its final super-equipartition value increases with magnetic Reynolds number. The large scale field generation can be identified as being due to nonlocal interactions originating from the forcing scale, which is characteristic of the alpha-effect. Both alpha and turbulent magnetic diffusivity eta_t are determined simultaneously using numerical experiments where the mean-field is modified artificially. Both quantities are quenched in a R_m-dependent fashion. The evolution of the energy of the mean field matches that predicted by an alpha^2 dynamo model with similar alpha and eta_t quenchings. For this model an analytic solution is given which matches the results of the simulations. The simulations are numerically robust in that the shape of the spectrum at large scales is unchanged when changing the resolution from 30^3 to 120^3 meshpoints, or when increasing the magnetic Prandtl number (viscosity/magnetic diffusivity) from 1 to 100. Increasing the forcing wavenumber to 30 (i.e. increasing the scale separation) makes the inverse cascade effect more pronounced, although it remains otherwise qualitatively unchanged.Comment: 21 pages, 26 figures, ApJ (accepted

Mean Field Dynamos with Algebraic and Dynamic alpha-Quenchings

Calculations for mean field dynamo models (in both full spheres and spherical shells), with both algebraic and dynamic $\alpha$--quenchings, show qualitative as well as quantitative differences and similarities in the dynamical behaviour of these models. We summarise and enhance recent results with extra examples. Overall, the effect of using a dynamic $\alpha$ appears to be complicated and is affected by the region of parameter space examined.Comment: 6 pages, 2 postscript figures, also available at http://www.maths.qmw.ac.uk/~eo

Competition of rotation and stratification in flux concentrations

In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to spontaneously form local flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism, called the negative effective magnetic pressure instability (NEMPI), is of interest in connection with dynamo scenarios where most of the magnetic field resides in the bulk of the convection zone, and not at the bottom. Recent work using the mean-field hydromagnetic equations has shown that NEMPI becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1.}{Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the effective magnetic pressure. We also quantify the kinetic helicity from direct numerical simulations (DNS) and compare with earlier work.}{To calculate the rotational effect on the effective magnetic pressure we consider both DNS and analytical studies using the $\tau$ approach. To study the effects of rotation on the development of NEMPI we use both DNS and mean-field calculations of the 3D hydromagnetic equations in a Cartesian domain.}{We find that the growth rates of NEMPI from earlier mean-field calculations are well reproduced with DNS, provided the Coriolis number is below about 0.06. In that case, kinetic and magnetic helicities are found to be weak. For faster rotation, dynamo action becomes possible. However, there is an intermediate range of rotation rates where dynamo action on its own is not yet possible, but the rotational suppression of NEMPI is being alleviated.}{Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours.}Comment: 13 pages, 13 figures submitted to A&

Magnetic flux concentrations from dynamo-generated fields

The mean-field theory of magnetized stellar convection gives rise to the two possibility of distinct instabilities: the large-scale dynamo instability, operating in the bulk of the convection zone, and a negative effective magnetic pressure instability (NEMPI) operating in the strongly stratified surface layers. The latter might be important in connection with magnetic spot formation, but the growth rate of NEMPI is suppressed with increasing rotation rates, although recent direct numerical simulations (DNS) have shown a subsequent increase in the growth rate. We examine quantitatively whether this increase in the growth rate of NEMPI can be explained by an alpha squared mean-field dynamo, and whether both NEMPI and the dynamo instability can operate at the same time. We use both DNS and mean-field simulations (MFS) to solve the underlying equations numerically either with or without an imposed horizontal field. We use the test-field method to compute relevant dynamo coefficients. DNS show that magnetic flux concentrations are still possible up to rotation rates above which the large-scale dynamo effect produces mean magnetic fields. The resulting DNS growth rates are quantitatively well reproduced with MFS. As expected, for weak or vanishing rotation, the growth rate of NEMPI increases with increasing gravity, but there is a correction term for strong gravity and large turbulent magnetic diffusivity. Magnetic flux concentrations are still possible for rotation rates above which dynamo action takes over. For the solar rotation rate, the corresponding turbulent turnover time is about 5 hours, with dynamo action commencing in the layers beneath.Comment: 10 pages, 10 figures, submitted to A&

Inertial range scaling in numerical turbulence with hyperviscosity

Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case--presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t^{-1.25} decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.Comment: Final version (7 pages

Rotational effects on the negative magnetic pressure instability

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.Comment: 8 pages, 13 figures, submitted to A&