새프먼-테일러 불안정

편집 역사 토론
유체역학 전산유체역학 시뮬레이션 마지막 수정: 2026-08-12 04:38:19

1. 개요[편집]

새프먼-테일러 불안정
Saffman–Taylor instability
다른 이름점성 손가락 (viscous fingering)
제시P. G. 새프먼 · G. I. 테일러, 1958
무대헬레-쇼 셀 / 다공성 매질
발생 조건저점성 유체가 고점성 유체를 밂
지배 법칙다르시 법칙 + 라플라스 압력 도약
안정화 요인표면장력 (짧은 파장 차단)
유명한 결과채널에서 손가락 폭 → 채널폭의 1/2

물로 기름을 밀어내면, 물은 기름을 밀지 않는다. 기름을 뚫고 지나간다.

새프먼-테일러 불안정은 점성이 낮은 유체가 점성이 높은 유체를 밀어낼 때 두 유체의 경계면이 불안정해져, 저점성 유체가 가늘고 긴 손가락 모양으로 침투하는 현상이다. 필립 새프먼과 제프리 잉그램 테일러가 1958년 논문에서 헬레-쇼 셀과 다공성 매질을 함께 다루며 이론과 실험을 동시에 제시했고, 이후 “점성 손가락”이라는 이름으로 비평형 패턴 형성 연구의 표준 문제가 됐다.

원리는 한 줄이다. 유체가 밀려 나가려면 자기 앞의 유체를 밀어야 하는데, 앞에 있는 놈이 점성이 커서 잘 안 밀린다. 경계면 어딘가가 우연히 조금 앞서 나가면 그쪽은 남은 고점성 구간이 짧아져 저항이 줄고, 따라서 더 빨리 나간다. 앞선 놈이 더 빨라지는 양의 되먹임 — 불안정의 교과서적 형태다. 레일리-테일러 불안정에서 중력이 하던 역할을 여기서는 압력 구배와 이동도 차이가 대신한다.

공학적 무대는 지하다. 석유 회수에서 물이나 가스를 주입해 원유를 밀어내는 순간 정확히 이 불안정이 발생하고, 손가락이 생산정까지 뚫고 나오면(breakthrough) 그 뒤로는 주입한 물만 나온다. 채굴 효율이 이 하나로 결정되므로 다공성 매질 유동 시뮬레이션에서는 가장 비싼 물리 현상이다.1

2. 헬레-쇼 셀은 왜 2차원 다르시인가[편집]

**헬레-쇼 셀**은 간격 bb 가 아주 좁은 평행 유리판 두 장이다. 간격 방향으로 유동이 완전 발달한 포아죄유 흐름이라 가정하고 두께 방향으로 평균 내면, 면내 평균 속도가 압력 구배에 정비례한다.

u=−b212μ ∇p\mathbf{u} = -\frac{b^{2}}{12\mu}\,\nabla p

형태가 다르시 법칙 u=−(K/μ)∇p\mathbf u = -(K/\mu)\nabla p 와 완전히 같고, 투과율이 K=b2/12K = b^2/12 로 기하만으로 결정된다. 다공성 매질의 미로 같은 유동을 유리판 두 장으로 눈에 보이게 재현할 수 있다는 것이 헬레-쇼 셀이 100년 넘게 사랑받는 이유다.

비압축 조건 ∇⋅u=0\nabla\cdot\mathbf u = 0 을 대입하면 각 유체 영역에서 압력이 라플라스 방정식을 만족한다. 그래서 이 문제는 자유경계를 가진 라플라스 문제가 된다 — 계면의 모양이 압력장을 정하고, 압력장이 계면을 움직이고, 다시 압력장이 바뀐다. 확산 제한 응집이나 유전파괴와 수학적 골격이 같은 “라플라스 성장” 부류의 대표 주자다.

경계에서의 두 조건은 이렇다. 계면을 가로질러 법선 속도가 연속이고, 압력은 표면장력 γ\gamma 와 계면 곡률 κ\kappa 에 의해 도약한다(라플라스 압력).

p2−p1=γ κp_{2} - p_{1} = \gamma\,\kappa

여기서 짚어야 할 근사가 하나 있다. 헬레-쇼 계면의 곡률에는 면내 곡률 말고 간격 방향 곡률 ∼2/b\sim 2/b 도 있는데, 표준 유도에서는 이것을 상수로 보고 무시한다. 실제로는 젖음막(wetting film)이 남으면서 이 성분이 속도에 의존하고(브레서턴형 Ca2/3Ca^{2/3} 보정), 이 사실이 뒤에 나올 손가락 폭 문제의 실험 편차를 상당 부분 설명한다.

