Shielding Calculation for the Target Station Design

Similar documents
Minimizing Activation, Size and Costs and yet Maximizing Efficiency

Assessment of Gamma Dose Rate for Hypothetical Radioactive Waste Container

On the neutron transparency of various materials proposed for the AIDA enclosure

Reduction of Backscattered Radiation in Enclosure X-ray Radiography

ESS Target Facility Design. John Haines for the ESS Target Project Team

The Pennsylvania State University. The Graduate School. College of Engineering SUPPORTING ANALYSIS FOR THE DEVELOPMENT OF A RECONFIGURABLE

Current Status of Pulsed Spallation Neutron Source of J PARC

IFMIF LITHIUM LOOP ACTIVATION AND ASSOCIATED OPERATOR DOSE RATES

Measurement of Radio-nuclides in Radioactive aerosols produced in a 120-GeV Proton Target Station

Planning for the Decommissioning of the ASTRA-Reactor

Available online at ScienceDirect. Procedia Materials Science 4 (2014 )

Backsplash study. Alex Moiseev 05/22/2001

production using Lighthouse accelerator J.W. Blokker

Neutron applications laboratory for ESS-BILBAO

Status Update 18MW ILC Beam Dump Design

High Performance Shielding Material to Prevent Radiation Leaks

Radiation at accelerator laboratories

LOS ALAMOS AQUEOUS TARGET/BLANKET SYSTEM DESIGN FOR THE ACCELERATOR TRANSMUTATION OF WASTE CONCEPT

CHARACTERISTICS OF SELF-POWERED NEUTRON DETECTORS USED IN POWER REACTORS

Estimation of the Dose Effects of Gamma-Ray Streaming Through Cracks in the MVST Storage/Disposal Cask

COLD NEUTRON SOURCE AT CMRR

SAMPLING AND RADIOLOGICAL ANALYSIS OF COMPONENTS OF THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER

Specification for Phase VII Benchmark

ANALYSIS OF FISSION RATIO DISTRIBUTION IN SPHERICAL LITHIUM METAL ASSEMBLY WITH A GRAPHITE REFLECTOR. August 1975

M. R. Gilbert, S. L. Dudarev, S. Zheng, L. W. Packer, and J.-Ch. Sublet. EURATOM/CCFE Fusion Association, UK

Journal of American Science 2014;10(2) Burn-up credit in criticality safety of PWR spent fuel.

KIPT ACCELERATOR DRIVEN SYSTEM DESIGN AND PERFORMANCE

The new material irradiation infrastructure at the BR2 reactor. Copyright 2017 SCK CEN

Analysis of Core Physics Test Data and Sodium Void Reactivity Worth Calculation for MONJU Core with ARCADIAN-FBR Computer Code System

RADIATION DOSE EVALUATION FOR HYPOTHETICAL ACCIDENT WITH TRANSPORT PACKAGE CONTAINING IRIDIUM-192 SOURCE

SHORT-TERM STORAGE CONSIDERATIONS FOR SPENT PLUTONIUM THORIUM FUEL BUNDLES

Production of Rhenium by Transmuting Tungsten Metal in Fast Reactors with Moderator

Neutron Transport and Material Activation in a Power Plant Based on the HCLL Blanket Concept

Energy Spectrum of Cosmic-Ray Photons Observed at the Ground Level in Japan

BINP accelerator based neutron source S. Taskaev

Burn up Analysis for Fuel Assembly Unit in a Pressurized Heavy Water CANDU Reactor

Wir schaffen Wissen heute für morgen

CONTAMINATED SCRAP METAL-DIFICULTIES TO DETECT. Paulo Fernando Lavalle Heilbron Filho Nerbe J. Ruperti Junior. Abstract

Conceptual Design of Temporally Storage Area in Hot Cell for Fusion DEMO Reactor )

Tools and applications for core design and shielding in fast reactors

RESEARCH REACTOR FRJ-1 (MERLIN) THE CORE STRUCTURES OF THE REACTOR BLOCK ARE DISMANTLED

WM2011 Conference, February 27 March 3, 2011, Phoenix, AZ

Soil, Groundwater and Cooling Water

Co γ -ray attenuation coefficient of barite concrete

REDUCING THE ACTIVATION OF THE IRIS REACTOR BUILDING USING THE SCALE/MAVRIC METHODOLOGY

NuMI/LBNF Radiation Protection. Jim Hylen NBI September 2017

Lecture (3) on. Nuclear Reactors. By Dr. Emad M. Saad. Mechanical Engineering Dept. Faculty of Engineering. Fayoum University

FUSION TECHNOLOGY INSTITUTE

ARIES-ACT-DCLL NWL Distribution and Revised Radial Build

Core Modification for the High Burn-up to Improve Irradiation Efficiency of the Experimental Fast Reactor Joyo

WM2011 Conference, February 27-March 3, 2011, Phoenix, AZ. The DN30 Package for the Transport of Enriched Reprocessed Uranium

Radiation and Enviromental studies for BDF

Ability of New Concept Passive-Safety Reactor "KAMADO" - Safety, Economy and Hydrogen Production -

Future technologies for the electrical isolation of III-V semiconductor devices. Final report on grant no. GR/N32969 (including GR/S52797/01)

Design and construction of a thermal neutron target for the RPI linac

DEMO Concept Development and Assessment of Relevant Technologies

SAFETY ANALYSIS AND ASSESSMENT FOR TEMPORARY STORAGE OF 106 IRRADIATED VVR-M2 FUEL TYPE ASSEMBLIES IN DALAT REACTOR

Anil Kumar a, *, Yujiro Ikeda b, Mahmoud Z. Youssef a, Mohamed A. Abdou a, Fujio Maekawa b, Yoshimi Kasugai b

