본문으로 건너뛰기
Park JiHo Research notes · Seoul
#frontier-one#ai-mathematics#formal-verification#navier-stokes

나비에–스토크스, 소용돌이 그림에서 수식까지

나비에–스토크스 방정식의 각 항부터 에너지 수지와 원문 Figure의 수축률까지 공대생 눈높이로 연결한다. 공식 형식화의 실제 좌표·속도식을 읽고, 차수 추정과 전체 증명을 구별한다.

#1. 이 글에서 이해할 질문

처음 접하는 주제라면 Blogger의 쉬운 해설: Navier–Stokes 방정식이란? AI가 난제를 풀었다는 말 쉽게 보기부터 읽어도 좋다. 그 글에서는 소용돌이와 특이점의 뜻을 먼저 살펴보고, 여기서는 같은 질문을 물질미분·에너지식·수축 좌표로 이어서 계산한다.

커피를 저으면 소용돌이가 생긴다. 그러나 이 글의 목표는 ‘소용돌이는 복잡하다’에서 끝나는 것이 아니다. 왜 전체 에너지가 유한하다는 사실만으로 모든 위치의 속도가 유한하다고 결론 내릴 수 없는지, 그리고 속도를 크게 만드는 것보다 매끄러운 외력을 남기는 일이 왜 더 어려운지를 수식으로 이해하는 것이 목표다.

고등학교의 지수법칙, 미분·적분, 대학 기초의 벡터와 연쇄법칙에서 시작한다. 먼저 방정식을 읽고, 그림을 해석한 다음, 공식 형식화에 기록된 좌표와 속도식을 따라간다. ‘쉽다’는 말은 수식을 없앤다는 뜻이 아니다. 수식이 등장하기 전에 필요한 도구를 준비한다는 뜻이다.

핵심 질문은 하나다. 유한한 에너지와 매끄러운 외력이 있어도, 흐름이 한 작은 영역에 집중되어 매끄러운 해를 더는 유지할 수 없을까?

#2. 유체에 적용한 F = ma

물 한 덩어리를 입자 하나로 취급하면 그 안에서 서로 다르게 움직이는 부분을 놓친다. 그래서 공간의 각 위치에 속도를 배정한다. 이 함수가 속도장(velocity field)이다.

u(x,t)=[ux(x,t)uy(x,t)uz(x,t)].\mathbf u(\mathbf x,t)=\begin{bmatrix}u_x(\mathbf x,t)\\u_y(\mathbf x,t)\\u_z(\mathbf x,t)\end{bmatrix}.

위치 x와 시각 t를 넣으면 그곳의 유체가 어느 방향으로 얼마나 빨리 움직이는지가 나온다. 압력도 위치와 시간의 함수다. 밀도와 점성계수가 일정한 비압축성 뉴턴 유체에서는 운동량 방정식을 다음과 같이 쓸 수 있다. [1]

ρ[∂tu+(u⋅∇)u]=−∇P+μΔu+ρb.\rho\left[\partial_t\mathbf u+(\mathbf u\cdot\nabla)\mathbf u\right] =-\nabla P+\mu\Delta\mathbf u+\rho\mathbf b.

왼쪽은 단위 부피의 질량에 가속도를 곱한 것이고, 오른쪽은 압력·점성·외부 힘의 합이다. ρ는 밀도, P는 물리적 압력, μ는 점성계수, b는 단위 질량당 외력이다.

밀도가 일정하므로 식 전체를 ρ로 나누고 p = P/ρ, ν = μ/ρ, f = b로 놓으면 아래 형태가 된다.

∂u∂t+(u ⁣⋅ ⁣∇)u=−∇p+νΔu+f,∇ ⁣⋅ ⁣u=0.\begin{aligned} \frac{\partial\mathbf u}{\partial t}+(\mathbf u\!\cdot\!\nabla)\mathbf u &=-\nabla p+\nu\Delta\mathbf u+\mathbf f,\\[0.65em] \nabla\!\cdot\!\mathbf u&=0. \end{aligned}

이것이 여기서 사용할 나비에–스토크스 방정식이다. 원 논문의 정리 1.1은 외력 항을 우변에 두고 같은 운동량 방정식을 사용한다. [9] 밀도를 흡수한 표기와 모든 변수를 무차원화한 것은 다르다. 이 단계에서는 물리 단위가 남아 있다. [1]

