2016 Fiscal Year Annual Research Report
Development of Neutrino Radiation Hydrodynamic Schemes toward the Era of exa-scale computing
| Project/Area Number |
26870823
|
|---|
| Research Institution | National Astronomical Observatory of Japan |
|---|
Principal Investigator |
滝脇 知也 国立天文台, 理論研究部, 助教 (50507837)
|
|---|
| Project Period (FY) |
2014-04-01 – 2017-03-31
|
| Keywords | 高精度計算 / ニュートリノ / 輻射輸送 |
|---|
| Outline of Annual Research Achievements |
本研究の主眼は次世代のスパコンの仕様に合ったニュートリノ輻射輸送の手法を開発することであった。昨年度、その主たる開発目標であるM1-Closureの論文が出版されており、本年度はそのコードを速く走らせるためのチューニングに力を注いだ。また、輻射輸送法を他グループと比較する計画も大いに発展した。一方Sn法はコード開発はほぼ終了しているが、サイエンティフィックランと論文執筆は今後行う予定であり急ぐ必要がある。
本研究で用いられているM1-Closureではニュートリノのエネルギーごとにグリッドを切り、それぞれ計算している。特にニュートリノと物質の相互作用を計算するパートでは、その個数と運動量の４変数について、一般に240x240逆行列を解く必要がある。最新のスパコンを用いても1日1msしか計算が進まないような状況であったため、逆行列の解き方等を工夫し、計算を数倍スピードアップさせた。現在、このコードで計算を進めている。
また、近年、”似たような”セットアップで計算しているのに、計算するグループごとで結果が違うということが超新星のモデル化で大きな課題になっている。この問題を解決するため、アメリカ・ドイツのグループの計算に対して慎重にセットアップを合わせた。結果、ほとんど同じ結果を得ることが分かった。つまり”似たような”といっても少し違うセットアップで計算していることが結果の違いを生んでいることがわかった。
Sn計算はこれまでラグランジュ的手法で解かれていたコードを我々のオイラー的手法で解くコードに移植し、テスト計算でおよその一致を得るところまで計算が進んでいる。
|