Subcritical Experiments in Uranium-Fueled Core with Central Test Zone of Tungsten

OVERVIEW OF THE NSNS TARGET STATION* To be presented at the American Nuclear Society 1997 Annual Meeting Orlando, Florida June 1-5, 1997

Measuring Ytterbium Neutron Cross Section at 2.5 MeV. Norbert Hugger

Benchmark Specification for HTGR Fuel Element Depletion. Mark D. DeHart Nuclear Science and Technology Division Oak Ridge National Laboratory

Experimental and Theoretical Investigations of Ferro Boron as In-vessel Shield Material in FBRs

Study of mass attenuation properties for cobalt-free maraging steels at photon energies from 200 to 3000keV

Irradiation Damage Mechanism of Reactor Pressure Vessel Materials in PWR Nuclear Power Plant

1. Requirements 2. MARS Simulation for Beam Dump 3. Beam Dump Cooling 4. Layout of Beam Dump & Muon Pit 5. Summary & Schedule

M. Short (MIT) F. A. Garner (REC) M. B. Toloczko (PNNL) L. Shao, T. Chen, J. Gigax, E. Aydogan, C.-C. Wei (TAMU) V. N.

REPORT for TASK of the EFDA Technology Programme Reference: Document: Level of confidentiality Author(s): Date: Distribution list: Abstract:

The use of Serpent code for investigations at the University of Tartu

David A. McClintock Maxim N. Gussev

Irradiation Testing of Structural Materials in Fast Breeder Test Reactor

Validation of the Monte Carlo Code MVP on the First Criticality of Indonesian Multipurpose Reactor

BORAL CONTROL BLADE THERMAL-MECHANICAL ANALYSIS

Radiation Environment and Shielding Design Optimization at MEIC

Geosynchrotron emission from air showers at SLAC. Konstantin Belov UCLA

BORAL CONTROL BLADE THERMAL-MECHANICAL ANALYSIS

Comparative Analysis of ENDF, JEF & JENDL Data Libraries by Modeling the Fusion-Fission Hybrid System for Minor Actinide Incineration

A feasible DEMO blanket concept based on water cooled solid breeder

KIPT ADS Facility. Yousry Gohar 1, Igor Bolshinsky 2, Ivan Karnaukhov 3. Argonne National Laboratory, USA 2. Idaho National Laboratory, USA 3

Study on experiment and simulation shielding properties of security work chamber based on MCNP 1

Development of Nuclear Power Reactor Shielding Using Two Different Types of Heavy Concrete

J. R. Venhuizen D. W. Nigg G. E. Tripard. September 9, th International Congress On Neutron Capture Therapy

Design of Transport Container

Shielding evaluation and acceptance testing of a prefabricated, modular, temporary radiation therapy treatment facility

MYRRHA FUEL TRANSIENT TESTS PROJECT AT THE TRIGA-ACPR REACTOR

Definition of Small Gram Quantity Contents for Type B Radioactive Material Transportation Packages: Activity-Based Content Limitations, Rev.

Prospects of new ODS steels

Combined effect of molten fluoride salt and irradiation on Ni-based alloys

Hydrogen isotope retention in W irradiated by heavy ions and helium plasma

Richard C. McCall Stanford Linear Accelerator Center Stanford University, Stanford, California 94305

High-energy x-ray production with pyroelectric crystals

Calculation of Water Activity in Point 7

High-Energy Resolution Microcalorimeter EDS System for Electron Beam Excitation

ABSTRACT. 1. Introduction

Neutron Detector development at the ILL Based on 3 He and alternative convertors

MIT NUCLEAR REACTOR LABORATORY. Technical Analysis and Administrative Issues of Criticality Study for Different MITR Facilities

I?. 0. BOX2008, MS-6363 Oak Ridge, Tennessee , USA -

Conversion of MNSR (PARR-2) from HEU to LEU Fuel

A POSITRON SOURCE USING CHANNELING IN CRYSTALS FOR LINEAR COLLIDERS

Fuel and material irradiation hosting systems in the Jules Horowitz reactor

Transcription:

Shielding Calculation for the Target Station Design Contents Calculation Code & Nuclear Data Libraries Biological Shield 1-Dimensional 3-Dimensional Topics in RED will be presented. Neutron Beam Line (NBL) Methodology Junction of Bulk Shield and NBL Shielding Calculation for NBL Beam Stop Neutron Beam Shutter To-Chopper, Activation & Blade Size Fujio MAEKAWA, Masaya TAMURA (JAERI) Masayoshi KAWAI (KEK)

Code & Data MCNPX NMTC/JAM JENDL (<20MeV), LA-150 (<150MeV) Cutoff-γ file MCNP Cutoff-n file (< E tr ) Activation σ FENDL/A-2.0 (modified) Nuclide yield data file n γ n (> E tr ) γ n (< E tr ) DCHAIN-SP E tr = 20 or 150 MeV Radioactivity

Biological Shield (1-D): Objectives Conceptual design for the biological shield Learn about transmission phenomena of neutron and photon in shield materials Select appropriate shielding materials Determine approximate shield thickness

Biological Shield (1-D): Step-1 Source Term Calculation Model 3-dimensional TMRA in detail Mercury target, Lead reflector, Moderators, Beam holes, etc. Impinge a proton beam of 3GeV & 1MW to the mercury target Calculate neutron energy spectra leaking from the TMRA Every 15 degree angle ranges Total 12 (=180/15) spectra 30-45 deg. TMRA 15-30 deg. (120cm 120cm) Protom Beam 0-15 deg.

