Search

The decay of kinetic helicity is studied in numerical models of forced turbulence using either an externally imposed forcing function as an inhomogeneous term in the equations or, alternatively, a term linear in the velocity giving rise to a linear instability. The externally imposed forcing function injects energy at the largest scales, giving rise to a turbulent inertial range with nearly constant energy flux while for linearly forced turbulence the spectral energy is maximum near the dissipation wavenumber. Kinetic helicity is injected once a statistically steady state is reached, but it is found to decay on a turbulent time scale regardless of the nature of the forcing and the value of the Reynolds number., Comment: 7 pages, 8 figures, Astron. Nachr. (accepted)

Turbulence and chaos play a fundamental role in stellar convective zones through the transportof particles, energy and momentum, and in fast dynamos, through the stretching, twisting and folding of magnetic flux tubes. A particularly revealing way to describe turbulent motions is through the analysis of Lagrangian coherent structures (LCS), which are material lines or surfaces that act as transport barriers in the fluid. We report the detection of Lagrangian coherent structures in helical MHD dynamo simulations with scale separation. In an ABC--flow, two dynamo regimes, a propagating coherent mean--field regime and an intermittent regime, are identified as the magnetic diffusivity is varied. The sharp contrast between the chaotic tangle of attracting and repelling LCS in both regimes permits a unique analysis of the impact of the magnetic field on the velocity field. In a second example, LCS reveal the link between the level of chaotic mixing of the velocity field and the saturation of a large--scale dynamo when the magnetic field exceeds the equipartition value.

We present a unified three-dimensional model of the convection zone and upper atmosphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity of the small--scale field is transported outwards and form large structures like magnetic clouds., Comment: 4 pages, 4 figures, proceedings of IAU Symp. 274, Advances in Plasma Astrophysics

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., Comment: consistent with the published version

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., Comment: 14 pages, 19 figures 3 tables (submitted to A & A)

Two setups with interlocked magnetic flux tubes are used to study the evolution of magnetic energy and helicity on magnetohydrodynamical (MHD) systems like plasmas. In one setup the initial helicity is zero while in the other it is finite. To see if it is the actual linking or merely the helicity content that influences the dynamics of the system we also consider a setup with unlinked field lines as well as a field configuration in the shape of a trefoil knot. For helical systems the decay of magnetic energy is slowed down by the helicity which decays slowly. It turns out that it is the helicity content, rather than the actual linking, that is significant for the dynamics., Comment: 3 pages, 6 figures, proceedings of IAU Symp. 274, Advances in Plasma Astrophysics, ed. A. Bonanno, E. de Gouveia dal Pino and A. Kosovichev

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential., Comment: 11 pages, 10 figures, submitted to Physics of Plasmas

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars., Comment: 4 pages, 1 figure, to appear in the proceedings of IAU Symposium 273 "Physics of Sun and Starspots"

Magnetic buoyancy is believed to drive the transport of magnetic flux tubes from the convection zone to the surface of the Sun. The magnetic fields form twisted loop-like structures in the solar atmosphere. In this paper we use helical forcing to produce a large-scale dynamo-generated magnetic field, which rises even without magnetic buoyancy. A two layer system is used as computational domain where the upper part represents the solar atmosphere. Here, the evolution of the magnetic field is solved with the stress--and--relax method. Below this region a magnetic field is produced by a helical forcing function in the momentum equation, which leads to dynamo action. We find twisted magnetic fields emerging frequently to the outer layer, forming arch-like structures. In addition, recurrent plasmoid ejections can be found by looking at space--time diagrams of the magnetic field. Recent simulations in spherical coordinates show similar results., Comment: 4 pages, 8 figures, To appear in the proceedings of the IAU273 "Physics of Sun and Star Spots"

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection., Comment: To appear in the proceedings of the IAU272 "Active OB stars: structure, evolution, mass loss and critical limits"

