Main menu

2D shock along x cooling test

Note: the plots are in log-linear so sharp dips are from negative values.

This plot shows the steady state solution in black for a shock with a Mach number of 8.56, where the shock is moving in the x direction. The cooling curve is a simple power law with an exponent of 2 and a coefficient of 1.23786e-34. The steady state is put into the code as initial conditions with a resolution of 100 zones across the 2.5e17 cm. The red curve shows the state of the shock after 0.2 mega years. The time step is a factor of 100 smaller than the cooling time and all epsilons were set to 0.1.