Multicontinuum Generalized Multiscale Finite Element Method (MC-GMsFEM). Theory and applications to upscaling of two-phase flow

기술 요약: 다중 연속체 일반화 다중 스케일 유한 요소법 (MC-GMsFEM)

문제 정의
불균질 매질에서의 다중 스케일 문제는 일반적으로 미세 스케일 정보를 다중 스케일 기저 함수를 통해 포함하는 조대 스케일(coarse-scale) 모델을 구축함으로써 해결된다. 기존 방법들(예: GMsFEM)은 조대 격자 위에서 미세 스케일 해를 효율적으로 근사하지만, 물리적으로 일관된 거시적 미분 방정식을 도출하기 위한 체계적인 프레임워크는 일반적으로 결여되어 있다. 이러한 방법들의 조대 스케일 표현은 공간적 연속성이 부족한 자유도를 포함하며, 이는 근본적인 연속체 모델과의 연결성을 모호하게 만든다. 또한, 전통적인 다중 연속체 균질화 기술은 상호작용하는 연속체의 수가 공간적, 동적으로 변하는 복잡한 시스템에서 어려움을 겪으며, 구현하기 까다로운 특수한 셀 문제(cell problems)의 해결을 요구하는 경우가 많다.

본 논문은 정합적인(conforming) 다중 스케일 기저 함수를 구축하는 동시에 거시적 지배 방정식을 체계적으로 도출할 수 있는 통합 프레임워크를 개발함으로써 이러한 한계를 해결한다. 주요 적용 사례로는 여러 상호작용하는 연속체(이동성 및 포획된 상을 포함)가 발생하는 불균질 다상 매질 내의 이상(two-phase) 비혼합 흐름에 대한 업스케일링(upscaling)이 고려된다.

방법론: MC-GMsFEM 프레임워크
제안된 다중 연속체 일반화 다중 스케일 유한 요소법(MC-GMsFEM)은 다중 스케일 유한 요소 구축과 다중 연속체 표현을 통합한다. 핵심 방법론은 다음과 같이 진행된다:

1. 연속체 식별: 방법론은 연속체를 전역적으로 식별(예: 특성 함수 $\psi_i$를 통해)한다고 가정한다. 이상 흐름의 맥락에서, 이러한 연속체는 서로 다른 유체 상 또는 포획된 갱글리아(ganglia)와 같은 특정 영역에 대응한다.
2. 보조 함수 구축:
   • 연속체 함수 ($N_i$): 적분 제약 조건을 가진 국부 문제를 풀어 각 $i$번째 연속체를 나타내는 국부 함수를 구축한다. 구체적으로, $N_i$는 조대 블록 $K$ 내에서 연속체 $i$에 대해 단위 평균을 가지며, 다른 모든 연속체에 대해서는 영(zero)의 평균을 갖도록 정의된다: $\int_K N_i \psi_j = \delta_{ij} \int_K \psi_j$.
   • 기울기 보정 함수 ($M_i$): 미세 스케일 특징과 기울기를 고려하기 위해 보조 함수 $M_i$(또는 이들의 이산 형태인 $z_i = M_i \cdot \nabla \phi_0$)를 구축한다. 이들은 조대 블록 경계에서 영 경계 조건을 가지며 연속체에 대한 직교 제약을 만족하는 국부 문제를 충족한다.
3. 기저 함수 표현: 다중 스케일 기저 함수는 다음과 같이 표현된다:
   $$\phi^j_i = N_i \phi^j_0 + M_i \cdot \nabla \phi^j_0$$
   여기서 $\phi^j_0$는 표준 단위 분할(partition of unity) 함수이다. 이 표현은 조대 변수와 그 기울기를 분리한다.
4. 거시적 방정식 도출: 조대 스케일 해는 $u_H = \sum_i (N_i U_i + M_i \cdot \nabla U_i)$로 분해되며, 여기서 $U_i$는 연속체 변수이다. 이 분해를 지배 방정식(예: 타원형 방정식 또는 스토크스 흐름)의 변분 형식에 대입함으로써 저자들은 거시적 방정식을 도출한다. 이 방정식에는 포획된 영역에서의 효과를 포착하는 데 중요한 연속체 변수의 기울기와 관련된 항(Brinkman 유사 항)이 자연스럽게 포함된다.
5. 이상 흐름에 대한 적용: 프레임워크는 결합된 스토크스 흐름 및 수송 방정식에 적용된다. 이 방법은 길게 늘어진 포획된 유체 영역을 별개의 연속체로 취급한다. 비압축성 및 운동량 방정식은 단일 압력 및 다중 연속체 압력 공식 모두에 대해 도출되며, 결과적으로 연속체 간의 전달 항과 비압축성 제약을 고려하는 거시적 모델을 생성한다.

