Phase transitions and 2D spherical convection in a large icy satellite

Three ice phase transitions, based on experimental and theoretical results, were incorporated into a 2D spherical convection model. These phase transitions are the exothermic Ice 1-2, the endothermic Ice 2-6, and the exothermic Ice 6-8 phase transitions. The fluid is assumed isoviscous and the anelastic liquid formulation was used. The object is an icy satellite whose physical characteristics resemble those of Ganymede or Callisto. Structural models, similar to those of Mueller and McKinnon, with different core sizes, i.e. with different degree of differentiation were studied. The shells, i.e. the mantle of the satellite, is heated both internally and from below to account for the decaying radiogenic heating and the heat flow from the solid core. The lower boundary of the mantle is rigid and isothermal, the upper boundary is isothermal. Calculations with rigid and shear stress free upper boundary are carried out in order to assess the role of the different boundary conditions. Two different Rayleigh numbers, depending on the assumed value of the viscosity, were used in the calculations and the thermal evolution of the satellite was studied. The suit of calculations presented demonstrates that phase transition cannot be ignored when the thermal evolution of a large icy satellite has to be studied. The importance of the construction of the realistic phase diagram for ice was also documented. On the other hand the degree of differentiation, simulated with the structures with different core's radius does not seem to drastically influence the overall behavior of the convective activity. Within the limits of our simulations, the post-accretional evolution seems to be independent on the degree of primordial differentiation.