We investigate whether tidal forcing can result in sound waves steepening into shocks at the surface of a star. To model the sound waves and shocks, we consider adiabatic non-spherical perturbations of a Newtonian perfect fluid star. Because tidal forcing of sounds waves is naturally treated with linear theory, but the formation of shocks is necessarily nonlinear, we consider the perturbations in two regimes. In most of the interior, where tidal forcing dominates, we treat the perturbations as linear, while in a thin layer near the surface we treat them in full nonlinearity but in the approximation of plane symmetry, fixed gravitational field and a barotropic equation of state. Using a hodograph transformation, this nonlinear regime is also described by a linear equation. We show that the two regimes can be matched to give rise to a single mode equation which is linear but models nonlinearity in the outer layers. This can then be used to obtain an estimate for the critical mode amplitude at which a shock forms near the surface. As an application, we consider the tidal waves raised by the companion in an irrotational binary system in circular orbit. We find that shocks form at the same orbital separation where Roche lobe overflow occurs, and so shock formation is unlikely to occur.

We investigate the run-up of a shock wave from inside to the surface of a perfect fluid star in equilibrium and bounded by vacuum. Near the surface we approximate the fluid motion as plane-symmetric and the gravitational field as constant. We consider the "hot" equation of state $P=(\Gamma-1)\rho e$ and its "cold" (fixed entropy, barotropic) form $P=K_0\rho^\Gamma$ (the latter does not allow for shock heating). We find numerically that the evolution of generic initial data approaches universal similarity solutions sufficiently near the surface, and we construct these similarity solutions explicitly. The two equations of state show very different behaviour, because shock heating becomes the dominant effect when it is allowed. In the barotropic case, the fluid velocity behind the shock approaches a constant value, while the density behind the shock approaches a power law in space, as the shock approaches the surface. In the hot case with shock heating,the density jumps by a constant factor through the shock, while the sound speed and fluid velocity behind the shock diverge in a whiplash effect. We tabulate the similarity exponents as a function of the equation of state parameter $\Gamma$and the stratification index $n_*$.

We are interested in the generic behaviour of nonlinear sound waves as they approach the surface of a star, here assumed to have the polytropic equation of state $P=K\rho^\Gamma$. Restricting to spherical symmetry, and considering only the region near the surface, we generalise the methods of Carrier and Greenspan (1958) for the shallow water equations on a sloping beach to this problem. We give a semi-quantitative criterion for a shock to form near the surface during the evolution of generic initial data with support away from the surface. We show that in smooth solutions the velocity and the square of the sound speed remain regular functions of Eulerian radius at the surface.

Problem with the figures should be corrected. Apparently a broken uuencoder was the cause.

Problem with the figures should be corrected. Apparently a broken uuencoder was the cause.