Posts for the month of January 2022

CEE

Computing

  • Extension on Frontera — no news, write an email to Chris Hempel?
  • Parallel HDF5 (currently only being used for post-processing)
    • Bug has been identified and debugging ongoing (Jonathan)

Jet paper

  • Waiting for the OK to submit
  • Do another run? Or two? If so, what?
    • WD (as Run J6) but with flat profile? 10 x higher Mdot? Accretion off for simplicity? — but WD is almost order of magnitude more expensive because of jet speed
    • MS star with flat profile? 10 x higher Mdot as with Run J5? Accretion off for simplicity?

EOS paper

New analysis

  • Did unbound mass calculation for production runs 277 (MESA EoS without radiation) and 282 (gamma 5/3 ideal gas)
    • Results are generally consistent with the previous results from Run 271 (higher ambient density, lower resolution MESA EOS) and Run 143 (old gamma 5/3 ideal gas run from our other papers)
    • However, the difference between MESA EOS and gamma 5/3 runs is a bit larger
    • Also, unbinding remains flatter and goes up more at the end — this must be due to the 7 times decrease in ambient density
    • Note that these are not the final graphs because they include ambient gas, but the overall conclusions are unlikely to change

https://www.pas.rochester.edu/~lchamandy/EoS/Figures/Mass_282_277.png https://www.pas.rochester.edu/~lchamandy/EoS/Figures/Mass_282_277_half.png

  • Did the energy conservation analysis for both runs — we are under ~2% percent change

https://www.pas.rochester.edu/~lchamandy/EoS/Figures/Energy_percent_change_new_282_277.png

  • This is what their separation curves look like

https://www.pas.rochester.edu/~lchamandy/EoS/Figures/p_mult_282_277.png

Status of runs

The plan is to compare 3 runs (277,282,283) which are identical except for the EOS:

  • Run 277: MESA EOS (completed up to frame 247 on Stampede2)
  • Run 282: Ideal gas gamma=5/3 EOS (completed up to frame 250 on Stampede2)
  • Run 283: MESA EOS with recombination energy removed from EOS tables (completed up to frame 109 on Frontera)
  • Run 276: MESA EOS with maxlevel increased by 1 compared to Run 277 (for convergence study — completed up to frame 47 on Frontera)
  • Run 28?: MESA EOS with maxlevel reduced by 1 compared to Run 277 (for convergence study — not yet started on Frontera)

Next steps

  • Mass unbinding analysis of Runs 277 and 282, as above but now with star material only (i.e. excluding the ambient)
  • Extend Runs 277 and 282 — but for how long?
    • Energy conservation
    • Available SUs
    • What more will we learn?
    • Outflow through boundaries
    • Do we trust it as much after reductions in softening radius?
  • Extend Run 283 (but needs Frontera extension)
  • Ionization/recombination analysis as done with Run 271, now for Run 277

Fermi Project Update January 2022

The journal article on overleaf is at: https://www.overleaf.com/3756433924fzhskcnhbkmz


What Else I Need to Do?

  • Results Section: fix labels on parameter sweep (Fig 2). Use the parameter names on the paper rather than in the code
  • Appendix B: currently has the relevant timescales and ratios written down, what I need to do is figure out which 4/5 parameters or ratios uniquely determine the output. Or influence the output the most.
  • Appendix C: currently has density calculations (as function of radius). I need to copy my results showing technology as function of density. Then I need to use both of those to figure out technology as a function of radius.

What Is Left?

  • send out the paper to collaborators for review
  • write introduction, conclusion, and abstract