HUSCAP logo Hokkaido Univ. logo

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

Files in This Item:
JCP_Fujii.pdf442.43 kBPDFView/Open
Please use this identifier to cite or link to this item:

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/
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
Type: article (author version)
Appears in Collections:工学院・工学研究院 (Graduate School of Engineering / Faculty of Engineering) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)

Submitter: 藤井 宏之

Export metadata:

OAI-PMH ( junii2 , jpcoar_1.0 )

MathJax is now OFF:


 - Hokkaido University