3. 선형 안정성과 분산관계[편집]

평평한 계면이 xx 방향으로 속도 UU 로 전진한다고 하고, 뒤쪽(미는 쪽) 점성을 μ1\mu_1, 앞쪽(밀리는 쪽) 점성을 μ2\mu_2 라 하자. 계면 교란을 η=η^ eiky+σt\eta = \hat\eta\,e^{iky + \sigma t} 로 놓고 각 영역에서 라플라스 방정식을 풀어 두 조건을 맞추면 다음이 나온다.

σ(k)  =  kμ1+μ2[(μ2−μ1) U  −  b2γ12 k2]\sigma(k) \;=\; \frac{k}{\mu_1 + \mu_2}\left[(\mu_2 - \mu_1)\,U \;-\; \frac{b^{2}\gamma}{12}\,k^{2}\right]

읽는 법은 이렇다.

  • 불안정 판정. 대괄호 안 첫 항의 부호가 전부다. μ2>μ1\mu_2 > \mu_1, 즉 점성이 큰 쪽을 작은 쪽이 밀 때만 불안정하다. 반대로 밀면 모든 파수가 감쇠한다. 실무에서는 이동도비 M=μ2/μ1M = \mu_2/\mu_1 로 쓰고 M>1M>1 이면 불리한 배치(unfavorable displacement)라 부른다.
  • 표면장력이 짧은 파장을 죽인다. k3k^3 항 때문에 차단 파수 kc=12(μ2−μ1)U/(b2γ)k_c = \sqrt{12(\mu_2-\mu_1)U/(b^2\gamma)} 위로는 전부 안정하다. 표면장력이 없으면 σ∝k\sigma \propto k 가 되어 짧은 파장일수록 무한히 빨리 자란다 — 문제 자체가 부적정(ill-posed)이 된다.
  • 최불안정 파장. dσ/dk=0d\sigma/dk = 0 에서 km=kc/3k_m = k_c/\sqrt3 이고, 파장으로 쓰면
λm  =  2πkm  =  πb γ(μ2−μ1) U  =  πbCa\lambda_m \;=\; \frac{2\pi}{k_m} \;=\; \pi b\,\sqrt{\frac{\gamma}{(\mu_2-\mu_1)\,U}} \;=\; \frac{\pi b}{\sqrt{Ca}}

마지막 등호는 모세관 수 Ca=(μ2−μ1)U/γCa = (\mu_2-\mu_1)U/\gamma 로 묶은 것이다. 결론이 깔끔하다 — 손가락 굵기는 셀 간격에 비례하고 속도의 제곱근에 반비례한다. 천천히 밀면 굵은 손가락 몇 개, 빨리 밀면 가는 손가락 여러 개. 실험에서 눈으로 확인되는 가장 유명한 스케일링이다.

다공성 매질에서도 같은 형태가 성립하는데, 공극 규모의 모세관 효과를 거시적으로 뭉뚱그린 유효 계면장력을 넣어야 한다(추크 등, 1959). 이 유효값은 물성이 아니라 모형 상수라, 저수지 시뮬레이션에서 손가락 폭이 격자 크기에 끌려다니는 근본 원인이 된다.

4. 손가락 폭 1/2 문제[편집]

여기서부터가 이 문제가 40년 가까이 사람을 괴롭힌 이유다.

새프먼과 테일러는 표면장력을 0으로 두고 폭 WW 인 채널을 정상 속도로 전진하는 손가락의 정확해를 등각사상으로 구했다. 손가락이 채널의 비율 λ\lambda 를 차지한다고 하면 손가락 끝을 원점에 둔 계면 모양은

x  =  W(1−λ)2π ln⁡ ⁣[12(1+cos⁡2πyλW)]x \;=\; \frac{W(1-\lambda)}{2\pi}\,\ln\!\left[\frac{1}{2}\left(1 + \cos\frac{2\pi y}{\lambda W}\right)\right]

이고, 질량 보존에서 손가락 전진 속도는 먼 상류 평균 속도의 1/λ1/\lambda 배다. 문제는 이 해가 0<λ<10 < \lambda < 1 인 모든 λ\lambda 에 대해 존재한다는 것이다. 방정식이 답을 하나로 고르지 못한다.

