The structure of a shock front in atomic hydrogen

We examine the problem of a shock wave propagating in a gravitational field in the presence of pressure and density gradients by attacking the non-linear equations of fluid flow. Our approach is analytical rather than numerical, and we analyze the characteristic equations of a fluid in the presence of gravity with radiative dissipation. Because the radiation field enters the fluid equations in the form of an integral, radiative dissipation may be considered an inhomogeneity which does not affect the characteristic directions. The fluid equations remain hyperbolic and thus are amenable to solution by the standard techniques of gas analysis.We give an equation of path for a shock wave and we enumerate the physical conditions which lead to stability or instability. We find that shock waves are generally unstable in most stellar atmospheres unless they are very weak. The form of the instability is that of a spicule deformation similar to that observed in the upper solar chromosphere.This work was carried out at the Smithsonian-Harvard Astrophysical Observatory and was presented in a thesis to Brandeis University, May 1963.