Source author record

M. Rieutord

M. Rieutord 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

23works
8topics
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

23 published item(s)

preprint2022arXiv

Spectroscopic detection of Altair's non-radial pulsations

Context: Rapid rotation is a common feature of early-type stars but which remains a challenge for the models. The understanding of its effect on stellar evolution is however imperative to interpret the observed properties of numerous stars. Aims: We wish to bring more observational constraints on the properties of fast rotating stars, especially on their oscillation modes. Methods: We focus on the nearby star Altair which is known as a very rapidly rotating star with an equatorial velocity estimated recently at 313 km/s. We observed this star with the high-resolution spectropolarimeter Neo-Narval during six nights, with one night of interruption, in September 2020. Results: We detect significant line profile variations on the mean line profile of the spectra. Their time-frequency analysis shows that these variations are induced by gravito-inertial waves propagating at Altair's surface with azimuthal wavenumbers of order $m=10-15$. With a preliminary computation of the eigenspectrum using the most recent concordance model of Altair we can give a first modelling of the observed waves. Conclusions: Altair was known as the brightest $δ$ Scuti star. We now see that it is the brightest hybrid oscillating star with excited gravito-inertial waves and acoustic waves. Clearly, more observations and more advanced models are needed to explain the observations in greater details

preprint2020arXiv

Oscillations of 2D ESTER models. I. The adiabatic case

Recent numerical and theoretical considerations have shown that low-degree acoustic modes in rapidly rotating stars follow an asymptotic formula and recent observations of pulsations in rapidly rotating delta Scuti stars seem to match these expectations. However, a key question is whether strong gradients or discontinuities can adversely affect this pattern to the point of hindering its identification. Other important questions are how rotational splittings are affected by the 2D rotation profiles expected from baroclinic effects and whether it is possible to probe the rotation profile using these splittings. Accordingly, we numerically calculate pulsation modes in continuous and discontinuous rapidly rotating models produced by the 2D ESTER (Evolution STEllaire en Rotation) code. This spectral multi-domain code self-consistently calculates the rotation profile based on baroclinic effects and allows us to introduce discontinuities without loss of numerical accuracy. Pulsations are calculated using an adiabatic version of the Two-dimensional Oscillation Program (TOP) code. The variational principle is used to confirm the high accuracy of the pulsation frequencies and to derive an integral formula that closely matches the generalised rotational splittings, except when modes are involved in avoided crossings. This potentially allows us to probe the the rotation profile using inverse theory. Acoustic glitch theory, applied along the island mode orbit deduced from ray dynamics, can correctly predict the periodicity of the glitch frequency pattern produced by a discontinuity or the Gamma1 dip related to the He II ionisation zone in some of the models. The asymptotic frequency pattern remains sufficiently well preserved to potentially allow its detection in observed stars.

preprint2019arXiv

A realistic two-dimensional model of Altair

Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-the-art stellar interior and atmosphere models together with interferometric, spectroscopic, and asteroseismic observations. We use ESTER 2D stellar models to produce the relevant surface parameters needed to create intensity maps from atmosphere models. Interferometric and spectroscopic observables are computed from these intensity maps and several stellar parameters are then adjusted using the MCMC algorithm Emcee. We determined Altair's equatorial radius to be 2.008 +/- 0.006 Rsun, the position angle 301.1 +/- 0.3 degrees, the inclination 50.7 +/- 1.2 degrees, and the equatorial angular velocity 0.74 +/- 0.01 times the Keplerian angular velocity. This angular velocity leads to a flattening of 0.220 +/- 0.003. We also deduce from the spectroscopically derived vsini ~ 243 km/s, a true equatorial velocity of ~314 km/s corresponding to a rotation period of 7h46m (~3 c/d). The data also impose a strong correlation between mass, metallicity, hydrogen abundance, and core evolution. Thanks to asteroseismic data, we constrain the mass of Altair to 1.86 +/- 0.03 Msun and further deduce its metallicity Z = 0.019 and its core hydrogen mass fraction Xc = 0.71, assuming an initial solar hydrogen mass fraction X = 0.739. These values suggest that Altair is ~100 Myrs old. Finally, the 2D ESTER model also gives the internal differential rotation of Altair, showing that its core rotates approximately 50% faster than the envelope, while the surface differential rotation does not exceed 6%.

preprint2016arXiv

An algorithm for computing the 2D structure of fast rotating stars

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.

preprint2016arXiv

Relation between trees of fragmenting granules and supergranulation evolution

