No Arabic abstract
We perform 2D, fully compressible, time-implicit simulations of convection in a solar-like model with the MUSIC code. Our main motivation is to explore the impact of a common tactic adopted in numerical simulations of convection that use realistic stellar conditions. This tactic is to artificially increase the luminosity and to modify the thermal diffusivity of the reference stellar model. This work focuses on the impact of these modifications on convective penetration (or overshooting) at the base of the convective envelope of a solar-like model. We explore a range of enhancement factors for the energy input and confirm the increase in the characteristic overshooting depth with the increase in the energy input. Our results highlight the importance of the impact of penetrative downflows on the thermal background below the convective boundary. This is a result of compression and shear which induce local heating and thermal mixing. The artificial increase in the energy flux intensifies the heating process by increasing the velocities in the convective zone and at the convective boundary, revealing a subtle connection between the local heating of the thermal background and the plume dynamics. This heating also increases the efficiency of heat transport by radiation which may counterbalance further heating and helps to establish a steady state. The modification of the thermal background by penetrative plumes impacts the width of the overshooting layer. Our results suggest that an artificial modification of the radiative diffusivity in the overshooting layer, rather than only accelerating the thermal relaxation, could also alter the dynamics of the penetrating plumes and thus the width of the overshooting layer. Results from simulations with an artificial modification of the energy flux and of the thermal diffusivity should be regarded with caution if used to determine an overshooting width.
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 $alphaOmega$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $alpha$ effect in the northern (southern) hemisphere and a negative radial gradient of $Omega$ 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.
We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large Peclet number regime. These results should contribute to the 321D link.
The last decade has seen a rapid development in asteroseismology thanks to the CoRoT and Kepler missions. With more detailed asteroseismic observations available, it is becoming possible to infer exactly how oscillations are driven and dissipated in solar-type stars. We have carried out three-dimensional (3D) stellar atmosphere simulations together with one-dimensional (1D) stellar structural models of key benchmark turn-off and subgiant stars to study this problem from a theoretical perspective. Mode excitation and damping rates are extracted from 3D and 1D stellar models based on analytical expressions. Mode velocity amplitudes are determined by the balance between stochastic excitation and linear damping, which then allows the estimation of the frequency of maximum oscillation power, $ u_{max}$, for the first time based on ab initio and parameter-free modelling. We have made detailed comparisons between our numerical results and observational data and achieved very encouraging agreement for all of our target stars. This opens the exciting prospect of using such realistic 3D hydrodynamical stellar models to predict solar-like oscillations across the HR-diagram, thereby enabling accurate estimates of stellar properties such as mass, radius and age.
Using a thin flux tube model in a rotating spherical shell of turbulent, solar-like convective flows, we find that the distribution of emerging flux tubes in our simulation is inhomogeneous in longitude, with properties similar to those of active longitudes on the Sun and other solar-like stars. The large-scale pattern of flux emergence our simulations produce exhibits preferred longitudinal modes of low order, drift with respect to a fixed reference system, and alignment across the Equator at low latitudes between 15 degrees. We suggest that these active-longitude-like emergence patterns are the result of columnar, rotationally aligned giant cells present in our convection simulation at low latitudes. If giant convecting cells exist in the bulk of the solar convection zone, this phenomenon, along with differential rotation, could in part provide an explanation for the behavior of active longitudes.
We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from $sim$ 50 Myr to $sim$ 4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long standing problem of lithium depletion in pre-main sequence and main sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.