주요 결과
본 논문은 불균질 및 천공된 영역에서의 타원형 방정식과 천공 매질에서의 이상 흐름 시뮬레이션에 대한 수치 실험을 제시한다.

• 타원형 문제: MC-GMsFEM은 불균질 및 천공된 미세 구조에 대해 테스트되었다. 이 방법은 높은 정확도를 입증하였으며, 다양한 테스트 케이스 전반에서 업스케일링 및 다운스케일링된 해의 상대 오차가 일반적으로 20% 미만으로 유지되었다(표 1). 시각화 결과는 다운스케일링된 다중 스케일 해가 참조 미세 격자 해와 밀접하게 일치함을 확인시켜 주었다.
• 이상 흐름:
  • 갱글리아 탐지: 포획된 유체 영역(갱글리아)을 식별하고 이를 별개의 연속체로 취급하기 위해 너비 우선 탐색(BFS) 알고리즘이 사용되었다.
  • 정확도: 다양한 핑거(finger) 실현을 포함하는 대상 영역 1 및 2에 대한 테스트(테스트 A–E)에서, MC-GMsFEM 방식(전역 및 국부 변형 모두)은 참조 미세 격자 평균과 비교하여 작은 상대 오차를 가진 조대 스케일 해를 생성하였다(표 2–16).
  • 기울기 항의 영향: 기울기 항($M_i \cdot \nabla U_i$)이 포함된 모델과 포함되지 않은 모델을 비교한 결정적인 실험(테스트 C)을 수행하였다. 결과(표 17)에 따르면, 이 항들을 생략할 경우(모델을 표준 다중 연속체 Darcy 모델로 축소할 경우), 특히 포획된 갱글리아의 속도에서 오차가 크게 증가하였다(오차가 약 16%에서 약 96%로 증가). 이는 기울기 항이 포획된 영역을 정확하게 모델링하는 데 필수적임을 나타낸다.
  • 연속체 분리: 갱글리아를 별개의 연속체로 취급하는 것이 모델의 물리적 일관성을 개선한다는 것을 연구를 통해 보여주었으나, 특정 테스트 케이스에서 농도 오차에 미치는 수치적 영향은 미미한 것으로 나타난 반면, 속도 표현은 크게 개선되었다.

의의 및 주장
본 논문은 MC-GMsFEM이 복잡한 다중 스케일 시스템을 모델링하기 위한 강력하고 일반적인 접근법을 제공한다고 주장한다. 본 연구의 주요 의의는 다음과 같다:

1. 통합 프레임워크: 다중 스케일 이산화로부터 거시적 지배 방정식을 직접 도출하는 체계적인 경로를 제공하여, 수치 근사와 연속체 균질화 사이의 간극을 메운다.
2. 유연성: 기존의 다중 연속체 접근 방식과 달리, 이 프레임워크는 모든 구성에 대해 특수한 셀 문제를 해결할 필요 없이, 상호작용하는 연속체의 수가 공간적, 동적으로 변하는 불균질 매질을 다룰 수 있다.
3. 물리적 일관성: 표준 유한 요소 설정 내에 이 방법을 자연스럽게 내장함으로써, 이상 흐름에서 포획된 상의 역학을 포착하는 데 결정적인 것으로 입증된 필요한 기울기 항을 포함하는 거시적 방정식을 도출할 수 있게 한다.
4. 계산 효율성: 미세 스케일 정보를 보존하면서도 효율적인 조대 격자 시뮬레이션을 가능하게 하여, 전체 미세 스케일 모델링에 대한 계산 효율적인 대안을 제공한다.

저자들은 이 방법이 기존 균질화 기술보다 장점을 제공하지만, 국부 문제를 풀지 않고도 쉽게 사전 계산할 수 있는 폐쇄형(closed form)으로 거시 방정식의 기울기 항 계수를 명시적으로 정의하지는 않는다는 점을 겸허히 언급하였으나, 변화하는 연속체와 일반적인 불균질성에 대한 이 방법의 적응성이 주요 강점임을 강조하였다.