Context: The determination of the underlying mechanisms of the magnetic elements diffusion over the solar surface is still a challenge. Understanding the formation and evolution of the solar network (NE) is a challenge, because it provides a magnetic flux over the solar surface comparable to the flux of active regions at solar maximum. Aims: We investigate the structure and evolution of interior cells of solar supergranulation. From Hinode observations, we explore the motions on solar surface at high spatial and temporal resolution. We derive the main organization of the flows inside supergranules and their effect on the magnetic elements. Method: To probe the superganule interior cell, we used the Trees of Fragmenting Granules (TFG) evolution and their relations to horizontal Results: Evolution of TFG and their mutual interactions result in cumulative effects able to build horizontal coherent flows with longer lifetime than granulation (1 to 2 hours) over a scale up to 12\arcsec. These flows clearly act on the diffusion of the intranetwork (IN) magnetic elements and also on the location and shape of the network. Conclusions: From our analysis during 24 hours, TFG appear as one of the major elements of the supergranules which diffuse and advect the magnetic field on the Sun's surface. The strongest supergranules contribute the most to magnetic flux diffusion in the solar photosphere.

preprint2016arXiv

Tidal inertial waves in the differentially rotating convective envelopes of low-mass stars - I. Free oscillation modes

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. In low-mass stars, their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. Turbulent convection can sustain differential rotation in their envelope, with an equatorial acceleration (as in the Sun) or deceleration, which can modify the waves' propagation properties. We explore in this first paper the general propagation properties of free linear inertial waves in a differentially rotating homogeneous fluid inside a spherical shell. We assume that the angular velocity background flow depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We use i) an analytical approach in the inviscid case to get the dispersion relation, from which we compute the characteristic trajectories along which energy propagates. This allows us to study the existence of attractor cycles and infer the different families of inertial modes; ii) high-resolution numerical calculations based on a spectral method for the viscous problem. We find that modes that propagate in the whole shell (D modes) behave the same way as with solid-body rotation. However, another family of inertial modes exists (DT modes), which can propagate only in a restricted part of the convective zone. Our study shows that they are less common than D modes and that the characteristic rays and shear layers often focus towards a wedge - or point-like attractor. More importantly, we find that for non-axisymmetric oscillation modes, shear layers may cross a corotation resonance with a local accumulation of kinetic energy. Their damping rate scales very differently from what we obtain for standard D modes and we show an example where it is independent of viscosity (Ekman number) in the astrophysical regime in which it is small.

preprint2015arXiv

Free inertial modes in differentially rotating convective envelopes of low-mass stars : numerical exploration

Tidally-excited inertial waves in stellar convective regions are a key mechanism for tidal dissipation in stars and therefore the evolution of close-in binary or planetary systems. As a first step, we explore here the impact of latitudinal differential rotation on the properties of free inertial modes and identify the different families of modes. We show that they differ from the case of solid-body rotation. Using an analytical approach as well as numerical calculations, we conclude that critical layers (where the Doppler-shifted frequency vanishes) could play a very important role for tidal dissipation.

preprint2014arXiv

Inertial waves in differentially rotating low-mass stars and tides

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. Their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. If the star is assumed to be rotating as a solid-body, the waves' Doppler-shifted frequency is restricted to $[-2 Ω, 2 Ω]$ ($Ω$ being the angular velocity of the star) and they can propagate in the entire convective region. However, turbulent convection can sustain differential rotation with an equatorial acceleration (as in the Sun) or deceleration that may modify waves propagation. We thus explore the properties of inertial modes of oscillation in a conically differentially rotating background flow whose angular velocity depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We find that their frequency range is broadened by differential rotation, and that they may propagate only in a restricted part of the envelope. In some cases, inertial waves form shear layers around short-period attractor cycles. In others, they exhibit a remarkable behavior when a turning surface or a corotation layer exists in the star. We discuss how all these cases can impact tidal dissipation in stars.

preprint2014arXiv

Structure and evolution of solar supergranulation using SDO/HMI data

Context: Studying the motions on the solar surface is fundamental for understanding how turbulent convection transports energy and how magnetic fields are distributed across the solar surface. Aims: From horizontal velocity measurements all over the visible disc of the Sun and using data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), we investigate the structure and evolution of solar supergranulation. Methods: Horizontal velocity fields were measured by following the proper motions of solar granules using a newly developed version of the coherent structure tracking (CST) code. With this tool, maps of horizontal divergence were computed. We then segmented and identified supergranular cells and followed their histories by using spatio-temporal labelling. With this dataset we derived the fundamental properties of supergranulation, including their motion. Results: We find values of the fundamental parameters of supergranulation similar to previous studies: a mean lifetime of 1.5 days and a mean diameter of 25~Mm. The tracking of individual supergranular cells reveals the solar differential rotation and a poleward circulation trend of the meridional flow. The shape of the derived differential rotation and meridional flow does not depend on the cell size. If there is a background magnetic field, the diverging flows in supergranules are weaker. Conclusions: This study confirms that supergranules are suitable tracers that may be used to investigate the large-scale flows of the solar convection as long as they are detectable enough on the surface.

