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.

Superhistory-based Differential Operator Method for Generalized Responses Sensitivity Calculations

GENERALIZED SENSITIVITY ANALYSIS CAPABILITY WITH THE DIFFERENTIAL OPERATOR METHOD IN RMC CODE

The Hydrodynamics of Low-Density Particles and Optimization of the EMMS Calculation Process

Improved Generalized Perturbation Theory Method for Sensitivity Analysis of Generalized Response Function

Nuclear Data Sensitivity and Uncertainty Analysis for Generalized Response with RMC Code

Improvement of Sensitivity and Uncertainty Analysis Capabilities of Generalized Response in Monte Carlo Code RMC

Computing Eigenvalue Sensitivity Coefficients to Nuclear Data by Adjoint Superhistory Method and Adjoint Wielandt Method Implemented in RMC Code

Asymptotic Wielandt Method and Superhistory Method for Source Convergence in Monte Carlo Criticality Calculation

Preliminary Study of Generalized Sensitivity Analysis with Continuous-Energy Monte Carlo Code RMC

Development of Sensitivity Analysis Capabilities of Generalized Responses to Nuclear Data in Monte Carlo Code RMC

Simultaneous Approach for Simulation of A High-Temperature Gas-Cooled Reactor

Generalized Sensitivity Analysis with Continuous-Energy Monte Carlo Code RMC

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

New Strategies of Sensitivity Analysis Capabilities in Continuous-Energy Monte Carlo Code RMC

Geometric Sensitivity Analysis of Generalized Response Function with RMC Code

The Super Equivalence Method in Monte Carlo Based Homogenization

Sensitivity Embedded Homotopy-Based Backtracking Method for Chemical Process Simulation

Optimal Batch Size Growth for Wielandt Method and Superhistory Method

Generalized Sensitivity Analysis of Responses to Nuclear Data in 2-D/1-D Transport Code KYADJ

Temperature Perturbation Method Using On-the-fly Treatment of the Cross-Sections in the Resolved Resonance Region

Development of an Improved Direct Kinetic Simulation Capability in RMC Code