기호·항의미단위 또는 읽는 방법
u속도장m/s
∂ₜu한 위치에서 시간에 따라 바뀌는 속도m/s²
(u·∇)u유체가 다른 위치로 이동하면서 겪는 속도 변화m/s²
p = P/ρ밀도로 나눈 압력m²/s²; 물리적 압력 P의 Pa와 구별
−∇p압력 변화에 따른 가속도m/s²
νΔu점성이 속도 차이를 퍼뜨리는 효과ν의 단위는 m²/s
f단위 질량당 외력m/s²

삼각형을 뒤집은 ∇는 공간 미분을 모은 기호다. ∇p는 압력이 공간에서 어떻게 변하는지, ∇·u는 작은 부피에서 순유출이 있는지, Δu는 각 속도 성분의 공간 2차 미분 합을 나타낸다.

∇p=[∂xp∂yp∂zp],Δux=∂xxux+∂yyux+∂zzux.\nabla p=\begin{bmatrix}\partial_xp\\\partial_yp\\\partial_zp\end{bmatrix},\qquad \Delta u_x=\partial_{xx}u_x+\partial_{yy}u_x+\partial_{zz}u_x.

이 수식의 각 항이 모두 가속도 단위로 맞아야 한다. 방정식의 크기를 비교하기 전 단위를 확인하는 습관은 어떤 공학 논문에서도 유용하다.

MOMENTUM BALANCE

흐름을 만드는 다섯 항

∂u∂t+(u ⁣⋅ ⁣∇)u=−∇p+νΔu+f\frac{\partial\mathbf u}{\partial t}+(\mathbf u\!\cdot\!\nabla)\mathbf u=-\nabla p+\nu\Delta\mathbf u+\mathbf f

비압축성 조건 ∇ ⁣⋅ ⁣u=0\nabla\!\cdot\!\mathbf u=0

시간 변화
∂u∂t\frac{\partial\mathbf u}{\partial t}

한 위치에서 관찰한 속도의 시간 변화

이동에 따른 변화
(u ⁣⋅ ⁣∇)u(\mathbf u\!\cdot\!\nabla)\mathbf u

유체 입자가 이동하며 겪는 가속도

압력
−∇p-\nabla p

압력 차이가 만드는 가속도

점성
νΔu\nu\Delta\mathbf u

속도 차이를 퍼뜨리는 효과

외력
f\mathbf f

단위 질량에 외부에서 가하는 힘

방정식을 쓰는 것과 모든 해의 매끄러움을 보장하는 것은 다른 문제다.

#3. 시간에 따라 변하지 않는 속도장에서도 유체는 가속한다

가장 낯선 항은 (u·∇)u일 것이다. 이것을 ‘어려운 비선형 항’이라고 외우기보다, 움직이는 물덩어리를 따라가 보자.

그 물덩어리의 위치를 X(t)라고 하면 속도는 u(X(t),t)다. 연쇄법칙으로 미분하면 다음과 같다.

ddtu(X(t),t)=∂tu+∑j=13dXjdt∂xju=∂tu+(u⋅∇)u.\frac{d}{dt}\mathbf u(\mathbf X(t),t) =\partial_t\mathbf u+\sum_{j=1}^{3}\frac{dX_j}{dt}\partial_{x_j}\mathbf u =\partial_t\mathbf u+(\mathbf u\cdot\nabla)\mathbf u.

위치 변화율 dX/dt가 바로 유체 속도 u이기 때문이다. 이 전체가 물질미분(material derivative)이다. 고정된 카메라가 보는 변화와 움직이는 물덩어리가 경험하는 변화가 다른 것이다.

이를 확인하는 국소적인 교육용 속도장을 놓아 보자.

u(x,y,z)=(ax,−ay,0),a>0.\mathbf u(x,y,z)=(ax,-ay,0),\qquad a>0.

시간 t가 없으므로 ∂ₜu = 0이다. 하지만 x축 방향으로 이동하는 물덩어리는 x가 더 큰 곳으로 가면서 더 큰 속도 ax를 경험한다.

[(u⋅∇)u]x=(ax)∂x(ax)+(−ay)∂y(ax)=a2x.[(\mathbf u\cdot\nabla)\mathbf u]_x=(ax)\partial_x(ax)+(-ay)\partial_y(ax)=a^2x.

시간에 따라 변하지 않는 속도장에서도 물덩어리는 가속할 수 있다. 다만 이 선형 속도장을 무한한 공간 전체에 놓으면 유한 에너지 조건을 만족하지 않는다. 여기서는 미분 항의 의미만 설명한 것이지, 난제의 해를 만든 것이 아니다.

#4. 비압축성은 ‘작아지는 구조가 없다’는 뜻이 아니다

비압축성 조건을 풀어 쓰면 다음과 같다.

