나비에-스토크스 방정식의 유도

모든 나비에-스토크스 유도는 같은 아이디어에서 출발합니다: 뉴턴의 제2법칙 — 힘은 질량 곱하기 가속도 — 을 움직이는 유체의 작은 덩어리에 적용하는 것입니다.

물, 공기, 또는 다른 어떤 유체의 작은 덩어리를 골라 보십시오. 그 덩어리에는 질량이 있습니다. 여러 힘이 작용합니다: 압력이 사방에서 짓누르고, 내부 마찰이 잡아당기며, 중력이 아래로 끌어당깁니다. 뉴턴의 법칙에 따르면, 알짜힘이 이 덩어리의 가속도를 결정합니다.

강체와 달리, 유체 덩어리는 움직이면서 변형될 수 있습니다. 흐름을 타면서 늘어나고 일그러집니다. 따라서 "가속도"는 단일 물체를 추적하는 것만큼 단순하지 않습니다. 덩어리가 흐름을 타고 이동하는 것을 따라가야 합니다.

일상적인 비유를 들면: 강 위에 카누를 타고 있다고 상상해 보십시오. 여러분의 가속도는 현재 위치에서 시간에 따른 물살의 변화뿐만 아니라, 흐름이 여러분을 더 빠르거나 느린 지역으로 데려가고 있다는 사실에도 좌우됩니다. 두 효과 모두 속도 변화에 기여합니다.

이 균형 — 모든 힘이 질량 곱하기 이 결합된 가속도와 같다는 것 — 을 유체의 모든 점에서 동시에 쓰는 것이 나비에-스토크스 방정식 유도의 출발점입니다.

나비에-스토크스 방정식의 유도는 코시(Cauchy) 운동량 방정식에서 시작합니다. 이것은 유체와 함께 움직이는 물질 체적 $\Omega(t)$에 대해 뉴턴의 제2법칙을 연속체역학 형태로 적용한 것입니다.

밀도 $\rho$와 속도장 $u(x,t)$를 갖는 유체에 대해, 선형 운동량 보존은 다음과 같이 서술됩니다.

$$\frac{d}{dt}\int_{\Omega(t)} \rho\, u\, dV = \int_{\partial\Omega(t)} T\, n\, dS + \int_{\Omega(t)} \rho\, f\, dV,$$

여기서 $T$는 코시 응력 텐서, $n$은 바깥 방향 단위 법선 벡터, $f$는 단위 질량당 체적력(일반적으로 중력)을 나타냅니다.

레이놀즈 수송 정리를 사용하여 국소화하면 미분 형태가 됩니다.

$$\rho\frac{Du}{Dt} = \nabla \cdot T + \rho f,$$

여기서 물질 도함수는

$$\frac{Du}{Dt} = \partial_t u + (u \cdot \nabla)u$$

로, 국소적 시간 변화율과 대류 가속도를 모두 포착합니다. 이것은 연속체의 각 점에서 $ma = F$의 정확한 대응물입니다.

나비에-스토크스 방정식의 유도에는 세 종류의 힘이 등장합니다:

1. 압력. 유체가 덩어리를 사방에서 밀어냅니다. 한쪽의 압력이 다른 쪽보다 높으면, 덩어리는 낮은 압력 쪽으로 밀려납니다. 이 불균형은 압력 기울기 $-\nabla p$로 표현됩니다: 높은 압력에서 낮은 압력 방향을 가리키며 유체가 어느 방향으로 가속할지를 알려주는 벡터입니다.

2. 점성력. 서로 다른 속도로 움직이는 인접한 유체 층이 서로를 끌어당깁니다 — 손을 비비는 것과 비슷하지만, 유체 전체에서 연속적으로 일어납니다. 빠른 층은 느린 이웃을 끌어당기고, 느린 층은 빠른 이웃을 잡아둡니다. 이 내부 마찰은 속도 차이를 매끄럽게 만듭니다. 물과 같은 단순한 유체에서 이 마찰의 세기는 하나의 숫자, 즉 점성도 $\mu$로 결정됩니다.

