A Finite Element Simulation Of Three-dimensional Nonlinear Stokes Model For Laohugou Glacier

The numerical modeling of glacier and ice sheet evolution is a subject of growing interest,in part because of the potential for moadels to inform estimates of global sea level.In this dissertation,we use nonlinear Stokes equation to make a finite element numerical simulation for Laohugou Glacier,focusing on the determination of velocity field and pressure field.The viscosity coefficient is related to the speed,pressure,temperature and other associated variables.In this dissertation,we use the finite element method to solve the nonlinear Stokes equation,the high resolution mesh subdivision is helpful to obtain higher accuracy of numerical solution of velocity and pressure for the finite element numerical simulation.According to this view,in this dissertation,the three dimensional region of the glacier is constructed on the basis of the geothermal map,then,receiveing arbitrary effective subdivision of triangulation unit data with subdivision software.For the solution of the nonlinear Stokes equation,this dissertation conducted the tetrahedron for the area,and construct the second interpolation basis function,in order to further improve the smoothness of the numerial solution.that is,the velocity field and pressure field function space are constructed by the second interpolation basis function.Then,using Picard iterative method,we get the iterative scheme.in the end,the velocity and pressure values of Laohugou Glacier are obtained.The validity and rationality of the method is proved by numerical experiments.