Update 02/11

PNe Simulation

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

  • Production runs: initial wind runs to 168/300 frames. Estimated wall time remaining is 10 days, but I've doubled the cores (120) and should take less than that.
  • Need to reserve nodes on Bluehive for the rest (high-/low-momentum, adiabatic and cooling for each).
  • Low-momentum, refinement issue solved by reducing mintracer to 0.0001 (0.1 originally)
  • cooling, issue with tracer density when triggering temperature protection. Should be a quick fix. ​plots
  • Inner boundary conditions, look for zero pressure gradient at the inner radius where outflows start.

Moving on

  • Begin working on Toro's book
  • Pick up recombination

ProtectOutfowTemp subroutine added on Jan 31

SUBROUTINE ProtectOutflowTemp(pos,t,dt,q,Ghost)
  USE GlobalDeclarations
  REAL(KIND=qPREC), DIMENSION(:), INTENT(IN) :: pos
  REAL(KIND=qPREC), DIMENSION(:), INTENT(INOUT) :: q
  REAL(KIND=qPREC), INTENT(IN) :: t, dt
  LOGICAL, INTENT(IN) :: Ghost
  if (q(particle%outflowobj%itracer) > 1e-4*particle%outflowobj%density) then
     q(iE)=max(q(1)*FloorTemp/TempScale, q(iE))
  end if
END SUBROUTINE ProtectOutflowTemp

Comments

No comments.