비압축성 vs. 압축성 나비에-스토크스

"비압축성 대 압축성"은 실제로 밀도에 관한 질문입니다. 유체가 흐를 때, 밀도가 점마다 일정하게 유지되는가, 아니면 변하는가?

주사기에 물을 채우고 피스톤을 밉니다. 물은 움직이지만 압축되지 않습니다 — 같은 공간에 더 많은 물을 넣을 수 없습니다. 물은 거의 비압축성입니다. 이제 주사기에 공기를 채우고 끝을 막습니다. 피스톤을 밀면 공기가 압축되는 것을 느낄 수 있습니다. 같은 질량의 공기가 더 작은 부피를 차지합니다. 이것이 압축성 유동입니다.

비압축성 유동에서 밀도 $\rho$는 유체 전체에서 일정합니다. 유체의 미소 입자는 각각 움직이면서 부피를 유지합니다. 압축성 유동에서 밀도는 변수가 됩니다 — 장소와 시간에 따라 변할 수 있습니다. 제트 엔진 주위의 공기, 폭발 속의 가스, 대규모 대기: 이것들은 밀도 변화가 물리를 주도하는 압축성 유동을 수반합니다.

이 구분이 중요한 이유는 나비에-스토크스 방정식의 형태, 예측, 풀이의 난이도를 바꾸기 때문입니다.

비압축성 가정은 밀도장 $\rho$를 유동 영역 전체에서 양의 상수로 설정합니다. 물리적으로, 이것은 모든 물질 체적 요소가 유동 사상(flow map) 아래에서 부피를 보존한다는 의미입니다. 압축성 설정은 $\rho(x,t)$를 고유한 발전 방정식에 의해 지배되는 완전한 미지수로 승격시킵니다.

이것은 동일한 시스템의 두 표기법이 아닙니다. 미지수의 수, 제약 방정식의 구조, 압력의 성격이 다릅니다. 비압축성 시스템은 $d+1$개의 스칼라 미지수($\mathbb{R}^d$에서 $u$와 $p$)를 가지며; 압축성 시스템은 통상 $d+2$개의 주 필드(속도, 밀도, 내부 에너지 또는 온도)를 가지며, 압력은 상태 방정식을 통해 결정됩니다.

이 구분은 압력의 PDE 유형도 바꿉니다: 비압축성의 경우 타원형(압력이 즉각 조정)인 반면, 압축성의 경우 전체 시스템은 혼합 쌍곡-포물형이며 압력 교란이 유한 속도로 전파됩니다. 이것은 기술적 사항이 아닙니다 — 유체를 통해 정보가 이동하는 방식을 제어합니다.

비압축성 나비에-스토크스 방정식은 밀도가 일정한 유체를 기술합니다. 이것이 클레이 밀레니엄 문제에 등장하는 버전이며 이 사이트가 집중하는 버전입니다.

시스템은 두 부분으로 구성됩니다. 운동량 방정식:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

그리고 비압축성 제약:

$$\nabla \cdot u = 0$$

제약 $\nabla \cdot u = 0$은 속도장이 발산 없음을 의미합니다: 유체가 어디에서든 쌓이거나 희박해지지 않습니다. 미소 영역으로 들어오는 것은 같은 속도로 나가야 합니다. 이 하나의 조건이 밀도 방정식 전체를 대체합니다 — 밀도가 변하지 않으므로 이를 추적하는 방정식이 필요 없습니다.

이 시스템에서 압력은 특별한 역할을 합니다. 열역학적 법칙(예: 이상 기체 법칙)에 의해 결정되는 것이 아닙니다. 대신, 유동을 발산 없이 유지하기 위해 모든 곳에서 즉각 조정됩니다. 수학적으로, $p$는 제약에서 도출된 포아송 방정식을 풉니다. 이것은 압력 변화가 무한히 빠르게 전파된다는 의미입니다 — 비압축성 유동에는 "음속"이 없습니다.

