3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.

Modeling of Heat Transfer in Microchannel Gas Flow

Due to the existence of a velocity slip and temperature jump on the solid walls, the heat transfer in microchannels significantly differs from the one in the macroscale. In our research, we have focused on the pressure driven gas flows in a simple finite microchannel geometry, with an entrance and an outlet, for low Reynolds (Re<200) and low Knudsen (Kn<0.01) numbers. For such a regime, the slip induced phenomena are strongly connected with the viscous effects. As a result, heat transfer is also significantly altered. For the optimization of flow conditions, we have investigated various temperature gradient configurations, additionally changing Reynolds and Knudsen numbers. The entrance effects, slip flow, and temperature jump lead to complex relations between flow behavior and heat transfer. We have shown that slip effects are generally insignificant for flow behavior. However, two configuration setups (hot wall cold gas and cold wall hot gas) are affected by slip in distinguishably different ways. For the first one, which concerns turbomachinery, the mass flow rate can increase by about 1% in relation to the no-slip case, depending on the wall-gas temperature difference. Heat transfer is more significantly altered. The Nusselt number between slip and no-slip cases at the outlet of the microchannel is increased by about 10%.

Observations of the northern seasonal polar cap on Mars: I. Spring sublimation activity and processes

Spring sublimation of the seasonal CO2 northern polar cap is a dynamic process in the current Mars climate. Phenomena include dark fans of dune material propelled out onto the seasonal ice layer, polygonal cracks in the seasonal ice, sand flow down slipfaces, and outbreaks of gas and sand around the dune margins. These phenomena are concentrated on the north polar erg that encircles the northern residual polar cap. The Mars Reconnaissance Orbiter has been in orbit for three Mars years, allowing us to observe three northern spring seasons. Activity is consistent with and well described by the Kieffer model of basal sublimation of the seasonal layer of ice applied originally in the southern hemisphere. Three typical weak spots have been identified on the dunes for escape of gas sublimed from the bottom of the seasonal ice layer: the crest of the dune, the interface of the dune with the interdune substrate, and through polygonal cracks in the ice. Pressurized gas flows through these vents and carries out material entrained from the dune. Furrows in the dunes channel gas to outbreak points and may be the northern equivalent of southern radially-organized channels ("araneiform" terrain), albeit not permanent. Properties of the seasonal CO2 ice layer are derived from timing of seasonal events such as when final sublimation occurs. Modification of dune morphology shows that landscape evolution is occurring on Mars today, driven by seasonal activity associated with sublimation of the seasonal CO2 polar cap.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Vapor deposition coating of fused silica tubes with amorphous selenium

A set of optimized deposition conditions for the inner wall coating of fused silica tubes with amorphous selenium was elaborated. The method is based on the vapor transport deposition of pure elemental selenium on a cooled substrate held at liquid nitrogen temperatures. Morphological and structural examination of the deposited layer was performed by optical microscopy and X-ray diffraction studies. Neutron activated selenium was used to monitor the deposition pattern and its stability under high gas flows. Monte Carlo simulations allowed the estimation of the different Se species composing the amorphous phase, at the given experimental deposition conditions. The versatility of the coating method presented in this work allows for the coating of tubes of different lengths and diameters, opening the way for several applications of amorphous selenium films in various fields.