A kinetic model of gas flow in a porous cometary mantle

A numerical study of gas flow through a porous cometary mantle is presented. A kinetic model based on the well-known Test Particle Monte Carlo Method for the solution of rarefied gas dynamics problems is proposed. The physical model consists of two spatial plane regions: the condensed ice phase and a porous dust mantle. The structure of the porous dust layer is described as a bundle of cylindrical inclined channels not crossing each other. A vertical temperature gradient may exist across the dust mantle. The aim is to investigate how the characteristics of molecular flow depend on the capillary length, inclination angle, and temperature gradient. Examples illustrating a significant deviation of some results from equilibrium values are shown. In particular, the gas velocity distribution at both ends of the pore is strongly non-Maxwellian if there is an important temperature contrast across the pore. The emergent gas flow rate is found to vary with the pore length/radius ratio in excellent agreement with Clausing's empirical formula. The degree of collimation of the flow is quantitatively studied as a function of the length/radius ratio, and consequences for the jet force of outgassing through a dust mantle or, indeed, a rough surface are estimated.     

Estimating the strength of the nucleus material of comet 67P Churyumov–Gerasimenko

Consideration is given to the estimates for the strength of the consolidated material forming the bulk of the nucleus of comet 67P Churyumov–Gerasimenko and those for the strength of the surface material overlying the consolidated material at the sites of the first and last contact of the Philae lander with the nucleus. The strength of the consolidated material was estimated by analyzing the terrain characteristics of the steep cliffs, where the material is exposed on the surface. Based on these estimates, the tensile strength of the material is in the range from 1.5 to 100 Pa; the shear strength, from ～13 to ?30 Pa; and the compressive strength, from 30 to 150 Pa, possibly up to 1.5 kPa. These are very low strength values. Given the dependence of the measurement results on the size of the measured object, they are similar to those of fresh dry snow at –10°C. The (compressive) strength of the surface material at the site of the first touchdown of Philae on the nucleus is estimated from the measurements of the dynamics of the surface impact by the spacecraft’s legs and the geometry of the impact pits as 1–3 kPa. For comparison with the measurement results for ice-containing materials in terrestrial laboratories, it needs to be taken into account that the rate of deformation by Philae’s legs is four orders of magnitude higher than that in typical terrestrial measurements, leading to a possible overestimation of the strength by roughly an order of magnitude. There was an attemp to put one of the MUPUS sensors into the surface material at the site of the last contact of Philae with the nucleus. Noticeable penetration of the tester probe was not achieved that led to estimation of the minimum compressive strength of the material to be ?4 MPa4 This fairly high strength appears to indicate the presence of highly porous ice with grains “frozen” at contacts.     

Observations of the long-lasting activity of the distant Comets 29P Schwassmann-Wachmann 1, C/2003 WT42 (LINEAR) and C/2002 VQ94 (LINEAR)

We investigated three comets, which are active at large heliocentric distances, using observations obtained at the 6-m BTA telescope (SAO RAS, Russia) in the photometric mode of the focal reducer SCORPIO. The three comets, 29P/Schwassmann-Wachmann 1, C/2003 WT42 (LINEAR), and C/2002 VQ94 (LINEAR), were observed after their perihelion passages at heliocentric distances between 5.5 and 7.08 AU. The dust production rates in terms of Afρ was measured for these comets. Using the retrieved values, an average dust production rate was derived under different model assumptions. A tentative calculation of the total mass loss of the comet nucleus within a certain observation period was executed. We calculated the corresponding thickness of the depleted uppermost layer where high-volatile ices completely sublimated. The results obtained in our study strongly support the idea that the observed activity of Comet SW1 requires a permanent demolition of the upper surface layers.     

Mass Transfer in the Near-Surface Layer of a Cometary Nucleus: A Gas-Kinetic Approach

A self-consistent model of the kinetically nonequilibrium near-surface layer of a cometary nucleus is developed on the basis of the gas-kinetic approach. The weight method of direct statistical simulation is used to model numerically the two-dimensional gas outflow from an ice sample subjected to radiative heating. The effective coefficient of water ice sublimation is estimated. Mass transfer in a porous ice and mineral (scattering) nonisothermal medium is investigated by the method of test particles, and the effective gas release is evaluated taking into account the proper rotation of the cometary nucleus for various model parameters. In these calculations, allowance is made for the kinetic character of the flow and volume sublimation and condensation of the volatile constituents of the material of the cometary nucleus.     

Activity of comets: Gas transport in the near-surface porous layers of a cometary nucleus

