No Arabic abstract
We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the structure models. The eigenfrequencies obtained are compared to those of 1-D models relying on the usual phenomenologies of convection and to observations of the MDI instrument aboard SoHO. We also investigate how the magnetic activity could change the eigenfrequencies and the solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects a 1.2 kG horizontal field. All models and observed eigenfrequencies are fairly close below 3 mHz. Above 3 mHz the eigenfrequencies of the phenomenological convection models are above the observed eigenfrequencies. The frequencies of the models based on the 3-D simulations are slightly below the observed frequencies. Their maximum deviation is ~ 3 mu Hz at 3 mHz but drops below 1 mu Hz at 4 mHz. Replacing the hydrodynamic by the magnetohydrodynamic simulation increases the eigenfrequencies. The shift is negligible below 2.2 mHz and then increases linearly with frequency to reach ~ 1.7 mu Hz at 4 mHz. The impact of the simulated activity is a 14 milliarcsecond shrinking of the solar layers near the optical depth unity.
The space-borne missions have provided a wealth of highly accurate data. However, our inability to properly model the upper-most region of solar-like stars prevents us from making the best of these observations. This problem is called surface effect and a key ingredient to solve it is turbulent pressure for the computation of both the equilibrium models and the oscillations. While 3D hydrodynamic simulations help to include properly the turbulent pressure in the equilibrium models, the way this surface effect is included in the computation of stellar oscillations is still subject to uncertainties. We aim at determining how to properly include the effect of turbulent pressure and its Lagrangian perturbation in the adiabatic computation of the oscillations. We also discuss the validity of the gas-gamma model (GGM) and reduced gamma model (RGM) approximations, which have been used to compute adiabatic oscillations of equilibrium models including turbulent pressure. We use a patched model of the Sun with an inner part constructed by a 1D stellar evolution code (CESTAM) and an outer part by the 3D hydrodynamical code (CO$^5$BOLD). Then, the adiabatic oscillations are computed using the ADIPLS code for the GGM and RGM and with the MAD code imposing the adiabatic condition on an existing time-dependent convection (TDC) formalism. We show that the computation of the oscillations using the TDC formalism in the adiabatic limit improves significantly the agreement with the observed frequencies compared to the GGM and RGM. Of the components of the turbulent pressure perturbation, the perturbation of the density and advection term is found to contribute most to the frequency shift. We propose a formalism to evaluate the frequency shift due to the inclusion of the term with the turbulent pressure perturbation in the variational principle in order to extrapolate our result to other stars.
Solar activity and helioseismology show the limitation of the standard solar model and call for the inclusion of dynamical processes in both convective and radiative zones. We concentrate here on the radiative zone and first show the sensitivity of boron neutrinos to the microscopic physics included in solar models. We confront the neutrino predictions of the seismic model to all the detected neutrino fluxes. Then we compute new models of the Sun including a detailed transport of angular momentum and chemicals due to internal rotation that includes meridional circulation and shear induced turbulence. We use two stellar evolution codes: CESAM and STAREVOL to estimate the different terms. We follow three temporal evolutions of the internal rotation differing by their initial conditions: very slow, moderate and fast rotation, with magnetic braking at the arrival on the main sequence for the last two. We find that the meridional velocity in the present solar radiative zone is extremely small in comparison with those of the convective zone, smaller than 10^-6 cm/s instead of m/s. All models lead to a radial differential rotation profile but with a significantly different contrast. We compare these profiles to the presumed solar internal rotation and show that if meridional circulation and shear turbulence were the only mechanisms transporting angular momentum within the Sun, a rather slow rotation in the young Sun is favored. The transport by rotation slightly influence the sound speed profile but its potential impact on the chemicals in the transition region between radiation and convective zones. This work pushes us to pursue the inclusion of the other dynamical processes to better reproduce the present observable and to describe the young active Sun. We also need to get a better knowledge of solar gravity mode splittings to use their constraints.
Several works have reported changes of the Suns subsurface stratification inferred from f-mode or p-mode observations. Recently a non-homologous variation of the subsurface layers with depth and time has been deduced from f-modes. Progress on this important transition zone between the solar interior and the external part supposes a good understanding of the interplay between the different processes which contribute to this variation. This paper is the first of a series where we aim to study these layers from the theoretical point of view. For this first paper, we use solar models obtained with the CESAM code, in its classical form, and analyze the properties of the computed theoretical f-modes. We examine how a pure variation in the calibrated radius influences the subsurface structure and we show also the impact of an additional change of composition on the same layers. Then we use an inversion procedure to quantify the corresponding f-mode variation and their capacity to infer the radius variation. We deduce an estimate of the amplitude of the 11-year cyclic photospheric radius variation.
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.
It is widely reported that the power spectra of magnetic field and velocity fluctuations in the solar wind have power law scalings with inertial-range spectral indices of -5/3 and -3/2 respectively. Studies of solar wind turbulence have repeatedly demonstrated the impact of discontinuities and coherent structures on the measured spectral index. Whether or not such discontinuities are self-generated by the turbulence or simply observations of advected structures from the inner heliosphere has been a matter of considerable debate. This work presents a statistical study of magnetic field and velocity spectral indices over 10 years of solar-wind observations; we find that anomalously steep magnetic spectra occur in magnetically dominated intervals with negative residual energy. However, this increase in negative residual energy has no noticeable impact on the spectral index of the velocity fluctuations, suggesting that these intervals with negative residual energy correspond to intermittent magnetic structures. We show statistically that the difference between magnetic and velocity spectral indices is a monotonic function of residual energy, consistent with previous work which suggests that intermittency in fluctuations causes spectral steepening. Additionally, a statistical analysis of cross helicity demonstrates that when the turbulence is balanced (low cross-helicity), the magnetic and velocity spectral indices are not equal, which suggests that our observations of negative residual energy and intermittent structures are related to non-linear turbulent interactions rather than the presence of advected pre-existing flux-tube structures.