Numerical analysis of 1-D compression wave propagation in saturated poroelastic media

A numerical procedure is presented for the simulation of 1-D compression wave propagation in saturated poroelastic media. The media are modelled as a two-phase system consisting of compressible fluid and solids. Viscous coupling forces resulting from the relative motion between phases are characterized as Darcy type. The numerical procedure can account for effects of axial strain, nonlinear material behaviour, and various drained and undrained boundary conditions. Time integration is carried out explicitly and isothermal conditions are assumed. The method is capable of modelling shock wave fronts without introducing artificial viscosity. Numerical results are in close agreement with analytical solutions for several simplified cases and indicate that mass coupling may have important effects on fluid velocity and wave speed. Corresponding effects on solid velocity and wave speed are much smaller. Numerical results also indicate that damping occurs in a saturated poroelastic column and is dependent on the value of hydraulic conductivity.