Source author record

A. S. Brun

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

27 published item(s)

preprint2022arXiv

MHD study of planetary magnetospheric response during extreme solar wind conditions: Earth and exoplanet magnetospheres applications

Context: The stellar wind and the interplanetary magnetic field modify the topology of planetary magnetospheres. Consequently, the hazardous effect of the direct exposition to the stellar wind, for example regarding the integrity of satellites orbiting the Earth or the habitability of exoplanets, depend upon the space weather conditions. Aims: The aim of the study is to analyze the response of an Earth-like magnetosphere for various space weather conditions and interplanetary coronal mass ejections. The magnetopause stand off distance, open-close field line boundary and plasma flows towards the planet surface are calculated. Methods: We use the MHD code PLUTO in spherical coordinates to perform a parametric study regarding the dynamic pressure and temperature of the stellar wind as well as the interplanetary magnetic field intensity and orientation. The range of the parameters analyzed extends from regular to extreme space weather conditions consistent with coronal mass ejections at the Earth orbit for the present and early periods of the Sun main sequence. In addition, implications of sub-Afvenic solar wind configurations for the Earth and exoplanet magnetospheres are analyzed. Results: The direct precipitation of the solar wind at the Earth day side in equatorial latitudes is extremely unlikely even during super coronal mass ejections. On the other hand, for early evolution phases along the Sun main sequence once the Sun rotation rate was at least $5$ times faster (< 440 Myr), the Earth surface was directly exposed to the solar wind during coronal mass ejections. Nowadays, satellites at High, Geosynchronous and Medium orbits are directly exposed to the solar wind during coronal mass ejections, because part of the orbit at the Earth day side is beyond the nose of the bow shock.

preprint2022arXiv

MOVES V. Modelling star-planet magnetic interactions of HD 189733

Magnetic interactions between stars and close-in planets may lead to a detectable signal on the stellar disk. HD 189733 is one of the key exosystems thought to harbor magnetic interactions, which may have been detected in August 2013. We present a set of twelve wind models at that period, covering the possible coronal states and coronal topologies of HD 189733 at that time. We assess the power available for the magnetic interaction and predict its temporal modulation. By comparing the predicted signal with the observed signal, we find that some models could be compatible with an interpretation based on star-planet magnetic interactions. We also find that the observed signal can be explained only with a stretch-and-break interaction mechanism, while that the Alfvén wings scenario cannot deliver enough power. We finally demonstrate that the past observational cadence of HD 189733 leads to a detection rate of only between 12 to 23%, which could explain why star-planet interactions have been hard to detect in past campaigns. We conclude that the firm confirmation of their detection will require dedicated spectroscopic observations covering densely the orbital and rotation period, combined with scarcer spectropolarimetric observations to assess the concomitant large-scale magnetic topology of the star.

preprint2022arXiv

Two-dimensional simulations of solar-like models with artificially enhanced luminosity. II. Impact on internal gravity waves

Artificially increasing the luminosity and the thermal diffusivity of a model is a common tactic adopted in hydrodynamical simulations of stellar convection. In this work, we analyse the impact of these artificial modifications on the physical properties of stellar interiors and specifically on internal gravity waves. We perform two-dimensional simulations of solar-like stars with the MUSIC code. We compare three models with different luminosity enhancement factors to a reference model. The results confirm that properties of the waves are impacted by the artificial enhancement of the luminosity and thermal diffusivity. We find that an increase in the stellar luminosity yields a decrease in the bulk convective turnover timescale and an increase in the characteristic frequency of excitation of the internal waves. We also show that a higher energy input in a model, corresponding to a larger luminosity, results in higher energy in high frequency waves. Across our tests with the luminosity and thermal diffusivity enhanced together by up to a factor of 104, our results are consistent with theoretical predictions of radiative damping. Increasing the luminosity also has an impact on the amplitude of oscillatory motions across the convective boundary. One must use caution when interpreting studies of internal gravity waves based on hydrodynamical simulations with artificially enhanced luminosity.

preprint2021arXiv

Flux ropes and dynamics of the heliospheric current sheet

