Multiscale formulation of frictional contact mechanics at the pore scale

JOURNAL OF COMPUTATIONAL PHYSICS

卷:430

Article

Mehmani, Yashar¹  Castelletto, Nicola²  Tchelepi, Hamdi A.¹ 

[1] Stanford Univ, Dept Energy Resources Engn, Stanford, CA 94305 USA

[2] Lawrence Livermore Natl Lab, Atmospher Earth & Energy Div, Livermore, CA USA

关键词: Porous media;    Pore-scale modeling;    Multiscale method;    Contact mechanics;    Linear elasticity;    Preconditioning;   

DOI  :  10.1016/j.jcp.2020.110092

来源: Elsevier

【 摘 要 】

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. This work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free