오어-조머펠트 방정식

1. 개요[편집]

오어-조머펠트 방정식(Orr-Sommerfeld equation)은 평행 전단류(parallel shear flow)에 가해진 미소 교란의 시간적 성장 여부를 지배하는 4계 상미분방정식으로, 점성 유동의 선형 안정성 이론(linear stability theory)의 핵심 방정식이다. 쉽게 말하면 “이 얌전해 보이는 층류가 언제 참지 못하고 난류로 폭발하는가”를 계산해주는 방정식이다. 층류→난류 천이 문제의 출발점이자, 유체역학에서 가장 유명한 고유값 문제 중 하나.

나비에-스토크스 방정식을 기저 유동 $U(y)$ 위에 얹힌 미소 교란에 대해 선형화하고, 교란을 파동 형태로 가정해 정리하면 튀어나온다. 이름은 윌리엄 맥패든 오어(William McFadden Orr, 1907)와 아르놀트 조머펠트(Arnold Sommerfeld, 1908)가 각자 독립적으로 유도한 데서 왔다.¹

2. 유도와 형태[편집]

기저 유동을 $x$ 방향으로만 흐르는 평행류 $\mathbf{U} = (U(y), 0, 0)$로 잡는다. 여기에 미소 교란 유동함수(stream function)를 정규모드(normal mode)로 가정한다.

$$\psi(x, y, t) = \phi(y)\, e^{i\alpha(x - ct)}$$

여기서 $\alpha$는 실수 파수(wavenumber), $c = c_r + i c_i$는 복소 위상속도다. 핵심은 $c_i$의 부호. $c_i > 0$이면 교란이 시간에 따라 $e^{\alpha c_i t}$로 지수 성장하여 불안정, $c_i < 0$이면 감쇠하여 안정, $c_i = 0$이면 중립(neutral)이다. 이 가정을 선형화된 방정식에 넣으면 진폭함수 $\phi(y)$에 대한 오어-조머펠트 방정식이 나온다.

