| 研究課題/領域番号 |
18K03438
|
|---|
| 研究機関 | 岡山理科大学 |
|---|
研究代表者 |
大江 貴司 岡山理科大学, 理学部, 教授 (90258210)
|
|---|
| 研究分担者 |
町田 学 近畿大学, 工学部, 准教授 (40396916)
|
|---|
| 研究期間 (年度) |
2018-04-01 – 2024-03-31
|
| キーワード | 異常拡散現象 / 数値解法 / 代用電荷法 / 基本解 / CQM |
|---|
| 研究実績の概要 |
令和4年度は主として、拡散方程式に対する代用電荷法のNaiveな実装において時間刻みを小さくした場合に生ずる数値的不安定性への対処について、研究を進めた。前年度末に、Convolution Quadrature Method (CQM) を適用することにより安定性の改善の可能性があるという着想を得た。この着想を確認すべく、予備的な数値実験を行った結果、数値的に安定に計算を行うことができたたため、その理由を探るべくいくつかの観点から検討を行った。まず、ステップ毎の仮想電荷の計算に必要となる連立方程式の性質、特に行列の条件数を調べたところ、大幅に小さくなっていることが確認できた。このことが数値的不安定性の除去につながったものと考えられる。一方でCQMを用いるためには、近似解の表現として時間方向についてはシンプルな重ね合わせではなく畳み込み積分を離散化する手法を用いている。そのため、得られた数値解が拡散方程式を厳密に満たすという性質を失っている可能性がある。また、メッシュレス的な方法が持っている、時間・空間に関する計算点の自由度がかなり失われており、これらの点について改善の余地があることが判明した。
なお先に述べたように、代用電荷法のNaiveな実装において仮想電荷の計算に必要となる連立方程式の行列は悪条件なものであるが、対角性は強いため、その性質を生かした安定化が図れないかについても検討中である。
以上の結果について、令和4年度中の発表はかなわなかったものの、5年度の早い時期に行われる学会で発表を行う。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
4: 遅れている
理由
令和4年度は、令和3年度の研究で判明した拡散方程式に対する代用電荷法のNaiveな実装において、短い時間刻みの場合に生ずる不安定性に対する対処を中心に研究を進めた。その結果、一定の成果を得ることはできたものの、まだまだ不十分なものであると考えている。
また、目標である層状領域に対する適用については、令和4年度は手を付けることができなかった。そのため、研究期間をさらに1年延長することで見込みを得たい。
以上を考慮し、当初の研究予定からは「遅れている」と判断したが、研究計画の段階では知りえなかった興味深い数理的現象が数多く発見された。これらの発見は、非定常現象に対する代用電荷法の適用の研究の今後の発展に寄与するものと考えている。
|
| 今後の研究の推進方策 |
令和4年度の研究で、時間刻みを小さくしたときに生ずる数値的不安定性に対し、Convolution Quadrature Method (CQM) が有効であることがわかった。また、線型安定解析の結果として、仮想電荷の計算で必要となる連立方程式にあらわれる行列について、その条件数を大幅に小さくすることができることが分かった。
これらの結果を踏まえ、令和5年度は最終年度であることを考慮し、次の２つのテーマに絞って研究を進める。(i) CQM を適用した手法の高精度化、特にCQMに陰的 Runge-Kutta 法を組み合わせることにより、線型多段解法より高精度な計算の可能性を追求する。(ii) 層状の領域に対する適用の検討、特に空間的な電荷点および拘束点の配置に関する分析を行う。
|
| 次年度使用額が生じた理由 |
令和4年度は、過去2年間に比べ新型コロナウイルスの感染拡大の影響も小さくなり、国内の研究集会・学会の一部は対面で実施されるようになった。しかし、このような状態になったのは年度の後半であったため、学会出席等のために用意した旅費が未使用で残ることとなった。
令和5年度は年度初めから、ほとんどの学会・研究集会が対面で行われる見込みである。この状況を鑑み、延長した費用については主として国内学会発表・出席の旅費として使用する予定である。また、分担者との研究打ち合わせ用にも使用する予定である。
|
| 備考 |
・日本応用数理学会環瀬戸内応用数理研究部会 第26回シンポジウム
https://sites.google.com/view/kanseto-jsiam-2022/
本研究課題に関連して、上記の研究集会を主催者・研究代表者として実施した。
|
