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

この論文は、クリーン数値シミュレーションを用いた計算結果から、初期条件の差が極めて微小（$10^{-40}$ 程度）であってもナビエ - ストークス方程式が異なる大域解を持つ可能性を示唆し、クレイ数学研究所のミレニアム懸賞問題に関連する解の一意性と存在性について新たな洞察を提供するものである。

要約

この論文は、流体力学の「聖杯」とも呼ばれるナビエ・ストークス方程式（空気や水の流れを記述する方程式）に関する、非常に驚くべき発見を報告したものです。

専門用語を避け、日常の例え話を使って、何が起きたのかをわかりやすく解説します。

1. 核心となる発見：「ほんの少しの違い」が「全く違う未来」を作る

この研究の結論は一言で言うと、**「ナノメートル単位（髪の毛の直径の何万分の 1）の差しかない、ほぼ同じ状態から出発しても、時間が経つと『全く異なる』流れ方になる」**という証拠を、スーパーコンピューターで見つけたということです。

通常、私たちは「同じような条件なら、結果も同じはずだ」と考えがちです。しかし、この研究は**「ナビエ・ストークス方程式の世界では、それは通用しない」**と示唆しています。

2. 使われた「魔法の道具」：クリーン・シミュレーション（CNS）

なぜこれまでこの発見ができなかったのでしょうか？それは、従来のコンピューター計算には「ノイズ（計算誤差）」がつきものだからです。

• 従来の計算（DNS）：
  想像してください。非常に敏感なバランスの取れた塔を積み上げようとしていますが、指先が少し震えていて、微細な揺れ（計算誤差）が常に起きている状態です。最初は安定して見えても、時間が経つにつれてその揺れが増幅され、塔は崩れてしまいます。
  従来の計算では、この「計算の揺れ」が、本当の流れ（物理現象）よりも大きくなってしまい、結果が信用できなくなってしまうのです。

• 今回の計算（CNS）：
  著者たちは**「クリーン・シミュレーション（CNS）」という新しい方法を使いました。これは、「計算の揺れを極限まで抑え込んだ、超精密な計算」です。
  例えるなら、「完全な静寂の中で、微動だにしない手つきで、何千時間もの間、塔を積み上げ続ける」**ようなものです。これにより、計算の誤差が本物の物理現象に混入するのを防ぎ、非常に長い時間（300 秒間）にわたって、本当の流れに極めて近い軌跡を追うことができました。

3. 実験の内容：「蝶の羽ばたき」の証明

実験では、2 次元の「コルモゴロフ流」と呼ばれる乱流（カオス的な流れ）をシミュレーションしました。

• 条件 A（基準）： 完全に対称な形をした初期状態。
• 条件 B（実験）： 条件 A に、**「10 の -40 乗」**という、もはや実在しないような微小な歪み（$\delta$）を加えた状態。
  • 例え話: 地球全体の重さに対して、**「砂粒 1 粒」**の重さ分だけ違う状態、いや、それよりもっと小さい違いです。

結果：
最初は、条件 A と条件 B の流れは全く同じように見えました。しかし、時間が経つにつれて、条件 B の「微小な歪み」が**「バタフライ効果（蝶の羽ばたきが嵐を引き起こす現象）」**のように急激に増幅されました。

• 条件 A の結果： 回転対称性を持った、ある特定のパターンで落ち着きました。
• 条件 B の結果： 対称性が崩れ、全く異なるパターンに変わってしまいました。

驚くべき点は：
この「微小な歪み」は、計算の誤差ではなく、**意図的に加えた「物理的な初期条件の違い」です。つまり、「同じ方程式から、ほぼ同じ出発点でも、2 つの全く異なる未来（解）が生まれる」**ことが数値的に証明されたのです。

4. なぜこれが重要なのか？

この発見は、数学界の「ミレニアム懸賞問題（100 万ドルの賞金がかかる難問）」の一つである**「ナビエ・ストークス方程式の解の一意性（解は一つだけか？）」**に大きな疑問を投げかけています。

