MVP/GMVPII: General Purpose Monte Carlo Cock for Neutron and Photon Transport Calculations Based on Continuous Energy and Multigroup Method
Total Page:16
File Type:pdf, Size:1020Kb
00 CO JP0550301 MVP/GMVPII: General Purpose Monte Carlo Cock for Neutron and Photon Transport Calculations based on Continuous Energy and Multigroup Method June 2005 0 ffl % PJr Japan Atomic Energy Research Institute 0 A^7 tlHgis mm) fflill ®£ (Sf 3fH &$ It IE LilT ?flf i^ffl Japan Atomic Energy Research Institute Board of Editors Shinzaburo MATSUDA (Chief Editor) Zenko YOSHIDA Toshi NAGAOKA Shigeru TANAKA Tomihiro KAMIYA Kiyonobu YAMASHITA Atsushi TANAKA Yasuyuki ARATONO Tsutomu ISHIGAMI Kazuo ARAKAWA Hiroshi AIKAWA Hidetoshi AMANO Motoji NIIMI Kimiaki SAITOH Shusaku KOUNO Yukio TACHIBANA Yujiro IKEDA Yasuji MORITA Ryutaro HINO Shin-ichi SHAMOTO Osamu NAITO Akira NAGASHIMA Takamasa MORI Nobuyuki HOSOGANE Taikan HARAMI Hiroshi IKEZOE Kiyoshi OKUNO Osamu NARITA *°- h \t 319-1195 hx, ft JAERI reports are reviewed by the Board of Editors and issued irregularly. Inquiries about availability of the reports should be addressed to Research Information Division, Department of Intellectual Resources, Japan Atomic Energy Research Institute, Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195, Japan. Japan Atomic Energy Research Institute, 2005 0 JAERI1348 MVP/GMVPII: General Purpose Monte Carlo Codes for Neutron and Photon Transport Calculations based on Continuous Energy and Multigroup Methods Yasunobu NAGAYA, Keisuke OKUMURA, Takamasa MORI and Masayuki NAKAGAWA Department of Nuclear Energy System Tokai Research Establishment Japan Atomic Energy Research Institute Tokai-mura, Naka-gun, Ibaraki-ken (Received September 13, 2004) Abstract In order to realize fast and accurate Monte Carlo simulation of neutron and photon transport prob- lems, two vectorized Monte Carlo codes MVP and GMVP have been developed at JAERI. MVP is based on the continuous energy model and GMVP is on the multigroup model. Compared with conventional scalar codes, these codes achieve higher computation speed by a factor of 10 or more on vector super- computers. Both codes have sufficient functions for production use by adopting accurate physics model, geometry description capability and variance reduction techniques. The first version of the codes was released in 1994. They have been extensively improved and new functions have been implemented. The major improvements and new functions are (1) capability to treat the scattering model expressed with File 6 of the ENDF-6 format, (2) time-dependent tallies, (3) reaction rate calculation with the pointwise response function, (4) flexible source specification, (5) continuous- energy calculation at arbitrary temperatures, (6) estimation of real variances in eigenvalue problems, (7) point detector and surface crossing estimators, (8) statistical geometry model, (9) function of reactor noise analysis (simulation of the Feynman-o; experiment), (10) arbitrary shaped lattice boundary, (11) periodic boundary condition, (12) parallelization with standard libraries (MPI, PVM), (13) supporting many platforms, etc. This report describes the physical model, geometry description method used in the codes, new func- tions and how to use them. Keywords: Monte Carlo Method, Transport Calculation, Neutron, Photon, Continuous Energy, Multigroup, Nuclear Data Library, Vectorization, Parallel Computation, User's Manual JAERI 1348 MVP/GMVP II: IE (2004 ^ - H MVP t MVP (6) mm, , (8) i (10) fii ^ (11) 319-H95 JAERI 1348 in Contents 1 Introduction 1 2 Basics of Monte Carlo Methods 3 2.1 Particle Transport Calculation 3 2.2 Particle Transport with Monte Carlo Methods 5 2.3 Spatial Description 7 2.3.1 Spatial Description by Combinatorial Geometry 7 2.3.2 Repeated Geometry (LATTICE Geometry) 8 2.4 Boundary Conditions 9 2.5 Treatment of Nuclear Reactions 10 2.6 Estimation of Physical Quantities 11 2.6.1 Collision Estimator 11 2.6.2 Track Length Estimator 12 2.6.3 Point Detector Estimator 12 2.6.4 Surface Crossing Estimator 14 2.7 Statistical Uncertainty 14 2.8 History and Batch 15 2.9 Variance Reduction Techniques 16 2.9.1 Russian Roulette 18 2.9.2 Splitting 19 2.9.3 Path Stretching (Exponential Transform) 19 2.9.4 Options for Variance Reduction 20 2.10 Fixed-source and Eigenvalue Problems 21 3 Outline of Codes 23 3.1 Functions 23 3.2 Vector Calculation 24 3.3 Parallel Calculation 25 4 Physics Model 28 4.1 Representation of Cross Sections 28 4.1.1 Neutron Cross Sections 28 4.1.2 Photon Cross Sections 29 4.1.3 Preparation of Neutron Cross Sections at Arbitrary Temperatures 30 4.2 Analysis of Neutron Collision 31 4.2.1 Selection of Collision Nuclide 31 4.2.2 Generation of Fission Neutrons 32 4.2.3 Absorption 33 4.2.4 Scattering 33 4.2.5 Thermal Neutron Scattering 34 4.2.6 Speed of Neutrons 35 4.3 Estimation of Eigenvalue 35 4.4 Estimation of Real Variance in Eigenvalue Problems 36 4.5 Reactor Noise Analysis 38 iv JAERI1348 4.6 Photon Reaction 39 4.6.1 Detailed and Simple Models 40 4.6.2 Photoelectric Effect 40 4.6.3 Electron Pair Production 41 4.6.4 Incoherent (Compton) Scattering 41 4.6.5 Coherent (Thomson) Scattering 42 4.6.6 Bremsstrahlung Photon 42 4.6.7 Speed of Photons 42 MVP/GMVP Geometry 43 5.1 Combinatorial Geometry 43 5.2 BODIes 44 5.3 Definition of ZONEs 44 5.3.1 Operation of BODY Combination 44 5.3.2 Unsupported BODY Operation 45 5.3.3 OR at the Top of BODY Operations 45 5.3.4 INPUT-ZONE and ZONE 47 5.3.5 Multiply Defined ZONE 47 5.3.6 Identification of ZONEs and Undefined Regions 48 5.3.7 Finiteness of ZONEs 48 5.4 Specification of Materials and Boundary Conditions 48 5.5 Repeated Geometry 50 5.5.1 Concept of LATTICE Geometry 50 5.5.2 Hierarchical Structure of Multiple LATTICES 52 5.5.3 Definition of CELLs 53 5.5.4 Definition of LATTICES and Coordinate Transformation 54 5.5.5 Specification of CELLs Placed at SUBFRAMEs 57 5.5.6 Allocation of LATTICE at FRAME ZONE 60 5.5.7 Shape of FRAME ZONE 61 5.6 Statistical Geometry Model 62 5.6.1 How to Input the Statistical Geometry Model 63 5.7 Hierarchical Expression of REGION Names 65 5.8 TALLY REGION 66 5.9 REGION Name and REGION Number 67 Input Instructions 68 6.1 Features 68 6.2 Structure of Input Data 68 6.3 Effective Columns and Comments in a Line 69 6.4 Option Block 70 6.5 Format of Input Data Description 70 6.5.1 Labeled Free Format Input 71 6.5.2 Data Block 71 6.5.3 Labeled Data and Arithmetic Operation 72 6.5.4 Repetition and Interpolation of Values 74 6.5.5 Operations and Functions Applicable to Data Series 75 6.5.6 Reference of Variables Defined in the Codes as Symbolic Parameters 76 6.6 Data Input by REGION 76 JAERI1348 6.6.1 REGION-dependent Data 76 6.6.2 Input Format - 1 77 6.6.3 Input Format - 2 78 6.6.4 Specification of REGIONS by Wild Card Characters "*" and "?" 78 7 Title and Options 80 7.1 Title for Input Data 80 7.2 Rule for Option Block 80 7.3 List of Available Options 80 7.4 Conflict between Options 84 7.5 Detailed Information for Each Option 86 8 Calculation Control 105 8.1 Control Data for Both MVP and GMVP 105 8.2 Control Data Only for MVP 112 8.3 Input Data for PICTURE Option 115 9 Cross Section and Material Data 117 9.1 Input Data for MVP 117 9.1.1 CROSS SECTION or XSEC Data Block 117 9.1.2 Specification of Compositions (IDMAT and Number Density) 117 9.1.3 Specification of Nuclide ID in Photon Problems 119 9.1.4 Use of Arbitrary-temperature Neutron Cross Sections 119 9.2 Options for MVP Cross Section Data 120 9.3 Input Data for GMVP 123 9.3.1 Cross Section Form 123 9.3.2 CROSS SECTION Data Block 123 9.3.3 Input Parameters of Multigroup Cross Sections 124 9.3.4 Specification of Material Compositions 133 10 Geometry Data 134 10.1 Structure of Geometry Data 134 10.2 Definition of BODIes 134 10.2.1 BODY RPP (Right Parallelepiped) 136 10.2.2 BODY CYL (Cylinder in z-direction) 137 10.2.3 BODY RCC (Right Circular Cylinder) 138 10.2.4 BODY SPH (Sphere) 139 10.2.5 BODY BOX (Parallelepiped) 139 10.2.6 BODY WED or RAW (Wedge) 140 10.2.7 BODY ARB (Arbitrary Polyhedron with 4,5 or 6 Faces) 141 10.2.8 BODY HAF (Half Space with a Planer Surface) 142 10.2.9 BODY RHP (Right Hexagonal Prism) 143 10.2.10 BODY HEX (Hexagonal Prism) 143 10.2.11 BODY RCL (Right Cylinder with Specified Radial Direction) 144 10.2.12 BODY TRC (Truncated Right Cone) 145 10.2.13 BODY TEC (Truncated Elliptic Cone) 146 10.2.14 BODY ELL (Ellipsoid by Rotation) 147 10.2.15 BODY GEL (General Ellipsoid) 148 10.2.16 BODY ELT (Elliptical Torus) 148 vi JAERI 1348 10.2.17 BODY GQS (General Quadratic Curved Surface) 149 10.2.18 BODY BBC (BODY-by-BODY Combination) 150 10.3 Input Data for ZONEs 150 10.3.1 ZONE Name 151 10.3.2 REGION Name 151 10.3.3 Material ID 151 10.3.4 Combinatorial Geometry 152 10.3.5 'INPUT-ZONE' and 'ZONE' 152 10.3.6 ZONE Specification and Calculation Time 152 10.3.7 Periodic Boundary Condition 153 10.4 Input Data for LATTICE Geometry 154 10.4.1 Specification of LATTICE Structure 154 10.4.2 Examples of Input Data for LATTICE Geometry 159 10.4.3 Arrangement of LATTICES 162 10.5 Definition of CELLs 164 10.5.1 Definition of ZONEs in a CELL 165 10.6 Examples for LATTICE Geometry 165 10.7 Tallying in the Same CELLs 172 10.7.1 Labeling FRAME ZONEs 172 10.7.2 Labeling SUBFRAMEs 172 10.8 Input Data for Statistical Geometry Model 173 10.8.1 Definition of the STGM Region 174 10.8.2 Allocation of STGM Regions 175 10.8.3 Definition of STG Particles 176 10.8.4 Dummy CELL for MATRIX 177 10.8.5 Example for Multiple Kinds of STG Particles 177 10.9 Definition of TALLY REGIONS 179 11 Source Data 182 11.1 General Rules for Source Input Data 182 11.1.1 SOURCE Data Block 182 11.1.2 Set of Source Information 182 11.2 Representation of Sampling Information 184 11.3 Values and Mathematical Expression 184 11.4 Sampling with Functions