The gas transport through non-volatile random porous media is investigated numerically. We extend our previous research of the transport of molecules inside the uppermost layer of a cometary surface ( [Skorov and Rickman, 1995] and [Skorov et al., 2001]). We assess the validity of the simplified capillary model and its assumptions to simulate the gas flux trough the porous dust mantle as it has been applied in cometary physics. A microphysical computational model for molecular transport in random porous media formed by packed spheres is presented. The main transport characteristics such as the mean free path distribution and the permeability are calculated for a wide range of model parameters and compared with those obtained by more idealized models. The focus in this comparison is on limitations inherent in the capillary model. Finally a practical way is suggested to adjust the algebraic Clausing formula taking into consideration the nonlinear dependence of permeability on layer porosity. The retrieved dependence allows us to accurately calculate the permeability of layers whose thickness and porosity vary in the range of values expected for the near-surface regions of a cometary nucleus.     

Direct statistical simulation of the near-surface layers of a cometary atmosphere. II: A nonspherical nucleus

The study presents the results of numerical simulations of mass-transfer processes in the near-surface layer of the cometary nucleus and in the inner part of the cometary atmosphere, which is formed under the action of solar radiation. The gas-kinetic model of the inner part of the cometary atmosphere surrounding a spherical nucleus (Skorov et al., 2004) is extended to the case of a nonspherical nucleus with axial symmetry. After high-resolution images of comets 19P/Borrelly and Wild 2 have been obtained by Deep Space 1 and Stardust spacecraft, such an extension seems to be vital and important. The nucleus and the inner part of the coma are closely related to each other because of the permanent exchange of energy and mass; therefore, they are modeled consistently. As in the first part of our study, the boundary conditions at the inner boundary of the simulation domain, which are necessary for gas-kinetic simulations, were determined from the self-consistent model of heat and mass transfer in a porous cometary nucleus that was developed earlier by the authors. The model took into account the volumetric character of the radiation absorption in a porous sublimating medium, the kinetic regime of the transport of sublimation products in the pores, and the backward gas fluxes from the coma due to intermolecular collisions. We considered different models of the nucleus structure that determined the effective gas production. Using the direct simulation Monte Carlo method, we computed the two-dimensional gas flow from a heterogeneous nonspherical cometary nucleus. The simulations were performed using the SMILE software. The parallel computer implementation of the software made it possible to calculate the spatial structure of the gas flow for the entire circumnucleus zone.     

Thermophysical Modelling of Comet P/Borrelly Effects of Volume Energy Absorption and Volume Sublimation

In this work, we continue revising the theoretical basis ofnumerical models describing the transport of matter andenergy inside a porous dust-ice mixture at low temperature. Amodel of a light-absorbing near-surface layer of a comet nucleus isinvestigated. Gas transport is considered simultaneously with thesolution of the general heat transfer equation. Thequasi-stationary temperature distribution and the H₂O massflux and sublimation rate are computed for a nucleus model ofcomet 19P/Borrelly at the Deep Space 1 (DS1) encounter. Theenergy is deposited in a layer of about 20 particle radii: Thiscorresponds to a solid-state greenhouse effect. The surfacetemperature of the layer-absorbing model as well as the gasproduction rate are significantly smaller than the ones in thesurface-absorbing model. An active fraction of 40–50% would berequired to explain the observed water production rate ofP/Borrelly with our layer-absorption model at the time of the DS1encounter.     

On the Light-Absorbing Surface Layer of Cometary Nuclei: I. Radiative Transfer

A small but increasing volume of observations of cometary nuclei has accumulated during the past two decades. This development is accelerating with upcoming space missions such as Stardust, Contour, and Rosetta. In response to the growing need for a theoretical understanding of optical properties of cometary nuclei, we have calculated synthetic reflectance spectra in the wavelength region 0.2-2.0 μm, photometric colors in the Johnson-Kron-Cousins UBVRI system, and visual geometric albedos for a large number of porous ice-dust mixtures with differing composition, regolith grain sizes, and grain morphologies, such as core-mantle grains, dense clusters of such grains, and large irregular particles with internal scatterers. The calculations are based on Mie theory, the discrete dipole approximation, Hapke theory, and a numerical solution to the equation of radiative transfer in particulate media. In addition, wavelength-integrated directional-hemispherical albedos and flux attenuation profiles in the regolith as functions of depth have been calculated in order to improve the energy budget and treatment of energy boundary conditions in thermal models of cometary nuclei.Our results are compared with spectra and colors of observed cometary nuclei. Our main conclusions are that only regolith consisting of relatively large core-mantle grains, or clusters of smaller core-mantle grains, is capable of reproducing the red colors seen in comets; that ice-dust mixtures actually can be darker than ice-free regolith in certain circumstances; and that solar radiation sometimes penetrates to a depth that is comparable to the region in which diurnal temperature variations occur.