• これまでの常識： 「同じ初期条件なら、解は一つしかないはずだ」と考えられてきた。
• この論文の示唆： 「もし初期条件が『無限に近い精度』で一致しない限り（あるいは、無限に小さな差でも）、解は複数存在する可能性がある」。

著者たちは、この数値実験が、数学者たちが「解の非一意性」を証明するための**「道しるべ（羅針盤）」**になることを期待しています。

まとめ

この論文は、**「超精密な計算機を使って、ナノレベルの『小さな違い』が、乱流の世界では『大きな違い』として現れることを突き止め、同じ方程式から複数の異なる未来が生まれる可能性があることを示した」**という画期的な研究です。

まるで、**「同じように見えた 2 つの川が、源流で砂粒 1 つの差があっただけで、下流では全く異なる海に流れ着く」**ような現象を、計算機の中で鮮明に捉えたようなものです。

技術概要

以下は、提示された論文「Non-uniqueness of smooth solutions of the Navier-Stokes equations from almost the same initial conditions（ナヴィエ - ストークス方程式の滑らかな解の非一意性：ほぼ同一の初期条件から）」に関する詳細な技術的サマリです。

1. 研究の背景と課題 (Problem)

ナヴィエ - ストークス（NS）方程式の滑らかな解の「存在と一意性」は、クレイ数学研究所が定めるミレニアム懸賞問題の一つです。

• 既存の知見: Coiculescu と Palasek は、臨界空間 $BMO^{-1}$ において、同じ初期データから 2 つの異なる滑らかな大域解が存在することを証明しましたが、その構成は無限エネルギーを必要とする簡略化されたメカニズムに依存していました。
• 数値計算の限界: 従来の数値シミュレーション（直接数値シミュレーション：DNS）では、丸め誤差や打ち切り誤差といった「数値ノイズ」が避けられません。カオス系である乱流において、この微小な数値ノイズは「バタフライ効果」により指数関数的に増幅され、短時間で巨視的な誤差となり、真の解からの乖離を引き起こします。そのため、有限エネルギーを持つ初期条件から、微小な違いが異なる大域解を生むかどうかを数値的に証明することは、従来の手法では不可能とされてきました。

2. 手法 (Methodology)

本研究では、従来の DNS の限界を克服するために、クリーン・ナメリカル・シミュレーション（CNS: Clean Numerical Simulation） を採用しました。

• CNS の特徴:
  • 高精度計算: 物理量・数値パラメータすべてに対して多倍長演算（multiple-precision）を使用し、丸め誤差を極めて微小なレベルまで抑制します。
  • 誤差制御: 打ち切り誤差と丸め誤差を必要な精度まで低減させることで、有限時間内において数値ノイズを厳密に無視可能（negligible）にします。
  • 検証プロセス: 異なる精度設定（より高い次数 $M$ と有効桁数 $N_s$）で計算を行い、結果を比較することで、対象とする時間区間内の解が「収束した基準解（benchmark solution）」として真の解に極めて近いことを確認します。
• 数値設定:
  • 方程式: 2 次元非圧縮性コルモゴロフ流（Kolmogorov flow）。
  • 領域: 周期境界条件を持つ正方形領域 $[0, 2\pi] \times [0, 2\pi]$。
  • パラメータ: レイノルズ数 $Re = 2000$、コルモゴロフ強制スケール$n_K = 16$。
  • 離散化: 空間方向はフーリエ擬スペクトル法（グリッド $1024 \times 1024$）、時間方向は $M$ 次テイラー展開（時間ステップ $\Delta t = 10^{-3}$）。
  • 計算リソース: 天河次世代スーパーコンピュータを使用。

3. 実験設定 (Experimental Setup)

ほぼ同一の初期条件から出発し、その微小な差異が時間発展にどう影響するかを調査しました。

