Documents and Papers

1.​Stone&Proga 2009

Simulation and Restuls

4.15.2015

1. Parameters

Check blog:afrank04152015 for details

Hot Jupiter

M_{p} = M_j

R_{p} = R_j

n_{p} = m_h

T_{p0} = 10^4

Stellar Wind

V_{s} = 100 km/s

Other parameters.

Resolution: 60 zone/planet radius -check

Allow planetary wind to fill grid for 2 of its own crossing times ( L_grid/V_pw(max) ) then turn on stellar wind and let it run for 10 crossing times ( L_grid/V_s )

2. Check lambda — no stellar wind

Rp=Outflow radius + thickness !! where the radius defines the location of boundary and thickness for wrapping up the velocity

ns/np=1e-4 gamma = 1.01;

lambda=2.6	Mp = 0.5 Mj	Rp = 2Rj	Tp=2e4	namb=1e-4	​low res movie	​high res movie	; ​zoomed high res movie
lambda=5.3	Mp = 0.5 Mj	Rp = 2Rj	Tp=1e4	namb=1e-4	​low res movie	​high res movie
lambda=5.3	Mp = 0.5 Mj	Rp = 2Rj	Tp=1e4	namb=1e-5		​high res movie	;;​zoomed in high res movie
lambda=5.3	Mp = 0.5 Mj	Rp = 2Rj	Tp=1e4	namb=1e-6		​high res movie
lambda=5.3	Mp = 0.5 Mj	Rp = 2Rj	Tp=1e4	namb=1e-8		​high res movie
lambda=5.3	Mp = 1 Mj	Rp = 4Rj	Tp=1e4	namb=1e-4	​low res movie
lambda=10.6	Mp = 1 Mj	Rp = 2Rj	Tp=1e4	namb=1e-4	​low res movie	​high res movie
lambda=10.6	Mp = 0.5 Mj	Rp = 1Rj	Tp=1e4	namb=1e-4	​low res movie
lambda=10.6	Mp = 0.5 Mj	Rp = 2Rj	Tp=5e3	namb=1e-6			Tamb=1000; ​zoomed in movie

3. Backflow check for different line width of line plot — positive value means back flow

Flux=\int_{0}^{L} 2\pi x\rho(x) v_{y}(x) dx Where L is the width of line plot for the integrate.

L= 0.5 Rp at y=1.5 Rp
L = 1 Rp at y=1.5 Rp

4. Runs and Results

Mp=0.5; Rp=2; lambda=5.3

#	gamma	Solar Wind Mach #	n_s/n_p
1	1.01	0.8	10-4	​low res 2.5D	​high res 2.5D
2	1.01	0.8	10-4
3	1.3	0.8	10-4	​low res 2.5D namb=1e-4;​low res 2.5D namb=1e-5;​low res 2.5D namb=1e-8;	transfering
4	1.3	0.8	10-4
5	1.01	5.0	10-4	​low res 2.5D	​high res 2.5D
6	1.01	5.0	10-4
7	1.01	5.0	10-5	​low res 2.5D;​low res 2.5D large	​high res 2.5D
8	1.01	5.0	10-3	​low res 2.5D namb=1e-5	to be analyzed
9	1.3	5.0	10-4	​low res 2.5D namb=1e-5;​low res 2.5D namb=1e-6	​high res 2.5D
10	1.01	0.8	10-3	​low res 2.5D namb=1e-5;	analyzing

5. New runs

blog:afrank04272015

Mp=0.5Mj; Vs=100km/s; namb=1e-6; lambda=5.3

#	gamma	Solar Wind Mach #	n_s/n_p	Tp	Ts	low-res results	high-res results
1	1.01	0.8	10-5	1E4 K	1.86E+6 K	​low res 2.5D Tamb=0.01K small box; ​low res 2.5D Tamb=0.01K large box;​low res 2.5D Tamb=1000K;	running on Stampede 138/200
2	1.01	5.0	10-5	1E4 K	47632.6 K		​high res 2.5D Tamb=0.01K;
3	1.3	0.8	10-5	1E4 K	1.45E+6 K	​low res 2.5D, Tamb=0.01K;​low res 2.5D, Tamb=1000K;	​high res 2.5D Tamb=1000K;
4	1.3	5.0	10-5	1E4 K	37006.8 K	​low res 2.5D Tamb=0.01K; ​low res 2.5D Tamb=1000K;	​Zoomed high res 2.5D Tamb=1000K;

6. 3D runs

Mp=0.5Mj; Vs=100km/s; namb=1e-6; lambda=5.3

#	gamma	Solar Wind Mach #	n_s/n_p	Tp	Ts	low-res results	high-res results
1	1.01	0.8	10-5	1E4 K	1.86E+6 K	​low res 3D rho;​low res 3D T; ​low res 3D vol rendering;	queued on Stampede
2	1.01	5.0	10-5	1E4 K	47632.6 K	​low res 3D rho;​low res 3D T; ​low res 3D vol rendering;	185/300; queued on stampede; Analyzing the 185 frames
3	1.3	0.8	10-5	1E4 K	1.45E+6 K	to be analyzed?	queued on stampede
4	1.3	5.0	10-5	1E4 K	37006.8 K	to be analyzed?	queued on stampede; cut the dimension to half to save time

7. Rotating runs

Mp=0.5Mj; Vs=100km/s; namb=1e-6; lambda=5.3; distance between planet and mass center ~0.1AU, orbital angular velocity of the planet 0.5

#	gamma	Solar Wind Mach #	n_s/n_p	Tp	Ts	low-res results	high-res results
1	1.01	0.8	10-5	1E4 K	1.86E+6 K
2	1.01	5.0	10-5	1E4 K	47632.6 K	running on standard queue 98/150
3	1.3	0.8	10-5	1E4 K	1.45E+6 K
4	1.3	5.0	10-5	1E4 K	37006.8 K