Nonlinear regime of radially spreading extensional flows. Part 1. Newtonian fluids

Ice shelves that spread into the ocean can develop rifts that can trigger iceberg calving and enhance ocean-induced melting. Fluid mechanically, this system is analogous to an extensionally dominated radial spreading of a non-Newtonian fluid into a relatively inviscid and denser ambient fluid. Laboratory experiments have shown that rift patterns can emerge when the spreading fluid is shear thinning. Linear stability analysis supports these findings, predicting that while the instability mechanism is active in Newtonian fluids, it is suppressed by stabilising secondary-flow cellular vortices. Here, we explore the fully nonlinear evolution of a radially spreading Newtonian fluid, assessing whether large-amplitude perturbations could drive an instability. We use a quasi-three-dimensional numerical simulation that solves the full nonlinear shallow-shelf approximation, tracing the evolving fluid front, and validate it with known axisymmetric solutions and predictions from linear-stability theory. We find that large-amplitude perturbations induce nonlinear effects that give rise to non-axisymmetric patterns, including cusp-like patterns along the fluid front and complex secondary-flow eddies, which have neither been predicted theoretically nor observed experimentally. However, despite these nonlinear effects, large-amplitude perturbations alone are insufficient to induce rift-like patterns in Newtonian fluids. Strain-rate peaks at the troughs of the fluid front suggest that shear-thinning fluids may become more mobile in these regions, potentially leading to rift formation. This coincides with the likely weakening of stabilising forces as the fluid becomes more shear-thinning. These findings elucidate the critical role of nonlinear viscosity on the formation of rift-like patterns, which is the focus of Part 2 of this study.