Melting probe technology for subsurface exploration of extraterrestrial ice – Critical refreezing length and the role of gravity K. Sch¨uller, J. Kowalski∗ AICES Graduate School, RWTH Aachen University, Schinkelstr. 2, 52062 Aachen, Germany. Abstract The ’Ocean Worlds’ of our Solar System are covered with ice, hence the water is not directly accessible. Using melting probe technology is one of the promising technological approaches to reach those scientifically interesting water reservoirs. Melting probes basically consist of a heated melting head on top of an elongated body that contains the scientific payload. The traditional engineering approach to design such melting probes starts from a global energy balance around the melting head and quantifies the power necessary to sustain a specific melting velocity while preventing the probe from refreezing and stall in the channel. Though this approach is sufficient to design simple melting probes for terrestrial applications, it is too simplistic to study the probe’s performance for environmental conditions found on some of the Ocean’s Worlds, e.g. a lower value of the gravitational acceleration. This will be important, however, when designing exploration technologies for extraterrestrial purposes. We tackle the problem by explicitly modeling the physical processes in the thin melt film between the probe and the underlying ice. Our model allows to study melting regimes on bodies of different gravitational acceleration, and we explicitly compare melting regimes on Europa, Enceladus and Mars. In addition to that, our model allows to quantify the heat losses due to convective transport around the melting probe. We discuss to which extent these heat losses can be utilized to avoid the necessity of a side wall heating system to prevent stall, and introduce the notion of the ’Critical Refreezing Length’. Our results allow to draw important conclusions towards the design of melting probe technology for future missions to icy bodies in our Solar System. Keywords: exploration technology, melting probe, contact melting, Europa, Enceladus Nomenclature W melting velocity (r, z) coordinates,seefigure1 A cross-sectional area (u, w) velocitycomponentsinthemeltfilm c heat capacity p Ste Stefan number F exerted force Re Reynoldsnumber F∗ buoyancy corrected exerted force α = k/(ρ c ) thermaldiffusivity g gravitational acceleration p Γ boundary h latent heat of melting m γ parameter introduced in equation (29) h∗ reduced latent heat of melting m δ melt film thickness J , Y Bessel functions of first and second kind 0 0 ǫ efficiency k thermal conductivity η power conversion efficiency L length of the melting probe µ dynamic viscosity of the water L∗ critical refreezing length ρ density m mass n, d fit constants Indices P Power arXiv:1803.04883v1 [physics.flu-dyn] 13 Mar 2018 p pressure L liquid phase (ρL, cp,L, kL, αL) p ambient pressure S solid phase (ρS , cp,S , kS , αS , TS ) atm Γ ˙ Q˙ heat flow rate C at the phase interface ( C, QC) Γ ˙ Q˙ minimum heat flow rate H at the melting head ( H, QH) min Γ ˙ q˙ heat flux E at the outflow boundary ( E, QE) R radius of the melting head T temperature 1. Introduction Tm melting temperature The presence of subglacial liquid water on the icy moons of ∗Corresponding author our Solar System [14] implies the possibility of habitable en- Email address: [email protected] (J. Kowalski) vironmental conditions. Especially the cryovolcanically active Saturnian moon Enceladus seems to be a promising candidate region bounded internally by a circular cylinder [6]. It reads [13] and there is some hope in the scientific community that an 2 L ∞ exp − αS b z exploration mission to Enceladus might unravel the existence ˙ 8kS TS W QL = dbdz, (2) of extraterrestrial life. Next generation mission concepts focus π Z0 Z0 2 2 b J0 (Rb)Y0 (Rb) on orbiting and sample-returning of plume material [13, 18]. Should these further strengthen any evidence for life, then the in which b is the integration argument of the Bessel functions natural next step is to sample and analyze the subglacial ocean J0 and Y0, and z is the spatial coordinate along the longintudinal directly [18, 10]. In order to access the extraterrestrial sub- axis of the melting probe. L and R are the radius and the length glacial water reservoirs, a thick ice layer must be penetrated. of the cylindrical melting probe and αS = kS /(ρS cp,S ) is the A very promising technological approach for this task is to thermal diffusivity of the ice, in which kS denotes its thermal use a thermal melting probe [10]. Melting probes enforce ice conductivity. penetration by heating, such that the ice in the vicinity of the Summing up the minimum heat flow rate to open the chan- probe melts and the probe eventually sinks down. Because nel Q˙ min and the lateral heat requirement Q˙ L provides a good the amount of necessary power roughly scales with the cross- approximation for the overall power necessary to sustain a spe- sectional area of the melting channel, a melting probe typically cific melting velocity while preventing the probe from refreez- looks like an elongated cylinder with a heated melting head. In ing and stall in the channel. This ’simple’ approach has been comparison to other ice penetration technologies, e.g. hot wa- used to design thermal melting probe robots both for terrestrial ter or mechanical ice drilling, the advantage of a melting probe field tests, e.g. in Antarctica and Greenland[11, 2], and for con- for space exploration purposes is its smaller, lighter, and me- ceptual studies to prepare extraterrestrial exploration missions chanically less complex design. Melting probes are not a novel [10, 24]. technology as they have already been applied successfully for The main technological issues for the latter have been sum- terrestrial research since the 1960’s [7, 16]. In recent years, marized in [22]. In that article, the authors find that a major however, more advanced melting probe designs have been pro- challenge for the design of melting probe technology for ex- posed and tested [24, 20, 11, 23, 9]. traterrestrial purposes, is the very low ice temperatures at the A very common and relevant engineering approach to design target bodies, which result in a very low efficiency of the melt- melting probes also dates back to the 1960’s and considers a ing process. Power efficiency is hence of major concern, espe- straight forward energy balance [1]: Knowing the power P im- cially, when facing restrictive power constraints during space ˙ ˙ plemented in a melting probe’s head as well as its conversion missions. Efficiency can be defined as ǫ = Qmin/(Qmin + ˙ ˙ efficiency η allows to infer on the heat flow rate at the melting i Qloss,i), in which i Qloss,i denotes the sum of all losses. One head’s surface Q˙ = ηP. The corresponding heat flux is given Ppotential loss is forP example given by convective losses within byq ˙ = Q˙/A, in which A stands for the cross-sectional area of the micro-scalemelt film between the melting probeand the ice. The efficiency associated with these losses can either be studied the melting head. The minimum heat flow rate Q˙ min required to operate the melting probe at a target melting velocity W is experimentally [8] or through advanced modeling techniques then given by the sum of the energy necessary to increase the [17] that go beyond the engineering design approach covered temperature of the ice in front of the probe and the energy that by (1). is eventually needed to melt the ice: Another great challenge is associated with initiating a melt- ing mission in a pressure regime that is below the triple pointof water (<6.1mbar), in which ice sublimates if heated. This com- ˙ − Qmin = WAρS hm + cp,S (Tm TS ) . (1) plicates the initial penetration phase of a melting probe, which h i is operated on bodies like Enceladus or Europa [8]. After the Here, ρS is the ice density, hm is the latent heat of melting, cp,S probe reaches a certain depth, however, the melting channel is is the heat capacity of the ice, TS is the ice temperature and Tm believed to refreeze and consequently the channel will sustain is the melting temperature of ice. Note, that from equation (1) a pressure above the triple point, such that the probe operates it is now evident that the melting velocity scales inversely with in a pure melting regime [21]. The initial low pressure regime the cross-sectional area of the melting probe. In orderto accom- can be further shortened by using a top cap as proposed in [5]. modate some scientific payload it is hence beneficial to increase Although a reliable and robust technological solution for the the probe’s length L rather than its radius R and often a probe initial low pressure phase is unarguably key to a successful mis- design is characterized by R << L. Such an elongated geom- sion, still the larger part of the melting transit through hundreds etry, however, poses another problem, namely the risk of stall of meters of ice will take place in a pressure regime above the due to refreezing [21]. Various concepts have been proposed triple point. Consequently, the vast amount of energy is spent to avoid stall during melting, e.g. by overheating the melting in a regime that is characterized by melting rather than sublima- head beyond the minimally required power Q˙ min given by (1) as tion. This motivates to further study efficiency and dynamics of proposed in [1], or by implementing a side wall heating system, melting probe technology in a pressure regime above the triple such as realized in [11].