나비에-스토크스 방정식의 엄밀해

나비에-스토크스 방정식은 악명 높은 비선형 방정식입니다. 그런데 어떻게 정확하게 풀 수 있을까요?

답은 대칭성입니다. 유동의 기하 구조가 충분히 단순할 때 — 직선 파이프, 두 평판, 무한 평면 — 속도는 한 방향만을 가리키며 하나 또는 두 개의 좌표에만 의존하게 됩니다. 이러한 고전적 대칭 설정에서 비선형 항은 사라지거나 충분히 단순해져서 방정식이 선형 PDE로 축소됩니다 — 정상 상태에서는 종이와 연필로 풀 수 있는 ODE가 되는 경우도 많습니다.

이렇게 생각하면 됩니다: 나비에-스토크스 방정식은 모든 가능한 유체 운동을 기술합니다. 그러나 유체를 매우 질서 정연한 상황 — 소용돌이 없이, 혼돈 없이, 모든 것이 한 방향으로 진행하는 — 으로 몰아넣으면, 방정식의 복잡성 대부분이 무관해집니다. 나비에-스토크스의 어려운 부분은 유체가 자기 자신을 밀어내는 되먹임 고리입니다. 이러한 대칭 유동에서는 밀어낼 대상이 없어 비선형 자기 이류 항이 탈락합니다.

이러한 엄밀해는 단순한 호기심거리가 아닙니다. 유체역학 교육의 토대이자, 수치 코드의 벤치마크이며, 실제 유동이 이상화된 상태에서 얼마나 벗어나는지를 이해하기 위한 출발점입니다.

나비에-스토크스 방정식의 엄밀해란 다음을 만족하는 속도장 $u$와 압력 $p$입니다.

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0$$

수치적 근사 없이 닫힌 형태로 주어집니다.

핵심 메커니즘은 비선형 이류 항 $(u \cdot \nabla)u$의 소멸 또는 단순화입니다. 평행 유동 — 직교 좌표에서 $u = (u(y,t),\, 0,\, 0)$이거나 원통 좌표에서 $u = (0,\, 0,\, u(r,t))$인 유동 — 의 경우, 속도장은 자동으로 비압축성을 만족하고 이류 항은 항등적으로 사라집니다:

$$(u \cdot \nabla)u = u\,\partial_x u\, \hat{e}_x = 0$$

왜냐하면 $u$가 유선 방향 좌표 $x$에 의존하지 않기 때문입니다. 운동량 방정식은 선형 PDE로 축소됩니다 — 정상 유동의 경우 ODE가 됩니다.

유동이 정상($\partial_t u = 0$)이기도 하면, 남는 것은

$$\nu \Delta u = \nabla p,$$

이는 푸아송 방정식으로, 그 해가 점성 유동의 고전적 엄밀해인 포아즈유 유동, 쿠에트 유동 및 그 변형들입니다.

비정상 평행 유동($\partial_t u \neq 0$이지만 이류는 여전히 사라지는)의 경우 방정식은 확산 방정식 $\partial_t u = \nu\, \partial_{yy} u$가 되어, 스토크스 제1, 제2 문제를 줍니다.

포아즈유 유동(하겐-포아즈유 유동이라고도 함)은 나비에-스토크스 방정식의 가장 중요한 엄밀해이며, 대부분의 공학자가 처음 접하는 해입니다.

물이 길고 곧은 파이프를 통해 꾸준히 흐르는 모습을 상상해 보십시오. 양쪽 끝 사이의 일정한 압력 차이가 흐름을 구동합니다. 파이프 벽은 정지해 있으므로, 벽에 닿는 유체는 속도 0에 고정됩니다(미끄러짐 없는 조건). 벽에서 멀리 떨어진 파이프 중심부의 유체가 가장 빠르게 움직입니다.

결과적인 속도 프로파일은 포물선입니다: 벽에서 0, 중심에서 최대, 그 사이에 매끄러운 곡선을 그립니다. 파이프를 잘라 단면을 보면, 속도 분포가 뒤집힌 그릇 모양으로 보일 것입니다.