In the interstellar medium the turbulence is believed to be forced mostly through supernova explosions. In a first approximation these flows can be written as a gradient of a potential being thus devoid of vorticity. There are several mechanisms that could lead to vorticity generation, like viscosity and baroclinic terms, rotation, shear and magnetic fields, but it is not clear how effective they are, neither is it clear whether the vorticity is essential in determining the turbulent diffusion acting in the ISM. Here we present a study of the role of rotation, shear and baroclinicity in the generation of vorticity in the ISM., Comment: 2 pages, 1 figure, to be published in Proceedings of IAU Symp. 271, Astrophysical Dynamics: from Stars to Galaxies, ed. N. Brummell and A.S. Brun, CUP

We investigate the emergence of a large-scale magnetic field. This field is dynamo-generated by turbulence driven with a helical forcing function. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Time series of the magnetic field structure show recurrent plasmoid ejections., Comment: 2 pages, 2 figures, conference proceeding

In the absence of rotation and shear, and under the assumption of constant temperature or specific entropy, purely potential forcing by localized expansion waves is known to produce irrotational flows that have no vorticity. Here we study the production of vorticity under idealized conditions when there is rotation, shear, or baroclinicity, to address the problem of vorticity generation in the interstellar medium in a systematic fashion. We use three-dimensional periodic box numerical simulations to investigate the various effects in isolation. We find that for slow rotation, vorticity production in an isothermal gas is small in the sense that the ratio of the root-mean-square values of vorticity and velocity is small compared with the wavenumber of the energy-carrying motions. For Coriolis numbers above a certain level, vorticity production saturates at a value where the aforementioned ratio becomes comparable with the wavenumber of the energy-carrying motions. Shear also raises the vorticity production, but no saturation is found. When the assumption of isothermality is dropped, there is significant vorticity production by the baroclinic term once the turbulence becomes supersonic. In galaxies, shear and rotation are estimated to be insufficient to produce significant amounts of vorticity, leaving therefore only the baroclinic term as the most favorable candidate. We also demonstrate vorticity production visually as a result of colliding shock fronts., Comment: 9 pages, 10 figures, Accepted for publication in A&A

Observations have shown that the Sun's magnetic field has helical structures. The helicity content in magnetic field configurations is a crucial constraint on the dynamical evolution of the system. Since helicity is connected with the number of links we investigate configurations with interlocked magnetic flux rings and one with unlinked rings. It turns out that it is not the linking of the tubes which affects the magnetic field decay, but the content of magnetic helicity., Comment: 2 pages, 3 figures, proceedings of IAU Symp. 271, Astrophysical Dynamics: from Stars to Galaxies, ed. N. Brummell and A.S. Brun, CUP

We study the ability of magnetic helicity expulsion to alleviate catastrophic $\alpha$-quenching in mean field dynamos in two--dimensional spherical wedge domains. Motivated by the physical state of the outer regions of the Sun, we consider $\alpha^2\Omega$ mean field models with a dynamical $\alpha$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow ("solar wind") and meridional circulation. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain parameter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field., Comment: consistent with the published version

Magnetic helicity has risen to be a major player in dynamo theory, with the helicity of the small-scale field being linked to the dynamo saturation process for the large-scale field. It is a nearly conserved quantity, which allows its evolution equation to be written in terms of production and flux terms. The flux term can be decomposed in a variety of fashions. One particular contribution that has been expected to play a significant role in dynamos in the presence of mean shear was isolated by Vishniac & Cho (2001, ApJ 550, 752). Magnetic helicity fluxes are explicitly gauge dependent however, and the correlations that have come to be called the Vishniac-Cho flux were determined in the Coulomb gauge, which turns out to be fraught with complications in shearing systems. While the fluxes of small-scale helicity are explicitly gauge dependent, their divergences can be gauge independent. We use this property to investigate magnetic helicity fluxes of small-scale field through direct numerical simulations in a shearing-box system and find that in a numerically usable gauge the divergence of the small-scale helicity flux vanishes, while the divergence of the Vishniac-Cho flux remains finite. We attribute this seeming contradiction to the existence of horizontal fluxes of small-scale magnetic helicity with finite divergences even in our shearing-periodic domain., Comment: 8 pages, 5 figures, Accepted, ApJ

