Updates 11/05

Planetary Nebula Simulation

Current status and things to be done

  1. Tuning AMR
    • Refinements triggered on density gradient.
    • Restarted from high-momentum run02 frame with jets out.
    • Ran 5-level AMR (figure below). However, the outflow was turned off.
    • Current run starts from the same frame with 7-level AMR, and outflow turned on. The job is running on 120 nodes with 20-gb memory per node. It has been running for 4 hours (as I'm writing this blog) to generate 1 frame. Need better tuning on qTolerence, or reduce AMR level.
    • Next, start with no jets and then turn on the outflow. See if AMR capture the jets as we want.
  1. Resolution on cavity walls
    • Issues on pressure and density protection. (tbd)
    • Turn on outflow after initial-wind fills the entire simulation domain.
  1. Gamma = 1 (Cooling)
    • Set gamma_7 = 1.001 in the code. (where exactly?)
    • Generate Frame0: Read in SPH results (pressure, kinetic and thermal energy, gamma = 5/3), keep pressure fixed, calculate new thermal energy with gamma = 1.
    • Criteria, say, apply cooling if v > 100 km/s ?

level 5 mesh on 5-level AMR

http://www.pas.rochester.edu/~yzou5/PNe_graphics/amr_mesh_level5_rho.png

temperature

Taken at the beginning of the simulation and just before/after turning on the fast wind.

High momentum run 01 High momentum run 02 Initial Wind
InterpOpts = 0 InterpOpts = 1 InterpOpts = 1 , and no fast wind
http://www.pas.rochester.edu/~yzou5/PNe_graphics/protection_issues/High_Run001/T_all.png http://www.pas.rochester.edu/~yzou5/PNe_graphics/protection_issues/High_Run002/T_all.png http://www.pas.rochester.edu/~yzou5/PNe_graphics/protection_issues/IW/T_all.png

​Movies and ​Graphics

7-level AMR run returns the following error

bhc0074.bluehive.circ.private:UCM:2a62:11a41700: 520618407 us(520618407 us!!!): dapl async_event: DEV ERR 10
bhc0074.bluehive.circ.private:UCM:2a5a:3db57700: 521043336 us(419038 us): dapl async_event: DEV ERR 9

density and pressure protections in physics.data

!================================================================================
! Density Protection Options
!================================================================================
lRestartOnDensityProtections = .false.  ! Do density protections trigger restarts?  [.false.]
lTrackDensityProtections = .false.      ! Track density protections [.false.]
iDensityProtect  = 0                    ! 0-Set to MinDensity, 1-Set to minimum nearby density, [2]-Average nearby densities
iMomentumProtect = 1                    ! 0-Conserve momentum, 1-Set to zero, [2]-Average nearby velocities
MinDensity       = 1d-15                ! Minimum computational density before protection is triggered [1d-10]

!================================================================================
! Pressure Protection Options
!================================================================================
lRestartOnPressureProtections = .false. ! Do pressure protections trigger restarts? [.false.]
lTrackPressureProtections = .false.     ! Track pressure protections [.false.]
iPressureProtect = 0                    ! 0-Set to MinTemp, 1-Set to minimim nearby pressure, [2]-Set to average nearby pressure, 3-Set to minimum nearby temperature, 4-Set to average nearby temperature
MinTemp          = 1d-10                ! [1d-10] minimum allowed temperature for the system in Kelvin before protection is triggered

!================================================================================


Comments

No comments.