simple_snII

simple_snII

This project computes simple models of the passage of a shock through the layers of a massive star and the attendant nucleosynthesis. This project was introduced in the paper Explosive nucleosynthesis of ¹⁵N in a massive-star model. Instructions for downloading, installing, and running the codes to make the calculations in that paper are here.

Installing the project.

If you have not previously done so, install the necessary [libraries]. Next, if you have not previously done so, create a projects directory:

** mkdir projects **

Change into the projects directory:

** cd projects **

Now you can checkout the project:

** svn checkout svn://svn.code.sf.net/p/nucnet-projects/code/trunk/projects/simple_snII**

Compiling the executables.

Change into the project directory:

** cd simple_snII **

and compile the codes:

** ./project_make **

The project compiled successfully if you see Make succeeded. Otherwise, please submit a Ticket or make a post to the Discussion board (see links above).

Running the simple explosion code.

Create a presupernova model xml file. For illustration, we will use the ../nucnet-tools-code/data_pub/s15a28a29_presn.xml created in the [rauscher_xml] tutorial. To run the explosion code for an explosion energy of 1.0 B (1 B = 1 Bethe = 1.0 x 10⁵¹ ergs), type:

** ./simple_snII ../nucnet-tools-code/data_pub/s15a28a29_presn.xml 1.0e51 expl_1.0e51.xml**

The file expl_1.0e51.xml is the output. Properties in this file may be examined with the codes in the ../nucnet-tools-code/examples/analysis directory. The abundances in this file are still the pre-supernova abundances.

Running the single-zone nucleosynthesis code.

Get the data to run the nucleosynthesis codes by typing

** make simple_snII_data**

Now run the single-zone code for a particular zone. For example, to run the code for the zone with label1 = 600, type:

** ./run_single_zone expl_1.0e51.xml ../nucnet-tools-code/data_pub/simple_snII__zone.xml "[@label1 = '600']" my_output_600.xml**

You can analyze the output file my_output_600.xml with the analysis codes as described in this blog post. The zones in this output file are the time dumps during the calculation.

Running the multi-zone nucleosynthesis code.

If you have not already done so, get the data for the nucleosynthesis calculations:

** make simple_snII_data**

Now run the multi-zone code for all zones by typing:

** ./run_multiple_zone_omp expl_1.0e51.xml ../nucnet-tools-code/data_pub/simple_snII__zone.xml my_output_1.e51.xml ""**

The zones in the output file my_output_1.0e51.xml contain the final mass fractions of each zone at the end of the calculation. You can use the nucnet-tools-code/examples/analysis codes to study the output, but remember the zones now represent the different spatial zones in the supernova, not time dumps.

You can run a subset of zones by choosing the appropriate XPath expression. For example, to run calculations only for zones with interior mass greater than 5.0 Solar masses and less than 7.0 Solar masses (one Solar mass is 1.9891e33 grams), type:

** ./run_multiple_zone_omp expl_1.0e51.xml ../nucnet-tools-code/data_pub/simple_snII__zone.xml my_output_1.e51.xml "[optional_properties[property[@name = 'mass below'] > (5.0 * 1.9891e33) and property[@name = 'mass below'] < (7.0 * 1.9891e33)]]"**

Project Members:

• Bradley S. Meyer (admin)