비압축성 나비에-스토크스 시스템은 두 개의 미지 필드를 가집니다: 속도장 $u$와 압력장 $p$. 이 상대적 단순함은 기만적입니다 — 비선형 항 $(u \cdot \nabla)u$가 여전히 3차원에서 시스템을 극히 어렵게 만듭니다.

동점성계수 $\nu > 0$을 가진 $\mathbb{R}^3$ 위의 비압축성 나비에-스토크스 시스템:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

발산 없는 조건 $\nabla \cdot u = 0$은 점별 제약이지 발전 방정식이 아닙니다. 이것은 국소적 체적 보존을 부호화합니다: 유동 사상 $\Phi_t$는 모든 $t$에 대해 $\det(D\Phi_t) = 1$을 만족합니다.

운동량 방정식에 발산 연산자를 적용하고 비압축성을 사용하면 압력 포아송 방정식을 얻습니다:

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

이것은 각 고정 시간에서 $p$에 대한 타원형 방정식입니다. 압력은 독립적인 열역학적 변수가 아니며; 속도장에 의해 전역적이고 즉각적으로 결정됩니다. 압력을 통한 정보는 무한 속도로 전파됩니다 — 압축성 시스템과의 근본적 구조적 차이입니다.

미지수는 $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$와 $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$입니다. 클레이 공식(Fefferman, 2000)에서, $u_0 \in C^\infty(\mathbb{R}^3)$는 발산 없이이며 질문은 $u$가 유계 에너지를 가지며 $C^\infty(\mathbb{R}^3 \times [0,\infty))$에 머무는지 여부입니다.

압축성 나비에-스토크스 방정식은 밀도가 변하는 유동을 지배합니다. 더 많은 미지수와 더 많은 방정식을 가진 더 크고 복잡한 시스템입니다.

여전히 운동량 방정식이 있지만, 이제 밀도 $\rho$가 명시적으로 나타납니다:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

제약 $\nabla \cdot u = 0$은 사라집니다. 그 자리에 밀도가 어떻게 변하는지를 추적하는 연속 방정식을 얻습니다:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

이것은 질량이 보존된다는 뜻입니다: 밀도는 유동이 유체 입자를 압축하거나 팽창시킴에 따라 변합니다.

시스템은 또한 에너지 방정식과 상태 방정식 — 압력을 밀도와 온도에 연결하는 열역학적 관계(예: $p = \rho R T$, 이상 기체 법칙) — 을 필요로 합니다. 압력은 더 이상 제약의 수동적 강제자가 아닙니다. 고유한 물리학, 고유한 역학을 가지며, 유한 속도 — 음속 — 로 전파됩니다.

압축성 시스템은 고속 공기역학, 천체물리학적 기체 역학, 연소, 그리고 밀도 변화가 중요한 모든 유동에 필수적입니다. 그러나 비압축성 방정식과는 수학적으로 진정 다른 대상입니다 — 더 많은 미지수, 더 많은 방정식, 다른 PDE 구조.

압축성 나비에-스토크스 시스템은 속도 $u(x,t)$, 밀도 $\rho(x,t)$, 압력 $p(x,t)$, 비내부 에너지 $e(x,t)$(또는 온도 $\theta$)를 결합합니다. 보존 형태:

연속: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

운동량: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

에너지: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

여기서 $E = e + \tfrac{1}{2}|u|^2$는 총 비에너지, $\tau$는 점성 응력 텐서(뉴턴 유체의 경우, $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$이며 체적 점성 $\lambda$), $\kappa$는 열전도율입니다.

폐합에는 상태 방정식이 필요합니다. 예: 단열 지수 $\gamma$를 가진 이상 기체의 경우 $p = (\gamma - 1)\rho e$.

압축성 시스템에서 음향 교란은 유한 속도로 전파되며, 특성 음속 $c = \sqrt{\partial p / \partial \rho |_s}$를 가집니다. 이것은 비압축성 모델의 타원형이고 즉각적인 압력 결정과 대조됩니다. 압축성 방정식은 비압축성 유동에서 대응물이 없는 충격파, 팽창파, 접촉 불연속을 지원합니다.

실제로 압축성 유동과 비압축성 유동을 나누는 기준선은 마하 수입니다:

