Changes between Version 35 and Version 36 of u/erica/2DShockedClumpsSNR


Ignore:
Timestamp:
01/24/18 06:21:40 (7 years ago)
Author:
Erica Kaminski
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/erica/2DShockedClumpsSNR

    v35 v36  
    11'''Shocked Clump (2D)'''
    22
    3 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:
     3Am 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 (note all values given are approximate; actual values are calculated to machine precision in the attached problem module):
    44
    55|| $T_{amb}$ || $10^4 ~K$ ||
    … …  
    99The sound speed of the ambient material for these parameters is (for $\gamma=5/3$):
    1010
    11 || $C_{amb}$ || $10 ~km~s^{-1}$ ||
     11|| $C_{amb}$ || $10 km~s^{-1}$ ||
    1212
    1313which for a injected wind speed of $v_s=2,000 ~km/s$ gives a Mach number of:
    1414
    15 || $M_{wind}$  || $200$ ||
     15|| $M_{wind}$  || $194$ ||
    1616
    17 The wind parameters are then given by the jump conditions for an adiabatic shock ($\gamma=5/3$):
     17The wind parameters are then given by the R-H jump conditions (see attached pdf for equations used):
    1818
    19 || $T_{wind}$ || $[2\gamma (\gamma-1)/(\gamma+1) ]M^2 = 33,333 ~K $ ||
    20 ||$n_{wind}$ || $4\times 10^4 ~cm^{-3}$ ||
    21 ||$v_{wind}$||$2,000 ~km~s^{-1} $||
     19|| $T_{wind}$ || $1e+8~ K $ ||
     20||$n_{wind}$ || $4e+4 ~cm^{-3}$ ||
     21||$v_{wind}$||$1.5e+8 ~cm~s^{-1} $||
    2222
    2323The clump is initially in pressure equilibrium with the ambient medium, and its properties are given by:
    … …  
    4545|| $dt$ || $10^7~s~=~1/50~t_{cc}$ ||
    4646
    47 === Run 1 - Pulse Case ===
     47=== Run 1 - Continuously injected supersonic wind case ===
    4848
    49 Supersonic wind is injected at boundary t = 0 only.
     49Supersonic wind is injected each timestep over the course of the simulation.
    5050
    5151Visualizations: