Meeting Update 04/24/2012 - Eddie

I successfully added MHD to the 1D radiative shock simulations. The magnetic field is only in the y-direction (perpendicular to the fluid flow). This was initially trickier than I thought it would be. Having magnetic fields changes the shock jump equations and the ODE within the cooling region. I made a lot of changes to my problem.f90 so that the module as a whole could easily switch between hydro and MHD. As expected, the presence of the magnetic field lengthens the cooling region. Here is a movie that compares the temperature profile for hydro and MHD (both simulations use DM cooling). DM_compare_temp.gif​

I'm still working out some kinks in using the next type of cooling (NEQ cooling). The tricky part here is now I have to create initialization profiles for hydrogen H and ionized hydrogen HII (I'm not even considering molecular hydrogen, or helium and its ionized species). The cooling routines calculate these ionization and recombination rates, and these can be used to update my initial quantities to create an initial profile. These rates have to be self-consistent with the cooling rate that NEQ cooling gives. So the problem essentially becomes solving two different ODEs simultaneously. I think I have this done correctly, and I am now just running into issues with the way the main cooling routine interfaces with NEQ cooling.

Attachments (1)

Download all attachments as: .zip

Comments

No comments.