A priori Assessment of Tensor-Network Encoding for Isotropic Turbulent Flows

技術要約：等方性乱流におけるテンソルネットワーク符号化の事前評価

問題提起
直接数値シミュレーション（DNS）による乱流の計算は、レイノルズ数の増加に伴うメモリおよびプロセッシング要求量の指数関数的な増大により、計算コストが極めて高い。エネルギー含有エディから散逸スケールに至る全スケールのスペクトルを解像するには、現在のスーパーコンピュータの能力を超える極めて細かい格子と小さなタイムステップが必要となる。次元削減（ROM）手法は存在するものの、それらは通常、モード分解（例：POD）や投影ベースの手法に依存しており、等方性乱流に固有のマルチスケール相関を完全に保存できない可能性がある。同時に、行列積状態（MPS）に代表されるテンソルネットワーク（TN）手法は、量子多体系やその他の分野において高次元データの圧縮に非常に効果的であることが証明されているが、発達した乱流の特性評価への適用可能性については、依然としてほとんど探求されていない。MPSベースの圧縮が、様々な空間スケールにおける乱流統計量にどのような影響を与えるかを定量化する、系統的な「事前評価（a priori assessment）」が不足している。

手法
本研究では、行列積状態（MPS）アンザッツを、低次化された表現を構築するための低次モデル化フレームワークとして採用し、等方性乱流データの簡約化を行う。手法の詳細は以下の通りである：

1. データ符号化： 2つのDNSデータセットを使用する：非圧縮3次元流体速度場（$Re_\lambda \approx 433$）および、同様の流れによって輸送される保存された受動的フィック型スカラー場（$Re_\lambda \approx 180$）。両方の場は$1024^3$の格子上に定義されている。
2. テンソル再整形： 3次元場を、量子ビットベースの符号化（$d=2$）を用い、各物理インデックスが特定の空間方向の長さスケールに対応するように、次数 $N = \log_2(n)$ の非構造化テンソルへと再整形する。
3. インデックスの順序付け： 相関の捕捉を最適化するために、空間テンソルインデックスを再配置する。本研究では、「split」、「staircase」、および「interleaved」の順序付けを比較検討する。interleaved（インターリーブ）順序（$x_1, y_1, z_1, x_2, y_2, z_2, \dots$）は、異なる空間方向からの同等の長さスケールのインデックスを隣接させることで、主要なテンソル間相関を短距離のボンド・インデックス内に局所化させるため、選択された。
4. SVDによる圧縮： テンソルを一連の特異値分解（SVD）を通じて、ボンド・インデックスで接続されたサイト・テンソルの連鎖へと分解する。小さな特異値は、ユーザー定義のカットオフ閾値（$10^{-C}$, ここで $C \in \{2, 3, 4, 5\}$）に基づいて破棄される。この切り捨ては$L_2$（フロベニウス）の意味で最適であり、誤差は破棄された特異値の二乗和によって抑えられる。
5. 再構成と検証： 縮退したMPSを元の$1024^3$格子へと逆変換（コントラクト）する。再構成された場を、インフィデリティ（不忠実度）、圧縮比、および包括的な統計指標のセットを用いて、元のDNSデータと比較する。

主な貢献

• 系統的な事前評価： 本論文は、低次および高次の統計量を幅広く検討することにより、高忠実度DNSデータに対するMPS符号化の厳密な評価を提供し、これまでの限定的な研究を超えたものとなっている。
• インデックス順序付けの最適化： インターリーブ順序が、等方性乱流においてsplit順序よりも大幅に優れていることを実証した。これは、MPSの連鎖構造を、流れの物理的な相関構造（$x, y, z$ 方向の同程度のサイズのスケール同士が統計的に関連していること）に適合させているためである。
• トレードオフの定量化： エネルギー・スペクトル、確率密度関数（PDF）、構造関数、および勾配依存の指標を含む、様々な統計量における圧縮比（ストレージ削減）と精度の間のトレードオフを系統的にマッピングした。

結果

• 圧縮効率： MPS表現は大幅なストレージ削減を実現する。速度場の場合、カットオフ $C=4$ において、約22の圧縮比（元のメモリの約5%を保持）に対し、インフィデリティは$0.18\%$であった。スカラー場の場合、$C=4$ で約6.5の圧縮比（元のメモリの約15%を保持）に対し、インフィデリティは$0.18\%$であった。スカラー場の最高忠実度設定（$C=5$）では、圧縮比は約1.7に低下し、元のメモリの約59%を保持するが、インフィデリティは$0.02\%$となる。スカラー場は勾配が急峻であるため、同等の忠実度を得るために速度場よりも多くのメモリを必要とする傾向がある。
• 統計的忠実度：
  • 低次統計量： 全運動エネルギーおよびスカラーエネルギーは、$C=4$ において相対誤差$0.2\%$未満で回収される。
  • 勾配依存統計量： 平均散逸やスカラー散逸などの空間微分に依存する量は、切り捨てに対して高い感度を示す。$C=4$ において、平均散逸の誤差は速度場では約$4\%$、スカラー場では約$10\%$である。これらの誤差は、$C=5$ ではそれぞれ約$2.6\%$および約$1.8\%$へと減少する。
  • スペクトルおよび構造的特徴： 再構成されたエネルギー・スペクトルは、カットオフ依存の散逸領域まで、コルモゴロフの$k^{-5/3}$慣性スケーリングを正確に追従する。速度勾配テンソルの不変量（$Q-R$）の典型的な「ティアドロップ（涙滴型）」分布、および第三次構造関数のコルモゴロフの4/5則も、$C=4$ および $C=5$ において良好に保持されている。
  • PDF： MPSは、速度成分の近ガウス分布および速度勾配や散逸のヘビーテイル（重い裾）分布を捉える。しかし、極端な事象（PDFの裾の部分）は、低いカットオフ（$C=3$）では過小評価され、正確な再現には$C=4$または$C=5$を必要とする。
• シュミット・スペクトルの解析： シュミット・スペクトルが非対称なピラミッド構造を示すことを明らかにした。小スケールの変動の広範囲な相関を捉えるためには、MPS連鎖の「微細スケール」側において、より大きなボンド次元が必要とされる一方、大スケールの構造はより少ない特異値で捉えられる。

意義と主張
著者らは、これらの知見が、複雑な乱流データセットのスケール可能な低次解析に対するMPSの適合性を支持していると主張している。本研究は、適切な切り捨て閾値が設定されている場合、TNベースの手法が、エネルギー含有構造のような大規模構造だけでなく、混合や間欠性に伴う微細なスケールをも保存できることを示している。MPS構築の適応的な性質（相関が強い箇所でより多くの特異値を保持する性質）により、スケールを横断した流れの物理現象のバランスの取れた保存が可能となる。

結論として、MPSは事前の圧縮および解析には有効であるが、時間発展を伴うシミュレーションや、より複雑な流れ（異方性、壁境界、あるいは反応性流など）への適用には、さらなる探求が必要であるとしている。今後の課題として、3次元乱流の特定の相関構造をより良く捉えられる可能性のある他のTNアーキテクチャ（例：ツリー・テンソル・ネットワーク、PEPS）の調査や、TNの切り捨てと、ラージ・エディ・シミュレーション（LES）におけるサブグリッド・スケール・モデリングとの類似性の探求を提案している。