Monte Carlo methods for radiation transport analysis on vector computers

Abstract

The development of advanced computers with special capabilities for vectorized or parallel calculations demands the development of new calculational methods. The very nature of the Monte Carlo process precludes direct conversion of old (scalar) codes to the new machines. Instead, major changes in global algorithms and careful selection of compatible physics treatments are required. Recent results for Monte Carlo in multigroup shielding applications and in continuous-energy reactor lattice analysis have demonstrated that Monte Carlo methods can be successfully vectorized. The significant effort required for stylized coding and major algorithmic changes is worthwhile, and significant gains in computational efficiency are realized. Speedups of at least twenty to forty times faster than CDC-7600 scalar calculations have been achieved on the CYBER-205 without sacrificing the accuracy of standard Monte Carlo methods. Speedups of this magnitude provide reductions in statistical uncertainties for a given amount of computing time, permit more detailed and realistic problems to be analyzed, and make the Monte Carlo method more accessible to nuclear analysts. Following overviews of the Monte Carlo method for particle transport analysis and of vector computer hardware and software characteristics, both general and specific aspects of the vectorization of Monte Carlo are discussed. Finally, numerical results obtained from vectorized Monte Carlo codes run on the CYBER-205 are presented.

Using Monte Carlo Method For Simulation Of Electromagnetic Radiation Penetration In Bodys Tissues

“…Using random number generation method in recent years has found a wide application in physics [4]. Different simulation software's has been designed and used.…”

Recent Advances and Future Prospects for Monte Carlo

“…The capabilities of these GPGPU processors have grown rapidly in the last few years, leading to significant gains in speed if a calculation can be cast in vector form. Since vectorization of Monte Carlo was proven 30 years ago, 4,5) it would appear straightforward. In practice, however, existing mature codes would have to be entirely restructured and largely rewritten to take advantage of vectorization; that is prohibitively expensive and would invalidate many years of verification and validation work.…”

“…), often involving drums, extended core memory devices, disks, and punched cards. In the 1980s, a few codes were developed for vector computers 4,5,6) and then extended to parallel/vector computers such as the Cray-XMP, -YMP, -C90, -J90 series. During the 1990s, parallel Monte Carlo became commonplace due to the arrival of workstation clusters and modestly priced parallel systems.…”

Parallel Fission Bank Algorithms in Monte Carlo Criticality Calculations

“…Also recognized was the fact that while providing an elegant means run on the first vector computer, the ILLIAC-IV [2]. The general principles from that work were later refined and extended greatly through the work of Forrest Brown in the 1980s [3]. However, as Brown's work shows, the single-instruction multiple-data (SIMD) parallel model inherent to vector processing does not lend itself to the parallelism on particles in Monte Carlo simulations.…”