FLUKA Monte Carlo Modelling of the CHARM Facility s Test Area: Update of the Radiation Field Assessment

Size: px

Start display at page:

Download "FLUKA Monte Carlo Modelling of the CHARM Facility s Test Area: Update of the Radiation Field Assessment"

Amberly Copeland
5 years ago
Views:

CERN-ACC-NOTE-2017-0059 14-11-2017 angelo.infantino@cern.

1 CERN-ACC-NOTE FLUKA Monte Carlo Modelling of the CHARM Facility s Test Area: Update of the Radiation Field Assessment Angelo Infantino CERN (EN-STI-FDA), Geneva, Switzerland Abstract The present Accelerator Note is a follow-up of the previous report CERN-ACC-NOTE In the present work, the FLUKA Monte Carlo model of CERN s CHARM facility has been improved to the most up-to-date configuration of the facility, including: new test positions, a global refinement of the FLUKA geometry, a careful review of the transport and physics parameters. Several configurations of the facility, in terms of target material and movable shielding configuration, have been simulated. The full set of results is reported in the following and can act as a reference guide to any potential user of the facility. Keywords: CHARM, FLUKA, R2E, Monte Carlo Simulation. This is an internal CERN publication and does not necessarily reflect the views of the CERN management.

2 Contents 1 Introduction FLUKA Monte Carlo model Geometry Physics and Transport Scoring Reference Spectra & Radiation Level Data References i

3 List of Figures 1 FLUKA model of the CHARM facility at target height (130 cm). Positions T0, PC0 and PCint are at a different height, not visible in the current section Physics and transport parameters of the 2016 FLUKA simulations Particle spectra for the configuration cpoooo - test location R Particle spectra for the configuration cpoooo - test location R Particle spectra for the configuration cpcioo - test location R Particle spectra for the configuration cpcioo - test location R Particle spectra for the configuration cpciic - test location R Particle spectra for the configuration cpciic - test location R Particle spectra for the configuration aloooo - test location R Particle spectra for the configuration aloooo - test location R Particle spectra for the configuration alcioo - test location R Particle spectra for the configuration alcioo - test location R Particle spectra for the configuration alciic - test location R Particle spectra for the configuration alciic - test location R Particle spectra for the configuration alhoooo - test location R Particle spectra for the configuration alhoooo - test location R Particle spectra for the configuration alhcioo - test location R Particle spectra for the configuration alhcioo - test location R Particle spectra for the configuration alhciic - test location R Particle spectra for the configuration alhciic - test location R This is an internal CERN publication and does not necessarily reflect the views of the CERN management.

4 List of Tables 1 Integral data for the configuration: cpoooo Hardness Factors for the configuration: cpoooo Integral data for the configuration: cpcioo Hardness Factors for the configuration: cpcioo Integral data for the configuration: cpciic Hardness Factors for the configuration: cpciic Integral data for the configuration: aloooo Hardness Factors for the configuration: aloooo Integral data for the configuration: alcioo Hardness Factors for the configuration: alcioo Integral data for the configuration: alciic Hardness Factors for the configuration: alciic Integral data for the configuration: alhoooo Hardness Factors for the configuration: alhoooo Integral data for the configuration: alhcioo Hardness Factors for the configuration: alhcioo Integral data for the configuration: alhciic Hardness Factors for the configuration: alhciic This is an internal CERN publication and does not necessarily reflect the views of the CERN management.

5 1 Introduction In this report, the latest update of the FLUKA [1, 2] Monte Carlo model of the CHARM [3] facility is reported. The update of the model involves mainly the geometry, the physics and transport parameters and the scoring. More information about the modelling and the CHARM facility can be found in the previous report [4], which the present Accelerator Note will refer to. 2 FLUKA Monte Carlo model 2.1 Geometry A detailed cross-check of the dimensions of the irradiation room has been conducted in order to review the full FLUKA Monte Carlo model. The review of the facility layout has been conduted through: a detailed cross check of the dimensions of the irradiaiton room through measurements (laser distance meter), a 3D model of the irradiaiton room obtained through a laser scan. Compared to the model reported in [4], the area of the patch panel, including the external corridoir, has been completely re-modeled as well as the shielded rest position of the target revolver for an accurate modelling of the back-scattered radiation from that area to a new test location, an Aluminum table called in the following T0. Several new test locations have been added to the model: a small Aluminum table, close to the target selector, for high dose rate tests (T0), test locations at the patch pannel area (G0 and G0old), test locations for dedicated power converter conveyer (PC0 1 and its itermediate position PCint). Figure 1 shows the up-to-date FLUKA model of the CHARM facility and all the test locations for which data are available. 2.2 Physics and Transport Physics and transport parameters have been modified accordingly to the latest studies conducted within the continuos CHARM calibration process. FLUKA simulations of the dose field within the irradiation room have been compared with experimental measurements performed by RadFET (RADiation sensing Field Effect Transistors) dosimeters in the standard test location (1-13). The deials of this study can be found in Ref. [5]. With regard to the purpose of the present report, the following parameters have been set for all the simulations: NEW-DEFA default transport threshold for Protons, Pions, Muons, Kaons set to 100 kev EMFCUT (production & transport) set to 100 kev for e, e +, γ Figure 2 shows the implementation of the physics and transport parameters. 2.3 Scoring To update the results already reported in up-to-date model of the CHARM facility: [4], the following quantities have been calcuted using the Dose High Energy Hadron fluence (HEH and HEH equivalent) Neutron Fluence 1 In other publication, PC0 might be referred as position G15. 1

