[논문 리뷰] Travelling wave solutions for gravity fingering in porous media flows

Koondanibha Mitra, Andreas Rätz외 1명|arXiv (Cornell University)|2020. 01. 01.

Fluid Dynamics and Thin Films참고 문헌 29인용 수 1

한 줄 요약

이 논문은 동적 캡일러리 압력과 히스테리시스 모델을 사용하여 불포화 다공성 매체에서 중력 패드링에 대한 진행파 해를 연구한다. 시간에 따라 변하는 문제를 이동 기준 프레임으로 변환함으로써, 일정한 속도로 진행하는 단일 패드링에 대한 자유 경계 문제를 수립하고, 변분 방법을 통해 해의 존재성을 증명하며, 시간에 따라 변하는 시뮬레이션과 실험 결과와 일치하는 패드링의 형태와 속도를 보여주는 수치 결과를 제시한다.

ABSTRACT

We study an imbibition problem for porous media. When a wetted layer is above a dry medium, gravity leads to the propagation of the water downwards into the medium. In experiments, the occurrence of fingers was observed, a phenomenon that can be described with models that include hysteresis. In the present paper, we describe a single finger in a moving frame and set up a free boundary problem to describe the shape and the motion of one finger that propagates with a constant speed. We show the existence of solutions to the travelling wave problem and investigate the system numerically.

연구 동기 및 목표

• 기본 리처즈 방정식이 패드링 흐름 패atters를 포착하지 못하는 건조한 다공성 매체에서 중력 구동 흡수 현상을 모델링하기 위해.
• 표준 모델의 한계를 극복하기 위해, 특히 p = pc(s) + τ∂ts 형태의 τ 보정 모델을 포함한 히스테리시스와 동적 캡일러리 압력 효과를 통합하기 위해.
• 일정한 속도로 진행하는 단일 패드링의 형태와 속도를 기술하기 위해, 공-moving 프레임에서의 진행파 수식을 분석하기 위해.
• 주어진 유량과 경계 조건을 갖는 반무한 영역에서 발생하는 자유 경계 문제의 해 존재성을 확립하기 위해.
• 수치적으로 패드링 프로파일과 전파 속도를 계산하여 시간에 따라 변하는 시뮬레이션과 실험 관측 결과와의 일치 여부를 검증하기 위해.

제안 방법

• 시간에 따라 변하는 흡수 문제를 좌표 (y, z + ct)를 사용하여 진행파 프레임으로 변환함으로써, PDE 시스템을 이동 기준 프레임에서 정적 문제로 축소시킨다.
• 연결된 PDE 시스템으로 수식화: c∂zs = ∇⋅(k(s)[∇p + gez]) 및 cτ∂zs = [p − pc(s)]+ 로, 질량 보존과 동적 캡일러리 히스테리시스를 나타낸다.
• 에너지 기능에 k(s)와 압력 기울기를 포함하는 변분 원리를 사용하여, 절단된 영역 ΩH = (0,L)×(0,H)에서 해의 존재성을 증명한다.
• 경계 조건으로 z = 0에서 s = s0 및 p = p0 (하단), z → ∞에서 F∞ 유량, 측면에 무유량 조건을 적용한다.
• 수치적 연속성과 유한 요소 방법을 사용하여 절단된 문제를 해결하고, 패드링 프로파일과 전파 속도를 추출한다.
• sH와 ∇pH에 대한 리프시츠 유계를 확보하기 위해 부호 함수의 정규화를 사용하여, H → ∞의 극한에서 안정성과 수렴성을 보장한다.

실험 결과

연구 질문

• RQ1 히스테리시스와 동적 캡일러리 압력이 있는 다공성 매체에서 단일 중력 패드링에 대해 진행파 해를 구성할 수 있는가?
• RQ2 이러한 해의 존재 조건은 이동 기준 프레임에서 어떤가? 그리고 이는 유량, 속도, 경계 자료에 어떻게 의존하는가?
• RQ3 절단된 문제의 수치적 해는 전체 무한 도메인 문제를 얼마나 잘 근사하는가? H → ∞의 극한에서 압력과 포화도의 행동은 어떠한가?
• RQ4 계산된 패드링 형태와 전파 속도는 시간에 따라 변하는 시뮬레이션과 실험 결과와 일치하는가?
• RQ5 동적 캡일러리 매개변수 τ는 패드링 프로파일과 안정성에 어떤 역할을 하는가?

주요 결과

• 에너지 기능을 최소화하는 변분 수식을 통해 절단된 영역 ΩH에서 진행파 문제의 해가 존재함을 입증하였다.
• H → ∞의 극한에서 압력가 유계이기 위해서는 고정된 높이 z0에서 k(sH)의 적분이 0으로부터 멀리 떨어져 있어야 하며, 그렇지 않으면 압력가 무한대가 된다.
• 수치 시뮬레이션은 시간에 따라 변하는 시뮬레이션과 실험 관측 결과와 일치하는 패드링 형태의 포화도 프로파일을 생성하여 진행파 접근법의 타당성을 확인하였다.
• 전파 속도 c는 시간에 따라 변하는 계산에서 얻은 값과 양호한 일치를 보였지만, 다양한 수치적 방법 간에 약간의 변동이 있었다.
• 해는 공간과 시간에 대해 리프시츠 연속성을 보이며, |∂zsH|와 |∂ysH|의 유계는 c, τ 및 초기 자료 s0, p0에 따라 달라진다.
• 분석 결과, 해의 비유일성이 계산된 c 값의 변동성을 설명할 수 있으며, 동일한 경계 조건에서도 여러 가능한 패드링 속도가 존재할 수 있음을 시사한다.