3. 외력. 외부에서 유체에 작용하는 모든 힘 — 가장 흔한 것은 중력입니다. 이것은 체적력으로, 표면이 아닌 체적 내 모든 유체에 작용합니다.

나비에-스토크스 방정식은 모든 점에서 "질량 × 가속도 = 압력 + 점성력 + 외력"을 쓴 결과물입니다.

코시 응력 텐서 $T$는 모든 내부 접촉력을 부호화합니다. 임의의 유체에 대해 등방성 압력 부분과 편향(점성) 부분으로 분해됩니다:

$$T = -pI + \tau,$$

여기서 $p$는 역학적 압력(압축성 유동의 경우 $-\tfrac{1}{3}\mathrm{tr}\,T$로 정의)이고, $\tau$는 점성 응력 텐서입니다.

압력 기울기. $-pI$에 발산 정리를 적용하면 단위 체적당 힘 $-\nabla p$를 얻습니다. 이것은 응력 발산의 등방성 부분입니다.

점성 응력. 텐서 $\tau$는 응력을 변형률에 연결하는 구성 관계에 따라 달라집니다. 그 발산 $\nabla \cdot \tau$가 단위 체적당 점성력을 줍니다. 이 단계에서 $\tau$의 형태는 아직 지정되지 않았으며, 다음 절의 뉴턴 유체 가정에서 결정됩니다.

체적력. 중력과 같은 외력은 단위 체적당 $\rho f$로 들어옵니다. 응력 분해를 코시 운동량 방정식에 대입하면 다음을 얻습니다.

$$\rho\frac{Du}{Dt} = -\nabla p + \nabla \cdot \tau + \rho f.$$

이것은 뉴턴 유체든 아니든 임의의 단순 유체에 대한 일반 운동량 방정식입니다. 시스템을 닫으려면 $\tau$를 지정하는 구성 관계가 필요합니다.

모든 유체가 응력에 같은 방식으로 반응하는 것은 아닙니다. 꿀은 물과 다르게 운동에 저항합니다. 케첩은 흔들면 더 잘 흐릅니다. 옥수수 전분을 물에 섞으면 치면 더 단단해집니다.

나비에-스토크스 방정식은 특정한 가정을 합니다: 유체가 뉴턴 유체라는 것입니다. 이는 내부 마찰(점성 응력)이 유체가 변형되는 속도에 정비례한다는 의미입니다. 변형률을 두 배로 하면 응력도 두 배가 됩니다. 구성 가정은 선형입니다: 점성 응력이 순간 변형률에 비례합니다.

물, 공기, 그리고 많은 일반적인 유체는 뉴턴 유체로 매우 잘 모델링됩니다. 그러나 이것은 진정한 가정이지, 뉴턴 법칙의 결과가 아닙니다. 유도에 이 가정이 필요합니다. 이것이 없으면 다른 종류의 방정식(비뉴턴 유체 모델)이 됩니다.

비례 상수가 동적 점성도 $\mu$입니다: 유체가 전단에 얼마나 저항하는지를 측정하는 물질 특성입니다. 물은 점성도가 낮고, 꿀은 점성도가 높습니다.

유체가 압축되거나 팽창할 때 중요한 두 번째 점성도 매개변수 $\lambda$도 있는데, 제2 점성 계수라고도 합니다. 흔히 사용되는 추가적 단순화인 스토크스 가설은 $\lambda = -\frac{2}{3}\mu$로 설정합니다. 이것은 정리가 아닌 추가 가정입니다.

뉴턴 유체에 대한 구성 관계는 점성 응력 $\tau$가 변형률 텐서 $D(u) = \tfrac{1}{2}(\nabla u + \nabla u^T)$의 선형, 등방성 함수라고 가정합니다. 등방성 텐서 함수에 대한 표현 정리에 의해, 가장 일반적인 형태는 다음과 같습니다.

$$\tau = 2\mu\, D(u) + \lambda (\nabla \cdot u)\, I,$$

이를 동치하게 쓰면

$$\tau = \mu(\nabla u + \nabla u^T) + \lambda(\nabla \cdot u)\, I,$$

여기서 $\mu > 0$은 동적(전단) 점성도이고, $\lambda$는 제2 점성 계수입니다.

