Volume 53, Number 2, April-June 2018
|Page(s)||145 - 148|
|Published online||16 May 2018|
Comparison of mass attenuation coefficients of concretes using FLUKA, XCOM and experiment results
Department of Physics, Karnatak University,
2 Department of Nuclear Energy Engineering, Faculty of Engineering, Sinop University, 57000, Sinop, Turkey
* Corresponding author: email@example.com
Accepted: 24 February 2018
The mass attenuation coefficients of seven different types of normal and heavy concretes like ordinary, hematite-serpentine, ilmenite-limonite, basalt-magnetite, ilmenite, steel-scrap and steel-magnetite concretes has been simulated using FLUKA Monte Carlo code at high energies 1.5, 2, 3, 4, 5 and 6 MeV. The mass attenuation coefficients and linear attenuation coefficient of the concretes were found dependent upon the chemical composition, density and gamma ray energy. FLUKA Monte Carlo code results were found in good agreement with experimental and theoretical XCOM data. Our investigations for high energy gamma-ray interaction validate the FLUKA Monte Carlo code for use where experimental gamma-ray interaction results are not available.
Key words: gamma / mass attenuation coefficients / FLUKA / XCOM
© EDP Sciences 2018
Rapid increase in investigations of radiation interaction processes in computer environment using the Monte Carlo simulation has made easy nuclear engineering and technology. The Monte Carlo simulation is used in radiation transportation, shielding, detector response, medical applications and radio-biology, etc. The simulation process is the method for study of high-energy radiation interaction process where experiments are not possible. The shielding efficiency of a compound or a mixture is characterized by mass attenuation coefficients. The mass attenuation coefficients are a principle parameter for gamma-ray interaction. Other gamma-ray interaction parameters (e.g. half-value layer, tenth-value layer, effective atomic number and effective electron density) are being derived by using the mass attenuation coefficients. Therefore, it is essential to investigate the mass attenuation coefficients of materials used in shielding applications for nuclear reactors, accelerators, medical facilities, radiation protection and radiation dosimetry, etc. The radiation shielding has been a thrust area for optimization of radiation protection. The shielding materials are chosen as combination of low- and high-Z elements for gamma-ray and neutron (Bashter, 1997; Singh and Badiger, 2012). Various types of normal and heavy concretes shielding materials have been developed to minimize the construction cost of reactor biological shielding and containments with improved shielding efficiency materials (Bashter, 1997; Shirmardi et al., 2013). The glasses as transparent shielding materials are also being used in nuclear reactors (Singh et al., 2014). Various investigations have been done for calculation of shielding parameters and effectiveness of the different types of concretes (Makarious et al., 1988; Bashter et al., 1996; Makarious et al., 1996; Bashter, 1997; Singh and Badiger, 2012; Akkurt and El-Khayatt, 2013; Shirmardi et al., 2013; Singh et al., 2014). Bashter (1997) experimentally investigated the attenuation coefficients (shielding parameters) of concretes. The simulations of gamma-ray for shielding application are found elsewhere in literatures. The shielding for reactor core and accelerator are designed using the Monte Carlo code simulation results.
The gamma-ray interaction is characterized by partial interactions namely, photoelectric absorption, Compton scattering and pair production depending upon the energy and atomic number of the material or effective atomic numbers of the compound or mixture. The theoretical values for mass attenuation coefficients and cross sections for different elements, compounds and mixtures have been tabulated by Berger and Hubbell and given in the form of XCOM program at energies 1 keV to 100 GeV (Berger et al., 2010). The new version of this software, called WinXcom (Gerward et al., 2004) is nowadays used as user-friendly software to generate the desired data in Microsoft excel file in windows operating system.
FLUKA is a general purpose Monte Carlo simulation package for calculations of particle transport and interactions with matter. It is being used in various applications such as proton, electron accelerators shielding design, activity, dosimetry, detector design, cosmic rays, neutron physics, radiotherapy, accelerator driven system, etc. Also, it is useful in many scientific areas (high energy experimental physics and engineering, detector and telescope design, medical physics and radio-biology). FLUKA can simulate with high accuracy the interaction of more than 60 different types of particles such as heavy-ions, electrons, neutrons, photons, neutrinos, muons, and their antiparticles in many types of research fields and applications (Ferrari et al., 2005; Battistoni et al., 2007; Mark et al., 2007). The FLUKA can be used for transport of synchrotron radiation and optical photons too. FLUKA simulation code has been used in low- and intermediate- energy gamma-ray for investigation of radiation characteristics of soil (Wielopolski et al., 2005), shielding materials (Agosteo et al., 2005), X-, gamma-ray or radiation protection (Nariyama et al., 2003; Beskrovnaia et al., 2008), neutron shielding characteristics (Korkut et al., 2010), and water, concrete and bakelite (Demir et al., 2013). FLUKA code has been used in low- and intermediate-energies for gamma ray calculation of mass attenuation coefficients and good agreement is observed with experimental results (Demir et al., 2013). However, the FLUKA has not been tested for high-energy gamma ray with experimental results. This has encouraged us to utilize the FLUKA for gamma-ray interaction parameters for high-energies.
The aim of the present study is investigation of mass attenuation coefficients for high energy gamma-ray for ordinary, hematite-serpentine, ilmenite-limonite, basalt-magnetite, ilmenite, steel-scrap and steel-magnetite concretes using FLUKA code. First of all, the mass attenuation coefficients of the selected concretes for photon energy 1.5, 2, 3, 4, 5 and 6 MeV were calculated by using FLUKA code, and the linear attenuation coefficients were estimated. Finally the simulated results of linear attenuation coefficients were compared with the experimental data provided in the literature (Bashter, 1997). Good agreement among FLUKA code, XCOM data and experimental results for high energy gamma-ray was observed. The variation of mass attenuation coefficients determined using FLUKA code is shown graphically.
Different types of normal and heavy concretes have been taken in the literature of Bashter (1997), whose elemental compositions and densities are given in Table 1. These concretes are ordinary (OR), hematite-serpentine (HS), ilmenite-limonite (IL), basalt-magnetite (BM), ilmenite (IT), steel-scrap (SS) and steel-magnetite (SM) used in various applications of the shielding.
The Monte Carlo method is based on random numbers and mathematical algorithms (Ramirez-Lopez et al., 2011). It can be applied for physical systems, especially in nuclear science and engineering (Ferrari et al., 2005; Battistoni et al., 2007) as a Monte Carlo simulator. It is a Monte Carlo package used in interactions between all subatomic particles and matter. It has many advantages in terms of wide energy range. There are several studies using this Monte Carlo simulation code (Korkut et al., 2010, 2011, 2012; Ramirez-Lopez et al., 2011; Singh et al., 2015).
In the simulations, the latest version of FLUKA (2011.2b.4) was used. We have obtained I/I0 photon transmission values at 1.5, 2, 3, 4, 5 and 6 MeV photon energies by means of FLUKA code. The simulation has been done for all types of concretes. Linear attenuation coefficient is calculated using Lambert Beer Law (I/I0 = exp(-μx)). In this law, I0 is photon transmission, I is linear attenuation coefficient and x is the thickness of the sample. After the simulation process gamma transmission values have been read from FLUKA output file.1
The transmission of gamma-ray (I = I0 exp(-μt)) is dependent upon the thickness, t of the interacting medium and linear attenuation coefficient, μ. The μ is calculated by multiplication of mass attenuation coefficients, μ/ρ and density. The μ/ρ of the concretes are calculated by the mixture rule where wi is the proportion by weight and (μ/ρ)i is mass attenuation coefficient of the ith element by using XCOM or WinXcom. The linear attenuation coefficient of the concrete is multiplication of μ/ρ and the density of the concrete. The atomic number and atomic mass of elements have been taken from atomic weight of elements 2011, IUPAC (Michael et al., 2013). The uncertainties in μ/ρ values is about 1% for low-Z (1 < Z < 8) in Compton region (30 keV to 100 MeV). Below 30 keV energy, the uncertainties are as much as 5–10% because of correction to experiments for high-Z impurities and departure of Compton cross section from Klein-Nishina theory.
The linear attenuation coefficients, µ and mass attenuation coefficients, µ/ρ of the concretes have been investigated for high-energy (1.5, 2, 3, 4, 5 and 6 MeV) gamma-rays. Comparison of linear attenuation coefficients by using FLUKA code and experiment are provided in Table 2. The variation of mass attenuation coefficients for the ordinary concrete is shown in Figure 1.
Mass attenuation coefficients of ordinary concretes using FLUKA Monte Carlo simulation codes.
The linear attenuation coefficient, μ of the concretes by FLUKA code, XCOM and experiment are shown for photon energies 1.5, 2, 3, 4, 5 and 6 MeV in Table 2. It is observed that the linear attenuation coefficients simulated using FLUKA, XCOM and the experiment are in very good agreement. Therefore, it is concluded that the FLUKA is a useful simulation code for high energy gamma-rays interactions where data may not be available, analogous to the experiment.
The mass attenuation coefficients, µ/ρ using FLUKA code for ordinary concrete (as an example) with gamma-ray energy is shown in Figure 1. These µ/ρ values of the concretes decrease with the increase of gamma energy. The similar variation of the µ/ρ values with XCOM program can be found. The variation of the µ/ρ can be explained based on the partial interaction process Compton scattering and pair-production in the high energy gamma-rays. The experimental data of µ/ρ of the concretes at energies 1.5, 2, 3, 4, 5 and 6 MeV can be calculated using the linear attenuation coefficients. The mass attenuation coefficients are parameterized by polynomial fitting of order two.
The linear attenuation coefficients of various types of high photon energies (1.5, 2, 3, 4, 5 and 6 MeV) concretes using FLUKA simulation in the present investigation shows that FLUKA simulation is a very effective and capable tool for simulation of shielding materials at low, medium  as well as high energies.
The mass attenuation coefficients, μ/ρ and linear attenuation coefficients, μ of different types of normal and heavy concretes, ordinary, hematite-serpentine, ilmenite-limonite, basalt-magnetite, ilmenite, steel-scrap and steel-magnetite are compared for FLUKA simulation, experimental and XCOM theoretical data for energies 1.5, 2, 3, 4, 5 and 6 MeV. It was observed that the experimental and theoretical results are in very good agreement for the mass attenuation coefficients and linear attenuation coefficients. It is concluded that the FLUKA Monte Carlo simulation code is a useful tool for high energy gamma-ray interactions where data may not be available, analogous to the experiment.
- Agosteo S, Cammi A, Garlati L, Lombardi C, Padovani E. 2005. Gamma dose from activation of internal shields in IRIS reactor. Radiat. Prot. Dosim. 115: 86–91. [CrossRef] [Google Scholar]
- Akkurt I, El-Khayatt AM. 2013. Effective atomic number and electron density of marble concrete. J. Radioanal. Nucl. Chem. 295: 633–638. [CrossRef] [Google Scholar]
- Bashter II. 1997. Calculation of radiation attenuation coefficients for shielding concretes. Ann. Nucl. Energy 24: 1389–1401. [CrossRef] [Google Scholar]
- Bashter II, Makarious AS, El-Sayed Abdo AA. 1996. Investigation of hematite-serpentine and ilmenite-limonite concretes for reactor radiation shielding. Ann. Nucl. Energy 23: 65–71. [CrossRef] [Google Scholar]
- Battistoni G, Muraro S, Sala PR, Cerutti F, Ferrari A, Roesler S, Fasso A, Ranft J. 2007. The FLUKA code: Description and benchmarking. AIP Conf. Proc. 896: 31–49. [CrossRef] [Google Scholar]
- Berger MJ, Hubbell JH, Seltzer SM, Chang J, Coursey JS, Sukumar R, Zucker DS, K. Olsen K. 2010. XCOM: photon cross sections database, NIST standard reference database (XGAM). http://www.nist.gov/pml/data/xcom/index.cfm. [Google Scholar]
- Beskrovnaia L, Florko B, Paraipan M, Sobolevsky N, Timoshenko G. 2008. Verification of Monte Carlo transport codes FLUKA, GEANT4 and SHIELD for radiation protection purposes at relativistic heavy ion accelerators. Nucl. Instr. Meth. Phys. Res. B 266: 4058–4060. [CrossRef] [Google Scholar]
- Demir N, Akar UT, Popovici MA, Demirci ZN, Gurler O, Akkurt I. 2013. Investigation of mass attenuation coefficients of water, concrete and Bakelite at different energies using the FLUKA Monte Carlo code. J. Radioanal. Nucl. Chem. 298: 1303–1307. [CrossRef] [Google Scholar]
- Ferrari A, Sala PR, Fasso A, Ranft J. 2005. FLUKA: A multi-particle transportcode, CERN-2005-010, INFNTC_05/11, SLAC-R-773. [Google Scholar]
- Gerward L, Guilbert N, Jensen KB, Levring H. 2004. WinXcom-a program for calculating X-ray attenuation coefficients. J. Radiat. Phys. Chem. 71: 653–654. [CrossRef] [Google Scholar]
- Korkut T, Karabulut A, Budak G, Korkut H. 2010. Investigation of fast neutron shielding characteristics depending on boron percentages of MgB2, NaBH4 and KBH4. J. Radioanal. Nucl. Chem. 286: 61–65. [CrossRef] [Google Scholar]
- Korkut T, Korkut H, Karabulut A, Budak G. 2011. A new radiation shielding material: amethyst ore. Ann. Nucl. Energy 38: 56–59. [CrossRef] [Google Scholar]
- Korkut T, Karabulut A, Budak G, Aygun B, Genel O, Hancerliogullari A. 2012. Investigation of neutron shielding properties depending on number of boron atoms for colemanite, ulexite and tincal ores by experiments and FLUKA Monte Carlo simulations. Appl. Radiat. Isot. 70: 341–345. [CrossRef] [PubMed] [Google Scholar]
- Makarious S, Bashter II, Kany AM. 1988. Radiative capture gamma rays arising from iron fibre additions to ilmenite concrete shields. Ann. Nucl. Energy 15(10/11): 513–521. [CrossRef] [Google Scholar]
- Makarious S, Bashter II, El-Sayed AA, M. Samir AA, Kansouh WA. 1996. On the utilization of heavy concrete for radiation shielding. Ann. Nucl. Energy 23: 195–206. [CrossRef] [Google Scholar]
- Mark S, Khomchenko S, Shifrin M, Haviv Y, Schwartz JR, Orion I. 2007. TVF-NMCRCó A powerful program for writing and executing simulation inputs for the FLUKA Monte Carlo code system. Nucl. Instrum. Meth. Phys. Res. A 572: 929–934. [CrossRef] [Google Scholar]
- Michael EW et al. 2013. Atomic weight of elements 2011 (IUPAC Technical Report). Pure Appl. Chem. 85(5): 1047–1078. [CrossRef] [Google Scholar]
- Nariyama N, Konnai A, Ohnishi S, Odano N. 2003. Calculation of dosimeter response for in-human-phantom measurement to low-energy photons. In: Proceedings of the eleventh EGS4 users meeting in Japan, KEK proceedings 15: 53–58. [Google Scholar]
- Ramirez-Lopez A, Soto-Cortes G, Gonzalez-Trejo J, Munoz-Negron D. 2011. Computational algorithms for simulating the grain structure formed on steel billets using cellular automaton and chaos theories. Int. J. Miner. Metall. Mater. 18: 24–34. [CrossRef] [Google Scholar]
- Shirmardi SP, Shamsaei M, Naserpour M. 2013. Comparison of photon attenuation coefficients of various barite concretes and lead by MCNP code, XCOM and experimental data. Ann. Nucl. Energy 55: 288–291. [CrossRef] [Google Scholar]
- Singh VP, Badiger NM. 2012. Comprehensive study of energy absorption and exposure buildup factor for concrete shielding in photon energy range 0.015-15 MeV upto 40 mfp penetration depth: dependency of density, chemical element, photon energy. Int. J. Nucl. Eng. Sci. Tech. 7: 75–99. [Google Scholar]
- Singh VP, Badiger NM, Chanthima N, Kaewkhao J. 2014. Evaluation of gamma-ray exposure buildup factors and neutron shielding for bismuth borosilicate glasses. Radiat. Phys. Chem. 98: 14–21. [CrossRef] [Google Scholar]
- Singh VP, Shirmardi SP, Medhat ME, Badiger NM. 2015. Determination of mass attenuation coefficient for some polymers using Monte Carlo simulation. J. Vacuum 119: 284–288. [CrossRef] [Google Scholar]
- Wielopolski L, Song Z, Orion I, Hanson AL, G. Hendrey G. 2005. Basic considerations for Monte Carlo calculations in soil. Appl. Radiat. Isot. 62: 97–107. [CrossRef] [PubMed] [Google Scholar]
Detailed information can be seen in the fluka web page www.fluka.org.
Cite this article as: Singh VP, Korkut T, Badiger NM. 2018. Comparison of mass attenuation coefficients of concretes using FLUKA, XCOM and experiment results. Radioprotection 53(2): 145–148
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.