Autocatalysis is a fundamental concept, used in a wide range of domains. From the most general definition of autocatalysis, that is a process in which a chemical compound is able to catalyze its own formation, several different systems can be described. We detail the different categories of autocatalyses, and compare them on the basis of their mechanistic, kinetic, and dynamic properties. It is shown how autocatalytic patterns can be generated by different systems of chemical reactions. The notion of autocatalysis covering a large variety of mechanistic realisations with very similar behaviors, it is proposed that the key signature of autocatalysis is its kinetic pattern expressed in a mathematical form., Comment: 10 pages, 5 figures, Submitted to J Phys Chem A

The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic alpha operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the alpha layer. We perform a systematic study of a simple axisymmetric two-layer alpha-omega dynamo in a spherical shell for Reynolds numbers in the range 1 < Rm < 10^5. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary., Comment: 8 pages, 5 figures (Accepted in Geophysical and Astrophysical Fluid Dynamics)

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to be converged for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough., Comment: 15 pages, 18 figures, 2 tables, ApJ, accepted

Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we want to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. We consider a model with just radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic alpha owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic alpha effect that feeds back on the kinetic alpha effect. In some cases we compare with results obtained using a simple algebraic alpha quenching formula. In agreement with earlier findings, the magnetic alpha effect in the dynamical alpha quenching formalism has the opposite sign compared with the kinetic alpha effect and leads to a catastrophic decrease of the saturation field strength with increasing magnetic Reynolds numbers. However, at high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward due to the opposite sign of alpha. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching., Comment: 9 pages, 14 figures, submitted to A & A

We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux., Comment: 12 pages, 12 figures. Submitted to MNRAS

We perform direct numerical simulations of forced and freely decaying 3D magnetohydrodynamic turbulence in order to model magnetic field evolution during cosmological phase transitions in the early Universe. Our approach assumes the existence of a magnetic field generated either by a process during inflation or shortly thereafter, or by bubble collisions during a phase transition. We show that the final configuration of the magnetic field depends on the initial conditions, while the velocity field is nearly independent of initial conditions., Comment: 10 pages, 6 figures, references added, PRD accepted

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number., Comment: 10 pages, 12 figures, published version. Version with higher resolution figures can be found at http://www.helsinki.fi/~kapyla/publ.html

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields., Comment: 88 pages, 26 figures, published in Reports on Progress in Physics

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges., Comment: 6 pages 5 figures

The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic helicity need to be considered to characterize the decay of the field., Comment: 7 pages, 10 figures

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large-scale magnetic structures. Using three-dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two-dimensional and three-dimensional mean-field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large-scale instability results in the formation of loop-like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean-field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena., Comment: 9 pages, 10 figures, Astron. Nachr. (submitted)

Using the test-field method for nearly irrotational turbulence driven by spherical expansion waves it is shown that the turbulent magnetic diffusivity increases with magnetic Reynolds numbers. Its value levels off at several times the rms velocity of the turbulence multiplied by the typical radius of the expansion waves. This result is discussed in the context of the galactic mean-field dynamo., Comment: 2 pages, 1 figure, to appear in "Magnetic Fields in Diffuse Media", proceedings of Joint Discussion at the 2009 XXVII IAU General Assembly in Rio de Janeiro from 12 to 14 August, 2009

The role of magnetic helicity in astrophysical large-scale dynamos is reviewed and compared with cases where there is no energy supply and an initial magnetic field can only decay. In both cases magnetic energy tends to get redistributed to larger scales. Depending on the efficiency of magnetic helicity fluxes, the decay of a helical field can speed up. Likewise, the saturation of a helical dynamo can speed up through magnetic helicity fluxes. The astrophysical importance of these processes is reviewed in the context of the solar dynamo and an estimated upper limit for the magnetic helicity flux of 10^46 Mx^2/cycle is given., Comment: 14 pages, 4 figure, Proceedings of 36th EPS Conference on Plasma Physics, to appear in Plasma Physics and Controlled Fusion

