A boundary-based net-exchange Monte Carlo method for absorbing and scattering thick media

A boundary-based net-exchange Monte Carlo method was introduced (JQSRT 74 (2001) 563) that allows to bypass the difficulties encountered by standard Monte Carlo algorithms in the limit of optically thick absorption (and/or for quasi-isothermal configurations). With the present paper, this method is extended to scattering media. Developments are fully 3D, but illustrations are presented for plane parallel configuration. Compared to standard Monte Carlo algorithms, convergence qualities have been improved over a wide range of absorption and scattering optical thicknesses. The proposed algorithm still encounters a convergence difficulty in the case of optically thick, highly scattering media.