Finite element models of non-Newtonian crater relaxation

The effects of a non-Newtonian rheology on the profiles of relaxing craters (such as those seen on the surfaces of the icy Galilean and Saturnian satellites) were studied. Two-dimensional finite element simulations of non-Newtonian viscous flow were performed, and the results were compared with those associated with Newtonian rheology. Viscous relaxation of craters in a non-Newtonian medium was significantly different from that in a Newtonian medium. Crater rims are observed to relax at a more rapid rate in a non-Newtonian region as a result of the movement of the low viscosity region to underneath the crater rim after the initial relaxation of the bowl. Significant differences are also found when central depth is plotted as a function of time. For a Newtonian medium, crater relaxation is exponential in form. In contrast, non-Newtonian crater relaxation is initially rapid, in response to the large initial stresses and small viscosities; however, as stresses decrease, this relaxation becomes extremely gradual.