∇⋅u=∂xux+∂yuy+∂zuz=0.\nabla\cdot\mathbf u=\partial_xu_x+\partial_yu_y+\partial_zu_z=0.

아주 작은 부피에서 한 방향으로 나가는 양과 다른 방향으로 들어오는 양이 균형을 이룬다는 뜻이다. 어떤 한 방향으로 수축하는 것 자체를 금지하지 않는다.

원통좌표에서 θ 방향 변화가 없는 흐름을 생각하면 조건은 더 직관적이다.

1r∂r(rur)+∂zuz=0.\frac{1}{r}\partial_r(ru_r)+\partial_zu_z=0.

국소적인 설명 예시로 uᵣ = −αr, u_z = 2αz를 넣자. 중심으로 들어오는 흐름과 위아래로 나가는 흐름이다.

1r∂r(−αr2)+∂z(2αz)=−2α+2α=0.\frac{1}{r}\partial_r(-\alpha r^2)+\partial_z(2\alpha z) =-2\alpha+2\alpha=0.

반지름이 작아지는 만큼 축 방향으로 늘어나 부피를 보존할 수 있다. 실제로 같은 유체 입자들을 따라가는 물질 영역 Ω(t)는 아래 조건을 만족한다.

ddtVol⁡(Ω(t))=∫Ω(t)∇⋅u dV=0.\frac{d}{dt}\operatorname{Vol}(\Omega(t)) =\int_{\Omega(t)}\nabla\cdot\mathbf u\,dV=0.

반면 그림에 표시된 ‘빠른 흐름이 집중되는 영역’은 같은 입자들의 집합일 필요가 없다. 물질 영역의 부피와 고속 영역의 크기를 구별해야 한다. 고속 영역이 작아진다는 말은 물 자체를 압축했다는 말이 아니다.

#5. Figure 1을 먼저 정확하게 읽자

OpenAI 원고 Figure 1: 수축하는 중심 코어와 시간에 따른 세장비 변화
원 논문 Figure 1 시각이 진행될수록 코어의 반지름과 높이가 모두 줄고 반지름이 더 빠르게 줄어든다. 수축률의 차이는 도식에서 과장됐으며 실제 실험의 측정 영상이 아니다. OpenAI · Finite Time Blowup for Navier–Stokes · Figure 1 · 출처 · 원문 Figure 인용 · 그림과 캡션을 함께 발췌 · 내용 변경 없음

그림 위의 시각은 t₁ < t₂ < t₃ < 1 순서다. 왼쪽에서 오른쪽으로 특이점 발생 시각에 접근한다. 여기서 1은 해당 구성에서 택한 시간 기준이지, 실제 물이 1초 만에 그렇게 된다는 뜻이 아니다.

캡션의 핵심은 세 가지다. 유체는 안쪽으로 말려 들어오고, z = 0 가까이에 있는 분리층의 양쪽에서 축 방향으로 흐른다. 시간이 흐르면서 강한 흐름의 영역은 작아지고 속도는 증가한다. 그리고 반지름이 높이보다 더 빨리 줄어든다.

따라서 그림을 ‘높이는 실제로 늘어나고 굵기만 줄어든다’고 읽으면 부정확하다. 높이도 줄지만 반지름이 더 빨리 줄어, 상대적인 세장비가 커지는 것이다. 캡션은 이 차이를 과장해 그렸다고 명시한다. 그림의 픽셀 길이로 수축 지수를 추정해서도 안 된다.

OpenAI의 국소 소용돌이 시각화. 안쪽으로 회전하고 축 방향으로 흐르는 구성을 나타낸다.
공식 발표의 소용돌이 시각화. 색상은 각속도를 나타내며, 아래의 수축 스케일을 읽기 위한 도식이다. 실제 특이점을 실험으로 관측한 사진이 아니다.OpenAI · 공식 연구 발표 · 원문과 그림 설명

공식 발표는 색을 각속도와 연결한다. 그런데 회전 반지름이 r, 각속도가 Ω이면 원주 방향 속도는 uθ = rΩ이다. 각속도가 빠른 곳이 반드시 선속도도 가장 빠른 곳은 아니다. 그림은 색상 범례가 무엇의 크기를 표시하는지부터 읽어야 한다. [2]

#6. 전체 에너지와 한 점의 속도는 다른 정보를 담는다

밀도를 1로 정규화한 운동에너지는 다음과 같다.

E(t)=12∫Ω∣u(x,t)∣2 dx=12∥u(t)∥L22.E(t)=\frac12\int_{\Omega}|\mathbf u(\mathbf x,t)|^2\,d\mathbf x =\frac12\|\mathbf u(t)\|_{L^2}^2.

