Exact Analytical Solution of the One-Dimensional Time-Dependent Radiative Transfer Equation with Linear Scattering

Vladimir Allaxwerdian,Dmitry V.Naumov
DOI: https://doi.org/10.1016/j.jqsrt.2023.108726
2024-01-17
Abstract:The radiative transfer equation (RTE) is a cornerstone for describing the propagation of electromagnetic radiation in a medium, with applications spanning atmospheric science, astrophysics, remote sensing, and biomedical optics. Despite its importance, an exact analytical solution to the RTE has remained elusive, necessitating the use of numerical approximations such as Monte Carlo, discrete ordinate, and spherical harmonics methods.
Optics,High Energy Physics - Experiment,Mathematical Physics,Classical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to derive the exact analytical solution of the one - dimensional time - dependent radiative transfer equation (RTE). Specifically, the authors aim to provide an accurate mathematical description for the propagation of light in random media (such as water or ice). The importance of this problem lies in the fact that although the radiative transfer equation has a wide range of applications in describing the propagation of electromagnetic radiation through media, in the one - dimensional time - dependent case, no exact analytical solution has been found previously. ### Problem Background The radiative transfer equation (RTE) is the fundamental equation for describing the propagation of electromagnetic radiation in media and has important applications in fields such as atmospheric science, astrophysics, remote sensing, and biomedical optics. However, due to its complexity, numerical approximation methods (such as the Monte Carlo method, the discrete - ordinate method, and the spherical - harmonic - function method) are usually required to solve it. ### Main Contributions of the Paper 1. **Exact Analytical Solution**: The paper presents the exact analytical solution of the one - dimensional time - dependent RTE and shows the transition to the diffusion region by analyzing the moments of the photon distribution. 2. **Equivalence and Causality**: It is proved that the one - dimensional RTE is equivalent to the Klein - Gordon equation with an imaginary mass term, and under appropriate boundary conditions, the solution is strictly causal rather than superluminal. 3. **Improvement of the Monte Carlo Method**: By adjusting the parameters of Monte Carlo sampling, a method with faster convergence is proposed, which significantly improves the efficiency of the Monte Carlo method. 4. **Verification and Benchmark Testing**: The correctness of the analytical solution is verified by Monte Carlo simulation, and the performance of the Monte Carlo method is benchmark - tested. ### Specific Contents - **Equation Form**: Starting from the one - dimensional time - dependent RTE, the paper finally obtains the exact solution containing the modified Bessel function by introducing new variables and transformations. - **Monte Carlo Method**: By changing the input parameters \((\mu_s, g)\) to \((\mu_s(1 - g)/2, - 1)\), a faster convergence rate is achieved. - **Result Discussion**: The moments of the photon flux are discussed in detail, showing the propagation characteristics of light in the medium and explaining why the propagation of light will transform into diffusion behavior under certain conditions. ### Conclusion This research not only makes important contributions to the theoretical understanding of the RTE but also provides new ideas for improving numerical methods (such as the Monte Carlo method) to solve the RTE more accurately. In addition, the research results show that the Klein - Gordon equation with an imaginary mass term can maintain causality under appropriate boundary conditions, which is contrary to the intuitive expectation, emphasizing the importance of careful analysis of light - transmission problems.