Poisson-Box Sampling algorithms for three-dimensional Markov binary mixtures. (February 2018)

• Record Type: Journal Article
• Authors: Larmier, Coline; Zoia, Andrea; Malvagi, Fausto; Dumonteil, Eric; Mazzolo, Alain
• Is Part Of: Journal of quantitative spectroscopy & radiative transfer. Volume 206(2018)
• Page Start: 70
• Page End: 82
• Publication Date: 2018-02
• Subjects: Chord Length Sampling -- Markov geometries -- Poisson -- Box -- Benchmark -- Monte Carlo -- Tripoli-4®; Spectrum analysis -- Periodicals; Radiation -- Periodicals; Analyse spectrale -- Périodiques; Rayonnement -- Périodiques; Radiation; Spectrum analysis; Periodicals; 543.0858
• DOI: 10.1016/j.jqsrt.2017.10.020
• Languages: English
• ISSNs: 0022-4073

Abstract

Highlights: We introduce a family of new algorithms to generalize the Chord Length Sampling method. We test these Poisson-Box Sampling (PBS) methods against CLS and reference solutions. We revisit the Adams, Larsen and Pomraning benchmark. PBS is further extended so to take into account partial memory effects. Abstract: Particle transport in Markov mixtures can be addressed by the so-called Chord Length Sampling (CLS) methods, a family of Monte Carlo algorithms taking into account the effects of stochastic media on particle propagation by generating on-the-fly the material interfaces crossed by the random walkers during their trajectories. Such methods enable a significant reduction of computational resources as opposed to reference solutions obtained by solving the Boltzmann equation for a large number of realizations of random media. CLS solutions, which neglect correlations induced by the spatial disorder, are faster albeit approximate, and might thus show discrepancies with respect to reference solutions. In this work we propose a new family of algorithms (called 'Poisson Box Sampling', PBS) aimed at improving the accuracy of the CLS approach for transport in d -dimensional binary Markov mixtures. In order to probe the features of PBS methods, we will focus on three-dimensional Markov media and revisit the benchmark problem originally proposed by Adams, Larsen and Pomraning [1] and extended by Brantley [2]: for these configurations we will compare reference solutions, standard CLS solutions and the new PBS solutions for scalar particle flux, transmission and reflection coefficients. PBS will be shown to perform better than CLS at the expense of a reasonable increase in computational time.