적분은 공간 전체를 합산한다. 반면 최대 속도를 나타내는 L∞ 노름은 가장 큰 값 쪽을 본다.

∥u(t)∥L∞=ess sup⁡x∈Ω∣u(x,t)∣.\|\mathbf u(t)\|_{L^\infty}=\operatorname*{ess\,sup}_{\mathbf x\in\Omega}|\mathbf u(\mathbf x,t)|.

ess sup는 측도 0인 예외를 무시하는 상한이다. 여기서 다루는 매끄러운 속도장에서는 ‘어느 위치에서 가장 큰가’를 보는 상한으로 이해하면 된다. 무한한 영역에서는 최대값이 실제 한 점에서 달성되지 않을 수 있어 max보다 sup가 정확하다.

여러 값을 합산하는 정보와 가장 높은 값을 보는 정보는 다르다. 따라서 유한한 E만으로 최대 속도의 상한을 얻었다고 말할 수 없다. [1]

이를 직접 보이는 수학 예시를 만들 수 있다. 매끄럽고 유한한 영역 안에서만 0이 아닌 벡터장 Φ를 고르고, ∇·Φ = 0이라고 하자. ℓ > 0에 대해 다음과 같이 크기를 조절한다.

uℓ(x)=ℓ−1Φ(x/ℓ).\mathbf u_\ell(\mathbf x)=\ell^{-1}\Phi(\mathbf x/\ell).

ℓ이 줄면 공간 폭은 좁아지고 높이는 커진다. 최대 속도는 다음처럼 변한다.

∥uℓ∥L∞=ℓ−1∥Φ∥L∞⟶∞.\|\mathbf u_\ell\|_{L^\infty}=\ell^{-1}\|\Phi\|_{L^\infty}\longrightarrow\infty.

반면 y = x/ℓ로 적분 변수를 바꾸면 3차원 부피요소는 dx = ℓ³dy가 된다.

Eℓ=12∫ℓ−2∣Φ(x/ℓ)∣2dx=ℓ2∫∣Φ(y)∣2dy⟶0.E_\ell=\frac12\int\ell^{-2}|\Phi(\mathbf x/\ell)|^2d\mathbf x =\frac{\ell}{2}\int|\Phi(\mathbf y)|^2d\mathbf y\longrightarrow0.

이 예시는 최대 속도 증가와 에너지 감소가 동시에 가능하다는 것을 보여준다. 각 ℓ의 속도장은 매끄럽고 발산이 0일 수도 있다. 하지만 이것은 아직 시간에 따라 나비에–스토크스를 만족하는 해가 아니다. ℓ을 시간의 함수로 만들었을 때 필요한 외력이 매끄러운지는 전혀 다른 문제다. 바로 그 간격이 증명에서 어렵다.

#7. 에너지식을 유도하면 점성이 하는 일이 드러난다

다음 계산은 t < T에서 충분히 매끄러운 해를 가정한다. 공간은 무한원에서 경계항이 사라질 만큼 감쇠하는 경우 또는 주기적 경계를 생각한다. 임의의 입출구가 있는 관에 경계항 없이 그대로 적용해서는 안 된다.

운동량 방정식 양쪽에 u를 내적하고 공간에서 적분한다. 시간 미분 항은 운동에너지의 변화율이 된다. 대류 항은 다음 항등식으로 정리된다.

u⋅[(u⋅∇)u]=∇⋅(12∣u∣2u)−12∣u∣2∇⋅u.\mathbf u\cdot[(\mathbf u\cdot\nabla)\mathbf u] =\nabla\cdot\left(\frac12|\mathbf u|^2\mathbf u\right) -\frac12|\mathbf u|^2\nabla\cdot\mathbf u.

발산이 0이고 경계항이 사라지므로 전체 적분은 0이다. 압력 항도 적분하면 사라진다. 점성 항은 부분적분으로 음수가 된다.

ν∫u⋅Δu dx=−ν∫∣∇u∣2dx.\nu\int\mathbf u\cdot\Delta\mathbf u\,d\mathbf x =-\nu\int|\nabla\mathbf u|^2d\mathbf x.

결국 다음을 얻는다.

dEdt+ν∫∣∇u∣2dx=∫f⋅u dx.\boxed{\frac{dE}{dt}+\nu\int|\nabla\mathbf u|^2d\mathbf x =\int\mathbf f\cdot\mathbf u\,d\mathbf x}.

