Variational Variance Reduction for Monte Carlo Eigenvalue and Eigenfunction Problems

A new variational variance reduction (VVR) technique is developed for improving the efficiency of Monte Carlo multigroup nuclear reactor eigenvalue and eigenfunction calculations. The VVR method employs a variational functional, which requires detailed estimates of both the forward and adjoint fluxe...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Nuclear science and engineering Ročník 146; číslo 2; s. 121 - 140
Hlavní autoři: Densmore, Jeffery D., Larsen, Edward W.
Médium: Journal Article
Jazyk:angličtina
Vydáno: La Grange Park, IL Taylor & Francis 01.02.2004
American Nuclear Society
Témata:
ISSN:0029-5639, 1943-748X
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:A new variational variance reduction (VVR) technique is developed for improving the efficiency of Monte Carlo multigroup nuclear reactor eigenvalue and eigenfunction calculations. The VVR method employs a variational functional, which requires detailed estimates of both the forward and adjoint fluxes. The direct functional, employed in standard Monte Carlo calculations, requires only limited information concerning the forward flux. The variational functional requires global information about the forward and adjoint fluxes and hence is more expensive to evaluate but is more accurate than the direct functional. In calculations, this increased accuracy outweighs the extra expense, resulting in a more efficient Monte Carlo simulation. In our work, we evaluate the variational functional using Monte Carlo-calculated forward flux estimates and deterministically calculated adjoint flux estimates. Also, we represent the adjoint flux as a low-order polynomial in space and angle, which is accurate for diffusive systems. (In such systems, which are common in reactor analysis problems, the angular flux is locally nearly linear in space and angle.) Using this adjoint representation, we develop specific VVR methods for eigenvalue problems, in which an estimate of the eigenvalue k in a criticality calculation is desired, and eigenfunction problems, in which an estimate of a detector response due to a fission neutron source during a criticality calculation is desired. The resulting VVR method is very efficient for the problems of interest. With a set of example problems, we demonstrate the increased efficiency of the VVR method over standard Monte Carlo.
Bibliografie:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0029-5639
1943-748X
DOI:10.13182/NSE04-A2398