The effects of uniform horizontal shear on a stably stratified layer of gas is studied. The system is initially destabilized by a magnetically buoyant flux tube pointing in the cross-stream direction. The shear amplifies the initial field to Lundquist numbers of about 200-400, but then its value drops to about 100-300, depending on the value of the sub-adiabatic gradient. The larger values correspond to cases where the stratification is strongly stable and nearly isothermal. At the end of the runs the magnetic field is nearly axisymmetric, i.e. uniform in the streamwise direction. In view of Cowling's theorem the sustainment of the field remains a puzzle and may be due to subtle numerical effects that have not yet been identified in detail. In the final state the strength of the magnetic field decreases with height in such a way that the field is expected to be unstable. Low amplitude oscillations are seen in the vertical velocity even at late times, suggesting that they might be persistent., Comment: 10 pages, 16 figures

The mechanisms for explaining how a stable asymmetric chemical system can be formed from a symmetric chemical system, in the absence of any asymmetric influence other than statistical fluctuations, have been developed during the last decades, focusing on the non-linear kinetic aspects. Besides the absolute necessity of self-amplification processes, the importance of energetic aspects is often underestimated. Going down to the most fundamental aspects, the distinction between a single object -- that can be intrinsically asymmetric -- and a collection of objects -- whose racemic state is the more stable one -- must be emphasized. A system of strongly interacting objects can be described as one single object retaining its individuality and a single asymmetry; weakly or non-interacting objects keep their own individuality, and are prone to racemize towards the equilibrium state. In the presence of energy fluxes, systems can be maintained in an asymmetric non-equilibrium steady-state. Such dynamical systems can retain their asymmetry for times longer than their racemization time., Comment: 8 pages, 7 figures, submitted to Origins of Life and Evolution of Biospheres

Using two- and three-dimensional hydromagnetic simulations for a range of different flows, including laminar and turbulent ones, it is shown that solutions expressing the field in terms of Euler potentials (EP) are in general incorrect if the EP are evolved with an artificial diffusion term. In three dimensions, standard methods using the magnetic vector potential are found to permit dynamo action when the EP give decaying solutions. With an imposed field, the EP method yields excessive power at small scales. This effect is more exaggerated in the dynamic case, suggesting an unrealistically reduced feedback from the Lorentz force. The EP approach agrees with standard methods only at early times when magnetic diffusivity did not have time to act. It is demonstrated that the usage of EP with even a small artificial magnetic diffusivity does not converge to a proper solution of hydromagnetic turbulence. The source of this disagreement is not connected with magnetic helicity or the three-dimensionality of the magnetic field, but is simply due to the fact that the nonlinear representation of the magnetic field in terms of EP that depend on the same coordinates is incompatible with the linear diffusion operator in the induction equation., Comment: 8 pages, 10 figures, MNRAS, accepted

The turbulent diffusivity tensor is determined for linear shear flow turbulence using numerical simulations. For moderately strong shear, the diagonal components are found to increase quadratically with Peclet and Reynolds numbers below about 10 and then become constant. The diffusivity tensor is found to have components proportional to the symmetric and antisymmetric parts of the velocity gradient matrix, as well as products of these. All components decrease with the wave number of the mean field in a Lorentzian fashion. The components of the diffusivity tensor are found not to depend significantly on the presence of helicity in the turbulence. The signs of the leading terms in the expression for the diffusion tensor are found to be in good agreement with estimates based on a simple closure assumption., Comment: 8 pages, 7 figures, accepted by Phys. Rev. E

Using mean-field models with a dynamical quenching formalism we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result either from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity concentration, or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong both for large-scale and small-scale fields. In the presence of shear there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behavior emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain., Comment: 9 pages, 10 figures, submitted to MNRAS