Context. Solar Orbiter and PSP jointly observed the solar wind for the first time in June 2020, capturing data from very different solar wind streams, calm and Alfvénic wind as well as many dynamic structures. Aims. The aim here is to understand the origin and characteristics of the highly dynamic solar wind observed by the two probes, in particular in the vicinity of the heliospheric current sheet (HCS). Methods. We analyse the plasma data obtained by PSP and Solar Orbiter in situ during the month of June 2020. We use the Alfvén-wave turbulence MHD solar wind model WindPredict-AW, and perform two 3D simulations based on ADAPT solar magnetograms for this period. Results. We show that the dynamic regions measured by both spacecraft are pervaded with flux ropes close to the HCS. These flux ropes are also present in the simulations, forming at the tip of helmet streamers, i.e. at the base of the heliospheric current sheet. The formation mechanism involves a pressure driven instability followed by a fast tearing reconnection process, consistent with the picture of Réville et al. (2020a). We further characterize the 3D spatial structure of helmet streamer born flux ropes, which seems, in the simulations, to be related to the network of quasi-separatrices.

preprint2016arXiv

Characterizing the feedback of magnetic field on the differential rotation of solar-like stars

The aim of this article is to study how the differential rotation of solar-like stars is influenced by rotation rate and mass in presence of magnetic fields generated by a convective dynamo. We use the ASH code to model the convective dynamo of solar-like stars at various rotation rates and masses, hence different effective Rossby numbers. We obtained models with either prograde (solar-like) or retrograde (anti-solar-like) differential rotation. The trends of differential rotation versus stellar rotation rate obtained for simulations including the effect of the magnetic field are weaker compared with hydro simulations ($ΔΩ\propto (Ω/Ω_{\odot})^{0.44}$ in the MHD case and $ΔΩ\propto (Ω/Ω_{\odot})^{0.89}$ in the hydro case), hence showing a better agreement with the observations. Analysis of angular momentum transport revealed that the simulations with retrograde and prograde differential rotation have opposite distribution of the viscous, turbulent Reynolds stresses and meridional circulation contributions. The thermal wind balance is achieved in the prograde cases. However, in retrograde cases Reynolds stresses are dominant for high latitudes and near the top of the convective layer. Baroclinic effects are stronger for faster rotating models.

preprint2016arXiv

Flux-tube geometry and solar wind speed during an activity cycle

The solar wind speed at 1 AU shows variations in latitude and in time which reflect the evolution of the global background magnetic field during the activity cycle. It is commonly accepted that the terminal wind speed in a magnetic flux-tube is anti-correlated with its expansion ratio, which motivated the definition of widely-used semi-empirical scaling laws relating one to the other. In practice, such scaling laws require ad-hoc corrections. A predictive law based solely on physical principles is still missing. We test whether the flux-tube expansion is the controlling factor of the wind speed at all phases of the cycle and at all latitudes using a very large sample of wind-carrying open magnetic flux-tubes. We furthermore search for additional physical parameters based on the geometry of the coronal magnetic field which have an influence on the terminal wind flow speed. We use MHD simulations of the corona and wind coupled to a dynamo model to provide a large statistical ensemble of open flux-tubes which we analyse conjointly in order to identify relations of dependence between the wind speed and geometrical parameters of the flux-tubes which are valid globally (for all latitudes and moments of the cycle). Our study confirms that the terminal speed of the solar wind depends very strongly on the geometry of the open magnetic flux-tubes through which it flows. The total flux-tube expansion is more clearly anti-correlated with the wind speed for fast rather than for slow wind flows, and effectively controls the locations of these flows during solar minima. Overall, the actual asymptotic wind speeds attained are also strongly dependent on field-line inclination and magnetic field amplitude at the foot-points. We suggest ways of including these parameters on future predictive scaling-laws for the solar wind speed.

preprint2016arXiv

Modelling turbulent stellar convection zones: sub-grid scales effects

The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were used to obtain these results. It appears today mandatory to assess their robustness with respect to the details of the numerics, and in particular to the treatment of turbulent sub-grid scales. We report on an ongoing comparison between two global models, the ASH and EULAG codes. In EULAG the sub-grid scales are treated implicitly by the numerical scheme, while in ASH their effect is generally modelled by using enhanced dissipation coefficients. We characterize the sub-grid scales effect in a turbulent convection simulation with EULAG. We assess their effect at each resolved scale with a detailed energy budget. We derive equivalent eddy-diffusion coefficients and use the derived diffusivities in twin ASH numerical simulations. We find a good agreement between the large-scale flows developing in the two codes in the hydrodynamic regime, which encourages further investigation in the magnetohydrodynamic regime for various dynamo solutions.