preprint2013arXiv

Comparison of solar horizontal velocity fields from SDO/HMI and Hinode data

The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ Helioseismic and Magnetic Imager (SDO/HMI) data. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured from HMI white light data. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the two components of the horizontal velocity Vx and Vy measured in the sky plane and the simultaneous line of sight component from SDO/HMI dopplergrams v_D, we derive the spherical velocity components (Vr, Vtheta, Vphi). The azimuthal component Vphi gives the solar differential rotation with a high precision (+-0.037km/s) from a temporal sequence of only three hours. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial and temporal resolutions from hours to months and years with the SDO data.

preprint2013arXiv

Comparison of solar surface flows inferred from time--distance helioseismology and coherent structure tracking using HMI/SDO observations

We compare measurements of horizontal flows on the surface of the Sun using helioseismic time--distance inversions and coherent structure tracking of solar granules. Tracking provides 2D horizontal flows on the solar surface, whereas the time--distance inversions estimate the full 3-D velocity flows in the shallow near-surface layers. Both techniques use HMI observations as an input. We find good correlations between the various measurements resulting from the two techniques. Further, we find a good agreement between these measurements and the time-averaged Doppler line-of-sight velocity, and also perform sanity checks on the vertical flow that resulted from the 3-D time--distance inversion.

preprint2013arXiv

Present status of two-dimensional ESTER models: Application to Be stars

ESTER two-dimensional models solve the steady state structure of fast rotating early-type stars including the large scale flows associated with the baroclinicity of the radiative zones. Models are compared successfully to the fundamental parameters of the two main components of the triple system $δ$ Velorum that have been derived from interferometric and orbit measurements. Testing the models on the Be star Achernar ($α$ Eri), we cannot reproduce the data and conclude that this star has left the main sequence and is likely crossing the Herzsprung gap. Computing main sequence evolution of fast rotating stars at constant angular momentum shows that their criticality increases with time suggesting that the Be phenomenon and the ensuing mass ejections is the result of evolution.

preprint2012arXiv

Acoustic Events in the Solar Atmosphere from Hinode/SOT NFI observations

We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet sun using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard \textit{Hinode}. Line profiles of Fe \textsc{i} 557.6 nm were recorded by the Narrow band Filter Imager (NFI) on a $82" \times 82"$ FOV during 75 min with a time step of 28.75 s and 0.08$"$ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130 and 180 km) using the bisector technique allowing the determination of energy flux in the range 3-10 mHz using two complementary methods (Hilbert transform and Fourier power spectra). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. The net energy flux is upward. In the range 3-10 mHz, a full FOV space and time averaged flux of 2700 W m$^{-2}$ (lower layer 80-130 km) and 2000 W m$^{-2}$ (upper layer 130-180 km) is concentrated in less than 1% of the solar surface in the form of narrow (0.3$"$) AE. Their total duration (including rise and decay) is of the order of $10^{3}$ s. Inside each AE, the mean flux is $1.6 10^{5}$ W m$^{-2}$ (lower layer) and $1.2 10^{5}$ W m$^{-2}$ (upper). Each event carries an average energy (flux integrated over space and time) of $2.5 10^{19}$ J (lower layer) to $1.9 10^{19}$ J (upper). More than $10^{6}$ events could exist permanently on the Sun, with a birth and decay rate of 3500 s$^{-1}$. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.

preprint2012arXiv

Inertial waves in a differentially rotating spherical shell