그런데 실험은 아주 분명한 답을 준다. 속도를 올릴수록 λ\lambda 가 0.5로 수렴한다. 채널 폭의 정확히 절반. 자연이 어떻게 이 값을 고르는가 — 이것이 선택 문제(selection problem)다.

답은 무시했던 표면장력에 있었다. 표면장력을 넣으면 방정식의 최고차 미분항이 바뀌므로 이것은 특이섭동이다. 맥린과 새프먼(1981)이 무차원 표면장력 매개변수

B  =  b2γ12 μ U W2B \;=\; \frac{b^{2}\gamma}{12\,\mu\,U\,W^{2}}

를 넣고 적분방정식을 수치적으로 풀어, 연속체였던 해의 족이 이산적인 가지들로 쪼개지고 그중 물리적으로 안정한 가지가 B→0B \to 0 에서 λ→1/2\lambda \to 1/2 로 간다는 것을 보였다. 즉 표면장력은 아무리 작아도 답을 고르는 일을 한다.

왜 정칙 섭동 전개로는 안 보이는가. 콩베스코·아킴·동브르·포모·퓌미르, 홍·랑거, 슈라이만이 1986년에 독립적으로 밝힌 바에 따르면, 선택 조건은 BB 의 모든 차수 전개에서 0인 지수적으로 작은 항 ∼e−c/B\sim e^{-c/\sqrt{B}} 의 균형에서 나온다. 이른바 모든 차수를 넘어선 점근(asymptotics beyond all orders) 또는 가해성 이론(solvability theory)이다. 테일러 급수를 아무리 오래 계산해도 절대 보이지 않는 항이 답을 정한다는 것 — 수치해석 하는 사람 입장에서는 소름 끼치는 결론이다.2

같은 구조가 수지상 성장의 가지 끝 반지름 선택에도 그대로 나타난다. 이방성이나 표면장력 같은 “작은” 항이 대역적 패턴 규모를 결정한다는 것이 이 계열 문제의 공통 교훈이고, 위상장 모형이 왜 계면 물리를 정확히 넣어야만 쓸모가 있는지의 이유이기도 하다.

5. 레일리-테일러와 무엇이 같고 무엇이 다른가[편집]

두 불안정은 자주 형제로 묶이는데, 정확히 어디서 같은지 짚어 두면 헷갈리지 않는다.

같은 것. 둘 다 계면을 가로지르는 어떤 “구동 대비”가 양수면 불안정하고 표면장력이 짧은 파장을 잘라 최불안정 파장을 만든다. 실제로 수직 헬레-쇼 셀에서 중력까지 넣으면 대괄호 안의 구동항이

(μ2−μ1)U⏟점성 대비  +  b212(ρ2−ρ1) g∥⏟밀도 대비\underbrace{(\mu_2-\mu_1)U}_{\text{점성 대비}} \;+\; \underbrace{\frac{b^{2}}{12}(\rho_2-\rho_1)\,g_{\parallel}}_{\text{밀도 대비}}

형태로 나란히 더해진다(g∥g_\parallel 는 전진 방향 중력 성분). 점성차 구동과 밀도차 구동이 문자 그대로 같은 괄호 안에 앉는다는 것 — 이것이 두 불안정이 구조적으로 같은 문제라는 말의 정확한 뜻이다.

다른 것. 성장률의 파수 의존성이 다르다.

계지배 동역학성장률 (표면장력 무시)
고전 레일리-테일러관성 (비점성)σ∼Agk\sigma \sim \sqrt{A g k}
새프먼-테일러점성 (다르시)σ∼k\sigma \sim k

레일리-테일러는 관성이 지배해 성장률이 k\sqrt{k} 로 자라고, 새프먼-테일러는 점성 저항이 지배해 kk 에 정비례한다. k\sqrt{k} 냐 kk 냐가 관성계와 다르시계를 가르는 지문이다. 헬레-쇼 셀에 무거운 유체를 위에 얹으면 레일리-테일러가 일어나되 성장률은 kk 에 비례한다 — 중력 구동인데도 다르시 지문을 갖는 것이다. 켈빈-헬름홀츠 불안정이 전단 구동이라는 제3의 축을 차지한다.

6. 손가락에서 프랙탈로[편집]

