Source author record

Juri Toomre

Juri Toomre 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

26works
6topics
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

26 published item(s)

preprint2024arXiv

Confinement of the Solar Tachocline by a Non-Axisymmetric Dynamo

We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find three solution regimes, corresponding to weak, medium, and strong dynamo magnetic field strengths. In the weak-field regime, the large-scale magnetic field is mostly axisymmetric with regular, periodic polarity reversals (reminiscent of the observed solar cycle), but fails to create a confined tachocline. In the strong-field regime, the large-scale field is mostly non-axisymmetric with irregular, quasi-periodic polarity reversals, and creates a confined tachocline. In the medium-field regime, the large-scale field resembles a strong-field dynamo for extended intervals, but intermittently weakens to allow temporary epochs of strong differential rotation. In all regimes, the amplitude of poloidal field strength in the RZ is very well explained by skin-depth arguments, wherein the oscillating field that gives rise to the skin depth (in the medium- and strong-field cases) is a non-axisymmetric field structure rotating with respect to the RZ. These new simulations reaffirm that tachocline confinement by the solar dynamo (the so-called fast magnetic confinement scenario) is possible, but suggest a new picture in which non-axisymmetric field components rotating with respect to the RZ play the primary role, instead of the regularly reversing axisymmetic field associated with the 22-year cycle.

preprint2022arXiv

Longitudinally Modulated Dynamo Action in Simulated M-Dwarf Stars

M-dwarf stars are well known for the intense magnetic activity that many of them exhibit. In cool stars with near-surface convection zones, this magnetic activity is thought to be driven largely by the interplay of convection and the large scale differential rotation and circulations it establishes. The highly nonlinear nature of these flows yields a fascinatingly sensitive and diverse parameter space, with a wide range of possible dynamics. We report here on a set of three global MHD simulations of rapidly rotating M2 (0.4 $M_\odot$) stars. Each of these three models established nests of vigorous convection that were highly modulated in longitude at low latitudes. Slight differences in their magnetic parameters led each model to disparate dynamo states, but the effect of the convective nest was a unifying feature. In each case, the action of longitudinally modulated convection led to localized (and in one case, global) reversals of the toroidal magnetic field, as well as the formation of an active longitude, with enhanced poloidal field amplitudes and flux emergence.

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.

preprint2020arXiv

Exploring Bistability in the Cycles of the Solar Dynamo through Global Simulations

The calling card of solar magnetism is the sunspot cycle, during which sunspots regularly reverse their polarity sense every 11 years. However, a number of more complicated time-dependent behaviors have also been identified. In particular, there are temporal modulations associated with active longitudes and hemispheric asymmetry, when sunspots appear at certain solar longitudes or else in one hemisphere preferentially. So far, a direct link between between this asymmetric temporal behavior and the underlying solar dynamo has remained elusive. In this work, we present results from global, 3D magnetohydrodynamic (MHD) simulations, which for the first time display both behavior reminiscent of the sunspot cycle (regular polarity reversals and equatorward migration of internal magnetic field) and asymmetric, irregular behavior that in the simulations we interpret as active longitudes and hemispheric asymmetry. The simulations are thus bistable, in that the turbulent convection can stably support two distinct flavors of magnetism at different times, in superposition, or with smooth transitions from one state to the other. We discuss this new family of dynamo models in the context of the extensive observations of the Sun's surface magnetic field with the Solar and Heliospheric Observatory (SOHO) and the Solar Dynamics Observatory (SDO), as well as earlier observations of sunspot number and synoptic maps. We suggest that the solar dynamo itself may be bistable in nature, exhibiting two types of temporal behavior in the magnetic field.

preprint2020arXiv

Probing the Influence of a Tachocline in Simulated M-Dwarf Dynamos

In recent years, observers have found that the fraction of M-stars demonstrating significant magnetic activity transitions sharply from roughly $10\%$ for main-sequence stars earlier (more massive) than spectral type M3.5 (0.35 M$_\odot$) to nearly $90\%$ for stars later than M3.5. The latter are typically rotating quite rapidly, suggesting differing spin-down histories. Tantalizingly, it is also later than M3.5 at which main-sequence stars become fully convective, and may no longer contain a tachocline, a layer of rotational shear revealed by helioseismology to separate the radiative zone (RZ) and convection zone (CZ). We turn here to the more massive M-stars to study the impact such a layer may have on their internal dynamics. Using the spherical 3D MHD simulation code Rayleigh, we compare the properties of convective dynamos generated within quickly rotating (1, 2, and 4 $Ω_\odot$) M2 (0.4 M$_\odot$) stars, with the computational domain either terminating at the base of the convection zone or permitting overshoot into the underlying stable region. We find that a tachocline is not necessary for the organization of strong toroidal wreaths of magnetism in these stars, though its presence can increase the coupling of mean field amplitudes to the stellar rotation rate. Additionally, in stars that undergo periodic cycles, we find that the presence of a tachocline tends to make the cycles both longer and more regular than they would have otherwise been. Finally, we find that the tachocline helps to enhance the surface poloidal fields and organize them into larger spatial scales, both of which provide favorable conditions for more rapid angular momentum loss through a magnetized stellar wind.

preprint2020arXiv

Revisiting the Sun's Strong Differential Rotation along Radial Lines

Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition.

preprint2015arXiv

Convective Dynamo Simulation with a Grand Minimum

The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. It is found that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals occurring roughly every 6.2 years. The magnetic energy also rises and falls with a regular period. The magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the processes leading to polarity reversals are more complex, appearing to arise from the interaction of convection with the mean toroidal fields. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation, and potentially to a nonlinear dynamo wave. This simulation also enters a grand minimum lasting roughly 20~years, after which the dynamo recovers its regular polarity cycles.

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.