이는 위 지배방정식에서 직접 유도한 에너지 수지다. 왼쪽 두 번째 항은 점성에 의한 운동에너지 소산이고, 오른쪽은 외력이 하는 일이다. 외력이 있으면 운동에너지가 늘어날 수도 있으므로 항상 에너지가 보존된다고 말하면 안 된다. 외력이 없을 때에는 이 조건 아래 운동에너지가 증가하지 않는다.

중요한 점은 대류가 전체 운동에너지를 만들어 내지 않는다는 것이다. 그러면서도 에너지가 공간에서 어떻게 배치되는지에는 관여한다. 전체 합계를 제어하는 것과 국소 집중을 막는 것은 다르다.

#8. 왜 3차원 소용돌이는 더 어려운가

속도장에 회전 연산을 적용한 와도(vorticity)를 ω = ∇×u로 정의한다. 속도 자체가 아니라 국소 회전의 방향과 강도를 나타내는 양이다. 상수 밀도·비압축성 조건에서 운동량 방정식의 curl을 취하면 다음 식을 얻는다. [7]

∂tω+(u⋅∇)ω=(ω⋅∇)u+νΔω+∇×f.\partial_t\boldsymbol\omega+(\mathbf u\cdot\nabla)\boldsymbol\omega =(\boldsymbol\omega\cdot\nabla)\mathbf u+\nu\Delta\boldsymbol\omega+\nabla\times\mathbf f.

오른쪽 첫 번째 항은 와도가 속도 기울기에 의해 늘어나거나 방향을 바꾸는 효과다. 와류 신장(vortex stretching)을 포함한다. 점성이 변화를 퍼뜨리는 동안 흐름의 변형은 특정 방향의 와도를 강화할 수 있다.

평면 2차원 흐름에서는 u_z = 0이고 z방향으로 변하지 않는다. 와도는 z축 방향인데 ∂zu = 0이므로 이 신장 항이 사라진다. 3차원에서는 일반적으로 사라지지 않는다. 이것이 복잡해지는 한 원인을 수식으로 보여준다. 그러나 신장 항이 있다는 것만으로 반드시 특이점이 생긴다는 결론은 나오지 않는다.

보완해서, 와도에 대한 제곱적분을 만들면 신장 항은 부호가 정해지지 않은 채 남는다.

12ddt∥ω∥22+ν∥∇ω∥22=∫ω⋅[(ω⋅∇)u]dx+∫ω⋅(∇×f)dx.\frac12\frac{d}{dt}\|\boldsymbol\omega\|_2^2+\nu\|\nabla\boldsymbol\omega\|_2^2 =\int\boldsymbol\omega\cdot[(\boldsymbol\omega\cdot\nabla)\mathbf u]d\mathbf x +\int\boldsymbol\omega\cdot(\nabla\times\mathbf f)d\mathbf x.

‘점성이 있으니 무조건 괜찮다’는 직관만으로 끝나지 않는 이유다. 이 수식은 방금 얻은 와도 방정식에서, 같은 정칙성과 경계 가정 아래 직접 유도한 것이다.

#9. 공식 코드의 좌표를 읽으면 그림의 수축률을 계산할 수 있다

여기부터는 일반 배경을 넘어 공식 형식화의 실제 정의를 읽는다. 원 논문 §3.1의 식 (3.2)는 다음 관계를 정의한다. 공식 형식화의 CoordinateAlgebra.lean과 SimilarityCoordinates.lean에도 같은 좌표 관계가 기록되어 있다. [3][4][9]

A=12+h,D=12−h,τ=1−t.A=\frac12+h,\qquad D=\frac12-h,\qquad\tau=1-t.
τ=q(1−η2),z=qDη,X=r22q.\tau=q(1-\eta^2),\qquad z=q^D\eta,\qquad X=\frac{r^2}{2q}.

τ는 특이점까지 남은 시간, r은 축에서의 거리, z는 축 방향 좌표다. q는 시간과 축 위치를 함께 반영하는 양의 스케일 변수다. X와 η는 수축을 보정한 좌표다. 여기의 D는 미분 연산자가 아니라 지수다.

η를 없애면 다음 식을 얻는다.

τ=q−z2q2h.\tau=q-z^2q^{2h}.

따라서 q를 모든 위치에서 무조건 τ로 놓으면 안 된다. z = 0에서는 q = τ이지만, 일반적으로는 위 관계를 풀어야 한다. 공식 소스는 이 양의 좌표가 존재하고 유일하다는 성질을 따로 다룬다. [4]

