Estimation of Wavenumber Spectra from Particle Tracking Velocimetry Data

技術要約：粒子追跡流速計（PTV）データからの波数スペクトルの推定

問題提起
粒子追跡流速計（PTV）は、3成分・3次元の速度場データを提供するが、シード粒子のランダムな分布により、空間サンプリングは本質的に不規則である。この不規則性により、中間的な補間を行うことなく、従来の高速フーリエ変換（FFT）を直接適用して波数スペクトルを計算することが妨げられる。補間はサンプリングの問題を解決するものの、得られるスペクトル推定値のバイアスと分散に関する不確理性を導入する。具体的には、不規則にサンプリングされたデータに対して、期待値（バイアス）と変動性（分散）のバランスを最も適切に保つ推定手法がどれであるかが不明確である。特に、高波数領域の散逸領域を解像しようとする場合において、中程度のシード密度で散逸領域のスペクトル推定を可能にするような、無偏な推定法が存在するかどうかという問いがある。本研究では、この問題に取り組む。

手法
著者らは、レイノルズ数2000の統計的に定常な2次元コルモゴロフ流（k-flow）の直接数値シミュレーション（DNS）から生成された合成データセットを利用した。この流れは、二重カスケード構造（逆エネルギー伝達および直接エンストロフィー伝達）を特徴としている。PTV測定をシミュレートするため、DNS速度場内に仮想トレーサー粒子を様々な数密度（$N = 5,000, 10,000, 20,000$）でランダムに配置し、速度はオイラー格子から補間され、5%のガウスノイズが付加された。

4つの異なるスペクトル推定法を、完全に解像されたDNSデータから導出された参照スペクトルと比較した：

1. 直接スペクトル推定（NUDFT）： 再構成を行わず、散布されたサンプルに直接適用される非一様離散フーリエ変換。
2. 補間 + FFT： 粒子データを均一なDNS格子上に自然隣接補間し、その後に標準的なFFTを適用する手法。
3. 物理情報に基づくニューラルネットワーク（PINN）+ FFT： 渦度・速度PINNを用いた速度場の再構成。研究では2つのデータ忠実度損失関数、すなわち平均二乗誤差（$L_2$）と平均絶対誤差（$L_1$）を比較した。疎なデータ領域における$L_2$損失の過度な平滑化傾向を軽減するために、$L_1$ノルムが選択された。
4. ファジースロット相関（FSC）： 自己相関関数に基づく推定法。これは「ファジースロッティング」（三角形の重み関数を用いて近傍のビンにペアの寄与を分配する）と、「局所正規化」（相関のない測定ノイズの影響を抑制する）を組み合わせたものである。

これらの推定法の性能は、粒子数密度（$N_p$）およびアンサンブルフレーム数（$N_f$）の関数として、期待値（バイアス）および分散の観点から評価された。

主な結果

• NUDFT： 原理的には無偏であるが、この手法は参照スペクトルのエネルギー含有領域よりも数桁低い平坦なノイズフロアを生じさせ、密度補正なしではスペクトル形状を回復させる上で効果的ではなかった。
• 補間： この手法はローパスフィルタとして機能した。大規模なエネルギーを正確に回復したが、局所的な粒子間隔よりも小さい構造を解像できないために、散逸スケール（$k_\eta$）付近の高波数において深刻な正のバイアス（オーバーシュート）を導入した。
• PINN：
  • $L_2$-PINNは$k^{-3}$の勾配を再現したが、過度な平滑化により高波数領域全体で持続的な負のバイアス（過小評価）を示した。
  • $L_1$-PINNは、このバイアスを大幅に軽減し、$k_\eta$付近のDNS参照値を正確に追跡し、中間波数においてわずかなバイアスを示すのみであった。極めて重要な点として、$L_1$-PINNは単一フレームから滑らかに収束したスペクトルを生成し、アンサンブルサイズに対して鈍感であった。
• ファジースロット相関（FSC）： この手法は、粒子数密度によって決定される波数限界まで、無偏な推定を提供できた唯一の手法であった。しかし、単一フレームのFSC推定は強い振動（高い分散）を示した。これらの変動はアンサンブル平均によって抑制され、不確実性は$N_f^{-1/2}$として減衰した。
• 解像限界： FSC推定量の解像帯域幅（$k_{max}$）は、粒子密度の増加とともに増加したが、$\sqrt{N_p}$に対して劣線形にスケールした。テストされた最高密度（$N_p=20,000$）において、FSCは約$0.92 k_\eta$に達した。

計算コスト
本研究は、精度と計算努力の間のトレードオフを強調した。格子ベースの手法（補間およびPINN）は計算効率が高かった（フレームあたり数秒）。対照的に、FSCおよびNUDFTは散布サンプルに直接作用するため、FSCのコストはペアの相関計算により$O(N_p^2)$でスケールする。$N_p=20,000$において、FSCは1フレームあたり約72秒を要し、PINNアプローチよりも大幅に長い時間を要した。

意義および主張
論文は、推定法の選択は、バイアス、分散、および計算リソースに関する具体的な要件に依存すると結論付けている：

• $L_1$-PINNは、小さな残留バイアスが許容される場合に、高速で滑らかな単一フレーム推定を得るための最適な手法として提示されている。
• 局所正規化を伴うFSCは、分散を抑制するために十分な計算リソースとアンサンブルサイズが利用可能であれば、無偏なアンサンブル平均スペクトルを散逸領域まで得るための優れた手法として特定されている。

著者らは、本研究は2次元の流れに対して行われたものであるが、これらの推定量のバイアス・分散特性に関する根本的な結論は3次元の流れにも保持されると主張している。ただし、FSCの計算コストは、粒子のペアの立方体的な増加に伴い大幅に増大することに注意が必要である。本研究は、不規則なPTVデータから無偏なスペクトル推定を行うことは可能であるが、バイアスと分散のトレードオフを管理するために、アンサンブル平均と適切な推定法の選択が必要であることを確立した。