Update 11/15

  • Removed extra refinement from MHD simulations and extended boundaries by 0.025 in each direction. Still have subsonic material leaving grid in x-y plane, as well as both subsonic and sub-alfvenic eventually leaving in the z direction. Also, alfven speed is becoming very large in a concentrated area around frame 24 and increasing computation time significantly - not sure if there's anything to do about this.

Definitely going to increase grid size more. Anything else to do about computation time and irregularity of contours?

  • WASP-12b simulations appear to have gone extremely well, on the other hand. Ran a test with new parameters for density, bow shock radius, etc.:
0.036180 ratio of planet radius to orbital separation
0.324198 ratio of stellar radius to orbital separation
0.069149 ratio of Hill radius to orbital separation (from planet)
0.125384 ratio of bow shock radius to orbital separation (from planet)
0.314214 ratio of planet sonic radius to orbital separation (from planet)
0.314214 ratio of coriolis radius to orbital separation (from planet)
87.507134 ratio of densities at bow shock
19.542018 stellar lambda
17.369379 planetary lambda
Orbital separation 0.022930 AU
Mass of Star 1.349999 solar masses
Mass of Planet 1.405393 Jupiter masses
Radius of Star 1.599000 solar radii
Radius of Planet 1.736000 Jupiter radii
Temperature of Star 999999.016096 Kelvin
Temperature of Planet 10009.919196 Kelvin
Density of Star 3.211200e-17 g/cc
Density of Planet 3.211200e-15 g/cc
Orbital period 1.090679 days
Mass loss from Star 1.002167e-18 solar masses per yr
Mass loss from Planet 5.460294e+07 g/s
Mach number of Stellar wind at shock 2.050641e-01
Mach number of Planetary wind at shock 4.368280e-01
predicted bow shock radius in units of lscale 1.255086e-01

Mass loss rates may be a bit low - can increase both densities to get higher, but sims look good:

​Frames

Still a slight spiral feature in the initial stellar wind, before planet is created, but doesn't seem to persist (may just be Coriolis effects - doesn't appear to originate from stellar surface). Ran 9 full orbits - reaches quasi-steady state after about 3 days:

​Movie

Attachments (4)

Comments

No comments.