Secondary Sonic Boom Ray Path Prediction using sBOOM

Sonic boom ray path equations are solved numerically to show the extent of secondary booms compared to primary sonic boom carpet. Secondary booms are shown to reach altitudes of ~125 km before turning toward the ground and impinging at lateral distances of ~50 km and beyond. Exact closed-form solutions of the ray paths and their partial derivatives with respect to certain relevant parameters are derived for the simplified case of windless atmospheres and compared against numerical solutions. Ray tube areas are computed using the Jacobian of a single ray as opposed to the traditional technique of using 4 rays that are temporally and azimuthally perturbed. The ray tube areas corresponding to secondary booms can get 1000 times larger than those from primary booms. The numerically computed ray tube areas are verified using exact expressions for the special case of windless atmospheres. Finally, pressure scaling due to ray tube area changes and atmospheric stratification is presented and the shortcomings associated with the current atmospheric interpolation techniques are discussed. It is determined that piece-wise linear interpolation of the atmospheric temperature may be inadequate and could lead to numerical noise in the prediction of secondary booms.