왜 복잡한 좌표를 만들까? 매 순간 작아지는 구조를 원래 좌표에서 보면 형태까지 계속 변한다. 구조와 함께 축척을 조정하면 그 모양을 고정된 범위의 함수로 다룰 수 있다. 이것이 이 좌표를 읽는 직관이다.

고정된 코어 구간에서 |η| ≤ η₀ < 1이라고 하자. 그러면 q는 τ와 상수배 범위로 비슷하다. X도 고정된 유한 범위에서 보면 r = √(2qX)이므로 아래 수축률이 나온다.

R(τ)∝τ1/2,H(τ)∝τ1/2−h,HR∝τ−h.R(\tau)\propto\tau^{1/2},\qquad H(\tau)\propto\tau^{1/2-h},\qquad \frac{H}{R}\propto\tau^{-h}.

h가 작고 양수이면 두 길이 모두 0으로 향하지만 세장비 H/R은 커진다. 이것이 전체 크기는 작아지면서 상대적으로 가늘어진다는 캡션의 내용을 수식으로 읽은 결과다. η가 ±1에 접근하는 끝부분까지 이 단순 비교를 그대로 확장해서는 안 된다.

원 논문 §3.1은 0 < h < 1/100을 두고 이후 작게 선택한다. NaturalAxisData.lean의 국소 구성에는 0 < h ≤ 1/1000, 0 < j ≤ 1/1000이라는 더 구체적인 작은 매개변수 범위가 기록되어 있다. 이는 국소 모듈의 조건이며, 그것만으로 전체 증명이 끝난다는 뜻은 아니다. 이처럼 작은 지수가 쓰이므로 그림의 세장비 변화는 특히 정량 축척으로 읽으면 안 된다. [5]

#10. 속도는 정확히 어떤 성분에서 커지는가

공식 코드 NaturalCore.lean의 coreVelocity_swirl은 원주 방향 속도를 아래처럼 표현한다. 코드에서 프로필 이름으로 쓰는 f는 외력 f와 혼동하지 않도록 여기서 F로 바꿔 썼다. [6]

uθ=q−A2X F(X,η).u_\theta=q^{-A}\sqrt{2X}\,F(X,\eta).

F는 수축을 보정한 좌표에서의 모양을 나타낸다. 핵심은 그 앞에 붙은 q⁻ᴬ다. A > 0이므로 q가 줄면 속도의 대표 크기는 커진다.

하지만 회전축에서는 원주 방향 기저가 특이해 보일 수 있다. 좌표가 특이해 보이는 것과 실제 벡터장이 특이한 것은 다르다. 위 식을 r로 쓰면 다음과 같다.

uθ=rq−(1+h)F(X,η).u_\theta=rq^{-(1+h)}F(X,\eta).

고정된 t < 1에서 원점으로 접근하면 q는 양수이고 r은 0으로 간다. 프로필이 매끄러우면 축에서 회전 속도가 처음부터 무한한 것이 아니다.

오히려 coreVelocity_at_origin은 원점의 축 방향 속도를 직접 보여준다. [6]

ucore(0,t)=j(1−t)−(1/2+h)ez.\mathbf u_{\mathrm{core}}(\mathbf 0,t)=j(1-t)^{-(1/2+h)}\mathbf e_z.

j는 양수이고 e_z는 축 방향 단위벡터다. 따라서 t가 1보다 작은 동안에는 유한하지만 1로 접근할수록 그 크기에 공통 상한을 둘 수 없다. 이 국소 코어와 최종 전체 흐름을 연결하는 일은 별도로 필요하다. 코드 역시 국소 코어 존재와 최종 외력의 정칙성을 구분한다.

이제 에너지와 연결해 보자. 고정된 재스케일 코어에서 프로필이 제어되고 대표 속도가 U에 비례한다고 가정한 차수 추정은 다음과 같다.

Ecore≈U2R2H∝τ−2(1/2+h)τ1τ1/2−h=τ1/2−3h.E_{\mathrm{core}}\approx U^2R^2H \propto\tau^{-2(1/2+h)}\tau^{1}\tau^{1/2-h} =\tau^{1/2-3h}.

h < 1/6이면 지수가 양수다. 따라서 이 스케일만 보면 속도는 증가하면서 코어의 에너지는 감소할 수 있다. 공식 국소 매개변수 범위는 이 조건 안에 있다.

이 계산은 에너지와 속도 스케일의 양립 가능성을 설명한다. 전체 해의 에너지 유계성을 증명한 것은 아니다. 코어 밖의 흐름, 보정항, 미분, 전역 적분까지 별도로 통제해야 한다. 단순한 차수 계산과 완성된 정리를 혼동하지 않는 것이 논문을 읽는 핵심이다.

