Version 3 (modified by 6 years ago) ( diff ) | ,
---|
Summary of HD209458b parameters
Simulation domain: 20 R_p cube (0.0146 AU)
rho_p = 1.625d-15 g/cm3
T_p = 3000 K (not really important as long as wind gets hotter - T_wind should be >5300 K for this planet)
M_p = 0.73 MJ
R_p = 1.529 RJ
a = 65 R_p = 0.0475 AU
M_s = 1.23 MS (for tidal forces)
F_UV = 2d13 phot/cm2/s (Sun-like for orbital radius)
F_Ly-alpha = 4.1d14 phot/cm2/s, 4.1d15 phot/cm2/s, 4.1d13 phot/cm2/s (HD209458b flux, 10x, 1/10x - would it be better to do 5x and 1/5x?)
(ramping up in approximately 1.5 CT, or 1/6 orbit)
HD209458b results
Steady state of ionizing-only simulation:
Fairly different from the steady state with our inflated planet, but still see that neutral tail from backside (wind is no longer swept laterally around the planet, however, so tail isn't shocked).
Radiation pressure first few frames:
High
Med
MHD planets
Results from larger local runs:
Beta of 0.1 corresponds to surface magnetic field of 0.23 G.
Late state of the beta=0.1 run:
Late state of Beta = 10 run:
Starting global simulation.