이것이 뉴턴 유체의 정의적 가정입니다. 제1원리에서 유도되는 것이 아니라, 많은 일반 유체에 대해 경험적으로 검증된 구성 가설입니다.

스토크스 가설은 $\lambda = -\frac{2}{3}\mu$로 놓아 체적 점성도 $\kappa = \lambda + \frac{2}{3}\mu$를 0으로 만듭니다. 이 단순화는 널리 사용되지만, 뉴턴 가설 자체나 열역학의 결과가 아닌 독립적인 가정입니다. 고전 기체 운동론은 이상화된 가정 하에 단원자 이상 기체의 체적 점성도가 0이라고 예측하지만, 많은 실제 기체와 액체에서 스토크스 가설은 근사적일 뿐입니다.

클라우지우스-뒤엠(Clausius-Duhem) 부등식으로부터의 열역학적 제약은 $\mu \geq 0$과 $3\lambda + 2\mu \geq 0$ (동치하게 $\kappa \geq 0$)만을 요구합니다.

이제 조각들을 하나로 합칩니다. 다음이 준비되어 있습니다:

• 좌변에 질량 곱하기 가속도 (속도의 물질 도함수)
• 우변에 압력, 점성 마찰, 중력
• 마찰과 변형률을 연결하는 뉴턴 유체 규칙

대입하고 정리하면 압축성 나비에-스토크스 운동량 방정식이 됩니다:

$$\rho\Big(\frac{\partial u}{\partial t} + (u \cdot \nabla) u\Big) = -\nabla p + \mu\, \Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u) + \rho\, f$$

각 항은 물리적 의미를 갖습니다:

• $\rho\,\partial_t u$ — 고정된 점에서 속도가 시간에 따라 변하는 것
• $\rho(u \cdot \nabla) u$ — 유체가 자신의 속도를 이곳저곳으로 실어 나르는 것 (이류)
• $-\nabla p$ — 높은 압력에서 낮은 압력 쪽으로 미는 것
• $\mu\,\Delta u$ — 점성이 속도 차이를 매끄럽게 만드는 것
• $(\mu + \lambda)\nabla(\nabla \cdot u)$ — 유체가 압축되거나 팽창할 때만 중요한 추가 점성 항
• $\rho f$ — 중력 같은 외력

이것이 완전한 운동량 방정식입니다. 질량 보존과 상태 방정식(압력과 밀도의 관계)을 함께 쓰면 압축성 나비에-스토크스 시스템이 됩니다.

뉴턴 구성 관계 $\tau = 2\mu D(u) + \lambda(\nabla \cdot u)I$를 운동량 방정식 $\rho \frac{Du}{Dt} = -\nabla p + \nabla \cdot \tau + \rho f$에 대입하고 $\tau$의 발산을 계산합니다:

$$\nabla \cdot \tau = \nabla \cdot \big[\mu(\nabla u + \nabla u^T)\big] + \nabla\big[\lambda(\nabla \cdot u)\big].$$

$\mu$와 $\lambda$가 상수이면(많은 유도에서의 표준 가정), 이것은 다음으로 단순화됩니다.

$$\nabla \cdot \tau = \mu\,\Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u),$$

벡터 항등식 $\nabla \cdot (\nabla u^T) = \nabla(\nabla \cdot u)$를 사용합니다. 압축성 나비에-스토크스 운동량 방정식은 다음과 같습니다.

$$\rho\big(\partial_t u + (u \cdot \nabla)u\big) = -\nabla p + \mu\,\Delta u + (\mu + \lambda)\,\nabla(\nabla \cdot u) + \rho f.$$

이 방정식은 질량 보존(연속 방정식)

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

과 상태 방정식 $p = p(\rho, \theta)$ (또는 에너지 방정식)과 결합되어야 시스템이 닫힙니다. 스토크스 가설 $\lambda = -\frac{2}{3}\mu$에서는 계수 $\mu + \lambda = \frac{1}{3}\mu$가 됩니다. 비압축성 시스템과의 상세한 비교는 비압축성 대 압축성 나비에-스토크스를 참조하십시오.