Biological Shield (1-D) : Source Spectra from TMRA Low-energy < 1 MeV almost isotropic High-energy > 10 MeV enhanced strongly to forward direction Leakage Neutron [n/ lethargy/ proton] 10 2 10 1 10 0 10-1 10-2 10-9 10-7 10-5 10-3 10-1 10 1 10 3 Neutron Energy [MeV] angle with respect to the proton beam direction 120-135 deg. 30-45 deg. 75-90 deg.

Biological Shield (1-D) : Step-2 Deep Penetration 1-Dimensional spherical model One of 12 source spectra Materials Pure iron Concrete Density: 7.2 g/cm 3 void included TMRA Source Thickness: 4.8m, 5.0m, 5.2m from the center Concrete Ordinary concrete (2.2 g/cm 3 ) Magnetite concrete (3.4 g/cm 3 ) Steel punched magnetite concrete (4.6 g/cm 3 ) Variance reduction by weight window CPU Time: 500 min. with Pentium-III (784 MHz) for 1 angle range Target dose: 0.1 µsv/h to achieve 12.5 µsv/h with considering a safety factor of ~10 a correction factor of ~10 for omission of streaming effects Iron 0 0.7 4.8 8m

Biological Shield (1-D) : Energy Spectra Flux [n/cm 2 /lethagy/s] @ 1 MW 10 14 10 12 10 10 10 8 10 6 10 4 10 2 10 0 10-2 10-4 10-6 Neutron Angle Range: 45-60 deg. Iron 0.8 m Iron Iron 1.8 m Iron 2.8 m Magnetite Concrete 5.3 m 3.8 m Magnetite Concrete 4.8 m 5.8 m Magnetite Concrete 6.8 m 10-9 10-7 10-5 10-3 10-1 10 1 10 3 Neutron Energy [MeV] Photon Angle Range: 45-60 deg. 0.8 m 1.8 m 2.8 m 3.8 m 4.8 m 5.3 m 5.8 m 6.8 m 10-2 10-1 10 0 10 1 10 2 Photon Energy [MeV] Monte Carlo simulation can be used for bulk shielding calculation.

Biological Shield (1-D) : Dose 10 12 Iron Magnetite Concrete Dose [µsv/h] 10 10 10 8 10 6 10 4 10 2 10 0 10-2 High-energy Neutron (> 1 MeV) Thermal neutron Low-energy Neutron (ev - kev) Photon 0.1 10-4 10-6 Angle Range: 45-60 deg. 0 100 200 300 400 500 600 700 800 Distance from the center [cm]

Biological Shield (1-D) : Results - Magnetite concrete is very effective to reduce total shield thickness when compared to ordinary concrete. Total shield thickness to attain 0.1 µsv/h (Iron shield up to 480 cm) 800 - No remarkable benefit is found in use of steel punched magnetite concrete. Magnetite concrete was selected for lateral outer region of the biological shield. Distance from the Center [cm] 750 700 650 600 550 500 Ordinary Concrete Magnetite Concrete Steel Punched Magnetite Concrete 0-15 15-30 30-45 45-60 60-75 75-90 90-105 105-120 120-135 135-150 150-165 165-180 Angle with Respect to the Proton Beam

Biological Shield (1-D) : Results Distance from the Center [cm] 800 750 700 650 600 550 Total shield thickness to attain 0.1 µsv/h (Magnetite Concrete) Iron Shield Thickness (Distance from the Center) < 480 cm < 500 cm < 520 cm Effective Dose [µsv/h] Breakdown 0.12 0.10 0.08 0.06 0.04 0.02 Iron: 480 cm from the center High-energy Neutron (> 1 MeV) Low-energy Neutron (ev - kev) Total Photon Thermal neutron 500 0.00 0-15 15-30 30-45 45-60 60-75 75-90 90-105 105-120 120-135 135-150 150-165 165-180 Angle with Respect to the Proton Beam 0-15 15-30 30-45 45-60 60-75 75-90 90-105 105-120 120-135 135-150 150-165 165-180 Angle with Respect to the Proton Beam

Biological Shield (1-D) : Optimization 10 13 10 11 10 9 Iron (7.2 g/cm 3 ) Low-energy Neutron (ev-kev) Ordinary (2.2 g/cm 3 ) 10 13 10 11 10 9 Iron (7.2 g/cm 3 ) Low-energy Neutron (ev-kev) Ordinary (2.2 g/cm 3 ) µ 10 7 10 5 10 3 10 1 10-1 Thermal Neutron Photon High-energy Neutron (> 1MeV) 1 Dose [µsv/hr] @ 1 MW 10 7 10 5 10 3 10 1 10-1 Thermal Neutron Photon High-energy Neutron (> 1MeV) 1 10-3 10-3 10-5 0 1 2 3 4 5 6 7 10-5 0 1 2 3 4 5 6 7 Distance from the Center [cm] - A thin concrete layer in the steel shield is very effective to attenuate low energy neutron fluxes. - This idea was adopted for the biological shield for the vertical direction.

Biological Shield (1-D) : Final Results Preliminary shield thickness was determined. Backward: Iron up to 4.8 m magnetite concrete up to 6.4 m Secondary gamma-rays are dominant to the total. Forward Iron up to 4.8 m magnetite concrete 6.5 ~ 8.0 m High-energy neutrons are dominant to the total. Plan View 45-60 deg. 60-75 deg. 75-90 deg. 90-105 deg. 105-120 deg. 30-45 deg. R7000 R6700 R6500 120-135 deg. R7300 R6400 15-30 deg. R6400 135-150 deg. R7600 R6400 150-165 deg. R8000 R6400 0-15 deg. Void Vessel R1500 TMRA R600 R6400 165-180 deg. R8000 Magnetite Concrete Density 3.4 g/cm 3 R6400 Iron Shield R4800 Proton Beam

