Hokkaido University Collection of Scholarly and Academic Papers >
Graduate School of Engineering / Faculty of Engineering >
Peer-reviewed Journal Articles, etc >
Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach
This item is licensed under:Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
Title: | Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach |
Authors: | Fujii, Hiroyuki Browse this author →KAKEN DB | Yamada, Yukio Browse this author →KAKEN DB | Chiba, Go Browse this author →KAKEN DB | Hoshi, Yoko Browse this author →KAKEN DB | Kobayashi, Kazumichi Browse this author →KAKEN DB | Watanabe, Masao Browse this author →KAKEN DB |
Keywords: | Radiative transfer equation | Discrete ordinates method | Highly forward-peaked scattering phase function | Renormalization approach | Double-exponential formula |
Issue Date: | 1-Dec-2018 |
Publisher: | Elsevier |
Journal Title: | Journal of Computational Physics |
Volume: | 374 |
Start Page: | 591 |
End Page: | 604 |
Publisher DOI: | 10.1016/j.jcp.2018.07.047 |
Abstract: | 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. |
Rights: | © 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/ | https://creativecommons.org/licenses/by-nc-nd/4.0/ |
Type: | article (author version) |
URI: | http://hdl.handle.net/2115/79910 |
Appears in Collections: | 工学院・工学研究院 (Graduate School of Engineering / Faculty of Engineering) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)
|
Submitter: 藤井 宏之
|