skip to main content

Title: SECOND-ORDER CROSS TERMS IN MONTE CARLO DIFFERENTIAL OPERATOR PERTURBATION ESTIMATES

Given some initial, unperturbed problem and a desired perturbation, a second-order accurate Taylor series perturbation estimate for a Monte Carlo tally that is a function of two or more perturbed variables can be obtained using an implementation of the differential operator method that ignores cross terms, such as in MCNP4C{trademark}. This requires running a base case defined to be halfway between the perturbed and unperturbed states of all of the perturbed variables and doubling the first-order estimate of the effect of perturbing from the ''midpoint'' base case to the desired perturbed case. The difference between such a midpoint perturbation estimate and the standard perturbation estimate (using the endpoints) is a second-order estimate of the sum of the second-order cross terms of the Taylor series expansion. This technique is demonstrated on an analytic fixed-source problem, a Godiva k{sub eff} eigenvalue problem, and a concrete shielding problem. The effect of ignoring the cross terms in all three problems is significant.
Authors:
;
Publication Date:
OSTI Identifier:
776530
Report Number(s):
LA-UR-01-1640
TRN: AH200123%%119
DOE Contract Number:
W-7405-ENG-36
Resource Type:
Conference
Resource Relation:
Conference: Conference title not supplied, Conference location not supplied, Conference dates not supplied; Other Information: PBD: 1 Mar 2001
Research Org:
Los Alamos National Lab., NM (US)
Sponsoring Org:
US Department of Energy (US)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; CONCRETES; EIGENVALUES; IMPLEMENTATION; SERIES EXPANSION; SHIELDING