Immiscible two-phase flow in porous media: a statistical mechanics approach

이 논문은 에드윈 제인스가 제안한 정보 이론적 엔트로피 기반의 통계역학 접근법을 활용하여, 다공성 매체 내 불혼화성 2 상 유동의 연속체 규모 현상을 미시적 물리학과 직접적으로 연결하면서도 관리 가능한 복잡도로 설명하는 새로운 열역학적 형식주의를 제시합니다.

핵심

이 논문은 불투과성 두 유체 (예: 기름과 물) 가 스펀지 같은 다공성 매체 ( porous media) 를 통과할 때 일어나는 복잡한 현상을 설명하는 새로운 이론을 제시합니다.

기존의 공학적 접근법이 너무 단순하거나 복잡해서 실용적이지 못했던 문제를, **'통계역학 (Statistical Mechanics)'**이라는 물리학의 강력한 도구를 빌려 해결하려는 시도입니다.

이 복잡한 내용을 일상적인 언어와 비유로 쉽게 풀어서 설명해 드리겠습니다.

---

1. 문제의 핵심: "스펀지 속의 혼란스러운 교통체증"

상상해 보세요. 거대한 스펀지 (다공성 매체) 가 있습니다. 이 스펀지 안에는 아주 작은 구멍들이 무수히 많이 뚫려 있습니다. 이제 이 스펀지를 통해 **물 (젖은 유체)**과 **기름 (젖지 않은 유체)**을 동시에 흘려보낸다고 칩시다.

• 기존의 문제점: 과학자들은 오랫동안 이 현상을 설명하기 위해 '상대 투과도 (Relative Permeability)'라는 공식을 썼습니다. 하지만 이는 마치 "차량이 많으면 속도가 느려진다"는 경험칙만 있을 뿐, 왜 그런지, 어떻게 상호작용하는지에 대한 깊은 이치 (미시적 물리) 를 설명하지 못했습니다. 마치 교통 체증을 설명할 때 "차량이 많으니까 막힌다"라고만 하고, 왜 차들이 서로 끼어들고, 어디로 흐르는지 설명하지 못하는 것과 비슷합니다.
• 새로운 접근: 이 논문은 이 현상을 **'통계역학'**이라는 렌즈로 다시 봅니다. 통계역학은 보통 원자나 분자처럼 아주 작은 입자들이 모여 거시적인 현상 (온도, 압력 등) 을 만들어내는 원리를 설명합니다. 저자들은 "스펀지 속의 유체 흐름도 마치 입자들이 모여 거시적인 흐름을 만드는 것과 같다"고 가정하고, 이를 설명하는 새로운 수학적 틀을 만들었습니다.

2. 새로운 도구: "정보 이론"과 "엔트로피"

이론의 핵심은 에드윈 제인스 (Edwin Jaynes) 가 제안한 **'정보 이론적 엔트로피'**를 사용하는 것입니다.

• 비유: 스펀지 속의 유체 배치를 알지 못할 때, 우리는 "모든 배치가 일어날 확률이 같다"고 가정합니다. 이때 우리가 '모르는 것'의 양을 수치화한 것이 엔트로피입니다.
• 적용: 저자들은 스펀지 속의 유체가 어떻게 흐르는지 (미시적 상태) 를 알지 못하는 상황에서, 우리가 알고 있는 것들 (전체 유량, 각 유체의 점유 면적 등) 을 제약 조건으로 두고, **가장 무질서한 상태 (최대 엔트로피)**를 찾아내어 거시적인 흐름 법칙을 유도했습니다.

3. 발견된 새로운 개념들: "유체 열역학"

이 과정을 통해 일반 열역학 (온도, 압력, 화학적 퍼텐셜) 과 매우 유사한 새로운 '유체 열역학'이 탄생했습니다. 여기서 등장하는 재미있는 개념들은 다음과 같습니다.

A. '아기처 (Agiture)' - 유체의 '흥분도'