많은 유동 — 파이프 속의 물, 느린 기류, 해양 순환 — 에서 유체의 밀도는 본질적으로 일정합니다. 이 단순화를 적용하면 일반 방정식이 더 깔끔하고 다루기 쉬운 시스템으로 변환됩니다.

일정한 밀도란 $\rho$가 어디서나 항상 같다는 의미입니다. 그러면 질량 보존에 의해 속도장이 $\nabla \cdot u = 0$을 만족해야 합니다: 유체가 어디서도 압축되거나 팽창하지 않습니다. 이것이 비압축성 조건입니다.

$\nabla \cdot u = 0$이면 추가 점성 항 $(\mu + \lambda)\nabla(\nabla \cdot u)$가 완전히 사라집니다. 제2 점성 계수 $\lambda$는 역할을 잃습니다 — 유체가 압축될 수 없을 때는 필요 없습니다. 운동량 방정식은 다음이 됩니다:

$$\rho\Big(\frac{\partial u}{\partial t} + (u \cdot \nabla)u\Big) = -\nabla p + \mu\,\Delta u + \rho f$$

양변을 $\rho$로 나누고, $\nu = \mu / \rho$ (동점성도)를 쓰며, 압력을 재정의하여 $1/\rho$ 인수를 흡수하면 표준 형태가 됩니다:

$$\frac{\partial u}{\partial t} + (u \cdot \nabla)u = -\nabla p + \nu\,\Delta u + f$$

$$\nabla \cdot u = 0$$

이것이 클레이 밀레니엄 문제에서 연구하는 시스템입니다. 이 사이트 전체와 대부분의 나비에-스토크스 수학적 논의에서 만나게 될 형태입니다. 이 형태와 압축성 시스템의 상세한 비교는 비압축성 대 압축성 나비에-스토크스를 참조하십시오.

$\nu = 0$ (점성이 전혀 없는 경우)을 설정하면 오일러 방정식이 됩니다 — 관련되지만 중요하게 다른 시스템입니다.

일반적인 비압축성 특수화는 유동 전체에서 밀도 $\rho$가 일정하다고 가정합니다. 연속 방정식 $\partial_t \rho + \nabla \cdot(\rho u) = 0$에서 $\nabla \cdot u = 0$이 강제됩니다. 더 일반적으로, 비압축성은 $\nabla \cdot u = 0$ (연속 방정식을 통해 $D\rho/Dt = 0$과 동치)으로 부호화되는데, 이는 원칙적으로 변밀도도 허용하지만, 일정 밀도 경우가 클레이 문제의 표준 설정입니다.

비압축성 하에서 항 $(\mu + \lambda)\nabla(\nabla \cdot u)$는 항등적으로 사라집니다. 제2 점성 계수 $\lambda$는 무관해집니다 — 스토크스 가설이나 $\lambda$에 대한 어떤 가정도 비압축성 시스템에서는 필요 없습니다. 운동량 방정식은 다음으로 축소됩니다.

$$\rho\big(\partial_t u + (u \cdot \nabla)u\big) = -\nabla p + \mu\,\Delta u + \rho f.$$

$\rho$로 나누고 동점성도 $\nu = \mu/\rho$를 정의하면 (그리고 $p$를 $1/\rho$ 인수를 흡수하도록 재정의하면) $\mathbb{R}^3$ 위의 비압축성 나비에-스토크스 시스템이 됩니다:

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

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

이것이 클레이 밀레니엄 문제 (Fefferman, 2000)에서 연구하는 시스템입니다. 공식 문제는 표준 클레이 정식화(Fefferman, 2000)에서 $\mathbb{R}^3$ (또는 $\mathbb{T}^3$) 위의 매끄러운 발산 없는 초기 데이터에 대한 전역 정칙성에 관한 것입니다.

압력 $p$는 독립적인 동역학 변수가 아닙니다: 운동량 방정식의 발산을 취하고 $\nabla \cdot u = 0$을 사용하면 푸아송 방정식 $\Delta p = -\nabla \cdot [(u \cdot \nabla)u] + \nabla \cdot f$ (또는 $f$가 발산 없거나 없을 때 $\Delta p = -\partial_i \partial_j(u_i u_j)$)이 되어 (상수까지) 결정됩니다. 이 타원형 결합은 비압축성 시스템의 특징적인 성질입니다.

