44,364 research outputs found
Hybrid-dual-fourier tomographic algorithm for a fast three-dimensionial optical image reconstruction in turbid media
A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model
Evaluation and development of new VAS remote sensing algorithms
An algorithm, developed to permit the simultaneous retrieval of temperature/moisture profiles and surface skin temperature by direct analytical solution of the radiative transfer equation, alleviates the problem associated with the interdependencies of water vapor retrieval on temperature retrieval and their dependencies on surface emissions. Simultaneous solution of all these quantities is achieved in one calculation using the available radiance observations. Since only a single matrix inversion is required for the specification of all parameters, the solution is computationally efficient. Ancillary observations of temperature and/or moisture from surface sensors on aircraft can be readily incorporated into the solution. Simulation tests of the method indicate improved performance over the previous iterative technique, particularly for the lower troposphere and for water vapor
Radiative transfer in disc galaxies I - A comparison of four methods to solve the transfer equation in plane-parallel geometry
Accurate photometric and kinematic modelling of disc galaxies requires the
inclusion of radiative transfer models. Due to the complexity of the radiative
transfer equation (RTE), sophisticated techniques are required. Various
techniques have been employed for the attenuation in disc galaxies, but a
quantitative comparison of them is difficult, because of the differing
assumptions, approximations and accuracy requirements which are adopted in the
literature. In this paper, we present an unbiased comparison of four methods to
solve the RTE, in terms of accuracy, efficiency and flexibility. We apply them
all on one problem that can serve as a first approximation of large portions of
disc galaxies: a one-dimensional plane-parallel geometry, with both absorption
and multiple scattering taken into account, with an arbitrary vertical
distributions of stars and dust and an arbitrary angular redistribution of the
scattering. We find that the spherical harmonics method is by far the most
efficient way to solve the RTE, whereas both Monte Carlo simulations and the
iteration method, which are straightforward to extend to more complex
geometries, have a cost which is about 170 times larger.Comment: 12 pages, 4 figures, accepted for publication in MNRA
A fast GPU Monte Carlo Radiative Heat Transfer Implementation for Coupling with Direct Numerical Simulation
We implemented a fast Reciprocal Monte Carlo algorithm, to accurately solve
radiative heat transfer in turbulent flows of non-grey participating media that
can be coupled to fully resolved turbulent flows, namely to Direct Numerical
Simulation (DNS). The spectrally varying absorption coefficient is treated in a
narrow-band fashion with a correlated-k distribution. The implementation is
verified with analytical solutions and validated with results from literature
and line-by-line Monte Carlo computations. The method is implemented on GPU
with a thorough attention to memory transfer and computational efficiency. The
bottlenecks that dominate the computational expenses are addressed and several
techniques are proposed to optimize the GPU execution. By implementing the
proposed algorithmic accelerations, a speed-up of up to 3 orders of magnitude
can be achieved, while maintaining the same accuracy
Unified Gas-kinetic Wave-Particle Methods III: Multiscale Photon Transport
In this paper, we extend the unified gas-kinetic wave-particle (UGKWP) method
to the multiscale photon transport. In this method, the photon free streaming
and scattering processes are treated in an un-splitting way. The duality
descriptions, namely the simulation particle and distribution function, are
utilized to describe the photon. By accurately recovering the governing
equations of the unified gas-kinetic scheme (UGKS), the UGKWP preserves the
multiscale dynamics of photon transport from optically thin to optically thick
regime. In the optically thin regime, the UGKWP becomes a Monte Carlo type
particle tracking method, while in the optically thick regime, the UGKWP
becomes a diffusion equation solver. The local photon dynamics of the UGKWP, as
well as the proportion of wave-described and particle-described photons are
automatically adapted according to the numerical resolution and transport
regime. Compared to the -type UGKS, the UGKWP requires less memory cost
and does not suffer ray effect. Compared to the implicit Monte Carlo (IMC)
method, the statistical noise of UGKWP is greatly reduced and computational
efficiency is significantly improved in the optically thick regime. Several
numerical examples covering all transport regimes from the optically thin to
optically thick are computed to validate the accuracy and efficiency of the
UGKWP method. In comparison to the -type UGKS and IMC method, the UGKWP
method may have several-order-of-magnitude reduction in computational cost and
memory requirement in solving some multsicale transport problems.Comment: 27 pages, 15 figures. arXiv admin note: text overlap with
arXiv:1810.0598
- …