A coupled multivalued model for ice streams and its numerical simulation

^** Email: nati{at}dma.uvigo.es^*** Email: durany{at}dma.uvigo.es^**** Email: anaisabel.munoz{at}urjc.es^***** Email: emanuele.schiavi{at}urjc.es^****** Email: carlosv{at}udc.es This paper deals with the numerical solution of a non-linearmodel describing a free-boundary problem arising in modern glaciology.Considering a shallow, viscous ice sheet flow along a soft,deformable bed, a coupled non-linear system of differentialequations can be obtained. Particularly, an obstacle problemis then deduced and solved in the framework of its complementarityformulation. We present the numerical solution of the resultingmultivalued system modelling the ice sheet non-Newtonian dynamicsdriven by the underlying drainage system. Our numerical resultsshow the existence of fast ice streams when positive wave-likeinitial conditions are considered. The solutions are numericallycomputed with a decoupling iterative method and finite-elementtechnique. A duality algorithm and a projected Gauss–Seidelmethod are the alternatives used to cope with the resultingvariational inequality while the explicit treatment, Newtonmethod or a duality method are proposed to deal with the non-linearsource term. Finally, the numerical solutions are physicallyinterpreted and some comparisons among the numerical methodsare then discussed.