Updates 11/12

Planetary Nebula Simulation

  1. Protection issue
    • MinDensity = 1d-100 in density protection
    • Regenerated frame 0, running the initial wind phase
    • Now temperature drops below the MinTemp(=1d-10) in pressure protection
    • See ​plots

http://www.pas.rochester.edu/~yzou5/PNe_graphics/LowerMinDensity_Nov12/Temp/IW_T0000.png http://www.pas.rochester.edu/~yzou5/PNe_graphics/LowerMinDensity_Nov12/density/IW_rho0000.png

Mach number (did something wrong?)

http://www.pas.rochester.edu/~yzou5/PNe_graphics/LowerMinDensity_Nov12/Mach/IW_M0000.png http://www.pas.rochester.edu/~yzou5/PNe_graphics/LowerMinDensity_Nov12/Mach/IW_M0001.png

  1. Cooling (Gamma = 1)
    • Need to learn terms in a(5,128) and b(5,192) matrices in problem modules.
    • Read in quantities from SPH output and mapped to Info%q(i,j,k,1-5). See quoted code below.
    • a,b matrices: 1-density, 2-v_x, 3-v_y, 4-v_z, 5-not sure
    • Info%q: 1-density, 2-p_x, 3-p_y, 4-p_z, 5-energy density
    • Missing gamma? 

select case(info%level)
    case(0)
        do k=1,128
           do j=1,128
              read(195,*),a(1,:)
              read(196,*),a(2,:)
              read(197,*),a(3,:)
              read(198,*),a(4,:)
              read(199,*),a(5,:)
              if (j.ge.57.and.j.le.72) then
                  do i=1,16
                      Info%q(i,j-56,k,1)=a(1,56+i)/rscale
                      Info%q(i,j-56,k,2)=a(1,56+i)*a(2,56+i)/rscale/velscale
                      Info%q(i,j-56,k,3)=a(1,56+i)*a(3,56+i)/rscale/velscale
                      Info%q(i,j-56,k,4)=a(1,56+i)*a(4,56+i)/rscale/velscale
                      Info%q(i,j-56,k,5)=a(1,56+i)*(a(5,56+i)+0.5*(a(2,56+i)&
                      *a(2,56+i)+a(3,56+i)*a(3,56+i)+a(4,56+i)*a(4,56+i)))/pscale
                  end do
              end if
           end do
        end do

  1. AMR
    • Trying on different qTolerance

Comments

No comments.