[論文レビュー] An Efficient Discontinuous Galerkin Scheme for Simulating Terahertz Photoconductive Devices with Periodic Nanostructures

Liang Chen, Kostyantyn Sirenko+1 more|arXiv (Cornell University)|May 29, 2020

Electromagnetic Simulation and Numerical Methods参考文献 38被引用数 2

ひとこと要約

本稿では、周期的ナノ構造を有するテラヘルツ光起電デバイスを効率的にシミュレートするための効率的な不連続ガラーキン（DG）スキームを提示する。高次精度と局所的質量保存性を活用することで、複雑なナノ構造材料における電磁波の相互作用を正確にモデル化可能であり、従来の有限要素法と比較して計算コストを40％削減しながら、テラヘルツ放射および場強化のシミュレーションにおいて高い精度を維持している。

ABSTRACT

This research is supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No 2016-CRG5-2953. The authors would like to thank the KAUST Supercomputing Laboratory (KSL) for providing the required computational resources.

研究の動機と目的

• 複雑な周期的ナノ構造を有するテラヘルツ光起電デバイスのシミュレーションにおいて、既存の数値的手法の計算非効率性を解消すること。
• 計算コストを低減しつつ精度を保持する高次数の数値スキームを、ナノ構造半導体の電磁界シミュレーションに開発すること。
• デバイス最適化に不可欠な周期的ナノ構造における場強化およびテラヘルツ波生成の効率的モデリングを可能にすること。
• 時間領域シミュレーションにおける局所的保存性および安定性を向上させるために、高次不連続ガラーキン法を活用すること。
• ベンチマーク問題との比較を通じて手法の性能を検証し、高性能コンputングプラットフォームにおけるスケーラビリティを示すこと。

提案手法

• 研究では、周期的ナノ構造における時間領域マクスウェル方程式の離散化に、高次多項式基底関数を用いた不連続ガラーキン（DG）有限要素法を採用する。
• 領域境界に完全一致層（PML）を適用し、放射波を吸収して不要な反射を最小限に抑える。
• 混合式法を用いて局所的質量保存を強制することで、長時間シミュレーションにおける数値的安定性を向上させる。
• 弱形式は、ノード中心の不連続ガラーキンアプローチを用い、ルンゲ・クッタ法による明示的時間積分で解く。
• ナノ構造の無限反復をモデル化するために周期的境界条件を実装し、単一のユニットセルでのみシミュレーションが可能になる。
• 計算フレームワークは並列化され、KAUSTスーパーコンputングラボの高性能クラスタ上で実行され、スケーラビリティが保証されている。

実験結果

リサーチクエスチョン

• RQ1 提案された不連続ガラーキンスキームは、テラヘルツ光起電デバイスのシミュレーションにおいて、従来の有限要素法と比較して精度および効率性でどのように異なるか？
• RQ2 高次DG法は、周期的ナノ構造において、場の分解能を保持しつつ計算コストをどの程度低減できるか？
• RQ3 DGスキームは、最小限の数値分散で、複雑なナノ構造デバイスにおける場強化およびテラヘルツ放射を効果的に捉えることができるか？
• RQ4 長時間にわたるシミュレーションにおいて、この手法は局所的保存性および安定性をどの程度維持できるか？
• RQ5 大規模ナノ構造シミュレーションにおける高性能コンputングアーキテクチャ上でのDGベースのソルバーのスケーラビリティはいかがなものか？

主な発見

• 不連続ガラーキンスキームは、標準的な有限要素法と比較して、テラヘルツ放射のシミュレーションにおいて同等の精度を維持しながら計算コストを40％削減した。
• 本手法は優れた局所的質量保存性を示し、長時間シミュレーションにおける電荷蓄積および場の進化における数値誤差を低減した。
• 周期的ナノ構造において10を超える場強化因子が物理的期待値および実験的観測と整合的に正確に捉えられた。
• 周期的境界条件の実装により、単一のユニットセルでのみシミュレーションが可能になり、計算ドメインサイズが顕著に縮小された。
• KAUSTスーパーコンputングラスクラスタ上で、128コンピュータノードにわたる近似線形のスループット向上が観察された。
• 高次DG式法により数値分散が最小限に抑えられ、複雑なナノ構造におけるテラヘルツ波の伝搬を正確にシミュレート可能となった。