Source author record

Allan Sacha Brun

Allan Sacha Brun appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

Exploring the Solar Poles: The Last Great Frontier of the Sun

Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere.

preprint2022arXiv

Powering Stellar Magnetism: Energy Transfers in Cyclic Dynamos of Sun-like Stars

We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ($ΔΩ\sim (|Ω|/Ω_{\odot})^{0.46}$) in the MHD solutions than in their HD counterpart $(|Ω|/Ω_{\odot})^{0.66})$, yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between $0.15 \lesssim Ro_f \lesssim 0.65$ the solutions possess long magnetic cycle, if $Ro_f \lesssim 0.42$ a short cycle and if $Ro_f \gtrsim 1$ (anti-solar-like differential rotation) a statistically steady state. We show that short-cycle dynamos follow the classical Parker-Yoshimura rule whereas the long-cycle period ones do not. We further demonstrate that the Rossby number dependency of the large-scale surface magnetic field in the simulation ($B_{L,surf} \sim Ro_{f}^{-1.26}$) agrees better with observations ($B_{V} \sim Ro_{s}^{-1.4 \pm 0.1}$) and differs from dynamo scaling based on the global magnetic energy ($B_{bulk} \sim Ro_{f}^{-0.5}$). We also show that up to few percents of the stellar luminosity can be channelled into the star's magnetism, hence providing a large energy reservoir for possible surface eruptive events.

preprint2022arXiv

Validation of a wave heated 3D MHD coronal-wind model using Polarized Brightness and EUV observations

The physical properties responsible for the formation and evolution of the corona and heliosphere are still not completely understood. 3D MHD global modeling is a powerful tool to investigate all the possible candidate processes. To fully understand the role of each of them, we need a validation process where the output from the simulations is quantitatively compared to the observational data. In this work, we present the results from our validation process applied to the wave turbulence driven 3D MHD corona-wind model WindPredict-AW. At this stage of the model development, we focus the work to the coronal regime in quiescent condition. We analyze three simulations results, which differ by the boundary values. We use the 3D distributions of density and temperature, output from the simulations at the time of around the First Parker Solar Probe perihelion (during minimum of the solar activity), to synthesize both extreme ultraviolet (EUV) and white light polarized (WL pB) images to reproduce the observed solar corona. For these tests, we selected AIA 193 A, 211 A and 171 A EUV emissions, MLSO K-Cor and LASCO C2 pB images obtained the 6 and 7 November 2018. We then make quantitative comparisons of the disk and off limb corona. We show that our model is able to produce synthetic images comparable to those of the observed corona.

preprint2020arXiv

Alfvén-wave driven magnetic rotator winds from low-mass stars I: rotation dependences of magnetic braking and mass-loss rate

Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfvén-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $τ_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( Ω_\ast / Ω_\odot \right)^{2.82}$, yielding a spin-down law $Ω_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfvén waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto Ω_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfvén radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfvén radius matches previous works which used thermally-driven winds. Our results suggest that the Alfvén-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence.

preprint2019arXiv

Exoplanet host-star properties: the active environment of exoplanets

The primary objectives of the ExoplANETS-A project are to: establish new knowledge on exoplanet atmospheres; establish new insight on influence of the host star on the planet atmosphere; disseminate knowledge, using online, web-based platforms. The project, funded under the EU's Horizon-2020 programme, started in January 2018 and has a duration ~3 years. We present an overview of the project, the activities concerning the host stars and some early results on the host stars.

preprint2016arXiv

Age dependence of wind properties for solar type stars: a 3d study

Young and rapidly rotating stars are known for intense, dynamo generated magnetic fields. Spectropolarimetric observations of those stars in precisely aged clusters are key input for gyrochronology and magnetochronology. We use ZDI maps of several young K-type stars of similar mass and radius but with various ages and rotational periods, to perform 3D numerical MHD simulations of their coronae and follow the evolution of their magnetic properties with age. Those simulations yield the coronal structure as well as the instant torque exerted by the magnetized, rotating wind on the star. As stars get older, we find that the angular momentum loss decreases with $Ω^3$, which is the reason for the convergence on the Skumanich law. For the youngest stars of our sample, the angular momentum loss show signs of saturation around $8Ω_{\odot}$, which is a common value used in spin evolution models for K-type stars. We compare these results to semi-analytical models and existing braking laws. We observe a complex wind speed distribution for the youngest stars with slow, intermediate and fast wind components, which are the result of the interaction with intense and non axisymmetric magnetic fields. Consequently, in our simulations, the stellar wind structure in the equatorial plane of young stars varies significantly from a solar configuration, delivering insight about the past of the solar system interplanetary medium.