$$\text{Ma} = \frac{|u|}{c}$$

여기서 $|u|$는 유동 속도이고 $c$는 음속입니다. 마하 수는 유동이 압력 교란의 전파 속도에 비해 얼마나 빠르게 이동하는지를 측정합니다.

경험적 규칙: $\text{Ma} < 0.3$이면 밀도 변화는 약 5% 미만이며, 비압축성 방정식이 훌륭한 근사입니다. 방 안의 공기, 파이프 속의 물, 건물 주위의 바람 — 이것들은 모두 낮은 마하 수 유동입니다. 공기는 기술적으로 압축성이지만, 마치 아닌 것처럼 행동합니다.

$\text{Ma} \approx 0.3$ 이상에서는 압축성 효과가 커집니다. $\text{Ma} \approx 1$ 부근의 천음속 영역에서는 국소적 초음속 영역이 나타나고 충격파가 형성되기 시작합니다. 전투기, 로켓 노즐, 재진입 우주선이 이 영역에 속합니다.

이것은 이진 스위치가 아닌 스펙트럼입니다. 낮은 마하 수 영역은 비압축성 나비에-스토크스 방정식이 물리적 모델로 적용되는 곳입니다. 대부분의 일상적 유체 유동 — 그리고 클레이 밀레니엄 문제 — 은 바로 이 영역에 속합니다.

마하 수 $\text{Ma} = |u|/c$는 압축성의 중요도를 매개변수화하며, 여기서 $c = \sqrt{\partial p / \partial \rho |_s}$는 등엔트로피 음속입니다. 형식적으로, 비압축성 방정식은 압축성 시스템의 낮은 마하 수 극한으로 발생합니다.

$\text{Ma}^2$의 거듭제곱에 대한 점근 전개(Klainerman & Majda, 1981, 1982; Schochet, 1986 참조)는 적절한 초기 데이터 하에서 $\text{Ma} \to 0$일 때 압축성 해가 비압축성 해로 수렴함을 보여줍니다. 잘 준비된 데이터를 가진 표준 낮은 마하 수 무차원화에서, 압력은 흔히 거의 공간적으로 균일한 열역학적 배경에 극한에서 비압축성 제약을 강제하는 더 작은 동적 보정을 더한 것으로 쓰입니다.

이것은 특이 극한입니다: 음속 $c \to \infty$이며 압축성 시스템의 쌍곡적 성격이 비압축성 유동의 타원형 압력 방정식으로 퇴화합니다. 음향 모드는 무한히 빨라지며 와류 역학으로부터 분리됩니다.

$\text{Ma} < 0.3$ 영역은 공학적 경험칙입니다. 이것은 이 범위에서 상대적 밀도 변화 $\delta\rho / \rho \sim \text{Ma}^2 / 2$가 약 5% 미만으로 유지된다는 경험적 관찰을 반영합니다. 수학적 정당화는 잘 준비된 초기 데이터와 양립하는 경계 조건을 필요로 하는 낮은 마하 수 극한의 수렴 정리입니다.

클레이 밀레니엄 문제는 정확한 질문을 합니다: 매끄러운 발산 없는 초기 속도가 $\mathbb{R}^3$에서 주어졌을 때, 비압축성 나비에-스토크스 시스템은 항상 모든 시간에 대해 존재하는 매끄러운 해를 산출하는가?

왜 구체적으로 비압축성인가? 세 가지 이유입니다.

첫째, 이미 충분히 어렵습니다. 비압축성 3차원 방정식은 1934년 르레의 기초 작업 이후 전역적 정칙성의 증명에 저항해 왔습니다. 가변 밀도, 열역학, 충격파의 복잡성을 추가하면 문제가 더 다루기 쉬워지는 것이 아니라 훨씬 어려워질 것입니다.

