Nonlinear finite volume discretization of geomechanical problem

Elliptic differential operators describe a wide range of processes in mechanics relevant to geo-energy applications. Extensively used in reservoir modeling, the Finite Volume Method with TPFA can be consistently applied to discretize only a specific type of application under severe assumptions. In this paper, we introduce a positivity preserving Nonlinear Two Point Stress Approximation (NTPSA) based on the recently developed collocated Finite Volume scheme for linear elastic mechanics. The gradient reconstruction is different from the one used in Nonlinear TPFA, but a similar form of weighting scheme is employed to reconstruct the traction vector at each interface. The convergence of the scheme is tested with a homogeneous anisotropic stiffness tensor. The motivation behind the implementation of a new discretization framework in mechanics is to develop a uniform discretization technique preserving monotonicity for generic poromechanics applications.