다공성 매질 유동

전산유체역학 유체역학 시뮬레이션 마지막 수정: 2026-07-16 04:26:41

1. 개요[편집]

다공성 매질 유동 Flow in Porous Media
기본 법칙	다르시 법칙 (Darcy, 1856)
핵심 물성	투과율 $K$, 다공도 $\varepsilon$
고속 보정	포체이머 항 (관성)
확장 방정식	브링크만 방정식
주 사용처	지하수, 석유, 연료전지, 열교환기

스펀지 속 구멍을 다 격자로 그릴 것인가, 아니면 스펀지를 그냥 “저항 있는 공기”로 칠 것인가. 후자를 택한 사람들이 만든 학문.

다공성 매질 유동은 고체 골격 사이의 미세한 공극(pore)을 통과하는 유체의 흐름을, 공극 하나하나를 해상하지 않고 평균화된 연속체로 기술하는 유동 해석 분야다. 모래층, 촉매층, 라디에이터 핀, 연료전지 확산층, 필터 — 안이 뻥 뚫린 게 아니라 복잡하게 막혀 있는 모든 것이 대상이다.

핵심 트릭은 스케일 분리다. 공극 규모의 실제 유동은 나비에-스토크스 방정식을 따르지만, 그걸 다 풀려면 자갈 하나 주위 격자를 수백만 개 깔아야 한다. 대신 여러 공극이 통계적으로 대표되는 부피에 대해 방정식을 평균 내면, 복잡한 벽면 마찰 전체가 투과율이라는 스칼라 하나로 압축된다. 미시적 지옥을 거시적 계수 하나로 갈아 넣는 것 — 다공성 매질 모델의 전부이자 한계다.

2. 다르시 법칙[편집]

1856년 프랑스 디종의 엔지니어 앙리 다르시(Henry Darcy)는 도시 급수 필터를 설계하다가, 모래층을 통과하는 물의 유량이 압력 강하에 정확히 비례한다는 것을 실험으로 찾아냈다.¹ 미분 형태로 쓰면 이렇다.

\mathbf{u}_D = -\frac{K}{\mu} \nabla p

$\mathbf{u}_D$ 는 다르시 속도(겉보기 속도, superficial velocity), $K$ 는 투과율 $[\mathrm{m}^2]$ , $\mu$ 는 점성계수다. 형태를 보면 알겠지만 이건 열전도의 푸리에 법칙, 확산의 픽 법칙과 완전히 같은 골격이다. 비압축 연속 방정식과 결합하면 압력에 대한 라플라스 방정식이 튀어나온다.

\nabla \cdot \left( \frac{K}{\mu} \nabla p \right) = 0

즉 균질 매질의 다르시 유동은 그냥 퍼텐셜 문제다. 비선형 대류항이 사라지므로 포텐셜 유동이나 정상 열전도와 수학적으로 형제 사이. 발산이라는 단어를 모르고 살 수 있는 몇 안 되는 CFD 영역이다.

여기서 짚고 넘어갈 두 가지 속도가 있다. 다르시 속도는 매질 전체 단면으로 나눈 값이고, 공극 안 실제 유속(간극 속도, interstitial velocity)은 이보다 빠르다. 둘은 두피-포크하이머(Dupuit–Forchheimer) 관계로 이어진다.

\mathbf{u}_D = \varepsilon \, \mathbf{u}_{\text{pore}}

이 $\varepsilon$ 을 헷갈려서 유속을 다공도만큼 틀리게 보고하는 것이 이 분야 신입의 통과의례다.

3. 투과율과 다공도[편집]

다공도(porosity) $\varepsilon$ 은 전체 부피 중 빈 공간의 비율이다. 무작위로 쌓은 구는 대략 0.36~0.40, 금속 폼은 0.9 이상, 잘 다져진 점토는 0.3 근처.

투과율 $K$ 는 이름과 달리 유체 물성이 아니라 기하 구조만의 물성이다. 같은 모래층이면 물이 흐르든 기름이 흐르든 $K$ 는 같다. 점성 영향은 이미 $\mu$ 로 따로 빠져 있기 때문. 단위가 면적( $\mathrm{m}^2$ )인 것이 힌트인데, 실무에서는 다르시(darcy, $\approx 0.987 \times 10^{-12}\,\mathrm{m}^2$ )라는 단위를 쓴다. 석유공학 논문에서 밀리다르시(mD)가 튀어나오면 이 얘기다.