• 初期条件 1 ($\psi_1$): 回転対称性と並進対称性の両方を持つ解。
  $$\psi_1(x, y, 0) = -\frac{1}{2}[\cos(x + y) + \cos(x - y)]$$
• 初期条件 2 ($\psi_2$): $\psi_1$ に微小な摂動 $\delta \sin(x + y)$ を加えたもの。
  $$\psi_2(x, y, 0) = \psi_1(x, y, 0) + \delta \sin(x + y)$$
  ここで、$\delta$ として $10^{-10}, 10^{-20}, 10^{-40}$ の 3 つの値を設定しました。
  • 注意: $\psi_2$ は $\delta$ が極めて小さくても、初期状態において $\psi_1$ とは異なる対称性（回転対称性の欠如）を持ちます。

4. 主要な結果 (Key Results)

CNS による計算結果（時間区間 $t \in [0, 300]$）から、以下の重要な発見が得られました。

1. 解の分岐:

   • 初期段階では、摂動を加えた 3 つのケース（Flow CNS'$_1$, CNS'$_2$, CNS'$_3$）は、基準となる Flow CNS とほぼ同じ挙動を示します。
   • しかし、ある時間経過（$t > 100$ 付近）後、これらは Flow CNS から明確に分岐します。
   • 空間平均の運動エネルギー散逸率 $\langle D \rangle_A$ において、$t > 100$ 以降、Flow CNS の値は他の 3 つのケースの約 2 倍の大きさになります。

2. 対称性の破れと多様性:

   • Flow CNS: 初期条件の持つ「回転対称性」と「並進対称性」を維持したまま時間発展します。
   • Flow CNS'$_i$: 微小な摂動 $\delta$ がカオス的な「ノイズ拡大カスケード（noise-expansion cascade）」により巨視的なレベルまで増幅され、元の回転対称性を破壊します。その結果、これらは初期条件 $\psi_2$ と同じ「並進対称性」のみを持つ異なる解となります。
   • 統計的差異: 時間・空間平均された統計量（確率密度関数など）も、Flow CNS と他の 3 つのケースの間で明確な差異を示します。

3. 有限エネルギーの解:

   • 得られたすべての解は有限のエネルギー密度を持ち、滑らかです。

5. 結論と意義 (Significance)

• 非一意性の数値的証拠:
  本研究は、「無限に小さな差（$\delta \to 0$）を持つほぼ同一の初期条件から、ナヴィエ - ストークス方程式が異なる大域滑らかな解（finite energy）を許容する」 ことを数値的に示しました。
• 仮説の提示:
  著者は以下の仮説を提唱しています。

  「任意に小さな $\delta$（あるいは $\delta \to 0$）に対して、$|U_1 - U_2| \le \delta$ となる初期条件 $U_1$ と $U_2$ から、$t > 0$ に対して異なる滑らかな大域解が存在する。」

• 数学的意義:
  これらの結果は、NS 方程式の解の一意性に関する数学的証明への道筋（ロードマップ）や、反例の構築への示唆を与える可能性があります。特に、従来の数値計算では「ノイズによる誤差」として扱われていた現象が、実際には「異なる解への分岐」を意味し得ることを示した点に革新性があります。
• CNS の有効性:
  乱流のようなカオス系において、高精度な数値シミュレーション（CNS）を用いることで、数値ノイズを排除し、真の物理的挙動（解の非一意性）を捉えることが可能であることを実証しました。

まとめ

この論文は、クリーン・ナメリカル・シミュレーション（CNS）という高精度手法を用いることで、ナヴィエ - ストークス方程式において、無限に微小な初期条件の差異が、時間経過とともに巨視的に異なる解（異なる対称性、異なる統計的性質を持つ解）へと成長することを示しました。これは、NS 方程式の解の一意性に関するミレニアム懸賞問題に対し、有限エネルギーの滑らかな解の文脈で新たな視点と数値的証拠を提供する重要な研究です。