Paper
29 September 2009 Fourier decomposition of sharply peaked phase functions: Legendre expansions versus trapezoidal rule
Author Affiliations +
Abstract
Most numerical methods used in radiation transfer problems decompose the radiation field in Fourier series and solve a reduced radiative transfer equation for each Fourier mode. Classically, Legendre polynomials expansions provide the kernels of these reduced transfer equations. For highly peaked phase functions, the Legendre series expansion converges very slowly. We show in this paper that this expansion can advantageously be replaced by a direct numerical evaluation using the trapezoidal rule. The improvement afforded by this direct evaluation yields highly accurate results with orders of magnitude fewer arithmetic operations than the Legendre series, avoids the very slow convergence of the Legendre series and exploits instead the rapid decay of the Fourier coefficients for exponential convergence, and finally bypasses the need for phase function truncations.
© (2009) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Alain Sei "Fourier decomposition of sharply peaked phase functions: Legendre expansions versus trapezoidal rule", Proc. SPIE 7475, Remote Sensing of Clouds and the Atmosphere XIV, 74750J (29 September 2009); https://doi.org/10.1117/12.828842
Lens.org Logo
CITATIONS
Cited by 1 scholarly publication.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Scattering

Radiative transfer

Clouds

Numerical analysis

MATLAB

Particles

Spherical lenses

Back to Top