Abstract:Making the conversion from the geometrical spatial scale to the optical depth spatial scale is useful in obtaining numerical solutions for the radiative transfer equation. This is because it allows for the use of exponential integrators, while enforcing numerical stability. Such a conversion involves the integration of the total opacity of the medium along the considered ray path. This is usually approximated by applying a piecewise quadrature in each spatial cell of the discretized medium. However, a rigorous analysis of this numerical step is lacking. This work is aimed at clearly assessing the performance of different optical depth conversion schemes with respect to the solution of the radiative transfer problem for polarized radiation, out of the local thermodynamic equilibrium. We analyzed different optical depth conversion schemes and their combinations with common formal solvers, both in terms of the rate of convergence as a function of the number of spatial points and the accuracy of the emergent Stokes profiles. The analysis was performed in a 1D semi-empirical model of the solar atmosphere, both in the absence and in the presence of a magnetic field. We solved the transfer problem of polarized radiation in different settings: the continuum, the photospheric Sr i AA modeled under the assumption of complete frequency redistribution, and the chromospheric Ca i AA taking the partial frequency redistribution effects into account during the modeling. High-order conversion schemes generally outperform low-order methods when a sufficiently high number of spatial grid points is considered. In the synthesis of the emergent Stokes profiles, the convergence rate, as a function of the number of spatial points, is impacted by both the conversion scheme and formal solver. The use of low-order conversion schemes significantly reduces the accuracy of high-order formal solvers. In practical applications, the use of low-order optical depth conversion schemes introduces large numerical errors in the formal solution. To fully exploit high-order formal solvers and obtain accurate synthetic emergent Stokes profiles, it is necessary to use high-order optical depth conversion schemes.

Radiant fluence from ray tracing in optical multipass systems

Quasi-Monte Carlo Ray Tracing Algorithm for Radiative Flux Distribution Simulation

Comparison of radiant intensity in aqueous media using experimental and numerical simulation techniques

GPU-based Monte Carlo Ray Tracing Simulation Considering Refraction for Central Receiver System

An Efficient 3D Tube-Level Flux-Thermal Model of External Tubular Receivers in Solar Power Tower Systems

Radiance Cascades: A Novel High-Resolution Formal Solution for Multidimensional Non-LTE Radiative Transfer

Impact of the numerical conversion to optical depth on the transfer of polarized radiation

Transient and Steady-State Combined Heat-Transfer in Semitransparent Materials Subjected to A Pulse or A Step Irradiation

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

In Situ Detailed Fluence Rate Distributions in a UV Reactor with Multiple Low-Pressure Lamps: Comparison of Experimental and Model Results

Inverse Diffusion Approximation for Extraction of Scattering and Absorption Coefficients in Highly Scattering Media

Numerical Characterization of a High Flux Solar Simulator Using Forward and Inverse Methods

Holistic opto-thermal simulation framework for high-brightness light sources based on fluorescent conversion

Optical Design of a High Flux Solar Simulator Configurable to Achieve Either a Line or a Point Focus

A three-dimensional finite elements approach for the coupled radiative transfer equation and diffusion approximation modeling in fluorescence imaging

Opto-thermal simulation model for optimizing laser-excited remote phosphor systems

Development of a generalized radiation calculation system

A finite element approach for modeling photon transport in tissue

Finite Element Simulation for Short Pulse Light Radiative Transfer in Homogeneous and Nonhomogeneous Media

A low-time complexity semi-analytic Monte Carlo radiative transfer model: Application to optical characteristics of complex spatial targets

Characterization of the Optical Properties of Photoluminescent Turbid Media Using an Integrating Sphere and Monte Carlo Simulations