채널이 아니라 중심에서 방사형으로 주입하면 이야기가 달라진다. 손가락 끝이 자라면서 국소 곡률이 작아지고, 어느 순간 표면장력이 끝을 안정시키지 못해 끝이 둘로 갈라진다(tip splitting). 이것이 반복되면 프랙탈에 가까운 가지 구조가 나오고, 측정된 프랙탈 차원이 확산 제한 응집 군집과 비슷하다는 것이 패터슨(1984)에 의해 지적됐다. 압력장이 라플라스 방정식을 만족하고 계면 성장률이 −∇p-\nabla p 에 비례한다는 점에서 두 문제는 같은 방정식이며, 유전파괴 모형이 이 계열을 하나의 지수 η\eta 로 묶는다.

반대로 이방성을 인위적으로 넣으면 — 유리판에 격자무늬 홈을 파거나, 손가락 끝에 작은 기포를 붙여 국소적으로 성질을 바꾸면 — 끝 갈라짐 대신 규칙적인 수지상 가지가 나온다. 잡음이 지배하면 DLA, 이방성이 지배하면 덴드라이트라는 이 대비는 결정 성장 이론 전체를 관통하는 구도이며, 헬레-쇼 셀이 그 실험실 판본이다.

액정, 점탄성 유체, 겔 등으로 밀리는 쪽을 바꾸면 손가락 대신 균열 같은 갈라짐이 나오기도 한다. 손가락의 모양은 결국 밀리는 물질의 구성 방정식이 정한다.

7. 수치해석에서[편집]

  • 표면장력을 빼면 부적정 문제다. γ=0\gamma = 0 인 헬레-쇼 문제는 짧은 파장이 무한히 빨리 자라고, 정확해가 유한 시간에 첨점(cusp) 특이점을 만든다는 것이 알려져 있다. 그래서 격자를 조밀하게 할수록 결과가 더 지저분해지는 현상이 물리가 아니라 수치 노이즈일 수 있다. 수렴 검사는 손가락 무늬가 아니라 통계량(침투 면적, 손가락 개수 밀도, 돌파 시각)으로 해야 한다.
  • 경계적분법이 고전적 무기다. 압력이 각 영역에서 조화함수이므로 미지수를 계면 위 밀도로 줄일 수 있다. 계면만 이산화하니 정확도가 높지만, 끝 갈라짐 이후 위상이 바뀌면 다시 꿰매야 한다는 것이 한계.
  • 확산계면 접근이 대안이다. 위상장 모형이나 레벨셋 방법·VOF 방법을 쓰면 위상 변화가 저절로 처리된다. 대가는 계면 두께 ε\varepsilon 이 유효 표면장력에 섞여 들어간다는 것 — 선택 문제가 지수적으로 작은 항에 걸려 있는 계에서 계면 두께를 물리값보다 두껍게 잡으면 손가락 폭이 통째로 틀린다.
  • 저수지 시뮬레이션의 격자 의존성. 산업용 다공성 매질 코드는 손가락을 해상하지 않고 격자 셀 안에서 뭉갠다. 이 경우 실제 손가락 폭이 격자보다 작으면 회수율이 심각하게 과대평가되므로, 코츠-토드 혼합 모형 같은 하위격자 모형으로 보정한다. 격자를 못 줄이면 물리를 모형화해야 한다는 CFD의 오랜 교훈이 여기서도 반복된다.3
  • 격자 볼츠만 방법이나 공극 규모 직접 해석은 반대 방향의 접근이다. 공극 하나하나를 해상해 손가락의 씨앗을 물리적으로 만들지만, 저수지 규모로는 여전히 갈 수 없다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 공저자인 G. I. 테일러는 이미 레일리-테일러 불안정에 이름을 올린 그 테일러다. 계면 불안정 두 개에 이름이 박힌 사람이라 “테일러 불안정”이라고만 쓰면 어느 쪽인지 알 수 없다. 논문 제목에 다공성 매질과 헬레-쇼 셀이 나란히 적혀 있는 것도 상징적인데, 유리판 두 장이 곧 지층의 축소판이라는 선언이었다. ↩

  2. 우리가 수치해석에서 “작은 항은 무시해도 된다”고 배운 감각이 정면으로 부서지는 사례다. e−c/Be^{-c/\sqrt{B}} 는 B=0.01B=0.01 이면 10−510^{-5} 수준인데, 그 항이 없으면 답이 아예 정해지지 않는다. 존재하지 않는 것과 지수적으로 작은 것은 다르다. ↩

  3. 석유공학에서 회수율을 과대평가한 시뮬레이션은 곧 투자 결정 오류다. “손가락은 격자보다 가늘어서 안 보였습니다”라는 사후 보고가 실제로 존재하며, 그 손가락은 대체로 현장에서 먼저 발견된다. ↩