A new mixed finite element method for poro‐elasticity |
| |
| Authors: | Maria Tchonkova John Peters Stein Sture |
| |
| Institution: | 1. 9009 Great Hills Trail, Suite 224, Austin, TX 78759, U.S.A.;2. US Army Engineer Research and Development Center, Vicksburg, MS 39190, U.S.A.;3. Department of Civil, Environmental and Architectural Engineering, Campus Box 428, University of Colorado at Boulder, Boulder, CO 80309, U.S.A. |
| |
| Abstract: | Development of robust numerical solutions for poro‐elasticity is an important and timely issue in modern computational geomechanics. Recently, research in this area has seen a surge in activity, not only because of increased interest in coupled problems relevant to the petroleum industry, but also due to emerging applications of poro‐elasticity for modelling problems in biomedical engineering and materials science. In this paper, an original mixed least‐squares method for solving Biot consolidation problems is developed. The solution is obtained via minimization of a least‐squares functional, based upon the equations of equilibrium, the equations of continuity and weak forms of the constitutive relationships for elasticity and Darcy flow. The formulation involves four separate categories of unknowns: displacements, stresses, fluid pressures and velocities. Each of these unknowns is approximated by linear continuous functions. The mathematical formulation is implemented in an original computer program, written from scratch and using object‐oriented logic. The performance of the method is tested on one‐ and two‐dimensional classical problems in poro‐elasticity. The numerical experiments suggest the same rates of convergence for all four types of variables, when the same interpolation spaces are used. The continuous linear triangles show the same rates of convergence for both compressible and entirely incompressible elastic solids. This mixed formulation results in non‐oscillating fluid pressures over entire domain for different moments of time. The method appears to be naturally stable, without any need of additional stabilization terms with mesh‐dependent parameters. Copyright © 2007 John Wiley & Sons, Ltd. |
| |
| Keywords: | finite element method Biot consolidation porous medium elasticity Darcy flow |
|
|