구 충전층의 투과율은 코제니-카르만(Kozeny–Carman) 식으로 추정한다.

K = \frac{\varepsilon^3 d_p^2}{180 (1-\varepsilon)^2}

$\varepsilon$ 의 3제곱이 눈에 띌 것이다. 다공도가 0.4에서 0.3으로 떨어지면 투과율은 4배 넘게 급감한다. 압축된 지층이 왜 그렇게 안 뚫리는지, 필터가 왜 갑자기 막히는지가 이 지수 하나로 설명된다.

이방성 매질(층리가 있는 지층, 직조 섬유)에서는 $K$ 가 스칼라가 아니라 2계 텐서다. 수평 투과율이 수직의 10배인 셰일층 같은 것.

4. 포체이머 보정[편집]

다르시 법칙은 느린 유동에서만 맞다. 유속이 올라가면 공극 안에서 유체가 굽이치며 관성 효과가 생기고, 압력 강하가 속도의 제곱에 비례하는 성분을 갖는다. 이를 보정한 것이 포체이머(Forchheimer) 방정식이다.

-\nabla p = \frac{\mu}{K}\mathbf{u}_D + \frac{C_F \rho}{\sqrt{K}} |\mathbf{u}_D| \mathbf{u}_D

앞이 점성(다르시) 항, 뒤가 관성(포체이머) 항. 판별 기준은 공극 규모 레이놀즈수로, 대략 $Re_p \gtrsim 1{\sim}10$ 부터 다르시가 깨지기 시작한다.² 여기서 흔한 오해 하나 — 관성항이 생겼다고 난류가 된 게 아니다. 이 영역은 여전히 층류이고, 진짜 공극 난류는 $Re_p$ 수백 이상에서 온다.

충전층에 대해서는 두 항을 한꺼번에 준 에르군(Ergun) 식이 사실상의 표준이다.

\frac{\Delta p}{L} = 150 \frac{(1-\varepsilon)^2}{\varepsilon^3} \frac{\mu u_D}{d_p^2} + 1.75 \frac{1-\varepsilon}{\varepsilon^3} \frac{\rho u_D^2}{d_p}

150과 1.75라는 마법의 상수는 이론이 아니라 실험 회귀에서 나온 값이다. 1952년 논문의 숫자 두 개가 아직도 전 세계 화학공장의 압력 강하를 계산하고 있다.

5. 브링크만 방정식[편집]

다르시 법칙에는 치명적 결함이 있다. 1계 방정식이라 점착 조건(no-slip)을 못 준다. 벽에서 속도가 0이 되어야 하는데 다르시 해는 그냥 벽까지 유한한 속도로 흘러간다. 자유 유동 영역과 다공 영역이 맞닿는 문제(강물 옆 대수층, 연료전지 채널과 확산층)에서는 이게 그냥 틀린 답이 된다.

브링크만(Brinkman) 방정식은 여기에 점성 확산항을 되살린다.

-\nabla p + \tilde{\mu} \nabla^2 \mathbf{u}_D - \frac{\mu}{K}\mathbf{u}_D = 0

$\tilde{\mu}$ 는 유효 점성계수. $K \to \infty$ 면 스토크스 유동으로, $K \to 0$ 이면 다르시로 매끄럽게 넘어간다. 하나의 방정식으로 자유 유동과 다공 유동을 동시에 푸는 통합 프레임을 주기 때문에, COMSOL Multiphysics 같은 코드가 좋아한다. 다만 브링크만 층의 두께는 $\sqrt{K}$ 스케일이라 보통 마이크론 단위 — 이걸 격자로 해상하려면 결국 격자를 지독하게 깔아야 한다.

6. 대표체적요소[편집]

다공성 매질 모델링 전체를 떠받치는 개념이 대표체적요소(Representative Elementary Volume, REV)다. 다공도를 재려고 측정 부피를 키워 나가면, 처음에는 값이 요동친다(공극 하나에 걸치느냐 골격에 걸치느냐에 따라 0이나 1). 부피가 커지면 값이 안정되는 평탄 구간이 나타나고, 더 키우면 매질 자체의 대규모 불균질성 때문에 다시 변한다. 그 중간 평탄 구간이 REV다.

\ell_{\text{pore}} \ll \ell_{\text{REV}} \ll L_{\text{macro}}