Fig. 1: FLUKA model of the CHARM facility at target height (130 cm). Positions T0, PC0 and PCint are at a different height, not visible in the current section. Fig.

6 Fig. 1: FLUKA model of the CHARM facility at target height (130 cm). Positions T0, PC0 and PCint are at a different height, not visible in the current section. Fig. 2: Physics and transport parameters of the 2016 FLUKA simulations. Proton Fluence Thermal Neutron equivalent Silicon 1 MeV-equivalent R-factor The above quantities 2 have been calculated for all the standard and the new test locations. In addition to the integral data, particle spectra 2 have been assessed in the same locations. With regard to the position T0, 11, 12, 13 a dedicated set of scores have been set for studing the radiation field gradient in those areas. The comparison between FLUKA simulations and experimental measurements has been accurately reported in two dedicated publications (a peer-reviewed paper in Radiation Physics and Chemistry and a CERN EDMS Test Report) which will be available sooner. 3 Reference Spectra & Radiation Level Data In this section both reference spectra and radiation level data, for each configuration, are reported. Particle lethargy spectra in position R1 and R13 are reported as representative spectra of the lateral (R1-2 All the quantites - except the R-factor - are normalized per Proton On Target (POT) 2

7 R10) and the downstream (R11- R13) test locations respectively. Fig. 3: Particle spectra for the configuration cpoooo - test location R1. 3

8 Fig. 4: Particle spectra for the configuration cpoooo - test location R13. 4

9 Fig. 5: Particle spectra for the configuration cpcioo - test location R1. 5

10 Fig. 6: Particle spectra for the configuration cpcioo - test location R13. 6

11 Fig. 7: Particle spectra for the configuration cpciic - test location R1. 7

12 Fig. 8: Particle spectra for the configuration cpciic - test location R13. 8

13 Fig. 9: Particle spectra for the configuration aloooo - test location R1. 9

14 Fig. 10: Particle spectra for the configuration aloooo - test location R13. 10

15 Fig. 11: Particle spectra for the configuration alcioo - test location R1. 11

16 Fig. 12: Particle spectra for the configuration alcioo - test location R13. 12

17 Fig. 13: Particle spectra for the configuration alciic - test location R1. 13

18 Fig. 14: Particle spectra for the configuration alciic - test location R13. 14

19 Fig. 15: Particle spectra for the configuration alhoooo - test location R1. 15

20 Fig. 16: Particle spectra for the configuration alhoooo - test location R13. 16

21 Fig. 17: Particle spectra for the configuration alhcioo - test location R1. 17

22 Fig. 18: Particle spectra for the configuration alhcioo - test location R13. 18

23 Fig. 19: Particle spectra for the configuration alhciic - test location R1. 19

24 Fig. 20: Particle spectra for the configuration alhciic - test location R13. 20

25 21 Table 1: Integral data for the configuration: cpoooo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-01 R e e e e e e e e+00 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 G e e e e e e e e+00 G0old 2.529e e e e e e e e+00 PC e e e e e e e e+00 PCint 1.478e e e e e e e e+00 T e e e e e e e e-02 Internal Note

26 Table 2: Hardness Factors for the configuration: cpoooo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e+00 R e e e+00 R e e e-01 G e e e-02 G0old 4.002e e e-02 PC e e e-01 PCint 1.045e e e-01 T e e e-02 22

27 23 Table 3: Integral data for the configuration: cpcioo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-01 R e e e e e e e e+00 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 G e e e e e e e e+01 G0old 4.520e e e e e e e e+01 PC e e e e e e e e+01 PCint 2.397e e e e e e e e+01 T e e e e e e e e-02 Internal Note

28 Table 4: Hardness Factors for the configuration: cpcioo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+00 R e e e-01 G e e e-02 G0old 3.161e e e-02 PC e e e-02 PCint 7.699e e e-02 T e e e-02 24

29 25 Table 5: Integral data for the configuration: cpciic POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-01 R e e e e e e e e+00 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 G e e e e e e e e+01 G0old 2.622e e e e e e e e+01 PC e e e e e e e e+01 PCint 2.083e e e e e e e e+01 T e e e e e e e e-02 Internal Note

