레일리-테일러 불안정

유체역학 전산유체역학 마지막 수정: 2026-07-20 04:22:44

1. 개요[편집]

레일리-테일러 불안정 Rayleigh–Taylor Instability
약칭	RTI
발생 조건	무거운 유체가 가벼운 유체 위에 놓임
구동력	중력(또는 가속도) × 밀도 차
특징 형상	버섯 구름형 손가락(finger)·스파이크
사촌 불안정	켈빈-헬름홀츠 불안정

무거운 물을 가벼운 공기 위에 얹으려는 순간, 자연이 그 배치를 거부한다.

레일리-테일러 불안정(Rayleigh–Taylor instability, RTI)은 밀도가 다른 두 유체가 겹쳐 있을 때, 무거운 유체가 가벼운 유체 위에 놓이면 경계면이 불안정해지며 격렬하게 뒤섞이는 현상이다. 물이 든 컵을 재빨리 뒤집으면 물이 쏟아지는 것도, 미시적으로 보면 무거운 물과 가벼운 공기의 경계가 이 불안정을 일으켜 무너지는 것이다. 자연은 무거운 것을 아래에 두고 싶어 하며, 그 반대 배치는 아주 작은 교란만 있어도 스스로 붕괴한다.

RTI는 유체 불안정의 교과서적 대표 사례로, 켈빈-헬름홀츠 불안정과 함께 자주 언급된다. 둘의 차이는 구동력에 있다 — 켈빈-헬름홀츠는 두 층의 속도 차(전단)가, 레일리-테일러는 밀도 차와 중력(가속도)이 불안정을 만든다. 초신성 폭발, 관성 핵융합, 대기·해양 대류, 그리고 다상유동 시뮬레이션의 벤치마크까지, RTI는 스케일을 가리지 않고 등장하는 우주적 단골손님이다.¹

2. 발생 조건과 물리[편집]

두 정지 유체가 수평 경계면을 맞대고 있다고 하자. 위쪽 밀도 $\rho_1$ , 아래쪽 밀도 $\rho_2$ 일 때, 불안정 여부는 단순하다.

\rho_1 > \rho_2 \quad \Rightarrow \quad \text{불안정}

무거운 쪽이 위에 있으면 계의 위치에너지가 배치를 뒤집을수록 낮아지므로, 자연은 그 방향으로 굴러떨어진다. 경계면에 미세한 물결이 생기면 무거운 유체는 아래로 파고들며 스파이크(spike) 를, 가벼운 유체는 위로 솟으며 버블(bubble) 을 형성한다. 이 상호 침투가 진행되면 특유의 버섯 구름 모양 손가락 구조가 자란다. 중력 $g$ 대신 유체가 겪는 임의의 가속도가 무거운 유체에서 가벼운 유체 쪽을 향하기만 하면 동일한 불안정이 발생하므로, 관성 핵융합 캡슐처럼 급가속되는 계에서도 나타난다.

3. 성장률과 아트우드 수[편집]

선형 단계에서 진폭이 작을 때, 파수 $k$ 인 교란은 지수적으로 성장한다. 표면장력과 점성을 무시한 이상적 경우 성장률 $n$ 은

n = \sqrt{A \, g \, k}, \qquad A = \frac{\rho_2 - \rho_1'}{\rho_2 + \rho_1'}

로 주어진다. 여기서 $A$ 가 아트우드 수(Atwood number)로, 두 유체의 밀도 대비를 나타내는 무차원수다.² $A$ 가 1에 가까울수록(예: 물 위의 공기) 성장이 빠르고, 0에 가까우면(비슷한 밀도) 완만하다. 성장률이 $\sqrt{k}$ 에 비례한다는 것은 짧은 파장일수록 빠르게 자란다는 뜻인데, 이 성질이 수치 시뮬레이션에서 골칫거리가 된다. 격자를 촘촘히 할수록 더 짧은 파장이 해상되어 더 빠른 불안정 모드가 깨어나므로, 초기 교란을 제어하지 않으면 격자 의존적인 결과가 나온다.

현실에서는 이 무한 성장을 막는 안정화 요인이 둘 있다. 표면장력은 짧은 파장의 곡률 에너지를 키워 매우 짧은 파장을 억제하고, 점성은 전 파장에 걸쳐 성장을 늦춘다. 그 결과 특정 파장에서 성장률이 최대가 되는 “가장 위험한 모드”가 존재한다.

4. 발달 단계[편집]

RTI는 시간에 따라 뚜렷한 국면을 거친다.

선형 단계: 진폭이 파장보다 작을 때. 각 파장 모드가 독립적으로 지수 성장한다. 이론이 잘 맞는 유일한 구간.
약비선형 단계: 스파이크와 버블이 뚜렷해지고 형상이 비대칭해진다. 무거운 스파이크는 좁고 날카롭게, 가벼운 버블은 넓고 둥글게 자란다.
비선형·혼돈 단계: 스파이크 끝에서 켈빈-헬름홀츠 불안정이 2차로 발생해 버섯 갓처럼 말려 올라가고, 여러 파장 모드가 병합하며 큰 구조로 성장한다.
난류 혼합 단계: 경계면이 완전히 붕괴해 두 유체가 뒤섞인 혼합층이 형성된다. 혼합층 두께가 $\propto A g t^2$ 로 자란다는 것이 실험·시뮬레이션의 공통 결론이다.

5. 수치 시뮬레이션[편집]

RTI는 전산유체역학에서 계면 추적(interface capturing) 기법의 성능을 재는 표준 벤치마크다. 두 유체의 경계면이 접혔다 말리며 극도로 복잡해지므로, 계면을 얼마나 날카롭게 유지하는지가 관건이다. 주로 쓰이는 도구는 다음과 같다.

VOF 방법: 체적 분율로 계면을 추적. 질량 보존이 강점.
레벨셋 방법: 부호 거리 함수로 계면을 매끄럽게 표현. 곡률 계산이 정확.

수치해석에서 RTI의 함정은 앞서 말한 격자 의존성이다. 물리적 초기 교란을 명시하지 않으면, 격자 크기·이산화 오차·수치 소산과 분산이 만들어내는 인공 교란이 씨앗 역할을 해 격자마다 다른 손가락 무늬가 나온다. 그래서 RTI 검증(검증 및 확인) 문제는 초기 섭동의 파장·진폭을 명확히 규정하는 것이 국룰이다.³

6. 관련 문서[편집]

RTI가 가장 극적으로 보이는 곳은 우주다. 게자리 성운(Crab Nebula)의 필라멘트 구조, 초신성 잔해의 뒤엉킨 손가락들이 모두 폭발로 밀려난 무거운 물질이 가벼운 물질을 뚫고 나온 RTI의 흔적이다. 컵의 물이 쏟아지는 것과 별이 터지는 것이 같은 방정식을 따른다는 게 유체역학의 낭만. ↩
조지 아트우드(George Atwood)의 이름을 딴 무차원수지만, 정작 그는 18세기 물리학자로 RTI와는 무관하다. 도르래 실험(아트우드 기계)으로 유명한 사람의 이름이 밀도비에 붙은 것은, 두 질량의 비를 다루는 형태가 닮아서다. ↩
“격자를 더 촘촘히 했더니 결과가 더 지저분해졌다”는 RTI 초심자의 흔한 당혹감이다. 이건 수렴 실패가 아니라 물리다 — 짧은 파장이 더 빨리 자라니 고해상 격자가 더 많은 손가락을 깨운다. 진짜 수렴을 보려면 통계량(혼합층 두께)을 봐야지, 손가락 무늬 하나하나를 맞추려 들면 안 된다. ↩