Modeling a class of thermal ice probes for accessing the solar system’s ocean worlds

The subsurface oceans of icy moons in the outer solar system are prime candidates for harboring extra-terrestrial life. Missions to melt through the ice shells of these worlds utilizing an ice probe or “cryobot” to access the oceans have been proposed. To design such a cryobot for a successful mission, the relationship between a cryobot’s heat budget and descent speed within any given environment must be understood. Large uncertainties in the ice shell thickness and temperature profile require that the problem be treated via probabilistic techniques, such as a Monte Carlo simulation. By contrast, the speed of a cryobot descending through ice with specified far-field temperature and material properties can, in principle, be quantified with little uncertainty. A model of cryobot performance that can be computed quickly and accurately will save significant computational time for the broader simulation, and allow for the mission’s basic thermal needs to be understood. The work reported here builds upon an existing analytic thermal model of ice probes that, while suitable for terrestrial conditions, begins to incur potentially problematic errors in extraterrestrial environments. The new model provides improvement to the previous model for flight-relevant probes in extra-terrestrial ice conditions, predicting the total heat budget to within 1.2% of high-fidelity numerical simulations while requiring a small fraction of the computational time. The model can be non-dimensionalized in a compact form, making it suitable for modeling a broad class of cylindrical ice probes in a probabilistic modeling framework for determining the time required to access the sub-surface ocean. Such models can then provide the high-level requirements for cryobot system architectures that will enable successful ocean-access missions.