Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach

Fujii, Hiroyuki; Yamada, Yukio; Chiba, Go; Hoshi, Yoko; Kobayashi, Kazumichi; Watanabe, Masao

Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach

Fujii, Hiroyuki; Yamada, Yukio; Chiba, Go et al.

Journal of Computational Physics, 2018, 374, 591-604

Permalink : https://hdl.handle.net/2115/79910

このアイテムのアクセス数:635件（2026-08-18 02:17 集計）

閲覧可能ファイル

ファイル	フォーマット	サイズ	閲覧回数	説明
JCP_Fujii	pdf	442 KB	630

論文情報

ファイル出力 EndNote Basic出力 Mendeley出力

アクセス権	open access
URI	https://hdl.handle.net/2115/79910
タイトル	Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach
著者	Fujii, Hiroyuki ; Yamada, Yukio ; Chiba, Go ; Hoshi, Yoko ; Kobayashi, Kazumichi ; Watanabe, Masao
e-Rad_Researcher 00632580 著者名 Fujii, Hiroyuki
e-Rad_Researcher 10334583 著者名 Yamada, Yukio
e-Rad_Researcher 50421524 著者名 Chiba, Go
e-Rad_Researcher 50332383 著者名 Hoshi, Yoko
e-Rad_Researcher 80453140 著者名 Kobayashi, Kazumichi
e-Rad_Researcher 30274484 著者名 Watanabe, Masao
言語	英語
キーワード	Radiative transfer equation
	Discrete ordinates method
	Highly forward-peaked scattering phase function
	Renormalization approach
	Double-exponential formula
発行日	2018-12-01
出版者	Elsevier
収録物名	Journal of Computational Physics
巻(号)	374
ページ	591-604
抄録	This paper considers an accurate and efficient numerical scheme for solving the radiative transfer equation (RTE) based on the discrete ordinates method (DOM) in highly forward-peaked scattering media. The DOM approximates the scattering integral of the RTE including the phase function into a quadrature sum with the total number of discrete angular directions (TND) by a quadrature set. Due to large numerical errors of the scattering integral based on the DOM in highly forward-peaked scattering, the phase function is renormalized to satisfy its normalization conditions. Although the renormalization approaches of the phase function improve the accuracy of the numerical results of the RTE, the computational efficiency of the RTE is still required. This paper develops the first order renormalization approach using the double exponential formula for three quadrature sets: level symmetric even, even and odd, and Lebedev sets in a wide range of the TND from 48 to 1454. Numerical errors of the three-dimensional time-dependent RTE are investigated by the analytical solutions of the RTE. The investigation shows that the level symmetric even set with the TND of 48 using the developed approach provides the most accurate results of the RTE among the quadrature sets in the range of the TND, while to obtain the same accuracy by the conventional zeroth order renormalization approach, the TND needs to be larger than 360. The results suggest the large reduction of computational loads by the developed approach to less than 10 percent from those in the conventional approach.
資源タイプ	学術雑誌論文
出版タイプ	Accepted Manuscript
権利情報	Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
権利情報	© 2018. This manuscript version is made available under the CC-BY-NC-ND 4.0 license https://creativecommons.org/licenses/by-nc-nd/4.0/
PISSN	0021-9991
NCID	AA00696013
関連情報 (isVersionOf)	出版社版DOI https://doi.org/10.1016/j.jcp.2018.07.047