Biological Shield (3-D): Introduction Objective Design a target station layout that satisfies the radiation dose regulation Approach Three-dimensional Monte Carlo calculation Consider various components in the target station in a model, especially major gaps and void spaces to treat streaming effects precisely Design items Materials & dimensions of shield blocks Location & dimensions of components Widths of gaps between neighboring components Shutter stroke, positions of top & bottom void for shutter movement Heating & activation in the shield region, etc.

Biological Shield (3-D): Model-1 Ceiling Shield Ordinary Concrete Steel Helium-Vessel Void Shutter (closed) Magnetite Concrete Shield Magnetite Concrete Shield

Biological Shield (3-D): Model-2 Hydrogen Transfer-Line (Void) Helium Vessel Concrete Layer Shutter Top Void Shield (Steel) Shield (Steel) Shutter (closed) Water Cooled Shield Neutron Beam Extraction Pipe Helium Vessel Support Cylinder

Biological Shield (3-D): Model-3 Outer-Plug Inner- Plug Gap (10 mm) Gap (50 mm) Helium Vessel Shield (Steel) Gap (50 mm) Shutter Moderators Target Reflector Void Void Water Cooled Shield Neutron Beam Extraction Pipe Gap (40 mm)

Helium Vessel Biological Shield (3-D): Model-4 Vertical Cut View Horizontal Cut View Hydrogen Transfer-Line Concrete Void Gap Gap Shutter (closed) Shutters (closed) Neutron Beam Extraction Pipe Moderators Target Reflector Water Cooled Shield Neutron Beam Extraction Pipe Moderator Reflector Water Cooled Shield

Biological Shield (3-D): Conditions Code & Data MCNPX 2.2.6 + LA-150 Density Usual steel: 7.7 g/cm 3 Steel shield blocks: 7.4 g/cm 3 (including unavoidable vacancies) Ordinary concrete: 2.2 g/cm 3 Magnetite concrete: 3.4 g/cm 3 Source 3 GeV, 1 MW proton beam on the Hg-target Design goal for dose rate: 1 µsv/h by neutrons > 10 MeV Correction for underestimation of high-energy neutron fluxes by MCNPX+LA150 for steel (underestimation of 20% / m, 4.8 m thickness, 0.8 4.8 ~ 1/3): 3 Correction for dose by neutrons < 10 MeV and photons: 2 A priori safety factor for MC calculation: 2 3 x 2 x 2 x 1 µsv/h = 12 µsv/h, not to exceed the dose limit of 12.5 µsv/h

Biological Shield (3-D): Results-1 Beam-Line #2 47.5 degrees with respect to the proton beam direction Streaming due to working area working area Streaming due to shutter top void Required shield thickness: 7.5 m Envelope of 1 µsv/h

Biological Shield (3-D): Results-2 Beam-Line #7 90.8 degrees with respect to the proton beam direction Streaming due to working area working area Streaming due to shutter top void Required shield thickness: 7.0 m Envelope of 1 µsv/h

Biological Shield (3-D): Results-3 Beam-Line #12 134.0 degrees with respect to the proton beam direction Streaming due to working area working area Streaming due to shutter top void Required shield thickness: 6.5 m Envelope of 1 µsv/h

Biological Shield (3-D): Horizontal Layout 7.0m 7.1m 7.2m 7.3m 7.4m 7.5m Magnetite Concrete Proton Beam Steel Steel Magnetite Concrete 7.4m 7.0m 7.1m 7.1m 7.2m 7.4m Thickness of 6.5m is required for backward angles, but increased to 7.0 m. - No strong request to reduce the thickness to 6.5 m from users - Save construction cost for future beam-line shielding

Biological Shield (3-D): Vertical Layout Shield (Ordinary Concrete) Shield (Iron) Outer Liner Shutter Drive Shield (Iron) Ordinary Concrete Layer (0.2m) 9.6m 3.5m 2.1m 0.5m Shutter 3m Hydrogen Moderator Mercury Target 5.1m FL=0 BL=0 Shield (Heavy Concrete) Vessel Support Cylinder Reflector Be+SS304 In-Vessel Shield Shield (Heavy Concrete) Neutron Beam Line 1.6m 1.6m BL 9.6m Base Plate

Biological Shield (3-D): Summary-1 Steel is basically used inside of the outer-liner (R < 4.8 m). A 200 mm thick horizontal layer of ordinary concrete is introduced in the steel shield at a height about B.L.+2500 mm to attenuate low-energy neutron flux to reduce radiation dose for shutter drives and activation levels of components in the working area around top of the helium-vessel by 2 ~ 3 orders of magnitude. For shutter bottom shield, steel is replaced by magnetite concrete (ρ = 3.4 g/cm 3 ) to reduce construction cost. For outside of the outer-liner, magnetite concrete is selected instead of ordinary concrete to reduce total shield thickness. (Use of the magnetite concrete is also required for adequate shielding against neutrons passing through neutron-beam-lines.) Total lateral shield thickness measured from the target station center ranges from 7.0 m to 7.5 m (4.2 ~ 4.4 kg/cm 2 ).

Biological Shield (3-D): Summary-2 Top of the target station is covered with fifteen ceiling shield blocks composed of 800 mm thick steel and 1400 mm thick ordinary concrete. Thick concrete shields for both upward and lateral directions surrounding the target station attenuate low-energy neutron fluxes below 1 MeV accumulated in the steel shields effectively. Height of the reflector plug is minimized to suppress its weight for easy remote handling under a condition that activation dose at top of the plug is lower than 100 µsv/h. Shutter stroke is determined as 400 mm (+100 mm margin). A beam-line is closed when a shutter goes down. Heights of shutter and shutter top shield are minimized for precise alignment and easy remote handling of shutters without increasing shield thickness. Finally, a target station layout that satisfied the radiation dose regulation was determined.

