wiki:PlanetaryAtmospheres

BackLinksMenu()

Defining the parameter space for stellar-planetary wind interactions

Assuming circular orbits, we can define the problem using these 11 primary variables

Mass of planet
Radius of planet
Temperature at planet surface
Density at planet surface
Magnetic field strength of planet
Mass of star
Radius of star
Temperature at stellar surface
Density at stellar surface
Magnetic field of star
orbital separation

Now the Euler equations with 2 body potential are as follows.

If we consider the dimensionless form for the equations, we find that the dimensionless equations are unchanged under arbitrary changes in density, length, and time, as long as the initial conditions have the same mach number, Alfven number, and as long as the strength of the potential scales with the velocity2, and the shape of the potential scales with the length scale.

This implies that must be constant, must scale with the length scale, and that must scale with the inverse time scale.

Time and length similarities allow us to fix the total mass and separation without loss of generality. In addition, the actual densities don't matter - just their ratios, so we can also fix the planetary density without loss of generality. So we can reduce the list of 11 primary variables to the following eight dimensionless variables that define the interaction

mass ratio
ratio of densities at surfaces
dimensionless planetary radius
dimensionless stellar radius
characterizes planetary wind
characterizes stellar wind
beta of planet
beta of star

Now instead of those 8, we may want to define the following 5 length scales, and density ratio, and magnetic/kinetic ratios at the bow shock

Ratio of Hill radius to orbital radius
Ratio bow shock radius to orbital radius
Ratio of planetary radius to orbital radius
Ratio of sonic radius of planetary wind to planetary orbital radius
Density ratio at bow shock.
Ratio of stellar radius to orbital radius
Ratio of planetary magnetic pressure to ram pressure (plus thermal) at bow shock
Ratio of stellar magnetic pressure to ram pressure (plus thermal) at bow shock

Using the following relations,

orbital angular velocity
Hill radius
stellar sound speed
planetary sound speed
bow shock standoff distance

and the dimensionless solution to the Parker Wind

which gives us the following solutions for the stellar and planetary winds

We can directly calculate

and numerically solve the following 2 equations for and

Matsakos et al, compare the ordering of the Hill radius (Rt), the wind bow radius Rw (from ram pressure only), and the magnetic bow radius RM (from magnetic pressure) which give 6 different possible orderings. They lump them into 4 different types - based on 3 length scales.

I
II
III
III
IV
IV

Assuming that the planetary radius is always the smallest, there are still 4 free parameters to the problem. They did not essentially explore

  1. Where the location of the stellar sonic surface is in relation to the other length scales
  2. Where the location of the planetary sonic surcface is in relation to the other length scales (they mention that the planetary thermal pressure is assumed to be negligible
  3. The density ratio at the shock
  4. How magnetized the stellar wind is

They also did not explore the role of eccentricity or of the non-isotropic wind from the planet.

A few comments

  1. The bow shock radius will always be larger than the magnetic radius and the ram pressure radius, and perhaps it makes more sense to talk about the bow shock radius and the degree of magnetization then to talk about the 'magnetic' and 'ram pressure' radii as independent constructs.
  2. Also, note that while looks like a function of two variables, it is really a function of one variable

And if we set the mass ratio , and the ratio of the sonic surfaces by setting and , we fix the ratio of sound speeds, and the ratio of velocities coming into the bow shock. But we are also setting the ratio of densities at the bow shock as well as the location of the bow shock . These leads to an overdetermined system with no way to solve for the actual speeds.

So Instead of setting we can set either , or .

Sample params

0.005000
0.200000
0.200000
0.250000
0.010000
1.000000
0.100000
0.100000
Orbital separation 0.100000 AU
Mass of Star 0.997333 solar masses
Mass of Planet 2.795926 Jupiter masses
Radius of Star 4.301941 solar radii
Radius of Planet 1.046271 Jupiter radii
Temperature of Star 431184.397104 Kelvin
Temperature of Planet 143479.010318 Kelvin
Density of Star 3.996761e-07
Density of Planet 1.000000e-08
Orbital period 0.031615 years
Magnetic field of Star 0.000036 Gauss
Magnetic field of Planet 0.084617 Gauss

Matlab Code

See the 3 .m files attached

Planetary Atmospheres

Profiles

  • Density

  • Enclosed Mass

  • Pressure - and rho*R*T mismatch

Module supports

  • Global simulation in a fixed frame
  • Global simulation in a rotating frame
  • Local simulation in a rotating frame
  • Spatial based Refinement of planet
  • Stellar envelope in HSE
  • Source code problem.f90
  • Data file problem.data
  • Uses Particles, Ambients, Clumps, and Refinement Objects

Still working on

  • Basic testing of hydrostatic equilibrium for planet
  • Line transfer for stellar heating
  • Second AMR implicit solve may need to be added later (ie Howell and Greenough 2002)

Results

Working on getting stable planetary atmosphere using profile without a core. Turning on characteristic limiting seems to cause numerical artifacts which lead to 'explosion'. See ticket #

Last modified 9 years ago Last modified on 09/02/15 12:53:16

Attachments (13)

Download all attachments as: .zip

Note: See TracWiki for help on using the wiki.