arXiv雑要約
数値解析 - 2026/03/25 公開
連続データ同化による高レイノルズ数ナビエ-ストークス方程式に対するUzawa法ソルバーの高速かつ堅牢な収束 [math.NA, cs.NA]目的：高レイノルズ数ナビエ-ストークス方程式に対するUzawa法ソルバーの収束性向上
- 流体シミュレーションは，気象，海洋，航空力学など広範な分野で不可欠であり，高精度な数値解法が求められる。
- ナビエ-ストークス方程式の数値解法は，高レイノルズ数において計算コストが高く，収束性が保証されない場合がある。
- 連続データ同化を用いることで，従来のUzawa法ソルバーの収束を加速し，高レイノルズ数における安定性を確保することを目指す。
- 連続データ同化(CDA)をUzawa法に組み込んだCDA-Uzawa法が，従来のUzawa法よりも収束が速くなることを理論的に証明した。
- 十分な部分解データを用いることで，CDA-Uzawa法は複数の解が存在する場合でも，任意の高レイノルズ数において収束することが示された。
- ノイズを含むデータに対しても，ノイズの大きさ以上の範囲で収束性が保たれ，CDA-Uzawa法の限界に達した時点でニュートン法に移行する戦略が提案された。
単純だが単純ではない：2つの楕円体間の最小距離を見つけるための表面スライディング法 [cs.CG, cs.NA, math.NA]目的：2つの楕円体間の最小分離距離の算出
- コンピュータグラフィックスやシミュレーション等の分野において，形状間の距離計算は重要な要素である。
- 既存手法では，計算コストが高い，または局所解に陥りやすいという問題点がある。
- 本研究は，効率的かつ高精度に最小距離を算出する新たな手法を提案することで，上記問題を解決する。
- 提案手法は，楕円体表面上の2点をパラメータ空間で更新し，連結線分の張力をガイドとして滑らせる。
- 連結線分が楕円体表面に垂直になった際に最小距離に到達し，表面点が移動しなくなる。
- 数値実験により，高い精度，安定性，ロバスト性が確認された。また，幾何学的な解釈が容易で拡張性も高い。
ハミルトニアンモンテカルロ法のためのアルゴリズム的ウォームスタート [cs.DS, cs.LG, cs.NA, math.NA, math.ST, stat.ML, stat.TH]目的：高次元確率分布からのサンプリング効率改善
- 統計学，工学，自然科学における重要な課題であり，効率的なサンプリング手法が求められている。
- ハミルトニアンモンテカルロ法は広く利用されているが，次元数が高くなると収束に要する反復回数が問題となる。
- ウォームスタートの計算ボトルネックを解消し，高精度サンプリングの効率を向上させる。
- 強ログ凹性および３次微分条件を満たす分布において，ノンメトロポライズドHMCによるウォームスタートが$\tilde{O}(d^{1/4})$で実現可能となった。
- ウォームスタートを利用したメトロポライズドHMCにより，$\tilde{O}(d^{1/4})$という最速のアルゴリズムが実現された。
- これまでの$\tilde{O}(d^{1/2})$を改善し，高精度サンプリングにおける次元複雑性の問題を解決した。
非線形放物型方程式に対する行列フリー安定BDFスキーム：最大値原理とエネルギー安定性の無条件保存 [cs.HC, cs.RO, cs.RO, math.NA, cs.NA]目的：非線形放物型方程式に対する安定BDFスキームの開発
- 偏微分方程式の数値解法は，物理現象のシミュレーションや工学設計において不可欠である。
- 高次の時間離散化法では，安定性や最大値原理の保存が困難な場合が多い。
- 行列フリー実装により，大規模問題への適用可能性を向上させる。
- 提案手法は，最大値原理とエネルギー安定性を無条件に保存する。
- 固定点反復と pointwise cut-off 戦略により，行列フリーの実装を実現している。
- 混合境界条件を持つ問題において，ETD法と比較して高いロバスト性と効率性を示す。
非圧縮性ナビエ・ストークス方程式に対する物理情報ニューラルネットワークの一般化誤差について [cs.LG, cs.NA, math.AP, math.NA]目的：非圧縮性ナビエ・ストークス方程式の解近似における物理情報ニューラルネットワークの一般化誤差の上界
- 流体シミュレーションは，気象予測，航空力学，バイオエンジニアリングなど，幅広い分野で不可欠である。
- 従来の数値解法は計算コストが高く，高次元問題への適用が困難であるという課題がある。
- 物理情報ニューラルネットワークは，少ない計算資源で高精度な解を得る可能性を秘めている。
- 本研究では，深さ2のニューラルネットワークを用いたPINNの一般化誤差に対する初めての上界を導出した。
- 導出された一般化誤差は，ネットワークの幅に明示的に依存せず，流体の運動粘性や損失正則化パラメータで特徴付けられる。
- 次元に依存しないサンプル複雑性の上界が得られ，新しい活性化関数の利用が示唆された。
データ依存型RBF近似 [math.NA, cs.NA]目的：不連続点近傍における振動低減を目的としたRBF補間手法の改良
- RBF補間は高次元データの近似に有効だが，不連続点近傍で振動が生じやすい。
- 従来のRBF補間では，不連続点近傍の振動を抑制する効果的な手法が不足していた。
- データ依存的に形状パラメータを変化させ，不連続点近傍での振動を最小化することを目指す。
- 提案手法は，形状パラメータをデータ依存的に変化させることで，局所的にKronecker delta関数に類似したカーネル関数を実現する。
- 滑らかさ指標を用いることで不連続点を検出し，形状パラメータの適応的な調整を可能にした。
- 数値実験により，提案手法が様々なカーネルタイプにおいて，不連続点近傍の振動を大幅に低減することが示された。
安定かつ費用対効果の高い行列多項式評価のための体系的フレームワーク [math.NA, cs.MS, cs.NA]目的：行列多項式評価における安定な係数集合の獲得と検証
- 行列計算は科学技術計算の根幹であり，効率的な手法が求められている
- 既存手法では計算コストが高く，大規模行列での実用性に課題があった
- 数値的に安定した係数集合を選定し，計算コスト削減と精度の両立を目指す
- 本研究では，Paterson-Stockmeyer法よりも行列積の回数を1つ削減しつつ，数値安定性を確保するフレームワークを提示した
- 提案手法により，行列指数関数や等比級数といった応用例において，既存法と同等の精度を維持した
- フレームワークは，数値精度を評価する信頼性指標を提供し，安定な係数集合を効率的に特定するツールとしてMATLABで実装された
非グローバルなリプシッツ定数を持つ確率偏微分方程式に対する有限要素全離散スキームの長時間の誤差解析 [cs.RO, math.NA, cs.NA]目的：非グローバルなリプシッツ定数を持つ確率偏微分方程式の長時間の近似
- 確率偏微分方程式は，流体力学や金融工学など，様々な分野で現れる重要な数理モデルである。
- 既存の数値解法では，時間ステップ幅や空間離散化の精度が，誤差の大きさに大きく影響する。
- 時間ステップ幅と空間離散化の精度に制限を受けない，安定した長時間の近似スキームを開発すること。
- 提案スキームは，バナッハ空間$L^{r}(D)$における一様時間境界を保存し，時間空間離散化ステップサイズの比率に制限がない。
- 滑らかでない初期値を持つ決定論的線形放物型方程式に対する誤差評価を新たに導き出し，数値偏微分方程式の分野に貢献する。
- 空間・時間ホワイトノイズおよびトレースクラスノイズの場合において，一様時間収束率を得ることで，低規則性の設定下での誤差を詳細に分析した。
多次元非線形双曲型問題に対する任意の次数定常性保存安定化有限要素法：重力を含むオイラー方程式への応用 [math.NA, cs.NA]目的：多次元非線形双曲型保存則の数値解法
- 流体やプラズマなどの非線形現象をシミュレーションする上で重要な役割を担う
- 高精度な計算を行う際に，数値的な不安定性や振動が発生しやすい
- 定常解を正確に近似し，数値的な安定性を向上させることを目指す
- 本研究では，任意の次数で定常性を保存する安定化有限要素法を開発した。
- この手法は，従来のSUPG法や有限体積法と比較して，ロバスト性と精度が向上することが数値実験により示された。
- 特に，圧縮性オイラー方程式への適用により，等温静水圧平衡状態の機械精度での保存が可能となった。
指数的に圧縮された非ガウス集約分布に対する高分解能テンソルネットワークフーリエ法 [math.CO, cs.DB, cs.DM, stat.ML, cs.LG, cs.NA, math.NA, quant-ph]目的：独立な確率変数の加重和の特性関数表現の効率的圧縮
- 金融工学等の分野において，複雑な確率分布の正確な計算が重要である。
- 従来の数値計算手法では，高次元の確率分布を扱う際にメモリや計算時間の制約が生じる。
- テンソルネットワークを用いることで，高次元確率分布を効率的に表現し，計算コストを削減することを目指す。
- 特性関数は量子化されたテンソルトレイン（QTT）表現において低ランク構造を示し，非ガウス分布の指数的な圧縮を可能にする。
- ベルヌーイ変数と対数正規分布の加重和において，QTT構造による効率的な計算が確認された。
- 本手法により，標準的なハードウェアで高分解能の確率分布計算（$N = 2^{30}$）が可能となり，VaRやESの計算に貢献する。
制約付きサンプリングのための運動ランジュバン分割スキーム [eess.SP, cs.ET, stat.ME, cs.NA, math.NA, stat.ML]目的：制約付きサンプリング手法の開発
- 統計計算において，特定の制約下でのサンプリングは不可欠な課題である。
- 既存手法は，計算効率や収束性において課題が残されている。
- 運動ランジュバン動力学に基づく分割スキームを用いて，効率的なサンプリングを実現する。
- 提案手法は，既存の非分割法と比較して，より有利な計算複雑性を示す。
- Wasserstein縮小性に関する理論的な知見と収束結果が得られた。
- 数値実験により，提案手法の有効性が，簡単なモデルやベイズ線形回帰で確認された。
円筒面上の拡散に対するスペクトルラプラシアン形式を用いた厳密解析PGSE信号 [physics.med-ph, cs.NA, math-ph, math.MP, math.NA]目的：円筒面上の拡散に対する厳密な解析的解
- 拡散MRIは，生体組織内の分子拡散を非侵襲的に評価する強力な手法である。
- 円筒面上の拡散の既存モデルは，高拡散重みにおいて精度が低下する問題がある。
- 本研究は，任意の勾配時間と間隔における厳密な解析的解を導出し，その精度を向上させる。
- Bloch-Torrey方程式の厳密解を導出し，円筒面上の拡散に対する拡散MRI信号の式を得た。
- スペクトル行列形式により，拡散プロパゲーターやスピン位相分布に対する近似を回避した。
- 数値計算戦略を開発し，高速な評価とモンテカルロシミュレーションによる検証を行った。
より高速なMDS行列検証のための巡回行列部分行列群の構築 [math.NA, cs.NA]目的：MDS（最大距離分離可能）行列であることの検証に必要となる巡回行列の部分行列群
- 符号理論や暗号理論において，MDS行列は重要な役割を担うため，その検証は必須である。
- 一般的なMDS行列検証は計算量が大きく，特に大規模行列においては現実的な時間で検証が困難となる。
- 巡回行列の構造を利用し，検証に必要な部分行列の数を削減することで，高速な検証を可能とすること。
- 提案手法により，検証に必要な部分行列群のサイズを約2次元分削減することができた。
- 従来の検証方法と比較して，計算量の削減が期待できる。
- 本研究は，大規模な巡回MDS行列の検証を効率化する上で貢献する。
Truncated Tensor SVD近似のための確率的ブロック・クリロフ法 [math.NA, cs.NA]目的：Truncated Tensor SVD近似の効率化
- 大規模データ解析において，テンソル分解は重要な役割を果たす。
- Truncated Tensor SVDの計算コストは，特に大規模データに対して課題となる。
- 確率的ブロック・クリロフ法を用いて，計算コストを削減し，効率的な近似を目指す。
- 提案手法は，合成データおよび実データを用いた実験で有効性が確認された。
- 本手法は，Truncated Tensor SVDの計算時間短縮に貢献しうる。
- データ補完やデータ圧縮への応用可能性も示された。
圧縮性ナビエ・ストークス方程式に対する陰解的二重時間ステップ陽性保存エントロピー安定解法 [math.NA, cs.NA]目的：圧縮性ナビエ・ストークス方程式の数値解法
- 高精度な数値シミュレーションは，航空力学やプラズマ物理学等の分野で不可欠である。
- 陽解法では時間ステップ幅が制限され，計算コストが増大する可能性がある。
- 陰解法により，時間ステップ幅の制限を緩和し，大規模シミュレーションを可能とする。
- 本研究で提案する二重時間ステップ法は，エントロピー安定性と陽性保存性を数学的に保証している。
- これにより，衝撃波や接触不連続を含む超音速粘性流のシミュレーションにおいて，安定かつ正確な結果が得られる。
- 特に，高レイノルズ数における定常・非定常流れの計算に適しており，計算効率の向上が期待される。
非線形反応拡散方程式の安定化のためのペナルティに基づくフィードバック制御と有限要素解析 [math.NA, cs.NA, math.OC]目的：非線形反応拡散方程式の安定化手法
- 反応拡散方程式は，物理，化学，生物学など広範な分野で現象を記述する上で重要である。
- 境界フィードバック制御における安定性の解析は困難であり，特に非線形問題では解析が複雑になる。
- ペナルティ法を用いることで，境界制御問題を近似的に解き，安定性を解析することを試みる。
- ペナルティパラメータに関する二階ノルムにおいて，相当するロビン問題の安定性が確立された。
- ペナルティパラメータをゼロに近づけることで，ペナルティ付き制御問題の解が，元のディリクレ境界フィードバック制御問題の解に収束することが示された。
- 数値実験により，理論的な結果が検証され，ペナルティ付き解が元の解に収束することが確認された。
クォータニオン行列の低ランク近似に対するパス効率の良いランダムアルゴリズム [eess.SY, cs.SY, math.NA, cs.NA]目的：クォータニオン行列の低ランク近似
- データ量の増大に伴い，大規模行列の効率的な処理が重要になっている。
- 既存手法は，入力行列に対するパス回数を削減するパス効率が考慮されていない。
- パス回数を減らしつつ，近似精度を調整可能なアルゴリズムを開発すること。
- 提案手法は，ユーザーが指定したパス回数内で低ランク近似を行うランダムアルゴリズム群である。
- スペクトルがゆっくりと減衰する行列に対し，ブロック・クリロフ部分空間法を拡張し，収束を加速する。
- パス回数が増加するにつれて，期待される近似誤差が指数関数的に減少することが理論的に示された。
疎なBスプライン近似空間のテンソル積に対する誤差評価 [math.NA, cs.NA]目的：疎なテンソル積Bスプライン近似空間の誤差評価
- 複雑な形状に対する高精度な数値解析が求められており，効率的な近似空間の構築が重要である。
- 従来のテンソル積空間では，高次元問題に対して自由度が指数的に増加する問題がある。
- 疎なテンソル積空間を用いることで，自由度を削減し，計算コストを抑制することを目指す。
- 疎なグリッド結合法と階層的サブ空間分解法は数学的に同等であることが示された。
- 適切な正則性条件の下で，疎なテンソル積空間は標準的なテンソル積空間と同程度の近似精度を，より少ない自由度で達成することが示された。
- 非テンソル積領域においては，最適な収束率を得るために，混合微分に関するより強い正則性条件が必要となることが示された。
非構造化三角網目上のDG定式化におけるロバストなPAMPAスキーム：境界保存，振動の除去，境界条件 [math.NA, cs.NA]目的：非構造化三角網目上でのDG定式化におけるPAMPAスキームの改良
- 流体計算などのシミュレーションにおいて，精度と安定性の高い数値解法が求められている。
- 従来の数値解法では，振動が発生したり，境界条件の取り扱いが難しい場合がある。
- 境界保存性と非振動性を両立した，よりロバストな数値解法を開発すること。
- 本研究で提案するスキームは，局所的に保存的であり，質量行列の逆演算を必要としない。
- PAMPAとDiscontinuous Galerkin法の関連性を利用し，厳密な境界条件の実装方法を定式化した。
- 数値実験により，提案スキームが境界保存性と非振動性を広範なベンチマークで実現することを確認した。
楕円型偏微分方程式における未知境界条件緩和：内部点計測による手法 [math.NA, cs.NA]目的：ポアソン方程式の未知境界データ回復に関する最適化問題
- 偏微分方程式は自然現象や工学問題を記述する基礎であり，その数値解法は重要。
- 境界条件が不明な場合，解の精度向上が課題となる。特に内部点での計測は困難。
- 内部点計測データを用いた，より精度の高い境界条件回復アルゴリズムを開発。
- 本研究では，計測点におけるRiesz表現子の点ごとの誤差評価を導出した。
- この誤差評価を応用することで，回復アルゴリズムの性能に関する改善された推定値が得られた。
- 解の滑らかさに関する制約を緩和し，より広範な問題への適用が可能となった。
音響逆散乱問題に対するコントラスト源反転型手法の収束解析 [math.NA, cs.NA]目的：音響逆散乱問題に対するコントラスト源反転型手法の収束性
- 逆散乱問題は，医療診断や非破壊検査など幅広い分野で重要である。
- コントラスト源反転型手法は効率的だが，厳密な収束性の証明が課題となっていた。
- 固定周波数における非線形逆散乱問題を解く反復手法の収束性を確立すること。
- 本研究では，新しい$\ell_1$近接項を用いた2つの反復正則化コントラスト源反転型(IRCSI型)手法を提案した。
- 提案手法(IRCSI法，IRSOM法)は，元のCSI法やSOM法と同程度の計算量を持ち，自然かつ弱い条件の下で大域的収束が証明された。
- 数値実験により，提案手法の収束性と性能が示された。
クープマン作用素と遷移行列近似のための動的に学習される基底 [math.DS, cs.LG, cs.NA, math.NA]目的：複雑な非線形動的システムの表現と解析
- 複雑系理解への貢献が期待され，制御や予測といった応用範囲が広い。
- 効率的なスペクトル推定が困難であり，計算コストが高いという課題がある。
- 動的システムに特化した基底関数を学習し，近似精度向上を目指す。
- 提案手法により，クープマン作用素と遷移行列のスペクトル特性を効果的に抽出可能となった。
- 学習された基底は，システムのダイナミクスに適応的に変化し，高精度な近似を実現する。
- 固有関数や不変測度の効率的な復元にも貢献し，動的システムの理解を深める。
$\mathbb{R}^n$における外微分と平均値の等式 [math.DG, cs.NA, math.NA]目的：外微分の一般化された捉え方によるベクトル解析と解析における古典的な結果の再検討
- ベクトル解析は物理学や工学における様々な現象の基礎となる重要な分野である。
- 微分形式の滑らかさの条件が厳しく，応用範囲が限定される場合がある。
- 微分形式に対する平均値の定理を導き，Stokesの定理のより一般的な定式化を目指す。
- 外微分を「無限小のフラックス」の尺度と解釈する視点から，高次元の平均値の定理を導出した。
- 微分形式の完全な$C^1$滑らかさを必要としない，Stokesの定理の自然な定式化を提案した。
- 微分形式へのブラックボックスアクセスのみを必要とする外微分計算アルゴリズムを提案した。
有限体上の多項式環の部分集合に存在する最小分母 [math.NT, cs.NA, math.NA, math.PR]目的：有限体上の多項式環の部分集合における最小分母の分布
- 数論におけるディオファントス近似の研究は，整数や多項式の有理数近似に関する重要な課題である。
- 最小分母の分布に関する研究は，近似の精度と効率を評価する上で不可欠だが，解析が困難である。
- 本研究は，有限体上での最小分母の分布を解析し，離散的近似と連続的近似の同値性を示すことを目指す。
- 無限部分集合において，任意の次元と次数に対して，2つの確率分布が一致することが示された。
- これは，実数におけるBalazard-Martinの結果よりも強い主張であり，離散近似と連続近似の間のより深い関係を示唆する。
- 本研究は，ディオファントス近似の有限体上における性質を理解するための新たな視点を提供する。
高次元フォッカー・プランク方程式に対する適応確率フロー残差最小化 [physics.comp-ph, cs.LG, cs.NA, math.NA]目的：高次元フォッカー・プランク方程式の効率的な解法
- 物理学や確率力学において高次元方程式の解法は不可欠である。
- 次元の呪いにより，従来の数値解法は高次元で計算コストが増大する。
- 高次元における計算コストを削減し，精度を維持する解法の開発。
- 適応確率フロー残差最小化(A-PFRM)法を提案し，ヘッセ行列の明示的な計算を回避した。
- 連続正規化フローとハッチンソン・トレース推定器を利用し，計算複雑度をO(d)に削減した。
- 様々なベンチマークテストで，次元の呪いを軽減し，高い精度と一定の計算時間を維持することを示した。
- 1