The statistical estimator concept, created in the nuclear engineering field, has been adapted to the elaboration of a new and fast semianalytical Monte Carlo numerical simulation for time-resolved light-scattering problems. This concept has also been generalized to the case of unmatched boundaries. The model, discussed in detail in this paper, contains two stages. The first stage is the information generator in which, for each scattering event, the contribution to the total reflectance and transmittance is evaluated and subtracted from the photon current energy. This procedure reduces the number of photons required to produce a given accuracy, which makes it possible to store all event positions and energies. In the second stage, called the information processor, the results of the first stage are used to calculate analytically any desired result. Examples are given for scattering slabs of isotropic or anisotropic scatterers when collimated-beam incidence is used. Reflections at the boundaries are taken into account. The results obtained either with this new method or with classical Monte Carlo methods are very similar. However, the convergence of our new model is much better and, because of the separation into two stages, any quantity related to the problem can be easily calculated afterward without recomputing the simulation.
© 1996 Optical Society of America
Original Manuscript: July 26, 1995
Revised Manuscript: April 1, 1996
Manuscript Accepted: February 28, 1996
Published: September 1, 1996
Eric Tinet, Sigrid Avrillier, and Jean Michel Tualle, "Fast semianalytical Monte Carlo simulation for time-resolved light propagation in turbid media," J. Opt. Soc. Am. A 13, 1903-1915 (1996)