[논문 리뷰] Upscaling of unsaturated flow in fractured porous media

Florian List, Kundan Kumar외 2명|arXiv (Cornell University)|2018. 07. 16.

Advanced Mathematical Modeling in Engineering참고 문헌 48인용 수 5

한 줄 요약

이 논문은 박층 다공성 매체에서 불포화 유동에 대한 효과적 모델을 두꺼운 균열(폭 ε → 0)의 극한을 취하여 엄밀하게 유도한다. 이로 인해 유도된 다각도의 효과적 모델은 밀도 비율과 투수도 비율이 ε에 따라 어떻게 척도화되는지에 따라 달라진다. 분석을 통해 다양한 효과적 방정식, 예를 들어 균열 인터페이스에서의 1차원 리처드스 방정식 등이 κ 및 λ의 값에 따라 도출되며, ε ≤ 10⁻³일 경우 수치적으로 검증되었고 오차가 10⁻⁵ 이하로 매우 높은 정확도를 보였다.

ABSTRACT

In this work, we consider a mathematical model for flow in a unsaturated porous medium containing a fracture. In all subdomains (the fracture and the adjacent matrix blocks) the flow is governed by Richards' equation. The submodels are coupled by physical transmission conditions expressing the continuity of the normal fluxes and of the pressures. We start by analyzing the case of a fracture having a fixed width-length ratio, called $\varepsilon > 0$. Then we take the limit $\varepsilon o 0$ and give a rigorous proof for the convergence towards effective models. This is done in different regimes, depending on how the ratio of porosities and permeabilities in the fracture, respectively matrix scale with respect to $\varepsilon$, and leads to a variety of effective models. Numerical simulations confirm the theoretical upscaling results.

연구 동기 및 목표

• 박층 다공성 매체에서의 불포화 유동 시뮬레이션에 하위차원 균열 모델을 사용하는 것의 수학적 정당성을 제시하기 위해 얇은 균열 기하학에서 이를 엄밀하게 유도하는 것.
• 균열 두께 ε에 따른 밀도 및 투수도 비율의 척도화 방식에 따라 다양한 효과적 모델이 어떻게 도출되는지 규명하는 것.
• 양의 두께를 가진 균열 문제(ε-문제)가 ε → 0 극한에서 효과적 모델로 수렴하는지 엄밀히 증명하는 것.
• 지하 유동 시뮬레이션에서 사용되는 이산 균열 모델링 접근법에 물리적으로 일관된 기반을 제공하는 것.
• 작은 ε에 대해 높은 정확도로 수치 시뮬레이션을 통해 이론적 결과를 검증하는 것

제안 방법

• 매트릭스 블록과 두께 ε를 가진 얇은 균열에서 리처드스 방정식을 사용하여 문제를 수립하고, 인터페이스에서 압력과 법선 유량이 연속임을 보장한다.
• 척도 가정 도입: 밀도 비율 ∝ ε^κ 및 투수도 비율 ∝ ε^λ, 여기서 κ > -1 및 λ < 1.
• ε → 0 극한에서 효과적 모델을 도출하기 위해 엄밀한 점근적 분석 및 균질화 기법을 적용한다.
• κ 및 λ의 값에 따라 다양한 효과적 모델을 도출하며, λ < 1 및 κ > -1일 경우 균열 인터페이스에서 1차원 리처드스 방정식이 도출됨.
• 약한 공식화 및 에너지 추정을 사용하여 ε-문제의 약한 해의 존재성과 유일성을 증명하고, 효과적 모델로의 수렴성을 입증한다.
• 반-겐وخ텐–무레이엠 매개변수를 사용하여 수치 시뮬레이션을 수행하고, ε-해와 효과적 모델을 비교하여 하위영역 전체에서 L² 오차를 계산한다.

실험 결과

연구 질문

• RQ1 균열 두께가 작은 양의 값을 가지는 경우, 밀도 및 투수도 비율이 ε에 따라 어떻게 척도화되는지에 따라 ε → 0 극한에서 어떤 효과적 모델이 도출되는가?
• RQ2 밀도 및 투수도 척도화에서의 지수 κ 및 λ로 정의된 다양한 영역이 스케일업된 방정식의 구조에 어떻게 영향을 미치는가?
• RQ3 리처드스 방정식에 의해 지배되는 비퇴적 불포화 유동에 대해 ε-문제에서 효과적 모델로의 수렴을 엄밀히 증명할 수 있는가?
• RQ4 작은 ε에 대해 효과적 모델이 원래 ε-문제를 얼마나 정확하게 근사하는가? 특히 압력 헤드와 수분 함량 측면에서.
• RQ5 효과적 모델의 수치 해가 전체 2차원 균열 문제의 핵심 특성, 특히 유동 역학과 수분 함량 분포를 잘 반영하는가?

주요 결과

• 효과적 모델은 척도 지수 κ 및 λ에 따라 크게 달라진다: λ < 1 및 κ > -1일 경우 균열은 감소된 리처드스 방정식에 의해 지배되는 1차원 인터페이스로 행동한다.
• ε ≤ 10⁻³일 경우 균열 및 매트릭스 블록에서의 L² 오차는 10⁻⁵ 이하이며, 이는 비선형 솔버의 정밀도 허용 오차와 일치하여 높은 수치 정확도를 나타낸다.
• ε = 10⁻²일 경우에도 L² 오차는 10⁻⁴ 이하로 유지되어, ε → 0 근처에서 효과적 모델의 수렴성이 뚜렷하게 확인된다.
• ε-해에서 균열 두께에 따른 압력 헤드의 x-평균값이 효과적 모델 해와 매우 유사하게 일치하여, 평균화 가정이 타당함을 검증한다.
• 압력 연속성에도 불구하고, 유체 이동성의 저하로 인해 균열은 매트릭스 블록보다 수분 함량이 낮게 유지되며, 이는 모델이 예측한 그대로이다.
• 이론적 프레임워크는 실용적 시뮬레이션에서 이산 균열 모델의 사용을 지지하며, 그 물리적 일관성에 대한 엄밀한 수학적 기반을 제공한다.