preprint2015arXiv

From solar to stellar corona: the role of wind, rotation and magnetism

Observations of surface magnetic fields are now within reach for many stellar types thanks to the development of Zeeman-Doppler Imaging. These observations are extremely useful for constraining rotational evolution models of stars, as well as for characterizing the generation of magnetic field. We recently demonstrated that the impact of coronal magnetic field topology on the rotational braking of a star can be parametrized with a scalar parameter: the open magnetic flux. However, without running costly numerical simulations of the stellar wind, reconstructing the coronal structure of the large scale magnetic field is not trivial. An alternative -broadly used in solar physics- is to extrapolate the surface magnetic field assuming a potential field in the corona, to describe the opening of the field lines by the magnetized wind. This technique relies on the definition of a so-called source surface radius, which is often fixed to the canonical value of 2.5Rsun. However this value likely varies from star to star. To resolve this issue, we use our extended set of 2.5D wind simulations published in 2015, to provide a criteria for the opening of field lines as well as a simple tool to assess the source surface radius and the open magnetic flux. This allows us to derive the magnetic torque applied to the star by the wind from any spectropolarimetric observation. We conclude by discussing some estimations of spin-down time scales made using our technique, and compare them to observational requirements.

preprint2015arXiv

Grand Minima and Equatorward Propagation in a Cycling Stellar Convective Dynamo

The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $α$-$Ω$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified.

preprint2014arXiv

Modelling the Corona of HD 189733 in 3D

The braking of main sequence stars originates mainly from their stellar wind. The efficiency of this angular momentum extraction depends on the rotation rate of the star, the acceleration profile of the wind and the coronal magnetic field. The derivation of scaling laws parametrizing the stellar wind torque is important for our understanding of gyro-chronology and the evolution of the rotation rates of stars. In order to understand the impact of complex magnetic topologies on the stellar wind torque, we present three- dimensional, dynamical simulations of the corona of HD 189733. Using the observed complex topology of the magnetic field, we estimate how the torque associated with the wind scales with model parameters and compare those trends to previously published scaling laws.

preprint2014arXiv

On the diversity of magnetic interactions in close-in star-planet systems

Magnetic interactions between close-in planets and their host star can play an important role in the secular orbital evolution of the planets, as well as the rotational evolution of their host. As long as the planet orbits inside the Alfvén surface of the stellar wind, the magnetic interaction between the star and the planet can modify the wind properties and also lead to direct angular momentum transfers between the two. We model these star-planet interactions using compressible magneto-hydrodynamic (MHD) simulations, and quantify the angular momentum transfers between the star, the planet, and the stellar wind. We study the cases of magnetized and non-magnetized planets and vary the orbital radius inside the Alfvén surface of the stellar wind. Based on a grid of numerical simulations, we propose general scaling laws for the modification of the stellar wind torque, for the torque between the star and the planet, and for the planet migration associated with the star-planet magnetic interactions. We show that when the coronal magnetic field is large enough and the star is rotating sufficiently slowly, the effect of the magnetic star-planet interaction is comparable to tidal effects and can lead to a rapid orbital decay.

preprint2014arXiv

The effect of magnetic topology on thermally-driven winds: towards a general formulation of the braking law

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486- 4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI).

preprint2014arXiv

The influence of the magnetic topology on the braking of sun-like stars

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486-4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI). Ideas about how to compute the open flux from ZDI Maps are discussed.

preprint2014arXiv

Theoretical seismology in 3D : nonlinear simulations of internal gravity waves in solar-like stars

