Changes between Version 24 and Version 25 of 1DPulsedJets
- Timestamp:
- 02/17/12 09:41:26 (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
1DPulsedJets
v24 v25 90 90 Global Data: 91 91 Domain length = 30 AU 92 Final time = 0.01(computational units)92 Final time = 1.7 (computational units) 93 93 }}} 94 Note that these runs feature a stationary ambient, and they do not include MHD or jet perturbations. 95 In order to make some predictions based on these initial parameters, we need some equations... 94 Note that these runs feature a stationary ambient, and they do not include MHD or jet perturbations. To get an idea of the form of the shock waves that are produced, here is a lineout of velocity: 95 96 [[Image(vel_sample.png,width=400)]] 97 98 Now in order to make some predictions based on these initial parameters, we need some equations... 96 99 [[BR]] 97 100 ==== Important Equations ==== … … 136 139 A rate of 4.7863x10^-22^ was used, which is fairly consistent with Dalgarno !McCray. So if we have a domain that is 30 AU and we want 10 cells per cooling length, then we would need a resolution of about 94 cells. That doesn't seem too bad, so I will round up to 100 and start from there. 137 140 [[BR]] 138 ==== Simulations ====141 ==== Results ==== 139 142 Here we look at what post-shock temperature the simulation gives as a function of the effective resolution. The effective resolution = (# cells) x 2^(# AMR levels)^. For Tps I will use the highest temperature that I get from the simulation. 140 143 ||= Effective Resolution =||= cells/lcool =||= Tps (10^3^ K) =|| … … 145 148 ||= 1600 =||= 170.85 =||= =|| 146 149 ||= 3200 =||= 341.70 =||= =|| 150 And in a graphical form, the table looks like this: 151 152 [[Image(temp_res.png,width=400)]] 153 154 So the big question is, why does it appear that we need way more resolution (about 100x more) than is predicted from the cooling length?
