Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Efficient Approach for Optical Radiative Transfer Tomography using the Spherical Harmonics Discrete Ordinates Method
94
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper introduces a method to preform optical tomography, using 3D radiative transfer as the forward model. We use an iterative approach predicated on the Spherical Harmonics Discrete Ordinates Method (SHDOM) to solve the optimization problem in a scalable manner. We illustrate with an application in remote sensing of a cloudy atmosphere.
rate research
Read More
We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with a ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed frame approach; the source function is evaluated in the comoving frame whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated with various test problems and comparisons with results of a relativistic Monte Carlo code. These validations confirm that the code correctly calculates intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.
Spherical Harmonics, $Y_ell^m(theta,phi)$, are derived and presented (in a Table) for half-odd-integer values of $ell$ and $m$. These functions are eigenfunctions of $L^2$ and $L_z$ written as differential operators in the spherical-polar angles, $theta$ and $phi$. The Fermion Spherical Harmonics are a new, scalar and angular-coordinate-dependent representation of fermion spin angular momentum. They have $4pi$ symmetry in the angle $phi$, and hence are not single-valued functions on the Euclidean unit sphere; they are double-valued functions on the sphere, or alternatively are interpreted as having a double-sphere as their domain.
The radiative transfer equation models the interaction of radiation with scattering and absorbing media and has important applications in various fields in science and engineering. It is an integro-differential equation involving time, space and angular variables and contains an integral term in angular directions while being hyperbolic in space. The challenges for its numerical solution include the needs to handle with its high dimensionality, the presence of the integral term, and the development of discontinuities and sharp layers in its solution along spatial directions. Its numerical solution is studied in this paper using an adaptive moving mesh discontinuous Galerkin method for spatial discretization together with the discrete ordinate method for angular discretization. The former employs a dynamic mesh adaptation strategy based on moving mesh partial differential equations to improve computational accuracy and efficiency. Its mesh adaptation ability, accuracy, and efficiency are demonstrated in a selection of one- and two-dimensional numerical examples.
Many applications in science call for the numerical simulation of systems on manifolds with spherical topology. Through use of integer spin weighted spherical harmonics we present a method which allows for the implementation of arbitrary tensorial evolution equations. Our method combines two numerical techniques that were originally developed with different applications in mind. The first is Huffenberger and Wandelts spectral decomposition algorithm to perform the mapping from physical to spectral space. The second is the application of Luscombe and Lubans method, to convert numerically divergent linear recursions into stable nonlinear recursions, to the calculation of reduced Wigner d-functions. We give a detailed discussion of the theory and numerical implementation of our algorithm. The properties of our method are investigated by solving the scalar and vectorial advection equation on the sphere, as well as the 2+1 Maxwell equations on a deformed sphere.
The derivation of spherical harmonics is the same in nearly every quantum mechanics textbook and classroom. It is found to be difficult to follow, hard to understand, and challenging to reproduce by most students. In this work, we show how one can determine spherical harmonics in a more natural way based on operators and a powerful identity called the exponential disentangling operator identity (known in quantum optics, but little used elsewhere). This new strategy follows naturally after one has introduced Dirac notation, computed the angular momentum algebra, and determined the action of the angular momentum raising and lowering operators on the simultaneous angular momentum eigenstates (under $hat L^2$ and $hat L_z$).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