Methodology for NBL Calculation NBL Source Term Neutronics calculation for neutron beam lines (NBL) starting with a proton beam bombardment of the mercury target is not practical because a great amount of computation time is needed, & modeling of the whole system, TMRA, bulk shield and NBL, is too complicated. If an NBL source term that represents accurately neutrons traveling through a beam line is obtained, NBL calculation will be very efficient. 2-Step calculation Objective To create a methodology for NBL calculation a suitable NBL source term a subroutine to generate neutrons in Monte Carlo calculations

Methodology: Source Term Calculation Model TMRA & beam lines Calculate neutron current J(E,R) in an NBL at many positions from the moderator by using the DXTRAN spheres of the MCNPX code (source term) Two Duct Configurations 10 x 10 cm Straight Converging from 10 x 10 cm at the moderator to 2 x 2 cm at 11 m Beam Duct 10 x 10 cm Surfaces on which neutron current is calculated DXTRAN Spheres DXTRAN Spheres TMRA TMRA Proton Beam

Methodology: Calculated Results Neutron Current [n/p] 10x10 cm Straight Beam Duct 10 0 10-1 10-2 10-3 10-4 10-5 10-6 Wall Load 130(1/R 3 ) 65(1/R) 2 Current 10-3 10-4 10-5 10-6 10-7 10-8 10-9 0 1000 2000 3000 4000 5000 Distance from the Moderator [cm] Neutron Wall Load [n/cm/p] Neutron Current [n/p] 10 0 10-1 10-2 10-3 10-4 10-5 10-6 2x2 cm @ 11m Converging Beam Duct c: convergence rate (-8/1100 cm) X : viewed moderator size 0 (10 cm) Current 65(1/R+c/X 0 ) 2 130(1/R 3 +c/x 0 /R 2 ) Wall Load 0 500 1000 1500 2000 Distance from the Moderator [cm] 10-3 10-4 10-5 10-6 10-7 10-8 10-9 Neutron Wall Load [n/cm/p] Current: J(R) Wall Load: -dj(r)/dr Symbols: Monte Carlo calculation Lines: Theoretical estimation A decrease of neutron current between both ends of an interval just corresponds to a neutron wall load in the interval.

Methodology: Sampling X Moderator Area Wall Source 1. sample a position on a wall according to J(E,Z) 2. sample a position in the viewed moderator area uniformly 3. determine a direction as illustrated here 4. sample an energy according to J(E,Z) beam collimator Wall Source - beam line shield - shutter Duct Source 1. sample a position on a specific surface in the beam duct uniformly 2. sample a position in the viewed moderator area uniformly 3. determine a direction as illustrated here 4. sample an energy according to J(E,Z) Z Duct Source - To-chopper - beam stop - shutter - sample A source subroutines for MCNPX that generate source neutrons for NBL calculations by using the NBL source term was produced.

Methodology: Summary A new method for NBL shielding calculation has been established. Feature: In the 2nd step calculation, source neutrons are not generated on a certain plane across the NBL, but on beam duct walls along the NBL. This method enables us to treat a neutron wall load distribution, that is, a source term, accurately in the NBL shielding calculation.

Junction: Objectives Bulk shield Iron + Concrete Weak point in terms of beam-line shielding The concrete layer is needed for the bulk shield to stop low-energy neutrons penetrating through the thick iron shield. The concrete layer is, however, a weak point in terms of beam-line shielding. Objectives To determine an appropriate structure for the junction between the bulk shield and NBL α

Junction: Calculation R-Z model Iron for bulk shield up to 4.8 m from the center Concrete 1.3 or 1.9 m thickness Beam-line shield Iron (1.8m) + Polyethylene (0.3m) The NBL source term along the beam line Material Iron (7.2 g/cm 3, including void spaces) Concrete Ordinary concrete (2.2 g/cm 3 ) Magnetite concrete (3.5 g/cm 3 ) Steel punched magnetite concrete (4.6 g/cm 3 ) Polyethylene (0.9 g/cm 3 ) Target dose: 0.5 µsv/h by neutrons > 1 MeV assuming that the total dose < 1 µsv/h can be attained under this condition

Junction: 2D-Map Dose by Neutrons above 1 MeV, d205 400 10 3 Distance from the Beam-Line [cm] 300 200 100 0 Iron Bulk Shield Neutron Beam Duct Magnetite Concrete Iron Beam Line Shield 0.5 µsv/h Polyethylene Iron 300 400 500 600 700 800 10 2 10 1 10 0 10 1 Dose [µsv/h] Distance from the Moderator [cm]

Junction: Results 0m Ordinary Concrete 4.8m 6.7m Too High 0m Ordinary Concrete 4.8m 6.7m 4.2m 3.0m 2.4m 2.1m 2.6m 2.3m 2.1m 0m -1.6m 0m -1.6m No Concrete 0m Magnetite Concrete 4.8m 6.7m 0m Magnetite Concrete 4.8m 6.1m 0m Magnetite Concrete 4.8m 6.1m Suitable! 3.3m 2.7m 2.3m 2.1m 2.5m 2.3m 2.2m 2.1m 3.3m 2.8m 2.4m 2.2m 2.1m 0m -1.6m Too Thick 0m -1.6m No Concrete 0m -1.6m Steel Punched Magnetite Concrete 0m 4.8m 6.7m Steel Punched Magnetite Concrete 0m 4.8m 6.1m Steel Punched Magnetite Concrete 0m 4.8m 6.1m Suitable! 2.9m 2.5m 2.2m 2.1m 2.4m 2.3m 2.2m 2.1m 2.9m 2.7m 2.3m 2.2m 2.1m 0m -1.6m Too Thick 0m -1.6m No Concrete 0m -1.6m

