Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
An efficient parallel computing scheme for Monte Carlo criticality calculations
KTH, School of Engineering Sciences (SCI), Physics, Reactor Physics.ORCID iD: 0000-0002-7943-7517
KTH, School of Engineering Sciences (SCI), Physics, Reactor Physics.
2009 (English)In: Annals of Nuclear Energy, ISSN 0306-4549, E-ISSN 1873-2100, Vol. 36, no 8, 1276-1279 p.Article in journal (Refereed) Published
Abstract [en]

The existing parallel computing schemes for Monte Carlo criticality calculations suffer from a low efficiency when applied on many processors. We suggest a new fission matrix based scheme for efficient parallel computing. The results are derived from the fission matrix that is combined from all parallel simulations. The scheme allows for a practically ideal parallel scaling as no communication among the parallel simulations is required, and inactive cycles are not needed. (C) 2009 Elsevier Ltd. All rights reserved.

Place, publisher, year, edition, pages
2009. Vol. 36, no 8, 1276-1279 p.
National Category
Physical Sciences
Identifiers
URN: urn:nbn:se:kth:diva-14049DOI: 10.1016/j.anucene.2009.04.017ISI: 000269419800034Scopus ID: 2-s2.0-67651160330OAI: oai:DiVA.org:kth-14049DiVA: diva2:329375
Note

QC 20100709

Available from: 2010-07-09 Created: 2010-07-09 Last updated: 2014-09-24Bibliographically approved
In thesis
1. Development of New Monte Carlo Methods in Reactor Physics: Criticality, Non-Linear Steady-State and Burnup Problems
Open this publication in new window or tab >>Development of New Monte Carlo Methods in Reactor Physics: Criticality, Non-Linear Steady-State and Burnup Problems
2009 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

The Monte Carlo method is, practically, the only approach capable of giving detail insight into complex neutron transport problems. In reactor physics, the method has been used mainly for determining the keff in criticality calculations. In the last decade, the continuously growing computer performance has allowed to apply the Monte Carlo method also on simple burnup simulations of nuclear systems. Nevertheless, due to its extensive computational demands the Monte Carlo method is still not used as commonly as deterministic methods.

One of the reasons for the large computational demands of Monte Carlo criticality calculations is the necessity to carry out a number of inactive cycles to converge the fission source. This thesis presents a new concept of fission matrix based Monte Carlo criticality calculations where inactive cycles are not required. It is shown that the fission matrix is not sensitive to the errors in the fission source, and can be thus calculated by a Monte Carlo calculation without inactive cycles. All required results, including keff, are then derived via the final fission matrix. The confidence interval for the estimated keff can be conservatively derived from the variance in the fission matrix. This was confirmed by numerical test calculations of Whitesides's ``keff of the world problem'' model where other Monte Carlo methods fail to estimate the confidence interval correctly unless a large number of inactive cycles is simulated.

 

Another problem is that the existing Monte Carlo criticality codes are not well shaped for parallel computations; they cannot fully utilise the processing power of modern multi-processor computers and computer clusters. This thesis presents a new parallel computing scheme for Monte Carlo criticality calculations based on the fission matrix. The fission matrix is combined over a number of independent parallel simulations, and the final results are derived by means of the fission matrix. This scheme allows for a practically ideal parallel scaling since no communication among the parallel simulations is required, and no inactive cycles need to be simulated.

 

When the Monte Carlo criticality calculations are sufficiently fast, they will be more commonly applied on complex reactor physics problems, like non-linear steady-state calculations and fuel cycle calculations. This thesis develops an efficient method that introduces thermal-hydraulic and other feedbacks into the numerical model of a power reactor, allowing to carry out a non-linear Monte Carlo analysis of the reactor with steady-state core conditions. The thesis also shows that the major existing Monte Carlo burnup codes use unstable algorithms for coupling the neutronic and burnup calculations; therefore, they cannot be used for fuel cycle calculations. Nevertheless, stable coupling algorithms are known and can be implemented into the future Monte Carlo burnup codes.

 

Place, publisher, year, edition, pages
Stockholm: Universitetsservice US AB, 2009. xi, 49 p.
Series
Trita-FYS, ISSN 0280-316X ; 2009:20
Keyword
Monte Carlo, reactor physics, fission source, inactive cycles, convergence, burnup, steady-state, criticality, eigenvalue
National Category
Subatomic Physics
Identifiers
urn:nbn:se:kth:diva-10602 (URN)978-91-7415-366-8 (ISBN)
Public defence
2009-06-11, Sal FA32, AlbaNova, Roslagstullsbacken 21, Stockholm, 13:15 (English)
Opponent
Supervisors
Note
QC 20100709Available from: 2009-06-04 Created: 2009-06-01 Last updated: 2010-07-16Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full textScopus

Authority records BETA

Dufek, Jan

Search in DiVA

By author/editor
Dufek, JanGudowski, Waclaw
By organisation
Reactor Physics
In the same journal
Annals of Nuclear Energy
Physical Sciences

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 75 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf