2023 Fiscal Year Research-status Report
能動的プラズモンデバイス設計向け熱―弾性波―電磁界時間領域複合物理計算手法の開発
| Project/Area Number |
21K17753
|
|---|
| Research Institution | Nihon University |
|---|
Principal Investigator |
岸本 誠也 日本大学, 理工学部, 助教 (90843053)
|
|---|
| Project Period (FY) |
2021-04-01 – 2025-03-31
|
| Keywords | 電磁界解析 / 複合物理解析 / 電磁波 / 有限差分法 / FILT / マルチフィジックス / 時間領域 |
|---|
| Outline of Annual Research Achievements |
昨年度までの検討にて明らかになった課題として、電磁界―熱―弾性波の複合物理解析を行うと、計算量が想定以上に増大することが見込まれている。この問題を解決するため、２０２３年度は、①磁化と電磁界の連成解析法の検討と高速化手法の提案、②FILT法による大規模電磁界解析法の開発、を中心に検討を行った。これらのテーマに関する研究実績の概要は以下の通りである。
①磁化と電磁界の連成解析法の検討と高速化手法の提案：研究計画立案時に、想定上の計算負荷が生じる場合は、熱、弾性波以外の他の物理現象と電磁界の相互作用を考慮することを策定していたため、磁化の運動と電磁界の複合物理解析について検討を行った。磁化の運動に対する支配方程式であるLLG方程式と電磁界の支配方程式であるMaxwell方程式を連立し時間発展を計算する方法を開発した。特に、昨年度までの検討で課題となった計算負荷の削減を行うため、電磁界と磁化の分割サイズを変えてシミュレーションを行うマルチスケール解析法を提案した。提案法により既存の電磁界―磁化複合物理解析法の計算時間を分割領域のサイズ比と同等に高速化できることを確認した。
②FILT法による大規模電磁界解析法の開発：本研究課題で開発する計算手法である数値逆ラプラス変換（FILT : Fast inverse Laplace transform）法を用いた計算時間の削減を検討するため、散乱体が入射波長に対し数百倍となる問題を解析する方法を開発した。この方法にて、FILT法の計算量を削減する方法を検討し、計算精度の点から、計算時間に関わる無限級数の打ち切り項数の上限があることを確認し必要以上に計算をすること必要がないことを示した。
|
| Current Status of Research Progress |
Current Status of Research Progress
4: Progress in research has been delayed.
Reason
２０２３年度においては、以下の得られた研究成果と取り組みにより、本研究は遅れていると判断する。
①FILT法による熱および弾性波の数値シミュレーション方法に関する研究：FILT法を用いた熱伝導の数値シミュレーション方法について、電磁界の逐次計算とFILT法を連成する方法の確立が進んでいない。FILT法を用いた弾性波の数値シミュレーション方法について、固体中を伝搬する波については計算妥当性を検証するための逐次計算のみ計算コード作成が完了している状況である。更にFILT法を用いずに逐次計算のみで電磁界―熱―弾性波のシミュレーションを実施したところ、想定以上に計算時間が増大し、十分な検討が困難であることがわかったため，逐次計算でも計算量の削減が課題となっている。
②熱―電磁界を連成する方法：電磁界の逐次計算を行わず、特定の観測時間までの消費エネルギー等を扱う方法とFILT法による熱および弾性波の数値シミュレーション方法に関する研究を連成する検討を実施する予定であったが、開発法を十分に検証するための方法が確立できていない。また、検討中に従来法の範疇である熱―電磁界の逐次計算部においても、計算時間削減に余地があることがわかってきたため、この検討を先に行う必要がある。
本研究課題のベースとなるFILT法による熱および弾性波の数値シミュレーション方法に関する研究は目標を達成できているが、項目①②それぞれにおいて開発法を十分に検証るための方法の確立が遅れている。このため、総合的に現在までの達成度は遅れていると判断する。
|
| Strategy for Future Research Activity |
今後の研究推進方策として、本研究の根幹である電磁界の逐次計算を行わず、特定の観測時間までの消費エネルギー等を扱う方法を用いた複合物理解析の検証及び逐次計算における計算量の削減の検討を中心に、以下の研究項目・計画に従い実施する。
①FILT法を用いた電磁界解析と熱解析法の連成シミュレーション方法の確立：検討が終了している特定の観測時間までの消費エネルギー等を扱う方法と熱伝導解析法の連成について検討を行う。定式化は完了しているため、まずは電磁界と熱伝導解析を逐次計算する方法について、特定時間で相互作用を計算することを中心に検討を行い、提案法の有効性を確認する。計算妥当性については電磁界が連続的に変化する場合をモチーフとして扱い、提案法の有効性については、電磁界が急峻に変化する場合や、パルス波の応答について検証を行う。
②電磁界解析法の計算負荷を削減する方法の検討：複合物理解析の計算時間が想定以上に増大しているため、電磁界解析法の計算時間を短縮する方法について検討を行う。FILT法を用いた時分割並列計算の併用、逐次計算負荷をマルチスケールに変更する計算法の検討、圧縮センシングに基づく行列計算の削減を中心に検討する。
|
| Causes of Carryover |
２０２１年度、２０２２年度に計上していた研究成果の周知と情報収集のために計画していた国際会議への出張を感染症対策として取りやめたため、出張旅費分に差額が生じている。また「その他」の枠に論文掲載費としてＩＥＥＥオープンアクセス誌の額を基に計上していたが、査読結果や、論文採択までの時間を鑑みて、学術雑誌を変更し論文掲載費が減額したことから差額が生じている。
次年度以降の予算使用計画として２０２３年度から投稿準備を行っている論文に対して、英文校閲費、論文投稿費を充てる。また、参加を見送ってきた国際会議の出張旅費および、所属研究室に在籍している研究協力者の学会発表における参加費として予算を使用予定である。
|