이 포물선 프로파일에는 놀라운 결과가 있습니다. 파이프를 통과하는 총 유량은 파이프 반지름의 네제곱에 비례합니다. 반지름을 두 배로 하면 유량은 16배가 됩니다. 이것이 하겐-포아즈유 법칙이며, 아주 약간의 동맥 협착만으로도 혈류가 크게 감소하는 이유를 설명합니다.

가정: 포아즈유 유동은 유체가 비압축성이고 뉴턴 유체(점성도 일정)이며, 유동이 정상이고 완전히 발달(입구에서의 가속이 끝난 상태)되었으며, 유동이 층류 — 매끄럽고 질서 정연한 — 라고 가정합니다. 실제로 관 유동은 레이놀즈 수 약 2,300에서 난류로 전이합니다. 이 수치는 이론에서 유도된 것이 아닌 경험적 관찰입니다.

반지름 $R$의 원형 파이프에서 축방향 $x$로 균일한 압력 기울기 $dp/dx < 0$에 의해 구동되는 정상, 완전 발달, 축대칭 유동을 고려합니다. 원통 좌표 $(r, \theta, x)$에서 속도장은 $u = (0,\, 0,\, u(r))$ 형태입니다.

비압축성 조건 $\nabla \cdot u = 0$은 자동으로 만족됩니다. $u$가 $x$에 의존하지 않으므로 이류 항은 사라집니다. 축방향 운동량 방정식은 다음으로 축소됩니다.

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

동치하게 쓰면

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

$dp/dx$가 상수이므로, 이것은 $r$에 대한 ODE입니다. 경계 조건 $u(R) = 0$ (미끄러짐 없음)과 $du/dr|_{r=0} = 0$ (대칭성)으로 두 번 적분하면 포물선 속도 프로파일을 얻습니다:

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

중심선에서의 최대 속도는 $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$이고, 평균 속도는 $\bar{u} = u_{\max}/2$입니다.

단면적에 대해 적분하면 체적 유량에 대한 하겐-포아즈유 법칙을 얻습니다:

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

여기서 $\Delta p > 0$은 파이프 길이 $L$에 걸친 압력 강하입니다.

$R^4$ 의존성이 특징적입니다: 반지름의 작은 변화가 유량의 큰 변화를 야기합니다. 혈류역학에서 이것은 동맥 협착에 대한 혈류의 민감도를 뒷받침합니다.

유효 범위: 이 해는 비압축성, 뉴턴, 정상, 완전 발달, 층류 유동에 적용됩니다. 난류로의 전이는 실험적으로 $\text{Re} = \bar{u} D / \nu \approx 2{,}300$에서 일어나며, 여기서 $D = 2R$은 파이프 직경입니다. 이 임계값은 경험적 관찰이며 — 나비에-스토크스 방정식으로부터 예측하는 정리는 없습니다. 레이놀즈 수와 층류-난류 전이에 대한 논의는 레이놀즈 수, 난류, 그리고 왜 작은 스케일이 중요한가를 참조하십시오.

쿠에트 유동은 두 평행 평판 사이에 갇힌 유체에서, 한 평판이 움직이고 다른 평판은 정지해 있을 때의 엄밀해입니다.

테이블 위에 놓인 카드 한 벌을 상상해 보십시오. 맨 위 카드를 옆으로 끌면 아래 카드들도 따라 움직입니다 — 각각은 바로 위의 카드보다 조금 덜 움직입니다. 속도가 바닥의 0부터 위 평판의 속도까지 선형으로 변합니다.

이것이 단순 쿠에트 유동입니다: 직선적인 속도 프로파일. 아래 평판의 유체는 정지해 있고, 위 평판의 유체는 평판과 함께 움직이며, 그 사이의 모든 것은 선형으로 보간됩니다. 압력 기울기가 필요 없습니다 — 운동은 순전히 움직이는 경계에 의해 구동됩니다.

