This is a classical example to demonstrate the effect of hydromechanical coupling in a poroelastic medium. For more details we refer to a textbook [1], in which also the analytical solution is derived. As domain we consider a sphere, by symmetry we need to simulate only an octant.
Mesh
The boundary conditions of hydraulics are atmospheric pressure on the sphere surface and impermeable elsewhere. The boundary conditions of mechanics are an overburden (Neumann) applied as step load on the sphere surface at initial time $t=0$. The remaining sides are fixed in normal direction (Dirichlet).
Boundary conditions
The material is isotropic, linear elastic. Solid and fluid are assumed to be incompressible. In its initial state the sphere is not deformed and there is ambient pressure everywhere. A sudden load increase on the surface is instantly transferred on the pore pressure, whereas the solid needs time to deform, until it carries the load. Since the pore fluid is squeezed out of the outer layers first, they act like a tightening belt and consequently the pressure in the center rises, it may even exceed the value of the applied load. Finally the pore pressure approaches to ambient pressure.
All parameters are concluded in the following tables.
Property  Value  Unit 

Fluid density  $10^3$  kg/m$^3$ 
Viscosity  $10^{3}$  Pa$\cdot$s 
Porosity  $0.2$   
Permeability  $10\cdot 10^{12}$  m$^2$ 
Young’s modulus (bulk)  $10\cdot 10^6$  Pa 
Poisson ratio (bulk)  $0.1$   
Biot coefficient  $1.0$   
Solid density  $2.5\cdot 10^3$  kg/m$^3$ 
Overburden  $1000$  Pa 
Atmospheric pressure  $0$  Pa 
Property  Value  Unit 

Radius  $0.4$  m 
Finite Elements  $8741$  TaylorHood tetrahedral elements 
Time step  $10^{2}$  s 
Coupling scheme parameter  $0.7774$   
As predicted, the pressure in the center exceeds the applied load and then levels out.
Pressure at center of sphere
For more details about the staggered scheme we refer to the user guide  conventions.
[1] Verruijt, A. (2009): An introduction to soil dynamics. Springer Science and Business Media, DOI:https://doi.org/10.1007/9789048134410 https://link.springer.com/book/10.1007/9789048134410
This article was written by Dominik Kern. If you are missing something or you find an error please let us know.
Generated with Hugo 0.101.0
in CI job 282302

Last revision: November 22, 2022
Commit: [PL/THM] Implement freezing for temperature eq. 68ebbec
 Edit this page on