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
Stochastic Approximation for Monte Carlo Calculation of Steady-State Conditions in Thermal Reactors
KTH, School of Engineering Sciences (SCI), Physics, Reactor Physics.ORCID iD: 0000-0002-7943-7517
KTH, School of Engineering Sciences (SCI), Physics, Reactor Physics.
2006 (English)In: Nuclear science and engineering, ISSN 0029-5639, E-ISSN 1943-748X, Vol. 152, 274-283 p.Article in journal (Refereed) Published
Abstract [en]

A new adaptive stochastic approximation method for an efficient Monte Carlo calculation of steady-state conditions in thermal reactor cores is described The core conditions that we consider are spatial distributions of power, neutron flux, coolant density, and strongly absorbing fission products like Xe-135. These distributions relate to each other; thus, the steady-state conditions are described by a system of nonlinear equations. When a Monte Carlo method is used to evaluate the power or neutron flux, then the task turns to a nonlinear stochastic root-finding problem that is usually solved in the iterative manner by stochastic optimization methods. One of those methods is stochastic approximation where efficiency depends on a sequence of stepsize and sample size parameters. The stepsize generation is often based on the well-known Robbins-Monro algorithm; however, the efficient generation of the sample size (number of neutrons simulated at each iteration step) was not published yet. The proposed method controls both the stepsize and the sample size in an efficient way; according to the results, the method reaches the highest possible convergence rate.

Place, publisher, year, edition, pages
2006. Vol. 152, 274-283 p.
Keyword [en]
ROBBINS-MONRO PROCEDURE
National Category
Physical Sciences
Identifiers
URN: urn:nbn:se:kth:diva-7363ISI: 000235833100003Scopus ID: 2-s2.0-33645138278OAI: oai:DiVA.org:kth-7363DiVA: diva2:12365
Note
QC 20100709Available from: 2007-07-10 Created: 2007-07-10 Last updated: 2017-12-14Bibliographically approved
In thesis
1. Advanced Monte Carlo Methods in Reactor Physics, Eigenvalue and Steady State Problems
Open this publication in new window or tab >>Advanced Monte Carlo Methods in Reactor Physics, Eigenvalue and Steady State Problems
2007 (English)Licentiate thesis, comprehensive summary (Other scientific)
Place, publisher, year, edition, pages
Stockholm: KTH, 2007. viii, 27 p.
Series
Trita-FYS, ISSN 0280-316X ; 2007:40
National Category
Subatomic Physics
Identifiers
urn:nbn:se:kth:diva-4458 (URN)978-91-7178-716-3 (ISBN)
Presentation
2007-06-07, FA32, AlvaNova, Roslagstullsbacken 21, Stockholm, 11:00
Note
QC 20101104Available from: 2007-07-10 Created: 2007-07-10 Last updated: 2010-11-04Bibliographically approved
2. 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

Scopushttp://epubs.ans.org/rip/index.cgi?ja=nse-152-3-274-283

Authority records BETA

Dufek, Jan

Search in DiVA

By author/editor
Dufek, JanGudowski, Waclaw
By organisation
Reactor Physics
In the same journal
Nuclear science and engineering
Physical Sciences

Search outside of DiVA

GoogleGoogle Scholar

urn-nbn

Altmetric score

urn-nbn
Total: 213 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