Internal gravity waves (hereafter IGWs) are studied for their impact on the angular momentum transport in stellar radiation zones and the information they provide about the structure and dynamics of deep stellar interiors. We here present the first 3D nonlinear numerical simulations of IGWs excitation and propagation in a solar-like star. The aim is to study the behavior of waves in a realistic 3D nonlinear time dependent model of the Sun and to characterize their properties. We compare our results with theoretical and 1D predictions. It allows us to point out the complementarity between theory and simulation and to highlight the convenience but also the limits of the asymptotic and linear theories. We show that a rich spectrum of IGWs is excited by the convection, representing about 0.4% of the total solar luminosity. We study the spatial and temporal properties of this spectrum, the effect of thermal damping and nonlinear interactions between waves. We give quantitative results about the modes frequencies, evolution with time and rotational splitting and we discuss the amplitude of IGWs considering different regimes of parameters. This work points out the importance of high performance simulation for its complementarity with observation and theory. It opens a large field of investigation concerning IGWs propagating nonlinearly in 3D spherical structures. The extension of this work to other types of stars, with different masses, structures and rotation rates will lead to a deeper and more accurate comprehension of IGWs in stars.

preprint2013arXiv

Cycling Dynamo in a Young Sun: Grand Minima and Equatorward Propagation

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.

preprint2011arXiv

Assimilating data into an αω dynamo model of the sun: A variational approach

We have developed a variational data assimilation technique for the Sun using a toy αΩ dynamo model. The purpose of this work is to apply modern data assimilation techniques to solar data using a physically based model. This work represents the first step toward a complete variational model of solar magnetism. We derive the adjoint αΩ dynamo code and use a minimization procedure to invert the spatial dependence of key physical ingredients of the model. We find that the variational technique is very powerful and leads to encouraging results that will be applied to a more realistic model of the solar dynamo. We also note that the continuous adjoint mean-field induction equation is derived in the appendix.

preprint2011arXiv

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.

preprint2011arXiv

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én surface. Wind sources and braking torque application zones also vary accordingly.

preprint2011arXiv

Effects of turbulent pumping on stellar activity cycles

Stellar magnetic activity of solar like stars is thought to be due to an internal dynamo. While the Sun has been the subject of intense research for refining dynamo models, observations of magnetic cyclic activity in solar type stars have become more and more available, opening a new path to understand the underlying physics behind stellar cycles. For instance, it is key to understand how stellar rotation rate influences magnetic cycle period $\Pcyc$. Recent numerical simulations of advection-dominated Babcock Leighton models have demonstrated that it is difficult to explain this observed trend given a) the strong influence of the cycle period to the meridional circulation amplitude and b) the fact that 3D models indicate that meridional flows become weaker as the rotation rate increases. In this paper, we introduce the turbulent pumping mechanism as another advective process capable also of transporting the magnetic fields. We found that this model is now able to reproduce the observations under the assumption that this effect increases as $Ω^2$. The turbulent pumping becomes indeed another major player able to circumvent the meridional circulation. However, for high rotation rates ($Ω\simeq 5 Ω_\odot$), its effects dominate those of the meridional circulation, entering a new class of regime dominated by the advection of turbulent pumping and thus leading to a cyclic activity qualitatively different from that of the sun.

preprint2011arXiv

Hunting for giant cells in deep stellar convective zones using wavelet analysis

We study the influence of stratification on stellar turbulent convection near the stellar surface and in depth by carrying out 3D high resolution hydrodynamic simulations with the ASH code. Four simulations with different radial density contrast corresponding to different aspect ratio for the same underlying 4 Myr $0.7 M_\odot$ pre-main sequence star model are thus performed. We highlight the existence of giant cells which are embedded in the complex surface convective patterns using a wavelet and time correlation analysis. Next, we study their properties such as their lifetime, aspect ratio and spatial extension in the different models both in latitude and depth according to the density contrast. We find that these giants cells have a lifetime larger than the stellar period with a typical longitudinal width of 490 Mm and a latitudinal extension increasing with the radial density contrast overpassing 50\degree of latitude for the thickest convective zone. Their rotation rate is much larger than the local differential rotation rate increasing also with the radial density contrast. However, their spatial coherence as a function of depth decrease with the density contrast due to the stronger shear present in these more stratified cases.

