Nonlinear spherical sound waves at the surface of a perfect fluid star

Current numerical relativity codes model neutron star matter as a perfect fluid, with an unphysical "atmosphere" surrounding the star to avoid the breakdown of the equations at the fluid-vacuum interface at the surface of the star. To design numerical methods that do not require an unphysical atmosphere, it is useful to know what a generic sound wave looks near the surface. After a review of relevant mathematical methods, I will present results for low (finite) amplitude waves that remain smooth and, perhaps, for high amplitude waves that form a shock.