A numerical study of shock wave diffraction by a circular cylinder

Abstract

The nonstationary shock wave diffraction patterns generated by a blast wave impinging on a circular cylinder are numerically simulated using a second-order hybrid upwind method for solving the two-dimensional inviscid compressible Euler equations of gasdynamics. The complete diffraction patterns, including the transition from regular to Mach reflection, trajectory of the Mach triple point and the complex shock-on-shock interaction at the wake region resulting from the Mach shocks collision behind the cylinder are reported in detail. Pressure-time history and various contour plots are also included. Comparison between the work of Bryson and Gross (1961) which included both experimental schlieren pictures and theoretical calculations using Whitham's ray-shock theory and results of the present finite difference computation indicate good agreement in every aspect except for some nonideal gas and viscous effects which are not accounted for by the Euler equations.