-
Small-scale and large-scale dynamos in global convection simulations of solar-like stars
Authors:
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Matthias Rheinhard,
Mariangela Viviani,
Ameya Prabhu
Abstract:
It has been recently shown that a small-scale dynamo (SSD) instability could be possible in solar-like low magnetic Prandtl number Pm plasmas. It has been proposed that the presence of SSD can potentially have a significant impact on the dynamics of the large-scale dynamo (LSD) in the stellar convection zones. Studying these two dynamos, SSD and LSD, together in a global magnetoconvection model re…
▽ More
It has been recently shown that a small-scale dynamo (SSD) instability could be possible in solar-like low magnetic Prandtl number Pm plasmas. It has been proposed that the presence of SSD can potentially have a significant impact on the dynamics of the large-scale dynamo (LSD) in the stellar convection zones. Studying these two dynamos, SSD and LSD, together in a global magnetoconvection model requires high-resolution simulations and large amounts of computational resources. Starting from a well-studied global convective dynamo model that produces cyclic magnetic fields, we systematically increased the resolution and lowered the diffusivities to enter the regime of Reynolds numbers that allow for the excitation of SSD on top of the LSD. We studied how the properties of convection, generated differential rotation profiles, and LSD solutions change with the presence of SSD. We performed convective dynamo simulations in a spherical wedge with the Pencil Code. The resolutions of the models were increased in 4 steps by a total factor of 16 to achieve maximal fluid and magnetic Reynolds numbers of over 500. We found that the differential rotation is strongly quenched by the presence of the LSD and SSD. Even though the small-scale magnetic field only mildly decreases increasing Re, the large-scale field strength decreases significantly. We do not find the SSD dynamo significantly quenching the convective flows as claimed recently by other authors; in contrast, the convective flows first grow and then saturate for increasing Re. Furthermore, the angular momentum transport is highly affected by the presence of small-scale magnetic fields, which are mostly generated by LSD. These fields not only change the Reynolds stresses, but also generate dynamically important Maxwell stresses. The LSD evolution in terms of its pattern and field distribution is rather independent of the increase in Rm.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Data-driven simulations of magnetic field evolution in Active Region 11429: Magneto-friction method using PENCIL code
Authors:
P. Vemareddy,
Jörn Warnecke,
Ph. A. Bourdin
Abstract:
Coronal magnetic fields evolve quasi statically over long time scales and dynamically over short time scales. As of now there exists no regular measurements of coronal magnetic fields, and therefore generating the coronal magnetic field evolution using the observations of the magnetic field at the photosphere is of fundamental requirement to understand the origin of the transient phenomena from th…
▽ More
Coronal magnetic fields evolve quasi statically over long time scales and dynamically over short time scales. As of now there exists no regular measurements of coronal magnetic fields, and therefore generating the coronal magnetic field evolution using the observations of the magnetic field at the photosphere is of fundamental requirement to understand the origin of the transient phenomena from the solar active regions. Using the magnetofriction (MF) approach, we aim to simulate the coronal field evolution in the solar active region 11429. The MF method is implemented in open source \PC along with a driver module to drive the initial field with different boundary conditions prescribed from observed vector magnetic fields at the photosphere. In order to work with vector potential and the observations, we prescribe three types of bottom boundary drivers with varying free-magnetic energy. The MF simulation reproduces the magnetic structure, which better matches to the sigmoidal morphology exhibited by AIA images at the pre-eruptive time. We found that the already sheared field further driven by the sheared magnetic field, will maintain and further build the highly sheared coronal magnetic configuration, as seen in AR 11429. Data-driven MF simulation is a viable tool to generate the coronal magnetic field evolution, capturing the formation of the twisted flux rope and its eruption.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
From convective stellar dynamo simulations to Zeeman-Doppler images
Authors:
Thomas Hackman,
Oleg Kochukhov,
Mariangela Viviani,
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Jyri J. Lehtinen
Abstract:
Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preser…
▽ More
Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preserved when applying the ZDI method on simulated data, i.e. how accurately the field topology is preserved and to what extent stellar parameters influence the reconstruction. We use published magnetic field data from direct numerical MHD simulations. These have variable rotation rates, and hence represent different levels of activity, of an otherwise Sun-like setup. Our ZDI reconstruction is based on spherical harmonics expansion. By comparing the original values to those of the reconstructed images, we study the ability to reconstruct the surface magnetic field in terms of various characteristics of the field. The main large-scale features are reasonably well recovered, but the strength of the recovered magnetic field is just a fraction of the original input. The quality of the reconstruction shows clear correlations with the data quality. Furthermore, there are some spurious dependencies between stellar parameters and the characteristics of the field. Our study uncovers some limits of ZDI. Firstly, the recovered field strength will generally be lower than the "real" value as smaller structures with opposite polarities will be blurred in the inversion. Secondly, the axisymmetry is overestimated. The poloidality vs. toroidality is better recovered. The reconstruction works better for a stronger field and faster rotation velocity. Still, the ZDI method works surprisingly well even for a weaker field and slow rotation, provided the data has a high signal-to-noise and good rotation phase coverage.
△ Less
Submitted 22 January, 2024; v1 submitted 13 June, 2023;
originally announced June 2023.
-
Numerical evidence for a small-scale dynamo approaching solar magnetic Prandtl numbers
Authors:
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Frederick A. Gent,
Matthias Rheinhardt
Abstract:
Magnetic fields on small scales are ubiquitous in the universe. Though they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers ($Pr_M$) such as are present in the Sun and othe…
▽ More
Magnetic fields on small scales are ubiquitous in the universe. Though they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers ($Pr_M$) such as are present in the Sun and other cool stars. We have performed high-resolution simulations of isothermal forced turbulence employing the lowest $Pr_M$ values so far achieved. Contrary to earlier findings, the SSD turns out to be not only possible for $Pr_M$ down to 0.0031, but even becomes increasingly easier to excite for $Pr_M$ below $\simeq\,$0.05. We relate this behaviour to the known hydrodynamic phenomenon referred to as the bottleneck effect. Extrapolating our results to solar values of $Pr_M$ indicates that an SSD would be possible under such conditions.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
Simulations of solar and stellar dynamos and their theoretical interpretation
Authors:
Petri J. Käpylä,
Matthew K. Browning,
Allan Sacha Brun,
Gustavo Guerrero,
Jörn Warnecke
Abstract:
We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scal…
▽ More
We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
Constraining Global Solar Models through Helioseismic Analysis
Authors:
Andrey M. Stejko,
Alexander G. Kosovichev,
Nicholas A. Featherstone,
Gustavo Guerrero,
Bradley W. Hindman,
Loren I. Matilsky,
Jörn Warnecke
Abstract:
Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this…
▽ More
Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise.
△ Less
Submitted 11 April, 2022;
originally announced April 2022.
-
Investigating global convective dynamos with mean-field models: full spectrum of turbulent effects required
Authors:
Jörn Warnecke,
Matthias Rheinhardt,
Mariangela Viviani,
Frederick Gent,
Simo Tuomisto,
Maarit J. Käpylä
Abstract:
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport…
▽ More
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport coefficients, based on the MF framework, from an exemplary GCD simulation using the test-field method. These coefficients are then used as an input into an MF model. We find a good agreement between the MF and GCD solutions, which validates our theoretical approach. This agreement requires all turbulent effects to be included, even those which have been regarded as unimportant so far. Our results suggest that simple dynamo models, as are commonly used in the solar and stellar community, relying on very few, precisely fine-tuned turbulent effects, may not be representative of the full dynamics of dynamos in global convective simulations and astronomical objects.
△ Less
Submitted 5 October, 2021; v1 submitted 17 May, 2021;
originally announced May 2021.
-
Hunting down the cause of solar magnetism
Authors:
M. Viviani,
A. Prabhu,
J. Warnecke,
L. Duarte,
J. Pekkilä,
M. Rheinhardt,
M. J. Käpylä
Abstract:
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions…
▽ More
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Stellar X-rays and magnetic activity in 3D MHD coronal models
Authors:
J. Zhuleku,
J. Warnecke,
H. Peter
Abstract:
Observations suggest a power-law relation between the coronal emission in X-rays, $L_{\rm{X}}$, and the total (unsigned) magnetic flux at the stellar surface, $Φ$. The physics basis for this relation is poorly understood. We use three-dimensional (3D) magnetohydrodynamics (MHD) numerical models of the coronae above active regions, that is, strong concentrations of magnetic field, to investigate th…
▽ More
Observations suggest a power-law relation between the coronal emission in X-rays, $L_{\rm{X}}$, and the total (unsigned) magnetic flux at the stellar surface, $Φ$. The physics basis for this relation is poorly understood. We use three-dimensional (3D) magnetohydrodynamics (MHD) numerical models of the coronae above active regions, that is, strong concentrations of magnetic field, to investigate the $L_{\rm{X}}$ versus $Φ$ relation and illustrate this relation with an analytical model based on simple well-established scaling relations. In the 3D MHD model horizontal (convective) motions near the surface induce currents in the coronal magnetic field that are dissipated and heat the plasma. This self-consistently creates a corona with a temperature of 1 MK. We run a series of models that differ in terms of the (unsigned) magnetic flux at the surface by changing the (peak) magnetic field strength while keeping all other parameters fixed. In the 3D MHD models we find that the energy input into the corona, characterized by either the Poynting flux or the total volumetric heating, scales roughly quadratically with the unsigned surface flux $Φ$. This is expected from heating through field-line braiding. Our central result is the nonlinear scaling of the X-ray emission as $L_{\rm{X}}\propto Φ^{3.44}$. This scaling is slightly steeper than found in recent observations that give power-law indices of up to only 2 or 3. Assuming that on a real star, not only the peak magnetic field strength in the active regions changes but also their number (or surface filling factor), our results are consistent with observations. Our model provides indications of what causes the steep increase in X-ray luminosity by four orders of magnitude from solar-type activity to fast rotating active stars.
△ Less
Submitted 5 July, 2021; v1 submitted 1 February, 2021;
originally announced February 2021.
-
The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
Authors:
A. Brandenburg,
A. Johansen,
P. A. Bourdin,
W. Dobler,
W. Lyra,
M. Rheinhardt,
S. Bingert,
N. E. L. Haugen,
A. Mee,
F. Gent,
N. Babkovskaia,
C. -C. Yang,
T. Heinemann,
B. Dintrans,
D. Mitra,
S. Candelaresi,
J. Warnecke,
P. J. Käpylä,
A. Schreiber,
P. Chatterjee,
M. J. Käpylä,
X. -Y. Li,
J. Krüger,
J. R. Aarnes,
G. R. Sarson
, et al. (12 additional authors not shown)
Abstract:
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the…
▽ More
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Stellar coronal X-ray emission and surface magnetic flux
Authors:
J. Zhuleku,
J. Warnecke,
H. Peter
Abstract:
Observations show that the coronal X-ray emission of the Sun and other stars depends on the surface magnetic field. Using power-law scaling relations between different physical parameters, we build an analytical model to connect the observed X-ray emission to the magnetic flux. The basis for our model are the scaling laws of Rosner, Tucker \& Vaiana (RTV) that connect the temperature and pressure…
▽ More
Observations show that the coronal X-ray emission of the Sun and other stars depends on the surface magnetic field. Using power-law scaling relations between different physical parameters, we build an analytical model to connect the observed X-ray emission to the magnetic flux. The basis for our model are the scaling laws of Rosner, Tucker \& Vaiana (RTV) that connect the temperature and pressure of a coronal loop to its length and energy input. To estimate the energy flux into the upper atmosphere, we use scalings derived for different heating mechanisms, e.g. for field-line braiding or Alfven-wave heating. We supplement this by observed relations between active region size and magnetic flux and derive scalings of how X-ray emissivity depends on temperature. Based on our analytical model, we find a power-law dependence of the X-ray emission on the magnetic flux, $L_{\rm X}\propto Φ^m$, with a power-law index $m$ being in the range from about 1 to 2. This finding is consistent with a wide range of observations, from individual features on the Sun, e.g. bright points or active regions, to stars of different types and varying levels of activity. The power-law index $m$ depends on the choice of the heating mechanism, and our results slightly favour the braiding and nanoflare scenarios over Alfvén wave heating. In addition, the choice of instrument will have an impact on the power-law index $m$, which is because of the sensitivity of the observed wavelength region to the temperature of the coronal plasma. Overall, our simple analytical model based on the RTV scaling laws gives a good representation of the observed X-ray emission. This underlines that we might be able to understand stellar coronal activity though a collection of basic building blocks, i.e. loops, that we can study in spatially resolved detail on the Sun.
△ Less
Submitted 24 June, 2020;
originally announced June 2020.
-
Shapes of stellar activity cycles
Authors:
T. Willamo,
T. Hackman,
J. J. Lehtinen,
M. J. Käpylä,
N. Olspert,
M. Viviani,
J. Warnecke
Abstract:
Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other ste…
▽ More
Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one.
△ Less
Submitted 29 April, 2020; v1 submitted 11 February, 2020;
originally announced February 2020.
-
On the influence of magnetic helicity on X-rays emission of solar and stellar coronae
Authors:
Jörn Warnecke,
Hardi Peter
Abstract:
Observation of solar-like stars show a clear relation between X-ray emission and their rotation. Higher stellar rotation can lead to a larger magnetic helicity production in stars. We aim to understand the relation between magnetic helicity on the surface of a star to their coronal X-ray emission. We use 3D MHD simulations to model the corona of the solar-like stars. We take an observed magnetogra…
▽ More
Observation of solar-like stars show a clear relation between X-ray emission and their rotation. Higher stellar rotation can lead to a larger magnetic helicity production in stars. We aim to understand the relation between magnetic helicity on the surface of a star to their coronal X-ray emission. We use 3D MHD simulations to model the corona of the solar-like stars. We take an observed magnetogram as in photospheric activity input, and inject different values of magnetic helicity. We use synthesis emission to calculate the X-ray emission flux of each simulation and investigate how this scales with injected magnetic helicity. We find that for larger injected magnetic helicities an increase in temperature and an increase in X-ray emission. The X-ray emission scaled cubicly with the injected helicity. We can related this to increase of horizontal magnetic field and therefore higher Poynting flux at the coronal base. Using typical scaling of magnetic helicity production with stellar rotation, we can explain the increase of X-ray emission with rotation only by an increase of magnetic helicity at the surface of a star.
△ Less
Submitted 16 October, 2019; v1 submitted 15 October, 2019;
originally announced October 2019.
-
Rotational dependence of turbulent transport coefficients in global convective dynamo simulations of solar-like stars
Authors:
J. Warnecke,
M. J. Käpylä
Abstract:
For moderate and slow rotation, magnetic activity of solar-like stars is observed to strongly depend on rotation. These observations do not yet have a solid explanation in terms of dynamo theory. We aim to find such an explanation by numerically investigated the rotational dependency of dynamo drivers in solar-like stars. We ran semi-global convection simulations of stars with rotation rates from…
▽ More
For moderate and slow rotation, magnetic activity of solar-like stars is observed to strongly depend on rotation. These observations do not yet have a solid explanation in terms of dynamo theory. We aim to find such an explanation by numerically investigated the rotational dependency of dynamo drivers in solar-like stars. We ran semi-global convection simulations of stars with rotation rates from 0 to 30 times the solar value, corresponding to Coriolis numbers, Co, of 0 to 110. We measured the turbulent transport coefficients describing the magnetic field evolution with the help of the test-field method, and compared with the dynamo effect arising from the differential rotation. The trace of the $α$ tensor increases for moderate rotation rates with Co$^{0.5}$ and levels off for rapid rotation. This behavior agrees with the kinetic $α$, if one considers the decrease of the convective scale with increasing rotation. The $α$ tensor becomes highly anisotropic for Co$\gtrsim 1$. Furthermore, $α_{rr}$ dominates for moderate rotation (1<Co<10), and $α_{φφ}$ for rapid rotation (Co $\gtrsim 10$). The turbulent pumping effect is dominating the meridional transport of the magnetic field. Taking all dynamo effects into account, we find three distinct regimes. For slow rotation, the $α$ and Rädler effects are dominating in the presence of anti-solar differential rotation. For moderate rotation, $α$ and $Ω$ effects are dominant, indicative of $αΩ$ or $α^2Ω$ dynamos in operation, producing equatorward-migrating dynamo waves with a qualitatively solar-like rotation profile. For rapid rotation, an $α^2$ mechanism, with an influence from the Rädler effect appears to be the most probable driver of the dynamo. Our study reveals the presence of a large variety of dynamo effects beyond the classical $αΩ$ mechanism.
△ Less
Submitted 22 July, 2020; v1 submitted 15 October, 2019;
originally announced October 2019.
-
Data-driven model of the solar corona above an active region
Authors:
Jörn Warnecke,
Hardi Peter
Abstract:
We aim to reproduce the structure of the corona above a solar active region as seen in the extreme ultraviolet (EUV) using a three-dimensional magnetohydrodynamic (3D MHD) model. The 3D MHD data-driven model solves the induction equation and the mass, momentum, and energy balance. To drive the system, we feed the observed evolution of the magnetic field in the photosphere of the active region AR 1…
▽ More
We aim to reproduce the structure of the corona above a solar active region as seen in the extreme ultraviolet (EUV) using a three-dimensional magnetohydrodynamic (3D MHD) model. The 3D MHD data-driven model solves the induction equation and the mass, momentum, and energy balance. To drive the system, we feed the observed evolution of the magnetic field in the photosphere of the active region AR 12139 into the bottom boundary. This creates a hot corona above the cool photosphere in a self-consistent way. We synthesize the coronal EUV emission from the densities and temperatures in the model and compare this to the actual coronal observations. We are able to reproduce the overall appearance and key features of the corona in this active region on a qualitative level. The model shows long loops, fan loops, compact loops, and diffuse emission forming at the same locations and at similar times as in the observation. Furthermore, the low-intensity contrast of the model loops in EUV matches the observations. In our model the energy input into the corona is similar as in the scenarios of fieldline-braiding or flux-tube tectonics, that is, energy is transported to the corona through the driving of the vertical magnetic field by horizontal photospheric motions. The success of our model shows the central role that this process plays for the structure, dynamics, and heating of the corona.
△ Less
Submitted 29 April, 2019; v1 submitted 1 March, 2019;
originally announced March 2019.
-
Stellar dynamos in the transition regime: multiple dynamo modes and anti-solar differential rotation
Authors:
M. Viviani,
M. J. Käpylä,
J. Warnecke,
P. J. Käpylä,
M. Rheinhardt
Abstract:
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions…
▽ More
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $αΩ$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $α$ effect is comparable or even larger than the $Ω$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $α^2Ω$ or $α^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping.
△ Less
Submitted 22 November, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
Non-Fourier description of heat flux evolution in 3D MHD simulations of the solar corona
Authors:
Jörn Warnecke,
Sven Bingert
Abstract:
The hot loop structures in the solar corona can be well modeled by three dimensional magnetohydrodynamic simulations, where the corona is heated by field line braiding driven at the photosphere. To be able to reproduce the emission comparable to observations, one has to use realistic values for the Spitzer heat conductivity, which puts a large constraint on the time step of these simulations and m…
▽ More
The hot loop structures in the solar corona can be well modeled by three dimensional magnetohydrodynamic simulations, where the corona is heated by field line braiding driven at the photosphere. To be able to reproduce the emission comparable to observations, one has to use realistic values for the Spitzer heat conductivity, which puts a large constraint on the time step of these simulations and make them therefore computationally expensive. Here, we present a non-Fourier description of the heat flux evolution, which allow us to speed up the simulations significantly. Together with the semi-relativistic Boris correction, we are able to limit the time step constraint of the Alfvén speed and speed up the simulations even further. We discuss the implementation of these two methods to the Pencil Code and present their implications on the time step, and the temperature structures, the ohmic heating rate and the emission in simulations of the solar corona. Using a non-Fourier description of the heat flux evolution together with the Boris correction, we can increase the time step of the simulation significantly without moving far away from the reference solution. However, for values of the Alfvén speed limit of 3 000 km/s and below, the simulation moves away from the reference solution und produces much higher temperatures and much structures with stronger emission.
△ Less
Submitted 27 November, 2019; v1 submitted 5 November, 2018;
originally announced November 2018.
-
Chromospheric activity catalogue of 4454 cool stars. Questioning the active branch of stellar activity cycles
Authors:
S. Boro Saikia,
C. J. Marvin,
S. V. Jeffers,
A. Reiners,
R. Cameron,
S. C. Marsden,
P. Petit,
J. Warnecke,
A. P. Yadav
Abstract:
Chromospheric activity monitoring of a wide range of cool stars can provide valuable information on stellar magnetic activity and its dependence on fundamental stellar parameters such as effective temperature and rotation.We compile a chromospheric activity catalogue of 4454 cool stars from a combination of archival HARPS spectra and multiple other surveys, including the Mount Wilson data that hav…
▽ More
Chromospheric activity monitoring of a wide range of cool stars can provide valuable information on stellar magnetic activity and its dependence on fundamental stellar parameters such as effective temperature and rotation.We compile a chromospheric activity catalogue of 4454 cool stars from a combination of archival HARPS spectra and multiple other surveys, including the Mount Wilson data that have recently been released by the NSO. We explore the variation in chromospheric activity of cool stars along the main sequence for stars with different effective temperatures. Additionally, we also perform an activity-cycle period search and investigate its relation with rotation. The chromospheric activity index, S-index, was measured for 304 main-sequence stars from archived high-resolution HARPS spectra. Additionally, the measured and archived S-indices were converted into the chromospheric flux ratio log R'HK. The activity-cycle periods were determined using the generalised Lomb-Scargle periodogram to study the active and inactive branches on the rotation-activity-cycle period plane. The global sample shows that the bimodality of chromospheric activity, known as the Vaughan-Preston gap, is not prominent, with a significant percentage of the stars at an intermediate-activity level around log R'HK = -4.75. Independently, the cycle period search shows that stars can lie in the region intermediate between the active and inactive branch, which means that the active branch is not as clearly distinct as previously thought. The weakening of the Vaughan-Preston gap indicates that cool stars spin down from a higher activity level and settle at a lower activity level without a sudden break at intermediate activity. Some cycle periods are close to the solar value between the active and inactive branch, which suggests that the solar dynamo is most likely a common case of the stellar dynamo.
△ Less
Submitted 29 March, 2018;
originally announced March 2018.
-
Magnetic bipoles in rotating turbulence with coronal envelope
Authors:
I. R. Losada,
J. Warnecke,
A. Brandenburg,
N. Kleeorin,
I. Rogachevskii
Abstract:
The formation of sunspots and starspots is not yet fully understood and is therefore one of the major open problems in solar and stellar physics. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the f…
▽ More
The formation of sunspots and starspots is not yet fully understood and is therefore one of the major open problems in solar and stellar physics. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the flux concentrations and make them more prominent. It also allows for the formation of bipolar regions (BRs). We want to know whether the presence of an outer coronal envelope also changes the excitation conditions and the rotational dependence of NEMPI. We use direct numerical simulations and mean-field simulations. We adopt a simple two-layer model of turbulence that mimics the jump between the convective turbulent and coronal layers below and above the surface of a star, respectively. The computational domain is Cartesian and located at a certain latitude of a rotating sphere. We investigate the effects of rotation on NEMPI by changing the Coriolis number, the latitude, and the box resolution. Rotation has a strong impact on the process of BR formation. Even rather slow rotation is found to suppress their formation. However, increasing the imposed magnetic field strength also makes the structures stronger and alleviates the rotational suppression somewhat. The presence of a coronal layer itself does not significantly alleviate the effects of rotational suppression.
△ Less
Submitted 27 August, 2018; v1 submitted 12 March, 2018;
originally announced March 2018.
-
Dynamo cycles in global convection simulations of solar-like stars
Authors:
Jörn Warnecke
Abstract:
Several solar-like stars exhibit cyclic magnetic activity similar to the Sun as found in photospheric and chromospheric emission. We want to understand the physical mechanism involved in rotational dependence of these activity cycle periods. We use three-dimensional magnetohydrodynamical simulations of global convective dynamos models of solar-like stars to investigate the rotational dependency of…
▽ More
Several solar-like stars exhibit cyclic magnetic activity similar to the Sun as found in photospheric and chromospheric emission. We want to understand the physical mechanism involved in rotational dependence of these activity cycle periods. We use three-dimensional magnetohydrodynamical simulations of global convective dynamos models of solar-like stars to investigate the rotational dependency of dynamos. We further apply the test-field method to determine the $α$ effect in these simulations. We find dynamo with clear oscillating mean magnetic fields for moderately and rapidly rotating runs. For slower rotation, the field is constant or exhibit irregular cycles. In the moderately and rapidly rotating regime the cycle periods increase weakly with rotation. This behavior can be well explained with a Parker-Yoshimura dynamo wave traveling equatorward. Even though the $α$ effect becomes stronger for increasing rotation, the shear decreases steeper, causing this weak dependence on rotation. Similar as other numerical studies, we find no indication of activity branches as suggested by Brandenburg et al. (1998). However, our simulation seems to agree more with the transitional branch suggested by Distefano et al. (2017) and Olspert et al. (2017). If the Sun exhibit a dynamo wave similar as we find in our simulations, it would operate deep inside the convection zone.
△ Less
Submitted 7 May, 2018; v1 submitted 4 December, 2017;
originally announced December 2017.
-
Transition from axi- to nonaxisymmetric dynamo modes in spherical convection models of solar-like stars
Authors:
Mariangela Viviani,
Jörn Warnecke,
Maarit J. Käpylä,
Petri J. Käpylä,
Nigul Olspert,
Elizabeth M. Cole-Kodikara,
Jyri J. Lehtinen,
Axel Brandenburg
Abstract:
We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition…
▽ More
We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.
△ Less
Submitted 7 May, 2018; v1 submitted 27 October, 2017;
originally announced October 2017.
-
Long-term variations of turbulent transport coefficients in a solar-like convective dynamo simulation
Authors:
Frederick A. Gent,
Maarit J. Käpylä,
Jörn Warnecke
Abstract:
The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dip…
▽ More
The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dipolar and quadrupolar modes. We compute the turbulent transport coefficients, describing the effects of the turbulent velocity field on the mean magnetic field, using the test-field method. The test-field analysis provides a plausible explanation of the missing cycle in terms of the reduction of $α_{φφ}$ in advance of the reduced surface activity, and enhanced downward turbulent pumping during the event to confine some of the magnetic field at the bottom of the convection zone, where local maximum of magnetic energy is observed during the event. At the same time, however, a quenching of the turbulent magnetic diffusivities is observed, albeit differently distributed in depth compared to the other transport coefficients. Therefore, dedicated mean-field modelling is required for verification.
△ Less
Submitted 6 October, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.
-
A new look at sunspot formation using theory and observations
Authors:
Illa R. Losada,
Jörn Warnecke,
Kolja Glogowski,
Markus Roth,
Axel Brandenburg,
Nathan Kleeorin,
Igor Rogachevskii
Abstract:
Sunspots are of basic interest in the study of the Sun. Their relevance ranges from them being an activity indicator of magnetic fields to being the place where coronal mass ejections and flares erupt. They are therefore also an important ingredient of space weather. Their formation, however, is still an unresolved problem in solar physics. Observations utilize just 2D surface information near the…
▽ More
Sunspots are of basic interest in the study of the Sun. Their relevance ranges from them being an activity indicator of magnetic fields to being the place where coronal mass ejections and flares erupt. They are therefore also an important ingredient of space weather. Their formation, however, is still an unresolved problem in solar physics. Observations utilize just 2D surface information near the spot, but it is debatable how to infer deep structures and properties from local helioseismology. For a long time, it was believed that flux tubes rising from the bottom of the convection zone are the origin of the bipolar sunspot structure seen on the solar surface. However, this theory has been challenged, in particular recently by new surface observation, helioseismic inversions, and numerical models of convective dynamos. In this article we discuss another theoretical approach to the formation of sunspots: the negative effective magnetic pressure instability. This is a large-scale instability, in which the total (kinetic plus magnetic) turbulent pressure can be suppressed in the presence of a weak large-scale magnetic field, leading to a converging downflow, which eventually concentrates the magnetic field within it. Numerical simulations of forced stratified turbulence have been able to produce strong super-equipartition flux concentrations, similar to sunspots at the solar surface. In this framework, sunspots would only form close to the surface due to the instability constraints on stratification and rotation. Additionally, we present some ideas from local helioseismology, where we plan to use the Hankel analysis to study the pre-emergence phase of a sunspot and to constrain its deep structure and formation mechanism.
△ Less
Submitted 13 April, 2017;
originally announced April 2017.
-
Extended subadiabatic layer in simulations of overshooting convection
Authors:
Petri J. Käpylä,
Matthias Rheinhardt,
Axel Brandenburg,
Rainer Arlt,
Maarit J. Käpylä,
Andreas Lagg,
Nigul Olspert,
Jörn Warnecke
Abstract:
We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kr…
▽ More
We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.
△ Less
Submitted 23 August, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Current systems of coronal loops in 3D MHD simulations
Authors:
Jörn Warnecke,
Feng Chen,
Sven Bingert,
Hardi Peter
Abstract:
We study the magnetic field and current structure associated with a coronal loop. Through this we investigate to what extent the assumptions of a force-free magnetic field break down and where they might be justified. We analyse a 3D MHD model of the solar corona in an emerging active region with the focus on the structure of the forming coronal loops. The lower boundary of this simulation is take…
▽ More
We study the magnetic field and current structure associated with a coronal loop. Through this we investigate to what extent the assumptions of a force-free magnetic field break down and where they might be justified. We analyse a 3D MHD model of the solar corona in an emerging active region with the focus on the structure of the forming coronal loops. The lower boundary of this simulation is taken from a model of an emerging active region. As a consequence of the emerging magnetic flux and the horizontal motions at the surface a coronal loop forms self-consistently. We investigate the current density along magnetic field inside (and outside) this loop and study the magnetic and plasma properties in and around it. We find that the total current along the loop changes its sign from being antiparallel to parallel to the magnetic field. This is caused by the inclination of the loop together with the footpoint motion. Around the loop the currents form a complex non-force-free helical structure. This is directly related to a bipolar current structure at the loop footpoints at the base of the corona and a local reduction of the background magnetic field (i.e. outside the loop) caused by the plasma flow into and along the loop. The locally reduced magnetic pressure in the loop allows the loop to sustain a higher density, which is crucial for the emission in extreme UV. The acting of the flow on the magnetic field hosting the loop turns out to be also responsible for the observed squashing of the loop. The complex magnetic field and current system surrounding it can be modeled only in 3D MHD models where the magnetic field has to balance the plasma pressure. A 1D coronal loop model or a force-free extrapolation can not capture the current system and the complex interaction of the plasma and the magnetic field in the coronal loop, despite the fact that the loop is under low-$β$ conditions.
△ Less
Submitted 4 December, 2017; v1 submitted 18 November, 2016;
originally announced November 2016.
-
Convection-driven spherical shell dynamos at varying Prandtl numbers
Authors:
P. J. Käpylä,
M. J. Käpylä,
N. Olspert,
J. Warnecke,
A. Brandenburg
Abstract:
(abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at v…
▽ More
(abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at various thermal and magnetic Prandtl numbers (from 0.25 to 2 and 5, respectively). Results: We find that the rotation profiles for high thermal diffusivity show a monotonically increasing angular velocity from the bottom of the convection zone to the top and from the poles toward the equator. For sufficiently rapid rotation, a region of negative radial shear develops at mid-latitudes as the thermal diffusivity is decreased. This coincides with a change in the dynamo mode from poleward propagating activity belts to equatorward propagating ones. Furthermore, the cyclic solutions disappear at the highest magnetic Reynolds numbers. The total magnetic energy increases with the magnetic Reynolds number in the range studied here ($5-151$), but the energies of the mean magnetic fields level off at high magnetic Reynolds numbers. The differential rotation is strongly affected by the magnetic fields and almost vanishes at the highest magnetic Reynolds numbers. In some of our most turbulent cases we find that two regimes are possible where either differential rotation is strong and mean magnetic fields relatively weak or vice versa. Conclusions: Our simulations indicate a strong non-linear feedback of magnetic fields on differential rotation, leading to qualitative changes in the behaviors of large-scale dynamos at high magnetic Reynolds numbers. Furthermore, we do not find indications of the simulations approaching an asymptotic regime where the results would be independent of diffusion coefficients.
△ Less
Submitted 12 September, 2016; v1 submitted 19 May, 2016;
originally announced May 2016.
-
Turbulent transport coefficients in spherical wedge dynamo simulations of solar-like stars
Authors:
Jörn Warnecke,
Matthias Rheinhardt,
S. Tuomisto,
Petri J. Käpylä,
Maarit J. Käpylä,
Axel Brandenburg
Abstract:
We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine t…
▽ More
We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $α$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $φφ$-component of the $α$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case.
△ Less
Submitted 16 January, 2018; v1 submitted 14 January, 2016;
originally announced January 2016.
-
Limitations of force-free magnetic field extrapolations: revisiting basic assumptions
Authors:
H. Peter,
J. Warnecke,
L. P. Chitta,
R. H. Cameron
Abstract:
Force-free extrapolations are widely used to study the magnetic field in the solar corona based on surface measurements. The extrapolations assume that the ratio of internal energy of the plasma to magnetic energy, the plasma-beta is negligible. Despite the widespread use of this assumption observations, models, and theoretical considerations show that beta is of the order of a few percent to more…
▽ More
Force-free extrapolations are widely used to study the magnetic field in the solar corona based on surface measurements. The extrapolations assume that the ratio of internal energy of the plasma to magnetic energy, the plasma-beta is negligible. Despite the widespread use of this assumption observations, models, and theoretical considerations show that beta is of the order of a few percent to more than 10%, and thus not small. We investigate what consequences this has for the reliability of extrapolation results. We use basic concepts starting with the force and the energy balance to infer relations between plasma-beta and free magnetic energy, to study the direction of currents in the corona with respect to the magnetic field, and to estimate the errors in the free magnetic energy by neglecting effects of the plasma (beta<<1). A comparison with a 3D MHD model supports our basic considerations. If plasma-beta is of the order of the relative free energy (the ratio of the free magnetic energy to the total magnetic energy) then the pressure gradient can balance the Lorentz force. This is the case in the solar corona, and therefore the currents are not properly described. In particular the error in terms of magnetic energy by neglecting the plasma is of the order of the free magnetic energy, so that the latter can not be reliably determined by an extrapolation. While a force-free extrapolation might capture the magnetic structure and connectivity of the coronal magnetic field, the derived currents and free magnetic energy are not reliable. Thus quantitative results of extrapolations on the location and amount of heating in the corona (through current dissipation) and on the energy storage of the magnetic field (e.g. for eruptive events) are limited.
△ Less
Submitted 15 October, 2015;
originally announced October 2015.
-
Multiple dynamo modes as a mechanism for long-term solar activity variations
Authors:
Maarit J. Käpylä,
Petri J. Käpylä,
Nigul Olspert,
Axel Brandenburg,
Jörn Warnecke,
Bidya B. Karak,
Jaan Pelt
Abstract:
Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years,…
▽ More
Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very complex exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum. We interpret the overall irregular behavior to be due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses.
△ Less
Submitted 22 January, 2016; v1 submitted 20 July, 2015;
originally announced July 2015.
-
Influence of a coronal envelope as a free boundary to global convective dynamo simulations
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit J. Käpylä,
Axel Brandenburg
Abstract:
We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above…
▽ More
We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above. The interface represents a free surface. We compare our model to models that have no coronal envelope. The presence of a coronal envelope is found to modify the Reynolds stress and the $Λ$ effect resulting in a weaker and non-cylindrical differential rotation. This is related to the reduced latitudinal temperature variations that are caused by and dependent on the angular velocity. Some simulations develop a near-surface shear layer that we can relate to a sign change in the meridional Reynolds stress term in the thermal wind balance equation. Furthermore, the presence of a free surface changes the magnetic field evolution since the toroidal field is concentrated closer to the surface. In all simulations, however, the migration direction of the mean magnetic field can be explained by the Parker--Yoshimura rule, which is consistent with earlier findings. A realistic treatment of the upper boundary in spherical dynamo simulations is crucial for the dynamics of the flow and magnetic field evolution.
△ Less
Submitted 19 December, 2016; v1 submitted 17 March, 2015;
originally announced March 2015.
-
Bipolar region formation in stratified two-layer turbulence
Authors:
Jörn Warnecke,
Illa R. Losada,
Axel Brandenburg,
Nathan Kleeorin,
Igor Rogachevskii
Abstract:
This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We interpret the formation process in terms of negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations, we use a Cartesian domain of isothermal stratified gas that is divided into two…
▽ More
This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We interpret the formation process in terms of negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations, we use a Cartesian domain of isothermal stratified gas that is divided into two layers. In the lower layer, turbulence is forced with transverse nonhelical random waves, whereas in the upper layer no flow is induced. A weak uniform magnetic field is imposed in the entire domain at all times. In this study we vary the stratification by changing the gravitational acceleration, magnetic Reynolds number, strength of the imposed magnetic field, and size of the domain to investigate their influence on the formation process. Bipolar magnetic structure formation takes place over a large range of parameters. The magnetic structures become more intense for higher stratification until the density contrast becomes around $100$ across the turbulent layer. For the Reynolds numbers considered, magnetic flux concentrations are generated at magnetic Prandtl number between 0.1 and 1. The magnetic field in bipolar regions increases with higher imposed field strength until the field becomes comparable to the equipartition field strength of the turbulence. A larger horizontal extent enables the flux concentrations to become stronger and more coherent. The size of the bipolar structures turns out to be independent of the domain size. In the case of bipolar region formation, we find an exponential growth of the large-scale magnetic field, which is indicative of a hydromagnetic instability. Additionally, the flux concentrations are correlated with strong large-scale downward and converging flows. These findings imply that NEMPI is responsible for magnetic flux concentrations.
△ Less
Submitted 26 April, 2016; v1 submitted 12 February, 2015;
originally announced February 2015.
-
On the cause of solar-like equatorward migration in global convective dynamo simulations
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit J. Käpylä,
Axel Brandenburg
Abstract:
We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $αΩ$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due…
▽ More
We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $αΩ$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $α$ effect in the northern (southern) hemisphere and a negative radial gradient of $Ω$ outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation.
△ Less
Submitted 28 November, 2014; v1 submitted 10 September, 2014;
originally announced September 2014.
-
Coronal influence on dynamos
Authors:
Jörn Warnecke,
Axel Brandenburg
Abstract:
We report on turbulent dynamo simulations in a spherical wedge with an outer coronal layer. We apply a two-layer model where the lower layer represents the convection zone and the upper layer the solar corona. This setup is used to study the coronal influence on the dynamo action beneath the surface. Increasing the radial coronal extent gradually to three times the solar radius and changing the ma…
▽ More
We report on turbulent dynamo simulations in a spherical wedge with an outer coronal layer. We apply a two-layer model where the lower layer represents the convection zone and the upper layer the solar corona. This setup is used to study the coronal influence on the dynamo action beneath the surface. Increasing the radial coronal extent gradually to three times the solar radius and changing the magnetic Reynolds number, we find that dynamo action benefits from the additional coronal extent in terms of higher magnetic energy in the saturated stage. The flux of magnetic helicity can play an important role in this context.
△ Less
Submitted 2 October, 2013;
originally announced October 2013.
-
Bipolar magnetic structures driven by stratified turbulence with a coronal envelope
Authors:
Jörn Warnecke,
Illa R. Losada,
Axel Brandenburg,
Nathan Kleeorin,
Igor Rogachevskii
Abstract:
We report the spontaneous formation of bipolar magnetic structures in direct numerical simulations of stratified forced turbulence with an outer coronal envelope. The turbulence is forced with transverse random waves only in the lower (turbulent) part of the domain. Our initial magnetic field is either uniform in the entire domain or confined to the turbulent layer. After about 1--2 turbulent diff…
▽ More
We report the spontaneous formation of bipolar magnetic structures in direct numerical simulations of stratified forced turbulence with an outer coronal envelope. The turbulence is forced with transverse random waves only in the lower (turbulent) part of the domain. Our initial magnetic field is either uniform in the entire domain or confined to the turbulent layer. After about 1--2 turbulent diffusion times, a bipolar magnetic region of vertical field develops with two coherent circular structures that live during one turbulent diffusion time, and then decay during 0.5 turbulent diffusion times. The resulting magnetic field strengths inside the bipolar region are comparable to the equipartition value with respect to the turbulent kinetic energy. The bipolar magnetic region forms a loop-like structure in the upper coronal layer. We associate the magnetic structure formation with the negative effective magnetic pressure instability in the two-layer model.
△ Less
Submitted 29 October, 2013; v1 submitted 5 August, 2013;
originally announced August 2013.
-
Effects of enhanced stratification on equatorward dynamo wave propagation
Authors:
Petri J. Käpylä,
Maarit. J. Mantere,
Elizabeth Cole,
Jörn Warnecke,
Axel Brandenburg
Abstract:
We present results from simulations of rotating magnetized turbulent convection in spherical wedge geometry representing parts of the latitudinal and longitudinal extents of a star. Here we consider a set of runs for which the density stratification is varied, keeping the Reynolds and Coriolis numbers at similar values. In the case of weak stratification, we find quasi-steady dynamo solutions for…
▽ More
We present results from simulations of rotating magnetized turbulent convection in spherical wedge geometry representing parts of the latitudinal and longitudinal extents of a star. Here we consider a set of runs for which the density stratification is varied, keeping the Reynolds and Coriolis numbers at similar values. In the case of weak stratification, we find quasi-steady dynamo solutions for moderate rotation and oscillatory ones with poleward migration of activity belts for more rapid rotation. For stronger stratification, the growth rate tends to become smaller. Furthermore, a transition from quasi-steady to oscillatory dynamos is found as the Coriolis number is increased, but now there is an equatorward migrating branch near the equator. The breakpoint where this happens corresponds to a rotation rate that is about 3-7 times the solar value. The phase relation of the magnetic field is such that the toroidal field lags behind the radial field by about $π/2$, which can be explained by an oscillatory $α^2$ dynamo caused by the sign change of the $α$-effect about the equator. We test the domain size dependence of our results for a rapidly rotating run with equatorward migration by varying the longitudinal extent of our wedge. The energy of the axisymmetric mean magnetic field decreases as the domain size increases and we find that an $m=1$ mode is excited for a full $2π$ azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.
△ Less
Submitted 18 November, 2013; v1 submitted 11 January, 2013;
originally announced January 2013.
-
Spoke-like differential rotation in a convective dynamo with a coronal envelope
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit. J. Mantere,
Axel Brandenburg
Abstract:
We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ time…
▽ More
We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ times the value at the bottom of the spherical shell. We associate this region with the solar corona. We find solar-like differential rotation with radial contours of constant rotation rate, together with a near-surface shear layer. This non-cylindrical rotation profile is caused by a non-zero latitudinal entropy gradient that offsets the Taylor--Proudman balance through the baroclinic term. The meridional circulation is multi-cellular with a solar-like poleward flow near the surface at low latitudes. In most of the cases, the mean magnetic field is oscillatory with equatorward migration in two cases. In other cases, the equatorward migration is overlaid by stationary or even poleward migrating mean fields.
△ Less
Submitted 13 November, 2013; v1 submitted 10 January, 2013;
originally announced January 2013.
-
Solar-like differential rotation and equatorward migration in a convective dynamo with a coronal envelope
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit. J. Mantere,
Axel Brandenburg
Abstract:
We present results of convective turbulent dynamo simulations including a coronal layer in a spherical wedge. We find an equatorward migration of the radial and azimuthal fields similar to the behavior of sunspots during the solar cycle. The migration of the field coexist with a spoke-like differential rotation and anti-solar (clockwise) meridional circulation. Even though the migration extends ov…
▽ More
We present results of convective turbulent dynamo simulations including a coronal layer in a spherical wedge. We find an equatorward migration of the radial and azimuthal fields similar to the behavior of sunspots during the solar cycle. The migration of the field coexist with a spoke-like differential rotation and anti-solar (clockwise) meridional circulation. Even though the migration extends over the whole convection zone, the mechanism causing this is not yet fully understood.
△ Less
Submitted 2 November, 2012;
originally announced November 2012.
-
Magnetic twist: a source and property of space weather
Authors:
Jörn Warnecke,
Axel Brandenburg,
Dhrubaditya Mitra
Abstract:
We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or indirectly using time series of the skew-symmetric components of the magnetic correlation tens…
▽ More
We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius in the northern hemisphere.
△ Less
Submitted 18 July, 2012; v1 submitted 5 March, 2012;
originally announced March 2012.
-
Ejections of magnetic structures above a spherical wedge driven by a convective dynamo with differential rotation
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit. J. Mantere,
Axel Brandenburg
Abstract:
We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists…
▽ More
We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists of a convection zone and an extended coronal region to $\approx1.5$ times the radius of the sphere. The wedge contains a quarter of the azimuthal extent of the sphere and $150\degr$ in latitude. The magnetic field is self-consistently generated by the turbulent motions due to convection beneath the surface. Magnetic fields are found to emerge at the surface and are ejected to the coronal part of the domain. These ejections occur at irregular intervals and are weaker than in earlier work. We tentatively associate these events with coronal mass ejections on the Sun, even though our model of the solar atmosphere is rather simplistic.
△ Less
Submitted 30 September, 2012; v1 submitted 2 December, 2011;
originally announced December 2011.
-
Coronal ejections from convective spherical shell dynamos
Authors:
Jörn Warnecke,
Petri J. Käpylä,
Maarit. J. Mantere,
Axel Brandenburg
Abstract:
We present a three-dimensional model of rotating convection combined with a simplified corona in spherical coordinates. The motions in the convection zone generate a large-scale magnetic field that gets sporadically ejected into the outer layers above. Our model corona is approximately isothermal, but it includes density stratification due to gravity.
We present a three-dimensional model of rotating convection combined with a simplified corona in spherical coordinates. The motions in the convection zone generate a large-scale magnetic field that gets sporadically ejected into the outer layers above. Our model corona is approximately isothermal, but it includes density stratification due to gravity.
△ Less
Submitted 7 November, 2011;
originally announced November 2011.
-
Dynamo-driven plasmoid ejections above a spherical surface
Authors:
Jörn Warnecke,
Axel Brandenburg,
Dhrubaditya Mitra
Abstract:
We extend earlier models of turbulent dynamos with an upper, nearly force-free exterior to spherical geometry, and study how flux emerges from lower layers to the upper ones without being driven by magnetic buoyancy. We also study how this affects the possibility of plasmoid ejection. A spherical wedge is used that includes northern and southern hemispheres up to mid-latitudes and a certain range…
▽ More
We extend earlier models of turbulent dynamos with an upper, nearly force-free exterior to spherical geometry, and study how flux emerges from lower layers to the upper ones without being driven by magnetic buoyancy. We also study how this affects the possibility of plasmoid ejection. A spherical wedge is used that includes northern and southern hemispheres up to mid-latitudes and a certain range in longitude of the Sun. In radius, we cover both the region that corresponds to the convection zone in the Sun and the immediate exterior up to twice the radius of the Sun. Turbulence is driven with a helical forcing function in the interior, where the sign changes at the equator between the two hemispheres. An oscillatory large-scale dynamo with equatorward migration is found to operate in the turbulence zone. Plasmoid ejections occur in regular intervals, similar to what is seen in earlier Cartesian models. These plasmoid ejections are tentatively associated with coronal mass ejections. The magnetic helicity is found to change sign outside the turbulence zone, which is in agreement with recent findings for the solar wind.
△ Less
Submitted 17 August, 2011; v1 submitted 4 April, 2011;
originally announced April 2011.
-
Plasmoid ejections driven by dynamo action underneath a spherical surface
Authors:
Jörn Warnecke,
Axel Brandenburg,
Dhrubaditya Mitra
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 18 November, 2010;
originally announced November 2010.
-
Dynamo generated field emergence through recurrent plasmoid ejections
Authors:
Jörn Warnecke,
Axel Brandenburg
Abstract:
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 whe…
▽ More
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.
△ Less
Submitted 1 October, 2010;
originally announced October 2010.
-
Recurrent flux emergence from dynamo-generated fields
Authors:
Jörn Warnecke,
Axel Brandenburg
Abstract:
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.
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.
△ Less
Submitted 31 August, 2010;
originally announced August 2010.
-
Surface appearance of dynamo-generated large-scale fields
Authors:
J. Warnecke,
A. Brandenburg
Abstract:
Twisted magnetic fields are frequently seen to emerge above the visible surface of the Sun. This emergence is usually associated with the rise of buoyant magnetic flux structures. Here we ask how magnetic fields from a turbulent large-scale dynamo appear above the surface if there is no magnetic buoyancy. The computational domain is split into two parts. In the lower part, which we refer to as t…
▽ More
Twisted magnetic fields are frequently seen to emerge above the visible surface of the Sun. This emergence is usually associated with the rise of buoyant magnetic flux structures. Here we ask how magnetic fields from a turbulent large-scale dynamo appear above the surface if there is no magnetic buoyancy. The computational domain is split into two parts. In the lower part, which we refer to as the turbulence zone, the flow is driven by an assumed helical forcing function leading to dynamo action. Above this region, which we refer to as the exterior, a nearly force-free magnetic field is computed at each time step using the stress-and-relax method. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Strong current sheets tend to form above the neutral line, where the vertical field component vanishes. Time series of the magnetic field structure show recurrent plasmoid ejections. The degree to which the exterior field is force free is estimated as the ratio of the dot product of current density and magnetic field strength to their respective rms values. This ratio reaches values of up to 95% in the exterior. A weak outward flow is driven by the residual Lorentz force.
△ Less
Submitted 7 July, 2010; v1 submitted 19 February, 2010;
originally announced February 2010.