Evaluation of PBMR control rod worth using full three-dimensional deterministic transport methods
|
|
- Brittany McCoy
- 6 years ago
- Views:
Transcription
1 Available online at annals of NUCLEAR ENERGY Annals of Nuclear Energy 35 (28) Evaluation of PBMR control rod worth using full three-dimensional deterministic transport methods Bismark Tyobeka a, Kostadin Ivanov a, *, Andreas Pautz b a Department of Mechanical and Nuclear Engineering, Pennsylvania State University, 23 Reber Building, University Park, State College, PA 682, United States b Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) mbh, Garching, Germany Received 5 October 27; accepted 6 November 27 Available online 8 January 28 Abstract It is a well known fact that the neutron diffusion theory fails in the vicinity of strongly absorbing media, such as the control rods. This failure is much more pronounced in the PBMR because the location of control rods in the side reflector adds a directional dependence to the flux, and this complicates the problem further. Reactor control and safety can only be ensured if control rod worths are accurately predicted. In this work a thorough evaluation of different control rod models is performed by constructing an approximate two-dimensional (2D) model adopting the so-called grey curtain and this is compared to full 3D deterministic transport model. Both these models are compared to an explicit MCNP model. It is shown in this study that it is possible to have an accurate 2D model of control rods, utilizing appropriate equivalent cross sections and applying them to a control rod grey curtain. Further, this paper also shows that it is possible to obtain reasonably accurate control rod worths using a 3D deterministic transport method with the same fidelity of a Monte Carlo reference calculation, provided that the correct cross sections are used. These results confirm that the deterministic transport models can be successfully used for PBMR transient safety analysis. Ó 27 Elsevier Ltd. All rights reserved.. Introduction The failure of diffusion theory to model highly absorbing regions is well known and numerous methods have been developed to overcome it by using the so-called equivalent diffusion parameters. In the PBMR design, the positioning of these highly absorbing regions in the side reflector, where the leakage out of the core adds a directional dependence to the flux, further complicates the problem. One of these methods is the Method of Equivalent Cross Sections (MECS), which was proposed by (Scherer and Neef, 976). The principle is to model the absorber and its environment in transport theory (S N ) and then extract cross sections and diffusion parameters from the transport solution that will represent the absorber region * Corresponding author. Tel.: ; fax: address: kni@psu.edu (K. Ivanov). accurately in subsequent 3D diffusion calculations. Another very simplified method available, called the equivalent boron concentration (EBC), consists of adding a certain amount of boron absorber homogeneously into a borated region representing the control rod (Reitsma and Naidoo, 23). The concentration of the boron is adjusted so that the control rod worth is conserved. The method assumes that the rod reactivity worth is known from experiment or from other methods such as MECS. In this work the 3D neutron transport S N code TORT (Johnson, 992), was used with the cross sections generated from MICROX-2 (Mathews et al., 997) to perform control rod worth calculations with control rods accurately and explicitly modeled in three-dimensions. The task here is to first get the accurate control cross sections and to apply them in the neutronics code, while at the same time coming up with the best possible geometric representation of the control rod in TORT. The main objective of these /$ - see front matter Ó 27 Elsevier Ltd. All rights reserved. doi:.6/j.anucene.27..8
2 B. Tyobeka et al. / Annals of Nuclear Energy 35 (28) studies is to obtain an optimum control rod representation in TORT, and use the differential control rod worth curve resulting from this configuration to evaluate and adjust the 2D control rod approximations (grey curtain) used in diffusion codes. In order to achieve this, it is important to first evaluate the accuracy of the grey curtain representation of control rods by comparing it to the control rod worths obtained using MCNP5 (RSIC Computer Code Collection, 23) as a reference and then proceed to develop an explicit 3D model in TORT and MCNP. The objective of this latter step is to evaluate the grey curtain representation in MCNP against the explicit MCNP control rod representation. 2. Core and control rod models developed In this study, three control rod models were developed and extensively evaluated for accuracy. The first one as mentioned above was the grey curtain representation. The grey curtain was used in DORT (r z geometry) and also used in TORT (r h z geometry) with a dummy theta dimension. This exercise was meant to verify TORT, because using the same cross sections and this dummy theta dimension, the TORT model is the same as the DORT model, hence one expects that the results obtained by the two codes should ideally be the same. This grey curtain model was also developed with MCNP, just for verification purposes. Secondly, an explicit model in MCNP was developed to represent the control rods explicitly, i.e. not as a grey curtain. Thirdly, a TORT model was developed and was made to be as close as possible to the MCNP to accurately model the control rod configuration of the PBMR 268 MW design (Reitsma et al., 24). Before proceeding to the models themselves, it is important to first discuss the PBMR control rod design. The control rods in the PBMR are located in holes in the side reflector. These sleeves are made of graphite and the control rod is an annulus of Boron carbide within a cladding, as shown in Fig.. Control Rod Channel Sleeve He There are 8 control rods and 7 shut down rods filled with small absorber sphere called KLAKS. For purposes of this study, the shut down elements were not modeled as these are not used during normal operation. Inner radius of control rod channel sleeve is 6.5 cm and the outer radius is 7.3 cm. In this study, the control rod can be inserted as far as 85 cm, which is the height of the active core. Again, in the actual PBMR design, control rods are divided into two sets, with nine inserted from the bottom and nine from the top, thereby creating a possibility of control rod overlap. For purposes of this work, all rods are assumed to be inserted from the top and no rod overlap is considered. In both TORT and MCNP, the design data used was obtained from the PBMR (Pty) Ltd Company in South Africa, and this is the same data used for the reactor design. Table gives a summary of the dimensions of the control elements while the number densities and material composition data are given in Table 2. From the data provided, an MCNP model was developed for Case N2 of the PBMR 268 MW Benchmark problem (Reitsma et al., 24). Since this case requires cross sections generated at 9 K and MCNP can perform calculations at room temperatures, further cross section processing was performed with the NJOY code (MacFarlane et al., 999) so as to finally produce the cross section file xsdir for MCNP at 9 K. As mentioned earlier, two MCNP models were developed, namely, MCNP-grey curtain and MCNPexplicit. These models are shown in Figs Table Control element design data Control rod and KLAK channels Other info Outer radius (cm) PCD 375 cm Number of 8 control rods Number of 7 KLAK channels Incoloy inner radius 4. cm Thickness (cm) Control rods Inner incoloy 4.2. B4C Outer incoloy 5.25 Fig.. PBMR control rod design. Inner Incoloy Outer Incoloy B4C Control Rod Front (Facing Reactor Core) Table 2 Incoloy (Inc8H) material composition of control rods Nuclide Weight percentage C-nat.7 Si-nat. P-3.3 S-32.2 Cr-nat 2. Mn-55.5 Fe-nat Co-59. Ni-nat 32.
3 52 B. Tyobeka et al. / Annals of Nuclear Energy 35 (28) 5 55 Control region In developing the TORT model of the PBMR, BOT3P, the preprocessor (and postprocessor) of DORT and TORT was used. This approach was chosen because it would eliminate the possibility of geometrical errors in the model. With BOT3P, the user can actually visualize the geometry of the model and one of the salient features of this tool is that it can make the necessary adjustment in case of simple mistakes like overlaps or mesh size. It also has features that ensure that the volume of the model is preserved in case of some approximations due to small mesh size. The resulting TORT model is shown in Fig Control rod worth calculations and results Fig. 2. Radial cut of the PBMR 268 MW MCNP model using grey curtain CR representation. Central reflector Mixing zone Fuel region Side Reflector Carbon Insulation Fig. 3. Axial cut of the PBMR 268 MW MCNP model using grey curtain CR representation. Fig. 4. Top view of the PBMR 268 MW MCNP (explicit) model showing positions of control rods. From the developed models above, core calculations were performed with a view of computing control rod worths for each model. Firstly, it was important to verify that with TORT and DORT, using the same cross sections and a dummy theta dimension in TORT, i.e., both using the grey curtain representation of the control region, the results will be the same as expected. Calculations were performed at various insertion depths of the grey curtain with steps of cm. The results obtained are shown in Fig. 6. In this figure, the ratio of rod reactivity worth at a given insertion step to rod reactivity worth when rods are fully inserted is plotted against where x is the insertion step or distance, and H is the height of the active core, in this case 85 cm. From the figure, it is clear that as expected, TORT in two-dimensions completely reproduced the DORT results. This helps to verify TORT for PBMR modeling, as there is no published literature that the authors are aware of in this regard. The same calculations were performed with MCNP using the grey curtain representation of the control region. The results of MCNP were compared with those of DORT and TORT and as shown in Fig. 7, other than the apparent shift as the control rods are inserted deeper into the core, there is a reasonable agreement between these results. Following from the reasonable agreement between the codes while using the grey curtain model of control rod representation, further calculations were performed with MCNP but this time, using the explicit model of control rods. The results were compared to the MCNP results with the grey curtain, to evaluate the accuracy of the grey curtain representation. The results, in Fig. 8, show some differences on the control rod worth between these two models, but the overall result is a reasonable agreement. One of the objectives of performing these calculations is to finally benchmark the DORT-TD code for PBMR and to later couple it with THERMIX-DIREKT. The coupled code will be used for transient analysis and one of the envisaged transient cases to be modeled with the coupled code is the rod ejection incident. It is therefore important that DORT-TD can model rod worth calculations accurately. In two-dimensions, the only feasible way to model control rods is by means of a grey curtain approximation. For this reason, rod worth calculations from DORT-TD
4 B. Tyobeka et al. / Annals of Nuclear Energy 35 (28) Control rods Fig. 5. Radial and axial cut of the PBMR 268 MW TORT model showing positions of control rods. are compared with MCNP-explicit model in Fig. 9. It can be concluded from the figure above, that the grey curtain model in DORT-TD agrees excellent with the results of MCNP using explicit representation of control rods for the PBMR 268 MW design. So far, all the results presented above have been obtained from symmetric rod movement, i.e. it is assumed that the whole bank of control rods is moved. It is also important to model an asymmetric movement of control rods and to be able to do that, a three-dimensional deterministic model is necessary. The TORT model presented in Fig. 5 was utilized for this purpose. Using the equivalent boron concentration method, calculations were initially performed with MCNP (explicit model) with all control PBMR 268MWth Control Rod Worth Curve - Grey Curtain Model (DORT-TORT Model Verification) TO RT DO RT-TD Fig. 6. Control rod worth results using the grey curtain model for DORT and TORT. PBMR 268MWth Control Rod Worth Curve - Grey Curtain Model (MCNP Reference Included).8 MCNP (Grey-Curtain) TORT DORT-TD Fig. 7. Control rod worth results comparing grey curtain representation with MCNP reference.
5 54 B. Tyobeka et al. / Annals of Nuclear Energy 35 (28) 5 55 PBMR 268MWth Control Rod Worth Curve (Grey-Curtain v/s Explicit Models).8 MCNP (Explicit) MCNP (Grey-Curtain) rho (x)/rho(h) Fig. 8. Control rod worth results evaluating the accuracy of the grey curtain representation. PBMR 268MWth Control Rod Worth Curve (Explicit v/s Grey Curtain Models).8 MCNP (Explicit) DORT-TD Fig. 9. Control rod worth results evaluating the accuracy of the DORT-TD grey curtain. rods inserted in the core and the resulting eigenvalue was noted. Subsequently, an iterative series of calculations was performed between MICROX, GIP and TORT, the goal being to find the eigenvalue in TORT which is reasonably close to the MCNP value. The Boron number density in MICROX was adjusted accordingly and cross sections were generated and arranged into a format readable by GIP. GIP was used to process the cross sections into the FIDO format used by TORT and core calculations were then performed in TORT. This process was repeated until the target k eff value was obtained in TORT. With the final cross sections obtained, the calculation is repeated with all rods withdrawn and the comparison between TORT and MCNP yielded good agreement. Finally, a series of calculations in both TORT and MCNP was performed where, a single rod and subsequently a combination of rods in an asymmetric arrangement were inserted whilst the rest were fully withdrawn and the corresponding reactivity worth was calculated. The results are shown in Table 3 and Table 3 Reactivity worth of a combination of control rods Combination MCNP (dk/k) TORT (dk/k) % Difference CR+CR4+CR CR+CR5+CR CR+CR CR+CR3+CR
6 B. Tyobeka et al. / Annals of Nuclear Energy 35 (28) PBMR 268MWth Single Rod Worth Calculations (CR).8 MCNP TORT Fig.. PBMR 268 MW single rod worth calculations for CR. Fig.. The results show a good agreement between TORT and MCNP for both single rod movement and for the case where a cluster of control rods moves in one part of the core. This demonstrates that deterministic transport methods when applied with the correctly computed cross sections can reproduce Monte Carlo results. 4. Conclusions Accurate prediction of control rod worths is imperative for the safety of all reactor types, especially for new designs where there exists little or no experimental or plant data, like the PBMR design. With the failure of diffusion theory in carrying out this task, it is traditional to make use of Monte Carlo calculations, or use approximations such as the Method of Equivalent Cross Sections or the EBC. In this work, it was shown that deterministic transport codes, which are readily available, can actually be used for this purpose with almost the same fidelity as the Monte Carlo calculations. It was shown, however, that these methods can only be useful if proper cross section generation tools are used. Furthermore, results from these models can be used to evaluate approximate models like the MECS and the grey curtain model of control rods, which in this case have shown to be quite accurate when compared to both MCNP and TORT. More importantly, the comparative analysis performed in the reported study given in this paper has indicated that the neutron transport models can be successfully utilized for both PBMR control rod worth calculation and transient safety analysis involving control rod movements. References Johnson, J.O., 992. A users manual for Mash., a Monte Carlo Adjoint Shielding Code System (Contains the documentation for DORT), Oak Ridge National Laboratory Report ORNL/TM-778. MacFarlane, R.E., NJOY 99.8, 999. Code system for producing pointwise and multi-group neutron and photon cross sections from ENDF/ B data, Los Alamos, New Mexico. Mathews, D., 997. An Improved Version of the MICROX-2 Code, Paul Scherrer Institute, Switzerland, PSI Bericht Nr Reitsma, F., Naidoo, D., 23. Evaluating the control rod modeling approach used in the South African PBMR: comparison of VSOP calculations with ASTRA experiments. Nuclear Engineering and Design 222 (2 3), Reitsma, F., de Haas, J.B., Tyobeka, B., Ivanov, K., De Cruz, F., 24. PBMR steady state and coupled kinetics core thermal-hydraulics test problems, in: Transactions of the ANS National Meeting, Pittsburg, PA, USA. RSIC Computer Code Collection, April 23, MCNP5 Monte Carlo N- Particle Transport Code System, Oak Ridge National Laboratory. Scherer, W., Neef, H.J., 976. Determination of equivalent cross sections for representation of control rod regions in diffusion calculations, July- 3.
Neutronics Analysis of TRIGA Mark II Research Reactor. R. Khan, S. Karimzadeh, H. Böck Vienna University of Technology Atominstitute
Neutronics Analysis of TRIGA Mark II Research Reactor R. Khan, S. Karimzadeh, H. Böck Vienna University of Technology Atominstitute 23-03-2010 TRIGA Mark II reactor MCNP radiation transport code MCNP model
More informationApplication of MCNP Code in Shielding Design for Radioactive Sources
Application of MCNP Code in Shielding Design for Radioactive Sources Ibrahim A. Alrammah Abstract This paper presents three tasks: Task 1 explores: the detected number of as a function of polythene moderator
More informationPSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code. Jaakko Leppänen
PSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code Jaakko Leppänen Outline Background History The Serpent code: Neutron tracking Physics and interaction data Burnup calculation Output
More informationBEAVRS benchmark calculations with Serpent-ARES code sequence
BEAVRS benchmark calculations with Serpent-ARES code sequence Jaakko Leppänen rd International Serpent User Group Meeting Berkeley, CA, Nov. 6-8, Outline Goal of the study The ARES nodal diffusion code
More informationCALCULATION OF THE ACTIVITY INVENTORY FOR THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER (MHH) IN PREPARATION FOR DISMANTLING THE FACILITY
CALCULATION OF THE ACTIVITY INVENTORY FOR THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER (MHH) IN PREPARATION FOR DISMANTLING THE FACILITY Gabriele Hampel, Friedemann Scheller, Medical University
More informationOPTIMIZATION OF MONTE CARLO TRANSPORT SIMULATIONS IN STOCHASTIC MEDIA
PHYSOR 2012 Advances in Reactor Physics Linking Research, Industry, and Education Knoxville, Tennessee, USA, April 15-20, 2012, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2010) OPTIMIZATION
More informationInstallation of a Second CLICIT Irradiation Facility at the Oregon State TRIGA Reactor
Installation of a Second CLICIT Irradiation Facility at the Oregon State TRIGA Reactor Robert Schickler and Steve Reese Oregon State University, 100 Radiation Center Corvallis, OR, 97330 USA Corresponding
More informationOPTIMIZATION OF MONTE CARLO TRANSPORT SIMULATIONS IN STOCHASTIC MEDIA
PHYSOR 2012 Advances in Reactor Physics Linking Research, Industry, and Education Knoxville, Tennessee, USA, April 15-20, 2012, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2012) OPTIMIZATION
More informationVerification of the Hexagonal Ray Tracing Module and the CMFD Acceleration in ntracer
KNS 2017 Autumn Gyeongju Verification of the Hexagonal Ray Tracing Module and the CMFD Acceleration in ntracer October 27, 2017 Seongchan Kim, Changhyun Lim, Young Suk Ban and Han Gyu Joo * Reactor Physics
More informationThe Pennsylvania State University. The Graduate School. Department of Mechanical and Nuclear Engineering
The Pennsylvania State University The Graduate School Department of Mechanical and Nuclear Engineering IMPROVED REFLECTOR MODELING FOR LIGHT WATER REACTOR ANALYSIS A Thesis in Nuclear Engineering by David
More informationMCNP Monte Carlo & Advanced Reactor Simulations. Forrest Brown. NEAMS Reactor Simulation Workshop ANL, 19 May Title: Author(s): Intended for:
LA-UR- 09-03055 Approved for public release; distribution is unlimited. Title: MCNP Monte Carlo & Advanced Reactor Simulations Author(s): Forrest Brown Intended for: NEAMS Reactor Simulation Workshop ANL,
More informationDaedeok-daero, Yuseong-gu, Daejeon , Republic of Korea b Argonne National Laboratory (ANL)
MC 2-3/TWODANT/DIF3D Analysis for the ZPPR-15 10 B(n, α) Reaction Rate Measurement Min Jae Lee a*, Donny Hartanto a, Sang Ji Kim a, and Changho Lee b a Korea Atomic Energy Research Institute (KAERI) 989-111
More informationClick to edit Master title style
Introduction to Serpent Code Fusion neutronics workshop, Cambridge, UK, June 11-12, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline style Serpent overview
More informationRadiological Characterization and Decommissioning of Research and Power Reactors 15602
Radiological Characterization and Decommissioning of Research and Power Reactors 15602 INTRODUCTION Faezeh Abbasi *, Bruno Thomauske *, Rahim Nabbi * RWTH University Aachen The production of the detailed
More informationClick to edit Master title style
New features in Serpent 2 for fusion neutronics 5th International Serpent UGM, Knoxville, TN, Oct. 13-16, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline
More informationSERPENT Cross Section Generation for the RBWR
SERPENT Cross Section Generation for the RBWR Andrew Hall Thomas Downar 9/19/2012 Outline RBWR Motivation and Design Why use Serpent Cross Sections? Modeling the RBWR Generating an Equilibrium Cycle RBWR
More information2-D Reflector Modelling for VENUS-2 MOX Core Benchmark
2-D Reflector Modelling for VENUS-2 MOX Core Benchmark Dušan Ćalić ZEL-EN d.o.o. Vrbina 18 8270, Krsko, Slovenia dusan.calic@zel-en.si ABSTRACT The choice of the reflector model is an important issue in
More informationABSTRACT. W. T. Urban', L. A. Crotzerl, K. B. Spinney', L. S. Waters', D. K. Parsons', R. J. Cacciapouti2, and R. E. Alcouffel. 1.
COMPARISON OF' THREE-DIMENSIONAL NEUTRON FLUX CALCULATIONS FOR MAINE YANKEE W. T. Urban', L. A. Crotzerl, K. B. Spinney', L. S. Waters', D. K. Parsons', R. J. Cacciapouti2, and R. E. Alcouffel ABSTRACT
More informationCPM-3 BENCHMARKING to the DOE/B&W CRITICAL EXPERIMENTS
CPM-3 BENCHMARKING to the DOE/B&W CRITICAL EXPERIMENTS Kenneth M. Smolinske and Rodney L. Grow Utility Resource Associates Corporation 1901 Research Boulevard, Suite 405 Rockville, Maryland 20850 ABSTRACT
More informationMethodology for spatial homogenization in Serpent 2
Methodology for spatial homogenization in erpent 2 Jaakko Leppänen Memo 204/05/26 Background patial homogenization has been one of the main motivations for developing erpent since the beginning of the
More informationDEVELOPMENT OF A GRAPHICAL USER INTERFACE FOR IN-CORE FUEL MANAGEMENT USING MCODE
Advances in Nuclear Fuel Management IV (ANFM 2009) Hilton Head Island, South Carolina, USA, April 12-15, 2009, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2009) DEVELOPMENT OF A GRAPHICAL USER
More informationA MULTI-PHYSICS ANALYSIS FOR THE ACTUATION OF THE SSS IN OPAL REACTOR
A MULTI-PHYSICS ANALYSIS FOR THE ACTUATION OF THE SSS IN OPAL REACTOR D. FERRARO Nuclear Engineering Department, INVAP S.E. Esmeralda 356 P.B. C1035ABH, C.A.B.A, Buenos Aires, Argentina P. ALBERTO, E.
More informationIMPROVEMENTS TO MONK & MCBEND ENABLING COUPLING & THE USE OF MONK CALCULATED ISOTOPIC COMPOSITIONS IN SHIELDING & CRITICALITY
IMPROVEMENTS TO MONK & MCBEND ENABLING COUPLING & THE USE OF MONK CALCULATED ISOTOPIC COMPOSITIONS IN SHIELDING & CRITICALITY N. Davies, M.J. Armishaw, S.D. Richards and G.P.Dobson Serco Technical Consulting
More informationOECD/NEA EXPERT GROUP ON UNCERTAINTY ANALYSIS FOR CRITICALITY SAFETY ASSESSMENT: CURRENT ACTIVITIES
OECD/NEA EXPERT GROUP ON UNCERTAINTY ANALYSIS FOR CRITICALITY SAFETY ASSESSMENT: CURRENT ACTIVITIES Tatiana Ivanova WPEC Subgroup 33 Meeting Issy-les-Moulineaux May 11, 2011 EG UACSA: Objectives Expert
More informationIMPLEMENTATION OF SALIVARY GLANDS IN THE BODYBUILDER ANTHROPOMORPHIC PHANTOMS
Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) IMPLEMENTATION
More informationTREAT Modeling & Simulation Using PROTEUS
TREAT Modeling & Simulation Using PROTEUS May 24, 2016 ChanghoLee Neutronics Methods and Codes Section Nuclear Engineering Division Argonne National Laboratory Historic TREAT Experiments: Minimum Critical
More informationELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS
Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) ELECTRON DOSE
More informationState of the art of Monte Carlo technics for reliable activated waste evaluations
State of the art of Monte Carlo technics for reliable activated waste evaluations Matthieu CULIOLI a*, Nicolas CHAPOUTIER a, Samuel BARBIER a, Sylvain JANSKI b a AREVA NP, 10-12 rue Juliette Récamier,
More informationWPEC - SG45: procedure for the validation of IRSN criticality input decks
WPEC - SG45: procedure for the validation of IRSN criticality input decks LECLAIRE Nicolas IRSN May, 14 th 2018 Contents 1. IRSN calculations with MC codes 2. Validation database 3. Procedure a) Construction
More informationImproved Detector Response Characterization Method in ISOCS and LabSOCS
P Improved Detector Response Characterization Method in ISOCS and LabSOCS *1 1 1 1 1 R. VenkataramanP P, F. BronsonP P, V. AtrashkevichP P, M. FieldP P, and B.M. YoungP P. 1 PCanberra Industries, 800 Research
More informationLimitations in the PHOTON Monte Carlo gamma transport code
Nuclear Instruments and Methods in Physics Research A 480 (2002) 729 733 Limitations in the PHOTON Monte Carlo gamma transport code I. Orion a, L. Wielopolski b, * a St. Luke s/roosevelt Hospital, Columbia
More informationExperience in Neutronic/Thermal-hydraulic Coupling in Ciemat
Madrid 2012 Experience in Neutronic/Thermal-hydraulic Coupling in Ciemat Miriam Vazquez (Ciemat) Francisco Martín-Fuertes (Ciemat) Aleksandar Ivanov (INR-KIT) Outline 1. Introduction 2. Coupling scheme
More informationGraphical User Interface for Simplified Neutron Transport Calculations
Graphical User Interface for Simplified Neutron Transport Calculations Phase 1 Final Report Instrument No: DE-SC0002321 July 20, 2009, through April 19, 2010 Recipient: Randolph Schwarz, Visual Editor
More informationStatus and development of multi-physics capabilities in Serpent 2
Status and development of multi-physics capabilities in Serpent 2 V. Valtavirta VTT Technical Research Centre of Finland ville.valtavirta@vtt.fi 2014 Serpent User Group Meeting Structure Click to of edit
More informationEXPERIENCE AND EVALUATION OF ADVANCED ON-LINE CORE MONITORING SYSTEM BEACON AT IKATA SITE
EXPERIENCE AND EVALUATION OF ADVANCED ON-LINE CORE MONITORING SYSTEM BEACON AT IKATA SITE Nobumichi Fujitsuka, Hideyuki Tanouchi, Yasuhiro Imamura, Daisuke MizobuchiI IKATA Power Station Shikoku Electric
More informationVerification of the 3D Method of characteristics solver in OpenMOC
Verification of the 3D Method of characteristics solver in OpenMOC The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published
More informationResearch Article Development and Application of MCNP5 and KENO-VI Monte Carlo Models for the Atucha-2 PHWR Analysis
Science and Technology of Nuclear Installations Volume 2, Article ID 68347, 7 pages doi:.55/2/68347 Research Article Development and Application of MCNP5 and KENO-VI Monte Carlo Models for the Atucha-2
More informationHELIOS CALCULATIONS FOR UO2 LATTICE BENCHMARKS
M-UR- 98-22. Title: Author@): Submitted to: HELOS CALCULATONS FOR UO2 LATTCE BENCHMARKS R. D. Mosteller nt'l Conf. on Physics of Nuclear Science & Technology slandia, Long sland, NY October 5-8, 1998 Los
More informationA FLEXIBLE COUPLING SCHEME FOR MONTE CARLO AND THERMAL-HYDRAULICS CODES
International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011) Rio de Janeiro, RJ, Brazil, May 8-12, 2011, on CD-ROM, Latin American Section (LAS)
More informationA premilinary study of the OECD/NEA 3D transport problem using the lattice code DRAGON
A premilinary study of the OECD/NEA 3D transport problem using the lattice code DRAGON Nicolas Martin, Guy Marleau, Alain Hébert Institut de Génie Nucléaire École Polytechnique de Montréal 28 CNS Symposium
More information1 st International Serpent User Group Meeting in Dresden, Germany, September 15 16, 2011
1 st International Serpent User Group Meeting in Dresden, Germany, September 15 16, 2011 Discussion notes The first international Serpent user group meeting was held at the Helmholtz Zentrum Dresden Rossendorf
More informationDose rate calculation at transport and storage casks for spent nuclear fuel
Dose rate calculation at transport and storage casks for spent nuclear fuel B. Gmal, U. Hesse, K. Hummelsheim, R. Kilger, M. Wagner Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) mbh, Forschungsinstitute,
More informationA Method for Estimating Criticality Lower Limit Multiplication Factor. Yoshitaka NAITO NAIS Co., Ltd.
A Method for Estimating Criticality Lower Limit Multiplication Factor Yoshitaka NAITO NAIS Co., Ltd. Progress Grade up of computer performance Sub-criticality becomes to be decided using computed results
More informationSubplane-based Control Rod Decusping Techniques for the 2D/1D Method in MPACT 1. Aaron M. Graham, Benjamin S. Collins, Thomas Downar
Subplane-based Control Rod Decusping Techniques for the 2D/1D Method in MPACT 1 Aaron M. Graham, Benjamin S. Collins, Thomas Downar Department of Nuclear Engineering and Radiological Sciences, University
More informationNUC E 521. Chapter 6: METHOD OF CHARACTERISTICS
NUC E 521 Chapter 6: METHOD OF CHARACTERISTICS K. Ivanov 206 Reber, 865-0040, kni1@psu.edu Introduction o Spatial three-dimensional (3D) and energy dependent modeling of neutron population in a reactor
More informationA COARSE MESH RADIATION TRANSPORT METHOD FOR PRISMATIC BLOCK THERMAL REACTORS IN TWO DIMENSIONS
A COARSE MESH RADIATION TRANSPORT METHOD FOR PRISMATIC BLOCK THERMAL REACTORS IN TWO DIMENSIONS A Thesis Presented to The Academic Faculty By Kevin John Connolly In Partial Fulfillment Of the Requirements
More informationAutomated ADVANTG Variance Reduction in a Proton Driven System. Kenneth A. Van Riper1 and Robert L. Metzger2
Automated ADVANTG Variance Reduction in a Proton Driven System Kenneth A. Van Riper1 and Robert L. Metzger2 1 White Rock Science, P. O. Box 4729, White Rock, NM 87547, kvr@rt66.com Radiation Safety Engineering,
More informationEvaluation of the Full Core VVER-440 Benchmarks Using the KARATE and MCNP Code Systems
NENE 2015 September 14-17 PORTOROŽ SLOVENIA 24th International Conference Nuclear Energy for New Europe Evaluation of the Full Core VVER-440 Benchmarks Using the KARATE and MCNP Code Systems György Hegyi
More informationAttila4MC. Software for Simplifying Monte Carlo. For more info contact or
Attila4MC Software for Simplifying Monte Carlo For more info contact attila@varian.com or Gregory.Failla@varian.com MCNP and MCNP6 are trademarks of Los Alamos National Security, LLC, Los Alamos National
More informationMCNP Variance Reduction technique application for the Development of the Citrusdal Irradiation Facility
IYNC 2008 Interlaken, Switzerland, 20-26 September 2008 Paper No. 376 MCNP Variance Reduction technique application for the Development of the Citrusdal Irradiation Facility R Makgae Pebble Bed Modular
More informationEvaluation of radiative power loading on WEST metallic in-vessel components
Evaluation of radiative power loading on WEST metallic in-vessel components M-H. Aumeunier 1, P. Moreau, J. Bucalossi, M. Firdaouss CEA/IRFM F-13108 Saint-Paul-Lez-Durance, France E-mail: marie-helene.aumeunier@cea.fr
More informationClick to edit Master title style
Fun stuff with the built-in response matrix solver 7th International Serpent UGM, Gainesville, FL, Nov. 6 9, 2017 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline
More informationMURE : MCNP Utility for Reactor Evolution - Description of the methods, first applications and results
MURE : MCNP Utility for Reactor Evolution - Description of the methods, first applications and results O. Méplan, A. Nuttin, O. Laulan, S. David, F. Michel-Sendis, J. Wilson To cite this version: O. Méplan,
More informationMCNP CLASS SERIES (SAMPLE MCNP INPUT) Jongsoon Kim
MCNP CLASS SERIES (SAMPLE MCNP INPUT) Jongsoon Kim Basic constants in MCNP Lengths in cm Energies in MeV Times in shakes (10-8 sec) Atomic densities in units of atoms/barn*-cm Mass densities in g/cm 3
More informationModeling the ORTEC EX-100 Detector using MCNP
Modeling the ORTEC EX-100 Detector using MCNP MCNP is a general-purpose Monte Carlo radiation transport code for modeling the interaction of radiation with materials based on composition and density. MCNP
More informationMonte Carlo Method for Medical & Health Physics
Med Phys 774 Monte Carlo Method for Medical & Health Physics Chapter 5. MCNP Monte Carlo Code 1 MCNP stands for: A general-purpose M C N-P code Particles that can be transported??? See the references:
More informationMicrowell Mixing with Surface Tension
Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical
More informationMultiphysics simulations of nuclear reactors and more
Multiphysics simulations of nuclear reactors and more Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg klasjareteg@chalmersse Division of Nuclear Engineering Department of Applied Physics Chalmers
More informationIndustrial Radiography Simulation by MCNPX for Pipeline Corrosion Detection
More Info at Open Access Database www.ndt.net/?id=18672 Industrial Radiography Simulation by MCNPX for Pipeline Corrosion Detection Mostafa KHODABANDELOU 1, Amir MOVAFEGHI 2 *, Peyman ROSTAMI 1 and Jahanghir
More informationCORE MONITORING EXPERIENCE WITH GARDEL
CORE MONITORING EXPERIENCE WITH GARDEL Axel Becker, Alejandro Noël Studsvik Scandpower GmbH Studsvik Scandpower Suisse GmbH Abstract The GARDEL core surveillance and analysis system is a standard, modular
More informationMonte Carlo simulation for adaptive optics
Monte Carlo simulation for adaptive optics R. Valicu 1, P. Böni 1,3, J. Stahn 2, U. Filges 2, T. Panzner 2, Y. Bodenthin 2, M. Schneider 2,3, C.Schanzer 3 1 Physik-Department E21, James-Franck-Strasse,
More informationSOFTWARE REQUIREMENTS SPECIFICATION FOR THE TRAC-M-SPECIFIC DATA MAP ROUTINE IN THE COUPLED TRAC-M/PARCS CODE. R. Matthew Miller, Thomas J.
PU/NE-00-9 SOFTWARE REQUIREMENTS SPECIFICATION FOR THE TRAC-M-SPECIFIC DATA MAP ROUTINE IN THE COUPLED TRAC-M/PARCS CODE R. Matthew Miller, Thomas J. Downar School of Nuclear Engineering Purdue University
More informationAnnals of Nuclear Energy
Annals of Nuclear Energy 55 (2013) 248 257 Contents lists available at SciVerse ScienceDirect Annals of Nuclear Energy journal homepage: www.elsevier.com/locate/anucene Visual numerical steering in 3D
More informationTHE SIGACE PACKAGE FOR GENERATING HIGH TEMPERATURE ACE FILES USER MANUAL
INTERNATIONAL ATOMIC ENERGY AGENCY NUCLEAR DATA SERVICES DOCUMENTATION SERIES OF THE IAEA NUCLEAR DATA SECTION IAEA-NDS-212 17 January 2005 THE SIGACE PACKAGE FOR GENERATING HIGH TEMPERATURE ACE FILES
More informationIdentification of Shielding Material Configurations Using NMIS Imaging
Identification of Shielding Material Configurations Using NMIS Imaging B. R. Grogan, J. T. Mihalczo, S. M. McConchie, and J. A. Mullens Oak Ridge National Laboratory, P.O. Box 2008, MS-6010, Oak Ridge,
More informationMedical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran.
DXRaySMCS First User Friendly Interface Developed for Prediction of Diagnostic Radiology X-Ray Spectra Produced by Monte Carlo (MCNP-4C) Simulation in Iran M.T. Bahreyni Toosi a*, H. Moradi b, H. Zare
More informationMiddle East Technical University Mechanical Engineering Department ME 413 Introduction to Finite Element Analysis Spring 2015 (Dr.
Middle East Technical University Mechanical Engineering Department ME 413 Introduction to Finite Element Analysis Spring 2015 (Dr. Sert) COMSOL 1 Tutorial 2 Problem Definition Hot combustion gases of a
More informationDevelopment of a Radiation Shielding Monte Carlo Code: RShieldMC
Development of a Radiation Shielding Monte Carlo Code: RShieldMC Shenshen GAO 1,2, Zhen WU 1,3, Xin WANG 1,2, Rui QIU 1,2, Chunyan LI 1,3, Wei LU 1,2, Junli LI 1,2*, 1.Department of Physics Engineering,
More informationAttenuation Coefficients for Layered Ceiling and Floor Shields in PET/CT Clinics
Attenuation Coefficients for Layered Ceiling and Floor Shields in PET/CT Clinics Robert L. Metzger and Kenneth A. Van Riper Radiation Safety Engineering, Inc 3245 North Washington Street, Chandler, AZ
More informationSimulation of Radiographic Testing for POD Assessment
4th European-American Workshop on Reliability of NDE - Th.5.B.2 Simulation of Radiographic Testing for POD Assessment Gerd-Rüdiger JAENISCH, Carsten BELLON, Uwe EWERT, BAM, Berlin, Germany Abstract. Efficient
More informationParticle track plotting in Visual MCNP6 Randy Schwarz 1,*
Particle track plotting in Visual MCNP6 Randy Schwarz 1,* 1 Visual Editor Consultants, PO Box 1308, Richland, WA 99352, USA Abstract. A visual interface for MCNP6 has been created to allow the plotting
More informationCombining Analytical and Monte Carlo Modelling for Industrial Radiology
19 th World Conference on Non-Destructive Testing 2016 Combining Analytical and Monte Carlo Modelling for Industrial Radiology Carsten BELLON, Gerd-Rüdiger JAENISCH, Andreas DERESCH BAM Bundesanstalt für
More informationInvestigation of mixing chamber for experimental FGD reactor
Investigation of mixing chamber for experimental FGD reactor Jan Novosád 1,a, Petra Danová 1 and Tomáš Vít 1 1 Department of Power Engineering Equipment, Faculty of Mechanical Engineering, Technical University
More informationPerformance Test of Core Protection and Monitoring System with DLL for SMART Simulator Implementation
2014 KNS Spring Meeting in Jeju Performance Test of Core Protection and Monitoring System with DLL for Simulator Implementation BonSeung Koo, DaeHyun Hwang and KeungKoo Kim May 29-30, 2014 Korea Atomic
More informationParallel computations for the auto-converted MCNP5 models of the ITER ECRH launcher
Forschungszentrum Karlsruhe in der Helmholtz-Gemeinschaft Parallel computations for the auto-converted MCNP5 models of the ITER ECRH launcher A. Serikov, U. Fischer, R. Heidunger, L. Obholz, P. Spaeh,
More informationFrequency-based method to optimize the number of projections for industrial computed tomography
More Info at Open Access Database www.ndt.net/?id=18739 Frequency-based method to optimize the number of projections for industrial computed tomography Andrea Buratti 1, Soufian Ben Achour 1, Christopher
More informationComparison of radiosity and ray-tracing techniques with a practical design procedure for the prediction of daylight levels in atria
Renewable Energy 28 (2003) 2157 2162 www.elsevier.com/locate/renene Technical note Comparison of radiosity and ray-tracing techniques with a practical design procedure for the prediction of daylight levels
More informationMATHEMATICAL MODELING OF SILICON SINGLE CRYSTAL INDUSTRIAL GROWTH
MATHEMATICAL MODELING OF SILICON SINGLE CRYSTAL INDUSTRIAL GROWTH Dr. Andris Muižnieks, Dr. Andis Rudevičs, Dr. Armands Krauze, BSc. Vadims Suškovs, BSc. Kirils Surovovs, BSc. Kārlis Janisels 1. Introduction
More informationGraphical User Interface Development for the MARS Code
GENES4/ANP2003, Sep. 15-19, 2003, Kyoto, JAPAN Paper 1093 Graphical User Interface Development for the MARS Code J.-J. Jeong, M. Hwang, Y. J. Lee, K. D. Kim, B. D. Chung Korea Atomic Energy Research Institute
More informationUsing three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model
Boundary Elements XXVII 245 Using three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model J. J. Rencis & S. R. Pisani Department of Mechanical Engineering,
More informationInitialize data. read input file. particle in system? START. find the collided. photon. neutron. nuclide and its in bank? in bank? interaction.
Application of EDF uclear Data for Testing a Monte-Carlo eutron and Photon Transport Code P. Siangsanan, W.Dharmavanij and S. Chongkum Physics Division, Office of Atomic Eneegy for Peace (OAEP), Ministry
More informationTHE SIMULATION OF THE 4 MV VARIAN LINAC WITH EXPERIMENTAL VALIDATION
2007 International Nuclear Atlantic Conference - INAC 2007 Santos, SP, Brazil, September 30 to October 5, 2007 ASSOCIAÇÃO BRASILEIRA DE ENERGIA NUCLEAR - ABEN ISBN: 978-85-99141-02-1 THE SIMULATION OF
More informationLarge Plastic Scintillation Detectors for the Nuclear Materials Identification System
Large Plastic Scintillation Detectors for the Nuclear Materials Identification System J.S. Neal, J.T. Mihalczo, M. T. Hiatt, J. D. Edwards Oak Ridge National Laboratory P. O. Box 2008, Oak Ridge, Tennessee
More informationModeling Radiation Transport Using MCNP6 and Abaqus/CAE Chelsea A. D Angelo, Steven S. McCready, Karen C. Kelley Los Alamos National Laboratory
Modeling Radiation Transport Using MCNP6 and Abaqus/CAE Chelsea A. D Angelo, Steven S. McCready, Karen C. Kelley Los Alamos National Laboratory Abstract: Los Alamos National Laboratory (LANL) has released
More informationISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry
ISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry Sasha A. Philips, Frazier Bronson, Ram Venkataraman, Brian M. Young Abstract--Activity measurements require knowledge of the
More informationTRANSX-2005 New Structure and Features R.E.MacFarlane Los Alamos National Laboratory
TRANSX-2005 New Structure and Features R.E.MacFarlane Los Alamos National Laboratory TRANSX-2005 is a translation of TRANSX to Fortran- 90/95 style with an extended code-management scheme. The new features
More informationElectron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations
Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations A. Al-Basheer, G. Sjoden, M. Ghita Computational Medical Physics Team Nuclear & Radiological Engineering University
More informationEvaluation of RAPID for a UNF cask benchmark problem
Evaluation of RAPID for a UNF cask benchmark problem Valerio Mascolino 1,a, Alireza Haghighat 1,b, and Nathan J. Roskoff 1,c 1 Nuclear Science & Engineering Lab (NSEL), Virginia Tech, 900 N Glebe Rd.,
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationWhite Paper 3D Geometry Visualization Capability for MCNP
White Paper 3D Geometry Visualization Capability for MCNP J. B. Spencer, J. A. Kulesza, A. Sood Los Alamos National Laboratory Monte Carlo Methods, Codes, and Applications Group June 12, 2017 1 Introduction
More informationModeling the White Sands Missile Range Fast Burst Reactor Using a Discrete Ordinates Code, PENTRAN
Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-23-207 Modeling the White Sands Missile Range Fast Burst Reactor Using a Discrete Ordinates Code, PENTRAN Taylor R. Schulmeister
More informationSALOME-CŒUR : une plate-forme pour des études neutroniques à EDF
SALOME-CŒUR : une plate-forme pour des études neutroniques à EDF Hadrien Leroyer, Renaud Barate 27 novembre 2014 Journée des Utilisateurs SALOME ENSTA - Saclay 2014 EDF. No partial distribution of information
More informationMulti-Disciplinary Design of an Aircraft Landing Gear with Altair HyperWorks
Multi-Disciplinary Design of an Aircraft Landing Gear with Altair HyperWorks Altair Engineering, October 2008 Introduction Task: To design an aircraft landing gear that meets design requirements of several
More informationDeliverable D10.2. WP10 JRA04 INDESYS Innovative solutions for nuclear physics detectors
MS116 Characterization of light production, propagation and collection for both organic and inorganic scintillators D10.2 R&D on new and existing scintillation materials: Report on the light production,
More informationEvaluation of RayXpert for shielding design of medical facilities
Evaluation of Raypert for shielding design of medical facilities Sylvie Derreumaux 1,*, Sophie Vecchiola 1, Thomas Geoffray 2, and Cécile Etard 1 1 Institut for radiation protection and nuclear safety,
More informationGeometric Templates for Improved Tracking Performance in Monte Carlo Codes
Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 2013 (SNA + MC 2013) La Cité des Sciences et de l Industrie, Paris, France, October 27-31, 2013 Geometric Templates
More informationMesh Human Phantoms with MCNP
LAUR-12-01659 Mesh Human Phantoms with MCNP Casey Anderson (casey_a@lanl.gov) Karen Kelley, Tim Goorley Los Alamos National Laboratory U N C L A S S I F I E D Slide 1 Summary Monte Carlo for Radiation
More informationRadiographic Simulator artist: Version 2
18th World Conference on Nondestructive Testing, 16-20 April 2012, Durban, South Africa Radiographic Simulator artist: Version 2 Carsten Bellon 1, Andreas Deresch 1, Christian Gollwitzer 1, Gerd-Rüdiger
More informationAdaptation of the Nagra Activation Analysis Methodology to Serpent
Adaptation of the Nagra Activation Analysis Methodology to Serpent Valentyn Bykov May 31, 2018, Serpent UGM 2018 Nuclear Power in Switzerland Reactor Type Net First MWe power Beznau 1 PWR 365 1969 Beznau
More informationSHIELDING DEPTH DETERMINATION OF COBALT PHOTON SHOWER THROUGH LEAD, ALUMINUM AND AIR USING MONTE CARLO SIMULATION
Research Article SHIELDING DEPTH DETERMINATION OF COBALT PHOTON SHOWER THROUGH LEAD, ALUMINUM AND AIR USING MONTE CARLO SIMULATION 1 Ngadda, Y. H., 2 Ewa, I. O. B. and 3 Chagok, N. M. D. 1 Physics Department,
More information