30 Table 6: Hardness Factors for the configuration: cpciic POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+00 R e e e-01 G e e e-02 G0old 3.122e e e-02 PC e e e-02 PCint 5.661e e e-02 T e e e-02 26

31 27 Table 7: Integral data for the configuration: aloooo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+00 G0old 9.156e e e e e e e e+00 PC e e e e e e e e+00 PCint 8.139e e e e e e e e+00 T e e e e e e e e-02 Internal Note

32 Table 8: Hardness Factors for the configuration: aloooo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 4.325e e e-02 PC e e e-01 PCint 1.070e e e-01 T e e e-01 28

33 29 Table 9: Integral data for the configuration: alcioo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+01 G0old 1.844e e e e e e e e+01 PC e e e e e e e e+00 PCint 1.932e e e e e e e e+00 T e e e e e e e e-02 Internal Note

34 Table 10: Hardness Factors for the configuration: alcioo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 3.136e e e-02 PC e e e-02 PCint 7.508e e e-02 T e e e-01 30

35 31 Table 11: Integral data for the configuration: alciic POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+01 G0old 1.084e e e e e e e e+01 PC e e e e e e e e+00 PCint 1.830e e e e e e e e+00 T e e e e e e e e-02 Internal Note

36 Table 12: Hardness Factors for the configuration: alciic POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 3.093e e e-02 PC e e e-02 PCint 5.489e e e-02 T e e e-01 32

37 33 Table 13: Integral data for the configuration: alhoooo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+00 G0old 5.275e e e e e e e e+00 PC e e e e e e e e+00 PCint 4.724e e e e e e e e+00 T e e e e e e e e-02 Internal Note

38 Table 14: Hardness Factors for the configuration: alhoooo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e-01 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 4.442e e e-02 PC e e e-01 PCint 1.071e e e-01 T e e e-01 34

39 35 Table 15: Integral data for the configuration: alhcioo POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+01 G0old 1.171e e e e e e e e+01 PC e e e e e e e e+00 PCint 1.169e e e e e e e e+00 T e e e e e e e e-02 Internal Note

40 Table 16: Hardness Factors for the configuration: alhcioo POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 3.094e e e-02 PC e e e-02 PCint 7.514e e e-02 T e e e-01 36

41 37 Table 17: Integral data for the configuration: alhciic POSITION DOSE HEH SI1MEV NEUTRON HEHEQ THNEUTRON PROTON RFACTOR [Gy P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] [cm 2 P OT 1 ] R e e e e e e e e+01 R e e e e e e e e+01 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e+00 R e e e e e e e e-01 R e e e e e e e e-02 R e e e e e e e e-01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+01 M e e e e e e e e+00 M e e e e e e e e+00 M e e e e e e e e-01 M e e e e e e e e-02 G e e e e e e e e+01 G0old 6.567e e e e e e e e+01 PC e e e e e e e e+00 PCint 1.104e e e e e e e e+00 T e e e e e e e e-02 Internal Note

42 Table 18: Hardness Factors for the configuration: alhciic POSITION H1 H10 H50 [GeV ] [GeV ] [GeV ] R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e-02 R e e e+00 R e e e+01 R e e e+00 G e e e-02 G0old 3.497e e e-02 PC e e e-02 PCint 5.212e e e-02 T e e e-01 4 References [1] A. Ferrari, P.R. Sala, A. Fasso, and J. Ranft, FLUKA: a multi-particle transport code, CERN (2005). [2] T.T. Bohlen, F. Cerutti, M.P.W. Chin, A. FassÃš, A. Ferrari, P.G. Ortega, A. Mairani, P.R. Sala, G. Smirnov and V. Vlachoudis, The FLUKA Code: Developments and Challenges for High Energy and Medical Applications, Nuclear Data Sheets 120, (2014). [3] J. Mekki, M. Brugger, R. G. Alia, A. Thornton, N. C. Dos Santos Mota, S. Danzeca, CHARM: A Mixed Field Facility at CERN for Radiation Tests in Ground, Atmospheric, Space and Accelerator Representative Environments, IEEE TRANSACTIONS ON NUCLEAR SCIENCE, 63(4), (2016). [4] A. Thornton, CHARM Facility Test Area Radiation Field Description, CERN-ACC-NOTE , (2015). [5] M. Marzo, S. Bonaldo, M. Brugger, S. Danzeca, R. Gacia Alia, A. Infantino, A. Thorton, RadFET dose response in the CHARM mixed-field: FLUKA MC simulations, Proceeding of the 13th International Conference on Radiation Shielding (ICRS-13) & 19th Topical Meeting of the Radiation Protection & Shielding Division of the American Nuclear Society (RPSD-2016) (2016) 38

Irradiation facilities at the INFN National Laboratory of Legnaro

Irradiation facilities at the INFN National Laboratory of Legnaro A. Candelori Istituto Nazionale di Fisica Nucleare and Dipartimento di Fisica, Padova 2/4/2004 OUTLINE The SIRAD irradiation facility at