Premium
A meshless method for axi‐symmetric poroelastic simulations: numerical study
Author(s) -
Ferronato Massimiliano,
Mazzia Annamaria,
Pini Giorgio,
Gambolati Giuseppe
Publication year - 2006
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1931
Subject(s) - solver , mathematics , finite element method , poromechanics , mathematical optimization , computer science , regularized meshless method , meshfree methods , subspace topology , projection (relational algebra) , parametric statistics , algorithm , singular boundary method , mathematical analysis , structural engineering , geology , engineering , statistics , geotechnical engineering , porous medium , porosity , boundary element method
A meshless model, based on the meshless local Petrov–Galerkin (MLPG) approach, is developed and implemented for the solution of axi‐symmetric poroelastic problems. The solution accuracy and the code performance are investigated on a realistic application concerning the prediction of land subsidence above a deep compacting reservoir. The analysis addresses several numerical issues, including the parametric selection of the optimal size of the local sub‐domains for the weak form and the nodal supports, the appropriate integration rule, and the linear system solver. The results show that MLPG can be more accurate than the standard finite element (FE) method on coarse discretizations, with its superiority decreasing as the nodal resolution increases. This is due to both a slower convergence rate and a progressively higher computational cost compared to FE. These drawbacks can be partially mitigated by improving the efficiency of the numerical integration and the system solver with the aid of projection techniques based on Krylov subspace methods. The outcome of the present analysis supports the development of coupled methods where a limited number of MLPG nodes are used to locally improve a FE solution. Copyright © 2006 John Wiley & Sons, Ltd.