Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method

doi:10.19596/j.cnki.1001-246x.7570

Abstract

Abstract: A study of continuous energy adjoint weighted Monte Carlo method is carried out in several kinetic parameters, including effective delayed neutron fraction, prompt neutron generation time and prompt neutron decay constant, which is based on Iterated Fission Probability (IFP) method and selection of adjoint flux is extented. Difference among tracklength, collision, absorption and tally estimator are compared. Covariance and variance weight are also adopted based on these estimators. Meanwhile, different distribution domain of kinetic parameters are provided. IFP method is studied in more detail. Next event estimate and IFP method with different latent generation number are compared and influence on result are obtained. Fitness number of latent generation number is obtained. Particle kind and distribution information are automatically listed in output file.

Key words: continuous energy, Monte Carlo, adjoint weighted, kinetic parameters

CLC Number:

TL327

PENG Gang. Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 35(1): 87-94.

References

[1] KEEPIN G R. Physics of nuclear kinetics[M]. Addison-Wesley,1965.
[2] X-5 MONTE CARLO TEAM. MCNP-A General N-particle transport code[CP/OL]. Version 5.LA-UR-03-1987, Los Alamos National Laboratory[Apr 2003]. https://www.lanl.gov.
[3] BRETSCHER E E. Evalution of reactor kinetic parameters without the need for perturbation codes[C]. Proc Int Mtg Reduced Enrichment for Research and Test Reactors,JacksonHole,Wyoming,Oct 7-10,1997.
[4] STANKOVSKIY E Y. Calculation of the effective delayed neutron fraction with MCNP[J]. Trans Am Nucl Soc, 2009, 100:634-635.
[5] MEULEKAMP R K, STEVEN C, van der MARCK. Calculating the effective delayed neutron fraction with Monte Carlo[J]. Nucl Sci Eng,2006, 152:142-148.
[6] van der MARCK SC, MEULEKAMP R K. Calculating the effective delayed neutron fraction using Monte Carlo techniques[C]. PHYSOR 2004,The Physics of Fuel Cycles and Advanced Nuclear Systems:Global Developments, Chicago, Illinois, April 25-29, 2004.
[7] YASUSHI Nauchi,TAKAMORI Kameyama. Proposal of direct calculation of kinetic parameters β_eff and Λ based on continuous energy Monte Carlo method[J]. J Nucl Sci Techn,2005, 42:503-514.
[8] BRIAN C, Kiedrowski, et al.Calculating kinetics parameters and reactivity changes with continuous energy Monte Carlo[C]. Proc PHYSOR2010, Pittsburgh,Pennsylvania,USA,May 9-14,2010.
[9] JUDITH F, Briesmeister. MCNP^TM-A general Monte Carlo N-particle transporrt code[CP/OL]. Version 4C, LA-13709-M,March, 2000. https://mcnp.lanl.gov.
[10] BRIAN Kiedrowski. Adjoint-weighted kinetics parameters with continuous energy Monte Carlo[J]. Trans Am Nucl Soc,2009, 100:297-299.
[11] YASUNOBU Nagaya,KEN Nakajima. Approximate estimation of effective delayed neutron fraction with correated sampling method[C]. Proc PHYSOR2006,Interlaken,Switzerland,Sep 14-19,2008.
[12] STEPHANIE C Frankle. A suite of criticality benchmarks for validating nuclear data[EB/OL]. LA-13594,April 1999, https://www.lanl.gov.
[13] 徐传效,等.高通量工程试验堆(HFETR)运行十年论文集[M]. 四川科学技术出版社,1990.
[14] RAUCK S. Science V2 nuclear code package-qualification report[EB/OL]. FRAMATOME ANP,Ver.A,2004. https://www.areva.com.
[15] HALPERIN M. Almost linearly optimum combination of unbiased estimates[J]. Amer Stat Ass J, 1961, 56:36-43.
[16] SCHEFFH. On solutions of the Behrens Fisher problem based on the t-distribution[J]. Annals of Mathematical Statistics,1943,14:35-44.