균열이 있는 다공성 매질 내 유동 해석을 위한 비정합 메쉬 모델의 비교 및 적용

초록

본 연구는 균열 다공성 매질 내 단상 유동을 모델링하기 위해 이중 라그랑주 승수법과 변분 전이 기법을 결합한 통합 프레임워크를 제시합니다. 복잡한 비정합 기하구조와 다양한 균열 모델(하이브리드 차원, 등차원)을 효율적으로 처리할 수 있으며, 2D 및 3D 벤치마크에서 우수한 성능을 입증했습니다. 복잡한 고투수성 균열망에는 내장형 하이브리드 차원 모델이, 두꺼운 차단 균열에는 등차원 모델이 더 적합함을 보여줍니다.

상세 분석

본 논문의 핵심 기술적 기여는 서로 다른 차원과 기하학을 가진 균열-매질 도메인 간의 정확하고 효율적인 커플링을 가능하게 하는 통합 수치 프레임워크를 제시한 점에 있습니다. 기존의 일치 메쉬(conforming mesh) 기반 방법은 복잡한 3D 균열망에서 메쉬 생성이 매우 어려운 문제가 있었습니다. 이를 해결하기 위해 저자들은 두 가지 주요 비정합(non-conforming) 기법, 즉 내장형(embedded) 접근법과 모르타르(mortar) 도메인 분해법을 이중 라그랑주 승수 공간을 통해 하나의 체계로 통합했습니다.

여기서 ‘이중(Dual)’ 라그랑주 승수의 사용이 중요한데, 이는 기존의 프라이멀(Primal) 공간을 사용할 때 발생하는 부정부호(indefinite)의 새들-포인트(saddle-point) 문제를 대칭 양정부호(symmetric positive-definite) 문제로 변환시킵니다. 이로 인해 켤레 기울기법(Conjugate Gradient)과 같은 효율적인 솔버와 멀티그리드 프리컨디셔너를 사용할 수 있어 대규모 계산의 성능과 안정성이 크게 향상됩니다.

변분 전이(Variational Transfer)의 핵심은 L2-투영을 기반으로 비정합 메쉬 간의 필드 정보를 투영하는 것입니다. 이는 기하학적으로 일치하지 않는 메쉬 간에 변분적으로 일관된(variationally consistent) 정보 전달을 보장하며, 병렬 계산에 매우 적합합니다. 또한 논문에서는 오차 제어와 자유도 감소를 위해 이 커플링 프레임워크 내에 비정합 적응 메쉬 정련(AMR)을 통합하는 방법을 제시합니다.

수치 예제를 통해 저자들은 얇은 균열(코차원 1)을 표현하는 하이브리드 차원 모델이 메쉬 생성 효율성과 계산 비용 측면에서 매우 우수하지만, 두꺼운 균열이나 투과도가 매우 낮은 ‘차단’ 균열의 경우 정확도가 떨어짐을 지적합니다. 후자의 경우 균열을 매질과 동일한 차원(등차원)의 볼륨으로 표현하고 모르타르 방법으로 커플링하는 것이 더 정확한 결과를 제공함을 입증합니다. 이는 실용적 관점에서 문제의 물리적 특성(균열의 투과도, 두께)에 따라 최적의 모델링 기법을 선택할 수 있는 지침을 제시한다는 점에서 의미가 큽니다.