Analysis and Finite Element Approximation of a Nonlinear Stationary Stokes Problem Arising in Glaciology英文文献资料.docVIP

HindawiPublishingCorporation AdvancesinNumericalAnalysis Volume2011,ArticleID164581,24pages doi:10.1155/2011/164581 ResearchArticle AnalysisandFiniteElement Approximation ofaNonlinearStationaryStokes ProblemArisinginGlaciology GuillaumeJouvet1,2andJacquesRappaz 2 1Department ofMathematics andComputerScience,FreeUniversity ofBerlin,14195Berlin,Germany 2ChairofNumerical AnalysisandSimulation, EPFL,1015Lausanne, Switzerland Correspondence shouldbeaddressedtoGuillaume Jouvet,guillaume.jouvet@fu-berlin.de Received13September 2011;Accepted 25October2011 Academic Editor:MicheleBenzi Copyright q2011G.JouvetandJ.Rappaz. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction inanymedium,providedtheoriginalworkisproperlycited. The aim of this paper is to study a nonlinear stationary Stokes problem with mixed boundary conditions that describes the ice velocity and pressure ?elds of grounded glaciers under Glen’s ?owlaw.Usingconvexanalysisarguments, weprovetheexistence andtheuniqueness ofaweak solution. A ?nite element method is applied with approximation spaces that satisfy the inf-sup condition, and a priori error estimates are established by using a quasinorm technique. Several algorithms including Newton’s method are proposed to solve the nonlinearity of the Stokes problem and areproved to be convergent. Our results aresupported by numerical convergence studies. 1.Introduction Inthispaperweconsideramodelproblemthatiscommonlyusedbyglaciologiststocompute themotionofglaciers.Iceisassumedtobeanincompressiblenon-Newtonian?uidgoverned byGlen’slaw1.Glen’slawandthemassmomentumequationleadtoanonlinearstationary Stokesproblemwithastrain-dependentviscosity. Glacier models based on Glen’s law have already been studied by several authors. However, all of them have considered a simpli?ed model, called ?rst-order approximation 2. This model is obtained by rewriting the S