preprint2010arXiv

Global magnetic cycles in rapidly rotating younger suns

Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D MHD simulations carried out at three times the current solar rotation rate reveal hydromagnetic dynamo action that yields wreaths of strong toroidal magnetic field at low latitudes, often with opposite polarity in the two hemispheres. Our recent simulations have explored behavior in systems with considerably lower diffusivities, achieved with sub-grid scale models including a dynamic Smagorinsky treatment of unresolved turbulence. The lower diffusion promotes the generation of magnetic wreaths that undergo prominent temporal variations in field strength, exhibiting global magnetic cycles that involve polarity reversals. In our least diffusive simulation, we find that magnetic buoyancy coupled with advection by convective giant cells can lead to the rise of coherent loops of magnetic field toward the top of the simulated domain.

preprint2010arXiv

Global-scale Magnetism (and Cycles) in Dynamo Simulations of Stellar Convection Zones

Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. We explore one case rotates five times faster than the Sun in detail. This dynamo simulation, called case D5, has repeated quasi-cyclic reversals of global-scale polarity. We compare this case D5 to the broader family of simulations we have been able to explore and discuss how future simulations and observations can advance our understanding of stellar dynamos and magnetism.

preprint2010arXiv

Magnetic Cycles and Meridional Circulation in Global Models of Solar Convection

We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and systematic variability on time scales of 6-15 years. The second development concerns the maintenance of the meridional circulation. Recent high-resolution simulations have captured the subtle nonlinear dynamical balances with more fidelity than previous, more laminar models, yielding more coherent circulation patterns. These patterns are dominated by a single cell in each hemisphere, with poleward and equatorward flow in the upper and lower convection zone respectively. We briefly address the implications of and future of these modeling efforts.

preprint2009arXiv

Three dimensional non-linear evolution of a magnetic flux tube in a spherical shell: Influence of turbulent convection and associated mean flows

We present the first 3D MHD study in spherical geometry of the non-linear dynamical evolution of magnetic flux tubes in a turbulent rotating convection zone. We study numerically the rise of magnetic toroidal flux ropes from the base of a modelled convection zone up to the top of our computational domain where bipolar patches are formed. We compare the dynamical behaviour of flux tubes in a fully convective shell possessing self-consistently generated mean flows such as meridional circulation and differential rotation, with reference calculations done in a quiet isentropic zone. We find that two parameters influence the tubes during their rise through the convection zone: the initial field strength and amount of twist, thus confirming previous findings in Cartesian geometry. Further, when the tube is sufficiently strong with respect to the equipartition field, it rises almost radially independently of the initial latitude (either low or high). By contrast, weaker field cases indicate that downflows and upflows control the rising velocity of particular regions of the rope and could in principle favour the emergence of flux through Omega-loop structures. For these latter cases, we focus on the orientation of bipolar patches and find that sufficiently arched structures are able to create bipolar regions with a predominantly East-West orientation. Meridional flow seems to determine the trajectory of the magnetic rope when the field strength has been significantly reduced near the top of the domain. Finally differential rotation makes it more difficult for tubes introduced at low latitudes to reach the top of the domain.

preprint2006arXiv

Dynamo Action in the Solar Convection Zone and Tachocline: Pumping and Organization of Toroidal Fields

We present the first results from three-dimensional spherical shell simulations of magnetic dynamo action realized by turbulent convection penetrating downward into a tachocline of rotational shear. This permits us to assess several dynamical elements believed to be crucial to the operation of the solar global dynamo, variously involving differential rotation resulting from convection, magnetic pumping, and amplification of fields by stretching within the tachocline. The simulations reveal that strong axisymmetric toroidal magnetic fields (about 3000 G in strength) are realized within the lower stable layer, unlike in the convection zone where fluctuating fields are predominant. The toroidal fields in the stable layer possess a striking persistent antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. The associated mean poloidal magnetic fields there have a clear dipolar geometry, but we have not yet observed any distinctive reversals or latitudinal propagation. The presence of these deep magnetic fields appears to stabilize the sense of mean fields produced by vigorous dynamo action in the bulk of the convection zone.