| 研究課題/領域番号 |
16K05302
|
|---|
| 研究機関 | 苫小牧工業高等専門学校 |
|---|
研究代表者 |
高橋 労太 苫小牧工業高等専門学校, 創造工学科, 准教授 (40513453)
|
|---|
| 研究分担者 |
梅村 雅之 筑波大学, 計算科学研究センター, 教授 (70183754)
|
|---|
| 研究期間 (年度) |
2016-04-01 – 2020-03-31
|
| キーワード | 輻射流体 / 一般相対論 / ブラックホール / ボルツマン方程式 / モンテカルロ計算 |
|---|
| 研究実績の概要 |
本研究の研究代表者と分担者が開発した一般相対論的講師ボルツマン方程式ソルバーであるARTIST法（R. Takahashi, M. Umemura, 2017, MNRAS, 464, 4567）では、湾曲次空中での輻射輸送方程式を直接数値的に解くことを試みた。このソルバーでは光子散乱の効果は形式的は取り入れることが出来ているが、光子の多重散乱の効果が利き、光子拡散の効果が無視できないような状況は正確には取り扱うことができない。特に、相対論的流体中での光子多重散乱の効果がどのようになっているのかは、よく把握されていなかった。今年度は、当初の予定通り、この問題を解決することに取り組んだ。その結果、相対論的流体中での光子多重散乱を記述する全く新しい解析解を発見することができた。
研究は、光子の任意回数の酔歩を記述する解析的計算とGPU並列計算機を用いたモンテカルロ・シミュレーションの結果を比較しながら進めた。数値シミュレーションでは、約10億個の光子を用いて計算を行い、流体速度を表すローレンツ因子は任意の値で実行できる数値コードを用いた。解析的に書かれた式は、任意の散乱回数、任意の流体速度に対する時空中の光子確率密度関数を表し、モンテカルロ計算の結果を完全に再現することが確認できた。また、モンテカルロ計算では、粒子数の限界から精度よく求めることができない領域の確率密度関数も厳密に計算できることがわかった。この解析式から全ての散乱回数の効果を含んだ光子確率密度関数を得ることもできた。昨年度の段階では、数値的に確率密度関数を記述することを行っていたが、今回、解析的な式が存在し、それが発見されたことは、予想外の結果であった。これらの計算結果の一部は、日本天文学会などで公表し、完全な結果は次の日本天文学会で公表する予定である。また、解析式を計算した論文は現在執筆中である。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
光子多重散乱の効果を記述する解析式が発見されたことは、当初の予想を大きく上回る成果と判断できるが、解析的計算にやや時間がかかったので、全体の進行状況が少し遅れている。これら両面を考え、研究はおおむね順調に進展していると判断した。
|
| 今後の研究の推進方策 |
研究の推進方策として、主に２つ考えている。一つは、光子多重散乱の効果を解析的に記述できる計算技術をより一般的もしくは多様な状況に発展させる方向である。この方向では、一般相対論的流体中での散乱、非等方散乱、非弾性散乱などへの拡張が可能であると考えられる。もう一つの方向は、新たに得られた光子多重散乱の効果を記述する解析式を一般相対論的光子ボルツマン方程式ソルバーであるARTIST法と組み合わせることで、任意の光学的深さでの一般相対論的輻射輸送計算を実行可能な数値コードを作成する方向である。現時点では、これらの両方向に研究を進めて行きたいと考えている。
|
| 次年度使用額が生じた理由 |
研究計画で購入を予定していたGPU並列計算機を購入したが、予算残が生じたので次年度の研究では物品費及び旅費の一部として使用することを予定している。
|
