Arxiv:1803.04883V1 [Physics.Flu-Dyn] 13 Mar 2018

arXiv:1803.04883v1 [physics.flu-dyn] 13 Mar 2018 oad h eino etn rb ehooyfrftr mi future Keywords: for technology probe melting of design the towards h oe eesr ossanaseicmligvlct wh velocity su melting is specific approach a this sustain f Though starts to probes necessary an melting of power top such the on design head to melting approach heated approache engineering a technological of promising consist basically the probes of one is technology ceeain hswl eiprat oee,we design when however, important, be on will This found acceleration. conditions environmental for performance probe’s ytmt rvn tl,aditouetento fte’Cri the of notion the introduce and stall, prevent to system Nomenclature Abstract rudtemligpoe edsust hc xetteehe these extent which to discuss We o that, probe. to melting addition the In di around of Mars. bodies and on Enceladus regimes Europa, melting on study regimes to allows model Our ice. c h g h k n m T T R q A F F L L P Q Q J p p ˙ ˙ ˙ h OenWrd’o u oa ytmaecvrdwt c,h ice, with covered are System Solar our of Worlds’ ’Ocean The etcl h rbe yepiil oeigtepyia pr physical the modeling explicitly by problem the tackle We p 0 ∗ m m ∗ , ∗ atm m ∗ min mi address: Email author Corresponding , d etn rb ehooyfrsbufc xlrto fext of exploration subsurface for technology probe Melting Y 0 rs-etoa area cross-sectional xre force exerted capacity heat rvttoa acceleration gravitational force exerted corrected buoyancy aetha fmelting of heat latent eue aetha fmelting of heat latent reduced egho h etn probe melting the of length kind conductivity second thermal and first of functions Bessel mass length refreezing critical pressure Power constants fit etn temperature melting temperature head melting the of radius flux heat iiu etflwrate flow heat minimum rate flow heat pressure ambient xlrto ehooy etn rb,cnatmlig E melting, contact probe, melting technology, exploration [email protected] IE rdaeSho,RT ahnUiest,Schinkelst University, Aachen RWTH School, Graduate AICES ffi in odsg ipemligpoe o ersra appli terrestrial for probes melting simple design to cient ereiglnt n h oeo gravity of role the and length refreezing J Kowalski) (J. .Shulr .Kowalski Schüller, J. K. ia ereigLnt’ u eut lo oda importa draw to allow results Our Length’. Refreezing tical n xlrto ehooisfretaersra purpos extraterrestrial for technologies exploration ing o lbleeg aac rudtemligha n quant and head melting the around balance energy global a rom sost c oisi u oa System. Solar our in bodies icy to ssions rmdlalw oqatf h etlse u oconvective to due losses heat the quantify to allows model ur oeo h ca’ ols ..alwrvleo h gravitat the of value lower a e.g. Worlds, Ocean’s the of some orahtoesiniclyitrsigwtrreservoir water interesting scientifically those reach to s ff tlse a euiie oaodtencsiyo iewall side a of necessity the avoid to utilized be can losses at lnae oyta otistesinicpyod h tra The payload. scientific the contains that body elongated l rvnigtepoefo ereigadsali h cha the in stall and refreezing from probe the preventing ile cse nteti etfimbtentepoeadteunderly the and probe the between film melt thin the in ocesses rn rvttoa ceeain n eepiil compa explicitly we and acceleration, gravitational erent .Introduction 1. u oa ytm[4 mle h osblt fhbtbee habitable cryovolcanically of the possibility Especially the conditions. implies vironmental [14] System Solar our Indices netewtri o ietyacsil.Uigmligpro melting Using accessible. directly not is water the ence C S ρ µ ǫ δ γ Γ η ( eRyod number Reynolds number Stefan α Re Ste ( W H E L u r rp,Enceladus uropa, h rsneo ugaillqi ae nteiymosof moons icy the on water liquid subglacial of presence The , , = z w oriae,sefiue1 figure see coordinates, ) k eoiycmoet nteml film melt the in components velocity ) ttepaeitrae( interface phase the at ttemligha ( head melting the at tteoto onay( boundary outflow the at oi hs ( phase solid iudpae( phase liquid / ( ρ ∗ .2 26 ahn Germany. Aachen, 52062 2, r. c p hra di thermal ) density water the of viscosity dynamic e thickness film melt (29) equation in introduced parameter boundary oe ovrine conversion power etn velocity melting ffi ρ ciency ρ S L , ain,i stosmlsi osuythe study to simplistic too is it cations, aersra c Critical – ice raterrestrial , c c p p , S , L , Γ ff , k H usivity k Γ S L , , C , Q α Γ , ˙ α H S Q E L ˙ , , ) ) C T Q ffi ) ˙ S E ciency ) ) es. tconclusions nt .Melting s. emelting re transport ditional heating active ional nnel. ifies ing be n- 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 potentialP 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 melting 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]. point, while neglecting low pressure effects for the time being. The amount of heat required to avoid stall, in the following One aspect, which has not been investigated so far even for referred to as the lateral heat requirement Q˙ L, has been deter- conditions above the triple point, is the effect of the gravity on mined in [1] based on quantifying heat conduction in an infinite the melting process. The value of the gravitational accelera- 2 tion is for example much smaller on Enceladus than it is on Europa, or even on Earth. A melting probe of a certain mass will hence exert much less force on the ice on Enceladus than it does on Europa, or on Earth. In this article, we will investigate this question by developing a physics-based model that explicitly accounts for the processes in the thin melt film between the probe and the underlying ice. This kind of theory is commonly referred to as close-contact melting theory [3]. In section 2, we will describe the physical situation for a cylindrical melting probe that heats at constant power. In section 3, we derive a model, which quantifies the conductive losses in the melt film as a function of the exerted contact force, which in itself depends on the gravitational acceleration. In contrast to the majority of other models that have been proposed for the design of melting probes, the additional consideration of exerted forces is the particular strength of our approach, which Melting Probe makes an investigation of the role of gravitational acceleration on the melting process possible for the first time. Section 4 is devoted to results, while discussing concrete melting scenarios. Water Specifically, we will apply the derived model to analyze melting regimes on Enceladus, Europa and Mars. A direct quantification of the conductive losses allows us to discuss the necessity Ice of wall heaters to prevent stall under different extraterrestrial environmental conditions. Finally, we are able to approximate Figure 1: Schematic of the physical problem. The figure shows a cylindrical the temperature at the surface of the melting head. We will con- melting probe, which melts downwards. Note that the melt film thickness is not clude with a discussion of our findings, and an outlook. to scale.

2. Physical model along the probe’s longitudinal axis. A constant force F is acting on the melting probe. In the absence of additional actuator sys- We consider a cylindrical melting probe of length L and ra- tems, it would simplify to the weight of the probe that acts in z dius R that melts through ideal ice at a constant melting velocity direction. The exerted force causes a squeezing of the melt film, W, hence in an equilibrium melting regime. With ideal ice we and consequently, induces a flow field within the melt film with refer to polycrystalline ice without pores or fractures. Although velocity components u and w in r- and z-direction, respectively. these artifacts may be present in reality, especially close to the surface of icy moons [12], considering pure ice yields a con- servative approximation for the necessary power of a melting 3. Mathematical model probe. The quasi-steady melting process is sketched in figure 1. A In this work, we extend straight forward energy balance ar- δ thin melt film of uniform thickness separates the melting head guments to the entire melt film, i.e. to the domain bounded by surface Γ , from the ice domain. The water ice interface is H Γ , Γ and Γ . In order to describe the flow within the melt Γ Γ H C E denoted by C. At H , heat is transferred from the melting head film, we use an extended lubrication theory [15, 4]. The lubri- Q˙ Γ into the melt film according to H. At the phase interface C, cation approximation is a valid simplification for the physical heat is transferred from the melt film into the ice according to problem, because the melt film is very thin compared to the Q˙ Γ C. At the lateral outflow boundary of the melt film E, we width of the melting head, i.e. δ/R << 1. Scaling analysis then observe convectivetransport of heat out of the melt film domain shows that inertia terms are negligible compared to the pres- Q˙ Q˙ according to E. Therefore, E denotes the convective losses sure gradient, and also ∂/∂r2 << ∂/∂z2 as long as δ/R << 1 in the melt film, namely heat that has been transferred into the and (δ/R)Re << 1, in which Re is the Reynolds number. melt film but that has not been utilized to increase the probe’s With these assumptions, conservation of mass, momentum melting velocity. Considering the melt film as a closed system, and energy for the melt film are given by [15] then yields the following energy balance 1 ∂(ru) ∂w Q˙ − Q˙ − Q˙ = 0 (3) + = 0 (4) H E C r ∂r ∂z The coordinate system with axes r and z is fixed with respect ∂2u dp µ = (5) to the melting probe. Its origin is located at the center of the ∂z2 dr melting head surface ΓH. Gravitational acceleration g acts in ∂T ∂T ∂2T u + w = α (6) positive z-direction. It induces a hydrostatic pressure gradient ∂r ∂z L ∂z2 3 in which p is the pressure, µ is the dynamic viscosity and αL = 3.2. Heat flow rate at the outflow boundary Q˙ E kL/(ρL cp,L) is the thermal diffusivity of water. To avoid unnecessary complications, we will assume that all At the melting head surface ΓH we have no-slip/zero-inflow thermo-physical properties are phase-wise constant. Integrating conditions and a prescribed heat flow rate Q˙ H, which yields the momentum balance (5) twice with respect to z and substituting in the boundary conditions (7) and (9) yields the horizontal u(r, 0) = w(r, 0) = 0 (7) melt film velocity ∂T Q˙ = − H (8) 1 dp 2 u = z(z − δ). (15) ∂z z=0 πR kL 2µ dr

At the phase interface ΓC, we again have no-slip in horizontal Substituting the velocity (15) into continuity equation (4), inte- direction, however this time, the inflow velocity is equal to the grating the resulting equation over the melt film thickness and melting velocity of the probe W scaled by a factor ρS /ρL that inserting the boundary conditions (7) and (9) yields accounts for the density change across the phase interface: ρ r r dp − S W rd˜ r˜ = δ3. (16) ρS ρL Z 12µ dr u(r,δ) = 0, w(r,δ) = − W (9) 0 ρL The pressure gradient is found by integrating equation (16) and making use of dp/dr = 0 at r = 0 due to axisymmetry: Furthermore, two temperature related boundary conditions hold ρS at the phase interface: dp 6µ ρ Wr = − L (17) dr δ3 T(r,δ) = Tm (10) ∂T ρ W The horizontal velocity component (15) is then evaluated to be = − S h∗ (11) m ρS − ∂z z=δ kL 3 ρ Wrz(z δ) u = − L . (18) δ3 Here, Tm is the melting temperature and equation (11) is the so-called Stefan condition that accounts for local energy bal- To derive an equation for the melt film thickness, we at first ance at the phase interface. It contains the reduced latent heat use continuity to write the temperature equation in conservation ∗ − of melting hm = hm + cp,S (Tm TS ), which is the sum of the form. Integration over the melt film thickness then yields latent heat h and the specific energy necessary to increase the m δ 1 ∂ (r u T) ∂ (w T) δ ∂2T temperature of the ice. + dz = αL dz (19) " r ∂r ∂z # ∂z2 For a vertically oriented melting probe in a water filled melt- Z0 Z0 ing channel, the pressure increases linearly along the probe’s This can be simplified by substituting in boundary conditions length, which causes the pressure at the outflow boundary (r = (7)–(11): R) tobe δ d − ρS − ρS W ∗ ∂T ruTdz WTm = αL hm + (20) p(R) = p + ρ gL (12) rdr Z0 ρL kL ∂z z=0! 0 L We make a quadratic polynomial Ansatz for the temperature in which p0 is the pressure at the back of the probe. field in z-direction, which satisfies the three boundary conditions (8), (10) and (11). The resulting polynomial is − ∗ 3.1. Heat flow rate at the phase interface Q˙C q˙H ρS Whm 2 2 q˙H T = z − δ + (δ − z) + Tm (21) 2δkL kL At the phase interface Γ , the heat flow rate equals the heat C in which we recall that flux multiplied by the probe’s cross-sectional area 2 q˙H = Q˙ H/(πR ) (22) ∂ ˙ − 2 T QC = πR kL (13) is the heat flux at the melting head. Substituting the temperature ∂z z=δ Ansatz (21) and the horizontal velocity (18) into the conserva-

Together with the Stefan condition (11) this yields tion of energy (20), finally yields an expressionfor the melt film thickness δ Q˙ = πR2ρ Wh∗ (14) − ∗ C S m 20αLρL q˙H ρS Whm δ = ∗ (23) WρS 3˙qH + 7ρS Wh Hence, equation (14) is formally equivalent to the straight- m forward energy balance (1) discussed in the introduction. The which can also be written in terms of the heat flow rates at the important innovation in this paper is that we consider the realis- melting head surface and at the phase interface: tic case, in which less power is utilized for forward motion than 20α ρ Q˙ − Q˙ is inputted into the system, hence Q˙ < Q˙ . The next section L L H C C H δ = (24) will be devoted to quantifying the heat losses Q˙ E = Q˙ H − Q˙C. WρS 3Q˙ H + 7Q˙C 4 Still δ cannot be computed as the melting velocity is not known. Table 1: Thermo-physical properties. Data from [10] and [22]. In order to close the system, we consider the force balance as Mars(polar) Enceladus Europa an additional relation. It states that the integral of the pressure g [m/s2] 3.7 0.1 1.3 in the melt film must balance the exerted force F (e.g. given by kL [W/(mK)] 0.6 0.6 0.6 the gravity proportional weight of the melting probe). This is 3 ρL [kg/m ] 1000 1000 1000 R 2π cp,L [J/(kgK)] 4200 4200 4200 F = (p − p0)rdφdr (25) µ [N s/m2] 0.0013 0.0013 0.0013 Z0 Z0 kS [W/(mK)] 2.5 3.0 3.5 3 An expression for the pressure is given by integrating the pres- ρS [kg/m ] 921.3 925.1 927.8 sure gradient (17) subject to the pressure boundary condition cp,S [J/(kgK)] 1877 1656 1476 (12), which yields hm [J/kg] 333700 333700 333700 Tm [K] 273 273 273 3µ ρS W(R2 − r2) ρL T [K] 210 150 100 p − p = + ρ gL. (26) S 0 δ3 L

Substitution into the force balance (25) again results in an equa- 4. Results and discussion tion for the melt film thickness δ In this section, we apply our model to analyze the melting πR4µ ρS W 1/3 3 ρL process of a cylindrical melting probe. We choose a length of δ = (27) 2 F∗  L = 1m, a radius of R = 0.06m and a mass of m = 25kg.     These dimensions correspond to both the IceMole flight model  ∗ 2 in which F = F − πR ρLgL is the buoyancy corrected force. [10] and the Cryobot [24], both of which have been designed Note, that equation (27) is a mechanical approximation to the for extraterrestrial application. We consider three possible ex- melt film thickness δ, while the previous equation (24) has been traterrestrial targets, namely the polar caps on Mars, the south a thermodynamic approximation. Equating both yields a rela- polar region on Enceladus and Europa. The relevant thermo- tion that allows to quantify the fraction of the input heat that is physical properties are summarized in table 1. Within the scope utilized for the forward motion of the probe of this study and as discussed in the introduction, we will focus on pure melting without sublimation. 1 − 3γ Q˙C = Q˙ H , (28) 7γ + 1 4.1. The role of gravitational acceleration in which In the absence of an actuator system, the buoyancy corrected ∗ 4/3 1/3 force of the melting probe F is given by its weight and a buoy- 1 ρS 3πµ ancy term. Both are proportional to the gravitational accelera- γ = WR ∗ . (29) 20αL ρL ! 2F ! tion by definition. Studying the role of gravitational acceleration on the melting process, hence translates into studying its 3.3. Model summary force dependency. The effect of F∗ on the melting process can Equation (28) finally allows us to determine the convective directly be observed by analyzing the melt film thickness δ (27) ˙ losses at the melt film outflow boundary Q˙ E and the heat flow rate at the outflow boundary QE (30). For heavy melting probes or environments that exhibit a large grav- 1 − 3γ itational acceleration, the force F∗ increases, and the melt film Q˙ E = Q˙ H − Q˙C = 1 − Q˙ H (30) 7γ + 1! thickness tends to zero. Consequently, convective losses Q˙ E get small, which means that melting gets efficient and Q˙ H ≈ Q˙C. Our results allow us to split the heat flow rate at the melting This corresponds to the straight-forward energy balance (1), head surface into the portion of heat that is used for the for- which is assumed to operate at optimal efficiency, and hence is ˙ ward motion of the melting probe QC, and additional convective a special case of our model. If on the other hand exerted forces ˙ losses given as a fraction of QC get small, either for very light melting probes or environments + that exhibit a small gravitational acceleration, then the convec- 7γ 1 ∗ Q˙ H = Q˙C + − 1 Q˙C (31) tive losses increase. The limiting case F → 0 deserves fur- 1 − 3γ ! ther attention, as it results in the intuitive solution of a stagnant This relation can be used in two ways: One can either determine melting probe (W → 0). While the melting head remains ther- the required heat flow rate at the melting head Q˙ H to sustain a mally active, the phase interface moves gradually away from the target melting velocity for a given exerted force. The factor melting probe without any motion of probe itself. This physical (1 − 3γ)/(7γ + 1)is then to be interpreted as an efficiency factor situation is commonly referred to as a Stefan problem. less than 100%, which has to be compensated. Alternatively, Figure 2 and 3 show how the melting velocity varies with the one can determine the effective melting velocity, given a certain buoyancy corrected force for environmental conditions found input power Q˙ H. Since both γ and Q˙C depend on the W, the on Mars, Europa and Enceladus (see table 1). Figure 2 assumes latter requires a numerial root finding strategy. a heat flow rate at the melting head of Q˙ H = 5kW, while figure 5 3 assumes Q˙ H = 1kW. The vertical lines mark the buoyancy corrected forces for a melting probe of 25kg, which are 50.65N on Mars, 1.37N on Enceladus and 17.8N on Europa. Equa- tion (31) was solved numerically in order to obtain the melting velocities in the figures. 3.5 The plots clearly indicate that the melting velocity monotoni- cally increases with increasing force. In the limit (F∗ →∞) the melting velocity converges asymptotically to a value that cor- 3 responds to the velocity obtained through the straight-forward energy balance (1). The latter observation is more pronounced 2.5 in the Q˙ H = 1kW case, as the influence of reasonable force variations (1N – 1000N) on the melting velocity is larger if the heat flow rate at the melting head is large. Increasing the force 2 spectrum considered in the plot for the Q˙ H = 5kW case would however result in a similar asymptotic behavior. An additional observation is also that the melting velocity increases with in- 1.5 creasing ice temperature. The highest melting velocity for each force value is consistently found for Mars environmental con- 1 -2 -1 0 1 2 3 ditions that are the warmest, the lowest for Europa conditions, 10 10 10 10 10 10 which is also the coldest of the three. An interesting finding of our model is the combined impact of ice temperature and grav- Figure 2: Results for a melting probe (L = 1m, R = 0.06 m) with a heat flow itational acceleration. Despite of melting through much colder rate at the melting head of Q˙ H = 5 kW. The plot shows the melting velocity over ∗ − 2 ice on Europa, our 25kg probe still operates at a higher veloc- the buoyancy corrected force F = F πR ρLgL on different extraterrestrial ity than on Enceladus, which is warmer but has a significantly ice environments. The vertical lines mark the buoyancy corrected force at the target. lower gravitational acceleration. Hence, convective losses due to a low value of the exerted contact force seem to dominate the effect of the warmer ice on Enceladus. Our results refute earlier findings reported in [22] and based on [19] that the force of gravity is necessary to set the direction of penetration, but does not (to first order) control the melting velocity. Increas- ing the melting velocity of a probe on Enceladus hence cannot be achieved by increasing its mass alone, as the high masses needed would not be feasible for a space mission that starts on Earth. One solution would be to introduce an additional force 0.8 by means of an actuator system, e.g. an ice screw at the melting head, as it is realized with the IceMole melting probe [11]. 0.75

0.7 4.2. Critical refreezing length 0.65 One key question during the design of melting probes is whether a side wall heating system must be implemented in or- 0.6 der to prevent refreezing at the probe’s hull. A wrong design could cause stall of the melting probe and poses a severe risk 0.55 to a melting probe mission. To mitigate this risk, a common 0.5 approach is to over design the melting probe by utilizing wall heaters, even if they would not be needed. 0.45

In this section, we will investigate to which extent the con- 0.4 vective melt film losses naturally contribute to prevent refreez- 10-2 10-1 100 101 102 103 ing. Our analysis is based on an approximation to the lateral heat requirement (2) derived by [1, 22] Figure 3: Results for a melting probe (L = 1m, R = 0.06 m) with a heat flow ˙ ˙ ∗ rate at the melting head of QH = 1 kW. The plot shows the melting velocity over QL d L ∗ − 2 = n σ σ = the buoyancy corrected force F = F πR ρLgL on different extraterrestrial 2 with 2 (32) R W(Tm − TS ) WR ice environments. The vertical lines mark the buoyancy corrected force at the target. with the fit constants n = 932 Ws/K/m3 and d = 0.726. Ac- cording to [22], equation (32) is valid for 5 × 104 s/m2 <σ< 108 s/m2. Note the difference between the actual length of the probe L and L∗ used in equation (32). L∗ stands for the distance 6 2.5 4 Critical refreezing length Melting velocity with convective losses Melting velocity without convective losses 3.5 2

3 1 3.5 1.5 2.5 Critical refreezing length 0.9 Melting velocity with convective losses Melting velocity without convective losses 3 2 0.8 1 0.7 2.5 1.5 0.6 0.5 1 0.5 2

0.4 0 0.5 1.5 1000 1500 2000 2500 3000 3500 4000 4500 5000 0.3

0.2 1 Figure 4: Critical refreezing length and melting velocity over the heat flow rate 0.1 at the melting head for different buoyancy corrected forces in 210 K cold ice on 0 0.5 Mars. 1000 1500 2000 2500 3000 3500 4000 4500 5000 along the melting channel after which the temperature has re- Figure 5: Critical refreezing length and melting velocity over the heat flow rate laxed to melting temperature, hence the distance from the melt- at the melting head for different buoyancy corrected forces in 150 K cold ice on ing head after which refreezing would occur. It will be referred Enceladus. to as the critical refreezing length in the ongoing of this article. Assuming that the convective losses of the melt film Q˙ E contribute towards keeping the lateral melt channel open, we can investigate L∗ for

Q˙ L = Q˙ E. (33) Substituting equations (30) and (32) into equation (33) yields

Q˙ Wd−1R2(d−1) 7γ + 1 1/d L∗ = C − 1 . (34) " n (Tm − TS ) 1 − 3γ !# Hence, wall heaters to prevent refreezing must be implemented 0.6 3 ∗ Critical refreezing length if the probe is shorter than L . Melting velocity with convective losses Figures 4, 5 and 6 show the critical refreezing length for 0.5 Melting velocity without convective losses 2.5 the melting probe using equation (34) in combination with the melting velocity for Mars polar, Enceladus and Europa environ- 0.4 2 mental conditions, respectively. Note, that for a constant heat ∗ flow rate at the melting head, the critical refreezing length L 0.3 1.5 increases if the exerted force decreases. This can be explained by looking at the melting velocity, which depends on the ex- 0.2 1 erted force. When the melting probe slows down, the melt film thickness increases and higher convective losses occur. For one 0.1 0.5 value of the buoyancy corrected force, the critical refreezing length increases with increasing heat flow rate at the melting 0 0 head. Again, this is a result of increasing convective losses. 1000 1500 2000 2500 3000 3500 4000 4500 5000 When comparing the results for 210K (figure 4), 150K (figure 5) and 100K (figure 6), it can be seen that the critical refreezing length is smaller in a cold environment. This is be- Figure 6: Critical refreezing length and melting velocity over the heat flow rate ˙ ∝ − at the melting head for different buoyancy corrected forces in 100 K cold ice on cause more heat (QL (Tm TS ), see equation(32)) is needed at Europa. the probe’s hull to prevent stall. Our melting probe has a length of 1m. Since the critical refreezing length is always smaller than 1 m for reasonable buoyancy corrected forces (F∗ > 10 N), the melting probe will stall if no additional wall heaters are provided. 7 4.3. Optimal melting probe length 101 In the previous subsection, the critical refreezing length L∗ has been introduced as a measure to check if a melting probe of given length, radius and mass needs wall heaters, and at which 0 minimal distance from the head these must be implemented. If 10 the critical refreezing length is smaller than the length of the melting probe, the natural next question is: What probe length -1 would make the use of a side wall heating system dispensable. 10 Note, that this probe length is not given by the aforementioned critical refreeze length, as length of the probe itself alters the ∗ model. Instead, the critical refreezing length L in equation 10-2 (34) is at replaced by the length of the probe L. The resulting equation is

10-3 Q˙ Wd−1R2(d−1) 7γ + 1 1/d 1000 1500 2000 2500 3000 3500 4000 4500 5000 0 = L − C − 1 (35) " n (Tm − TS ) 1 − 3γ !# Note that γ and the melting velocity are now functions of the Figure 7: Maximum melting probe length for a 25 kg and R = 0.06 m melting ∗ probe that operates without a side wall heating system in different ice environ- buoyancy corrected force F and hence depend on the melting ments. The dashed line is used to mark the 1 m level in the plot. probe length L. Therefore, equation (35) has to be solved numerically for L. If the mass of a melting probe is given and no actuators are Therefore, more heat is required at the walls to keep the water used to increase the exerted force of the melting probe, the defi- from refreezing. Interestingly both effects cancel out for Mars nition of the buoyancycorrectedforce can be exploited to derive and Enceladus conditions. The maximum melting probe length an upper limit for the length of the melting probe as a function for Europa is in the order of 0.1m for Q˙ H = 5000W. A probe of its mass and radius. Motion is only possible if the buoyancy of this dimension and a mass of 25kg including scientific in- corrected force is larger than zero, which yields the relation struments is technically not feasible. We therefore conclude that melting probes for future Europa missions must be utilized m L < (36) with wall heaters, while it might be possible to design a system πR2ρ L for Enceladus and Mars without additional side wall heaters. For the melting probe of 25kg and R = 0.06m, equation (36) yields an upper limit for the length of approximately 2.21m, if 4.4. Melting head temperature and Stefan number 3 ρL = 1000kg/m . Figure 7 shows the maximum length of a melting probe (m = The temperature at the melting head as predicted by our 25kg, R = 0.06m) that can operatewithout side wall heaters for model is found by evaluating the quadratic temperature pro- different ice environments. As expected, the probe’s length can file (21) at z = 0, while making use of the computed melt film be significantly larger for warm ice. Therefore, most melting thickness (24). This yields probesfortemperate ice on Earth (i.e. for ice with a temperature ˙2 − ˙2 close of 273.15K) do not need wall heaters. However, there 10 QH QC are also colder ices on Earth with temperatures in the order of TH = + Tm (37) c ρ πR2W 3Q˙ + 7Q˙ 256K, e.g. in Antarctica. If the heat flow rate at the melting p,L S H C head is smaller than approximately 3250W, the melting probe Figure 8 shows the melting head temperature for three force needs to be smaller than 1m. Otherwise a wall heating system scenarios as expected on Enceladus conditions. Knowing the must be utilized to prevent stall of the probe. If the heat flow surface temperature T also allows us to quantify the Stefan rate is larger than 3250W, the resulting heat flow rate at the H number regime for our melting process. The Stefan number is outflow boundary Q˙ is sufficient to prevent stall of the probe E typically defined as Ste = c (T − T )/h∗ and relates the sen- in 256K cold ice. p H m m sible heat to the reduced latent heat of the considered system. Even for the high power scenario Q˙ =5kW, we observe that H Since the task of a melting probe is to enforce phase change the maximum melting probe length for all extraterrestrial en- of the surrounding ice, sensible heat in the system is undesir- vironments is less than 0.5m. Interestingly, the melting probe able and can hence be interpreted as a measure of inefficiency. length for Enceladus is very close to the one for Mars. This Therefore, small Stefan numbers indicate efficient melting. is again due to the combined effect of the low gravitational acceleration and the difference in ice temperature. A small grav- Based on the evaluated surface temperature of the melting head (37), the Stefan number can be translated in itational acceleration, as it is present on Enceladus compared to Mars, will decrease the buoyancy corrected force and hence 10 Q˙ 2 − Q˙ 2 increase the convective losses Q˙ . On the other hand, the ice H C E Ste = (38) temperature on Enceladus is significantly less than on Mars. Q˙C 3Q˙ H + 7Q˙C 8 of the convective losses in the melt film, and hence provides a 3 10 10 more accurate description of the efficiency of the melting process than a so far often used straight-forward energy balance. Furthermore, we introduced the notion of the critical refreezing length and used it to discuss to which extent convective losses 10 2 1 during the melting process can naturally be utilized to avoid otherwise over designed side wall heating systems. The con- cept of the critical refreezing length is not only an important measure for melting probes, but also important, when mitigat- 10 1 0.1 ing the risk due to radioisotope power sources in the context of planetary protection on the icy bodies in our Solar System as described in [12]. The three major findings of our work are:

10 0 0.01 1000 1500 2000 2500 3000 3500 4000 4500 5000 • Other than previous literature has stated, we find that the exerted force, hence the gravitational acceleration of the Figure 8: Melting head temperatures and Stefan numbers as a function of heat environment can have a significant effect on the melting flow rate at the melting head for Enceladus conditions. The three lines corre- velocity. In particular this is the case, if the power in- spond to three different values of the buoyancy corrected force, namely 0.1 N, put, hence the heat flow rate at the melting head is large. 10 N and 1000 N. The figure shows the results using equation (38). Small mass budgets for extraterrestrial exploration missions hence motivate the use of actuator systems, e.g. an In addition to the melting head surface temperatures, figure 8 ice screw as proposed in [11] in order to superpose the also shows the Stefan numbers over the the heat flow rate at probe’s intrinsic weight with an additional force. the melting head for three different buoyancy corrected forces. The curves have been produced using equation (38) with the • Wall heaters must be considered for the design of extrater- thermo-physical properties for water ice on Enceladus. Since restrial melting probes, at least if less power than 5kW the Stefan number Ste is just a different scaling of the melting is provided and the probe’s length is in the order of one head temperature, the curves are identical. meter. It can be seen that the temperature at the melting head TH increases with increasing heat flow rate Q˙ H and with decreasing • For very light melting probes without an additional actua- ∗ ∗ buoyancy corrected force F . For F = 0.1N, the temperature tor system, or if melting on icy bodies that are character- ◦ exceeds 100 C at Q˙ H ≈ 1500W, which at an ambient pres- ized by a low value of the gravitational acceleration, the sure of 1bar would induce boiling. Note, that a liquid-gaseous melting head can become dangerously hot, and counter- phase change is beyond the capabilities of the presented model measures have to be implemented. and therefore this curve should be interpreted with care. The melting head temperature for F∗ = 10N is between 15 ◦C at Our model is applicable only to very simple probe geome- Q˙ = 1000W and 94 ◦C at Q˙ = 4500W. Such high temper- H H tries, i.e. probe’s that have a flat melting head. In the future, atures can be dangerous if the maximum survival temperature it would be very interesting to also investigate the performance of the utilized electronic components in the melting head is ex- of parabolic melting heads, as it is not obvious, how changing ceeded. To overcome this risk, a temperature feedback control the head geometry interferes with the efficiency results we find might be implemented to switch off the melting head heaters in our paper. Though we believe that the adaption of our model if a critical temperature is exceeded. This in return, would to curved melting head geometries should be straight forward of course reduce the melting velocity and blow up the over- this remains to be future work. Another next step of our work all penetration time. The lowest temperatures are obtained for will be a rigorous assessment of the initial phase of the melting F∗ = 1000 N. For this relatively high buoyancy corrected force, mission, in which low pressure effects have to be taken into ac- the temperature varies between 3 ◦C at Q˙ = 1000W and28 ◦C H count and further complicate the physical process close to the at Q˙ = 5000 W. H surface.

5. Conclusions Acknowledgement We applied close-contact melting theory to develop an im- proved model for the design of melting probes that allows to study the role of gravitational acceleration as well as the melt- The project is supported by the Federal Ministry for Eco- ing heads surface temperature. Our model explicitly accounts nomic Affairs and Energy, Germany, on the basis of a decision for the physical processes in the melt film between the melting by the German Bundestag (50 NA 1502). It is part of the Ence- head and the water ice interface. It allows a direct quantification ladus Explorer initiative of the DLR Space Administration. 9 References

[1] Haldor WC Aamot. Heat transfer and performance analysis of a thermal probe for glaciers. Technical report, Cold Regions Research and Engi- neering Lab Hanover NH, 1967. [2] HALDOR WC Aamot. Self-contained thermal probes for remote mea- surements within an ice sheet. In International Symposium on Antarctic Glaciological Exploration (ISAGE), Hanover, NH, pages 63–68, 1968. [3] Adrian Bejan. Contact melting heat transfer and lubrication. In Advances in heat transfer, volume 24, pages 1–38. Elsevier, 1994. [4] Bernard J Hamrock et al. Fundamentals of fluid film lubrication. CRC press, 2004. [5] Mera Horne. Thermal probe design for europa sample acquisition. Acta Astronaut., 2017. [6] JC Jaeger. Conduction of heat in an infinite region bounded internally by a circular cylinder of a perfect conductor. Aust. J. Phys., 9(2):167–179, 1956. [7] Peter Kasser. Ein leichter thermischer eisbohrer als hilfsgerät zur instal- lation von ablationsstangen auf gletschern. Pure Appl. Geophys., 45(1): 97–114, 1960. [8] Norbert I. Kömle, Patrick Tiefenbacher, Peter Weiss, and Anastasiia Bendiukova. Melting probes revisited – ice penetration experiments under mars surface pressure conditions. Icarus, 2017. ISSN 0019-1035. doi: https://doi.org/10.1016/j.icarus.2017.08.006. [9] Norbert I Kömle et al. Melting probes revisited–ice penetration experiments under mars surface pressure conditions. Icarus, 2017. [10] Konstantinos Konstantinidis et al. A lander mission to probe subglacial water on saturn’s moon enceladus for life. Acta Astronaut., 106:63–89, 2015. [11] Julia Kowalski et al. Navigation technology for exploration of glacier ice with maneuverable melting probes. Cold Reg. Sci. Technol., 123:53–70, 2016. [12] Ralph D Lorenz. Thermal drilling in planetary ices: An analytic solution with application to planetary protection problems of radioisotope power sources. Astrobiology, 12(8):799–802, 2012. [13] Jonathan Lunine et al. Enceladus life finder: the search for life in a habitable moon. In EGU General Assembly Conference Abstracts, volume 17, 2015. [14] Jonathan I Lunine. Ocean worlds exploration. Acta Astronaut., 131:123– 130, 2017. [15] MK Moallemi, BW Webb, and R Viskanta. An experimental and ana- lytical study of close-contact melting. J. Heat Transfer, 108(4):894–899, 1986. [16] K Philberth. Geophysique-une methode pour mesurer les temperatures a linterieur dun inlandsis. Comptes rendus de l’Acadmie des sciences, 254 (22):3881, 1962. [17] K. Schüller and J. Kowalski. Spatially varying heat flux driven close- contact melting a lagrangian approach. Int. J. Heat Mass Transfer, 115 (Part B):1276 – 1287, 2017. ISSN 0017-9310. doi: https://doi.org/10. 1016/j.ijheatmasstransfer.2017.08.092. [18] Brent Sherwood. Strategic map for exploring the ocean-world enceladus. Acta Astronaut., 126:52–58, 2016. [19] RL Shreve. Theory of performance of isothermal solid-nose hotpoints boring in temperate ice. Journal of Glaciology, 4(32):151–160, 1962. [20] WC Stone et al. Progress towards an optically powered cryobot. Ann. Glaciol., 55(65):1–13, 2014. [21] Melanie Treffer et al. Preliminary studies concerning subsurface probes for the exploration of icy planetary bodies. Planet. Space Sci., 54(6): 621–634, 2006. [22] Stephan Ulamec et al. Access to glacial and subglacial environments in the solar system by melting probe technology. Rev. Environ. Sci. Biotech- nol., 6(1-3):71–94, 2007. [23] Dale P Winebrenner et al. Clean, logistically light access to explore the closest places on earth to europa and enceladus. In AGU Fall Meeting Abstracts, 2016. [24] Wayne Zimmerman et al. Cryobot: an ice penetrating robotic vehicle for mars and europa. In Aerospace Conference, 2001, IEEE Proceedings., volume 1, pages 1–311. IEEE, 2001.