We investigate the properties of small-amplitude inertial waves propagating in a differentially rotating incompressible fluid contained in a spherical shell. For cylindrical and shellular rotation profiles and in the inviscid limit, inertial waves obey a second-order partial differential equation of mixed type. Two kinds of inertial modes therefore exist, depending on whether the hyperbolic domain where characteristics propagate covers the whole shell or not. The occurrence of these two kinds of inertial modes is examined, and we show that the range of frequencies at which inertial waves may propagate is broader than with solid-body rotation. Using high-resolution calculations based on a spectral method, we show that, as with solid-body rotation, singular modes with thin shear layers following short-period attractors still exist with differential rotation. They exist even in the case of a full sphere. In the limit of vanishing viscosities, the width of the shear layers seems to weakly depend on the global background shear, showing a scaling in E^{1/3} with the Ekman number E, as in the solid-body rotation case. There also exist modes with thin detached layers of width scaling with E^{1/2} as Ekman boundary layers. The behavior of inertial waves with a corotation resonance within the shell is also considered. For cylindrical rotation, waves get dramatically absorbed at corotation. In contrast, for shellular rotation, waves may cross a critical layer without visible absorption, and such modes can be unstable for small enough Ekman numbers.

preprint2012arXiv

Quasi full-disk maps of solar horizontal velocities using SDO/HMI data