Junction: Summary Magnetite concrete is suitable for the outer layer of the bulk shield to reduce additional shield on the NBL. 0m Bulk Shield 4.8m 6.5m Magnetite Concrete Iron Additional Shield Beam Line Shield 3.5m Iron 1.6m 0.0m 6.3m

NBL: Objectives To provide rough estimates of beam line shield thickness for arranging beam lines in the experimental hall Straight duct / Converging duct Dependence on duct size Dependence on distance from the moderator Empirical formula Detailed shielding design for each beam line --> Future work

NBL: Calculation R-Z model IRON + POLYETHYLENE (30 cm) The NBL source term along the beam line Material Iron (7.2 g/cm 3, including void spaces) Polyethylene (0.9 g/cm 3 ) Target dose: 0.2 µsv/h by neutrons > 10 MeV assuming that the total dose < 1 µsv/h can be attained under this condition An empirical formula was derived from the calculated results.

NBL: Empirical Formula T = + P X 2 0 1 530 21 ln 130 L + 3 10 T: NBL shield thickness [cm] P: proton beam power [MW] X0: viewed moderator size [cm] X1: converged beam size at the sample [cm] L: distance between the moderator and the sample [cm] C: convergence rate [ =(X1-X0)/L ] Moderator X0 C XL 0 2 Sample X1 This formula was used to estimate approximate dimensions of NBL shields for arranging NBLs in the experimental hall. L

NBL: Estimated Thickness 250 Straight Duct Continuously Converging Duct (Viewed Moderator Size: 10x10 cm) 250 Shield Thickness [cm] 200 150 100 8 x 8 cm 10 x 10 cm Shield Thickness [cm] 200 150 100 2x2cm @20m 2x2cm @10m 2x2cm @100m 50 50 2x2cm @50m 0 0 20 40 60 80 100 120 140 Distance from the Moderator [m] 0 0 20 40 60 80 100 Distance from the Moderator [m]

NBL: Validation The empirical formula predicts adequately the NBL shield thickness. Dose by Neutrons above 10 MeV Distance from the Beam-Line [cm] 200 100 Neutron Beam Duct Envelope of NBL Shield Determined with the Empirical Formula 0.2 µsv/h 0 500 1000 1500 2000 2500 3000 3500 4000 4500 Distance from the Moderator [cm] 10 3 10 2 10 1 10 0 10 1 Dose [µsv/h]

NBL: Validation == LASCE Experiment == The methodology will be validated by using the LANSCE neutron beam-line shielding experiment. Neutron Beam Duct

Beam Stop: Objectives To provide rough estimates of neutron beam stop size for arranging beam lines in the experimental hall Dependence on incident neutron flux Empirical formula Optimization of a neutron beam stop for each beam line --> Future work

Beam Stop: Calculation R-Z model IRON only Increase of shield thickness due to hydrogen-containing plates is assumed to be 15 % of iron thickness. The NBL source term along the beam line Material Iron (7.2 g/cm 3, including void spaces) Target dose: 0.2 µsv/h by neutrons > 10 MeV assuming that the total dose < 1 µsv/h can be attained under this condition An empirical formula was derived from the calculated results.

Beam Stop: Calculated Results Distance from the Beam-Line [cm] 200 180 160 140 120 100 80 60 40 20 Dose by Neutrons above 10 MeV @ 1.0 MW 10 7 200 1.0 MW 0.2 µsv/h 10 6 180 160 0.1 MW 10 5 140 10 4 120 neutron beam 0 1450 1500 1550 1600 1650 1700 1750 1800 Distance from the Moderator [cm] 10 3 10 2 10 1 10 0 10 1 Dose [µsv/h] Distance from the Beam-Line [cm] 100 80 60 40 20 Dose by Neutrons above 10 MeV @ 0.1 MW neutron beam 0.2 µsv/h 0 1450 1500 1550 1600 1650 1700 1750 1800 Distance from the Moderator [cm] 10 7 10 6 10 5 10 4 10 3 10 2 10 1 10 0 10 1 Dose [µsv/h] Distance from the Beam-Line [cm] 200 180 160 140 120 100 80 60 40 20 Dose by Neutrons above 10 MeV @ 0.01 MW 10 7 0.01 MW 10 Beam stop size depends on incident 6 neutron beam 0.2 µsv/h 0 1450 1500 1550 1600 1650 1700 1750 1800 Distance from the Moderator [cm] 10 5 10 4 10 3 10 2 10 1 10 0 10 1 Dose [µsv/h] neutron flux. Change of incident neutron flux to beam stops is simulated by changing proton beam power in the calculation.

Beam Stop: Empirical Formula Z: along the neutron beam line Shield Thickness = 322 + 20 ln (P) R: direction perpendicular to the neutron beam line Shield Thickness = 200 + 18 ln (P) P: Incident neutron flux 1.0 for 10x10cm straight duct and a beam stop at 15 m from the moderator Shield Thickness [cm] 300 250 200 150 100 R Z Neutron Beam Stop 50 y = 173.33 + 36log(x) R= 0.99883 y = 280 + 40log(x) R= 1 0 10-3 10-2 10-1 10 0 Proton Beam Power [MW] Example: 10x10cm straight duct @ 15 m from the moderator, proton beam power 1 MW Length: {322 + 20 ln(1.0)} x 1.5 = 483 cm (1.5=1.0+0.5, 1.0 for forward & 0.5 for backward ) Horizontal width: {200 + 18 ln(1.0)} x 2 = 400 cm (2 for left & right) Height: {200 + 18 ln(1.0)} + 175 = 375 cm (175 cm is the beam line height)

