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 Poincar\'e 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.

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).

An efficient way of solving 2D stability problems in fluid mechanics is to use, after discretization of the equations that cast the problem in the form of a generalized eigenvalue problem, the incomplete Arnoldi-Chebyshev method. This method preserves the banded structure sparsity of matrices of the algebraic eigenvalue problem and thus decreases memory use and CPU-time consumption. The errors that affect computed eigenvalues and eigenvectors are due to the truncation in the discretization and to finite precision in the computation of the discretized problem. In this paper we analyze those two errors and the interplay between them. We use as a test case the two-dimensional eigenvalue problem yielded by the computation of inertial modes in a spherical shell. This problem contains many difficulties that make it a very good test case. It turns out that that single modes (especially most-damped modes i.e. with high spatial frequency) can be very sensitive to round-off errors, even when apparently good spectral convergence is achieved. The influence of round-off errors is analyzed using the spectral portrait technique and by comparison of double precision and extended precision computations. Through the analysis we give practical recipes to control the truncation and round-off errors on eigenvalues and eigenvectors.

Oscillation modes of rapidly rotating stars have not yet been calculated with precision, rotational effects being generally approximated by perturbation methods. We are developing a numerical method able to account for the deformation of the star by the centrifugal force and, as a first step, we determined the acoustic modes of a uniform density spheroid and studied how the spheroid flatness affects these modes.

In the limit of low viscosity, we show that the amplitude of the modes of oscillation of a rotating fluid, namely inertial modes, concentrate along an attractor formed by a periodic orbit of characteristics of the underlying hyperbolic Poincar\'e equation. The dynamics of characteristics is used to elaborate a scenario for the asymptotic behaviour of the eigenmodes and eigenspectrum in the physically relevant r\'egime of very low viscosities which are out of reach numerically. This problem offers a canonical ill-posed Cauchy problem which has applications in other fields.

We investigate the asymptotic properties of inertial modes confined in a spherical shell when viscosity tends to zero. We first consider the mapping made by the characteristics of the hyperbolic equation (Poincar\'e's equation) satisfied by inviscid solutions. Characteristics are straight lines in a meridional section of the shell, and the mapping shows that, generically, these lines converge towards a periodic orbit which acts like an attractor. We then examine the relation between this characteristic path and eigensolutions of the inviscid problem and show that in a purely two-dimensional problem, convergence towards an attractor means that the associated velocity field is not square-integrable. We give arguments which generalize this result to three dimensions. We then consider the viscous problem and show how viscosity transforms singularities into internal shear layers which in general betray an attractor expected at the eigenfrequency of the mode. We find that there are nested layers, the thinnest and most internal layer scaling with $E^{1/3}$-scale, $E$ being the Ekman number. Using an inertial wave packet traveling around an attractor, we give a lower bound on the thickness of shear layers and show how eigenfrequencies can be computed in principle. Finally, we show that as viscosity decreases, eigenfrequencies tend towards a set of values which is not dense in $[0,2\Omega]$, contrary to the case of the full sphere ($\Omega$ is the angular velocity of the system). Hence, our geometrical approach opens the possibility of describing the eigenmodes and eigenvalues for astrophysical/geophysical Ekman numbers ($10^{-10}-10^{-20}$), which are out of reach numerically, and this for a wide class of containers.