채널을 따라 압력 기울기도 함께 적용하면, 전단 구동과 압력 구동 유동의 조합이 됩니다. 속도 프로파일은 선형 프로파일 위에 포물선이 겹쳐진 형태가 됩니다. 이 결합 유동은 때때로 평면 포아즈유-쿠에트 유동이라 불립니다.

순수 압력 구동 — 두 평판 모두 정지, 압력 차이로 구동 — 의 경우가 평면 포아즈유 유동이며, 관 유동의 평판 대응물입니다. 프로파일은 포물선이고, 중간에서 가장 빠르며, 양쪽 벽에서 0입니다.

간격 $h$로 분리된 두 무한 평행 평판 사이의 정상 유동을 고려합니다. $y$를 평판에 수직인 좌표로, 아래 평판이 $y = 0$, 위 평판이 $y = h$에 있다고 합니다. 유동은 평행합니다: $u = (u(y),\, 0,\, 0)$.

단순 쿠에트 유동. 위 평판은 $x$ 방향으로 속도 $U$로 움직이고, 아래 평판은 정지; $dp/dx = 0$. 운동량 방정식은 $d^2u/dy^2 = 0$으로 축소되어 다음을 줍니다.

$$u(y) = U\frac{y}{h}.$$

이것은 나비에-스토크스 방정식의 가장 단순한 비자명 엄밀해입니다: 경계 조건에 의해서만 구동되는 선형 속도 프로파일.

평면 포아즈유 유동. 두 평판 모두 정지; 일정한 압력 기울기 $dp/dx < 0$이 유동을 구동합니다. 운동량 방정식은 다음이 됩니다.

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

$u(0) = u(h) = 0$. 해는 다음입니다.

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

중심선 $y = h/2$ 주위로 대칭인 포물선 프로파일입니다.

일반 쿠에트-포아즈유 유동. 움직이는 위 평판과 압력 기울기를 결합하면 중첩이 됩니다.

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

첫째 항은 전단 구동 기여이고, 둘째 항은 압력 구동 기여입니다. $dp/dx$의 부호와 크기에 따라 $\mu U / h^2$에 대해 프로파일은 단조, 내부 극대를 갖거나, 한쪽 벽 근처에서 역류까지 나타날 수 있습니다.

포아즈유와 쿠에트 유동은 정상 — 시간에 따라 변하지 않는 — 입니다. 스토크스 문제는 가장 단순한 비정상 엄밀해이며, 근본적인 사실을 드러냅니다: 점성은 운동량을 유체 속으로 확산시킵니다, 마치 열이 고체 속으로 확산하는 것과 같습니다.

스토크스 제1문제(레일리(Rayleigh) 문제라고도 함): 평평한 판 위에 광대한 정지 유체가 있다고 상상해 보십시오. 시간 0에서 판이 갑자기 옆으로 일정한 속도로 미끄러지기 시작합니다. 판 바로 옆의 유체는 즉시 끌려가지만, 더 먼 유체는 알아차리는 데 시간이 걸립니다. 매끄러운 경계층이 판에서 바깥쪽으로 점점 두꺼워지며 성장합니다.

판 위 임의의 높이에서의 속도는 그 높이와 특성 확산 길이 $\sqrt{\nu t}$ (여기서 $\nu$는 점성도, $t$는 경과 시간)의 비에 의존합니다. 유체가 점성이 클수록 운동이 위쪽으로 더 빨리 퍼집니다.

스토크스 제2문제: 같은 설정이지만, 이번에는 판이 일정한 속도가 아닌 사인파 형태로 앞뒤로 진동합니다. 진동은 유체 속으로 유한한 거리까지만 침투합니다 — 더 위에서는 유체가 거의 알아차리지 못합니다. 운동의 진폭이 높이에 따라 지수적으로 감소하여, 얇은 진동 경계층을 만듭니다. 이것이 깊은 수영장의 표면을 젓는 것이 표면 근처의 얕은 영역만 교란하는 이유입니다.

두 스토크스 문제 모두 평평한 판 $y = 0$ 위의 반무한 유체($y > 0$)를 다룹니다. 유동은 평행합니다: $u = (u(y,t),\, 0,\, 0)$. 이류 항은 사라지며, 지배 방정식은 1차원 확산 방정식입니다:

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

