When lifting the assumption of spatially-independent scattering centers in classical linear transport theory, collision rate is no longer proportional to vector flux radiance because the macroscopic cross-section @math depends on the distance @math to the previous collision or boundary. This creates a nonlocal relationship between collision rate and flux and requires revising a number of familiar deterministic and Monte Carlo methods. We generalize collision and track-length estimators to support unbiased estimation of either flux integrals or collision rates in generalized radiative transfer (GRT). To provide benchmark solutions for the Monte Carlo estimators, we derive the four Green's functions for the isotropic point source in infinite media with isotropic scattering. Additionally, new moment-preserving diffusion approximations for these Green's functions are derived and reduce to algebraic expressions involving the first four moments of the free-path lengths between collisions.
The Green's function derivations and related moment-preserving diffusion approximations in this paper expand upon our previous work @cite_137 by including the case of uncorrelated emission, and by performing a completely general moment analysis in the Fourier domain that avoids the previous need for Fourier inversions. In addition to proving conjectures from this previous work, we also prove new results about the diffusion coefficient's dependence on spatial dimension @math . The relationship to other diffusion asymptotics derived from a parameter-of-smallness approach @cite_104 @cite_39 is also discussed.
@cite_137 @cite_39 @cite_104
2049719606 2085186189 1967564273 2116625471
The Green’s function of the time dependent radiative transfer equation for the semi-infinite medium is derived for the first time by a heuristic approach based on the extrapolated boundary condition and on an almost exact solution for the infinite medium. Monte Carlo simulations performed both in the simple case of isotropic scattering and of an isotropic point-like source, and in the more realistic case of anisotropic scattering and pencil beam source, are used to validate the heuristic Green’s function. Except for the very early times, the proposed solution has an excellent accuracy (>98 for the isotropic case, and >97 for the anisotropic case) significantly better than the diffusion equation. The use of this solution could be extremely useful in the biomedical optics field where it can be directly employed in conditions where the use of the diffusion equation is limited, e.g. small volume samples, high absorption and or low scattering media, short source-receiver distances and early times. Also it represents a first step to derive tools for other geometries (e.g. slab and slab with inhomogeneities inside) of practical interest for noninvasive spectroscopy and diffuse optical imaging. Moreover the proposed solution can be useful to several research fields where the study of a transport process is fundamental. We consider a random, uniformly elliptic coefficient field (a(x) ) on the (d )-dimensional integer lattice ( Z ^d ). We are interested in the spatial decay of the quenched elliptic Green function (G(a;x,y) ). Next to stationarity, we assume that the spatial correlation of the coefficient field decays sufficiently fast to the effect that a logarithmic Sobolev inequality holds for the ensemble ( ). We prove that all stochastic moments of the first and second mixed derivatives of the Green function, that is, ( | _x G(x,y)|^p ) and ( | _x _y G(x,y)|^p ), have the same decay rates in (|x-y| 1 ) as for the constant coefficient Green function, respectively. This result relies on and substantially extends the one by Delmotte and Deuschel (Probab Theory Relat Fields 133:358–390, 2005), which optimally controls second moments for the first derivatives and first moments of the second mixed derivatives of (G ), that is, ( | _x G(x,y)|^2 ) and ( | _x _y G(x,y)| ). As an application, we are able to obtain optimal estimates on the random part of the homogenization error even for large ellipticity contrast. ABSTRACTMultiarrival Green’s functions are essential in seismic modeling, migration, and inversion. Huygens-Kirchhoff (HK) integrals provide a bridge to integrate locally valid first-arrival Green’s functions into a globally valid multiarrival Green’s function. We have designed robust and accurate finite-difference methods to compute first-arrival traveltimes and amplitudes, so that first-arrival Green’s functions can be constructed rapidly. We adapted a fast butterfly algorithm to evaluate discretized HK integrals. The resulting fast Huygens’ sweeping method has the following unique features: (1) it precomputes a set of local traveltime and amplitude tables, (2) it automatically takes care of caustics, (3) it constructs Green’s functions of the Helmholtz equation for arbitrary frequencies and for many point sources, and (4) for a fixed number of points per wavelength, it constructs each Green’s function in nearly optimal complexity O(N log N) in terms of the total number of mesh points N, where the prefa... In this study, Green’s function of the two-dimensional radiative transfer equation is derived for an infinitely extended anisotropically scattering medium, which is illuminated by a unidirectional source distribution. In the steady-state domain, the final results, which are based on eigenvalues and eigenvectors, are given analytically apart from the eigenvalues. For the time-dependent case an additional numerical inverse Fourier transform is required. The obtained solutions were successfully validated with another exact analytical solution in the time domain for isotropically scattering and with the Monte Carlo method for anisotropically scattering media.