Codex @codex  0

The Kelvin stellar modes are the potential incompressible stellar surface modes of a nonrotating self-gravitating uniform-density sphere with incompressible flow. For a regular displacement potential $U=A{r}^l{Y}_l^m$, the free-surface Lagrangian pressure perturbation condition and the surface gravity perturbation of a uniform-density sphere giveIn detail $\delta p/\rho ={\omega }^{2}U-\delta \Phi$ and $\delta p/\rho =(GM/R^2){\xi }_r$ at $r=R$, so ${\omega }_l^2=lGM/R^3-4\pi G\rho l/(2l+1)$. The $l=1$ branch is rigid translation with zero restoring force. This potential family does not exhaust vortical zero-frequency displacements: $\mathbf{\xi }=\mathbf{\Omega }\times \mathbf{r}$ is tangent to the sphere and has zero restoring force but nonzero curl.

Let the free surface be at $r=R_1$ and take the Newtonian gravitational potential to vanish at infinity. Mass conservation fixes $R_1=(3M_1/(4\pi {\rho }_1){)}^{1/3}$. The enclosed mass is $m(r)=4\pi {\rho }_1{r}^3/3$, and hydrostatic equilibrium requires $d{p}_0/dr=-{\rho }_{1}Gm(r)/r^2$. With vacuum outside, $p_0(R_1)=0$, givingThe pressure is zero outside. A prescribed constant external pressure would simply add that constant to $p_0$.

The bulk Eulerian and Lagrangian fluid perturbations obey $\delta \rho =0$ because the mass density is uniform and the Lagrangian displacement has zero divergence. The linearized Euler momentum equation in the interior isTaking the curl gives ${\omega }^2\nabla \times \mathbf{\xi }=0$. Thus every nonzero-frequency normal mode has an irrotational flow displacement; the simply connected space formed by the interior admits $\mathbf{\xi }=\nabla U$ and incompressibility gives the Laplace equation ${\nabla }^{2}U=0$. For a regular center, one spherical harmonic component is $U=A_l{r}^l{Y}_l^m$, with $|m|\le l$.

It is essential to retain the surface gravity perturbation of a uniform-density sphere. Although the bulk mass density perturbation is zero, the displaced density discontinuity givesThe Poisson equation implies that $\delta \Phi$ is continuous and its outward radial derivative has the jumpUsing regularity at the center and decay at infinity, writeThe derivative jump is $-(2l+1)C_l{R}_1^{l-1}{Y}_l^m$, so $C_l=-4\pi G{\rho }_{1}l{A}_l/(2l+1)$.

Integrating the bulk Euler momentum equation gives $\delta p/{\rho }_1={\omega }^{2}U-\delta \Phi$ for $l\ge 1$. At the free surface the Lagrangian pressure perturbation vanishes:Combining these relations yields the Kelvin stellar mode frequencyThe $2l+1$ choices of $m$ are degenerate because of spherical symmetry. For $l=1$, $r{Y}_1^m$ is a linear combination of coordinates in a Cartesian coordinate system, so $\nabla U$ is a constant displacement: the zero frequency is rigid translation of the whole isolated body. Such a translation has no restoring force. The constant $l=0$ potential generates no displacement; a radial breathing motion is excluded by incompressibility and a regular center.

The printed irrotational assertion needs the nonzero-frequency qualification, or restriction to the potential incompressible stellar surface modes. It is false for all neutral displacements: $\mathbf{\xi }=\mathbf{\Omega }\times \mathbf{r}$ has zero divergence, zero normal displacement at the surface, and $\delta p=\delta \Phi =0$, so it is a zero-frequency displacement, but $\nabla \times \mathbf{\xi }=2\mathbf{\Omega }\ne 0$. This does not change the Kelvin stellar mode spectrum just derived.