The discovery of vast subsurface oceans hidden under kilometers of ices on icy moons in our Solar System has sparked worldwide interests in ascertaining their potential habitability. In the case of Saturn’s moon Enceladus, supersonic plumes of water vapour and icy grains have been observed by the Cassini mission spewing from the surface, giving us indirect knowledge of the composition of this subsurface ocean. The exact mechanisms of the plumes however, and their effect on the composition of the ejected matter has yet to be clearly understood. The focus of this study is to experimentally investigate physical characteristics of the plumes located at the South Polar Terrain (SPT) of Enceladus. Using facilities at TU Delft faculty, we simulate experimentally the topology of the ice crevasses and the subsurface ocean with a narrow channel mounted atop a liquid water reservoir placed inside a vacuum chamber. We inquire upon the dependence of the channel walls temperature on the plume’s exhaust velocity. Using a straight channel, our results show that colder wall temperatures enable a saturated water vapour flow with a minima 1.5-3 % solid fraction and vent velocities reaching around 400-500 m/s. The data ranges for velocities and solid fraction extrapolated from the Cassini data (550-2000 m/s and 7-70 %) thus cannot be explained by straight channel models. Using a channel with an expansion ratio of 1.73 however, the measured supersonic plume velocity becomes comparable to some of the in situ Mach number determined at Enceladus. Using a method based on the energy conservation law, it is possible to extrapolate from our experimental data some plausible geometries of the ice crevasses of Enceladus. This work lays the ground work for a coming comprehensive parametric study of the channel geometry and its effect on exhaust Mach number, temperature and solid fraction.@en

Supersonic plumes of water vapour and icy particles have been observed by the Cassini spacecraft during several flybys over Enceladus. These plumes originate from the Tiger Stripes located in the South Polar Terrain (SPT), and indicate the presence of a subsurface ocean under the icy crust which is salty and contains complex organic molecules. Other characteristics of the plumes, such as the vent temperature, mass flow rate, velocity and mass fraction of icy particles can be used to determine the conditions in the channel, linking the subsurface ocean to the icy surface. In this paper, we developed a fluid dynamics model that accounts for nucleation, particle growth, wall accretion and sublimation. The channel behaves similarly to a converging–diverging nozzle, which forms supersonic plumes due to a pressure difference between the reservoir where the subsurface ocean is located and the exosphere. The geometry of the channel and its evolution with accretion of gas and sublimation of ice are studied to reproduce the characteristics of the plumes observed by Cassini. We first performed a parameter study on the channel geometry to determine how it influences the plumes’ velocity, solid fraction and exit temperature. Our results show that the size of the icy particles is primarily dependent on the length of the channel, indicating that large particles (∼75μm) must originate from within a kilometer below the surface, while smaller particles (∼3μm) can originate from only hundreds of meters below the surface. We further show that the velocity of the flow, exit temperature and nucleation depend directly on the exit-to-throat size ratio. We find that the channel geometry evolves within a few tens of hours until an equilibrium is reached, when considering the accretion of gas to the walls, or sublimation of ice from the walls. As the channel closes due to accretion, the flow becomes thinner, which in turn reduces accretion. After around 70 h, the accretion is sufficiently slowed such that the geometry does not evolve anymore. This equilibrium geometry produces higher Mach numbers and a larger particle size and solid fraction compared to the initial geometry.@en