2022 Fiscal Year Research-status Report
偏微分方程式の解に対する数値的検証法の新たな高度化の研究
| Project/Area Number |
21K03378
|
|---|
| Research Institution | Waseda University |
|---|
Principal Investigator |
中尾 充宏 早稲田大学, 理工学術院, その他(招聘研究員) (10136418)
|
|---|
| Project Period (FY) |
2021-04-01 – 2024-03-31
|
| Keywords | 数値解析 / 精度保証付き数値計算 / 有限要素法の構成的誤差評価 / 解の数値的検証法 / 計算機援用証明 |
|---|
| Outline of Annual Research Achievements |
本研究課題と関連した研究協力者との情報交換および研究連携のもとにを検討を行った。本年度得られた研究結果は主として次の通りである。
（１）発展型方程式の基本形である熱方程式に対し、その時間周期解に対する従来の全離散解の構成的誤差評価を改良した。具体的なモデル問題に対する数値例を実装し、十分な改良効果が実現されていることを確認した。（２）3次元Stokes 方程式の有限要素解に対するa priori誤差評価定数を用いて、3 次元定常Navier-Stokes 方程式の解に対する数値的検証法を定式化し、非凸多面体領域上の問題に適用して、その検証数値例により、有効性を実証した。これは世界的に見ても未だ例のない画期的な成果であり、国際学術誌に論文として公表した。（３）藤田型の発展方程式に対し、解の爆発可能性および時刻の特定に関して、それらを精度保証付きで求める手法を定式化し、それを用いて空間1次元の藤田型問題に対する爆発時刻包み込みの実例を示し、その有効性を立証した。これは従来の数値的近似とは異なり、解が有限時間爆発することおよびその時刻を、数学的に厳密に保証するものである。（４）線形楕円型作用素に対する近似逆作用素ノルムの厳密な逆作用素ノルムへの収束について考察しその条件を明らかにした。さらに収束オーダーについてもいくつかの可能性を示唆した。楕円型作用素は無限次元であり、それを有限次元で近似した作用素が作用素として収束することは期待できないが、そのノルムは収束することを明らかにするものである。（５）熱方程式の全離散解に対して、その構成的a priori誤差評価定数を計算機援用証明を用いずに導出した。従来のa priori誤差評価定数の算定では、行列固有値問題の解を精度保証付きで求める必要があり、定数の決定が空間領域の形状や個々の離散化スキームに依存し、非効率的であったのに対して、この結果は本質的な改善をもたらすものである。
|
| Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
熱方程式の全離散解に対し、その構成的a priori誤差評価定数を行列固有値問題の精度保証計算を用いずに導出する方法を与えたことは、従来の算定法に対して画期的な改良をもたらすものである。この結果によって、誤差評価の算定が空間領域の形状や個々の離散化スキームに依存せず容易に手計算での算定が可能となった。さらに、従来手法では空間高次元の場合、対応する固有値問題に現れる行列の次数は空間次元のべき乗に比例するため、空間多次元では、その計算コストは莫大なものとなる。したがって本結果は、放物型問題の計算機援用証明にも大きなブレークスルーを与えるものと期待される。また、３次元定常Navier-Stokes 方程式の解に対する数値的検証の実現や、藤田型方程式の解が有限時間爆発し、その時刻を数学的に厳密に特定できることを示した結果なども、前年度の成果をさらに発展させるものである。さらに線形楕円型作用素に対する近似逆作用素ノルムの厳密な逆作用素ノルムへの収束性証明とその収束オーダーの評価も、これまで解明されていない成果である。これらの研究成果は、本研究課題が当初の計画通り順調に進展していることを示すものである。
|
| Strategy for Future Research Activity |
これまでと同様に、本課題の関係研究者との緊密な情報交換と研究協力の下に、次の課題に対して恒常的に検討を進める。
（１）非線形楕円型方程式の解のNewton法による数値的検証の効率化に向けて、線形楕円型作用素に対する逆作用素の存在検証とそのノルム評価のさらなる改良を図る。（２）熱方程式の全離散近似に対する構成的誤差評価法の改良と、それを用いた非線形発展方程式の新たな数値的検証法の定式化とその実装を図る。（３）楕円型作用素の近似逆作用素ノルムの厳密な逆作用素ノルムへの収束オーダーに関する検討（４）無限領域における解の検証条件とその数値的検証法の定式化（５）高Reynolds領域における、3次元Navier-Stokes 方程式の解に対する数値的検証の実現（６）発展方程式の解の爆発に関して、初期値との関連性や爆発領域の特定を精度保証する手法の開発（７）熱方程式の周期解に対し、空間方向有限要素、時間方向スペクトル法にもとづく全離散近似に対する構成的誤差評価を与え、それを用いた非線形放物型方程式の周期解の存在検証を実現する（８）非線形楕円型問題に対する全解探索法を定式化し、その実用化を図る。
|
| Causes of Carryover |
新型コロナウィルス感染症による影響のため、予定した研究出張が取りやめとなったことが主な理由である。今年度は延期となっていた国内外の研究集会の開催等により、出張費が大幅に増大するため、前年度分を含めて全額使用の予定である。
|
Research Products
(11 results)