스토크스 제1문제 (레일리 문제). 초기 조건: $y > 0$에서 $u(y,0) = 0$. 경계 조건: $t > 0$에서 $u(0,t) = U$. 원방: $y \to \infty$에서 $u \to 0$.

유사 변수 $\eta = y / (2\sqrt{\nu t})$가 PDE를 ODE로 축소합니다. 해는 다음입니다.

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

여기서 $\operatorname{erfc}$는 여오차 함수입니다. 경계층 두께는 $\delta \sim \sqrt{\nu t}$로 성장하며, 확산 퍼짐의 특징입니다.

스토크스 제2문제. 판이 진동합니다: $u(0,t) = U\cos(\omega t)$. $u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$ 형태의 해를 구하면 다음을 얻습니다.

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

따라서

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

진폭은 침투 깊이 $\delta_s$에 따라 지수적으로 감쇠하고, 위상은 선형으로 지연됩니다 — 에너지가 점성에 의해 완전히 소산되는 횡파입니다. 주파수가 높을수록 침투 깊이가 얕아집니다.

포아즈유, 쿠에트, 스토크스 유동이 가장 유명하지만, 유일한 엄밀해는 아닙니다. 문헌에서 자주 등장하는 몇 가지가 더 있습니다:

• 테일러-그린 와류 — 2차원에서 감쇠하는 소용돌이 패턴. 진정한 와류 구조를 가진 유명한 엄밀해이며, 전산유체역학(CFD) 코드의 벤치마크로 널리 사용됩니다.
• 제프리-하멜(Jeffery-Hamel) 유동 — 수렴하거나 발산하는 쐐기형 채널에서의 유동. 유체가 좁아지는 틈으로 가속되거나 넓어지는 곳에서 감속되는 현상을 포착합니다.
• 히멘츠(Hiemenz) 정체점 유동 — 유체가 평평한 벽에 정면으로 충돌하는 유동. 유체가 벽에서 속도 0으로 감속되고 옆으로 비껴갑니다. 바람이 건물 정면에 부딪히거나 제트가 표면에 부딪힐 때의 현상을 모델링합니다.

이 각각은 특정한 기하학적 대칭을 활용하여 방정식을 다루기 쉽게 만듭니다. 이들은 특수한 맥락에서 중요하지만, 포아즈유와 쿠에트 유동이 입문 유체역학의 주력 해로 남아 있습니다.

평행 유동 외에도, 특정 대칭이나 자기 유사 구조를 활용하는 다른 엄밀해 계열이 있습니다:

