We describe the implementation of the complete Sea Level Equation (SLE) in a Finite Element (FE) self-gravitating 3D model. The procedure, originally proposed by Wu (2004), consists of iterating the solution of the SLE starting from a non self- gravitating model. At each iteration, the perturbation to the gravitational potential due to the deformation at the density interfaces is determined, and the boundary conditions for the following iteration are modified accordingly. We implemented the computation of the additional loads corresponding to the perturbations induced by glacial and oceanic forcings at the same iteration at which such forcings are applied. This implies an acceleration of the convergence of the iterative process that occurs actually in three to four iterations so that the complete procedure, for a 6,800 elements FE grid, can be run in about two hours of computing time, on a four-core 2.2 GHz Linux workstation. This spherical and self-gravitating FE model can be employed to simulate the deformation of the Earth induced by any kind of load (non necessarily of glacial origin) acting on the surface and/or internally.
Implementation of the Complete Sea Level Equation in a 3D Finite Elements Scheme: A Validation Study
SPADA, GIORGIO
2012
Abstract
We describe the implementation of the complete Sea Level Equation (SLE) in a Finite Element (FE) self-gravitating 3D model. The procedure, originally proposed by Wu (2004), consists of iterating the solution of the SLE starting from a non self- gravitating model. At each iteration, the perturbation to the gravitational potential due to the deformation at the density interfaces is determined, and the boundary conditions for the following iteration are modified accordingly. We implemented the computation of the additional loads corresponding to the perturbations induced by glacial and oceanic forcings at the same iteration at which such forcings are applied. This implies an acceleration of the convergence of the iterative process that occurs actually in three to four iterations so that the complete procedure, for a 6,800 elements FE grid, can be run in about two hours of computing time, on a four-core 2.2 GHz Linux workstation. This spherical and self-gravitating FE model can be employed to simulate the deformation of the Earth induced by any kind of load (non necessarily of glacial origin) acting on the surface and/or internally.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
