Initialize data. read input file. particle in system? START. find the collided. photon. neutron. nuclide and its in bank? in bank? interaction.
|
|
- Clinton Clark
- 5 years ago
- Views:
Transcription
1 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 of Science Technology and Environment, Thailand A Monte-Carlo photon and neutron transport code was developed at OAEP. The code was written in C and C++ languages in an object-oriented programming style. Constructive solid geometry (CSG), rather than combinatorial, was used such that making its input file more readable and recognizable. As the first stage of code validation, data from some EDF files, in the MCP's specific format, were used and compared with experimental data. The neutron (from a 00 mci Am/Be source) attenuation by water was chosen to compare the results. The agreement of the quantity 1=± among the calculation from SIPHO and MCP, and the experiment which are 10.9 cm, 9.71 cm and cm respectively was satisfactorily well within the experimental uncertainties. These results also agree with the 10.8 cm result of.m., Mirza, et al. 1. Intrduction We have developed our own Monte-Carlo particle transport code. The particles of our interest are neutron and photon, because they are available in our laboratory. We have used the reputed MCP code, including its ACE nuclear data libraries, as our benchmark. Obviously, the MCP is very long and complicated code because it was written in FORTRA and has an old-style programming. It could have taken us much effort to study MCP's source code directly. Instead, we have chosen to study the calculation method from other sources (including MCP manual, of course) and developed our own code using MCP as a guide. We chose C/C++ language because it has a neat structure, especially the OOP styleofc++. ot only the enhanced speed of code development, its source code reuseability is also important. We can add new features to the code with much less effort comparing to FORTRA. Some people may worry about speed penalty to choosing C/C++ over FORTRA, but we do not, because there are many good numerical libraries written in C/C++ that is comparable to that in FORTRA. In our case wehavechosen BLITZ++ (from because the vector operations are extensively used in our code. 2. The Code The name of the code is Simple Implementation of PHOton and eutron transport code (SIPHO). It is not final yet, but it has some major useful features. We have used the concept of Constructive Solid Geometry (CSG)[8] which define the geometry as a solid object. There are four basic shapes in SIPHO: cuboid, cylinder, sphere and cone. Objects of other shape could be obtained through Boolean operation between them. Some features of SIPHO are: ffl User interface via text input file (an example input file in this benchmark is shown below) ffl Repeated structure ffl Fixed particle sources ffl eutron and photon transport ffl Thermal neutron scattering treatment ffl Photon energy deposition in a cell ffl Track length estimate of particle flux in a cell The flowchart of the code is shown in fig.1. 1
2 Table 1: An example of input file objectf #declare srcwall = cylinderf<0,0,-.9>, <0,0,.9>,.87g mat<0> g #declare srcvial1 = cylinderf <0,0,-1.1>, <0,0,1.97>,1.1 g #declare srcvial2 = cylinderf <0,0,-.9>, <0,0,1.77>,.9 g objectf srcvial1 ~ srcvial2 mat<7> g #declare srctube = objectf #declare srcvoid1 = cylinderf <0,0,-15>, <0,0,2>, 2.g ~ #declare srcvoid2 = cylinderf <0,0,-14.8>, <0,0,2>, 2.1 g mat<4> g objectf srcvoid2 ~ srcvial1 mat<6> g objectf srcvial2 ~ srcwall mat<6> g #declare detvshroud = cylinderf <16.2,0,-1.5>, <16.2,0,1.5>, 1.1 g #declare uppercd1 = objectf #declare ucd1 = cylinderf <16.2,0,-1.5>, <16.2,0,-1>, 1.49 g ~ cylinderf <16.2,0,-1.5>, <16.2,0,-1>, 1.1 g mat<2> g #declare lowercd1 = objectf #declare lcd1 = cylinderf <16.2,0,1>, <16.2,0,1.5>, 1.49g ~ cylinderf <16.2,0,1>, <16.2,0,1.5>, 1.1g mat<2> g #declare knifetool1 = cylinderf <16.25,2,-14>, <16.25,2,14>, 2.1g #declare knifetool2 = cylinderf <16.25,-2,-14>, <16.25,-2,14>, 2.1g #declare knifetool = cylinderf <14.25,0,-14>, <14.25,0,14>, 2g #declare patchcd1 = objectf cylinderf <16.2,0,-1.5>, <16.2,0,-1>, 1.58g ~ ucd1 ~ knifetool1 ~ knifetool2 ~ knifetool mat<2> g #declare patchcd2 = objectf cylinderf <16.2,0,1>, <16.2,0,1.5>, 1.58 g ~ lcd1 ~ knifetool1 ~ knifetool2 ~ knifetool mat<2> g #declare dettube = objectf #declare detvoid1 = cylinderf <16.2,0,-15>, <16.2,0,2>, 1.9g ~ #declare detvoid2 = cylinderf <16.2,0,-14.8>, <16.2,0,2>, 1.7g mat<4> g objectf detvoid2 ~ uppercd1 ~ lowercd1 ~ patchcd1 ~ patchcd2 ~ detvshroud mat<6> g #declare detwall = objectf detvshroud ~ #declare detel = cylinderf <16.2,0, >, <16.2,0,10.795>, g mat<8> g #declare detector = objectf detel mat<5> g #declare splittank1 = cuboidf <-55,-55,-55>, <5.25,55,55> g #declare splittank2 = cuboidf < 5.25,-55,-55>, <8.25,55,55> g #declare splittank = cuboidf < 8.25,-55,-55>, <11.25,55,55> g #declare splittank4 = cuboidf < 11.25,-55,-55>, <55,55,55> g #declare tank = cylinderf <0,0,-16>, <0, 0,16>, 5g #declare watertank1 = objectf (tank & splittank1) ~ srcvoid1 mat<1> g #declare watertank2 = objectf tank & splittank2 mat<1> g #declare watertank = objectf tank & splittank mat<1> g #declare watertank4 = objectf (tank & splittank4) ~ detvoid1 mat<1> g objectf #declare otank = cylinderf tank mat<8> g objectf #declare env = spheref <0,0,0>,50 g ~ otank ~ srcvoid1 ~ detvoid1 mat<6> g objectf spheref <0,0,0>,55 g ~ env mat<0> g //+++ Last(data) Section mode = neutron source=fvolume<1>, vect=<1,0,0>, dir=<0>, erg=<d 1> g //Am/Be neutron energy distribution. cdf 1fi=[ ], p=[
3 -cont ]g mat 1[ lwtr.01t] //H2O mat 2[ ] //Cd mat [1.26e ] //B-10 mat 4[ ] // Fe mat 5[1.2875e ] //BF- mat 6[ ] // Air mat 7[ ] // Polyethylene mat 8[ ] // ss04 tallyf neutron flux within cell<detector> g tallyf neutron flux within cell <detector> multiplier<1,,107> g run particles imp[ ] START read input file Initialize data particle in system? find the collided photon neutron nuclide and its in bank? in bank? interaction photon running 1 calculate neutron production? mode? parameters then store in the bank running 1 neutron more than mode? production? 1 neutron? calculate photon parameters then scoring store in the bank Figure 1: The flowchart of SIPHO. Experiment The neutron attenuation in water was chosen to be the first validation scheme of the code. The experimental setup composed of a 00 mci Am/Be neutron source installed in the middle of a 70-cm-diameter water tank. A cadmium-covered BF detector was installed along side
4 of the source to measure the neutron activity. The distance of the detector to the source is allowed to vary from 2 cm to 26 cm with 2-cm-step increment. The thickness of cadmium is 1.8 mm and a 2-cm-wide window was provided at the middle of the detector. The window of the detector was aligned to the center of the source. The setup is shown in fig.2 and the detailed description of the detector is shown in fig. 70 cm 4.6 cm diam. 2 mm thick iron tube r Preamplifier 70 cm.8 cm diam. 2 mm thick iron tube HV supply 00 mci Am/Be neutron source MCA Cannberra S 10 distilled water Cd covered BF detector Stainless steel tank Figure 2: The experimental setup BF detector Detector Properties: 1.8 mm thick cadmium 2 cm window Overall length cm Sensitive length cm Outer shell...ss04, cm Fill pressure...40 cm Hg Fill gas...bf enriched to 96% B 10 Figure : The detector description 4. Calculation The model of the calculation at r = 16.2 cm is shown as the input file in table 1. In this calculation we have treat the implicit neutron capture, thermal neutron scattering and geometry splitting. At other r's the models are similar. 5. Results The results from the experiment and the calculation are shown in table 2. The plot of C B r 2 against r, where C B is the activation of boron within the detector and can be 4
5 calculated from C B = R 1 0 ± (n;ff) (E)ffi(E)dE, is also shown in fig.4. The source emission rate used was 96% of 6: n/s. r counting rate (n/s) C B = R 1 ± 0 (n;ff) (E)ffi(E)dE (cm) SIPHO MCP Experiment C B (rel.error) fom C B (rel.error) fom mean(rel.error) (0.0015) (0.024) (0.00) (0.01) (0.015) (0.006) (0.01) (0.018) (0.004) (0.011) (0.021) (0.0051) (0.01) (0.017) (0.0061) (0.0078) (0.015) (0.0075) (0.014) (0.016) (0.009) (0.016) (0.012) (0.01) (0.015) (0.01) (0.01) (0.024) (0.016) (0.016) (0.02) (0.02) (0.018) (0.02) (0.02).4 6.5(0.02) (0.026) (0.029) (0.027) Table 2: The counting rate from measurements and calculations Experiment, 1/Σ = cm SIPHO, 1/Σ = 10.9 cm MCP, 1/Σ = 9.71 cm C B r 2 (arbitrary unit) r (cm) Figure 4: Plot of C B r 2 against r to determine and compare 1/± 6. Conclusion It was seen from the plot in fig.4 that all curve have the same shape, though they are not exactly aligned. The determination of 1=± using the relation C B ο e ±r =r 2 or 5
6 r 2 C B ο e ±r at r > 1 cm, was found to be 10.9 cm for SIPHO, 9.71 cm for MCP and cm for the experiment. These values agree with 10.8 cm from the work of Mirza,.M., et al.[]. It should be noted that the Monte-Carlo efficiency of SIPHO is rather good, as shown by the figure of merit (fom) in the third column of table 2. However, there are somewhat difference between the result of SIPHO and MCP, which, unfortunately, we have not analysed in detail yet. Acknowledgements We would like to thank the free software community who contributed to every software we have used, especially BLITZ++, in our work. References [1] Lux, I. and L. Koblinger, Monte Carlo Particle Transport Methods: eutron and Photon calculations", CRC Press, Inc., Boca Raton, Florida. 517 p. (1990) [2] Cashwell, E. D. and C. J. Everett, Monte Carlo Method for Random Walk Problems", Pergamon Press, London. 15 p. (1959) [] Mirza,. M., S. M. Mirza and M. Iqbal, Determination of Mean Squared Slowing-down Distance for Am/Be eutrons in Water using BF -Detector", Radiat. Phys. Chem. vol. 48 no. 4 pp (1995) [4] Kludge, H. and K. Weise, The eutron Energy Spectrum of a 241 Am-Be(ff,n) Source and Resulting Mean Fluence to Dose EquivalentConversion Factors", Radiation Protection Dosimetry, vol. 2 o. 2 pp (1982) [5] Briesmeister, J. F., MCP A General Monte Carlo -Particle Transport Code Version 4B", Radiation Information Computational Center (RSICC), Oak Ridge, U.S.A., 709 p. (1997) [6] Eckel, B., Thinking in C++", Prentice Hall, ew Jersey, 81p. (1995) [7] Goldstein, H., Classical Mechanics", Addison-Wesley, U.S.A., pp (1980 2nd Ed.) [8] Foley, J. D., et al., Computer Graphics: Priciples and Practice", Addison-Wesley, U.S.A., pp (1992) 6
Application 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 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 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 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 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 informationValidation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition
Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition Bernd Schweizer, Andreas Goedicke Philips Technology Research Laboratories, Aachen, Germany bernd.schweizer@philips.com Abstract.
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 informationOutline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction
Monte Carlo Radiation Transport Modeling Overview () Lecture 7 Special Topics: Device Modeling Outline Principles of Monte Carlo modeling Radiation transport modeling with Utilizing Visual Editor (VisEd)
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 informationShielding factors for traditional safety glasses
Shielding factors for traditional safety glasses Malcolm McEwen, Hong Shen and Ernesto Mainegra-Hing Ionizing Radiation Standards, National Research Council Canada Alan DuSautoy, Radiation and Health Sciences
More informationArtifact Mitigation in High Energy CT via Monte Carlo Simulation
PIERS ONLINE, VOL. 7, NO. 8, 11 791 Artifact Mitigation in High Energy CT via Monte Carlo Simulation Xuemin Jin and Robert Y. Levine Spectral Sciences, Inc., USA Abstract The high energy (< 15 MeV) incident
More informationBreaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport
Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport GTC 2018 Jeremy Sweezy Scientist Monte Carlo Methods, Codes and Applications Group 3/28/2018 Operated by Los Alamos National
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 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 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 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 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 informationAnalysis of Radiation Transport through Multileaf Collimators Using BEAMnrc Code
American Journal of Biomedical Engineering 216, 6(4): 124-131 DOI: 1.5923/j.ajbe.21664.3 Analysis of Radiation Transport through Multileaf Collimators Using BEAMnrc Code Ankit Kajaria 1,*, Neeraj Sharma
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 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 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 informationMonte Carlo simulation of photon and electron transport
First Barcelona Techno Week Course on semiconductor detectors ICCUB, 11-15th July 2016 Monte Carlo simulation of photon and electron transport Francesc Salvat Monte Carlo 1 Simulations performed with the
More informationSuitability Study of MCNP Monte Carlo Program for Use in Medical Physics
Nuclear Energy in Central Europe '98 Terme Catez, September 7 to 10, 1998 SI0100092 Suitability Study of MCNP Monte Carlo Program for Use in Medical Physics R. Jeraj Reactor Physics Division, Jozef Stefan
More informationMonte Carlo Simulation for Neptun 10 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters
Monte Carlo Simulation for Neptun 1 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters M.T. Bahreyni Toossi a, M. Momen Nezhad b, S.M. Hashemi a a Medical Physics Research Center,
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 informationThe Monte Carlo simulation of a Package formed by the combination of three scintillators: Brillance380, Brillance350, and Prelude420.
EURONS I3 506065 JRA9 RHIB Report made during stay IEM-CSIC Madrid december 2006 MINISTERIO DE ASUNTOS EXTERIORES Y DE COOPERACIÓN AECI VICESECRETARÍA GENERAL The Monte Carlo simulation of a Package formed
More informationMCNP4C3-BASED SIMULATION OF A MEDICAL LINEAR ACCELERATOR
Computational Medical Physics Working Group Workshop II, Sep 3 Oct 3, 7 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (7) MCNP4C3-BASED SIMULATION
More informationMeasurement of depth-dose of linear accelerator and simulation by use of Geant4 computer code
reports of practical oncology and radiotherapy 1 5 (2 0 1 0) 64 68 available at www.sciencedirect.com journal homepage: http://www.rpor.eu/ Original article Measurement of depth-dose of linear accelerator
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 informationImprovement in the calibration techniques of approved dosimetry services in Spanish nuclear power plants of whole body counters with NaI detectors.
Improvement in the calibration techniques of approved dosimetry services in Spanish nuclear power plants of whole body counters with NaI detectors. E.Sollet (*), A.Felipe(**), JM.Perez (**), S. De Maria
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 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 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 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 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 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 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 informationDevelopment of a Variance Reduction Scheme in the Serpent 2 Monte Carlo Code Jaakko Leppänen, Tuomas Viitanen, Olli Hyvönen
Development of a Variance Reduction Scheme in the Serpent 2 Monte Carlo Code Jaakko Leppänen, Tuomas Viitanen, Olli Hyvönen VTT Technical Research Centre of Finland, Ltd., P.O Box 1000, FI-02044 VTT, Finland
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 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 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 informationMonte Carlo simulations
MC simulations Monte Carlo simulations Eirik Malinen Simulations of stochastic processes Interactions are stochastic: the path of a single ioniing particle may not be predicted Interactions are quantified
More informationProton dose calculation algorithms and configuration data
Proton dose calculation algorithms and configuration data Barbara Schaffner PTCOG 46 Educational workshop in Wanjie, 20. May 2007 VARIAN Medical Systems Agenda Broad beam algorithms Concept of pencil beam
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 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 informationIndrin Chetty Henry Ford Hospital Detroit, MI. AAPM Annual Meeting Houston 7:30-8:25 Mon 08/07/28 1/30
Review of TG105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning D. W. O. Rogers, Carleton Laboratory for Radiotherapy Physics,
More informationVIEWING MORSE-CG RADIATION TRANSPORT WITH. 3-D COLOR GRAPHICS* Abstract
SLAC-PUB-5170 January 1990 (1) VIEWING MORSE-CG RADIATION TRANSPORT WITH. 3-D COLOR GRAPHICS* YOSHIHITO NAMITO~, T.M. JENKINS AND W.R. NELSON Stanford Linear Accelerator Center Stanford University, Stanford,
More informationAccelerating koblinger's method of compton scattering on GPU
Available online at www.sciencedirect.com Procedia Engineering 24 (211) 242 246 211 International Conference on Advances in Engineering Accelerating koblingers method of compton scattering on GPU Jing
More informationEGS5 Monte Carlo Code: Installation
EGS5 Monte Carlo Code: Installation What do you need to install? EGS5 Cygwin CGView gvim (or any editor) EGS5 Monte Carlo Code: Installation Installing CYGWIN Installing EGS5 Running an EGS5 TEST Any trouble
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 informationI. INTRODUCTION. Figure 1. Radiation room model at Dongnai General Hospital
International Journal of Computational Engineering Research Vol, 04 Issue, 4 Simulation of Photon and Electron dose distributions 5 code for the treatment area using the linear electron accelerator (LINAC)
More informationStudy of scattered photons from the collimator system of Leksell Gamma Knife using the EGS4 Monte Carlo Code
Study of scattered photons from the collimator system of Leksell Gamma Knife using the EGS4 Monte Carlo Code Joel Y. C. Cheung Gamma Knife Centre, Canossa Hospital, 1 Old Peak Road, Hong Kong K. N. Yu
More informationImplementation of the EGSnrc / BEAMnrc Monte Carlo code - Application to medical accelerator SATURNE43
International Journal of Innovation and Applied Studies ISSN 2028-9324 Vol. 6 No. 3 July 2014, pp. 635-641 2014 Innovative Space of Scientific Research Journals http://www.ijias.issr-journals.org/ Implementation
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 informationMonte Carlo simulations. Lesson FYSKJM4710 Eirik Malinen
Monte Carlo simulations Lesson FYSKJM4710 Eirik Malinen MC simulations 1 Simulations of stochastic processes Interactions are stochastic: the path of a single ionizing particle may not be predicted Interactions
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 informationCALCULATION OF DOSE DISTRIBUTION AND DOSE RATE AT EACH DWELL POSITION IN 60 CO IRRADIATOR SVST-CO-60/B IN VIETNAM
CALCULATION OF DOSE DISTRIBUTION AND DOSE RATE AT EACH DWELL POSITION IN 60 CO IRRADIATOR SVST-CO-60/B IN VIETNAM Tran Van Hung Research and Development Center for Radiation Technology 202A, 11 Str., Linh
More informationCALCULATION AND MEASUREMENT OF EXPOSURE BUILDUP FACTORS FOR X-RAY RADIOGRAPHY
CALCULATION AND MEASUREMENT OF EXPOSURE BUILDUP FACTORS FOR X-RAY RADIOGRAPHY Taher 1. Aljundi and Joseph N. Gray Center for Nondestructive Evaluation Iowa State University Ames, Iowa 50011 INTRODUCTION
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 informationImproved Convergence Rates in Implicit Monte Carlo Simulations Through Stratified Sampling
Improved Convergence Rates in Implicit Monte Carlo Simulations Through Stratified Sampling ANS Winter Conference 2013 Alex Long and Ryan McClarren Texas A&M University Alex Long (Texas A&M) ANS Winter
More informationNeutronics 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 informationMonte Carlo methods in proton beam radiation therapy. Harald Paganetti
Monte Carlo methods in proton beam radiation therapy Harald Paganetti Introduction: Proton Physics Electromagnetic energy loss of protons Distal distribution Dose [%] 120 100 80 60 40 p e p Ionization
More informationCode characteristics
The PENELOPE Computer code M.J. Anagnostakis Nuclear Engineering Department National Technical University of Athens The PENELOPE code system PENetration and Energy LOss of Positrons and Electrons in matter
More informationGraphical User Interface for High Energy Multi-Particle Transport
Graphical User Interface for High Energy Multi-Particle Transport Phase I Final Report PREPARED BY: P.O. Box 1308 Richland, WA 99352-1308 PHONE: (509) 539-8621 FAX: (509) 946-2001 Email: randyschwarz@mcnpvised.com
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 informationSurface Area and Volume
Surface Area and Volume Level 1 2 1. Calculate the surface area and volume of each shape. Use metres for all lengths. Write your answers to 4 decimal places: a) 0.8 m Surface Area: Volume: b) 1 m 0.2 m
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 informationPreface. Med. Phys. 35(9), , Mechanical QA. Radiation Survey Mechanical tests Light radiation Table, Collimator, Gantry Jaws.
AAPM-SAM-2012-Das (1) Beam Data Collection and Commissioning for Linear Accelerators: Technical Considerations and Recommendations Preface Indra J. Das, PhD, FAAPM, FACR, FASTRO Department of Radiation
More informationDosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System
Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System A. Al-Basheer, M. Ghita, G. Sjoden, W. Bolch, C. Lee, and the ALRADS Group Computational Medical Physics Team Nuclear
More informationimass: Computational NRF Spectra Signal from Geant4
imass: Computational NRF Spectra Signal from Geant4 John Perry, Shanjie Xiao, Tatjana Jevremovic School of Nuclear Engineering, Purdue University, West Lafayette, IN 47907 {joperry, xiaosj, tatjanaj} @purdue.edu
More information11. Mensuration. Q 2 Find the altitude of a trapezium, the sum of the lengths of whose bases is 6.5 cm and whose area is 26 cm 2.
11. Mensuration Q 1 Find the volume of a cuboid whose length is 8 cm, breadth 6 cm and height 3.5 cm. Q 2 Find the altitude of a trapezium, the sum of the lengths of whose bases is 6.5 cm and whose area
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 informationImprovements in Monte Carlo simulation of large electron fields
SLAC-PUB-12908 October 2007 Improvements in Monte Carlo simulation of large electron fields Bruce A Faddegon 1, Joseph Perl 2 and Makoto Asai 2 1 University of California San Francisco Comprehensive Cancer
More informationDetailed analysis of scatter contribution from different simulated geometries of X-ray detectors
Detailed analysis of scatter contribution from different simulated geometries of X-ray detectors Elena Marimon 1,2*, Hammadi Nait-Charif 1, Asmar Khan 2, Philip A. Marsden 3, Oliver Diaz 4 1 Centre for
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 informationEvaluation of PBMR control rod worth using full three-dimensional deterministic transport methods
Available online at www.sciencedirect.com annals of NUCLEAR ENERGY Annals of Nuclear Energy 35 (28) 5 55 www.elsevier.com/locate/anucene Evaluation of PBMR control rod worth using full three-dimensional
More informationProduction of neutrons in laminated barriers of radiotherapy rooms: comparison between the analytical methodology and Monte Carlo simulations
JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, VOLUME 15, NUMBER 6, 2014 Production of neutrons in laminated barriers of radiotherapy rooms: comparison between the analytical methodology and Monte Carlo
More informationAssesing multileaf collimator effect on the build-up region using Monte Carlo method
Pol J Med Phys Eng. 2008;14(3):163-182. PL ISSN 1425-4689 doi: 10.2478/v10013-008-0014-0 website: http://www.pjmpe.waw.pl M. Zarza Moreno 1, 2, N. Teixeira 3, 4, A. P. Jesus 1, 2, G. Mora 1 Assesing multileaf
More informationX-Ray fluorescence and Raman spectroscopy
X-Ray fluorescence and Raman spectroscopy Advanced physics laboratory (nd part) 4CFU Catalini Letizia, De Angelis Giulia Vittoria, Piselli Verdiana Abstract In this paper we report about two different
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Monte Carlo Methods 1 Monte Carlo! Most accurate at predicting dose distributions! Based on
More informationSURFACE AREAS AND VOLUMES
CHAPTER 1 SURFACE AREAS AND VOLUMES (A) Main Concepts and Results Cuboid whose length l, breadth b and height h (a) Volume of cuboid lbh (b) Total surface area of cuboid 2 ( lb + bh + hl ) (c) Lateral
More informationChapter - 13 (Surface areas and Volumes)
Chapter - 13 (Surface areas and Volumes) Key Concepts SN. Name Figure Lateral/curved surface area 1 Cuboid Total surface area TSA Volume (V) Symbols use for b = breadth 2. Cube 4s 2 6s 2 s 3 s = side 3.
More informationFurther Volume and Surface Area
1 Further Volume and Surface Area Objectives * To find the volume and surface area of spheres, cones, pyramids and cylinders. * To solve problems involving volume and surface area of spheres, cones, pyramids
More informationJournal of American Science 2014;10(2)
LABVIEW Real Time Control for Pneumatic Transfer System in TRIGA Research Reactor Adel A.M. Abdelrahman Metallurgy Department, Nuclear Research Center, Atomic Energy Authority, Egypt elthakeel@yahoo.com
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 informationDEVELOPMENT AND CHARACTERISATION OF A HEAD CALIBRATION PHANTOM FOR IN VIVO MEASUREMENTS OF ACTINIDES
DEVELOPMENT AND CHARACTERISATION OF A HEAD CALIBRATION PHANTOM FOR IN VIVO MEASUREMENTS OF ACTINIDES 1 G. Gualdrini, 1 P. Battisti, 2 R. Biagini, 2 P. De Felice, 2 A. Fazio, 3 P. Ferrari 1; Ente Nazionale
More informationSupercomputing the Cascade Processes of Radiation Transport
19 th World Conference on Non-Destructive Testing 2016 Supercomputing the Cascade Processes of Radiation Transport Mikhail ZHUKOVSKIY 1, Mikhail MARKOV 1, Sergey PODOLYAKO 1, Roman USKOV 1, Carsten BELLON
More information9. Volume NOTES.notebook January 19, 2017
Starter - NO Calculators ) Find 20% of 248 2) Find - 5 4 Today's Learning: ) Find the highest common factor of 4 and 49. To revise volume of cubes, cuboids and prisms. 4) Find 5% of 80. 5) Calculate 256
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 informationCHAPTER 10: TALLYING IN MCNP
_or_.e_sa_m_h_us_se_in 6:..:...7 ---"M.=o:.c.;;nte.:-C-"-=arlo Particle Transport with MCNP CHAPTER 10: TALLYING IN MCNP Tallying is the process of scoring the parameters of interest, Le. providing the
More informationAn introduction to plotting data
An introduction to plotting data Eric D. Black California Institute of Technology February 25, 2014 1 Introduction Plotting data is one of the essential skills every scientist must have. We use it on a
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 informationGUIMesh: a tool to import STEP geometries into Geant4 via GDML
GUIMesh: a tool to import STEP geometries into Geant4 via GDML Abstract M. Pinto a *, P. Gonçalves a a LIP-Lisboa, Av. Gama Pinto, n.2, piso 3, 1649-003 Lisboa, Portugal Detailed radiation analysis of
More informationIMPROVING COMPUTATIONAL EFFICIENCY OF MONTE-CARLO SIMULATIONS WITH VARIANCE REDUCTION
International Conference on Mathematics and Computational Methods Applied to Nuclear Science & Engineering (M&C 2013), Sun Valley, Idaho, USA, May 5-9, 2013, on CD-ROM, American Nuclear Society, LaGrange
More informationA PRACTICAL LOOK AT MONTE CARLO VARIANCE REDUCTION METHODS IN RADIATION SHIELDING
A PRACTICAL LOOK AT MONTE CARLO VARIANCE REDUCTION METHODS IN RADIATION SHIELDING RICHARD H. OLSHER Health Physics Measurements Group, Los Alamos National Laboratory MS J573, P.O. Box 1663, Los Alamos,
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 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 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 informationA Monte Carlo shielding model for PET/CT clinics
A Monte Carlo shielding model for PET/CT clinics Robert L. Metzger and Kenneth A. Van Riper Radiation Safety Engineering, Inc Chandler, AZ Abstract Modern PET/CT clinics consist of a scanner room housing
More informationEffects of the difference in tube voltage of the CT scanner on. dose calculation
Effects of the difference in tube voltage of the CT scanner on dose calculation Dong Joo Rhee, Sung-woo Kim, Dong Hyeok Jeong Medical and Radiological Physics Laboratory, Dongnam Institute of Radiological
More information