圧縮性なしの多孔質媒体、ブシネスク、および圧縮性なしオイラー方程式における滑らかな強制項を用いた有限時間爆発

非圧縮性多孔質媒体、Boussinesq、および非圧縮性Euler方程式における滑らかな強制項を用いた有限時間爆発

AlpögeとBuckmasterによるエキサイティングなごく最近の研究は、以前のCórdovaとMartínez-Zoroaによる研究を基盤として、悪名高い非圧縮性3次元Navier-Stokes方程式のグローバル正則性問題に関する一般的なトピックを扱っています。現在、滑らかな初期データと滑らかな強制項を用意することで、これらの方程式が有限時間で特異性を発生させることが可能であると広く予想されており、強制項なしでもそれが可能であるはずです。これらの著者たちはまだこれらの目標を完全に達成してはいませんが、近い将来にこれらの目標を達成することが非常に現実的になるようなブレークスルーを遂げました。特に、彼らは3つのより単純なモデル方程式、すなわち非圧縮性多孔質媒体（IPM）方程式、2次元Boussinesq方程式、および3次元非圧縮性Euler方程式で、そのような有限時間爆発を実証しました。（最初のIPM方程式はCórdovaとMartínez-Zoroaによって既に扱われていましたが、AlpögeとBuckmasterは彼らの手法の変種を見つけ、それが他の2つの方程式にも拡張可能であり、Navier-Stokes方程式にも拡張される可能性が高いです。）現代の自動形式化エージェントの時代において驚くほど実現可能であるように、彼らの研究はLeanでも形式化されています。

現在予想されるように、これらの議論はAIによって大きく支援されていますが、著者たちは過去数週間、数学の歴史上最悪の記述と彼らが呼ぶものから、はるかに読みやすくプロフェッショナルな品質の証明へと簡略化・書き直しに取り組んでいます。これはまだ進行中の作業です。残念ながら、上記リンクに記載されている外部イベントのため、彼らは予備的なプレプリントを完全に消化・洗練する前に公開せざるを得ませんでした。それにもかかわらず、Boussinesq論文の序論は少なくともかなり良い状態にあり、初期の出発点として役立つことができます。私もTristan Buckmaster氏から、30分間の電話での会話で論文の主要なアイデアを説明してもらう幸運に恵まれましたが、それをより深く理解するには、おそらく黒板と最新のAIツールの組み合わせで、さらに作業する必要があります。今のところ、そのような仕事の主な価値と見なしているハイライトである、いくつかの主要なアイデアの簡単な要約を書いてみようと思います。これらの問題の実際の解決は、数学的理解と洞察を深めるという主要な目標の単なる代理目標にすぎません。そのような理解なしには、たとえ大衆メディアで宣伝されるよりも数学にとって本質的な重要性がはるかに低いとしても、Navier-Stokesの正則性問題のような問題でさえ、その意味合いは限定的です。

CórdovaとMartínez-Zoroaによる基本的な戦略は、そのような方程式の解を段階的に構築し、以前の（強制された）解に小さな高周波補正を繰り返し追加していく方法です。この方法は、解が爆発時間に向かってより特異になるようにしながら、強制項を良好に保ちます。特定の方程式を扱うのではなく、完全に抽象的な方程式を考えます。

［画像:
(\displaystyle N(u) = f)
］

ここで、［画像: u］は解、［画像: N］は運動方程式を表す非線形微分演算子、［画像: f］は強制項です。もちろん、この設定は完全な解析を行うにはあまりにも一般的ですが、この短い投稿には十分でしょう。

低周波解

［画像:
(\displaystyle N(u_{lo}) = f_{lo})
］

を既に構築できたと仮定し、それを摂動して新しい解

［画像:
(\displaystyle N(u_{lo} + u_{hi}) = f_{lo} + f_{hi})
］

を作成したいとします。この新しい解は、解に高周波補正［画像: u_{hi}］を追加し、強制項に（おそらく高周波の）補正［画像: f_{hi}］を伴います。解の補正［画像: u_{hi}］の振幅を比較的大きく保ちつつ、補正［画像: f_{hi}］の振幅を非常に小さく保つことができ、かつ補正の周波数が各反復で急速に増加する場合、この手順を反復して極限をとることで、解

［画像:
(\displaystyle N(u) = f)
］

を得ることができると期待できます。ここで、［画像: u］は有限時間での爆発を示し、［画像: f］は滑らかなままです。

この戦略を機能させるためには、［画像: f_{hi}］を小さく保つために、［画像: u_{hi}］は差分方程式

［画像:
(\displaystyle N(u_{lo} + u_{hi}) - N(u_{lo}) \approx 0)
］

を近似的に満たす必要があります。もし非線形効果を無視できるなら、これは基本的に線形化された方程式

［画像:
(\displaystyle N'(u_{lo}) u_{hi} \approx 0.)
］

を解くことに相当します。ゲームは、背景解［画像: u_{lo}］を、進化方程式［画像: N'(u_{lo}) u_{hi} = 0］がある種の利用可能な不安定性を示すように設計することです。この方程式の解［画像: u_{hi}］は、初期には指数関数的に小さく開始できますが、爆発時間近くで大きくなります。この時点で、非線形効果が現れて解の解析が非常に困難になることが予想されます。しかし、大きな振幅の出現をちょうど良いタイミングで調整できれば、爆発時間までに振幅が滑らかさを破壊するのに十分大きくなるが、解析を不安定にするほど大きくはない、というスイートスポットに到達できると期待できます。

少なくともBoussinesq方程式の場合、関連論文の序論に記載されている明示的なansatzがあり、爆発に近い時間と原点に近い場所では、［画像: u_{lo}］は空間的に線形に振る舞い、［画像: u_{hi}］は高周波平面波のように振る舞います。驚くべきことに、このansatzは（非線形補正項なしで）正確に解くことができ、必要な不安定性特性を持つ常微分方程式（ODE）変調方程式の明示的なシステムにつながります。これが爆発の基本的なメカニズムです。しかし、完全な議論を厳密にするためには、空間的なカットオフに関連する、膨大な数の技術的な複雑さが存在します。Alpöge–Buckmasterの構成は、古いCórdova–Martínez-Zoroaの構成に比べていくつかの技術的な改善があり、より一般的な流体方程式を扱えるようになっています。私はまだ正確な違いを消化していませんが、ODEはより不安定であり、高周波補正はより良い空間局在性を持つようです。

うまくいけば、将来的に著者たちによるこの素晴らしい結果についての、より良い解説や講演があるでしょう。もし時間があり、結果をより深く理解できれば、フォローアップのブログ記事でさらに詳細を提供できるかもしれません。

編集：また、有限時間爆発へのもう一つの主要なアプローチである、まず数値的または機械学習ツールを使用して近似的に自己相似な爆発プロファイルansatzを特定し、次にそれが実際の解に摂動できるほど安定していることを実証するというアプローチに関する、Ganeshram、Duruisseaux、Anandkumarによる独立したプレプリントもあります。Euler方程式（強制項や境界なし）については、物理情報ニューラルネットワーク（PINN）を使用して数値的に安定な候補解を特定しましたが、解の残差誤差の許容誤差内でその安定性を実際に確立することは、有限時間爆発の完全な厳密な実証に至るまでの主要な挑戦的なタスクです。

第二編集：Euler方程式の特異点形成に向けた最近の進展に関するCórdovaによるサーベイも参照してください。（サーベイはスペイン語ですが、現代の技術を使えば翻訳は容易です。）

第三編集：Javier Gomez-Serranoによる最近の論文に関する講演も参照してください。