Convective Babcock-Leighton Dynamo Models

We present the first global, three-dimensional simulations of solar/stellar convection that take into account the influence of magnetic flux emergence by means of the Babcock-Leighton (BL) mechanism. We have shown that the inclusion of a BL poloidal source term in a convection simulation can promote cyclic activity in an otherwise steady dynamo. Some cycle properties are reminiscent of solar observations, such as the equatorward propagation of toroidal flux near the base of the convection zone. However, the cycle period in this young sun (rotating three times faster than the solar rate) is very short ($\sim$ 6 months) and it is unclear whether much longer cycles may be achieved within this modeling framework, given the high efficiency of field generation and transport by the convection. Even so, the incorporation of mean-field parameterizations in 3D convection simulations to account for elusive processes such as flux emergence may well prove useful in the future modeling of solar and stellar activity cycles.Comment: 5 pages, 4 figures, submitted to the Astrophysical Journal Letter

Cyclic magnetic activity due to turbulent convection in spherical wedge geometry

We report on simulations of turbulent, rotating, stratified, magnetohydrodynamic convection in spherical wedge geometry. An initially small-scale, random, weak-amplitude magnetic field is amplified by several orders of magnitude in the course of the simulation to form oscillatory large-scale fields in the saturated state of the dynamo. The differential rotation is solar-like (fast equator), but neither coherent meridional poleward circulation nor near-surface shear layer develop in these runs. In addition to a poleward branch of magnetic activity beyond 50 degrees latitude, we find for the first time a pronounced equatorward branch at around 20 degrees latitude, reminiscent of the solar cycle.Comment: 5 pages, 5 figures, version published in ApJL, version with higher resolution figures http://www.helsinki.fi/~kapyla/publ.htm

Buoyant Magnetic Loops in a Global Dynamo Simulation of a Young Sun

The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and dynamo action in a spherical shell with solar stratification, but rotating three times faster than the current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic field in these wreaths to rise toward the top of the convection zone via a combination of magnetic buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding buoyant loops represents a significant step forward in combining numerical models of dynamo action and flux emergence.Comment: 6 pages, 3 figures, accepted in ApJ Letter

Toroidal Field Reversals and the Axisymmetric Tayler Instability

We present axisymmetric numerical simulations of the solar interior, including the convection zone and an extended radiative interior. We find that differential rotation in the convection zone induces a toroidal field from an initially purely poloidal field. This toroidal field becomes unstable to the axisymmetric Tayler instability and undergoes equatorward propagating toroidal field reversals. These reversals occur in the absence of a dynamo and without accompanying poloidal field reversals. The nature and presence of such reversals depends sensitively on the initial poloidal field strength imposed, with north-south symmetric reversals only seen at a particular initial field strength. Coupled with a dynamo mechanism which regenerates the poloidal field this could be one ingredient in the sunspot cycle.Comment: accepted at Ap

Magnetic Cycles in a Convective Dynamo Simulation of a Young Solar-type Star

Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate. Striking global-scale magnetic wreaths appear in the midst of the turbulent convection zone and show rich time-dependence. The dynamo exhibits cyclic activity and undergoes quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow and are more organized temporally and spatially than those realized in our previous simulations of the solar dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily through an $\Omega$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $\alpha$-effect. During a reversal the magnetic wreaths propagate towards the polar regions, and this appears to arise from a poleward propagating dynamo wave. The primary response in the convective flows involves the axisymmetric differential rotation which shows variations associated with the poleward propagating magnetic wreaths. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. [abridged]Comment: 20 pages, 14 figures, emulateapj format; accepted for publication in ApJ. Expanded and published version of sections 5-6 from http://arxiv.org/abs/0906.240

Effects of stratification in spherical shell convection

Peer reviewe

Meridional circulation dynamics in a cyclic convective dynamo

Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to slowing down the axisymmetric poleward horizontal component. The behavior of this large-scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large-scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid-latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.Centra-ISTGRPS-UdeMNatural Sciences and Engineering Research Council of CanadaNational Science FoundationUniversity of the Algarveinfo:eu-repo/semantics/publishedVersio

Characteristics of magnetic solar-like cycles in a 3D MHD simulation of solar convection

We analyse the statistical properties of the stable magnetic cycle unfolding in an extended 3D magnetohydroclynamic simulation of solar convection produced with the EULAG-MHD code. The millennium,simulation spans over 1650 years, in the course of which forty polarity reversals take place on a regular similar to 40yr cadence, remaining well-synchronized across solar hemispheres. In order to characterize this cycle and facilitate its comparison with measures typically used to represent solar activity, we build two proxies for the magnetic field in the simulation mimicking the solar toroidal field and the polar radial field. Several quantities that characterize the cycle are measured (period, amplitudes, etc.) and correlations between them are computed, These are then compared with their observational analogs. From the typical Gnevyshesv-Ohl pattern, to hints of Gleissberg modulation the simulated cycles share many of the characteristics of their observational analogs even though the simulation lacks poloidal field regeneration through active region decay, a mechanism nowadays often considered an essential component of the solar dynamo. Some significant discrepancies are, also identified, most notably the in-phase variation of the simulated poloidal and toroidal large-scale magnetic components, and the low degree of hemispheric coupling at the level of hemispheric cycle amplitudes. Possible causes underlying these discrepancies are discussed.Fundacao para a Ciencia e Tecnologia (FCT) [SFRH/BPD/68409/2010]; CENTRA-IST; Natural Sciences and Engineering Research Council of Canadainfo:eu-repo/semantics/publishedVersio

Coupling the solar dynamo and the corona: wind properties, mass and momentum losses during an activity cycle

We study the connections between the sun's convection zone and the evolution of the solar wind and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code (STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in good agreement with the available observations. The magnetic polarity reversal happens at different paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal distribution of the sources of open flux and solar wind and the geometry of the Alfv\'en surface. Wind sources and braking torque application zones also vary accordingly

Reynolds stress and heat flux in spherical shell convection

14 pages, 15 figures, submitted to Astron. Astrophys. Version with higher resolution figures can be found at http://www.helsinki.fi/~kapyla/publ.htmlContext. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.Peer reviewe