• 일반 열역학: 온도는 분자들이 얼마나 활발하게 움직이는지를 나타냅니다.
• 유체 열역학: 저자들은 **'아기처 (Agiture)'**라는 새로운 변수를 만들었습니다. 이는 유체가 얼마나 '흥분'되어 있는지, 즉 압력 구배 (압력 차이) 가 얼마나 큰지를 나타내는 척도입니다.
• 비유: 스펀지 속의 물과 기름이 미친 듯이 뛰어다니는 상태라면 '아기처'가 높은 것입니다. 압력을 강하게 가하면 유체들이 더 활발해지고, 이는 마치 온도가 오르는 것과 같은 효과를 줍니다.

B. '공이동 속도 (Co-moving Velocity)' - 함께 움직이는 속도

이 논문에서 가장 중요한 발견 중 하나입니다.

• 현상: 물과 기름이 스펀지 속을 흐를 때, 두 유체는 서로를 밀어내거나 끌어당기며 흐릅니다. 단순히 각 유체의 속도를 더하는 것만으로는 설명이 안 됩니다.
• 비유: 혼잡한 지하철을 생각해 보세요.
  • 일반 속도: 각 승객이 걷는 속도.
  • 공이동 속도: 승객들이 서로 부딪히거나, 옆사람을 밀어내며 함께 움직일 때 생기는 '공동의 흐름' 속도입니다.
  • 이 논문은 두 유체가 서로를 밀어내며 흐를 때, 마치 **한 덩어리의 유체처럼 움직이는 '공유된 속도'**가 존재한다고 말합니다. 이 속도를 알면 유체의 흐름을 훨씬 정확하게 예측할 수 있습니다.

C. '유체 유리 (Glassy Flow)' - 얼어붙은 상태

• 현상: 유체의 흐름 속도가 매우 느리면, 물과 기름의 경계면이 스펀지 구멍에 걸려 움직이지 않는 상태가 됩니다.
• 비유: **유리 (Glass)**처럼 딱딱하게 굳어버린 상태입니다. 액체처럼 보이지만 움직이지 않는 '유동성 없는' 상태입니다.
• 의미: 이 상태에서는 유체의 흐름이 예측 불가능하게 변하고 (히스테리시스 현상), 마치 유리가 깨지기 직전의 불안정한 상태와 비슷합니다. 저자들은 이 현상을 '유리 전이 (Glass Transition)'라고 부르며, 이것이 왜 유체 흐름이 갑자기 막히거나 터지는지 설명해 줍니다.

4. 결론: 왜 이것이 중요한가?

이 논문은 단순히 수식을 바꾼 것이 아닙니다.

1. 단순화: 복잡한 스펀지 속 유체 흐름을, 마치 '온도'와 '압력'을 다루듯 체계적인 '열역학'으로 정리했습니다.
2. 새로운 변수: '공이동 속도'와 같은 새로운 변수를 도입하여, 기존에는 설명하지 못했던 유체 간의 복잡한 상호작용을 설명할 수 있게 되었습니다.
3. 실용성: 석유 회수, 지하수 관리, 이산화탄소 저장 등 실제 공학 분야에서 유체 흐름을 더 정확하게 예측하고 제어할 수 있는 길을 열었습니다.

한 줄 요약:

"스펀지 속의 물과 기름 흐름을 설명할 때, 단순한 공식을 버리고 **'유체들의 집단 심리 (통계역학)'**를 분석하여, 마치 온도와 압력처럼 새로운 **'유체 열역학'**을 만들어낸 혁신적인 연구입니다."

이 연구는 복잡한 자연 현상을 이해할 때, 우리가 가진 직관적인 물리 법칙 (열역학) 을 확장하여 적용하면 놀라운 통찰을 얻을 수 있음을 보여줍니다.

기술 요약

논문 요약: 불혼화성 2 상 유동의 통계역학적 접근

1. 연구 배경 및 문제 제기 (Problem)

