COMMON ENVELOPE SIMULATIONS

New Work

  • Advanced Run 183 and produced figures and movies
  • Progress on deriving expressions for total energy conservation prescription
  • Plotted convergence test 185 to confirm that changing refinement region center from particle 1 to particle CM has negligible effect on energy conservation (confirmed)
  • Encountered (and solved) particle creation issue that was crashing code (code is now running on stampede2)

Results

  • Rate of energy gain for Run 183 is now quite small and apparently on track for staying below 10% of the old AGB run (Run 164).
  • Run 183 crashed because of particle creation….I created a switch that turns off particle creation.

Figures

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run183/p_mult_164_177_183.png
http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run183/Energy_percent_change_164_177_183.png

Movies

​Run 183 Density
​Run 183 Density with mesh
​Run 183 Pressure profiles (sliced along x-axis in simulation frame of reference)

Next steps

  • Continue Run 183 (and increase speed by using more nodes)
  • Further modify new energy conserving routines to take into account particle-gas and particle-particle energy terms
  • Further automate and improve simulation:
    • predefine refinement radius as a function of time (should vary smoothly rather than in sudden steps)
    • center of refinement shape should transition smoothly in time from particle 1 to particle CM (instead of suddenly)
    • experiment with using true AMR, e.g. refine on density or pressure gradients within refinement shape

Comments

No comments.