둘째, 난제는 순수한 유체역학입니다. 비압축성 시스템은 핵심 수학적 도전 — 비선형 이류 $(u \cdot \nabla)u$와 점성 소산 $\nu \Delta u$ 사이의 상호작용 — 을 열역학적 또는 음향적 복잡성 없이 고립시킵니다. 비압축성 경우는 이미 충분히 어려우므로 압축성, 열역학, 충격을 추가하기 전에 정칙성 문제를 고립시키기에 자연스러운 장소입니다.

셋째, 물리학이 깔끔합니다. 비압축성 방정식은 가장 흔한 일상적 유동을 모델링합니다. 매끄러운 데이터로부터 특이점을 생성할 수 있는지 이해하는 것은 고전적 유체역학의 수학적 일관성에 관한 근본적 질문입니다.

압축성 시스템은 자체적으로 깊은 미해결 문제를 가지고 있습니다 — 대규모 데이터에 대한 전역 해의 존재성, 충격의 형성과 상호작용 — 그러나 이것들은 다른 구조를 가진 다른 문제입니다. 클레이 상금은 비압축성 경우를 대상으로 하는데, 이것이 클레이가 3차원 나비에-스토크스에 대해 공식화한 특정 정칙성 문제이기 때문입니다.

공식 클레이 공식(Fefferman, 2000)은 $\mathbb{R}^3$ 위의 비압축성 시스템을 명시합니다:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

$u_0 \in C^\infty(\mathbb{R}^3)$가 발산 없이이며 질문은 $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$이면서 모든 $t \geq 0$에 대해 $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$가 유계인지 여부입니다.

비압축성 시스템의 선택은 수학적으로 동기가 부여됩니다. 핵심 미해결 문제 — 전역적으로 존재하지만 유일하거나 매끄럽지 않을 수 있는 르레-호프 약한 해와 국소적으로 존재하지만 폭발할 수 있는 고전적 매끄러운 해 사이의 간극 — 은 비압축성 3차원 방정식에 특유합니다. 2차원에서는 매끄러운 비압축성 나비에-스토크스 해의 전역적 정칙성이 알려져 있으며; 3차원은 미해결입니다.

압축성 시스템은 질적으로 다른 난제를 도입합니다: 충격 형성(매끄러운 데이터에 대해서도 오일러 방정식에서 발생), 진공 상태($\rho \to 0$), 그리고 와도와 음향 모드 사이의 결합. 이것들은 중요한 미해결 문제이지만, 비압축성 정칙성 문제와는 구조적으로 구별됩니다.

비압축성 문제는 에너지 초임계 비선형 항과 점성 소산 사이의 경쟁을 고립시킵니다. 자연스러운 에너지 추정은 $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$를 주며, 이것은 3차원에서 스케일링 임계 공간 $L^\infty_t \dot{H}^{1/2}_x$에 반 도함수만큼 미치지 못합니다. 이 간극을 좁히는 것 — 또는 좁힐 수 없음을 증명하는 것 — 이 밀레니엄 문제의 핵심입니다.

완전한 비압축성 시스템과 각 항의 의미에 대해서는 나비에-스토크스 방정식이란?을 읽으십시오.

방정식이 뉴턴의 제2법칙과 응력 텐서로부터 어떻게 구성되는지에 대해서는 나비에-스토크스 방정식의 유도를 참고하십시오.

점성을 완전히 제거하고 오일러 방정식을 얻을 때 어떤 일이 일어나는지에 대해서는 오일러 vs. 나비에-스토크스를 참고하십시오.

밀레니엄 문제의 정확한 명제와 "존재성과 매끄러움"이 실제로 의미하는 바에 대해서는 나비에-스토크스 존재성과 매끄러움 문제를 참고하십시오.

추천 경로:

• 나비에-스토크스 방정식이란? — 항별 분석과 함께하는 완전한 비압축성 시스템
• 나비에-스토크스 방정식의 유도 — 연속체 역학 공리에서 PDE 시스템까지
• 오일러 vs. 나비에-스토크스 — 비점성 극한 $\nu \to 0$과 정칙성에서의 점성의 역할
• 나비에-스토크스 존재성과 매끄러움 문제 — 정확한 클레이 공식, 르레-호프 이론, 그리고 약한 해와 강한 해 사이의 간극