[논문 리뷰] Simulation of Deformation and Flow in Fractured, Poroelastic Materials

Katja K. Hanowski, Oliver Sander|arXiv (Cornell University)|2016. 06. 18.

Advanced Numerical Methods in Computational Mathematics참고 문헌 3인용 수 8

한 줄 요약

이 논문은 단일 유한 길이의 균열을 가진 균열이 있는 다공성 매체에서 유체 흐름과 다공성 탄성 변형을 동시에 시뮬레이션하기 위한 새로운 변분 형식과 수치 방법을 제시한다. 매트릭스 기계학에 대해 XFEM을, 균열 흐름에 대해 저차원 Darcy 모델을 사용하며, 균열 끝부분에서 균열 두께가 0에 수렴함에도 불구하고 가중 치수 공간에서 해의 존재성을 증명하고, 2D 및 3D 시뮬레이션에서 최적 수렴 속도를 보이며, 빠른 서브구조화 솔버 수렴 성능을 입증한다.

ABSTRACT

We introduce a coupled system of PDEs for the modeling of the fluid-fluid and fluid-solid interaction in a poroelastic material with a single static fracture. The fluid flow in the fracture is modeled by a lower-dimensional Darcy equation, which interacts with the surrounding rock matrix and the fluid it contains. We explicitly allow the fracture to end within the domain, and the fracture width is an unknown of the problem. The resulting weak problem is nonlinear, elliptic and symmetric, and can be given the structure of a fixed-point problem. We show that the coupled fluid-fluid problem has a solution in a specially crafted Sobolev space, even though the fracture width cannot be bounded away from zero near the crack tip. For numerical simulations, we combine XFEM discretizations for the rock matrix deformation and pore pressure with a standard lower-dimensional finite element method for the fracture flow problem. The resulting coupled discrete system consists of linear subdomain problems coupled by nonlinear coupling conditions. We solve the coupled system with a substructuring solver and observe very fast convergence. We also observe optimal mesh dependence of the discretization errors even in the presence of crack tips.

연구 동기 및 목표

• 균열 내 유체 흐름과 주변 다공성 탄성 매트릭스의 변형 사이의 비선형 결합을 모델링하고, 유체-유체 및 유체-고체 상호작용을 포함한다.
• 도메인 내에 끝나는 균열, 즉 균열 끝부분에서 균열 두께가 0이 되어 해의 특이성이 발생하는 문제에 대응한다.
• 균열 두께가 0으로부터 유한하게 떨어져 있지 않은 경우에도 특별히 구성된 가중 치수 공간에서 연립 방정식의 약한 해의 존재성을 확립한다.
• 매트릭스 변위와 압력에 대해 XFEM을, 균열 흐름에 대해 표준 FEM을 사용하는 단일화된, 일치하는 이산화를 개발하여 특이성에도 불구하고 최적 수렴 속도를 보장한다.
• 서브구조화 솔버를 구현하고 검증하여 2D 및 3D 설정 모두에서 비선형 결합 시스템에 대해 빠른 수렴 성능을 확보한다.

제안 방법

• 매트릭스 내 3D 다공성 탄성성과 균열 내 1D Darcy 모델을 조합한 PDE 연립방정식을 설정하며, 균열 두께를 미지 변수로 포함한다.
• 비선형 결합을 다루기 위해 고정점 반복 프레임워크를 사용하며, 균열 투과도는 미지의 균열 두께에 반비례한다.
• 균열 끝부분에서 응력과 유체 압력 특이성을 정확히 포착하기 위해 특이 함수(예: $ r^{1/2} $)를 사용한 XFEM 보정을 적용하며, 메esh 정렬이 필요하지 않다.
• 매트릭스 문제를 XFEM 보정을 사용한 일치하는 유한요소로 이산화하며, 변위와 기포 압력을 다루고, 균열은 저차원 표준 FEM으로 처리한다.
• 결과로 생긴 연립 방정식을 서브도메인 문제를 별도로 다루고 비선형 결합 조건을 효율적으로 강제하는 서브구조화 솔버로 해결한다.
• 균열 두께가 균열 끝부분에서 0에 수렴함에도 불구하고 해의 존재성을 엄밀히 증명하기 위해 가중 치수 공간을 구성한다.

실험 결과

연구 질문

• RQ1 균열 두께가 균열 끝부분에서 0에 수렴하고 0으로부터 유한하게 떨어져 있지 않은 경우에도, 다공성 탄성-균열 흐름 연립계에 대해 약한 해가 존재할 수 있는가?
• RQ2 균열 끝부분 특이성 존재 시, 역투과도를 통한 균열 흐름과 매트릭스 변형 사이의 비선형 결합을 정확히 이산화할 수 있는가?
• RQ3 특이성이 존재하고 균열이 도메인 내에 끝나는 경우, 수치적 스킴의 수렴 거동은 어떠한가?
• RQ4 서브구조화 솔버는 단일화된 형식으로 유도된 비선형으로 결합된 시스템을 효율적으로 처리할 수 있는가?
• RQ5 적절한 보정 함수를 사용한 XFEM의 적용이 균열 끝부근에서 변위, 압력 및 균열 흐름에 대해 최적 수렴 속도를 유지하는가?

주요 결과

• 균열 두께가 균열 끝부분에서 0에 수렴함에도 불구하고 특별히 구성된 가중 치수 공간에서 연립 다공성 탄성-균열 흐름 문제에 대해 약한 해가 존재한다.
• 매트릭스 압력의 $ L^2 $ 오차는 $ ext{O}(h^3) $ 속도로 감소하며, 이는 초수렴성을 시사하며, 주로 압력 기울기 특이성을 정확히 재현하기 때문일 것이다.
• 균열 압력의 $ L^2 $ 오차는 $ ext{O}(h^2) $ 속도로 감소하며, 예상보다 우수한 성능을 보이며, 이는 매트릭스 압력에서 기인한 결합 항의 향상된 정규성 때문일 수 있다.
• 변위와 체적 압력의 $ L^2 $ 및 $ H^1 $ 오차는 각각 최적 수렴 속도 $ ext{O}(h^2) $ 및 $ ext{O}(h) $ 를 나타내어, XFEM 이산화의 강건성을 확인한다.
• 서브구조화 솔버는 모든 테스트 케이스에서 매우 적은 반복 수로 수렴하여 비선형으로 결합된 시스템에 대해 높은 효율성을 보였다.
• 3D 시뮬레이션은 이 방법이 실제 기하 구조에 적용 가능함을 보이며, 물리적으로 일관된 결과를 보여주는데, 유체 유출이 균열 개방과 균열 앞선에서 최대 von Mises 응력 발생을 유도한다.