$\nu = 0$으로 설정하면 비압축성 오일러 방정식을 회복합니다. 이 두 시스템의 관계는 수학적 유체역학의 많은 미해결 문제의 핵심입니다.

위의 유도는 나비에-스토크스 방정식을 제공합니다. 방정식이 무엇인지, 왜 그런 형태를 취하는지 알려줍니다. 모든 항은 물리적 원리나 명시적 가정으로 추적됩니다.

그러나 방정식을 유도하는 것은 그 해를 이해하는 것과 같지 않습니다. 유도는 다음에 답하지 않습니다:

• 해가 항상 모든 시간에 존재하는가?
• 매끄럽게 시작하면 매끄러운 상태를 유지하는가?
• 속도가 유한 시간 안에 무한대로 폭발할 수 있는가?

이것은 방정식의 물리적 기원이 아닌 수학적 행동에 대한 질문입니다. 그리고 3차원에서 이 질문들은 미해결 상태입니다. 이 간극이 클레이 밀레니엄 문제로 이어집니다: 매끄러운 3차원 비압축성 나비에-스토크스 데이터가 항상 전역 매끄러운 해를 만들어 내는지, 아니면 유한 시간 안에 특이점이 형성될 수 있는지.

클레이 수학연구소는 그 해결에 백만 달러를 현상금으로 내걸고 있습니다. 175년이 넘도록 이 문제는 풀리지 않고 있습니다.

유도는 나비에-스토크스 방정식이 운동량 보존과 뉴턴 구성 관계에 기반한, 물리적으로 잘 동기부여된 PDE 시스템임을 확립합니다. 그러나 수학적 적정성(well-posedness)의 핵심 질문에는 답하지 않습니다.

구체적으로, 유도는 다음에 대해 침묵합니다:

• 전역 존재: 매끄러운 초기 데이터가 주어졌을 때 매끄러운 해가 모든 $t > 0$에서 지속하는지 여부.
• 정칙성: 해가 $C^\infty(\mathbb{R}^3 \times [0,\infty))$에 남아 있는지 아니면 특이점이 발생할 수 있는지 여부.
• 유일성: 다양한 함수 공간 설정에서 해가 유일한지 여부.

르레(Leray) 이론(1934)은 $L^2$에서 약해의 전역 존재를 보장하지만, 3차원에서 르레-홉프(Leray-Hopf) 해의 유일성과 정칙성은 여전히 증명되지 않았습니다. 가용한 에너지 추정값과 비선형성의 스케일링 사이의 간극이 이 어려움의 해석학적 핵심입니다.

클레이 밀레니엄 문제는 $\mathbb{R}^3$에서의 전역 정칙성에 대한 증명 또는 반증을 정확히 요구합니다. 유도는 PDE 시스템에 동기를 부여하지만, 3차원에서의 존재, 정칙성, 유일성을 해결하지는 않습니다.

방정식의 유래를 보셨으니, 다음으로 갈 곳은 여기입니다:

• 나비에-스토크스 방정식이란? — 항별 설명을 포함한 방정식 자체
• 비압축성 대 압축성 나비에-스토크스 — 밀도가 변할 때 무엇이 달라지는지, 그리고 비압축성 경우가 왜 특별한지
• 오일러 방정식 대 나비에-스토크스 방정식 — 점성을 완전히 제거하면 무슨 일이 일어나는지
• 밀레니엄 문제 — 이 방정식을 유명하게 만든 미해결 문제

추가 학습을 위한 관련 페이지:

• 나비에-스토크스 방정식이란? — 각 항에 대한 상세한 논의를 포함한 표준 형태
• 비압축성 대 압축성 나비에-스토크스 — 완전한 압축성 시스템, 순압 닫힘, 저마하수 극한
• 오일러 방정식 대 나비에-스토크스 방정식 — 비점성 극한, $\nu$의 역할, 점성 소멸 문제
• 나비에-스토크스 존재와 매끄러움 — 정확한 클레이 정식화, 알려진 결과, 미해결 문제