preprint2016arXiv

Radio emission in Mercury magnetosphere

Context: Active stars possess magnetized wind that has a direct impact on planets that can lead to radio emission. Mercury is a good test case to study the effect of the solar wind and interplanetary magnetic field on radio emission driven in the planet magnetosphere. Such studies could be used as proxies to characterize the magnetic field topology and intensity of exoplanets. Aims: The aim of this study is to quantify the radio emission in the Hermean magnetosphere. Methods: We use the MHD code PLUTO in spherical coordinates with an axisymmetric multipolar expansion for the Hermean magnetic field, to analyze the effect of the interplanetary magnetic field (IMF) orientation and intensity, as well as the hydrodynamic parameters of the solar wind (velocity, density and temperature), on the net power dissipated on the Hermean day and night side. We apply the formalism derived by Zarka [2001, 2007] to infer the radio emission level from the net dissipated power. We perform a set of simulations with different hydrodynamic parameters of the solar wind, IMF orientations and intensities, that allow us to calculate the dissipated power distribution and infer the existence of radio emission hot spots on the planet day side, and to calculate the integrated radio emission of the Hermean magnetosphere. Results: The obtained radio emission distribution of dissipated power is determined by the IMF orientation (associated with the reconnection regions in the magnetosphere), although the radio emission strength is dependent on the IMF intensity and solar wind hydro parameters. The calculated total radio emission level is in agreement with the one estimated in [Zarka, 2001], between $5 \times 10^{5}$ and $2 \times 10^{6}$ W.

preprint2015arXiv

Characterizing the propagation of gravity waves in 3D nonlinear simulations of solar-like stars

The revolution of helio- and asteroseismology provides access to the detailed properties of stellar interiors by studying the star's oscillation modes. Among them, gravity (g) modes are formed by constructive interferences between progressive internal gravity waves (IGWs), propagating in stellar radiative zones. Our new 3D nonlinear simulations of the interior of a solar-like star allows us to study the excitation, propagation, and dissipation of these waves. The aim of this article is to clarify our understanding of the behavior of IGWs in a 3D radiative zone and to provide a clear overview of their properties. We use a method of frequency filtering that reveals the path of {individual} gravity waves of different frequencies in the radiative zone. We are able to identify the region of propagation of different waves in 2D and 3D, to compare them to the linear raytracing theory and to distinguish between propagative and standing waves (g modes). We also show that the energy carried by waves is distributed in different planes in the sphere, depending on their azimuthal wave number. We are able to isolate individual IGWs from a complex spectrum and to study their propagation in space and time. In particular, we highlight in this paper the necessity of studying the propagation of waves in 3D spherical geometry, since the distribution of their energy is not equipartitioned in the sphere.

preprint2015arXiv

Polar cap magnetic field reversals during solar grand minima: could pores play a role?

We study the magnetic flux carried by pores located outside active regions with sunspots and investigate their possible contribution to the reversal of the global magnetic field of the Sun. We find that they contain a total flux of comparable amplitude to the total magnetic flux contained in polar caps. The pores located at distances of 40--100~Mm from the closest active region have systematically the correct sign to contribute to the polar cap reversal. These pores can predominantly be found in bipolar magnetic regions. We propose that during grand minima of solar activity, such a systematic polarity trend, akin to a weak magnetic (Babcock-Leighton-like) source term could still be operating but was missed by the contemporary observers due to the limited resolving power of their telescopes.

preprint2015arXiv

The Solar-Stellar Connection

We discuss how recent advances in observations, theory and numerical simulations have allowed the stellar community to progress in its understanding of stellar convection, rotation and magnetism and to assess the degree to which the Sun and other stars share similar dynamical properties. Ensemble asteroseismology has become a reality with the advent of large time domain studies, especially from space missions. This new capability has provided improved constraints on stellar rotation and activity, over and above that obtained via traditional techniques such as spectropolarimetry or CaII H&K observations. New data and surveys covering large mass and age ranges have provided a wide parameter space to confront theories of stellar magnetism. These new empirical databases are complemented by theoretical advances and improved multi-D simulations of stellar dynamos. We trace these pathways through which a lucid and more detailed picture of magnetohydrodynamics of solar-like stars is beginning to emerge and discuss future prospects.

preprint2014arXiv

Numerical aspects of 3D stellar winds