• 핵심 문제: 다공성 매질 내 불혼화성 2 상 유동 (예: 물과 기름) 의 물리학에서 가장 중요한 과제는 스케일 업 (Scale-up) 문제입니다. 즉, 미세한 기공 (pore) 스케일의 물리 현상을 다공성 매질을 연속체로 간주할 수 있는 대규모 (Darcy 스케일) 로 어떻게 적절히 기술할 것인가 하는 문제입니다.
• 기존 접근법의 한계: 현재까지 실용적으로 사용되고 있는 유일한 방법은 상대 투과도 (Relative Permeability) 이론입니다. 이는 경험적 방정식 (Darcy 법칙의 일반화) 에 기반하며, 모세관 압력 ($p_c$) 과 상대 투과도 ($k_{rw}, k_{rn}$) 가 포화도 ($S_w$) 의 함수라고 가정합니다.
  • 약점: 이 이론은 현상론적 (phenomenological) 이며, 기공 스케일의 물리학과 명확한 연결고리가 부족합니다. 또한, 복잡도가 급증하여 기존 이론을 넘어서는 시도는 실용적인 적용에 실패해 왔습니다.
• 목표: 기공 스케일 물리학과 Darcy 스케일 연속체 역학을 연결하면서도 복잡도를 관리 가능한 수준으로 유지하는 새로운 이론적 틀을 구축하는 것입니다.

2. 방법론 (Methodology)

• 제인스 (Jaynes) 의 일반화된 통계역학 적용:
  • 저자들은 1950 년대 Edwin T. Jaynes 가 제안한 정보 이론적 엔트로피 (Shannon entropy) 기반의 통계역학 일반화 기법을 차용했습니다.
  • 기존의 열역학이 분자 스케일에서 거시적 열역학으로 스케일 업되는 것과 유사하게, 이 연구는 기공 스케일의 유체역학을 Darcy 스케일의 '열역학과 유사한 형식 (Thermodynamics-like formalism)'으로 스케일 업합니다.
• 대표 기본 면적 (REA, Representative Elementary Area) 설정:
  • 원통형 다공성 시료 내에서 유동 방향에 수직인 단면을 정의하고, 통계적 특성을 대표할 수 있으면서도 유동 변동을 포함하는 REA 를 설정합니다.
  • REA 를 통과하는 체적 유량 ($Q_p$), 습윤 유체 유량 ($Q_w$), 비습윤 유체 유량 ($Q_n$) 및 각 유체가 차지하는 면적 ($A_w, A_n$) 을 확률 변수로 취급합니다.
• 최대 엔트로피 원리 적용:
  • 알려진 평균값 (유량, 면적 등) 에 대한 제약 조건 하에서 **구성 엔트로피 (Configurational Entropy, $S$)**를 최대화하여 확률 분포 $p(X)$를 유도합니다.
  • 이를 통해 분배 함수 (Partition function) 와 새로운 열역학적 변수들을 도출합니다.

3. 주요 기여 및 이론적 결과 (Key Contributions & Results)

가. 비열적 열역학 (Non-thermal Thermodynamics) 의 정립

• 기존 열역학의 변수들을 대체하는 새로운 **집적 변수 (Intensive variables)**가 등장합니다:
  1. Agiture ($\theta$): '교란 온도 (Agitation temperature)'로 불리며, 압력 구배 ($|\nabla p|$) 에 비례합니다. 이는 시스템의 운동 에너지를 나타냅니다.
  2. Flow Derivative ($\mu$): 습윤 면적 ($A_w$) 또는 포화도 ($S_w$) 의 켤레 변수로, 화학 퍼텐셜과 유사한 역할을 합니다.
  3. Flow Pressure ($\pi$): 기공 면적 ($A_p$) 의 켤레 변수로, 기존 열역학에는 존재하지 않는 새로운 변수입니다.
• 이 형식을 통해 Darcy 스케일에서 열역학의 모든 관계식 (예: Legendre 변환, 오일러 정리 등) 을 적용할 수 있게 됩니다.

나. 공동 이동 속도 (Co-moving Velocity, $v_m$) 의 발견

