Meeting Update 02/11/2013 - Eddie

Pulsed Jets

I altered my 2.5D jet module a bit and added forced refinement of the injection region (it's at 16 cells per jet radius). I completed 4 runs: beta = inf (hydro), 5, 1, 0.4…need to visualize the data now.

Rad Shocks

Spent almost all weekend practically rewriting my 1D rad shock module because I found that at longer run times the shocks were not steady. I rewrote the jump conditions in a different way, and was much more explicit about tracking the electron number density. I also started accounting for electron mass. For a long time I couldn't get it to work, but I just figured out this morning that it was due to the velocity scale. My density was being scaled much more accurately now, but in physics_control.f90 VelScale ~ SQRT(1/mu). In all other scale factors (except rScale), mu cancels out so I actually don't need to multiply the temperature by mu like I originally thought in cooling.f90. I still need to run some tests, but the shock looks rock steady now.

Comments

No comments.