이 부등식이 성립할 때만 “연속체로서의 다공성 매질”이라는 개념이 성립한다. 깨지는 순간 — 채널 폭이 입자 3개짜리인 미세유체, 균열 하나가 전체 유량을 지배하는 파쇄암 — 다르시 법칙은 그냥 무의미해지고 공극 규모 직접 모사(격자 볼츠만 방법이 자주 쓰인다)나 이산요소법 기반 CFD-DEM으로 내려가야 한다. REV를 확인하지 않고 다르시를 쓰는 것은 검증 및 확인을 건너뛰는 것과 같다.³

7. CFD에서의 다공성 존[편집]

상용 코드에서 다공성 매질을 다루는 방식은 놀랍도록 실용적이다. 격자는 그냥 빈 공간처럼 깔고, 해당 셀 영역(porous zone)의 운동량 방정식에 저항 소스항을 더한다.

S_i = -\left( \frac{\mu}{K_i} u_i + C_{2,i} \frac{1}{2}\rho |u| u_i \right)

사용자가 입력하는 건 방향별 점성 저항 계수 $1/K_i$ 와 관성 저항 계수 $C_2$ 딱 두 세트. ANSYS Fluent, STAR-CCM+, OpenFOAM의 explicitPorositySource가 전부 이 구조다.

실무 적용 예:

열교환기 / 라디에이터: 핀 수천 장을 격자로 그리는 대신 다공성 블록 하나로. 자동차 공력해석에서 라디에이터를 이렇게 처리하는 게 표준이다.
촉매 컨버터 / 필터: 벌집 채널을 방향성 저항(축방향 $K$ 는 크게, 횡방향은 0에 가깝게)으로.
팬·다공판·메시 스크린: 압력 강하 실험 곡선을 그대로 회귀해서 $1/K$ , $C_2$ 로 환산.

계수를 얻는 정석은 실험 $\Delta p$ – $u$ 곡선을 $\Delta p/L = a u + b u^2$ 로 피팅해 $a = \mu/K$ , $b = C_2 \rho / 2$ 로 역산하는 것이다. 실험 데이터가 없으면? 벤더 카탈로그를 뒤지거나, 대표 셀 하나만 상세 격자로 풀어서 곡선을 뽑아낸다.

주의할 점: 다공성 존 안의 속도는 다르시 속도(겉보기 속도)로 보고되는 것이 기본값이라, 온도장이나 체류 시간을 볼 때 다공도 보정을 안 하면 조용히 틀린다. 그리고 다공성 존 안에서는 난류 모델링이 물리적 의미를 거의 잃는다 — 셀 크기가 공극보다 훨씬 크므로 그 안의 와류를 대변할 방법이 없다. 그래서 다공 영역의 난류 생성을 그냥 끄는 옵션이 존재한다.

8. 활용 분야[편집]

지하수 / 수문학: 대수층 모델링, 오염물 이송. MODFLOW가 이 바닥의 왕이다.
석유·가스: 저류층 시뮬레이션. 다상(물-기름-가스) 다르시에 상대투과율 곡선을 얹는다. 여기서는 다상유동과 완전히 결합한다.
연료전지 / 전지: 기체확산층(GDL)과 촉매층 안의 반응물 수송. 다공성 유동 + 전기화학 + 열이 얽힌 다물리 연성해석의 교과서적 사례.
CO2 지중 저장: 초임계 CO2의 부력 상승과 트래핑. 요즘 이 분야 논문의 절반.
생체: 뼈 속 골수 유동, 종양 조직의 약물 침투.

9. 관련 문서[편집]

논문 제목이 Les Fontaines Publiques de la Ville de Dijon(디종 시의 공공 분수, 1856)이다. 상하수도 설계 보고서 부록에 실린 실험 하나가 저류층 공학과 수문학이라는 학문 두 개를 낳았다. 요즘 기준으로는 “부록 D”에서 노벨상이 나온 셈. ↩
공극 레이놀즈수는 정의부터가 문헌마다 다르다. 특성 길이를 입자 지름으로 잡느냐 $\sqrt{K}$ 로 잡느냐, 속도를 다르시 속도로 잡느냐 간극 속도로 잡느냐에 따라 값이 10배씩 왔다 갔다 한다. “다르시 한계는 $Re \approx 1$ “이라는 문장을 볼 때마다 어느 $Re$ 인지 확인해야 하는 이유. ↩
파쇄암 대수층에서 다르시로 오염물 이동을 예측했다가, 실제로는 균열 하나를 타고 예측보다 훨씬 빨리 도달한 사례가 지하수 오염 소송사에 여럿 있다. 평균화는 평균이 의미 있을 때만 의미가 있다. ↩