This paper explores and compares the pitfalls of modelling the three-dimensional wind of a spherical star with a cartesian grid. Several numerical methods are compared, using either uniform and stretched grid or adaptative mesh refinement (AMR). An additional numerical complication is added, when an orbiting planet is considered. In this case a rotating frame is added to the model such that the orbiting planet is at rest in the frame of work. The three-dimensional simulations are systematically compared to an equivalent two-dimensional, axisymmetric simulation. The comparative study presented here suggests to limit the rotation rate of the rotating frame below the rotating frame of the star and provides guidelines for further three-dimensional modelling of stellar winds in the context of close-in star-planet interactions.

preprint2014arXiv

Soft X-ray emission in kink-unstable coronal loops

Solar flares are associated with intense soft X-ray emission generated by the hot flaring plasma. Kink unstable twisted flux-ropes provide a source of magnetic energy which can be released impulsively and account for the flare plasma heating. We compute the temporal evolution of the thermal X-ray emission in kink-unstable coronal loops using MHD simulations and discuss the results of with respect to solar flare observations. The model consists of a highly twisted loop embedded in a region of uniform and untwisted coronal magnetic field. We let the kink instability develop, compute the evolution of the plasma properties in the loop (density, temperature) without accounting for mass exchange with the chromosphere. We then deduce the X-ray emission properties of the plasma during the whole flaring episode. During the initial phase of the instability plasma heating is mostly adiabatic. Ohmic diffusion takes over as the instability saturates, leading to strong and impulsive heating (> 20 MK), to a quick enhancement of X-ray emission and to the hardening of the thermal X-ray spectrum. The temperature distribution of the plasma becomes broad, with the emission measure depending strongly on temperature. Significant emission measures arise for plasma at temperatures T > 9 MK. The magnetic flux-rope then relaxes progressively towards a lower energy state as it reconnects with the background flux. The loop plasma suffers smaller sporadic heating events but cools down conductively. The total thermal X-ray emission slowly fades away during this phase, and the high temperature component of emission measure distribution converges to the power-law distribution $EM\propto T^{-4.2}$. The amount of twist deduced directly from the X-ray emission patterns is considerably lower than the maximum magnetic twist in the simulated flux-ropes.

preprint2013arXiv

Flux emergence in a magnetized convection zone

We study the influence of a dynamo magnetic field on the buoyant rise and emergence of twisted magnetic flux-ropes, and their influence on the global external magnetic field. We ran 3D MHD numerical simulations using the ASH code and analysed the dynamical evolution of such buoyant flux-ropes from the bottom of the convection zone until the post-emergence phases. The global nature of this model represents very crudely and inaccurately the local dynamics of the buoyant rise, but allows to study the influence of global effects such as self-consistently generated differential rotation, meridional circulation and Coriolis forces. Although motivated by the solar context, this model cannot be thought of as a realistic model of the rise of magnetic structures and their emergence in the Sun where the local dynamics are completely different. The properties of initial phases of the buoyant rise in good agreement with previous studies. However, the effects of the interaction of the background dynamo field become increase as the flux-ropes evolve. During the buoyant rise across the CZ, the flux-rope's magnetic field strength and scales as $rho^a$, with $a\lesssim 1$. An increase of velocity, density and current precedes flux emergence at all longitudes. The geometry, latitude and relative orientation of the flux-ropes with respect to the background magnetic field influences the rise speeds, zonal flow amplitudes (which develop within the flux-ropes) and the corresponding surface signatures. This influences the morphology, duration and amplitude of the associated surface shearing and Poynting flux. The emerged flux influences the system's global polarity, leading in some cases to a polarity reversal while inhibiting background dynamo from doing so in some others. The emerged magnetic flux is slowly advected poleward, while being diffused and assimilated by the background dynamo field.

preprint2013arXiv

Magnetic energy cascade in spherical geometry: I. The stellar convective dynamo case

