Finite-amplitude realization of prescribed Reynolds stress inthree-dimensional viscous flow: projection, local forcing, and gridconvergence

技術要約：3次元粘性流における規定されたレイノルズ応力の有限振幅実現

問題の定義
本論文は、非圧縮性流体の構築という文脈における特定の数値的課題、すなわち、有限振幅・有限解像度の粘性ソルバーにおいて、規定されたレイノルズ応力テンソルとその誘起される力がどの程度正確に実現可能かという問題に取り組んでいる。近年の数学的構成（例：OpenAI 2026; Alp¨oge & Buckmaster 2026）は、爆発的シナリオに必要な運動量フラックスを供給するために高周波の振動速度場を利用しているが、これらの研究は存在証明に焦点を当てている。本研究では、「応力合成」という構成ステップを分離し、非圧縮性射影、粘性拡散、非線形自己相互作用、および不均一な背景流との相互作用によって課される実用的な限界を特定する。中心的な問いは、設計された振動場が、ソレノイダル部分空間への数値的射影を経て、目標とする応力と力を局所的に提供するために、著しい歪みを生じることなく生存できるかどうかである。

手法
本研究では、6次コンパクト空間微分と3段階ルンゲ・クッタ時間進展を用いた周期的な3次元非圧縮性ソルバー（Xcompact3d）を採用している。研究は以下の明確なフェーズを経て進行する：

1. 合成と射影による障害： 著者はまず、射影前に目標共分散に一致する速度場を直接コレスキー分解によって構築する方法は、圧力ポアソン射影後に大きく失敗することを示す。この射影は、波長ベクトルに平行な振幅ベクトルの成分を除去し、応力テンソルをオーダー1の係数（例：対角成分が最大2.03倍の係数まで再スケーリングされる）で変化させてしまう。
2. 発散ゼロ基底の構築： 射影による障害を克服するため、6モードのソレノイダル基底を構築する。これには、$\nabla \cdot u' \equiv 0$ が構成上保証されるように、座標軸に沿った波長ベクトルと2つの直交する偏極（正弦および余弦位相）を組み合わせる作業が含まれる。この最小基底を用いた目標テンソル $T$ の実現可能性は、テンソル成分のペアごとの相関に関する3つの不等式によって制約される。
3. 数値検証プロトコル： 本研究では「マッチド・パルス」制御戦略を用いる。正のパルス振幅（$+\varepsilon$）、負の振幅（$-\varepsilon$）、およびゼロ振幅（パルスなし）のシミュレーションを実行する。生成された場の偶成分を $(Q_{+\varepsilon} + Q_{-\varepsilon})/2 - Q_0$ によって抽出することで、奇数項の交差項および背景の進化から応力と力を分離する。
4. 診断的進化： 単独パルス、ガウス渦背景におけるパルス、および平面せん断層におけるパルスのテストを行う。決定的なことに、第2の実験では、非ゼロのダイバージェンス（力）を持つ、空間的に変化する正定値の目標テンソルを規定し、応力と力の両方の局所的な供給をテストする。
5. 格子収束性： 実験は、固定されたスペクトルカットオフ（$K=4, 6, 8$）の下で、$N=128, 192, 256$（および計画されている $N=512$）の解像度で行われる。これは、誘起される力の場の収束を評価するためである。

主要な結果

• 射影と実現可能性： 6モードのソレノイダル構成は、目標とする応力を射影後も保持することに成功した（テストされた一定の目標に対して変化は0.01%未満であった）。これは、直接的なコレスキー法によるオーダー1の誤差と比較して顕著である。しかし、最小基底はすべての正定値テンソルを実現できるわけではなく、ランダムな相関行列に対する実現可能性は44.5%に低下し、ランク1（完全な相関）の目標に対しては0%に崩壊する。
• 粘性減衰と供給の忠実度： 単独パルスの場合、提供される応力の忠実度 $\eta(t)$ は、理論的な形式 $1 - \exp(-2\nu_{eff}k^2t)$ に従い、$\nu_{eff}/\nu = 0.9991$ である。これにより、供給時間に関する厳格な運用上の制約 $t_{deliver} \ll (2\nu_{eff}k^2)^{-1}$ が確立される。
• 非線形自己相互作用： 顕著な自己相互作用は $\varepsilon k t \approx 1$ 付近で始まる。テストされた範囲内では、粘性減衰からの偏差は $\varepsilon k t \approx 1$ で有意（オーダー $10^{-2}$）になり、$\varepsilon k t \approx 4.32$ で大幅なエネルギー損失（21.9%）をもたらす。
• 局所的な力の供給： 空間的に変化する目標 $T(x)$ と非ゼロの力 $F(x) = -\nabla \cdot T(x)$ を用いた実験において、$N=128$ かつ $t=0.01$ における局所的な応力および全力の誤差は、それぞれ0.316%および1.796%であり、事前宣言された閾値を通過した。
• 格子収束性： $t=1.9$ において、誘起される3成分の力場の解像度間の差は、明確な収束を示した。$K=6$ における $N=192$ と $N=256$ の間の相対的な $L_2$ 差は、全力の成分に対して0.0114%であり、$N=128$ と $N=192$ の間の差（0.166%）よりも有意に小さかった。
• 数値的・物理的散逸： 本研究は、数値的散逸と物理的な非線形減衰を区別している。観察された「過剰な減衰」は、数値誤差ではなく二次非線形性（自己相互作用）であると特定された。なぜなら、それは振幅の2乗に比例し、かつ格子解像度に依存しなかったためである。6次スキームの実際の数値的散達は、粗い解像度（$k\Delta x \approx 0.98$）においても無視できる程度（$<0.02\%$）であることが判明した。

意義と主張
本論文は、粘性非圧縮性流において、規定された非定常目標に対して局所的な応力と力の供給が行われることを示す、検証された数値的事例を提供すると主張している。特に、進化するガウス渦背景における非定常な目標について扱っている。

• 限定的な範囲： 著者は、これらの結果がすべての目標テンソルに対する一般的な実現可能性定理を構成するものではなく、また、フィルタリングされていない全スペクトルの収束や、点別最大ノルムの予算閉鎖を確立するものでもないことを明示している。知見は、特定のテンソル、単一の背景パラメータ点、コア体積RMSノルム、および低モードスペクトルカットオフに限定されている。
• 運用上の制約： 本研究は、スケール分離、粘性供給限界、弱い非線形性（$\varepsilon k t \lesssim 0.5$）、および解像度（$k\Delta x \lesssim 1.2$）によって境界付けられる、局所的な供給が可能となるパラメータの「実行可能なウェッジ（くさび形領域）」を定義している。テストされたパラメータにおいて、この窓は $Re=10^5$ で存在するが、より低いレイノルズ数では空となる。
• 限界： 論文は、導入部で引用された特定の構成（OpenAI 2026; Alp¨oge & Buckmaster 2026）が、それらの明示的なプロファイルなしには、この実行可能な窓内で動作していることを確認できないことを認めている。さらに、計画されている $N=512$ のパラメータスキャンおよび第二のソルバーによる独立した再現実験はまだ実施されていないため、現在の主張は完了したシミュレーションに限定される。

要約すると、本研究は、規定されたレイブル応力を粘性ソルバー内で合成するために必要な代数的および数値的なメカニズムを分離し、検証しており、適切なパラメータ選択と発散ゼロの基底を用いることで、弱非線形性と十分な解像度の範囲内に留まる限り、局所的な応力と力の供給が高精度に達成できることを示している。