研究課題/領域番号 |
23K28103
|
補助金の研究課題番号 |
23H03413 (2023)
|
研究種目 |
基盤研究(B)
|
配分区分 | 基金 (2024) 補助金 (2023) |
応募区分 | 一般 |
審査区分 |
小区分60100:計算科学関連
|
研究機関 | 慶應義塾大学 |
研究代表者 |
飯盛 浩司 慶應義塾大学, 理工学部(矢上), 准教授 (50638773)
|
研究分担者 |
松島 慶 東京大学, 大学院工学系研究科(工学部), 助教 (60943537)
|
研究期間 (年度) |
2023-04-01 – 2027-03-31
|
研究課題ステータス |
交付 (2024年度)
|
配分額 *注記 |
16,900千円 (直接経費: 13,000千円、間接経費: 3,900千円)
2026年度: 1,430千円 (直接経費: 1,100千円、間接経費: 330千円)
2025年度: 1,430千円 (直接経費: 1,100千円、間接経費: 330千円)
2024年度: 1,170千円 (直接経費: 900千円、間接経費: 270千円)
2023年度: 12,870千円 (直接経費: 9,900千円、間接経費: 2,970千円)
|
キーワード | 波動散乱問題 / 散乱共振 / 高速周波数応答掃引 / メタデバイス / トポロジー最適化 / 境界要素法 |
研究開始時の研究の概要 |
無限領域で定義された波動散乱問題は、振動、音、光などの伝播を記述する数理モデルであり、広い分野における古くからの研究対象である。現在ではその数値解法と関連する最適設計法は成熟し、これらを利用したものづくりが行われている。一方で、解の高階周波数微分はこれまでほとんど注目されておらず、対象の周波数応答を高速に掃引する目的の研究が散発的に行われたにすぎない。本研究では、散乱問題の解の高階周波数微分に関する理論的な研究を行うとともに、その高精度かつ高速な数値計算法を開発する。構築した方法論を、高速周波数応答掃引のみならず、複素固有値問題や時間域波動散乱問題の高速解法、最適設計法へと応用する。
|
研究実績の概要 |
研究計画調書に記した計画に従い、開発済みであった2次元ノイマン問題に対する数値計算プログラムを3次元ノイマン問題へと拡張した。当初の予想通り、3次元問題に対しては階層型行列法に基づく直接解法の性能が十分でなかったため、これを反復法へと切り替えた。ここで現れるブロックテプリッツ型の疎なブロック三角行列に対するカルデロンの前処理を新規に開発し、その性能が十分であることを数値実験を通して確認した。この内容については国内ならびに国際学会で発表予定であり、その後、多重極法の実装を行う予定である。また、2次元問題に対しては階層型行列法も有効であったため、当初令和7年度に実施予定であった周期散乱場の高速周波数掃引法の開発を計画を前倒して実施した。平面波入射に対するフォノニック構造の周波数応答をいくつか計算したところ、提案法が従来法の性能を凌駕するとともに実用上有用であることが確認できたため、内容をまとめ国際誌に投稿した(現在査読中)。また、当初計画していなかったが重要なクラスの問題としてトランスミッション問題が挙げられるが、これに対する実装が容易で性能の高い境界要素法とその前処理の定式化を着想したため、計画を変更して研究を実施した。結果として、通常の境界要素法のみならず、高階周波数微分された境界積分積分方程式の解法としても非常に有効であることが明らかとなった。この成果についても既に内容をまとめ、国際誌に投稿・査読中である。
|
現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
当初予定していた3次元問題に対する散乱場の高階周波数微分のための境界要素法を実装できたため、おおむね順調に進展していると評価できると考える。さらに、当初予定していなかったトランスミッション問題に対する解法も構築できたことは計算科学の基盤研究として意義深いと考えている。一方、3次元問題に対する多重極法の実装は未だ手付かずであることから最上位の区分は選択しなかった。
|
今後の研究の推進方策 |
まず、前年度に開発した3次元問題に対する新しい境界要素法に対し、高速多重極法を実装する。開発した数値計算プログラムを用い、例えば全可聴周波数域における音響散乱解析を実施し、成果を国際会議・国際誌において発表する。その後は研究計画調書に沿って研究を進める。すなわち、散乱共鳴周波数とその感度解析について、ヘルムホルツ方程式の解の周波数微分を利用した新規の手法を開発するとともに、高速時間域境界要素法に関する基礎的な研究を実施する予定である。また、これまでに開発したスカラー波動場に対する方法論を、弾性波動方程式やマクスウェル方程式へと拡張することを試みる。
|