Abstract:Geometric sensitivity analyses of the k-eigenvalue have many applications in analyses of geometric uncertainty, calculations of differential control rod worth, and searches for critical geometry. The adjoint-weighted first-order geometric sensitivity theory is widely used and has continuously evolved with the Monte Carlo methods. However, the existing adjoint-weighted algorithm can do only uniform isotropic expansions or contractions of surfaces. The adjoint-weighted algorithm also requires computation of adjoint-weighted scattering and fission reaction rates exactly at material interfaces, which has an infinitesimal probability in reality. This paper presents an improved geometry adjoint-weighted perturbation algorithm that is incorporated into the continuous-energy Reactor Monte Carlo (RMC) code. The improvement of the adjoint-weighted algorithm is decomposed into three steps for constructing a cross-section function of geometric parameters using logical expressions, calculating the derivative of the cross-section function, and estimating the adjoint-weighted surface reaction rates. The improved algorithm can accommodate common one-parameter geometric perturbations of internal interfaces or boundary surfaces as well as those of cells as long as the perturbed cells can be described by logical expressions of spatial surface equations. The perturbation algorithm is compared with a direct difference method, the linear least-squares fitting method with central differences, for several typical geometric perturbations including translation, fixed-axis rotation, and uniform isotropic/anisotropic expansion transformations of planar, spherical, cylindrical, and conical surfaces. The differences between the two methods are not more than 3% and not more than 3 for the majority of the test examples. Even though the perturbation algorithm has higher figures of merit than the direct difference method for the majority of the test examples, there is no guarantee that the former can always be more efficient than the latter. The limitation in the efficiency of the perturbation algorithm was demonstrated by the totally reflecting light water reactor pin model.

Development of General Next Event Surface Crossing Estimators for Monte Carlo Method

Next event estimation for neutron transport Monte Carlo simulations in moving objects

The method of CADIS-NEE and CADIS-DXTRAN in NECP-MCX and its applications

Linked calculation method based on next event estimation for directional shielding problems

Quasi-Monte Carlo Ray Tracing Algorithm for Radiative Flux Distribution Simulation

Monte-Carlo Methods for the Neutron Transport Equation

Application of the First Collision Source Method to CSNS Target Station Shielding Calculation

Generalized Integral Method for Solving the Neutron Transport Equation with the Hybridized Discontinuous Galerkin Method

Calculating the K-Eigenvalue Sensitivity to Typical Geometric Perturbations with the Adjoint-Weighted Method in the Continuous-Energy Reactor Monte Carlo Code RMC

Calculation of Ex-Core Detector 3D Space Response Function via Monte Carlo Adjoint Transport Method

A Universal Adjoint-Weighted Algorithm for Geometric Sensitivity Analysis of K-Eigenvalue Based on Continuous-Energy Monte Carlo Method

Development of a Radiation Shielding Monte Carlo Code: RShieldMC

Calculation and Characteristic Analysis of Ex-Core Detector Spatial Response Function for High Temperature Gas-Cooled Reactor

Calculation of eigenfunction fluxes in nuclear systems

A new Excel-based Monte Carlo transport code for simulating energy response spectra and efficiencies in gamma-ray detection materials

Verification and Analysis of the Problem-Dependent Multigroup Macroscopic Cross-Sections for Shielding Calculations

A conservative first-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Generalized Sensitivity Analysis with Continuous-Energy Monte Carlo Code RMC

An extension of NASCA for SDDR analysis of mobile activated components with global variance reduction

Convergence Characteristics and Acceleration of the Transient Fixed Source Equation Solved by Monte Carlo Method

How Efficient Are Monte Carlo Calculations Together With the Q-System to Determine Radioactive Transport Limits? Case Study on Medical Radionuclides