We present a method to characterize the spectral transfers of magnetic energy between scales in simulations of stellar convective dynamos. The full triadic transfer functions are computed thanks to analytical coupling relations of spherical harmonics based on the Clebsch-Gordan coefficients. The method is applied to mean field $αΩ$ dynamo models as benchmark tests. From the physical standpoint, the decomposition of the dynamo field into primary and secondary dynamo families proves very instructive in the $αΩ$ case. The same method is then applied to a fully turbulent dynamo in a solar convection zone, modeled with the 3D MHD ASH code. The initial growth of the magnetic energy spectrum is shown to be non-local. It mainly reproduces the kinetic energy spectrum of convection at intermediate scales. During the saturation phase, two kinds of direct magnetic energy cascades are observed in regions encompassing the smallest scales involved in the simulation. The first cascade is obtained through the shearing of magnetic field by the large scale differential rotation that effectively cascades magnetic energy. The second is a generalized cascade that involves a range of local magnetic and velocity scales. Non-local transfers appear to be significant, such that the net transfers cannot be reduced to the dynamics of a small set of modes. The saturation of the large scale axisymmetric dipole and quadrupole are detailed. In particular, the dipole is saturated by a non-local interaction involving the most energetic scale of the magnetic energy spectrum, which points out the importance of the magnetic Prandtl number for large-scale dynamos.

preprint2013arXiv

On close-in magnetized star-planet interactions

We present 2D magnetohydrodynamic simulations performed with the PLUTO code to model magnetized star-planet interactions. We study two simple scenarios of magnetized star-planet interactions: the unipolar and dipolar interactions. Despite the simplified hypotheses we consider in the model, the qualitative behavior of the interactions is very well recovered. These encouraging results promote further developments of the model to obtain predictions on the effect and the physical manifestation of magnetized star--close-in planet interactions.

preprint2013arXiv

On dynamo action in the giant star Pollux : first results

We present preliminary results of a 3D MHD simulation of the convective envelope of the giant star Pollux for which the rotation period and the magnetic field intensity have been measured from spectroscopic and spectropolarimetric observations. This giant is one of the first single giants with a detected magnetic field, and the one with the weakest field so far. Our aim is to understand the development and the action of the dynamo in its extended convective envelope.

preprint2013arXiv

World-leading science with SPIRou - the nIR spectropolarimeter / high-precision velocimeter for CFHT

SPIRou is a near-infrared (nIR) spectropolarimeter / velocimeter proposed as a new-generation instrument for CFHT. SPIRou aims in particular at becoming world-leader on two forefront science topics, (i) the quest for habitable Earth-like planets around very- low-mass stars, and (ii) the study of low-mass star and planet formation in the presence of magnetic fields. In addition to these two main goals, SPIRou will be able to tackle many key programs, from weather patterns on brown dwarf to solar-system planet atmospheres, to dynamo processes in fully-convective bodies and planet habitability. The science programs that SPIRou proposes to tackle are forefront (identified as first priorities by most research agencies worldwide), ambitious (competitive and complementary with science programs carried out on much larger facilities, such as ALMA and JWST) and timely (ideally phased with complementary space missions like TESS and CHEOPS). SPIRou is designed to carry out its science mission with maximum efficiency and optimum precision. More specifically, SPIRou will be able to cover a very wide single-shot nIR spectral domain (0.98-2.35 μm) at a resolving power of 73.5K, providing unpolarized and polarized spectra of low-mass stars with a ~15% average throughput and a radial velocity (RV) precision of 1 m/s.

preprint2012arXiv

Solar Magnetic Field Reversals and the Role of Dynamo Families

The variable magnetic field of the solar photosphere exhibits periodic reversals as a result of dynamo activity occurring within the solar interior. We decompose the surface field as observed by both the Wilcox Solar Observatory and the Michelson Doppler Imager into its harmonic constituents, and present the time evolution of the mode coefficients for the past three sunspot cycles. The interplay between the various modes is then interpreted from the perspective of general dynamo theory, where the coupling between the primary and secondary families of modes is found to correlate with large-scale polarity reversals for many examples of cyclic dynamos. Mean-field dynamos based on the solar parameter regime are then used to explore how such couplings may result in the various long-term trends in the surface magnetic field observed to occur in the solar case.

preprint2011arXiv

Convection and Differential Rotation Properties of G & K Stars Computed with the ASH Code