For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the entire visible surface of the Sun, using SDO/HMI white-light data. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at different scales with its derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~0.25 km/s for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which range between 0.3 km/s and 1.8 km/s. A high correlation between velocities computed from Hinode and SDO/HMI has been found (85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar differential rotation, and the meridional velocity.

preprint2011arXiv

2D computations of g modes

We present complete 2D computations of g modes in distorted polytropic models of stars performed with the Two-dimensional Oscillation Program (TOP). We computed low-degree modes (l=1 modes with radial order n=-1...-14, and l=2,3 modes with n=-1...-5 and -16...-20) of a nonrotating model and followed them by slowly increasing the rotation rate up to 70 % of the Keplerian break-up velocity. We use these computations to determine the domain of validity of perturbative methods up to the 3rd order. We study the evolution of the regularities of the spectrum and show quantitative agreement with the traditional approximation for not too large values of the ratio of the rotation rate to the pulsation frequency. We also show the appearance of new types of modes, called "rosette" modes due to their spatial structure. Thanks to the ray theory for gravito-inertial waves that we developed, we can associate these modes with stable periodic rays.

preprint2011arXiv

A $r$-mode in a magnetic rotating spherical layer: application to neutron stars

The impact of the combination of rotation and magnetic fields on oscillations of stellar fluids is still not well known theoretically. It mixes Alfvén and inertial waves. Neutron stars are a place where both effects may be at work. We wish to decipher the solution of this problem in the context of $r$-modes instability in neutron stars, as it appears when these modes are coupled to gravitational radiation. We consider a rotating spherical shell filled with a viscous fluid but of infinite electrical conductivity and analyze propagation of modal perturbations when a dipolar magnetic field is bathing the fluid layer. We perform an extensive numerical analysis and find that the $m=2$ $r$-mode oscillation is influenced by the magnetic field when the Lehnert number (ratio of Alfvén speed to rotation speed) exceeds a value proportional to the one-fourth power of the Ekman number (non-dimensional measure of viscosity). This scaling is interpreted as the coincidence of the width of internal shear layers of inertial modes and the wavelength of the Alfvén waves. Applied to the case of rotating magnetic neutron stars, we find that dipolar magnetic fields above $10^{14}$ G are necessary to perturb the $r$-modes instability.

preprint2010arXiv

Gravity modes in rapidly rotating stars. Limits of perturbative methods

CoRoT and Kepler missions are now providing high-quality asteroseismic data for a large number of stars. Among intermediate-mass and massive stars, fast rotators are common objects. Taking the rotation effects into account is needed to correctly understand, identify, and interpret the observed oscillation frequencies of these stars. A classical approach is to consider the rotation as a perturbation. In this paper, we focus on gravity modes, such as those occurring in gamma Doradus, slowly pulsating B (SPB), or Be stars. We aim to define the suitability of perturbative methods. With the two-dimensional oscillation program (TOP), we performed complete computations of gravity modes -including the Coriolis force, the centrifugal distortion, and compressible effects- in 2-D distorted polytropic models of stars. We started with the modes l=1, n=1-14, and l=2-3, n=1-5,16-20 of a nonrotating star, and followed these modes by increasing the rotation rate up to 70% of the break-up rotation rate. We then derived perturbative coefficients and determined the domains of validity of the perturbative methods. Second-order perturbative methods are suited to computing low-order, low-degree mode frequencies up to rotation speeds ~100 km/s for typical gamma Dor stars or ~150 km/s for B stars. The domains of validity can be extended by a few tens of km/s thanks to the third-order terms. For higher order modes, the domains of validity are noticeably reduced. Moreover, perturbative methods are inefficient for modes with frequencies lower than the Coriolis frequency 2Omega. We interpret this failure as a consequence of a modification in the shape of the resonant cavity that is not taken into account in the perturbative approach.

preprint2010arXiv

On the interpretation of echelle diagrams for solar-like oscillations. Effect of centrifugal distortion

This work aims at determining the impact of slow to moderate rotation on the regular patterns often present in solar-like oscillation spectra. We focus on the well-known asteroseismic diagnostic echelle diagrams, examining how rotation may modify the estimates of the large and small spacings, as well as the identification of modes. We illustrate the work with a real case: the solar-like star $η$ Bootis. The modeling takes into account rotation effects on the equilibrium models through an effective gravity and on the oscillation frequencies through both perturbative and non-perturbative calculations. We compare the results of both type of calculations in the context of the regular spacings (like the small spacings and the scaled small spacings) and echelle diagrams. We show that for echelle diagrams the perturbative approach remains valid for rotational velocities up to 40-50 km/s. We show that for the rotational velocities measured in solar-like stars, theoretical oscillation frequencies must be corrected up to the second-order in terms of rotation rate, including near degeneracy effects. For rotational velocities of about 16 km/S and higher, diagnostics on large spacings and on modal identification through echelle diagrams can be significantly altered by the presence of the $m\neq0$ components of the rotationally split modes. We found these effects to be detectable in the observed frequency range. Analysis of the effects of rotation on small spacings and scaled small spacings reveals that these can be of the order of, or even larger than surface effects, typically turbulence, microscopic diffusion, etc. Furthermore, we show that scaled spacings are significantly affected by stellar distortion even for small stellar rotational velocities (from 10-15 km/s) and therefore some care must be taken when using them as indicators for probing deep stellar interiors.

preprint2010arXiv

On the power spectrum of solar surface flows

The aim of this work is to give new observational constraints on solar surface flows by determining the horizontal scale dependence of the velocity and intensity fields, as represented by their power spectra, and to offer some theoretical guidelines to interpret these spectra. We use long time series of images taken by SOT/Hinode and reconstruct both horizontal (by granule tracking) and vertical (by Doppler effect) velocity fields in a field of view 75x75Mm^2. At small sub-granulation scales, the kinetic energy spectral density associated with vertical motions exhibits a k^{-10/3}-like spectrum, while the intensity fluctuation spectrum follows a k^{-3} or k^{-17/3}-like spectrum at the two continuum levels investigated (525 and 450 nm respectively). We discuss the physical origin of these scalings and argue that they provide a direct observational signature of buoyancy-driven turbulent dynamics in a strongly thermally diffusive regime. In the mesogranulation range and up to a scale of 25Mm, we find that the vertical velocity field amplitude decreases like L^{-3/2} with the horizontal scale L. This behaviour corresponds to a k^2 spectral power law. Still in the mesoscale range, [2.5, 10]Mm, we find that intensity fluctuations in the blue continuum follow a power law in k^2. We show that granule tracking cannot sample scales below 2.5Mm. We locate the supergranulation energy peak around 30Mm and show that the emergence of a pore erases this spectral peak. Thanks to a scale height estimate, we find that supergranular flows are shallow. (abridged abstract).

preprint2009arXiv

Viscous dissipation by tidally forced inertial modes in a rotating spherical shell

We investigate the properties of forced inertial modes of a rotating fluid inside a spherical shell. Our forcing is tidal like, but its main property is that it is on the large scales. Our solutions first confirm some analytical results obtained on a two-dimensional model by Ogilvie (2005). We also note that as the frequency of the forcing varies, the dissipation varies drastically if the Ekman number E is low (as is usually the case). We then investigate the three-dimensional case and compare the results to the foregoing model. These solutions show, like their 2D counterpart, a spiky dissipation curve when the frequency of the forcing is varied; they also display small frequency intervals where the viscous dissipation is independent of viscosity. However, we show that the response of the fluid in these frequency intervals is crucially dominated by the shear layer that is emitted at the critical latitude on the inner sphere. The asymptotic regime is reached when an attractor has been excited by this shear layer. This property is not shared by the two-dimensional model. Finally, resonances of the three-dimensional model correspond to some selected least-damped eigenmodes. Unlike their two-dimensional counter parts these modes are not associated with simple attractors; instead, they show up in frequency intervals with a weakly contracting web of characteristics. Besides, we show that the inner core is negligible when its relative radius is less than the critical value 0.4E^{1/5}. For these spherical shells, the full sphere solutions give a good approximation of the flows (abridged abstract).

preprint2002arXiv

A simulation of solar convection at supergranulation scale

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$ω$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations.