켈빈-헬름홀츠 불안정

유체역학 전산유체역학 물리 마지막 수정: 2026-07-19 04:37:52

1. 개요[편집]

켈빈-헬름홀츠 불안정 Kelvin–Helmholtz Instability
약칭	KH 불안정
명명	Lord Kelvin, Hermann von Helmholtz (1868~71)
원인	경계면을 가로지르는 속도 전단
결과 구조	KH 롤업(고양이 눈 와류)
안정화 요인	표면장력, 안정 성층, 점성

켈빈-헬름홀츠 불안정(Kelvin–Helmholtz instability)은 유체 내부에 속도 전단이 존재하는 경계면에서 미소 교란이 지수적으로 증폭되어, 경계면이 파도처럼 말려 올라가 소용돌이 열로 발전하는 유동 불안정 현상이다. 바람이 수면을 스치면 파도가 서고, 담배 연기 기둥이 어느 높이에서 갑자기 흐트러지고, 목성 대적점 가장자리에 톱니 무늬가 생기는 것 — 전부 같은 놈이다.

시뮬레이션 하는 입장에서 KH 불안정이 중요한 이유는 두 가지다. 첫째, 층류에서 난류로 넘어가는 천이의 대표적 경로라서 전단층·혼합층·제트가 있는 모든 해석에 개입한다. 둘째, 수치 기법의 실력을 가장 잔인하게 드러내는 벤치마크다. 격자와 차분 도식이 조금만 부실해도 없어야 할 소용돌이가 생기거나, 있어야 할 소용돌이가 사라진다.¹

2. 물리적 메커니즘[편집]

경계면이 살짝 물결치면 볼록한 마루 쪽은 유선이 좁아져 유속이 빨라지고, 베르누이 방정식에 따라 그 지점의 압력이 낮아진다. 낮아진 압력은 마루를 더 끌어올린다. 골 쪽은 정반대로 더 눌린다. 교란이 스스로를 키우는 양의 되먹임이 성립하는 것이고, 이것이 불안정의 정체다.

와도 관점에서 보면 더 깔끔하다. 전단층은 와도가 얇게 뭉쳐 있는 와류층(vortex sheet)이며, 이 층은 균일 상태를 유지하지 못하고 일정 간격으로 뭉쳐(clumping) 이산 와류들로 재배열된다. 뭉친 와류들이 서로를 유도해 회전시키면서 생기는 것이 그 유명한 고양이 눈(cat’s eye) 모양의 KH 롤업이다.

3. 선형 안정성 해석[편집]

두 비점성·비압축성 층(밀도 $\rho_1$ , $\rho_2$ , 수평 속도 $U_1$ , $U_2$ )이 맞닿아 있고 중력 $g$ 와 표면장력 $\sigma$ 가 작용할 때, 파수 $k$ 의 교란에 대한 불안정 조건은 다음과 같다.

(U_1 - U_2)^2 > \frac{\rho_1 + \rho_2}{\rho_1 \rho_2} \left[ \frac{g(\rho_2 - \rho_1)}{k} + \sigma k \right]

( $\rho_2$ 가 아래쪽 무거운 유체) 우변을 $k$ 에 대해 최소화하면 임계 속도차가 나오는데, 이것이 켈빈이 계산한 바람에 의한 파도 발생 임계 풍속이다. 공기-물 물성으로 넣으면 대략 6~7 m/s 수준의 값이 나온다.²

우변의 두 항이 안정화 요인이라는 점을 눈여겨볼 만하다.

중력 항( $g/k$ ) — 무거운 유체가 아래에 있으면 긴 파장(작은 $k$ )을 억제한다. 성층이 반대로 뒤집히면 이 항이 부호를 바꿔 레일리-테일러 불안정이 된다.
표면장력 항( $\sigma k$ ) — 짧은 파장(큰 $k$ )을 억제한다.

즉 중간 파장대가 가장 먼저 터진다. 실제 관측되는 KH 파의 파장이 특정 대역에 몰리는 이유다.

4. 무한 성장률과 문제의 성질[편집]

중력도 표면장력도 없는 순수 와류층에서는 성장률이

\gamma = \frac{\sqrt{\rho_1 \rho_2}}{\rho_1 + \rho_2} \, k \, |U_1 - U_2|

로 $k$ 에 정비례한다. 파장이 짧을수록 무한히 빨리 자란다는 뜻이고, 이는 문제 자체가 비정상적(ill-posed)이라는 신호다. 현실에서는 점성이나 표면장력, 혹은 전단층의 유한한 두께가 짧은 파장을 잘라내면서 문제를 정상화한다.

수치해석적으로도 이건 남 일이 아니다. 격자가 무한 파수를 잘라주는 사실상의 정규화 장치로 작동하기 때문에, 격자를 바꾸면 물리적 답이 아니라 격자가 정한 답이 바뀐다. 격자 수렴을 확인하지 않은 KH 계산은 아무것도 증명하지 못한다. 검증 및 확인이 유난히 목소리를 높이는 영역.

5. 안정 성층과 리처드슨 수[편집]

대기·해양처럼 밀도가 연속적으로 성층된 유동에서는 국소 리처드슨 수로 판단한다.

\mathrm{Ri} = \frac{N^2}{(\mathrm{d}U/\mathrm{d}z)^2}, \qquad N^2 = -\frac{g}{\rho}\frac{\mathrm{d}\rho}{\mathrm{d}z}

