Non-uniqueness of smooth solutions of the Navier-Stokes equations from almost the same initial conditions

이 논문은 정밀 수치 시뮬레이션 (CNS) 을 통해 초기 조건의 차이가 $10^{-40}$정도로 미미함에도 불구하고 나비어 - 스토크스 방정식이 서로 다른 전역 해를 가질 수 있음을 수치적으로 증명하여, 클레이 수학연구소의 밀레니엄 문제 중 하나인 나비어 - 스토크스 방정식의 존재성과 유일성에 대한 통찰을 제공합니다.

핵심

이 논문은 수학의 거대한 난제 중 하나인 **'나비에-스토크스 (Navier-Stokes) 방정식'**에 대한 흥미로운 발견을 담고 있습니다. 어렵게 들릴 수 있는 이 내용을 일상적인 비유와 함께 쉽게 설명해 드리겠습니다.

🌊 핵심 주제: "미세한 차이가 거대한 파도를 만든다"

이 연구의 핵심은 **"거의 똑같은 조건에서 시작했는데, 전혀 다른 결과가 나올 수 있다"**는 것입니다.

상상해 보세요. 두 개의 아주 똑같은 컵에 물을 담았습니다. 그런데 한 컵에 원자 하나만큼의 먼지만 더 넣었다고 칩시다. 보통 우리는 "아, 거의 똑같으니 결과도 비슷하겠지?"라고 생각합니다. 하지만 이 연구는 **"아니, 그 아주 작은 먼지 때문에 두 컵의 물결이 완전히 다른 방향으로 흐르게 된다"**고 주장합니다.

🔍 연구의 배경: 왜 기존 컴퓨터로는 못 봤을까?

우리는 컴퓨터로 날씨나 물의 흐름을 시뮬레이션할 때, 아주 작은 오차 (숫자 계산의 반올림 오차 등) 를 피할 수 없습니다.

• 기존 방식 (DNS): 마치 아주 미세한 진동에도 흔들리는 저울 위에 설탕을 올리는 것과 같습니다. 아주 작은 오차가 '나비 효과 (Butterfly Effect)'처럼 증폭되어, 시간이 지나면 계산 결과가 실제 현상과 완전히 달라져버립니다. 그래서 "처음이 비슷했는데 결과가 왜 달라졌지?"라고 해도, "아마 계산 오류 때문이겠지"라고 치부해 왔습니다.

🚀 새로운 방법: '청소 Numerical Simulation (CNS)'

저자들은 기존 방법의 한계를 극복하기 위해 **'청소 (Clean) 수치 시뮬레이션 (CNS)'**이라는 새로운 기술을 개발했습니다.

• 비유: 마치 고층 빌딩을 짓기 전에, 땅의 미세한 진동까지 완벽하게 차단하고 아주 정밀한 자로 측정한 뒤 공사를 시작하는 것과 같습니다.
• 이 방법은 오차를 극도로 줄여서, 오류가 아닌 진짜 물리 법칙이 만들어내는 결과를 아주 오랫동안 정확하게 보여줍니다.

🧪 실험 내용: 10의 -40 제곱의 차이

연구진은 두 가지 아주 비슷한 초기 조건을 설정했습니다.

1. 조건 A: 완벽한 대칭을 가진 물의 흐름.
2. 조건 B: 조건 A 와 거의 똑같지만, **10 의 -40 제곱 (1 뒤에 0 이 40 개)**이라는, 상상할 수 없을 정도로 미세한 수치를 더한 조건.

이 차이는 원자보다도 훨씬 작고, 우리가 눈으로 볼 수 있는 한도보다 훨씬 작습니다.

🌪️ 결과: 완전히 다른 두 가지 세상

시간이 지나자 놀라운 일이 벌어졌습니다.

• 초반: 두 흐름은 거의 똑같았습니다.
• 중반: 아주 미세한 차이가 '나비 효과'로 인해 폭발적으로 커졌습니다.
• 결국: 두 흐름은 완전히 다른 모양이 되었습니다.
  • 조건 A는 대칭적인 아름다운 무늬를 유지했습니다.
  • 조건 B는 그 대칭이 깨진, 완전히 다른 난기류를 만들었습니다.
  • 심지어 통계적 성질 (에너지가 얼마나 소모되는지 등) 도 서로 달랐습니다.

💡 이 발견이 중요한 이유