추가로 U R/ν라는 국소 레이놀즈 수를 같은 스케일로 계산하면 τ⁻ʰ/ν에 비례한다. 작은 영역으로 줄어든다고 해서 이 비교량이 반드시 작아지는 것은 아니다. 속도가 함께 증가하기 때문이다. 이 역시 스케일 관찰이며 점성을 이겼다는 완전한 증명은 아니다.

#11. 진짜 난점은 ‘빠른 속도’가 아니라 ‘매끄럽게 남는 외력’이다

속도장 U와 압력 p₀를 임의로 제안했다고 하자. 이것을 방정식에 넣고 남는 값을 잔차(residual)라고 정의한다.

R(U,p0)=∂tU+(U⋅∇)U+∇p0−νΔU.\mathcal R(\mathbf U,p_0)=\partial_t\mathbf U+(\mathbf U\cdot\nabla)\mathbf U +\nabla p_0-\nu\Delta\mathbf U.

f를 이 잔차로 정하면 식은 형식적으로 맞는다. 하지만 그렇게 만든 f가 매끄럽다는 보장은 없다. U가 끝없이 커지면서 그 미분도 커지면 f가 같이 발산할 수 있다. 무한한 힘을 넣고 속도가 발산했다고 말하는 것은 문제의 조건을 만족시키지 못한다.

심지어 f의 값만 유한해도 부족하다. 미분이 망가지는 힘은 매끄러운 외력이 아니다. 공식 발표는 가속도·운동량 전달·압력·점성 항들이 커지면서도 정밀하게 상쇄되어 매끄러운 외력을 남기는 것이 핵심 과제라고 설명한다. [2]

그 구조를 읽기 위한 일반적인 보정식을 직접 전개해 보자. u = U + w, p = p₀ + π로 놓으면 다음이 성립한다.

R(U+w,p0+π)=R(U,p0)+∂tw+(U⋅∇)w+(w⋅∇)U+(w⋅∇)w+∇π−νΔw.\begin{aligned} \mathcal R(\mathbf U+\mathbf w,p_0+\pi) ={}&\mathcal R(\mathbf U,p_0)+\partial_t\mathbf w +(\mathbf U\cdot\nabla)\mathbf w+(\mathbf w\cdot\nabla)\mathbf U\\ &+(\mathbf w\cdot\nabla)\mathbf w+\nabla\pi-\nu\Delta\mathbf w. \end{aligned}

이 식은 해당 원고의 식 번호를 인용한 것이 아니라, 지배방정식을 대수적으로 전개한 교육용 설명이다. 보정 흐름을 추가하면 좋은 항 하나만 생기는 것이 아니라 상호작용, 시간 변화, 점성 비용이 함께 생긴다.

∇·w = 0이면 보정항의 비선형 부분은 다음처럼 쓸 수 있다.

(w⋅∇)w=∇⋅(w⊗w).(\mathbf w\cdot\nabla)\mathbf w=\nabla\cdot(\mathbf w\otimes\mathbf w).

⊗는 외적 곱으로, i,j 성분이 wᵢwⱼ인 행렬을 만든다. 여기서 벡터의 cross product와는 다르다. 이런 2차 항 때문에 평균이 0인 작은 진동도 평균 흐름에 영향을 줄 수 있다.

이를 설명하는 스칼라 예시로 한 주기 평균을 보면 다음과 같다.

⟨asin⁡θ⟩=0,⟨a2sin⁡2θ⟩=a22.\langle a\sin\theta\rangle=0,\qquad \langle a^2\sin^2\theta\rangle=\frac{a^2}{2}.

위아래로 진동해 평균은 0이어도 제곱의 평균은 0이 아니다. 이것이 진동과 평균 응력의 연결을 이해하는 출발점이다. 하지만 이 사인함수를 곧바로 해당 연구의 3차원 보정해라고 부를 수는 없다. 공간 의존성, 발산 조건, 경계, 압력, 점성까지 맞아야 한다.

특히 주파수가 큰 진동은 공간 2차 미분도 커져 점성 항의 비용을 키운다. ‘빠르게 흔들면 상쇄된다’는 직관만으로는 끝나지 않는다. 만들어 낸 모든 항을 같은 식 안에서 통제하는 것이 어려운 부분이다.

#12. 무엇을 증명했다고 주장하는지와 무엇을 읽어냈는지