$N$ 은 부력 진동수(Brunt–Väisälä 진동수)다. 마일스–하워드 정리에 따르면 유동장 전역에서 $\mathrm{Ri} > 1/4$ 이면 안정이 보장된다. 다만 이건 충분조건일 뿐이라, $\mathrm{Ri} < 1/4$ 라고 반드시 불안정한 것은 아니다. 이 지점을 헷갈려서 논문 리뷰에서 지적당하는 일이 잊을 만하면 반복된다.

대기 중 KH 파는 권층운이 규칙적인 파도 무늬로 나타나는 이른바 “빌로우 클라우드”로 눈에 보이고, 청천난류(CAT)의 주요 원인으로도 지목된다. 비행기가 구름 한 점 없는 하늘에서 갑자기 덜컹거리면 대개 이 친구 소행이다.

6. 시뮬레이션에서[편집]

전단이 있는 곳이면 어디든 나온다. 유동박리 후 형성되는 전단층, 카르만 와열의 초기 와류 형성, 제트 혼합층, 다상유동의 계면 파괴와 액적 생성(VOF 방법으로 계면을 추적하는 그 문제), 관성 밀폐 핵융합의 캡슐 혼합(하이드로코드 영역), 초신성 잔해와 은하 가스 흐름까지.

수치적으로 특히 악명 높은 지점들.

수치 점성이 롤업을 뭉갠다. 저차 상류 도식은 KH 소용돌이를 그냥 없애버린다. 차분 도식의 차수와 소산 특성을 반드시 확인해야 한다.
역으로, 미해상(under-resolved) 계산은 가짜 소용돌이를 만든다. 격자 스케일에서 발생한 수치 잡음이 KH 메커니즘에 의해 증폭되어, 물리처럼 보이는 구조로 자라난다. 그림은 그럴듯한데 격자를 두 배로 늘리면 완전히 다른 그림이 나온다.³
레이놀즈수가 명시되지 않은 KH 계산은 사실상 무의미하다. 점성이 없으면 소산 스케일이 없고, 결국 격자 크기가 콜모고로프 스케일 역할을 대신하게 되므로 결과가 격자에 종속된다.
그래서 KH는 스펙트럴 방법이나 고차 불연속 갤러킨 계열을 검증하는 단골 시험 문제이며, 대와류 모사나 직접수치모사 코드의 통과의례이기도 하다.

7. 선형 성장 이후[편집]

선형 안정성 해석이 답해주는 것은 딱 초기 단계, 교란이 아직 작을 때까지다. 그다음부터는 비선형 영역이고, 대략 다음 순서를 밟는다.

롤업 — 와도가 뭉치며 고양이 눈 와류가 형성된다. 여기까지는 2차원적이고 규칙적이다.
와류 병합(pairing) — 이웃한 와류 둘이 서로를 감고 돌다 하나로 합쳐진다. 병합이 일어날 때마다 특성 파장이 두 배가 되고, 이것이 혼합층이 하류로 갈수록 두꺼워지는 주된 성장 메커니즘이다.
3차원 2차 불안정 — 주 와류를 잇는 축방향 갈비뼈 모양의 와류(rib vortex)가 발달한다. 여기서부터 유동은 본격적으로 3차원이 된다.
난류 천이 — 스케일이 잘게 쪼개지며 콜모고로프 스케일까지 에너지가 흘러 내려간다.

2차원 계산에서는 3단계 이후가 원천적으로 불가능하다. 2차원에서는 와도 신장(vortex stretching)이 존재하지 않아 에너지가 역캐스케이드로 큰 스케일로 올라가 버리기 때문이다. 예쁜 고양이 눈 그림이 2차원 계산으로 넘쳐나는 것과 달리, 그 그림들이 실제 난류 혼합을 말해주지는 않는다는 뜻. 혼합률이나 연행(entrainment)을 논하려면 3차원이 필수다.

8. 사촌 격의 불안정들[편집]

레일리-테일러 불안정 — 무거운 유체가 가벼운 유체 위에 얹힌 경우. KH와 함께 나타나는 일이 많다.
리히트마이어-메시코프 불안정 — 충격파가 밀도 경계면을 통과할 때. RT의 충격파 버전.
홈스-슐리팅 파 — 벽 경계층 내부의 점성 천이 경로. KH가 비점성 메커니즘인 것과 대비된다.

9. 관련 문서[편집]

그래서 새 솔버를 만들면 다들 KH부터 돌려본다. 예쁜 고양이 눈이 나오면 자랑스럽게 스크린샷을 올리는데, 정작 격자 수렴 결과를 같이 올리는 사람은 드물다. ↩
켈빈의 계산은 물이 정지해 있고 공기가 균일하게 부는 이상화 모델이라, 실제 잔물결은 그보다 훨씬 약한 바람에서도 생긴다. 실제 파도 생성에는 마일스의 임계층 공명 등 다른 메커니즘이 함께 작동한다. 고전 결과의 값은 어디까지나 기준점으로 받아들이는 게 맞다. ↩
이 문제를 정면으로 파고든 연구들이 밝힌 바에 따르면, 계면을 매끄럽게 처리하지 않은 KH 초기조건은 애초에 격자 독립적인 답을 가질 수 없다. 그림이 예뻐도 물리가 아닐 수 있다는 뜻. 무지개 컨투어의 배신은 어디에나 있다. ↩