arXiv雑要約
数値解析 - 2026/06/16 公開
実数次数の複素引数に対するベッセル関数の効率的な多倍長計算 -- 第2部: 第二種修正ベッセル関数 $K_\nu(z)$ [math.NA, cs.NA, math-ph, math.MP]目的：複素引数と実数次数に対する第二種修正ベッセル関数 $K_\nu(z)$ の計算アルゴリズム
- ベッセル関数は，物理学，工学における様々な問題の解に現れるため，高精度な計算が重要である。
- 既存の方法では，パラメータ範囲や精度において制限があり，計算の信頼性が課題となっていた。
- 広範なパラメータ範囲において，高精度かつ効率的なベッセル関数計算を可能とする。
- 本研究では，冪級数，漸近展開，安定な漸化式などを適応的に選択するアルゴリズムを開発した。
- Fortranによる実装は，倍精度および四倍精度をサポートし，信頼性の高い計算領域を拡大している。
- ベンチマークテストの結果，既存の方法と比較して大幅に高速であり，数値的な不安定性を回避できることが示された。
実数次数と複素引数に対するベッセル関数の効率的な多倍長計算 - 第3部：第1種および第2種のベッセル関数 $J_{\nu}(z)$ および $Y_{\nu}(z)$ [math.NA, cs.NA, math-ph, math.MP]目的：実数次数と複素引数を持つベッセル関数の数値計算アルゴリズムと多倍長Fortran実装
- ベッセル関数は物理学，工学等の幅広い分野で頻繁に用いられる基礎的な特殊関数である。
- ベッセル関数の高精度計算は，計算コストが高く，数値安定性の問題が伴う。
- 計算精度と効率を両立したベッセル関数のアルゴリズムを開発し，適用範囲を拡大する。
- 提案アルゴリズムは，既存のAlgorithm 644と比較して，二重精度において実行時間が35-67%向上し，信頼性の高い結果を得る。
- Algorithm 912と比較して，二重精度では同等の精度，四重精度では大幅に高い精度を実現し，計算コストは大幅に低い。
- 本アルゴリズムは，安定性と精度を維持しながら，より広い領域で信頼性の高い計算が可能である。
定数係数時間調和波伝播問題に対する伝送変数を用いたHDG法 [math.NA, cs.NA]目的：時間調和波問題の高速反復解法
- 波伝播問題は，工学や科学の多くの分野で重要であり，正確な数値解析が求められる。
- 時間調和波問題の反復解法は，固有の性質により収束が遅いことが課題であった。
- 本研究は，標準HDG法の改良版であるCHDG法を一般化し，高速な反復解を可能にすることを目指す。
- 提案手法は，定数係数を持つ波伝播問題に対して適用可能であり，多様な波の種類を統一的に扱える。
- 得られたハイブリダイズ系は，固定点反復に適しており，良好な解法性を示す。
- 数値実験の結果，音波および渦波に対して，いくつかの反復スキームの収束性が確認された。
対流熱伝達を伴う最適断熱問題の有限要素近似 [math.NA, cs.NA]目的：対流熱伝達を伴う最適断熱問題
- 熱設計において断熱は重要であり，エネルギー効率の向上に不可欠である。
- 断熱材の配置最適化は非線形で複雑であり，効率的な解法が求められている。
- 有限要素法による効率的な数値解法を開発し，断熱設計の最適化を目指す。
- 本研究では，有限要素離散化による最適断熱問題のモデルを構築した。
- 連続設定の構造を保持する質量凝縮求積法を用いて，表面積分を近似した。
- 離散解の適切な性質と連続解への収束性，そしてブロック座標降下法の線形収束性が示された。
エネルギー保存型非オフセット電磁ポテンシャルPIC法 - 第I部：非相対論的一般化運動量定式化 [math.NA, cs.NA]目的：非相対論的Vlasov-Maxwell系の数値シミュレーション手法の開発
- プラズマ物理学において，高精度かつ効率的な数値シミュレーションは不可欠である。
- 従来のPIC法では，エネルギー保存性が十分でなく，長時間のシミュレーションで誤差が蓄積する。
- エネルギー保存性を厳密に保証する新しいPIC法を構築し，シミュレーションの精度と安定性を向上させる。
- 本研究では，電磁ポテンシャルに基づく非オフセットPIC法を提案し，Crank-Nicolson法を用いて場の更新を行う。
- 軌道平均化された散乱，集約，粒子移動スキームを導入することで，エネルギー保存性を実現した。
- 3次元のコールド二重流不安定性シミュレーションにおいて，本手法が厳密なエネルギー保存性を示すことを確認した。
物理法則適合な潜在的二重表現 [cs.LG, cs.NA, math.NA]目的：物理法則を尊重する高速予測モデルの学習
- 複雑な物理システムの予測・シミュレーションには，計算コストが低い代替モデルが不可欠である。
- 時間依存問題において，既存の代替モデルは物理法則や構造を必ずしも保存しないという課題がある。
- 物理法則を設計段階で組み込み，より信頼性の高い代替モデルを構築することを目的とする。
- 潜在的二重表現の枠組みを拡張し，潜在空間でのダイナミクスに物理構造を保存または消散させる制約を導入した。
- 制約伝達の視点から，元の状態空間の物理構造と潜在空間の制約との関連性を示し，潜在空間での制約強化が復号後の物理的欠陥の制御を改善することを示した。
- 線形・二次不変量保存や消散不等式の強制に関する代数的条件を導出し，数値実験で制約充足度，構造的忠実度，長期挙動の改善を確認した。
逆熱点源問題における Hankel 行列の零点に基づく源の数え上げ原理 [math.NA, cs.NA]目的：逆熱点源問題における点源の数の特定
- 熱伝導現象の解析は，工学，医学，地球科学など広範な分野で不可欠である。
- 熱源の数を事前に知る必要があるため，未知の熱源数の特定が困難であった。
- 境界での熱流束の測定値から熱源の数を直接推定する手法を確立すること。
- Hankel 行列の零点の消失次数が，真の熱源数を超えたときに変化することを示した。
- 熱源数は，零点特性の最初の非零の等高線数によって特徴付けられる。
- 測定ノイズや境界離散化などの摂動に対して，零点数の安定性が確認された。
高次元ウィグナー・ポアソン系の構造を保存する適応的ランクアプローチ [math.NA, cs.NA]目的：高次元ウィグナー・ポアソン系の数値解法
- 量子力学的電子ダイナミクスのシミュレーションにおいて，フェーズ空間モデルは重要である。
- 高次元フェーズ空間と非局所ウィグナーポテンシャルがシミュレーションの限界となる。
- 階層的タッカー形式を用いた適応的ランクソルバーにより，計算効率と精度を向上させる。
- 提案手法は，フーリエ・エルミート対称性を維持しつつ，離散的な保存則を保証する。
- 2次元2速度および3次元3速度のテストにおいて，丸め誤差レベルでの保存性と実数値逆変換の維持が確認された。
- フェーズ空間テンソル全体を組み立てることなく，長時間の3次元3速度ウィグナー・ポアソンシミュレーションが可能となった。
空間におけるハイブリッド高次法と時間における陰解法による双曲波動方程式の解法 [math.NA, cs.NA]目的：双曲波動方程式の数値解析
- 構造工学や振動解析において，高精度な波動伝播のシミュレーションが重要である。
- 既存の数値解法では，精度と計算コストのバランスが課題となっていた。
- ハイブリッド高次法と陰解法の組み合わせにより，高精度かつ安定な解法を確立すること。
- 空間および時間において，最適な収束オーダーが理論的に確認された。
- Newmark法とCrank-Nicolson法，いずれの陰解法も有効であることが数値実験により示された。
- 双曲波動方程式に対するハイブリッド高次法と陰解法の組み合わせに関する最初の研究である。
時間一様弱収束性を持つ減衰Langevin方程式の明示的数値解法 [math.NA, cs.NA, math.PR]目的：減衰Langevin方程式に対する時間一様弱収束性を持つ明示的数値解法の開発
- 統計力学におけるギブス測度のサンプリングや分子動力学シミュレーションの基礎モデルである。
- 明示的解法においては，多項式的に増加するポテンシャル下での時間一様弱収束性が未確立であった。
- 多項式ポテンシャル下での長時間の統計的観測量を正確に再現する効率的な解法を提供すること。
- 提案されたDSAV法は，スカラー補助変数に消散性を導入することで時間一様弱収束性を実現した。
- 精度εに対して，計算コストをO(ε⁻¹ log(ε⁻¹))に大幅に削減可能である。
- Malliavin微積分を用いずに，数値解の密度関数の存在と正定性を確立した。
多角形メッシュ上の線形弾性問題に対する時間差分不連続Galerkin法 [math.NA, cs.NA]目的：線形弾性問題に対する時間差分不連続Galerkin法の開発
- 構造解析において，より複雑な形状を扱う需要が高まっているため，汎用的なメッシュへの対応が重要である。
- 複雑な形状を扱う場合，従来の有限要素法ではメッシュ生成が困難であり，計算コストが増大する可能性がある。
- 多角形メッシュ上で安定かつ高精度な解を得るための数値解法の確立を目指す。
- 本研究で提案する時間差分不連続Galerkin法は，多角形メッシュ上で対称な応力空間を構築し，要素境界における法線連続性を実現する。
- この手法は，ロッキング現象を回避し，線形運動量および角運動量の局所的なバランスを満足することを理論的に証明した。
- 数値実験の結果，本手法の有効性と理論的結果との整合性が確認された。
有限体積法環境におけるメッシュ適応のための等方性回復型誤差推定アルゴリズム - 大気流れへの応用 [math.NA, cs.NA]目的：有限体積法におけるメッシュ適応のための等方性回復型誤差推定
- 大気流れなどの複雑な現象のシミュレーションにおいて，計算精度の向上と計算コスト削減が重要である。
- 有限体積法は広く利用されているが，誤差推定に基づいた効果的なメッシュ適応は未だ十分に研究されていない。
- 本研究は，有限体積法における誤差推定に基づいたメッシュ適応アルゴリズムを開発し，その有効性を検証する。
- 提案する等方性回復型誤差推定アルゴリズムは，数値的不安定性を抑制する効果が認められた。
- 従来のメッシュ適応アルゴリズムと比較して，シミュレーションの大部分において高い精度を示した。
- 適応メッシュを用いることで，固定メッシュと比較して35%から94%の計算コスト削減が実現された。
バイハーモニック・シュテクロフ固有値問題に対する適合および非適合仮想要素法 [cs.RO, math.NA, cs.NA, math.FA]目的：バイハーモニック・シュテクロフ固有値問題に対する仮想要素法の収束解析
- 偏微分方程式の数値解法は，工学や科学における様々な現象のシミュレーションに不可欠である。
- 複雑な形状の領域に対する高精度な数値解法は，依然として課題である。
- 既存手法では固有値のスペクトルに不要な固有値が生じる場合がある。
- 本研究では，離散$H^2$半ノルム，および$H^1$と$L^2$ノルムを用いた収束解析を行った。
- バブシュカ・オズボーンのスペクトル理論により，数値スキームがスペクトルを近似することが示された。
- 固有関数に対する最適な収束次数，および固有値に対する二重の収束次数が得られた。
テンソルネットワーク次元削減による指数規模線形代数 [math.NA, cond-mat.str-el, cs.DS, cs.NA, quant-ph]目的：指数規模線形代数問題の効率的な解法
- 現代科学計算において，次元の呪いに対処することは重要である。問題の規模が大きくなるにつれ，計算量が指数関数的に増加する。
- テンソルネットワークは指数規模のデータを扱う有用なツールだが，必ずしも安定したアルゴリズムに繋がるとは限らない。
- テンソルネットワークデータのランダム次元削減技術により，指数規模線形代数問題を解く。
- 提案手法は，トレース推定や固有値近似といった問題に対し，有効なアルゴリズムを提供する。
- 数値実験では，量子多体系物理学における最大$2^{200}$次元のデータに対して，手法の有効性を確認した。
- ランダム次元削減により，テンソルネットワークの表現を安定化させ，計算効率を向上させることが示された。
正値性保存微分リアプノフ方程式に対する強安定性保存積分因子ルンゲ・クッタ法 [math.NA, cs.NA]目的：微分リアプノフ方程式を解くための積分因子強安定性保存ルンゲ・クッタ法の開発
- 安定性解析や制御理論への応用において，数値解の重要な特性である対称性と正定値半正定値性の維持が不可欠である。
- 従来の数値解法では，対称性と正定値半正定値性の維持と，計算精度向上がトレードオフとなる場合が多い。
- 本研究は，対称性と正定値半正定値性を厳密に維持しつつ，従来の計算精度を上回る解法を開発することを目的とする。
- 提案手法は，数値解の対称性と正定値半正定値を厳密に保存しながら，従来の計算順序の限界を超える。
- 2次スキームに対する厳密な誤差評価を提供し，その精度を理論的に保証する。
- 数値実験により，提案手法の精度と正定値半正定値性の維持能力が検証された。
退化性拡散反応問題に対する最適化された高次IMEX-RKスキーム：移動波現象への応用 [math.NA, cs.NA]目的：退化性拡散反応問題の数値近似
- 生物学や物理学の多様な応用において，反応拡散モデルは重要な役割を果たす。
- 移動波解を持つモデルの数値解法では，急峻な波面を精度良く捉えることが難しい。
- 高次のIMEX-RKスキームにより，移動波の正確な数値シミュレーションを可能にする。
- 提案されたIMEX-RKスキームは，移動波の伝播を正確に捉えることが示された。
- 反応項の分割を回避するSI-IMEX-RKという新しい半陰解法が提案され，剛性の適切な処理が可能となった。
- 高次のポリゴナル不連続ガレルキン法との組み合わせにより，複雑な形状における多スケールダイナミクスの解析を可能にする柔軟かつ堅牢なフレームワークが提供された。
インピーダンス半空間における音響散乱に対するPML-BIE法の収束性 [math.NA, cs.NA]目的：音響散乱問題におけるPML-BIE法の収束性
- 無限領域における波動散乱問題の厳密解は困難であり，数値解析手法の重要性が高い。
- PML-BIE法は有効だが，物理的に重要な設定における厳密な収束理論が未完である。
- PML-BIE法の収束性を確立し，より一般的な散乱問題への応用を可能とすること。
- PML-BIE法の収束性のための一般的なフレームワークを提示した。
- インピーダンス半空間問題において，厳密なPML-BIEと計算可能なPML-BIEの収束性を証明した。
- 最終的な誤差境界はPML吸収電力の増加に伴い指数関数的に減衰する。
PyTorchにおける有限次元I型フォン・ノイマン代数：ランダムブロック対角演算子のGPUアクセラレーションフレームワーク [cs.MS, cs.NA, math.NA, math.OA, math.QA, math.SP, quant-ph]目的：有限次元I型フォン・ノイマン代数に関する数値実験のためのPythonライブラリ
- 量子情報科学や量子統計力学において，フォン・ノイマン代数は重要な数学的基盤である。
- 大規模な演算を行う場合，既存のライブラリではメモリ効率や計算速度が課題となる場合がある。
- GPUを活用し，メモリ効率と計算速度を向上させることで，大規模な数値実験を可能にすること。
- 本ライブラリは，テンソルを用いたコンパクトな表現と遅延評価により，メモリ使用量を削減している。
- 様々な固有値分布やユニタリーアンサンブルに対応したランダム演算子の生成機能を備えている。
- GPUアクセラレーションにより，大規模なモンテカルロスタディ（例：100x100演算子の2万サンプル）を実現している。
2次元非線形弾性体における一次互換ひずみ混合四角形有限要素 [math.NA, cs.NA]目的：2次元非線形弾性体の解析のための，一次互換ひずみ混合四角形有限要素の構築
- 構造解析において，複雑な形状や境界条件を扱う上で，有限要素法は不可欠な手法である。
- 既存の互換ひずみ混合有限要素は三角形要素に限定されており，四角形要素への拡張が課題であった。
- 四角形要素を用いた互換ひずみ混合有限要素を開発し，三角形要素と同等の精度を実現すること。
- 本研究で開発した四角形要素は，圧縮性と非圧縮性の両方の非線形弾性問題に適用可能である。
- 従来の三角形要素と同程度の精度を，より少ない自由度で達成できることが示された。
- 圧力場の要素レベルでの縮約により，非圧縮性解析における追加の自由度を回避することに成功した。
線形偏微分方程式データ同化のためのランダム特徴量カルマンフィルタリング [math.NA, cs.NA]目的：時間依存偏微分方程式のデータ同化におけるベイズ更新
- 気象予測や流体シミュレーションなど，多くの科学技術分野で，時間発展する現象の予測精度向上が求められている。
- 観測データが疎でノイズを含み，状態空間が高次元であるため，効率的なデータ同化手法が課題となっている。
- ランダム特徴量を用いたカルマンフィルタリングにより，高次元状態空間を効率的に扱い，予測精度を向上させる。
- ランダム特徴量を固定化し，線形偏微分方程式をGalerkin離散化することで，有限次元線形ガウス状態空間モデルが得られる。
- 質量正規化された座標系を導入することで，不要な方向を除去し，事後分布の次元を特定できる。
- 熱方程式において，事後分布の収縮と偏微分方程式の一貫性に関する確率的な定理を証明し，不確か量評価の保証を提供した。
ランダウ-リフシッツ-ギルバート方程式における構造保存ソルバーを用いた質量凝縮中心有限要素法の事前誤差解析 [math.NA, cs.NA]目的：ランダウ-リフシッツ-ギルバート方程式に対する質量凝縮中心有限要素スキームの誤差評価
- 強磁性体の磁化ダイナミクスを記述する基礎モデルであり，スピントロニクスなど応用範囲が広い。
- 特にカイラル磁性体におけるDMIの影響は重要だが，数値解法での扱いが難しい。
- 本研究は，DMIを含むLLG方程式に対する質量凝縮中心有限要素法の誤差解析を通して，高精度な数値シミュレーションを可能にする。
- 質量凝縮中心有限要素法がノードの単位長制約を正確に保存し，離散的なエネルギー則を満たすことが示された。
- 適切な正則性仮定の下，エネルギーノルムにおける最適な収束率（空間で一次，時間で二次）が証明された。
- 構造を保存する固定点反復法を提案し，ソルバーの許容誤差の影響を考慮しても同じ収束率を維持できることを示した。
格子グリーン関数を用いた無限格子離散カルデロン射影：アクティブノイズシールドおよび閉じ込め [math.NA, cs.NA, physics.comp-ph]目的：ヘльмホルツ方程式に対する無限格子離散カルデロン射影の構築
- アクティブノイズ制御は，音響環境の改善に不可欠であり，様々な応用分野で重要性が増している。
- 従来の数値手法では，境界条件の設定や計算コストが課題となり，複雑な形状への適用が困難であった。
- 格子グリーン関数を用いることで，これらの問題を解決し，効率的なノイズ制御を実現することを目指す。
- 格子グリーン関数に基づく離散カルデロン射影が，任意の形状の格子上でアクティブノイズシールドと閉じ込めを可能にすることを示した。
- この射影は，発射放射条件を正確に満たし，補助ヘльмホルツ問題を回避することで，数値計算の効率化に貢献する。
- 実験結果は，格子グリーン関数が整合性のある光源に対して機械的精度でのキャンセルを実現し，解析的平面波および点光源に対して2次収束することを確認した。
トポロジースペクトルギャップ，多重スケール・ウェイルの法則，および高コントラスト偏微分方程式の均質化 [cs.DB, math.NA, cs.NA, math.AP]目的：高コントラスト多重スケール偏微分方程式におけるスペクトルギャップと固有値分布の特性
- 高コントラスト問題は，現実世界の様々な物理現象を記述する上で重要である。
- 高コントラストな多重スケール問題では，スペクトルギャップの厳密な位置決定が困難である。
- 高コントラストな設定におけるスペクトルギャップの位置と固有値分布を解析的に決定すること。
- スペクトルギャップの位置は，高コントラストなインクルージョンに関連する局所的な零空間の次元によって普遍的に決定されることが示された。
- 無限次元カーネルを持つ系では，スペクトルギャップの位置はトポロジー的ベッチ数によって厳密に決定される。
- コントラストが無限大に近づくと，多重スケールスペクトルは背景行列のディリクレスペクトルとインクルージョンの内部ノイマンスペクトルに分岐する。
多孔弾性媒体における明示的波動伝播のためのアイソジオメトリック解析 [math.NA, cs.NA]目的：多孔弾性モデルにおけるアイソジオメトリック解析の利点に関する検討
- 波動伝播問題の正確な数値シミュレーションは，地盤工学や石油工学等の分野で不可欠である。
- 従来の有限要素法では，高次解離において偽のモードが発生し，計算コストが増大する問題がある。
- アイソジオメトリック解析を用いることで，偽のモードを抑制し，より効率的な波動伝播シミュレーションを実現すること。
- アイソジオメトリック解析は，多孔弾性体の問題においても，従来の有限要素法が抱える偽のモードの問題を回避できることが示された。
- 適切なスプライン空間を用いることで，安定性と正確な解が得られ，計算効率が向上することが確認された。
- アイソジオメトリック解析を用いることで，時間ステップ幅が解離次数に依存しにくくなり，大規模なシミュレーションが可能となる。
不均質井戸システムの破断解析：ケーシング偏心率とセメント-地層界面強度への影響 [cs.CE, cond-mat.mtrl-sci, cs.NA, math.NA, physics.geo-ph]目的：不均質井戸システムにおける破断の発生と伝播に関する解析
- 油田開発において，井戸の完全性が重要な課題であり，破断予測は安全な生産を確保する上で不可欠である。
- 複雑な地層条件下での破断予測は困難であり，特に不均質性や界面強度の影響評価が課題となっている。
- ケーシング偏心率やセメント-地層界面強度が生じる破断のメカニズムを解明し，井戸の信頼性向上に貢献すること。
- ケーシング偏心率は破断開始圧力に大きく影響し，同心円配置と比較して最大30%低減することが示された。
- 偏心率が50%を超えると，放射状の破断に加えて傾斜した破断が発生し，新たな破壊モードが確認された。
- 弱い界面では，セメント中の破断が地層界面に沿って逸れることで応力緩和が遅延し，追加の破断発生リスクが増加する。
グラスマン多様体上の測地線補間：GLERPと再帰的GIDER補間 [math.NA, cs.NA]目的：グラスマン多様体上のデータに対する測地線補間手法
- 多様体上のデータ処理は，データの幾何学的構造を考慮する上で重要である。特に，部分空間を扱う場合，その構造を適切に捉える必要がある。
- 既存の補間手法は，基底の曖昧性に対して不変ではなく，データの本質的な構造を反映しない可能性がある。
- 基底の曖昧性に対して不変な，より正確なグラスマン多様体上の補間手法を開発すること。
- GLERPは，グラスマン多様体上の球面線形補間の類似手法であり，指数写像と対数写像を用いて定義される。
- 再帰的な高次構成であるGIDER_nは，ネヴィルのアルゴリズムの各アフィン補間ステップをGLERPで置き換えることで得られる。
- 数値実験の結果，GIDER_nと接空間補間は滑らかな局所領域でほぼ区別がつかず，投影行列補間は外的な基準となることが示された。
特異摂動問題に対するPetrov-Galerkin変分物理情報ニューラルネットワークフレームワーク [math.NA, cs.LG, cs.NA]目的：二次元特異摂動問題の効率的な解法
- 工学・科学における境界層問題の正確な解析は，現象の理解と予測に不可欠である。
- 特異摂動問題は，多重解尺度を持つため，従来の数値解法では高精度な解を得るのが困難である。
- 本研究は，多重解尺度問題を克服し，高精度な数値解を効率的に得ることを目指す。
- 提案手法は，Petrov-Galerkin変分法と物理情報ニューラルネットワークを組み合わせることで，高い精度を実現した。
- 二次元特異摂動問題において，最大誤差およびL_2ノルムの両方で高い精度が確認された。
- 本手法は，多重解尺度特徴を捉える上で効率性とロバスト性を示すことが明らかになった。
完全エリックセン-レスリーモデルに対する幾何学的再構成と構造保存型回転離散勾配スキーム [math.NA, cs.NA]目的：完全エリックセン-レスリーモデルに対する構造保存型回転離散勾配スキームの構築
- 液晶材料は表示デバイス等に応用され，その数値シミュレーションは重要である。
- 液晶流れの数値計算では，ディレクターの単位長制約とエネルギー消散の維持が課題である。
- ディレクターの単位長制約を保ち，離散レベルでのエネルギー消散を確実に実現すること。
- 本研究では，エリックセン-レスリー系を回転形式に再構成することで，ディレクターの進化が単位球に接するようにした。
- 提案手法はディレクター長を保存し，無条件の離散エネルギー法則を満たすことが示された。
- 数値実験により，提案手法の精度，単位長制約の保持，エネルギー安定性が検証された。
放熱星上歩行法：ノイマン境界条件を持つ放物型熱方程式 [cs.CL, cs.GR, cs.NA, math.NA]目的：放物型熱方程式に対するモンテカルロ解法
- 偏微分方程式の数値解法は，科学技術計算の根幹であり，効率的な手法が求められている。
- 時間発展を伴う放物型方程式に対するモンテカルロ法の適用は，空間と時間の結合により困難であった。
- 時間変化する領域を扱うための枠組みを構築し，モンテカルロ法を放物型方程式に適用可能とする。
- 新たなモンテカルロ解法「放熱星上歩行法」を提案し，境界積分の枠組みを放物型方程式に拡張した。
- 時間座標と空間方向の対数パラメータ化により，二重層カーネルが独立なガンマ関数と一様成分に分解されることを示した。
- 空間微分を境界積分で表現する解離型勾配推定器を導出し，空間時間領域にヘテロスカダスティック回帰ベースのノイズ除去を適用した。
PhysGuard：フィッシャー情報に基づく勾配投影によるシミュレーションから現実世界へのニューラル偏微分方程式サロゲート [cs.LG, cs.NA, math.NA, physics.comp-ph]目的：シミュレーションデータで訓練されたニューラル演算モデルの，現実の実験測定値への適用における精度損失の軽減
- シミュレーションデータで学習したモデルを現実世界に適用する際，性能劣化が課題となる
- 限られた現実データでの微調整は有効だが，学習済みの物理法則に関する表現を損なう可能性がある
- 物理法則を保持しつつ，シミュレーションから現実世界への適応を可能にする手法の開発
- PhysGuardは，シミュレーションデータ上のフィッシャー情報行列を用いて，物理的に重要なパラメータ方向を特定する
- 微調整時の更新を，これらの重要な方向に干渉しない方向に制限することで，物理法則の保持を実現する
- ベンチマーク実験により，特にドメインシフトが激しい状況下で，標準的な微調整と比較して低周波誤差を最大32％削減できることが示された
カーネルに基づく多水準補間のノード表現 [cs.IR, math.NA, cs.NA]目的：散在データからの多変数関数の近似または学習
- 疎グリッド法など，様々な応用分野で重要性が増している。
- 補間演算子の範囲や基底関数の性質が十分解明されていない。
- 多水準補間演算子の範囲を特徴付け，基底関数を構築すること。
- ネストされた点集合に対して，ターゲット関数の値に依存するノード表現を導出した。
- 得られた基底関数は指数関数的な減衰を示すことが証明された。
- 非ネストされた点集合に対しても，一般化されたノード表現を導出した。
CREATORプロジェクト：計算電気機械研究所へ向けて [cs.CE, cs.NA, math.NA]目的：電気機械のシミュレーション駆動設計のための新たなパラダイムの確立
- 効率，電力密度，持続可能性への要求が高まり，電気機械の設計における重要性が増している。
- 従来の逐次的な設計ワークフローでは，多物理現象，高度な材料，複雑な形状の統合に対応できない。
- 電気工学，応用数学，流体工学，材料科学の専門知識を統合し，設計課題を解決することを目指す。
- 本研究は，シミュレーション駆動設計のための統合モデリング，シミュレーション，最適化手法の確立を目指す。
- 第一期（2022-2026年）の研究成果と，計算電気機械研究所に向けた現状の開発状況を概説する。
- ドイツとオーストリアの国家資金提供機関からの大規模プロジェクトとして，学際的な専門知識を組み合わせている。
混合境界条件を持つ楕円型問題に対するFEMにおける数値積分効果 [math.NA, cs.NA]目的：混合境界条件を持つ楕円型問題におけるFEMの数値積分の精度に関する研究
- 有限要素法は，工学や科学における様々な問題を数値的に解くための強力な手法である。
- 高次の要素を用いる場合，数値積分精度が計算結果に大きく影響することが知られている。
- 適切な数値積分精度を決定することで，高次要素の性能を最大限に引き出すことを目指す。
- 体積積分では少なくとも次数$2p-2$，表面積分では少なくとも次数$2p-1$の代数的精度を持つ数値積分を用いることで，$H^1$収束が最適に保たれる。
- $L^2$誤差に関しては，体積積分で次数$\max\{p,2p-2\}$以上，表面積分で次数$2p-1$以上の数値積分を用いることで，$O(h^{p+1})$の最適収束が達成される。
- 線形要素（$p=1$）の場合，得られた$L^2$誤差に関する結果が鋭いことが示され，ディリクレ境界条件に限定すると，既存の結果と比較してデータ依存性が改善される。
拘束されたヘルフレッヒ流れに対する幾何学的構造を保存するパラメータ有限要素近似 [cs.RO, math.NA, cs.NA]目的：閉曲線および曲面の拘束されたヘルフレッヒ流れに対する構造保存パラメータ有限要素法
- 細胞膜や液晶など，物理現象のモデル化に重要である。形状変化を精度良く捉える必要性が高い。
- 既存手法では，エネルギー保存や形状保存が十分でなく，計算不安定や精度低下を引き起こすことがある。
- エネルギー散逸と幾何学的保存を両立させ，より安定かつ高精度な数値解法を開発することを目指す。
- 提案手法は，エネルギーの減衰と正確な幾何学的保存を同時に達成する。
- 二段階の速度分割戦略を用いることで，線形システムと非線形システムを効率的に解くことが可能である。
- 数値実験により，初期データや自発的曲率に対する精度，ロバスト性，構造保存特性が確認された。
任意断面を持つ梁に対する，零空間法を用いた可変オフセット関節定式化 [math.NA, cs.NA]目的：任意断面を持つ梁の接続に関する定式化
- 構造設計において，複雑な形状の梁構造の正確な解析は不可欠である。
- 従来の定式化では，回転自由度への明示的なインターフェースが必要で，汎用性に欠ける。
- 異なる梁の相対位置を拘束・解放する，柔軟な接続方法を提案する。
- 本研究では，梁同士を接続する局所的な構成制約の変分定式化を提示した。
- 零空間法を用いることで，システム行列のサイズ縮小と条件数の改善を実現した。
- 提案手法を拡張可能なディレクターに基づく梁要素と，非線形弾塑性問題に適用し，数値例で検証した。
PF-DIC：相場法デジタル画像相関を用いた統合的な全場変位，ひずみ，損傷測定 [math.NA, cond-mat.mtrl-sci, cs.NA, physics.optics]目的：変位，ひずみ，損傷の全場測定のための新規デジタル画像相関フレームワーク
- 材料の強度や寿命を予測する上で，変位，ひずみ，損傷の正確な測定が不可欠である。
- 従来のDIC法では，複雑な損傷やき裂に対して正確な測定が困難であり，マスク処理に手間がかかる。
- 相場法をDICに導入することで，損傷の自動測定と測定精度の向上を目指す。
- 提案手法PF-DICは，様々な種類のき裂を捉える能力を示し，変位場の測定精度を向上させる。
- PF-DICは，特定の荷重条件における重要き裂を選択的に識別し，損傷評価に活用できる。
- 本手法は，相場法シミュレーションと実験測定の統合を可能にし，材料欠陥の特性評価や構造健全性監視に貢献する。
ラプラス固有値問題に対するNitsche法に基づく有限要素法：スペクトル近似と事後誤差解析 [math.NA, cs.NA]目的：ラプラス固有値問題の数値解析
- 工学・科学における様々な現象のモデル化に不可欠であり，高精度な解析が求められる。
- 境界条件の扱いが難しく，数値計算の安定性や精度を損なう可能性がある。
- Nitsche法を用いて境界条件を弱く課すことによる数値解の精度向上を目指す。
- コンパクト作用素理論の枠組みで離散固有値問題を解析し，離散解作用素のノルム収束性を示す。
- 固有値および固有関数の誤差評価を導出し，Nitsche法の変種に応じた収束率を明らかにする。
- 事後誤差解析を開発し，適応的なメッシュ細分化に適した残差ベースの誤差評価子を提案する。
ハイゼンベルクスピンガラスハミルトニアンに対する量子位相推定のための初期状態準備としてのVQE [quant-ph, cond-mat.dis-nn, cs.NA, math.NA]目的：ハイゼンベルクスピンガラス量子ハミルトニアンの初期状態準備におけるVQEの有効性
- 量子コンピュータによる量子ハミルトニアンの基底状態エネルギー計算は，科学的・商業的価値が期待される最も有望な応用分野である。
- 量子位相推定（QPE）は，高品質な初期状態を必要とするため，その準備がボトルネックとなっている。
- VQEが初期状態準備アルゴリズムとして有効か検証し，QPEへの応用可能性を探る。
- VQEは，研究対象のハミルトニアンに対して一般的に基底状態へ効率的に収束しないことが確認された。
- 低エネルギー状態と基底状態との重複度は必ずしも相関しないものの，ある程度の相関は見られた。
- VQEの変分回路の層数を増やしても，重複度やエネルギーの改善は見られず，VQEがQPEのためのヒューリスティックな初期状態準備アルゴリズムとして潜在的に利用できる可能性が示唆された。
確率論，数論，力学系の交差点における数値積分 [math.DS, cs.NA, math.NA]目的：一様分布数列の性質
- 確率論，数論，力学系は数学の基礎であり，応用範囲も広いため重要である。
- 一般的な非コンパクトトポロジー空間における一様分布の判定が難しい。
- ポーランド空間における一様分布数列の収束判定クラスの存在を示す。
- ポーランド空間において，正則な確率測度を持つ場合，一様分布数列に対する可算な収束判定クラスの存在が証明された。
- 可算な収束判定クラスが存在する条件下で，独立同一分布数列とエルゴード変換の軌道数列が一様分布を示すことが示された。
時間依存偏微分方程式に対する液体ランダム特徴手法 [physics.comp-ph, cs.NA, math.NA]目的：時間依存偏微分方程式のメッシュフリー時空間近似における，時間スケールの進化表現と計算コストの抑制
- 偏微分方程式は自然現象や工学問題を記述する基盤技術であり，その高精度な数値解法は不可欠である。
- 従来のメッシュフリー手法では，時間スケールの変化に対応できず，計算精度が制限される場合がある。
- 液体ランダム特徴手法を用いて，時間スケールの表現能力を高め，計算精度を向上させる。
- 液体ランダム特徴手法は，時間活性化関数を液体時間定数応答に置き換えることで，時間スケールの情報を埋め込む。
- この手法は，剛性，分散性，多重スケール性を持つ問題において，従来の近似手法よりも高い精度を示す。
- 時間スケールを凍結された試関数に直接組み込むことで，線形最小二乗ソルバーの簡便性を維持しつつ，高精度な時空間代替モデルを実現する。
傾斜入射によるインピーダンス円筒の結合境界積分方程式に対する高次ナイストローム法 [math-ph, cs.NA, math.MP, math.NA]目的：インピーダンス円筒による電磁波散乱の数値解法
- 電磁波散乱問題は，レーダー技術，アンテナ設計，電磁両立性など広範な分野で重要である。
- インピーダンス円筒の結合境界積分方程式は，特異積分や導関数を含むため，高精度な数値解法が困難である。
- 高次ナイストローム法を用いることで，インピーダンス円筒の電磁波散乱問題を安定かつ高精度に解くことを目指す。
- Kress型対数カーネル分解，周期積積分，フーリエ微分，ブロック対角事前条件付け子を用いた高次ナイストローム法を構築した。
- 連続解の一意性および離散安定性仮定の下，境界密度と遠場パターンの高次収束フレームワークを定式化した。
- 数値実験の結果，本方法は結合インピーダンスシステムに対する安定した高精度な前方ソルバーとして機能することが示された。
テンソルトレインランダムベクトルによる確率的トレース推定 [stat.ML, cs.LG, cs.NA, math.NA]目的：大規模行列のトレース近似手法
- 機械学習や科学計算において，大規模行列のトレースは重要な計算量
- テンソル構造を持つ場合，既存手法では計算コストやサンプル数が課題
- テンソルトレインベクトルを用いて，効率的なトレース推定を目指す
- テンソルトレインベクトルは，適切なランク選択により，古典的な推定器と同等の次元に依存しない保証を提供する。
- ランク$r \geq d-1$のテンソルトレインベクトルを用いた中央値平均法は，非構造化ガウスベクトルに基づく推定器と同等の精度と信頼性を持つ。
- テンソルトレインベクトルを用いたスケッチは，Nystr\"{o}m++フレームワークにおいて，$\mathcal{O}(\varepsilon^{-1})$のサンプル複雑度を達成可能である。
確率微分方程式に対するDLRAの指数的収束 [math.PR, cs.NA, math.NA]目的：確率微分方程式の動的直交近似の長時間の振る舞い
- 確率過程の数値解法は，科学技術におけるシミュレーションの基盤である。
- 高次元の確率微分方程式の解法は計算コストが高く，困難を伴う。
- 低ランク近似を用いることで，計算量を削減し，長時間の統計的性質を近似する。
- 強DOシステムがwell-posedであり，Wasserstein距離における誤差評価が確立された。
- 適切な条件の下で，強DOシステムの不変確率測度が存在することが証明された。
- 低ランク過程も不変確率測度を持つことが示され，不変測度の構造が例を通して議論された。
対称性分解されたエルミート行列の固有値厳密包含 [math.CO, cs.DM, cond-mat.str-el, cs.NA, math.NA, physics.comp-ph]目的：対称性分解されたエルミート行列の固有値の厳密な包含範囲の算出
- 量子多体系計算において，固有値は系の性質を決定する重要な情報である。その精度保証は不可欠。
- 浮動小数点演算の誤差や計算アルゴリズムの特性により，計算された固有値の精度が保証されない場合がある。
- 数値誤差を解析的に評価し，真の固有値との距離を厳密に評価できる固有値包含法を提供する。
- 本研究で開発された\texttt{symveig}は，NumPy/SciPyのみを用いて，エルミート行列のすべての固有値の厳密な包含範囲を計算する。
- 対称性を利用した分割を行うことで，計算速度を最大130倍，包含範囲の精度を3～9倍向上させることに成功した。
- 1次元および2次元ハイスンベルク模型での検証により，計算された固有値が常に包含範囲内に収まることが確認された。
ヘシアン変動下におけるチェビシェフ加速法：サイン・ヤコビ法 [math.OC, cs.NA, math.NA]目的：ヘシアン変動下における有限水平線一次勾配法の実装
- 最適化問題において，収束速度の向上は常に重要な課題である。
- ヘシアンが時間とともに変化する場合，最適化アルゴリズムの安定性や効率が低下する。
- ヘシアン変動下におけるチェビシェフ加速法の性能向上を目指す。
- チェビシェフ再帰法における前置精度でのゲインの鋭さを示す。
- ヤコビ座標を用いることで，終端のみが正確な手法を実現可能とする。
- 実験結果から，有限水平線効果やより低い再起動コストが確認された。
線形演算子の最適マルチスケール学習 [math.ST, cs.NA, math.NA, stat.ML, stat.TH]目的：線形演算子の学習における統計的・計算的限界
- 機械学習において，演算子の学習は重要な課題であり，様々な応用が期待されている。
- 高次元データに対する効率的な学習方法が確立されておらず，計算コストが課題となっている。
- マルチスケール解析を利用し，最適な学習コストで高精度な演算子推定を実現する。
- ウェーブレット座標系を用いて，問題をマルチスケール構造を持つ無限次元回帰問題として再定式化した。
- ソボレフ演算子ノルム損失の下でミニマックスレートを確立し，それを達成する有限解像度ブロックごとの最小二乗推定法を構築した。
- スケール適応的サンプルサイズを用いることで，最適な計算コストを実現できることを示した。
偏微分方程式に対する多段階深層学習：バーガース方程式への応用 [math.NA, cs.AI, cs.NA]目的：偏微分方程式の解法における多段階深層学習法の開発
- 偏微分方程式は科学技術の様々な分野で現れる基礎方程式であり，その効率的な解法が重要である。
- 深層ニューラルネットワークを用いた解法は最適化が困難であり，特に非線形な偏微分方程式では課題が多い。
- 本研究は，最適化の複雑さを軽減し，安定した階層的な精度向上を実現する手法を提案する。
- 提案手法であるTS-MGDLは，まず浅いネットワークを段階的に学習させ，低周波から高周波成分へと近似していく。
- 次に，学習済みのネットワークを初期値として，選択された層を再学習することで，解の精度を向上させる。
- 数値実験の結果，TS-MGDLは従来の単一段階学習よりも最大60倍高い精度でバーガース方程式を解くことができた。
分数ラプラシアン問題に対する直交貪欲アルゴリズムfOGA [math.NA, cs.NA]目的：分数ラプラシアン方程式の数値解法
- 偏微分方程式は，物理現象や工学問題を記述する上で不可欠である。
- 分数ラプラシアン方程式の効率的な数値解法は未だ確立されていない。
- 直交貪欲アルゴリズムを用いて，高精度かつ高速な数値解法を構築する。
- 有限差分法と浅いニューラルネットワーク近似を組み合わせた数値解法を提案した。
- 分数ラプラシアン演算子は，リーマン-リウヴィル型の方向表現を用いて離散化された。
- 二次元問題では，数値積分法と補助点を用いることで，演算子の評価を容易にした。
行列近接問題と固有値最適化 [math.NA, cs.NA, math.DS, math.OC]目的：行列近接問題の解法
- 力学系やグラフ理論など，多様な分野で現れる重要な問題である。
- 固有値や特異値を目的の位置に移動させるための最適な摂動が不明である。
- ランク1の摂動を用いて，効率的な解法を開発する。
- 最適な摂動行列は，通常ランク1であるか，ランク1行列の射影であることが示された。
- 内側の反復では，勾配法に基づくランク1行列微分方程式によって固有値最適化問題を解く。
- 外側の反復では，スカラー非線形方程式を解くことで最適な摂動サイズを決定する。
ひずみ勾配弾性問題に対する最適かつロバストな非適合有限要素法 [math.NA, cs.NA]目的：ひずみ勾配弾性モデルの解析手法
- 材料の微視的な構造がマクロな挙動に影響を与える解析において重要である。
- ひずみ勾配弾性問題の数値解析には，高次の要素が必要とされ，計算コストが高い。
- 低次要素を用いた，最適かつロバストな数値解法を提案することで，計算コストを削減する。
- 提案された非適合有限要素法は，ラメ定数やサイズパラメータに対して最適かつロバストであることが数値実験により確認された。
- 二次元および三次元の滑らかなストークス複体の非適合有限要素離散化も新たに考案された。
- 1
- 2