The degree of intermittency of the magnetic field of a large-scale dynamo is considered. Based on simulations it is argued that there is a tendency for the field to become more diffuse and non-intermittent as the dynamo saturates. The simulations are idealized in that the turbulence is strongly helical and shear is strong, so the tendency for the field to become more diffuse is somewhat exaggerated. Earlier results concerning the effects of magnetic buoyancy are discussed. It is emphasized that the resulting magnetic buoyancy is weak compared with the stronger effects of simultaneous downward pumping. These findings are used to support the notion that the solar dynamo might operate in a distributed fashion where the near-surface shear layer could play an important role., Comment: 10 pages, 7 figures, to be published in Solar-stellar dynamos as revealed by helio- and asteroseismology, ed. M. Dikpati et al

The generation of mean magnetic fields is studied for a simple non-helical flow where a net cross helicity of either sign can emerge. This flow, which is also known as the Archontis flow, is a generalization of the Arnold--Beltrami--Childress flow, but with the cosine terms omitted. The presence of cross helicity leads to a mean-field dynamo effect that is known as the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that exceed a certain critical value. Unlike $\alpha$ effect dynamos, it is found that the Yoshizawa effect is not noticeably constrained by the presence of a conservation law. It is argued that this is due to the presence of a forcing term in the momentum equation which leads to a nonzero correlation with the magnetic field. Finally, the application to energy convergence in solar wind turbulence is discussed., Comment: Accepted for publication in MNRAS

Selected topics in solar dynamo theory are being highlighted. The possible relevance of the near-surface shear layer is discussed. The role of turbulent downward pumping is mentioned in connection with earlier concerns that a dynamo-generated magnetic field would be rapidly lost from the convection zone by magnetic buoyancy. It is argued that shear-mediated small-scale magnetic helicity fluxes are responsible for the success of some of the recent large-scale dynamo simulations. These fluxes help in disposing of excess small-scale magnetic helicity. This small-scale magnetic helicity, in turn, is generated in response to the production of an overall tilt in each Parker loop. Some preliminary calculations of this helicity flux are presented for a system with uniform shear. In the Sun the effects of magnetic helicity fluxes may be seen in coronal mass ejections shedding large amounts of magnetic helicity., Comment: 7 pages, 4 figures, final version

The turbulent diffusion tensor describing the evolution of the mean concentration of a passive scalar is investigated for non-helically forced turbulence in the presence of rotation or a magnetic field. With rotation, the Coriolis force causes a sideways deflection of the flux of mean concentration. Within the magnetohydrodynamics approximation there is no analogous effect from the magnetic field because the effects on the flow do not depend on the sign of the field. rotation and magnetic fields tend to suppress turbulent transport, but this suppression is weaker in the direction along the magnetic field. Turbulent transport along the rotation axis is not strongly affected by rotation, except on shorter length scales, i.e., when the scale of the variation of the mean field becomes comparable with the scale of the energy-carrying eddies. These results are discussed in the context of anisotropic convective energy transport in the Sun., Comment: 8 page, 9 figures, compatible with published version

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo., Comment: Consistent with the published verstion (emulateapj)

Recent analytical and computational advances in the theory of large-scale dynamos are reviewed. The importance of the magnetic helicity constraint is apparent even without invoking mean-field theory. The tau approximation yields expressions that show how the magnetic helicity gets incorporated into mean-field theory. The test-field method allows an accurate numerical determination of turbulent transport coefficients in linear and nonlinear regimes. Finally, some critical views on the solar dynamo are being offered and targets for future research are highlighted., Comment: 19 pages, 2 figures, to appear in Space Sci. Rev