• 테일러-그린 와류. 2차원에서 $u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$이며, $aA = bB$ (비압축성)이고 시간에 따라 지수적 감쇠 $e^{-\nu(a^2+b^2)t}$합니다. 이것은 비선형 항을 포함한 완전한 2차원 나비에-스토크스 방정식의 엄밀해입니다(비선형 항이 기울기가 되어 압력에 흡수됩니다). DNS 코드의 표준 검증 사례로 사용됩니다.
• 제프리-하멜 유동. 반각 $\alpha$인 쐐기에서의 정상, 순방사 유동 $u_r(r,\theta)$. 유동 함수 $\psi = \nu f(\theta)$가 $\theta$에 대한 3차 비선형 ODE를 만족합니다. 수렴 및 발산 채널 모두에서 해가 존재하며, 높은 레이놀즈 수에서 풍부한 분기 구조를 보입니다.
• 히멘츠 정체점 유동. 평면 벽에 충돌하는 유동에 대한 2차원 유사해. 유동 함수 $\psi = \sqrt{a\nu}\, x\, f(\eta)$, $\eta = y\sqrt{a/\nu}$가 나비에-스토크스를 히멘츠 ODE로 축소합니다: $f''' + ff'' - f'^2 + 1 = 0$, $f(0)=f'(0)=0$, $f'(\infty)=1$.

이 모든 경우에서 나비에-스토크스 방정식을 정확히 풀 수 있다면, 왜 아직 백만 달러짜리 미해결 문제가 있을까요?

모든 엄밀해가 특별한 기교에 의존하기 때문입니다. 기하 구조가 너무나 신중하게 선택되어 방정식의 가장 어려운 부분 — 비선형 항 — 이 완전히 사라지거나 다루기 쉬운 것으로 축소됩니다. 관 유동은 사실상 1차원입니다. 쿠에트 유동은 직선입니다. 테일러-그린 와류조차 비선형성을 압력 안에 숨깁니다.

밀레니엄 상금 문제는 어떤 특별한 대칭도 없는 일반적인 3차원 유동에 대해 묻습니다. 임의의 초기 조건, 단순화하는 기하 구조 없이, 모든 스케일에서 완전한 비선형 상호작용이 일어나는 상황입니다. 이 설정에서 해가 항상 매끄럽게 유지된다는 것도, 그렇지 않다는 것도 아직 증명되지 않았습니다.

이러한 엄밀해들은 강한 대칭 가정 아래에서 방정식이 명시적인 매끄러운 해를 허용한다는 것을 보여줍니다. 미해결 문제는 이것이 항상 작동하는지 — 아니면 3차원 난류의 완전한 일반성에서 무언가가 재앙적으로 잘못될 수 있는지입니다.

위에서 논의한 모든 엄밀해는 비선형 이류 항 $(u \cdot \nabla)u$를 제거하거나 자명하게 만들어 다루기 쉬움을 달성합니다. 평행 유동에서는 항등적으로 사라집니다. 테일러-그린 와류에서는 기울기가 되어 압력에 흡수됩니다. 유사해에서는 변수 변환을 통해 저차원 문제로 축소됩니다.

클레이 밀레니엄 상금 문제는 임의의 매끄럽고 빠르게 감쇠하는 초기 데이터에 대한 3차원 비압축성 나비에-스토크스 방정식의 초기값 문제에 관한 것입니다 — 이러한 단순화가 전혀 적용되지 않는 바로 그 영역입니다.

문제는 다음의 해가

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

$\mathbb{R}^3$ 위에서 매끄러운 초기 데이터 $u_0 \in C_c^\infty(\mathbb{R}^3)$가 주어졌을 때 모든 시간 동안 $C^\infty(\mathbb{R}^3 \times [0,\infty))$에 남아 있는지 여부입니다. 엄밀해는 완전한 비선형 역학이 결코 관여하지 않는 해 공간의 부분공간에 존재하기 때문에 이 질문에 답하지 못합니다.

요약하면: 엄밀해는 고도로 대칭적인 영역에서 명시적인 매끄러운 해를 제공합니다. 미해결 문제는 적정성이 제한 없는 3차원 유동으로 확장되는지 여부입니다. 정확한 정식화는 나비에-스토크스 존재와 매끄러움을 참조하십시오.

방정식 자체와 각 항의 의미를 이해하려면 나비에-스토크스 방정식이란?에서 시작하십시오.

이 방정식이 제1원리에서 어떻게 구축되는지 보려면 나비에-스토크스 방정식의 유도를 읽으십시오.

포아즈유와 쿠에트 같은 층류가 언제 난류로 붕괴되는지 이해하려면 레이놀즈 수, 난류, 그리고 왜 작은 스케일이 중요한가를 읽으십시오.

엄밀해가 답할 수 없는 백만 달러짜리 질문을 이해하려면 밀레니엄 상금 문제: 존재와 매끄러움을 읽으십시오.

이 사이트의 추천 다음 페이지:

• 나비에-스토크스 방정식이란? — 지배 시스템과 그 항들
• 나비에-스토크스 방정식의 유도 — 운동량 균형과 뉴턴 구성 관계에서 PDE까지
• 레이놀즈 수, 난류, 그리고 왜 작은 스케일이 중요한가 — 여기서 기술한 층류 영역에서 완전 난류 유동으로의 전이
• 밀레니엄 상금 문제: 존재와 매끄러움 — 엄밀해가 해결하지 못하는 미해결 문제의 정확한 정식화