Changes between Version 21 and Version 22 of ShockedClumps


Ignore:
Timestamp:
05/24/11 11:31:19 (14 years ago)
Author:
Jonathan
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • ShockedClumps

    v21 v22  
    1111== Run Params ==
    1212
    13 || Run || v,,wind,, (km/s) || [[latex($\chi$)]] || final time (yr) ||  T,,bow,, (K)|| T,,trans,, (K)|| L,,bow,, (R,,c,,) || L,,trans,, (R,,c,,) ||
    14 || A || 50  ||  10   || 568.9 yr || 5.64e+04 || 5.64e+03 || 2.73e-01 || 7.86e+00 ||
    15 || B || 100 ||  10   || 248.4 yr || || || || ||
    16 || C || 200 ||  10   || 142.2 yr || 9.02e+05 || || || ||
    17 || D || 400 ||  10   || 71.1 yr  || 3.61e+06 || 3.61e+05 || 1.88e+03 || 2.11e-01 ||
    18 || E || 800 ||  10   || 35.6 yr  || 1.44e+07 || 1.44e+06 || 1.40e+04 || 2.30e+01 ||
    19 || F || 200 ||  100  || 142.2 yr || 9.02e+05 || 9.02e+04 || 3.93e+01 || 1.41e+00 ||
    20 || G || 200 ||  10   || 142.2 yr ||  9.02e+05 ||
     13|| Run || v,,wind,, (km/s) || [[latex($\chi$)]] || t,,cc,, (yr) || final time (yr) ||  T,,bow,, (K)|| T,,trans,, (K)|| L,,bow,, (R,,c,,) || L,,trans,, (R,,c,,) ||
     14|| A || 50  ||  10   ||  3.00 yr  || 568.9 yr || 5.64e+04 || 5.64e+03 || 2.73e-01 || 7.86e+00 ||
     15|| B || 100 ||  10   ||  1.50 yr || 248.4 yr || 2.26e+04 || 2.26e+03 || 8.10e-01 || 1.16e-01 ||
     16|| C || 200 ||  10   ||  .750 yr || 142.2 yr || 9.02e+05 || 9.02e+04 || 3.93e+01 || 3.34e-02 ||
     17|| D || 400 ||  10   || .375 yr  || 71.1 yr  || 3.61e+06 || 3.61e+05 || 1.88e+03 || 2.11e-01 ||
     18|| E || 800 ||  10   || .1875 yr  || 35.5 yr || 1.44e+07 || 1.44e+06 || 1.40e+04 || 2.30e+01 ||
     19|| F || 200 ||  100  || .750 yr  || 142.2 yr || 9.02e+05 || 9.02e+04 || 3.93e+01 || 1.41e+00 ||
     20|| G || 200 ||  10   || .750 yr  || 142.2 yr ||  9.02e+05 || 9.02e+04 || inf || inf ||
     21
    2122== Shock Structure ==
    2223
    … …  
    4950[[latex($T_{trans}=T_{bow}/\chi$)]]
    5051
     52 === Cloud crushing time ===
     53  In the strongly cooling case, the forward and reverse wind shock collapse to the center of the shock which is travelling at [[latex($v_{wind}/2$)]].  It makes sense to use this as the wind speed of the 'shock' for calculating the cloud crushing time etc...  Under these assumptions, the shock is first incident on the clump at [[latex($t_0=\frac{100 AU}{v_{wind}/2} = 18.96 yr \left(\frac{50 km/s}{v_wind}\right) = 3.333 \mbox{ frames}$)]].  Furthermore the cloud crushing time is [[latex($t_cc=\frac{\sqrt{\chi}100 AU}{v_{wind}} = \frac{\sqrt{\chi}}{2}t_0 = 6.666 \mbox{ frames } = 1.667 \sqrt{\chi} \mbox{ frames}$)]]
    5154
    5255 == Temperature lineout for non-cooling case ==
    … …  
    5659 == Cooling Lengths ==
    5760
    58 A naive calculation of L,,shock,, using the instantaneous cooling rate at T,,shock,, would implies that run A is bow only, run D is trans only, run E is neither, and that run F is marginally trans only.  However it does appear that run F has both a cooling bow shock and transmitted shock.  To understand the discrepancy we should consider the following figure where we have plotted the cooling lengths in units of clump radius as a function of wind velocity for the bow shock as well as for the transmitted shocks with [[latex($\chi$)]] = 10 & 100
     61[[Image(DMZonePlot250100.png, width=800)]]
     62
     63A naive calculation of L,,shock,, using the instantaneous cooling rate at T,,shock,, would implies that run A is bow only, run B is both, runs C & D are trans only, run E is neither, and that run F is marginally transonly only.  However it does appear that run F has both a cooling bow shock and transmitted shock.  To understand the discrepancy we should consider the following figure where we have plotted the cooling lengths in units of clump radius as a function of wind velocity for the bow shock as well as for the transmitted shocks with [[latex($\chi$)]] = 10 & 100
    5964
    6065Not here we have calculated the cooling length by integrating along the cooling curve until the temperature has dropped by half.  To see the same curve but using the instantaneous cooling rate see [attachment:CoolingLengths.png this image]
    … …  
    97102We know have data from a run at 200 km/s with a density constrast of 10 with and without cooling.  I've replaced the 800 km/s run with the 200 km/s without cooling and the run with a density contrast of 100 with one with a density constrast of 10. 
    98103
     104Below are schlieriens of runs G, A, C, D at frames 7, 11, and 15 at t = .6957, 1.4546, 2.2136 t,,cc,,
     105
    99106[[Image(200nc_50_100_200_7_11_15.png, width=800)]]
    100107