We use three-dimensional direct numerical simulations of the helically forced magnetohydrodynamic equations in spherical shell segments in order to study the effects of changes in the geometrical shape and size of the domain on the growth and saturation of large-scale magnetic fields. We inject kinetic energy along with kinetic helicity in spherical domains via helical forcing using Chandrasekhar-Kendall functions. We take perfect conductor boundary conditions for the magnetic field to ensure that no magnetic helicity escapes the domain boundaries. We find dynamo action giving rise to magnetic fields at scales larger than the characteristic scale of the forcing. The magnetic energy exceeds the kinetic energy over dissipative time scales, similar to that seen earlier in Cartesian simulations in periodic boxes. As we increase the size of the domain in the azimuthal direction we find that the nonlinearly saturated magnetic field organizes itself in long-lived cellular structures with aspect ratios close to unity. These structures tile the domain along the azimuthal direction, thus resulting in very small longitudinally averaged magnetic fields for large domain sizes. The scales of these structures are determined by the smallest scales of the domain, which in our simulations is usually the radial scale. We also find that increasing the meridional extent of the domains produces little qualitative change, except a marginal increase in the large-scale field. We obtain qualitatively similar results in Cartesian domains with similar aspect ratios., Comment: Updated with minor changes to be consistent with the published version

In the mean-field theory of magnetic fields, turbulent transport, i.e. the turbulent electromotive force, is described by a combination of the alpha effect and turbulent magnetic diffusion, which are usually assumed to be proportional respectively to the mean field and its spatial derivatives. For a passive scalar there is just turbulent diffusion, where the mean flux of concentration depends on the gradient of the mean concentration. However, these proportionalities are approximations that are valid only if the mean field or the mean concentration vary slowly in time. Examples are presented where turbulent transport possesses memory, i.e. where it depends crucially on the past history of the mean field. Such effects are captured by replacing turbulent transport coefficients with time integral kernels, resulting in transport coefficients that depend effectively on the frequency or the growth rate of the mean field itself. In this paper we perform numerical experiments to find the characteristic timescale (or memory length) of this effect as well as simple analytical models of the integral kernels in the case of passive scalar concentrations and kinematic dynamos. The integral kernels can then be used to find self-consistent growth or decay rates of the mean fields. In mean-field dynamos the growth rates and cycle periods based on steady state values of alpha effect and turbulent diffusivity can be quite different from the actual values., Comment: 16 pages, 16 figures, published version

Simulations of stochastically forced shear-flow turbulence in a shearing-periodic domain are used to study the spontaneous generation of large-scale flow patterns in the direction perpendicular to the plane of the shear. Based on an analysis of the resulting large-scale velocity correlations it is argued that the mechanism behind this phenomenon could be the mean-vorticity dynamo effect pioneered by Elperin, Kleeorin, and Rogachevskii in 2003 (Phys. Rev. E 68, 016311). This effect is based on the anisotropy of the eddy viscosity tensor. One of its components may be able to replenish cross-stream mean flows by acting upon the streamwise component of the mean flow. Shear, in turn, closes the loop by acting upon the cross-stream mean flow to produce stronger streamwise mean flows. The diagonal component of the eddy viscosity is found to be of the order of the rms turbulent velocity divided by the wavenumber of the energy-carrying eddies., Comment: 8 pages, 11 figures, Phys. Rev. E (in press)

The coefficients defining the mean electromotive force in a Galloway-Proctor flow are determined. This flow shows a two-dimensional pattern and is helical. The pattern wobbles in its plane. Apart from one exception a circular motion of the flow pattern is assumed. This corresponds to one of the cases considered recently by Courvoisier, Hughes and Tobias (2006, Phys. Rev. Lett., 96, 034503). An analytic theory of the alpha effect and related effects in this flow is developed within the second-order correlation approximation and a corresponding fourth-order approximation. In the validity range of these approximations there is an alpha effect but no gamma effect, or pumping effect. Numerical results obtained with the test-field method, which are independent of these approximations, confirm the results for alpha and show that gamma is in general nonzero. Both alpha and gamma show a complex dependency on the magnetic Reynolds number and other parameters that define the flow, that is, amplitude and frequency of the wobbling motion. Some results for the magnetic diffusivity eta_t and a related quantity are given, too. Finally a result for alpha in the case of a randomly varying flow without the aforementioned circular motion is presented. This flow may be a more appropriate model for studying the alpha effect and related effects in flows that are statistical isotropic in a plane., Comment: 12 pages, 14 figures, submitted to MNRAS