다상 흐름의 정상상태 유효 투과성 스케일링 연구

초록

본 논문은 2차원 모형 다공성 매체에서 공기와 물‑글리세린 용액의 두 상 흐름을 실험과 네트워크 시뮬레이션으로 조사한다. 정상상태에서 유효 투과성은 캡릴러리 수에 대한 거듭 제곱법칙을 따르며, 점도 대비가 큰 경우와 점도가 일치하는 경우 모두에서 동일한 경향이 관찰된다. 정체 클러스터가 존재하는 상태에서는 지수 0.500.54, 클러스터가 없을 때는 포화도에 따라 0.670.80 사이의 값이 나타난다.

상세 분석

이 연구는 두 상 흐름이 다공성 매체 내부에서 어떻게 상호작용하며 전체 투과성을 결정하는지를 정량적으로 규명하고자 한다. 실험에서는 2차원 유리 마이크로채널에 공기와 물‑글리세린 혼합액을 주입해 정상상태에 도달하도록 하였으며, 캡릴러리 수(Ca)와 유효 투과성(k_eff)의 관계를 측정하였다. 실험 결과는 k_eff ∝ Ca^γ 형태의 거듭 제곱법칙을 보였으며, γ는 약 0.54(점도 대비 큰 경우)와 0.50(점도 일치 경우)로 수렴하였다. 이를 뒷받침하기 위해 저자들은 2차원 정규 격자 네트워크 모델을 구축하고, 각 결절에 랜덤한 모공 크기와 접촉각을 부여해 실제 매체의 이질성을 모사하였다. 네트워크 시뮬레이션에서는 두 가지 전이 상태가 확인되었다. 첫 번째는 ‘정체 클러스터(stagnant cluster)’ 상태로, 일부 포어가 흐름에 거의 참여하지 못하고 고정된 채 남아 있다. 이 경우 흐름 경로가 제한되어 전체 투과성이 감소하고, γ는 0.50~0.54 범위에 머문다. 두 번째는 ‘클러스터 비존재’ 상태로, 모든 포어가 일정 수준 이상의 흐름을 유지한다. 여기서는 포화도(S)와 Ca에 따라 γ가 변동하며, 0.67에서 0.80 사이의 값을 보였다. 특히 포화도가 중간값에 가까울 때 γ가 최대에 도달한다는 점이 주목할 만하다. 점도 대비가 큰 경우와 점도 일치 경우 모두에서 동일한 스케일링 거동이 나타난 것은, 캡릴러리 힘이 점도 차이에 비해 지배적인 역할을 함을 시사한다. 또한, 정체 클러스터가 형성되는 메커니즘은 낮은 Ca에서 표면 장력에 의해 포어 간 압력 차가 크게 발생하고, 이로 인해 일부 포어가 ‘잠금’되는 현상으로 해석된다. 이러한 현상은 기존의 비정상 상태(전이 현상)와는 구별되며, 정상상태에서도 복잡한 구조적 비균일성이 투과성에 큰 영향을 미친다. 논문은 실험과 시뮬레이션 결과를 정량적으로 비교함으로써, 네트워크 모델이 실제 다공성 매체의 흐름 특성을 충분히 재현할 수 있음을 입증한다. 최종적으로, 캡릴러리 수에 대한 유효 투과성의 스케일링 법칙은 다상 흐름의 설계와 예측에 중요한 기준이 될 수 있다.