Using MATLAB, the development of a jet flow downstream of a nozzle is to be plotted. It is approached as an axi-symmetric problem, meaning that the problem can be plotted in 2D. The pressure at the nozzle is higher than the atmosphere, creating expansion fans. These expansion fans propagate and eventually develop into a shock wave. This shockwave formation is shown using the Method of Characteristics. Expansion happens when the characteristic lines diverge and compression happens when characteristic lines converge.
Initial and boundary conditions
M_exit = 2 (Mach number)
P_exit = 2P_atm (Pressure)
Increasing number of characteristics plot
It can be observed that the increasing number of characteristics, creates a more refined separation line (red) between the atmosphere and the jet.
Increasing inlet Mach number plot
It can be observed that the shock point, where all the characteristics overlap, keeps becoming further downstream the higher the mach number becomes.
(M = 1.5, M = 2, M = 2.5, respectively)
A more comprehensive and detailed view on the workout of the problem can be found in the full report: