변형 가능한 다공성 매체에서의 2상 흐름을 위한 순차 불연속 Galerkin 방법

초록

본 논문은 변형 가능한 다공성 매체 내 이소열 2상 흐름을 모델링하는 연속 방정식에 대해, 내부 페널티 불연속 Galerkin(DG) 공간과 순차적 시간 전진 방식을 결합한 수치 스키마를 제안한다. 압력(p_w, p_o)과 변위(u)를 기본 미지변수로 두고, 각 시간 단계에서 질량 보존식과 운동량 방정식을 차례대로 풀어 계산 비용을 크게 낮추면서도 해의 존재성 및 유일성을 증명한다. 3차원 실험을 통해 최적 수렴률과 이질성(불연속 모세관압, 투과성 변동) 및 외부 하중에 대한 강인성을 확인하였다.

상세 분석

이 연구는 Biot 이론에 기반한 2상 흐름‑포엘라스틱 연동 모델을 출발점으로 삼는다. 습윤상(p_w)과 비습윤상(p_o)의 압력과 포화도(s_w)를 변수로 두고, 모세관압 p_c(s_w)와 상대 투과성 k_r_i(s_w)를 비선형 함수로 설정한다. 연속 방정식(2)·(3)과 변위 방정식(4)은 각각 질량 보존과 quasi‑static 탄성 균형을 나타내며, C_i(p_o,p_w) 계수는 포화도와 모세관압에 의해 결정되는 복합 비선형 함수를 포함한다. 논문은 C_1·C_3이 비음(≥0)임을 이용해 순차적 스키마를 설계한다.

시간 이산화는 첫 단계 τ_0를 매우 작게 잡아 초기값을 정확히 초기화하고, 이후 일정 시간 간격 τ로 진행한다. 각 단계는 세 개의 서브스텝으로 구성된다. ① 습윤상 압력 P_{n+1}^w를 구할 때, C_1·ΔP_w/τ + C_2·ΔP_o/τ 형태의 질량 항과 내부 페널티 DG 형태의 확산 항 a(λ_n^w K;·,·)을 결합하고, 변위 변화에 대한 결합항 α b_u(S_n^w;·,·)을 포함한다. ② 비습윤상 압력 P_{n+1}^o는 유사한 구조로, C_3·ΔP_o/τ + C_4·ΔP_w/τ와 a(λ_n^o K;·,·), 그리고 (1‑S_n^w)에 대한 결합항을 사용한다. ③ 변위 U_{n+1}는 탄성 연산자 c(·,·)와 압력‑변위 결합 b_p(·,·)에 더해, γ·(U_{n+1}‑U_n)/τ 형태의 인공 점성(stabilization) 항을 도입한다. γ는 충분히 큰 양수로 선택해 DG 형태의 강제성을 확보한다.

DG 공간은 3차원 테트라hedral 메쉬에 대해 P1(선형) 다항식으로 정의되며, 내부 페널티 파라미터 σ_p, σ_u는 대칭 형태(ε = –1)일 경우 충분히 크게 설정한다. 점프·평균 연산자를 이용한 a(·)와 c(·)는 표준 interior‑penalty DG 스킴을 따르며, b_u와 b_p는 각각 발산·그라디언트 연산에 대한 DG 변형이다.

수학적 분석에서는 a와 c가 coercive함을 보이는 Lemma 1을 제시하고, λ_w, λ_o가 양의 하한을 갖는 경우(즉, 상대 투과성이 0이 아닌 경우) Proposition 1을 통해 각 서브스텝의 선형 시스템이 유일해짐을 증명한다. 존재성은 초기값을 L² 투영으로 정의하고, 시간 전진 과정에서 선형성에 의존해 귀납적으로 확보한다.

수치 실험에서는 (i) 제조해(solution of manufactured) 문제를 통해 L² 및 broken‑gradient 노름에서 2차 수렴률을 확인하고, (ii) 전통적인 McWhorter 문제(1‑D 역류)와 3‑D 이질 매질(불연속 모세관압, 급격한 투과성 변동) 사례를 적용했다. 특히, γ=10, σ_p=20, σ_u=14와 같은 파라미터 선택이 수렴과 안정성에 결정적임을 보여준다. 복합 매질에 대한 하중 적용 실험에서는 변위와 압력 분포가 물리적으로 기대되는 형태로 전파되는 것을 확인했으며, GMRES‑LU 프리컨디셔너를 이용한 선형 해석이 1~2 회 반복으로 수렴하는 효율성을 입증했다. 전체적으로 제안된 순차 DG 스키마는 완전 암시적(fully‑implicit) 방법에 비해 연산 비용을 크게 절감하면서도, 비선형 결합계수를 포함한 실제 지질공학 문제에 적용 가능한 정확도와 강인성을 제공한다.