수학계에는 **"나비에-스토크스 방정식은 초기 조건이 같으면 결과도 반드시 하나뿐이다 (해의 유일성)"**라는 믿음이 있었습니다. 이는 '밀레니엄 문제' 중 하나로, 증명되면 100 만 달러의 상금이 걸려 있습니다.

하지만 이 논문은 **"아니요, 초기 조건이 거의 같아도 (10 의 -40 제곱 차이), 해가 두 개 이상 나올 수 있습니다"**라고 보여주는 강력한 증거를 제시했습니다.

📝 결론: 수학자들에게 보내는 메시지

이 연구는 "우리가 계산한 게 틀린 게 아니라, 자연의 법칙 자체가 이렇게 복잡하고 예측 불가능할 수 있다"는 것을 보여줍니다.

• 비유: 같은 씨앗을 심었는데, 흙 한 알의 차이 때문에 하나는 장미가 되고 다른 하나는 선인장이 된다면, 우리는 그 씨앗이 '유일한' 결과를 내지 않는다고 말할 수 있습니다.

이 논문은 수학자들이 **"왜 해가 유일하지 않은지"**를 증명하는 데 중요한 단서 (로드맵) 를 제공하며, 유체 역학과 수학의 새로운 장을 열 수 있는 가능성을 보여줍니다.

---

한 줄 요약: "원자 하나만큼의 미세한 차이도 시간이 지나면 완전히 다른 세상을 만들 수 있으며, 나비에-스토크스 방정식은 그런 '비유일성'을 가질 수 있다"는 놀라운 발견을 정밀한 계산으로 증명했습니다.

기술 요약

논문 요약: 나비에 - 스토크스 방정식의 거의 동일한 초기 조건에서 발생하는 비유일성

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 난류 (Turbulence) 는 나비에 - 스토크스 (NS) 방정식으로 설명되며, NS 방정식의 매끄러운 해의 존재성과 유일성은 클레이 수학연구소 (Clay Institute) 의 밀레니엄 문제 중 하나입니다.
• 현재 상황: Coiculescu 와 Palasek 은 무한한 에너지를 가진 특수한 공간 (BMO−1) 에서 두 개의 서로 다른 전역 해가 존재할 수 있음을 증명했으나, 유한한 에너지를 가지면서도 명확한 수치적 예시를 찾는 것은 전통적인 수치 방법으로는 불가능했습니다.
• 핵심 문제: 전통적인 직접 수치 시뮬레이션 (DNS) 은 수치적 잡음 (round-off error, truncation error) 을 피할 수 없으며, 난류의 혼돈적 특성 (나비 효과) 으로 인해 이 작은 잡음이 기하급수적으로 증폭되어 정확한 해와 완전히 다른 궤적을 그리게 됩니다. 따라서, 초기 조건이 극도로 유사함에도 불구하고 해가 유일하지 않을 수 있는지에 대한 명확한 수치적 증거를 제시하는 것이 난제였습니다.

2. 방법론 (Methodology)

• 청정 수치 시뮬레이션 (Clean Numerical Simulation, CNS) 도입:
  • 기존 DNS 의 한계를 극복하기 위해 Liao 가 제안한 CNS 기법을 사용했습니다.
  • CNS 는 절단 오차 (truncation error) 와 반올림 오차 (round-off error) 를 극도로 낮추기 위해 다중 정밀도 (multiple-precision) 연산 (예: 260 자리의 유효 숫자) 과 고차 테일러 전개 (M=140 등) 를 적용합니다.
  • 이를 통해 수치 잡음을 엄밀하게 무시할 수 있는 수준으로 낮추어, 유한하지만 통계적으로 충분한 시간 구간 ($t \in [0, 300]$) 동안 NS 방정식의 정확한 시공간 궤적을 얻을 수 있습니다.
• 수치 모델 설정:
  • 방정식: 2 차원 비압축성 콜모고로프 흐름 (Kolmogorov flow) 을 2D 도메인 $[0, 2\pi] \times [0, 2\pi]$에서 고려합니다.
  • 조건: 레이놀즈 수 $Re = 2000$, 콜모고로프 힘 스케일$n_K = 16$.
  • 초기 조건 비교:
    1. Flow CNS: 대칭성을 가진 초기 조건 $\psi_1(x, y, 0) = -\frac{1}{2}[\cos(x+y) + \cos(x-y)]$.
    2. Flow CNS'$_{1,2,3}$: $\psi_1$에 미세한 교란 $\delta \sin(x+y)$를 추가한 조건. 여기서 $\delta$는 $10^{-10}, 10^{-20}, 10^{-40}$로 설정되어 초기 조건이 거의 동일함 ($\delta \to 0$) 을 보여줍니다.
