Rarefied gas dynamic simulation of transfer and escape in the Pluto–Charon system

Abstract

We apply the direct simulation Monte Carlo rarefied gas dynamic technique to simulations of Pluto's rarefied upper atmosphere motivated by the need to better understand New Horizons (NH) data. We present a novel three-dimensional DSMC model of the atmosphere that spans from several hundred km below the exobase – where continuum flow transitions to the rarefied regime – to fully free-molecular flow hundreds of thousands of km from Pluto's center. We find molecular collisions in Pluto's upper atmosphere to be significant in shaping the flowfield, both by promoting flux from the plutonian exobase to Charon and by increasing the proportion of that flux generated on the exobase's anti-Charon hemisphere. Our model accounts for the gravitational fields of both Pluto and Charon, the centripetal and Coriolis forces due to the rotation of Pluto in our reference frame, and the presence of Charon as a temporary sink for impacting particles. Using this model, we analyze the escape processes of N2 and CH4 from Pluto across different solar heating conditions, and evaluate the three-dimensional structure of the upper plutonian atmosphere, including gas transfer to and deposition on Charon.We find results consistent with the NH-determined escape rate, upper atmospheric temperature, and lack of a detectable Charon atmosphere. Gas-transfer structures are noted in a binary atmospheric configuration, including preferential deposition of material from Pluto's escaping atmosphere onto Charon's leading hemisphere that peaks at 315° E on the equator. As the moon gravitationally focuses incident flow, a high density structure forms in its wake. If molecules are permitted to escape from Charon in diffuse reflections from its surface, a returning flux forms to Pluto's exobase, preferentially directed toward its trailing hemisphere. Charon is capable of supporting a thin atmosphere at column densities as high as 1.5 ×1017m−2 in simulations with a plutonian exobase condition similar to the NH encounter. Results computed from a fit to the NH encounter exobase (Gladstone et al., 2016) predict a system escape rate of 7 ×1025 CH4 s−1 in close agreement with those reported by NH (Bagenal et al., 2016; Gladstone et al., 2016), and a net depositional flux to Charon of 2 ×1024s−1, of which ∼98% is methane.