(group 8)

Members: Rashi (Input Person), Raditya (Analysis Person), Alex B. (Code Person)

1. 1d: steady cloud of gas
2. 1d: hit it with supernova shock / plain wave
3. 1d: add source → turn cloud into wind
4. move to 2d
5. parametric study: What is the kill radius at given explosion energy and solar wind strength?
6. production of the “Battle of the Winds” movie blockbuster

Date: 06/02/2014

1. 1d version of steady cloud of gas created using the following parameters:
   • Pressure : pcloud = 1
   • Density : rcloud = 1
   • Velocity : ucloud = 0
   • Ambient density : ramb = 1e-04
2. Then we created a 1d shock wave and let it hit the cloud using the following parameters:
   • Pressure : pshock = 10000
   • Density : rshock = 0.5
   • Velocity : vshock = 0
3. Having succesfully performed the above two steps, we replced the steady cloud of gas with a source generating outgoing mass flow which basically represents 1d version of the solar wind. For the same, we introduced a velocity of 5 at the center in two oppisite directions and a density of 2.
4. Finally we ended the day with hitting the schock wave onto the solar wind to see the interaction beteen the two.

Date : 06/03/2014

1. All 3 of us installed netcdf and got the Serial version of VH1 running. Alex worked on setting up the 2D code for the solar wind and the supernova shock.
2. Myself and Raditya, figured out the input parameters for the problem using Brian Field's paper on “Supernova Collisions with the Heliosphere”. We tested those parameters on the 1D set-up for the problem, basically we replaced the random numbers we used on 06/02 with real numbers describing the problem in 1D format.
3. The resulting density, pressure and velocity plots in 1D looked reasonably good. So we gave those parameters to Alex to be used into the 2D code.
4. Changed the 2D inputs to more natural units ([L]=1 AU, [T]=1 yr, [M]=10^12 g) because VH-1 would crash when using the cgs system for our problem.
5. Fixed an error in our input densities.
6. Fixed the bug of density accumulating on the symmetry axis (we use cylindrical coordinates) by setting the density of the two innermost gridpoints in each column to the value of the third. This is a quick and dirty fix and unphysical! However, we think the problem is caused by some boundary condition/ghost zone issue that we were not able to resolve after spending several hours on it. We thank Richard for suggesting this fix!
7. Finally running a realistic simulation of the collision, now we are all set for producing a high-res video and a low-res parameter study tomorrow.