Monte Carlo solution convergence analysis via COMET in a stylized integral inherently safe LWR benchmark problem
Autor: | Dingkang Zhang, Farzad Rahnema |
---|---|
Rok vydání: | 2017 |
Předmět: |
Mathematical optimization
020209 energy Monte Carlo method 02 engineering and technology Solver 01 natural sciences Standard deviation 010305 fluids & plasmas Hybrid Monte Carlo Nuclear Energy and Engineering 0103 physical sciences Convergence (routing) 0202 electrical engineering electronic engineering information engineering Dynamic Monte Carlo method Benchmark (computing) Applied mathematics Monte Carlo integration Mathematics |
Zdroj: | Annals of Nuclear Energy. 100:3-11 |
ISSN: | 0306-4549 |
DOI: | 10.1016/j.anucene.2016.05.009 |
Popis: | In this paper, the coarse mesh radiation transport (COMET) method is used to study the convergence of the Monte Carlo solution in a stylized benchmark problem based on the beginning of the equilibrium cycle of the Integral Inherently Safe Light Water Reactor concept (I 2 S-LWR) with silicide U 3 Si 2 fuel. The benchmark problem is comprised of 121 fuel assemblies with radial, top and bottom reflectors. COMET consists of a stochastic local solver (response function generator) for generating a library of incident flux response expansion coefficients for the unique coarse meshes (e.g., all fuel assembly types and reflector/non-fuel blocks) and a deterministic iterative solver for whole-core calculations using this response library. The Monte Carlo method is regarded as the gold standard for benchmarking other numerical transport methods. However, for eigenvalue problems, Monte Carlo results require careful analysis to ensure the convergence of the fission source distribution. It is known that this issue can be of concern, particularly for large eigenvalue problems typical of realistic (i.e., operating) reactors. This became evident when the radially-integrated axial fission density distributions predicted by 50 independent MCNP runs in the I 2 S-LWR benchmark problem were compared to each other and were found to be significantly different, with the difference ranging from −9% to +9%. However, the average of the fission density distributions from these 50 runs agreed very well with that predicted by COMET, with a relative difference varying from −0.18% to 0.13%, being within one standard deviation of the average. This agreement confirms that the averaged Monte Carlo solution is converged. |
Databáze: | OpenAIRE |
Externí odkaz: |