Parameter Person - Sherwood Richers

Code Person - Justin Lietz

Analysis Person - Kaitlin Cook

Project Considered: Post-Processing nucleosynthesis

Extract a Lagrangian profile from the blast wave simulation and use the network to calculate the nucleosynthesis.

The plan: First use the Lane-Emden equation to model the density distribution. Use this as input for XNet and evolve at the right temperature for the initial NSE stage.

Modify VH1 for Lagrangian dynamics, then apply the Sedov blast on the initial density distribution (from the same solution to the Lane-Emden equation). This solution will then be used as input to XNet, using a 150 element network.

Interesting outputs: Abundance of elements as a function of time, and after the blast, as a function of radius. It is intended that this will be displayed in video form. Another output of interest is the total energy released/required in nucleosynthesis - is this energy significant enough to effect the hydrodynamics? To what extent?

——Logbook starts here ——-

2/6/14 3pm KJC: Lane-Emden Equation: We solved this in Mathematica, because, well, it only takes two lines. Far more lines to actually give a good output for FORTRAN, but that's ok. The notebook may be found on the git here. Output is three columns - radius (cm), density (g/cm^3) and pressure (dyn/cm^2). A sample output is in the github

Parameters for the Lane-Emden Equation:

ρ_c = 1.6*10^5

μ = 1/2

n = 3/2

γ = n +1/n

2/6/14 4pm KJC: JGL Got VH1 up and running for the Sedov blast as written for the problems done last week. We now have to combine the solution from the Mathematica output into the VH1 scripts. SR is working on this.

2/6/14 4pm JGL: Have netCDF installed correctly, and compiled serial version of VH1 using a modified makefile (github upload to follow). The vh1-serial code outputs an *.nc file, also installed ncview to view this file. With this, we can view the 2D Sod tube solution via the convenient ncview program. For the starter version, I added a python script written by Nathan Parzuchowski to animate the 1D Sedov blast solution by turning each time step data file into a frame for the animation. We will probably stick with 1D solutions, since we want to do our final simulation with the Lagrangian remap which is supposedly more well behaved in 1D. The python script requires python-numpy and python-matplotlib.

2/6/14 6pm JGL: Just added the reaction network code to the git repository. Should be up and running for the use of 150 isotopes now. We plan to eventually run the reaction network code for both just the alpha network as well as the 150 isotope network to compare the results. We may also try other networks if we have the time.

GitHub