The stellar luminosity and depth of the convective envelope vary rapidly with mass for G- and K-type main sequence stars. In order to understand how these properties influence the convective turbulence, differential rotation, and meridional circulation, we have carried out 3D dynamical simulations of the interiors of rotating main sequence stars, using the anelastic spherical harmonic (ASH) code. The stars in our simulations have masses of 0.5, 0.7, 0.9, and 1.1 M_sun, corresponding to spectral types K7 through G0, and rotate at the same angular speed as the sun. We identify several trends of convection zone properties with stellar mass, exhibited by the simulations. The convective velocities, temperature contrast between up- and down-flows, and meridional circulation velocities all increase with stellar luminosity. As a consequence of the trend in convective velocity, the Rossby number (at a fixed rotation rate) increases and the convective turnover timescales decrease significantly with increasing stellar mass. The 3 lowest mass cases exhibit solar-like differential rotation, in a sense that they show a maximum rotation at the equator and minimum at higher latitudes, but the 1.1 M_sun case exhibits anti-solar rotation. At low mass, the meridional circulation is multi-cellular and aligned with the rotation axis; as the mass increases, the circulation pattern tends toward a unicellular structure covering each hemisphere in the convection zone.

preprint2011arXiv

Fast Rotating solar-like stars using asteroseismic datasets

The NASA Kepler mission is providing an unprecedented set of asteroseismic data. In particular, short-cadence lightcurves (~60s samplings), allow us to study solar-like stars covering a wide range of masses, spectral types and evolutionary stages. Oscillations have been observed in around 600 out of 2000 stars observed for one month during the survey phase of the Kepler mission. The measured light curves can present features related to the surface magnetic activity (starspots) and, thus we are able to obtain a good estimation of the surface (differential) rotation. In this work we establish the basis of such research and we show a potential method to find stars with fast surface rotations.

preprint2011arXiv

Magnetic confinement of the solar tachocline: II. Coupling to a convection zone

The reason for the observed thinness of the solar tachocline is still not well understood. One of the explanations that have been proposed is that a primordial magnetic field renders the rotation uniform in the radiation zone. We test here the validity of this magnetic scenario through 3D numerical MHD simulations that encompass both the radiation zone and the convection zone. The numerical simulations are performed with the anelastic spherical harmonics (ASH) code. The computational domain extends from $0.07\;R_\odot$ to $0.97\;R_\odot$. In the parameter regime we explored, a dipolar fossil field aligned with the rotation axis can not remain confined in the radiation zone. When the field lines are allowed to interact with turbulent unstationary convective motions at the base of the convection zone, 3D effects prevent the field confinement. In agreement with previous work, we find that a dipolar fossil field, even when it is initially buried deep inside the radiation zone, will spread into the convective zone. According to Ferraro's law of iso-rotation, it then imprints on the radiation zone the latitudinal differential rotation of the convection zone, which is not observed.

preprint2009arXiv

Exploring the $P_{cyc}$ vs $P_{rot}$ relation with flux transport dynamo models of solar-like stars

Aims: To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field dynamo models with stellar activity observations Methods: 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational data. The novelty is that these 2-D mean field dynamo models incorporate scaling laws deduced from 3-D hydrodynamical simulations for the influence of rotation rate on the amplitude and profile of the meridional circulation. These models make also use of observational scaling laws for the variation of differential rotation with rotation rate. Results: We find that Babcock-Leighton flux transport dynamo models are able to reproduce the change in topology of the magnetic field (i.e. toward being more toroidal with increasing rotation rate) but seem to have difficulty reproducing the cycle period vs activity period correlation observed in solar-like stars if a monolithic single cell meridional flow is assumed. It may however be possible to recover the $P_{cyc}$ vs $P_{rot}$ relation with more complex meridional flows, if the profile changes in a particular assumed manner with rotation rate. Conclusions: The Babcock-Leighton flux transport dynamo model based on single cell meridional circulation does not reproduce the $P_{cyc}$ vs $P_{rot}$ relation unless the amplitude of the meridional circulation is assumed to increase with rotation rate which seems to be in contradiction with recent results obtained with 3-D global simulations.

preprint2009arXiv

Is the solar convection zone in strict thermal wind balance?

