OutflowWind: new set up for Co-rotating

small window with new parameters:

http://www.pas.rochester.edu/~bliu/OutflowWind/ForAdam/newSetup/newSetup.png

Tiny stellar wind density to check the planetary wind close to the planet surface.

​movie for 40 zones per radii

2. No rotation omega=0

​density omega=0; 20 zones per radii

​temperature omega=0; 20 zones per radii

​density omega=0; lambda=10

​temperature omega=0; lambda=10

  1. No rotation omega=0 High Res longer

Set stellar wind with very low density and run with high resolution. To check if it's possible to reproduce the planetary wind in 2.5D.

http://www.pas.rochester.edu/~bliu/OutflowWind/ForAdam/newSetup/rhoV_rot0000.png http://www.pas.rochester.edu/~bliu/OutflowWind/ForAdam//Finals/lambda5.3_gamma1.01_ns1e-5_d4.png

​density omega=0; 40-80 zones per radii
​Temperature omega=0; 40-80 zones per radii
​pressure omega=0; 40-80 zones per radii


  1. Line plots of density, temperature, Mach number and on dayside (from left boundary to critical radius)


Same run data as 3 above.

http://www.pas.rochester.edu/~bliu/OutflowWind/ForAdam/newSetup/plotCompare2dn3d_lambda5_highres.png

3D: 

gamma=1.01
3dpressure=(gamma-1)*(E-0.5*(px^2+py^2+pz^2)/rho)
vScale=sqrt(vx^2+vy^2+vz^2)
soundSpeed=sqrt(gamma*3dpressure/rho)
mach=vScale/soundSpeed

2D: 

gamma=1.01
2dpressure=(gamma-1)*(E-0.5*(px^2+py^2)/rho)
soundSpeed=sqrt(gamma*2dpressure/rho)
mach=sqrt(vx^2+vy^2)/2dsoundSpeed

Comments

No comments.