On a Doubly Nonlinear Parabolic Obstacle Problem Modelling Ice Sheet Dynamics

This paper deals with the weak formulation of a free (moving) boundary problem arising in theoretical glaciology. Considering shallow ice sheet flow, we present the mathematical analysis and the numerical solution of the second order nonlinear degenerate parabolic equation modelling, in the isothermal case, the ice sheet non-Newtonian dynamics. An obstacle problem is then deduced and analyzed. The existence of a free boundary generated by the support of the solution is proved and its location and evolution are qualitatively described by using a comparison principle and an energy method. Then the solutions are numerically computed with a method of characteristics and a duality algorithm to deal with the resulting variational inequalities. The weak framework we introduce and its analysis (both qualitative and numerical) are not restricted to the simple physics of the ice sheet model we consider nor to the model dimension; they can be successfully applied to more realistic and sophisticated models related to other geophysical settings.

[1]  A. Bermúdez,et al.  Duality methods for solving variational inequalities , 1981 .

[2]  Alfredo Bermúdez,et al.  Numerical solution of a three-dimensional thermoelectric problem taking place in an aluminium electrolytic cell , 1993 .

[3]  Petra Wittbold,et al.  On mild and weak solutions of elliptic-parabolic problems , 1996, Advances in Differential Equations.

[4]  Jesús Ildefonso Díaz Díaz,et al.  On a nonlinear parabolic problem arising in some models related to turbulent flows , 1994 .

[5]  Lawrence C. Evans,et al.  Instantaneous shrinking of the support of nonnegative solutions to certain nonlinear parabolic equations and variational inequalities , 1979 .

[6]  L. Lliboutry.  Very Slow Flows of Solids , 1987 .

[7]  O. Pironneau.  On the transport-diffusion algorithm and its applications to the Navier-Stokes equations , 1982 .

[8]  N. Calvo,et al.  Numerical approach of temperature distribution in a free boundary model for polythermal ice sheets , 1999, Numerische Mathematik.

[9]  W. Peltier,et al.  Terminating the 100 kyr ice age cycle , 1997 .

[10]  P. Lions,et al.  On the Cauchy problem for Boltzmann equations: global existence and weak stability , 1989 .

[11]  M. Chambat,et al.  Characteristics method for the formulation and computation of a free boundary cavitation problem , 1998 .

[12]  Ricardo H. Nochetto,et al.  A characteristic Galerkin method with adaptive error control for the continuous casting problem , 2000 .

[13]  Ricardo H. Nochetto,et al.  Approximation of Degenerate Parabolic Problems Using Numerical Integration , 1988 .

[14]  J. Durany,et al.  An elastohydrodynamic coupled problem between a piezoviscous Reynolds equation and a hinged plate model , 1997 .

[15]  Jesús Ildefonso Díaz Díaz,et al.  On a degenerate parabolic/hyperbolic system in glaciology giving rise to a free boundary , 1999 .

[16]  J. Durany,et al.  Numerical computation of ice sheet profiles with free boundary models 1 Partially supported by Resea , 2000 .

[17]  A. Fowler,et al.  A theory of ice-sheet surges , 1998, Journal of Glaciology.

[18]  O. Pironneau,et al.  Finite Elements and characteristics for some parabolic-hyperbolic problems , 1983 .

[19]  Jesús Ildefonso Díaz Díaz,et al.  On the Boussinesq system with nonlinear thermal diffusion , 1997 .

[20]  S. Antontsev,et al.  The support shrinking properties for solutions of quasilinear parabolic equations with strong absorption terms , 1995 .

[21]  Rodolfo Bermejo,et al.  A Multigrid Algorithm for the p-Laplacian , 1999, SIAM J. Sci. Comput..

[22]  Andrew C. Fowler,et al.  Mathematical Models in the Applied Sciences , 1997 .

[23]  Pekka Neittaanmäki.  Numerical Methods for Free Boundary Problems , 1991 .

[24]  J. Carrillo,et al.  Uniqueness of Renormalized Solutions of Degenerate Elliptic-Parabolic Problems , 1999 .

[25]  L. Véron,et al.  Local vanishing properties of solutions of elliptic and parabolic quasilinear equations , 1985 .

[26]  Andrew C. Fowler,et al.  Modelling ice sheet dynamics , 1992 .

[27]  J. Durany,et al.  Numerical simulation of a lubricated hertzian contact problem under imposed load , 2002 .

[28]  J. Lions,et al.  Les inéquations en mécanique et en physique , 1973 .