preprint2015arXiv

Helioseismic Imaging of Fast Convective Flows Throughout the Near-Surface Shear Layer

Using a new implementation of ring-diagram helioseismology, we ascertain the strength and spatial scale of convective flows throughout the near-surface shear layer. Our ring-diagram technique employs highly overlapped analysis regions and an efficient method of 3D inversion to measure convective motions with a resolution that ranges from $3 \ \mathrm{Mm}$ at the surface to $80 \ \mathrm{Mm}$ at the base of the layer. We find the rms horizontal flow speed to peak at $427 \ \mathrm{m \ s^{-1}}$ at the photosphere and fall to a minimum of $124 \ \mathrm{m \ s^{-1}}$ between $20 \ \mathrm{Mm}$ and $30 \ \mathrm{Mm}$. From the velocity amplitude and the dominant horizontal scales seen at each depth, we infer the level of rotational influence on convection to be low near the surface, but transition to a significant level at the base of the near-surface shear layer with a Rossby number varying between 2.2 to as low as 0.1.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

preprint2014arXiv

Multi-Ridge Fitting for Ring-Diagram Helioseismology

Inferences of sub-surface flow velocities using local domain ring-diagram helioseismology depend on measuring the frequency splittings of oscillation modes seen in acoustic power spectra. Current methods for making these measurements utilize maximum-likelihood fitting techniques to match a model of modal power to the spectra. The model typically describes a single oscillation mode, and each mode in a given power spectrum is fit independently. We present a new method that produces measurements with greater reliability and accuracy by fitting multiple modes simultaneously. We demonstrate how this method permits measurements of sub-surface flows deeper into the Sun while providing higher uniformity in data coverage and velocity response closer to the limb of the solar disk. While the previous fitting method performs better for some measurements of low-phase-speed modes, we find this new method to be particularly useful for high phase-speed modes and small spatial areas.

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.

preprint2012arXiv

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three times the current solar rate. These buoyant loops originate within sections of the magnetic wreaths in which turbulent flows amplify the fields to much larger values than is possible through laminar processes. These amplified portions can rise through the convective layer by a combination of magnetic buoyancy and advection by convective giant cells, forming buoyant loops. We measure statistical trends in the polarity, twist, and tilt of these loops. Loops are shown to preferentially arise in longitudinal patches somewhat reminiscent of active longitudes in the Sun, although broader in extent. We show that the strength of the axisymmetric toroidal field is not a good predictor of the production rate for buoyant loops or the amount of magnetic flux in the loops that are produced.

preprint2012arXiv

Effects of Granulation upon Larger-Scale Convection

We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergranular scales are realized, along with a differential rotation that becomes increasingly solar-like.

preprint2012arXiv

Magnetic Wreaths and Cycles in Convective Dynamos

Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried out at a nominal rotation rate of three times the solar value (3$Ω_\odot$), but similar dynamics may also apply to the Sun. Our previous simulations demonstrated that convective dynamos can build persistent toroidal flux structures (magnetic wreaths) in the midst of a turbulent convection zone and that high rotation rates promote the cyclic reversal of these wreaths. Here we demonstrate that magnetic cycles can also be achieved by reducing the diffusion, thus increasing the Reynolds and magnetic Reynolds numbers. In these more turbulent models, diffusive processes no longer play a significant role in the key dynamical balances that establish and maintain the differential rotation and magnetic wreaths. Magnetic reversals are attributed to an imbalance in the poloidal magnetic induction by convective motions that is stabilized at higher diffusion levels. Additionally, the enhanced levels of turbulence lead to greater intermittency in the toroidal magnetic wreaths, promoting the generation of buoyant magnetic loops that rise from the deep interior to the upper regions of our simulated domain. The implications of such turbulence-induced magnetic buoyancy for solar and stellar flux emergence are also discussed.

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

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 $Ω$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $α$-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]

preprint2010arXiv

Convection and dynamo action in B stars

Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields arise within the turbulence of the core and penetrate into the stably stratified radiative zone. These fields exhibit complex, time-dependent behavior including reversals in magnetic polarity and shifts between which hemisphere dominates the total magnetic energy.

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

Global-scale wreath-building dynamos in stellar convection zones

When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with 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. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.

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.

preprint2010arXiv

Modeling the Near-Surface Shear Layer: Diffusion Schemes Studied With CSS

As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.

preprint2010arXiv

Persistent Magnetic Wreaths in a Rapidly Rotating Sun

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field $α$-effect explains the regeneration of poloidal field.

preprint2009arXiv

Subsurface Circulations within Active Regions

Using high-resolution ring analysis (HRRA) we deduce subsurface flows within magnetic active regions and within quiet sun. With this procedure we are capable of measuring flows with a horizontal spatial resolution of 2 degrees in heliographic angle (or roughly 20 Mm). From the resulting flow fields we deduce mean inflow rates into active regions, mean circulation speeds around active regions, and probability density functions (PDFs) of properties of the flow field. These analyses indicate that active regions have a zonal velocity that exceeds that of quiet sun at the same latitude by 20 m/s, yet active regions advect poleward at the same rate as quiet sun. We also find that almost all active regions possess a mean inflow (20-30 m/s) and a cyclonic circulation (~5 m/s) at their peripheries, whereas their cores, where the sunspots are located, are zones of strong anticyclonic outflow (~50 m/s). From the PDFs, we find that active regions modify the structure of convection with a scale greater than that of supergranulation. Instead of possessing an asymmetry between inflows and outflows (with a larger percentage of the surface occupied by outflows), as is seen in quiet sun, active regions possess symmetric distributions.

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.