• 열역학적 속도 ($\hat{v}_w, \hat{v}_n$) 와 실제 유동 속도 ($v_w, v_n$) 의 분리:
  • 열역학적 속도는 엔트로피와 면적의 변화율로 정의되지만, 실험적으로 측정되는 실제 속도와는 다릅니다.
  • 두 속도 사이의 관계를 연결하는 핵심 변수로 **공동 이동 속도 ($v_m$)**가 등장합니다.
  • 관계식: $v_w = \hat{v}_w - S_n v_m$, $v_n = \hat{v}_n + S_w v_m$.
• $v_m$의 특성:
  • $v_m$은 포화도 변화에 따른 평균 유속 변화율과 실제 유속의 차이에서 기인합니다.
  • 실험 데이터 및 시뮬레이션을 통해 $v_m$이 $\mu$ (Flow Derivative) 에 대해 선형적으로 변하는 간단한 함수 형태를 가진다는 것이 확인되었습니다.
  • 이는 기존 상대 투과도 이론이 간과해 왔던 중요한 물리량으로, 2 상 유동의 비선형 거동을 설명하는 열쇠가 됩니다.

다. 위상도 및 유리 전이 (Phase Diagram & Glass Transition)

• 유리상 (Glassy Phase) 의 존재:
  • Darcy 스케일에서 유동은 모세관 수 (Capillary number, $Ca$) 에 따라 3 가지 위상 (I, II, III) 으로 나뉩니다.
  • 특히 낮은 $Ca$ 영역에서 **유리 전이 (Glass transition)**가 발생함이 발견되었습니다. 이는 인터페이스의 운동이 동적으로 멈추거나 (Kinetic arrest) 매우 느려지는 상태입니다.
  • Edwards-Anderson 질서 매개변수와 스핀 글래스 감수성 (Spin-glass susceptibility) 을 통해 이 전이가 열역학적 위상 전이와 유사한 임계선 (Critical line) 으로 확인되었습니다.
• 히스테리시스 설명: 유리상에서의 거동은 다공성 매질 유동에서 관찰되는 히스테리시스 현상과 큰 변동을 설명할 수 있는 새로운 틀을 제공합니다.

라. 일반 열역학과의 유사성 (Co-molar Volume)

• 이 연구에서 도출된 '공동 이동 속도 ($v_m$)' 개념은 일반 열역학의 **공 몰 부피 (Co-molar volume, $\nu_m$)**와 수학적으로 동일한 구조를 가집니다.
• 이는 불혼화성 2 상 유동이라는 비열적 시스템에서 도출된 개념이 일반 열역학 (액체 혼합물 등) 으로도 확장 가능함을 시사하며, 통계역학적 접근의 보편성을 보여줍니다.

4. 의의 및 결론 (Significance & Conclusion)

• 이론적 통합: 다공성 매질 내 불혼화성 2 상 유동이 단순히 공학적 현상이 아니라, 통계역학적 구조를 가진 비열적 시스템임을 규명했습니다.
• 실용적 함의:
  • 기존 상대 투과도 이론의 한계를 극복하고, 기공 스케일 물리와 대규모 거동을 연결하는 체계적인 이론적 틀을 제시했습니다.
  • 특히 **공동 이동 속도 ($v_m$)**는 새로운 관측 가능 변수로, 유동 모델의 정확도를 높이는 데 필수적입니다.
• 미래 전망:
  • 현재 연구는 정상 상태 (Steady-state) 유동에 국한되어 있으나, 이를 비정상 상태 (Non-equilibrium) 로 확장하고, 모세관 압력 함수 ($p_c$) 없이 Agiture ($\theta$) 와 Flow Derivative ($\mu$) 의 구배를 구동력으로 사용하는 비평형 열역학으로 발전시킬 수 있습니다.
  • 유리상 유동 (Glassy flow) 의 특성을 더 깊이 연구함으로써 실제 석유 회수, 지하수 오염 관리 등 공학적 응용 분야에 기여할 수 있습니다.

이 논문은 다공성 매질 유동 물리학에 통계역학의 강력한 도구를 도입하여, 복잡한 유동 현상을 체계적이고 예측 가능한 열역학적 형식으로 재해석한 획기적인 연구입니다.