Changes between Version 14 and Version 15 of u/adebrech


Ignore:
Timestamp:
07/06/16 11:47:41 (9 years ago)
Author:
adebrech
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/adebrech

    v14 v15  
    44= [[adebrech/PlanetSims|Planet Simulations]] =
    55
    6 =
     6= [[adebrech/Papers|Planet Papers]] =
    77
    88= Fluid Approximation of Charge Exchange =
    9 In [http://iopscience.iop.org/article/10.3847/0004-637X/820/1/3/pdf Christie et al.], they calculate the mean free path (sound speed over rate of reaction) for charge exchange and note that it is less than the planetary radius (in regions dominated by the planetary wind) to justify using the fluid approximation, with mixing and exchange primarily at turbulent boundaries. This carries an assumption that the hot and cold hydrogen have the same temperature, but changing the temp only changes the mean free path by a factor of sqrt(2) at most. In [http://iopscience.iop.org/article/10.1088/0004-637X/693/1/23/pdf Murray-Clay et al.], they justify the fluid approximation in general by comparing the scale height H to the mean free path of a particle, with the fluid approximation holding for H > lambda,,mfp,, at the sonic point (in other words, the exobase is above the sonic point). Thus, H > R,,p,, as well, and therefore greater than mfp of charge exchange near the planet.
     9In [http://iopscience.iop.org/article/10.3847/0004-637X/820/1/3/pdf Christie et al.], they calculate the mean free path (sound speed over rate of reaction) for charge exchange and note that it is less than the planetary radius (in regions dominated by the planetary wind) to justify using the fluid approximation, with mixing and exchange primarily at turbulent boundaries. This carries an assumption that the hot and cold hydrogen have the same temperature, but changing the temp only changes the mean free path by a factor of sqrt(2) at most. In [http://iopscience.iop.org/article/10.1088/0004-637X/693/1/23/pdf Murray-Clay et al.], they justify the fluid approximation in general by comparing the scale height H to the mean free path of a particle, with the fluid approximation holding for H > lambda,,mfp,, at the sonic point (in other words, when the exobase is above the sonic point).
    1010
    1111[[Image(ChristieChrgExch.png, width=400)]]
    … …  
    1313= Change in Bow Shock with Magnetic Field =
    1414If sigma,,*,, and sigma,,p,, are equal, the bow shock radius is unchanged with or without magnetic field - ratio of radius to orbital separation, chi,,bow,, = 0.240468. With sigma,,*,, = 1, sigma,,p,, = 0.1, chi,,bow,, = 0.148204; sigma,,*,, = 0.5, sigma,,p,, = 0.1, chi,,bow,, = 0.187300; sigma,,*,, = 0.1, sigma,,p,, = 0.5, chi,,bow,, = 0.302483; sigma,,*,, = 0.1, sigma,,p,, = 1, chi,,bow,, = 0.363674; and with sigma,,*,, = 0.5 and sigma,,p,, = 1, chi,,bow,, = 0.297793 ≈ chi,,Coriolis,,.
    15 
    16 = [http://arxiv.org/pdf/1206.5003v3.pdf Tremblin & Chiang], Computational Charge Exchange =
    17 Followup to 2008 and 2010 studies of charge exchange between planetary and stellar winds, which used Monte Carlo simulations of 'meta-particles' that were computationally obstructed by bow shocks. Here they use the hydrodynamic equations (no magnetism, Coriolis and centrifugal forces, or tidal gravity). A slow stellar wind (130 km/s) was chosen to approximate the solar wind, and the isothermal planetary wind was initialized as 80% ionized, following Murray-Clay et al. The planetary wind incorporated photoionization/recombination and advection. To incorporate charge exchange, the hydrodynamic code was augmented with chemical reaction solvers - 4 equations relate x,,i,,, i=1-4 representing each possible combination of hot or cold and neutral or ionized hydrogen, and n,,H,, with beta, the reaction rate. These equations take reverse exchange into account, so as not to overestimate neutral hydrogen too greatly (still slightly overestimated). x,,i,, is also included in the hydrodynamic equations. The simulations appear to reproduce the observed absorption curves well, with asymmetry between the two sides of the Doppler shift.
    18 
    19 [[Image(simabs.png, width=500)]]
    20 
    21 [[Image(TremblinEqns.png, width=400)]]
    22 
    23 Compare these with the equations used by Christie et al (see below) for incorporating charge exchange:
    24 
    25 [[Image(ChristieEqns.png,width=400)]]
    26 
    27 
    28 = [http://iopscience.iop.org/article/10.3847/0004-637X/820/1/3/pdf Christie] paper =
    29 2.5D spherical simulations of planetary and stellar wind interactions, including charge exchange, were performed. Hydrodynamic simulations were performed, with density fixed at the base of the planetary wind and an inflow boundary condition on one half of the simulation serving to emulate the stellar wind. In addition to charge exchange, advection, photoionization and recombination, and collisional ionization were included. The escape parameter lambda was used to categorize the models; it was found that there were two distinct regimes, with a transition region between. With lambda <= 4 (high planetary temp), the planetary wind becomes transonic before colliding with the stellar wind, creating a large tail that takes a significant amount of time to mix. With lambda >= 6 (low planetary temp), the planetary wind has no chance to become transonic before it encounters the stellar wind, and the winds mix turbulently rather than collide, resulting in a well-mixed, barely evident tail. The transition region between these is also shown clearly in the calculated mass-loss rates of the simulations.
    30 
    31 = Schneiter paper (2016) =
    32 Paper makes synthetic observations of Lyman-alpha absorption in tails created by interacting solar and planetary winds. Simulations are performed in 3d spherical coordinates, using the Guacho hydrodynamics code, with photoionization of hydrogen included (no magnetism). They have nineteen models of varying stellar UV flux (photoionization rate), stellar wind conditions, and the mass-loss rate of the planet. The planetary wind is initialized self-consistently in order to obtain the desired mass-loss rate, at 3R,,p,,. Both the stellar and planetary winds are isotropic, ignoring effects of tidal locking and atmospheric mixing. They approximate the radiation pressure from the star by reducing stellar gravity.
    33 
    34 They find that by including photoionization, a smaller neutral tail is formed, leading to less absorption; they also find a lower time to a stationary state than in previous models without photoionization. By numerically integrating to determine the optical depth, it is seen that the most absorption in in the blue-shifted side, between -130 and -40 km/s. This absorption is most dependent on the mass loss rate of the planet (with more material, there is more absorption) and on the ionizing flux (more ionization, less absorption, in an approximately linear relationship). By comparing their models to observation, the heat efficiency of HD 209458b (the planet modelled for simulations) can be predicted to be less than 50%. In addition, it can be seen that the observed Lyman-alpha absorption does not necessarily require charge exchange to accelerate the neutral hydrogen sufficiently.
    35 
    36 = Murray-Clay paper =
    37 Authors seek to numerically determine validity of hypothesis that hot Jupiters could be evaporated down to their rocky cores over the planetary lifetime. They use a one-dimensional model that includes heating/cooling terms, tidal gravity, and the effects of ionization on the mass-loss rate, and focus on the substellar point, at which tidal gravity and UV flux are greatest, thereby putting an upper limit on the possible mass-loss rate of the planet (by extending the one-dimensional result over the surface of the planet). They assume a planet of 1.4 R,,J,, and 0.7 M,,J,, and ignore the Coriolis force, under the assumption that the Lyman-alpha radiation from excited H is the only significant cooling term. Numerically, they use a relaxation solver, and find solutions iteratively by removing simplifying conditions one at a time.
    38 
    39 They find that, for main-sequence stars, about 20% of H is still neutral at the sonic point, and place an upper bound of ~3.3*10^10^ g/s on the mass loss rate. For hotter (T Tauri) stars, they find an upper bound of ~6.4*10^12^ g/s. At low flux, the mass loss is energy-limited, while for higher flux, the mass loss is radiation/recombination-limited. The assumption of a hydrodynamic wind is shown to be self-consistent, and they estimate that, due to the directionality of the tidal gravity and the UV irradiation, the maximum rate is an overestimate by ~4x. By reducing the wind speed to subsonic values and including a stellar wind, the day-side wind may be reduced or completely suppressed - they hypothesize that this may lead to night-side outflows, due to circulation of hot gases from the day-side.
    40 
    41 They compare observations to estimates from their model, and note that the disagreement in Lyman-alpha lines could be due to a variety of factors, including some missing physics or a cause unrelated to absorption by the planetary wind. A promising candidate is cited as acceleration of neutral hydrogen due to charge exchange. They note that modelled spectrally-unresolved measurements appear to be in agreement with observation.
    42 
    43 = Stone-Proga paper =
    44 Paper is a comparison of 2D simulations (of close-in hot Jupiters) to spherically symmetric simulations run by others. They characterize escape from the planet with the hydrodynamic escape parameter (lambda), ratio of gravitational potential to thermal potential - a small lambda suggests a thermal wind. They use no magnetism (hydrodynamic rather than MHD), and ignore heating and small-scale effects at the base of the wind. The wind is made self-consistent by fixing the density and energy at the base.
    45 
    46 In models with no stellar wind and an isothermal outflow, they find that the sonic surface of the wind is closer to the planet, with a slower radial velocity on the night side of the planet, and a very evident shock in the |v|/v,,s,, plot at various angles. For larger values of gamma, this shock produces a delta T, and the sonic surface is farther from the planet (but still nearer than in the spherically symmetric models). Introducing a stellar wind creates a back-swept profile, but has little effect on the sonic surfaces or mass-loss rate.