A Unified Non-Local Fluid Transport Model for Heterogeneous Saturated Porous Media

Yuqi Sun,Haitao Yu

DOI: https://doi.org/10.1016/j.cma.2021.114294

IF: 6.588

2021-01-01

Computer Methods in Applied Mechanics and Engineering

Abstract:A unified non-local fluid transport model is proposed for heterogeneous saturated porous media. First, a non-local governing equation of fluid transport is presented based on the Non-ordinary State-based Peridynamics (NOSB-PD), and a novel flux vector state and a non-local Darcy’s law is proposed to illustrate non-local fluid transport effects in heterogeneous saturated porous media. The proposed flux state vector ensures the non-local governing equation of fluid transport model in heterogeneous saturated porous media degenerates into the local one when the horizon of each material point approaches to zero. A non-local pressure gradient at each material point can be obtained by averaging all the flux state vectors in its horizon, which unifies the weak and strong discontinuities of pressure in heterogeneous saturated porous media in a consistent way. Second, a variational formulation of the unified non-local fluid transport model is developed. This variational formulation enables the proposed model to deal with complex boundary conditions, which is difficult or even impossible for current fluid transport models based on Bond-based Peridynamics (BB-PD) and Ordinary State-based Peridynamics (OSB-PD). Third, a fully implicit algorithm combined with the Newton–Raphson method is introduced to solve the general non-linear problems in the heterogeneous saturated porous media. The penalty function method is employed to eliminate the zero-energy mode oscillation in NOSB-PD. Further, numerical examples demonstrate that the proposed model is accurate and can well capture the weak and strong discontinuities of pressure at the material interfaces and crack surfaces in the heterogeneous saturated porous media. Numerical examples also indicate that the model can effectively eliminate the zero-energy mode oscillation inherently rooted in the NOSB-PD, and thus reaches a stable numerical solution.