$$(U - c)(\phi'' - \alpha^2 \phi) - U'' \phi = \frac{1}{i\alpha Re}(\phi'''' - 2\alpha^2 \phi'' + \alpha^4 \phi)$$

여기서 $Re$는 레이놀즈수다. 우변이 점성항인데, $Re \to \infty$로 보내 우변을 통째로 날리면 2계 방정식인 레일리 방정식(Rayleigh equation)이 된다. 즉 오어-조머펠트는 레일리 방정식에 점성이라는 4계 살을 붙인 버전인 셈. 재밌는 건 점성이 보통 안정화 요인일 것 같지만, 실제로는 점성이 오히려 불안정을 유발하는 경우(Tollmien-Schlichting 파)가 있다는 점이다.²

3. 고유값 문제와 임계 레이놀즈수[편집]

경계에서 교란이 사라진다는 조건($\phi = \phi' = 0$)을 벽면에 부여하면, 오어-조머펠트 방정식은 주어진 $(\alpha, Re)$에 대해 특정한 $c$에서만 자명하지 않은 해를 갖는 고유값 문제가 된다. 즉 $c = c(\alpha, Re)$. 이 관계를 분산 관계(dispersion relation)라 부른다.

$c_i = 0$이 되는 경계면을 $\alpha$–$Re$ 평면에 그리면 중립 안정 곡선(neutral stability curve)이 얻어지는데, 이 곡선이 감싸는 영역 안쪽이 불안정 영역이다. 이 곡선의 코, 즉 불안정이 시작되는 가장 낮은 레이놀즈수를 임계 레이놀즈수(critical Reynolds number) $Re_c$라 한다.

• 평면 푸아죄유 유동(두 평판 사이 압력구배 유동): $Re_c \approx 5772$ (Orszag, 1971)³
• 블라시우스 경계층: $Re_c \approx 520$ (변위두께 기준)
• 평면 쿠에트 유동(두 평판 상대운동): $Re_c = \infty$ — 선형적으로는 영원히 안정

마지막 항목이 유체역학의 오랜 골칫거리다. 쿠에트 유동은 오어-조머펠트로는 어떤 $Re$에서도 불안정하지 않은데, 실험에서는 $Re \sim 350$ 근처에서 멀쩡히 난류로 천이한다. 파이프 유동(하겐-푸아죄유)도 마찬가지로 선형 안정한데 실제로는 $Re \sim 2000$에서 난류가 된다.

4. 아임계 천이의 역설[편집]

이 불일치의 원인은 오어-조머펠트가 다루는 것이 어디까지나 선형·점근적(asymptotic) 성장이라는 데 있다. 실제 유동은 유한한 교란에 대해, 그리고 유한한 시간 동안 교란이 대수적으로 폭증하는 비정규(non-normal) 성장을 겪는다. 오어-조머펠트 연산자는 자기수반(self-adjoint)이 아니어서 고유함수들이 직교하지 않고, 이 비직교성 때문에 각 고유모드는 감쇠해도 그 합은 일시적으로 수백 배까지 커질 수 있다(transient growth). 이 순간적 증폭이 비선형 항을 깨워 천이를 촉발한다는 것이 1990년대 이후 정립된 관점이다.⁴ 그래서 요즘은 고유값 대신 최적 교란(optimal perturbation)과 유사스펙트럼(pseudospectrum)을 본다.

위 시뮬레이션은 원기둥 뒤 전단층 이야기이지만, 전단이 있는 유동이 임계 조건을 넘으면 교란이 스스로 자라 규칙적인 와류 구조로 발전한다는 오어-조머펠트의 직관을 눈으로 확인할 수 있다. 카르만 와열 역시 원기둥 후류 전단층의 불안정이 만드는 대표적 구조다.

5. 수치 해법[편집]

오어-조머펠트 방정식은 손으로는 거의 못 푼다. 표준 접근은 스펙트럴 방법이다. 진폭함수 $\phi(y)$를 체비쇼프 다항식으로 전개하고 체비쇼프-가우스-로바토 격자에서 미분 연산자를 이산화하면, 방정식이 일반화 고유값 문제 $A\mathbf{v} = c B\mathbf{v}$로 바뀐다. 이걸 QZ 알고리즘 같은 표준 고유값 솔버로 때리면 스펙트럼 전체가 나온다.

Orszag가 1971년에 이 체비쇼프-타우(tau) 방법으로 $Re_c = 5772.22$를 6자리까지 계산해낸 것이 기념비적 사례. 다만 고레이놀즈수에서는 행렬이 심하게 병조건(ill-conditioned)이 되어 가짜 고유값(spurious eigenvalue)이 스펙트럼에 섞여 나온다. 이걸 걸러내는 것이 실무의 관건이며, 유사스펙트럼 분석이 특히 예민하다. 관련해서 발생하는 수치 잡음은 절단오차와 부동소수점 정밀도 문제와도 얽힌다.

6. 관련 문서[편집]

• 나비에-스토크스 방정식
• 난류
• 카르만 와열
• 스펙트럴 방법
• 고유값
• 절단오차
• 무차원수
• 전산유체역학

7. Footnotes[편집]

1. 두 사람이 1년 간격으로 따로 유도했는데도 이름이 사이좋게 둘 다 붙었다. 요즘 같으면 우선권 분쟁으로 학회지 댓글창이 불탔겠지만, 20세기 초의 신사들은 그러지 않았다. 참고로 조머펠트는 노벨상을 네 번이나 놓친 것으로도 유명한 이론물리학자다. ↩

2. 점성이 불안정을 만든다니 직관에 반하는데, 벽 근처에서 점성이 위상을 살짝 어긋나게 만들어 교란이 기저 유동에서 에너지를 뽑아내도록 돕기 때문이다. “안정화의 대명사가 배신자였다”는 반전. ↩

3. Orszag, S. A. (1971). Accurate solution of the Orr-Sommerfeld stability equation. J. Fluid Mech. 이 숫자를 소수점 둘째 자리까지 외우고 다니면 유체 안정성 전공자로 인정받는다. 반대로 5772를 5722로 잘못 적으면 리뷰어가 바로 알아챈다. ↩

4. Trefethen et al. (1993), Science에 실린 “Hydrodynamic stability without eigenvalues”가 이 관점 전환의 상징. 제목부터 “고유값 없는 안정성 이론”이라며 대놓고 도발적이다. 결론은 오어-조머펠트가 틀렸다기보다, 우리가 고유값만 봤던 게 순진했다는 것. ↩