No Arabic abstract
Oscillations have been detected in a variety of stars, including intermediate- and high-mass main sequence stars. While many of these stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study on axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. These modes probe the deep layers of the star around its convective core. We consider a simplified model where the radiative zone of a star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral simulations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We characterise these attractors by the dependence of their Lyapunov exponent with the BV frequency of the background and the oscillation frequency of the mode. Finally, we note that, as modes associated with short-period attractors form dissipative structures, they could play an important role for tidal interactions but should be dismissed in the interpretation of observed oscillation frequencies.
While many intermediate- and high-mass main sequence stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study of axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. We consider a simplified model where the radiative zone of the star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral computations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We describe various geometries for the propagation domains, conditioning the surface visibility of the corresponding modes.
The gravito-inertial waves propagating over a shellular baroclinic flow inside a rotating spherical shell are analysed using the Boussinesq approximation. The wave properties are examined by computing paths of characteristics in the non-dissipative limit, and by solving the full dissipative eigenvalue problem using a high-resolution spectral method. Gravito-inertial waves are found to obey a mixed-type second-order operator and to be often focused around short-period attractors of characteristics or trapped in a wedge formed by turning surfaces and boundaries. We also find eigenmodes that show a weak dependence with respect to viscosity and heat diffusion just like truly regular modes. Some axisymmetric modes are found unstable and likely destabilized by baroclinic instabilities. Similarly, some non-axisymmetric modes that meet a critical layer (or corotation resonance) can turn unstable at sufficiently low diffusivities. In all cases, the instability is driven by the differential rotation. For many modes of the spectrum, neat power laws are found for the dependence of the damping rates with diffusion coefficients, but the theoretical explanation for the exponent values remains elusive in general. The eigenvalue spectrum turns out to be very rich and complex, which lets us suppose an even richer and more complex spectrum for rotating stars or planets that own a differential rotation driven by baroclinicity.
We investigate the asymptotic properties of axisymmetric inertial modes propagating in a spherical shell when viscosity tends to zero. We identify three kinds of eigenmodes whose eigenvalues follow very different laws as the Ekman number $E$ becomes very small. First are modes associated with attractors of characteristics that are made of thin shear layers closely following the periodic orbit traced by the characteristic attractor. Second are modes made of shear layers that connect the critical latitude singularities of the two hemispheres of the inner boundary of the spherical shell. Third are quasi-regular modes associated with the frequency of neutral periodic orbits of characteristics. We thoroughly analyse a subset of attractor modes for which numerical solutions point to an asymptotic law governing the eigenvalues. We show that three length scales proportional to $E^{1/6}$, $E^{1/4}$ and $E^{1/3}$ control the shape of the shear layers that are associated with these modes. These scales point out the key role of the small parameter $E^{1/12}$ in these oscillatory flows. With a simplified model of the viscous Poincare equation, we can give an approximate analytical formula that reproduces the velocity field in such shear layers. Finally, we also present an analysis of the quasi-regular modes whose frequencies are close to $sin(pi/4)$ and explain why a fluid inside a spherical shell cannot respond to any periodic forcing at this frequency when viscosity vanishes.
(abbreviated) In this paper we develop a consistent WKBJ formalism, together with a formal first order perturbation theory for calculating the properties of the inertial modes of a uniformly rotating coreless body (modelled as a polytrope and referred hereafter to as a planet) under the assumption of a spherically symmetric structure. The eigenfrequencies, spatial form of the associated eigenfunctions and other properties we obtained analytically using the WKBJ eigenfunctions are in good agreement with corresponding results obtained by numerical means for a variety of planet models even for global modes with a large scale distribution of perturbed quantities. This indicates that even though they are embedded in a dense spectrum, such modes can be identified and followed as model parameters changed and that first order perturbation theory can be applied. This is used to estimate corrections to the eigenfrequencies as a consequence of the anelastic approximation, which we argue here to be small when the rotation frequency is small. These are compared with simulation results in an accompanying paper with a good agreement between theoretical and numerical results. The results reported here may provide a basis of theoretical investigations of inertial waves in many astrophysical and other applications, where a rotating body can be modelled as a uniformly rotating barotropic object, for which the density has, close to its surface, an approximately power law dependence on distance from the surface.
We propose a numerical method to compute the inertial modes of a container with near-spherical geometry based on the fully spectral discretisation of the angular and radial directions using spherical harmonics and Gegenbauer polynomial expansion respectively. This allows to solve simultaneously the Poincare equation and the no penetration condition as an algebraic polynomial eigenvalue problem. The inertial modes of an exact oblate spheroid are recovered to machine precision using an appropriate set of spheroidal coordinates. We show how other boundaries that deviate slightly from a sphere can be accommodated for with the technique of equivalent spherical boundary and we demonstrate the convergence properties of this approach for the triaxial ellipsoid.