• 검증: 더 작은 수치 잡음을 가진 또 다른 CNS 결과를 '기준 (benchmark)'으로 사용하여, $t \in [0, 300]$ 구간에서 수치 오차가 해에 미치는 영향이 무시할 수 있음을 확인했습니다.

3. 주요 결과 (Key Results)

• 궤적의 분기 (Divergence of Trajectories):
  • 초기에는 모든 시뮬레이션 결과가 거의 동일하게 시작되지만, 시간이 지남에 따라 (특히 $t > 100$ 이후) 미세한 초기 차이 ($\delta$) 가 기하급수적으로 증폭되어 서로 완전히 다른 궤적을 보입니다.
  • 공간 평균 운동 에너지 소산율 ($\langle D \rangle_A$): $t > 100$ 이후, 대칭성을 가진 Flow CNS 의 소산율은 다른 세 가지 흐름 (Flow CNS'$_{1,2,3}$) 보다 약 2 배 더 큰 값을 보입니다.
• 대칭성의 파괴 (Symmetry Breaking):
  • Flow CNS: 초기 조건이 가진 회전 및 병진 대칭성을 유지합니다.
  • Flow CNS'$_{1,2,3}$: 초기에 존재했던 미세한 교란 ($\delta$) 이 '잡음 확장 연쇄 (noise-expansion cascade)'를 통해 거시적 수준으로 증폭되어 원래의 회전 대칭성을 파괴하고, 오직 병진 대칭성만 남깁니다.
  • 이는 $\delta$가 $10^{-40}$처럼 극도로 작아도 최종 해의 거시적 성질 (대칭성 포함) 을 완전히 바꿀 수 있음을 의미합니다.
• 통계적 차이:
  • 서로 다른 초기 조건에서 출발한 해들은 운동 에너지 소산율의 확률 밀도 함수 (PDF) 와 공간 - 시간 평균값에서 뚜렷한 통계적 차이를 보입니다.
• 유한 에너지: 모든 시뮬레이션된 해는 유한한 에너지 밀도를 가지며, 이는 물리적으로 타당한 조건입니다.

4. 주요 기여 (Key Contributions)

• 수치적 비유일성 증거 제시: 거의 동일한 초기 조건 (차이 크기 10^{-40 수준) 에서 출발하여, 유한한 에너지를 가진 매끄러운 전역 해가 서로 다르게 발산할 수 있음을 수치적으로 명확히 증명했습니다.
• CNS 기법의 유효성 입증: 전통적인 DNS 가 놓칠 수 있는 난류의 미세한 초기 민감도를 포착할 수 있는 CNS 의 능력을 보여주었습니다.
• 잡음 확장 연쇄 (Noise-Expansion Cascade) 발견: NS 난류에서 미세한 거시적 차이가 순차적으로 증폭되어 전체 유동의 성질을 근본적으로 변화시킨다는 메커니즘을 확인했습니다.

5. 의의 및 결론 (Significance)

• 밀레니엄 문제에 대한 시사점: 이 연구는 NS 방정식의 해의 유일성에 대한 새로운 관점을 제공합니다. 즉, 초기 조건이 임의로 작은 차이 ($\delta \to 0$) 를 가질지라도, $t > 0$에 대해 서로 다른 매끄러운 전역 해가 존재할 수 있다는 강력한 수치적 반례 (counter-example) 를 제시합니다.
• 수학적 증명에 대한 로드맵: 이 연구는 수학자들이 NS 방정식의 비유일성을 증명하기 위한 구체적인 방향성 (roadmap) 을 제시할 수 있는 영감을 줍니다.
• 결론: 저자들은 "NS 방정식은 임의로 작은 $\delta$에 대해 거의 동일한 초기 조건 $U_1$과 $U_2$ ($|U_1 - U_2| \le \delta$) 에서 서로 다른 매끄러운 전역 해를 허용한다"는 가설을 제시하며, 이는 NS 방정식의 해의 유일성 문제에 대한 중요한 통찰을 제공합니다.

---

참고: 이 연구는 중국 천허 (Tianhe) 슈퍼컴퓨터를 사용하여 수행되었으며, 관련 코드는 GitHub 를 통해 공개되어 있습니다.