Increasing the efficiency of radiation shielding calculations by using antithetic variates

John H. Halton and Pradip K. Sarkar

Mathematics and Computers in Simulation (MATCOM), 1998, vol. 47, issue 2, 309-318

Abstract: This paper considers the passage of a parallel, uniform stream of particles through a parallel-sided, homogeneous slab of material capable of scattering and absorption. Particle histories are followed until the particle escapes from the near or far surface of the slab. Each particle carries a weight (initially 1, at the near surface of the slab) and the total weight of particles escaping through the far surface is compared with the total number of particles initiated originally. To increase efficiency of the computation, (i) we disregard absorption, by multiplying the weight of each particle by its non-absorption probability at each collision, and (ii) we select N parallel geometric surfaces inside the slab (parallel to the sides of the slab) and, as a particle crosses any such surface – if it is crossing it in a backward direction – we play “Russian roulette” with probability 1/m (multiplying its weight by m if it survives), while – if it is crossing the surface in a forward direction – we “split” the particle into m identical particles (all initially moving with the same velocity vector), each with weight 1/m times that of the particle we split. In the case of splitting, the paper compares two alternative techniques. Either (1) we select m independent random numbers with which to compute the m free path lengths to the next (scattering) collision, or (2) we select one random number and form m “antithetically” balanced values with which to compute the m free path lengths (with appropriate weight-adjustments). The results obtained are very promising. The use of “antithetic splitting” increases the efficiency (in terms of the work required to obtain a given variance in the resulting estimate of the transmission coefficient) by a factor ranging from about 1 to about 30 (the improvement appears to be most marked for scattering probabilities around 0.5, rather than near 1.0).

Keywords: Radiation; Shielding; Monte Carlo; Antithetic variates (search for similar items in EconPapers)
Date: 1998
References: View complete reference list from CitEc
Citations: View citations in EconPapers (1)

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378475498001074
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:matcom:v:47:y:1998:i:2:p:309-318

DOI: 10.1016/S0378-4754(98)00107-4

Access Statistics for this article

Mathematics and Computers in Simulation (MATCOM) is currently edited by Robert Beauwens

More articles in Mathematics and Computers in Simulation (MATCOM) from Elsevier
Bibliographic data for series maintained by Catherine Liu ().