Shocked Clump (2D)

Am investigating the interaction between a cold clump in a supernovae remnant and a shock. To model this, am injecting a supersonic wind into the boundary, whose properties are specified by the jump conditions across an adiabatic shock. The jump conditions are calculated using the following parameters in the ambient medium:

T_{amb}	104 K
n_{amb}	104 cm-3
v_{amb}	0

The sound speed of the ambient material for these parameters is (for \gamma=5/3):

C_{amb} (for \gamma=5/3)

which for a injected wind speed of 2,000 km/s gives a Mach number of:

M (for v_{wind}= 2,000 ~km~s^{-1})

The wind parameters are then given by the jump conditions for an adiabatic shock:

T_{wind}	[2\gamma (\gamma-1)/(\gamma+1) ]M^2 =
n_{wind}	4\times 10^4 ~cm^{-3}
v_{wind}	2,000 ~km~s^{-1}

Q - shocked gas speed = shock speed?

The clump is initially in pressure equilibrium with the ambient medium, and its properties are given by:

T_{clump}	102 K
n_{clump}	106 cm-3
v_{clump}	0

The thermodynamics in the simulation are governed by an ideal EOS, and an adiabatic index of gamma = 5/3. The abundances are taken to be solar (\mu=1.3). Other parameters of interest are the clump radius, clump crushing time, sound crossing time (ambient and clump medium to travel a clump radius), wind crossing time (to travel clump radius):

R_{clump}	10^{16} ~cm
t_{cc}	5\times 10^8~s
t_{sc} (for ambient)	5\times 10^8~s
t_{sc} (for )	5\times 10^8~s
t_{wc} (for domain of L=10~R_c)	5\times 10^8~s

The domain dimensions, simulation time & framerate, are given by:

L_{x} = L_{y}	10~R_c
t_{sim}	5\times 10^8~s
dt	5\times 10^8~s