Shutter: Introduction Objectives To optimize shutter configurations, i.e., material, size and location To estimate radiation dose at a sample position Two Components to be Considered Beam duct component Uncollided part Collided part Bulk penetration component Beam Duct Component Bulk Penetration Component

Shutter: Need Hydrogen According to a rough estimation, radiation dose at a sample position will be ~ 100,000 µsv/h mainly due to low-energy neutrons when steel is used solely for a shutter core material. Combination of steel and hydrogen-contained material, such as polyethylene, is required for the shutter core.

Shutter: Bulk Penetration Component Calculation Source neutron for the 1-D bulk shielding calculation Monte Carlo calculation A spherical geometry with reflecting surface Point detector at the sample position at 11 m from the moderator Neutron Source Iron Concrete Shutter Region Point Detector @ 11 m Reflecting Surface (Dashed Red Lines)

Shutter: Bulk Penetration Component Total Shutter Length: 2.0 m (Iron 1.6 m + MAT-A 0.2 m + MAT-B 0.2 m) unit: µsv/h MAT-A MAT-B Thermal ev - kev > 1 MeV Photon Total PE Iron 0.6 385.0 0.6 6.3 392.5 B-PE Iron 0.8 388.8 0.6 5.1 395.3 Iron PE 35.5 78.1 0.6 69.4 183.5 Iron B-PE 1.1 79.1 0.6 19.0 99.7 PE SS-316 0.6 418.7 0.5 8.1 427.9 B-PE SS-316 0.6 434.5 0.6 6.2 441.9 PE Tungsten 0.4 297.3 0.3 8.0 305.9 B-PE Tungsten 0.5 305.1 0.4 7.6 313.5 Iron Tungsten 0.4 331.5 0.2 8.6 340.7 Tungsten Tungsten 0.5 325.6 0.2 8.8 335.0 Iron Ord. Conc. 60.1 257.2 0.5 33.4 351.2 Iron Mag. Conc. 20.7 390.0 0.5 44.5 455.7 Ord. Conc. Ord. Conc. 52.6 214.2 1.0 38.0 305.8 Mag. Conc. Mag. Conc. 18.0 280.4 0.8 47.4 346.7 PE PE 34.5 72.3 1.3 87.1 195.2 B-PE B-PE 1.0 59.6 1.4 18.4 80.3

Shutter: Bulk Penetration Component Results Low-energy neutrons and photons are the main contributor to the total dose. Since these low-energy neutrons travel all around the steel shield region, enhancement of the shutter core material is not effective to reduce the dose rate at a sample position. When polyethylene is used at the core end, the total dose is ~ 200 µsv/h. Hydrogen atoms in polyethylene attenuate effectively lowenergy neutrons penetrating through the steel shield blocks. Adding boron in polyethylene is further effective to reduce the dose due to thermal neutrons and photons, and the total dose is decreased to ~ 100 µsv/h.

Shutter: Beam Duct Component Uncollided Part The NBL source term Attenuation of neutron flux in the shutter due to total cross section was calculated by a simple exponential formula Collided Part The NBL source term Monte Carlo calculation with an R-Z model Point detector at the sample position at 11 m from the moderator Neutron Beam Iron Shield Neutron Beam Duct Shutter Region Point Detector @ 11 m

Shutter: Beam Duct Component unit: µsv/h Iron 180 cm + PE 20 cm Iron 220 cm + PE 20 cm Iron 180 cm + B-PE 20 cm Iron 220 cm + B-PE 20 cm Thermal ev - kev > 1 MeV Photon Total Uncollided 0.0 0.9 1.1 0.0 2.0 Collided 0.4 10.5 1.9 1.1 13.9 Total 0.4 11.4 3.0 1.1 15.9 Thermal ev - kev > 1 MeV Photon Total Uncollided 0.0 0.3 0.0 0.0 0.3 Collided 0.1 8.3 0.1 0.4 8.9 Total 0.1 8.6 0.1 0.4 9.2 Thermal ev - kev > 1 MeV Photon Total Uncollided 0.0 0.9 1.1 0.0 2.0 Collided 0.0 8.6 1.6 0.3 10.5 Total 0.0 9.5 2.7 0.3 12.5 Thermal ev - kev > 1 MeV Photon Total Uncollided 0.0 0.3 0.0 0.0 0.3 Collided 0.0 7.3 0.1 0.1 7.5 Total 0.0 7.6 0.1 0.1 7.8

Shutter: Beam Duct Component Results Dose rate due to the uncollided part is as small as 2 µsv/h for the 2.0 m length shutters. Collided low-energy neutron in ev-kev region is the main contributor to the dose rate at the sample position. Use of polyethylene at the core end is effective to reduce the total dose rate. When the total shutter length is 2.0 m (1.8 m steel + 0.2 m polyethylene), total dose rate is at most 20 µsv/h. Dose rate due to the beam duct component is not significant when it is compared with that due to the bulk penetration component.

Shutter: Summary A combination of steel (1.8 m) & polyethylene (0.2 m) was selected finally for the shutter core. No high density material like tungsten will be used. Although borated polyethylene is effective to reduce the dose rate at a sample position, it is not adopted because the material swells a little. A dose rate at a sample position will be at most 200 µsv/h. The dose rate can be reduced easily to 10 µsv/h by adding an auxiliary shield because low-energy neutrons and photons are dominant in the total dose rate.

