Gravito-inertial modes in a differentially rotating spherical shell

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.

Gravito-inertial modes in a differentially rotating spherical shell

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.

Numerical convergence in self-gravitating disc simulations: initial conditions and edge effects

We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $Omega^{-1}$, where $Omega$ is the local angular velocity.