Update on CE project

New Work

  1. Now able to run code on bluehive 2.5.
  2. Three runs in progress.
  3. How to optimize refinement?
  4. Edited code to be able to write some output to files.
  5. What other data to output?
  6. Reading and discussions conerning accretion onto secondary in CEE (presentation?).

Summary of New Results

  1. Code is running on bh2.5 and seems to be faster than regular bh (which is nice).
  2. New runs will be analyzed next week. Stampede queue time has increased (will eventually go to infinity as it is phased out)! Bluestreak queue is variable. No queue for Bluehive 2.5 as of now!
  3. There are a few possibilities:
    • Refine to maxlevel up to refinement radius of min(5e12cm,1.5*separation) about primary point particle (trying now)?
    • Choose refinement radius by hand at each restart (trying now)?
    • Set some refinement criteria within the refinement radius and force to lowest level outside refinement radius?
    • Refine to maxlevel within refinement radius but make refinement radius very small and use large buffer zones?
  4. Can now write accreted mass & momentum as well as time and acceleration on particles from gas gravity each time step.
  5. Aside from those, other output is also desired, so need to decide what exactly, then spend time implementing that.
    • angular momentum
    • particle separation
    • integrals of damping force out to a given radius around secondary
    • gas potential and kinetic energies
  6. Should I give a presentation at a CE meeting of what has been done so far on this topic of accretion in CEE, and also about what output needs to be generated?

Ongoing runs

Run to test making refinement zone decrease with time:

Binary run 133 similar to run 125 but now the refinement zone changes with time
Relaxation run: 096
First frame: 75 (5 RG freefall times, when velocity damping ended)
Last frame:
Total simulation time:
Machine and partition: Bluestreak standard (running, completed up to frame 181 in 4 days walltime)
Number of cores: 8192
Total wall time:
Hydro BCs: extrapolated
Poisson BCs: multipole expansion
Box size: L=4e13 cm
Base resolution: 9.0 Rsun (643 cells)
Highest resolution: 0.29 Rsun (20483 cells, 5 levels AMR)
AMR implementation: set by hand to have max level around point particles
Max resolution zone: within min(5e12cm, 1.5*particle_separation)
Buffer zones: 2 cells
Softening length: 2.4 Rsun
Ambient density: 6.7e-9 g/cc
Ambient pressure: 106 dyne/cm2
DeafaultAccretionRoutine=2 (Krumholz)

Run to test case where no accretion onto secondary is permitted:

Binary run 135 similar to run 133 but now DefaultAccretionRoutine=0 instead of 2. Also suppress generation of new sink particles.
Relaxation run: 096
First frame: 75 (5 RG freefall times, when velocity damping ended)
Last frame:
Total simulation time:
Machine and partition: bluehive2.5 standard (running, completed up to frame 196 in 2 days wall time)
Number of cores: 120
Total wall time: about
Hydro BCs: extrapolated
Poisson BCs: multipole expansion
Box size: L=4e13 cm
Base resolution: 9.0 Rsun (643 cells)
Highest resolution: 0.29 Rsun (20483 cells, 5 levels AMR)
AMR implementation: set by hand to have max level around point particles
Max resolution zone: within min(5e12cm, 1.5*particle_separation)
Buffer zones: 2 cells
Softening length: 2.4 Rsun
Ambient density: 6.7e-9 g/cc
Ambient pressure: 106 dyne/cm2
DeafaultAccretionRoutine=0 (No accretion)

Run that uses twice as high resolution and 10x lower ambient pressure as previous runs:

Binary run 132 with double max resolution, lower resolution in ambient medium, 10x smaller ambient pressure than run 116
Relaxation run: 129
First frame: 75 (5 RG freefall times, when velocity damping ended)
Last frame:
Total simulation time:
Machine and partition: Stampde 1 normal (running, completed up to frame 153, or 18 sim-days)
Number of cores: 1024
Total wall time: 4 days so far (starting from frame 75 of relaxation run)
Hydro BCs: extrapolated
Poisson BCs: multipole expansion
Box size: L=4e13 cm
Base resolution: 2.25 Rsun (2563 cells)
Highest resolution: 0.14 Rsun (40963 cells, 4 levels AMR)
AMR implementation: set by hand to have max level around point particles
Max resolution zone: within 5e12 cm (71.87 Rsun) of primary center and within a cylinder of radius 20 Rsun and height 20 Rsun around secondary center. After t~9.6e5s (~frame 123) refinement radius around primary was halved to 2.5e12cm.
Buffer zones: 0 cells (no buffer zones)
Softening length: 2.4 Rsun
Ambient density: 6.7e-9 g/cc
Ambient pressure: 105 dyne/cm2

Next steps

  • Analysis of runs 133, 135, 132 and comparison with previous runs.
  • Decide what data is needed to output and implement this in the code.
  • Improve style of code in preparation for more serious runs.
  • Run combining the merits of 132 and 133 that also outputs the relevant data.

Comments

No comments.