Context: The solar rotation profile is conical rather than cylindrical as one could expect from classical rotating fluid dynamics (e.g. Taylor-Proudman theorem). Thermal coupling to the tachocline, baroclinic effects and latitudinal transport of heat have been advocated to explain this peculiar state of rotation. Aims: To test the validity of thermal wind balance in the solar convection zone using helioseismic inversions for both the angular velocity and fluctuations in entropy and temperature. Methods: Entropy and temperature fluctuations obtained from 3-D hydrodynamical numerical simulations of the solar convection zone are compared with solar profiles obtained from helioseismic inversions. Results: The temperature and entropy fluctuations in 3-D numerical simulations have smaller amplitude in the bulk of the solar convection zone than those found from seismic inversions. Seismic inversion find variations of temperature from about 1 K at the surface up to 100 K at the base of the convection zone while in 3-D simulations they are of order 10 K throughout the convection zone up to 0.96 $R_{\odot}$. In 3-D simulations, baroclinic effects are found to be important to tilt the isocontours of $Ω$ away from a cylindrical profile in most of the convection zone helped by Reynolds and viscous stresses at some locations. By contrast the baroclinic effect inverted by helioseismology are much larger than what is required to yield the observed angular velocity profile. Conclusion: The solar convection does not appear to be in strict thermal wind balance, Reynolds stresses must play a dominant role in setting not only the equatorial acceleration but also the observed conical angular velocity profile.

preprint2006arXiv

On the interactions of turbulent convection and rotation in RGB stars

We have performed the first three-dimensional non-linear simulation of the turbulent convective envelope of a rotating 0.8 Msun RGB star using the ASH code. Adopting a global typical rotation rate of a tenth of the solar rate, we have analyzed the dynamical properties of the convection and the transport of angular momentum within the inner 50% in radius of the convective envelope. The convective patterns consist of a small number of large cell, associated with fast flows (about 3000 m/s) and large temperature fluctuations (about 300 K) in order to carry outward the large luminosity (L* = 400 Lsun) of the star. The interactions between convection and rotation give rise to a large radial differential rotation and a meridional circulation possessing one cell per hemisphere, the flow being poleward in both hemisphere. By analysing the redistribution of angular momentum, we find that the meridional circulation transports the angular momentum outward in the radial direction, and poleward in the latitudinal direction, and that the transport by Reynolds stresses acts in the opposite direction. From this 3-D simulation, we have derived an average radial rotation profile, that we will ultimately introduce back into 1-D stellar evolution code.

preprint2000arXiv

Influence of the Tachocline on Solar Evolution

Recently helioseismic observations have revealed the presence of a shear layer at the base of the convective zone related to the transition from differential rotation in the convection zone to almost uniform rotation in the radiative interior, the tachocline. At present, this layer extends only over a few percent of the solar radius and no definitive explanations have been given for this thiness. Following Spiegel and Zahn (1992, Astron. Astrophys.), who invoke anisotropic turbulence to stop the spread of the tachocline deeper in the radiative zone as the Sun evolves, we give some justifications for their hypothesis by taking into account recent results on rotating shear instability (Richard and Zahn 1999, Astron. Astrophys.). We study the impact of the macroscopic motions present in this layer on the Sun's structure and evolution by introducing a macroscopic diffusivity $D_T$ in updated solar models. We find that a time dependent treatment of the tachocline significantly improves the agreement between computed and observed surface chemical species, such as the $^7$Li and modify the internal structure of the Sun (Brun, Turck-Chièze and Zahn, 1999, in Astrophys. J.).

preprint1999arXiv

Standard Solar models in the Light of New Helioseismic Constraints II. Mixing Below the Convective Zone

In previous work, we have shown that recent updated standard solar models cannot reproduce the radial profile of the sound speed at the base of the convective zone (CZ) and fail to predict the Li7 depletion. In parallel, helioseismology has shown that the transition from differential rotation in the CZ to almost uniform rotation in the radiative solar interior occurs in a shallow layer called the tachocline. This layer is presumably the seat of large scale circulation and of turbulent motions. Here, we introduce a macroscopic transport term in the structure equations, which is based on a hydrodynamical description of the tachocline proposed by Spiegel and Zahn, and we calculate the mixing induced within this layer. We discuss the influence of different parameters that represent the tachocline thickness, the Brunt-Vaissala frequency at the base of the CZ, and the time dependence of this mixing process along the Sun's evolution. We show that the introduction of such a process inhibits the microscopic diffusion by about 25%. Starting from models including a pre-main sequence evolution, we obtain: a) a good agreement with the observed photospheric chemical abundance of light elements such as He3, He4, Li7 and Be9, b) a smooth composition gradient at the base of the CZ, and c) a significant improvement of the sound speed square difference between the seismic sun and the models in this transition region, when we allow the phostospheric heavy element abundance to adjust, within the observational incertitude, due to the action of this mixing process. The impact on neutrino predictions is also discussed.