To-Chopper: Activation Objectives To estimate radioactivity inventory in a chopper blade for maintenance To give a guideline for selecting a chopper blade material Calculation Conditions Blade size: 7 cm x 7 cm, 30 cm length Material: Nimonic-90 or Inconel X-750 Position: 6.2-6.5 m from the moderator Neutron beam size:5.64 cm x 5.64 cm Neutron flux: 4.1 x 10 9 n/s/cm 2 Operation: 1 MW, 5000 hours/year, 30 years

To-Chopper: Activation (Nimonic-90) 10 12 10 11 10 10 Nimonic-90 (Ni:53, Cr:20, Fe:5, Co:18, Ti:2.5, Al:1.5 wt%) 10 12 10 11 10 10 Nimonic-90 (Ni:53, Cr:20, Fe:5, Co:18, Ti:2.5, Al:1.5 wt%) without thermal neutron Total Activity in Blade [Bq] 10 9 10 8 10 7 10 6 Total Co-56 10 5 H-3 Co-57 V-49 Co-58 Cr-51 Co-60 10 4 Mn-54 Co-60m Mn-56 Ni-59 Fe-55 Ni-63 10 3 10-3 10-2 10-1 10 0 10 1 10 2 10 3 10 4 10 5 Cooling Time [Days] Total Activity in Blade [Bq] 10 9 10 8 10 7 10 6 Total Co-56 10 5 H-3 Co-57 V-49 Co-58 Cr-51 Co-60 10 4 Mn-54 Co-60m Mn-56 Ni-59 Fe-55 Ni-63 10 3 10-3 10-2 10-1 10 0 10 1 10 2 10 3 10 4 10 5 Cooling Time [Days]

To-Chopper: Activation (Inconel X-750) 10 12 10 11 10 10 Inconel X-750 (Ni:73, Cr:15.5, Fe:7, Ti:2.5, Al:0.7, Nb:1.0 wt%) 10 12 10 11 10 10 Inconel X-750 (Ni:73, Cr:15.5, Fe:7, Ti:2.5, Al:0.7, Nb:1.0 wt%) without thermal neutron Total Activity in Blade [Bq] 10 9 10 8 10 7 10 6 Total Co-58m 10 5 H-3 Co-60 Cr-51 Ni-59 Fe-55 Ni-63 10 4 Co-56 Nb-93m Co-57 Nb-94 Co-58 10 3 10-3 10-2 10-1 10 0 10 1 10 2 10 3 10 4 10 5 Cooling Time [Days] Total Activity in Blade [Bq] 10 9 10 8 10 7 10 6 Total Co-58m 10 5 H-3 Co-60 Cr-51 Ni-59 Fe-55 Ni-63 10 4 Co-56 Nb-93m Co-57 Nb-94 Co-58 10 3 10-3 10-2 10-1 10 0 10 1 10 2 10 3 10 4 10 5 Cooling Time [Days]

To-Chopper: Activation Summary Dominant radioactive nuclides are produced via the (n,γ) reactions. 59 Co(n,γ) 60 Co, 50 Cr(n,γ) 52 Cr, 62 Ni(n,γ) 63 Ni, 54 Fe(n,γ) 55 Fe Selecting a non-cobalt material is very effective to reduce total activity. The total activity in Inconel X-750 is about 1/10 to that in Nimonic-90. Since main nuclides in Inconel X-750 do not emit gamma-rays as intensely and/or energetically as 60 Co, radiation dose from the Inconel X-750 blade is further less than that from the Nimonic-90 blade. When a blade is made of Inconel X-750, total activity after 10 days cooling is approximately 10 9 Bq. This value is about 100 times less than the SNS estimation in a similar condition.

To-Chopper: Blade Size Objectives To give a guideline for a chopper blade size Large enough in terms of shielding for high-energy neutrons Small enough for short rise-time of neutron pulses and light weight Calculation Conditions Blade size: L x L x 30 cm, L = 0, 10, 10.5, 11, 11.5, 12, 14, 16 Material: Inconel X-750, 8.5 g/cm 3

To-Chopper: Blade Size Total Neutron Flux Neutron Flux > 1 MeV Neutron Flux > 10 MeV 10 10 10 10 8 10 7 10 9 Y-axis [cm] 5 0 5 10 8 Neutron Flux @ 1 MW [n/cm2/s] Y-axis [cm] 5 0 5 10 7 10 6 Neutron Flux @ 1 MW [n/cm2/s] Y-axis [cm] 5 0 5 10 6 10 5 Neutron Flux @ 1 MW [n/cm2/s] 10 7 10 10 5 0 5 10 10 10 5 0 5 10 10 10 5 0 5 10 X-axis [cm] X-axis [cm] X-axis [cm] 10 Neutron Flux > 100 MeV Y-axis [cm] 5 0 5 10 5 10 4 Neutron Flux @ 1 MW [n/cm2/s] Neutron flux distribution at the front end of the blade Low-energy neutron flux is well collimated while high-energy neutrons are scattered significantly. 10 3 10 10 5 0 5 10 X-axis [cm]

To-Chopper: Blade Size Results When the blade size is large enough, high-energy neutron flux attenuates more than 10 times, and low-energy neutron flux 10 3 ~ 10 4 times. When the blade size is just the same as the inner beam duct size, low-energy neutron flux increases more than 10 times. The blade size of 11 cm, i.e., 10 cm surrounded by an additional 0.5 cm layer, is enough. Neutron Flux / cm 2 / Lethargy 10 9 10 8 10 7 10 6 10 5 10 4 10 3 0.0 cm 10.0 cm 10.5 cm 11.0 cm 11.5 cm 16.0 cm 10 2 10-6 10-4 10-2 10 0 10 2 10 4 Neutron Energy [MeV] Neutron flux spectra on the 10 x 10 cm plane In the beam duct 1 m from the blade end as a function of blade size