Clay의 공식 문제는 A·B에서 외력이 없는 경우를, C·D에서 허용된 매끄러운 외력을 포함하는 붕괴 구성을 다룬다. 공식 형식화는 C·D에 대응하는 결과를 기록한다. 이번 구성의 설명을 외력 없는 경우의 해법으로 바꾸어 말하면 안 된다. [1][8]

여기까지의 수식은 지배방정식, 좌표 정의와 국소 코어의 속도를 연결한다. 에너지 수지와 보정식의 대수 전개가 맞더라도, 전체 해의 존재와 모든 경계·정칙성 조건은 별도로 증명해야 한다.

공식 메타데이터의 sorry_count: 0은 주요 형식화에 미완성 증명 표식을 사용하지 않았다는 자체 보고다. 검토 상태는 self-assessed로 기재돼 있다. 어떤 정리와 가정을 기계가 검사했는지, 원래 문제와 같은 정의를 사용하는지에 대한 독립 검토도 필요하다. [8]

유도한 식 하나가 맞는 것, 국소 코어가 존재하는 것, 전체 흐름을 구성하는 것, 최종 문제의 조건을 모두 충족하는 것은 구별해야 한다. 특히 NaturalCore.lean 자체도 최종 외력의 정칙성을 이 모듈에서 주장하지 않는다고 명시한다. [6]

수학 주장 공개부터 형식화, 독립 검토, 공식 인정까지의 네 단계
자체 제작 설명 도식형식 검증 자료는 중요한 증거지만 독립 검토와 공식 인정까지 자동으로 대신하지 않는다.JJo · 기사에서 확인한 원문을 바탕으로 재구성 · 실제 측정 Figure가 아님

#13. 이해를 확인하는 세 질문

첫째. 속도장에 t가 없으면 가속도도 0인가? 아니다. 물질미분에는 위치를 이동하면서 생기는 변화가 포함된다. 앞의 예시에서 ∂ₜu = 0이지만 x방향 가속도는 a²x였다.

둘째. 고속 영역의 크기가 줄면 비압축성과 모순인가? 아니다. 고속 영역과 같은 유체 입자들을 따라가는 물질 영역은 다르다. 같은 입자들의 부피 보존은 별도의 조건이다.

셋째. 에너지 차수 τ^(1/2−3h)가 양수이면 증명이 끝나는가? 아니다. 코어 밖의 흐름과 보정항, 외력의 모든 필요한 미분, 전체 공간에서의 조건이 남아 있다. 이 계산은 스케일의 양립성을 보여줄 뿐이다.

이 질문에 답할 수 있다면 ‘AI가 어려운 문제를 풀었다’는 제목을 넘어, 논문이 실제로 넘어야 하는 장벽을 설명할 수 있다.

#14. 핵심을 한 문장으로

이 연구를 이해하는 핵심은 ‘소용돌이가 빨라진다’가 아니라, 수축하는 공간 스케일과 증가하는 속도 스케일이 유한 에너지와 양립하면서도 나비에–스토크스 잔차가 매끄러운 외력으로 남아야 한다는 것이다.

#출처

[1] Charles L. Fefferman / Clay Mathematics Institute, 공식 문제 정의, pp. 1–2, 식 (1)–(7)과 조건 A–D.

[2] OpenAI, On the Navier–Stokes Millennium Prize Problem, 2026-09-08. 공식 발표의 색상 설명과 정밀 상쇄에 관한 서술. 기업의 연구 주장과 독립 검증을 구분한다.

[3] OpenAI, CoordinateAlgebra.lean, A, D, forward_coordinate_identity와 좌표 관련 주석.

[4] OpenAI, SimilarityCoordinates.lean, 양의 q 좌표 정의와 존재·유일성 정리.

[5] OpenAI, NaturalAxisData.lean, SmallParameters, A_bounds, D_bounds.

[6] OpenAI, NaturalCore.lean, similarityPoint, coreVelocity_swirl, coreVelocity_at_origin, extension_speedUnbounded 및 모듈 설명. 자연 코어의 방사 좌표는 AxisymmetricFields.lean의 radialEnergy를 따른다.

[7] MIT OpenCourseWare, Vorticity Theorem for a Viscous Fluid, 1.63, Fall 2002. 와도 방정식의 배경 강의 자료.

[8] OpenAI, formalization.yaml. 원고와 형식화의 대응 및 자체 보고 검증 상태.

[9] OpenAI, Finite Time Blowup for Navier–Stokes, 정리 1.1, §2.1–2.4, §3.1 식 (3.1)–(3.3), §3.3. 2026-09-08 공개.

관련 글: AI 연구 자동화의 지표 읽기 · NEWS

Connect