Premium
A new monolithic Newton‐multigrid‐based FEM solution scheme for large strain dynamic poroelasticity problems
Author(s) -
Obaid A.,
Turek St.,
Heider Y.,
Markert B.
Publication year - 2016
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.5315
Subject(s) - solver , finite element method , multigrid method , mathematics , nonlinear system , tangent , mathematical optimization , computation , computer science , partial differential equation , mathematical analysis , algorithm , geometry , engineering , structural engineering , physics , quantum mechanics
Summary This paper presents a new efficient monolithic finite element solution scheme to treat the set of PDEs governing a 2D, biphasic, saturated theory of porous media model with intrinsically coupled and incompressible solid and fluid constituents for infinitesimal and large elastic deformation. Our approach, which inherits some of its techniques from CFD, is characterized by the following aspects: (1) it only performs operator evaluation with no additional Gateaux derivatives. In particular, the computations of the time‐consuming material tangent matrix are not involved here; (2) it solves the non‐linear dynamic problem with no restriction on the strength of coupling; (3) it is more efficient than the linear u v p solver discussed in previous works; (4) it requires weaker derivatives, and hence, lower‐order FE can be tested; and (5) the boundary conditions are reduced, solution independent and more convenient to apply than in the old u v p formulation. For the purpose of validation and comparison, prototypical simulations including analytical solutions are carried out, and at the end, an adaptive time stepping procedure is introduced to handle the rapid change in the numbers of nonlinear iterations that may occur. Copyright © 2016 John Wiley & Sons, Ltd.