Consistent Adjoint Driven Importance Sampling using Space, Energy, and Angle

For challenging radiation transport problems, hybrid methods combine the accuracy of Monte Carlo methods with the global information present in deterministic methods. One of the most successful hybrid methods is CADIS Consistent Adjoint Driven Importance Sampling. This method uses a deterministic adjoint solution to construct a biased source distribution and consistent weight windows to optimize a specific tally in a Monte Carlo calculation. The method has been implemented into transport codes using just the spatial and energy information from the deterministic adjoint and has been used in many applications to compute tallies with much higher figures-of-merit than analog calculations. CADIS also outperforms user-supplied importance values, which usually take long periods of user time to develop. This work extends CADIS to develop weight windows that are a function of the position, energy, and direction of the Monte Carlo particle. Two types of consistent source biasing are presented: one method that biases the source in space and energy while preserving the original directional distribution and one method that biases the source in space, energy, and direction. Seven simple example problems are presented which compare the use of the standard space/energy CADIS with the new space/energy/angle treatments.

[1]  T. L. Becker,et al.  Hybrid Monte Carlo/Deterministic Methods for Radiation Shielding Problems. , 2009 .

[2]  Douglas E. Peplow,et al.  Monte Carlo Shielding Analysis Capabilities with MAVRIC , 2011 .

[3]  Douglas E. Peplow,et al.  Hybrid Monte Carlo/Deterministic Methods for Accelerating Active Interrogation Modeling , 2013 .

[4]  Thomas M. Evans,et al.  Automated Weight-Window Generation for Threat Detection Applications Using ADVANTG , 2009 .

[5]  Thomas M. Evans,et al.  Simultaneous Optimization of Tallies in Difficult Shielding Problems , 2009 .

[6]  R. P. Gardner,et al.  Monte Carlo nuclear well logging benchmark problems with preliminary intercomparison results , 1991 .

[7]  Alireza Haghighat,et al.  MONTE CARLO VARIANCE REDUCTION WITH DETERMINISTIC IMPORTANCE FUNCTIONS , 2003 .

[8]  Alireza Haghighat,et al.  Automated variance reduction of Monte Carlo shielding calculations using the discrete ordinates adjoint function , 1998 .

[9]  Thomas M. Evans,et al.  Automated Variance Reduction Applied to Nuclear Well-Logging Problems , 2008 .

[10]  T. M. Evans,et al.  An enhanced geometry-independent mesh weight window generator for MCNP , 1997 .

[11]  J Sweezy,et al.  Automated variance reduction for MCNP using deterministic methods. , 2005, Radiation protection dosimetry.

[12]  P. D. Soran,et al.  AVATAR -- Automatic variance reduction in Monte Carlo calculations , 1997 .

[13]  Kevin T. Clarno,et al.  Denovo: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE , 2010 .

[14]  N. Nariyama,et al.  Systematic Evaluation of Neutron Shielding Effects for Materials , 1996 .