Changes between Version 54 and Version 55 of 1DPulsedJets
- Timestamp:
- 06/19/12 14:36:04 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
1DPulsedJets
v54 v55 73 73 [[BR]] 74 74 75 == Connecting WithNEQ Cooling ==75 == Connecting To NEQ Cooling == 76 76 One component of the cooling rate calculated by NEQ cooling is called metal_loss. This is the energy loss due to metal excitation, and it is based on the ionization fraction and DM cooling curve. Z cooling passes the parameter iCooling to the subroutine Cool_Derivatives. Now, if iCooling==ZCool, metal_loss will be subtracted from the cooling rate, and it is going to be replaced by a different rate. 77 77 … … 278 278 279 279 [[BR]] 280 [[CollapsibleStart(Results)]] 280 281 === Results === 281 282 A run with no cooling was done to make sure the post-shock values were correct. With no cooling (adiabatic), the hydrodynamic quantities should jump discontinuously at the shock, and then remain constant while the shock itself remains stationary. Only the first 100 cells of the simulation are shown here. … … 312 313 || [[Image(Z_temp0010.png, width=400)]] || 313 314 ||= [attachment:Z_temp.gif Temperature Movie] =|| 315 316 [[CollpasibleEnd]] 317 318 [[BR]] 319 === Including MHD === 320 Making these simulations work with MHD was no simple task. Having a magnetic field adds another jump condition, and it makes the equations within the cooling region more complicated. Here are what the equations look like now: