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

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Allaxwerdian, Vladimir, Naumov, Dmitry V.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866929213192273920
author Allaxwerdian, Vladimir
Naumov, Dmitry V.
author_facet Allaxwerdian, Vladimir
Naumov, Dmitry V.
contents 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. In this paper, we present an exact solution to the one-dimensional time-dependent RTE. We delve into the moments of the photon distribution, providing a clear view of the transition to the diffusion regime. This analysis offers a deeper understanding of light propagation in the medium. Furthermore, we demonstrate that the one-dimensional RTE is equivalent to the Klein-Gordon equation with an imaginary mass term determined by the inverse reduced scattering length. Contrary to naive expectations of superluminal solutions, we find that our solution is strictly causal under appropriate boundary conditions, determined by the light transport problem. We validate the found solution using Monte Carlo simulations and benchmark the performance of the latter. Our analysis reveals that even for highly forward scattering, dozens of random light scatterings are required for an accurate estimate, underscoring the complexity of the problem. Moreover, we propose a method for faster convergence by adjusting the parameters of Monte Carlo sampling. We show that a Monte Carlo method sampling photon scatterings with input parameters $(μ_s,g)$, where $μ_s$ is the inverse scattering length and $g$ is the scattering anisotropy parameter, is equivalent to that with $(μ_s(1-g)/2,-1)$. This equivalence leads to a significantly faster convergence to the exact solution, offering a substantial improvement of the Monte Carlo method for the one-dimensional RTE.
format Preprint
id arxiv_https___arxiv_org_abs_2401_09511
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Exact Analytical Solution of the One-Dimensional Time-Dependent Radiative Transfer Equation with Linear Scattering
Allaxwerdian, Vladimir
Naumov, Dmitry V.
Optics
High Energy Physics - Experiment
Mathematical Physics
Classical Physics
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. In this paper, we present an exact solution to the one-dimensional time-dependent RTE. We delve into the moments of the photon distribution, providing a clear view of the transition to the diffusion regime. This analysis offers a deeper understanding of light propagation in the medium. Furthermore, we demonstrate that the one-dimensional RTE is equivalent to the Klein-Gordon equation with an imaginary mass term determined by the inverse reduced scattering length. Contrary to naive expectations of superluminal solutions, we find that our solution is strictly causal under appropriate boundary conditions, determined by the light transport problem. We validate the found solution using Monte Carlo simulations and benchmark the performance of the latter. Our analysis reveals that even for highly forward scattering, dozens of random light scatterings are required for an accurate estimate, underscoring the complexity of the problem. Moreover, we propose a method for faster convergence by adjusting the parameters of Monte Carlo sampling. We show that a Monte Carlo method sampling photon scatterings with input parameters $(μ_s,g)$, where $μ_s$ is the inverse scattering length and $g$ is the scattering anisotropy parameter, is equivalent to that with $(μ_s(1-g)/2,-1)$. This equivalence leads to a significantly faster convergence to the exact solution, offering a substantial improvement of the Monte Carlo method for the one-dimensional RTE.
title Exact Analytical Solution of the One-Dimensional Time-Dependent Radiative Transfer Equation with Linear Scattering
topic Optics
High Energy Physics - Experiment
Mathematical Physics
Classical Physics
url https://arxiv.org/abs/2401.09511