Research Article  Open Access
Calculation of LeadIron DoubleLayer Thickness for GammaRay Shielding by MATLAB Program
Abstract
This contribution is aimed at designing the optimal thickness of leadiron doublelayer container to store a radioactive waste releasing the photon energy at 1.3325 MeV and initial radiation intensity at 100 mSv/hr using the optimization design by MATLAB software. This design consisted of three parts of calculations to achieve 1000 times the radiation attenuation of container. The first was the logarithmic interpolation for the mass attenuation coefficient. The second was the bilogarithmic interpolation for the exposure buildup factor. The third was the contourplotting analytical technique for the optimal thickness of radiation container. The values of mass attenuation coefficient and exposure buildup factor were exactly validated as compared with the standard reference database. Furthermore, we have found that the optimal thickness was 3.2 cm for lead (1st layer) and 17.0 cm for iron (2nd layer). Container weight was 994.30 kg, whilst container cost was 167.30 USD. The benefit of our design can quickly and precisely apply for the radiation safety assessment of the occupational radiation workers who always work in the nuclear reactor area.
1. Introduction
Although nuclear technology is very useful for researches and industries, maintenance of nuclear reactor is necessary for the safety reason. International Commission on Radiological Protection (ICRP) has suggested that the limitation of radiation exposure dose for occupational radiation worker is 20 mSv/year [1]. During the maintenance period, Thailand Institute of Nuclear Technology (TINT) has informed that the gamma radiation intensity of a bolt taken from the nuclear reactor was about 100 mSv/hr, whereas photon energy was 1.3325 MeV that was similar to Co60 source. Consequently, this is considered as a radioactive waste. To evaluate the radiation safety due to the nuclear interactions, one of the best ways is the use of Monte Carlo simulation [2–4]. According to Monte Carlo approach, statistical processing time is prolonged at 13 hrs partly due to the complex geometry design [5, 6]. To avoid wasting of time, MATLAB software is certainly introduced to calculate the thickness of materials. This thickness will be used as an input variable in Monte Carlo simulation to reduce wasting of time. There are four advantages for MATLAB. First, this software is able to plot the contour graph to analyze the optimal thickness of materials. This contourplotting technique is easy and more suitable than 3D surface plotting technique for finding the optimum solution [7]. Second, there are a variety of commands used in forloop and ifclause processes [8] which are important components for the optimization model. Third, output data stored in matrix format after mathematical processing is very simple for the next mathematical process. Last, the command patterns in calculation for processing data are easily written.
Here we have performed the calculations for shielding materials of lead and iron layers for a 1.3325 MeV and 100 mSv/hr gammaray by MATLAB software. This design is comprised of three parts. The first and second parts are the logarithmic interpolation for mass attenuation coefficient and the bilogarithmic interpolation for exposure buildup factor, respectively. These parts are compared with the standard reference database from National Institute of Standards and Technology (NIST) [9] and American National Standard (ANS6.4.3) for checking a validation [10]. The third part is to apply this optimization model to determine the suitable doublelayer thickness (leadiron) for attenuating radiation intensity of 0.1 mSv/hr via the contourplotting analytical technique.
2. Simulation Procedures
2.1. Mathematical Format for Optimization
Objective function was defined as the gamma attenuation of doublelayer materials in a narrow (Figure 1(a)) and board beam (Figure 1(b)) [11].
(a)
(b)
To achieve 1000 times the radiation attenuation of container as informed by TINT and to follow the radiation dose limitation as informed by ICRP, our objective function was limited at 0.1 mSv/hr orwhere and were the initial and transmitted radiation intensity (mSv/hr), respectively. and were the linear attenuation coefficient (cm^{−1}) of lead and iron, respectively. and were the exposure buildup factor of lead and iron, respectively. Mass attenuation coefficient was calculated by (2) according to logarithmic interpolation (LI) as shown in Figure 2(a) [12, 13]. Exposure buildup factor was computed by (3)–(5) according to bilogarithmic interpolation (BI) modified from bilinear interpolation as seen in Figure 2(b) [14, 15]. where Constrained functions were divided into three parts: material thickness, weight, and cost. 3D and 2D configuration of containers were illustrated in Figures 3 and 4, respectively.
(a)
(b)
(a)
(b)
For economical reason, constrained function of the doublelayer thickness of cylindrical container was limited at 30 cm orFor technical reason, constrained function of the container weight was limited at 1,000 kg orwhere were the density of air, lead, and iron (g/cm^{3}), respectively. was the radius of the gap inside a container.
For economical reason, constrained function of the container’s cost was limited at 800 USD orwhere , were the cost of lead at 0.97 USD/lb as of April 21, 2017 [16] and iron at 0.028 USD/lb as of January 31, 2017 [17], respectively.
2.2. Contour Plot Technique
Photon energy (E = 1.3325 MeV), initial radiation intensity ( = 100 mSv/hr), material #1 (Pb), material #2 (Fe), material #1 thickness (21.0 cm), and material #2 thickness (2.0 cm) were imported as the input variables into MATLAB. Hence, linear attenuation coefficient, penetration depth, and exposure buildup factor were calculated as shown in Figure 5 in order to obtain the contour plot.
In detail, the mass attenuation coefficient from (2) was calculated as shown in Figure 6. After getting the penetration depth from (6), the exposure buildup factor from (3)–(5) was preceded as shown in Figure 7.
3. Results and Discussion
3.1. Mass Attenuation Coefficient and Exposure Buildup Factor
To investigate a validation of MATLAB, we firstly check according to condition 1 (red frame from Figure 6). For instance, when is 0.15 MeV, we get = 0.143600 cm^{2}/g from MATLAB that is the same number as informed by NIST [9]. On the other hand, we check according to condition 1 (red frame from Figure 7). For example, when we have = 5 mfp and = 20 MeV, we get = 5.22 from MATLAB that is the same number as ANS6.4.3 [10]. This number also gives better results than invariant embedding method (IEM) [18].
In our study, if we have = 1.3325 MeV according to C60 and = 100 mSv/hr, the mass attenuation coefficient and exposure buildup factor are calculated in Table 1 and they are summarized in Figure 8.

(a)
(b)
3.2. Optimal DoubleLayer Thickness
To understand the effects of container thickness on the container weight and cost, the space or void inside the container calculated by (7) is varied from 0 to 30 cm. The results are shown in Figure 9. It is observed that container’s weight and cost calculated by (8)(9) tend to increase if we increase the gap inside the container.
(a)
(b)
However, when we assume = 7.0 cm that is possibly suitable for keeping the radioactive waste, the optimal thickness analysis on the contour graph is presented in Figure 10. This optimal solution is identified as an intersection between the objective function and the all constrained functions [7]. Coordinate points on this identified line (red arrow in Figure 10) are therefore taken to calculate the weight and cost of container as shown in Table 2.

Finally, the optimal doublelayer thickness selected from the identified line is = 3.2 cm (lead) and = 17.0 cm (iron). Thus, container’s weight and cost are 994.30 kg and 167.30 USD, respectively.
To sum up, this optimization model concerns two merits. One is that the users are able to modify whatever the parameters (i.e., photon energy, initial intensity, material types and thicknesses, and material weight and cost) are to obtain the optimal thickness of materials of interests. Another one is that the simulation time is very fast just 10 sec.
4. Conclusions
From optimization model of the doublelayer shielding design and selection of lead and iron cylindrical container by MATLAB software to store the radioactive waste at 1.3325 MeV and 100 mSv/hr, we have found that the mass attenuation coefficient and exposure buildup factor were 0.056601 cm^{2}/g for lead, 0.051862 cm^{2}/g for iron and 4.7316 for lead, 1.6681 for iron, respectively. These numbers were the same as the standard reference database. The doublelayer thickness selected from the analysis on contourplotting graph was 3.2 cm for lead and 17.0 cm for iron to achieve 1000 times the radiation attenuation of container (0.1 mSv/hr). The total container weight and cost from these designed materials were 994.30 kg and 167.30 USD, respectively.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This work was supported by King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand. The authors also thank Kittiphot Songkaitiwong for helpful discussions.
References
 HPS Specialists in Radiation Protection, 2017. https://hps.org/publicinformation/ate/q8900.html.
 N. V. Hai and N. V. Hung, “Using weight window in MCNPX to calculate and simulate dose distribution at outside of treatment room in Dong Nai General Hospital to radiation safety assessment for this area,” International Journal of Science and Research, vol. 3, no. 12, pp. 445–448, 2014. View at: Google Scholar
 V. P. Singh, A. M. Ali, N. M. Badiger, and A. M. ElKhayatt, “Monte Carlo simulation of gamma ray shielding parameters of concretes,” Nuclear Engineering and Design, vol. 265, pp. 1071–1077, 2013. View at: Publisher Site  Google Scholar
 R. Merk, H. Kröger, L. EdelhäuserHornung, and B. Hoffmann, “PENELOPE2008 Monte Carlo simulation of gamma exposure induced by 60Co and NORMradionuclides in closed geometries,” Applied Radiation and Isotopes, vol. 82, pp. 20–27, 2013. View at: Publisher Site  Google Scholar
 Los Alamos Laboratory, MCNP – A General Monte Carlo NParticle Transport Code (Version 5) Volume 1: Overview and Theory, USA, 2003.
 Los Alamos Laboratory, MCNP – A General Monte Carlo NParticle Transport Code (Version 5) Volume 2: User’s Guide, USA, 2003.
 O. Venkataraman, Applied Optimization with MATLAB Programming, Wiley, USA, 2009.
 S. J. Chapman, MATLAB Programming with Applications for Engineers, Cengage Learning, USA, 2012.
 Physical Measurement Laboratory, National Institute of Standards and Technology, Xray Mass Attenuation Coefficients, USA, 1996.
 Oak Ridge National Laboratory, American National Standard, New GammaRay Buildup Factor Data, USA, 1988.
 M. AlArif and D. Kakil, “CalculatedExperimental Model for Multilayer Shield,” Aro, The Scientific Journal of Koya University, vol. 3, no. 1, pp. 23–27, 2015. View at: Publisher Site  Google Scholar
 G. S. Sidhu, P. S. Singh, and G. S. Mudahar, “Energy absorption buildup factor studies in biological samples,” Radiation Protection Dosimetry, vol. 86, no. 3, pp. 207–216, 1999. View at: Publisher Site  Google Scholar
 G. S. Sidhu, P. S. Singh, and G. S. Mudahar, “A study of energy and effective atomic number dependence of the exposure buildup factors in biological samples,” Journal of Radiological Protection, vol. 20, no. 1, pp. 53–68, 2000. View at: Publisher Site  Google Scholar
 K.T. Chang, Computation for Bilinear Interpolation: Introduction to Geographic Information Systems, McGrawHill, New York, NY, USA, 7th edition, 2014.
 Wikipedia, 2017, Bilinear Interpolation. https://en.wikipedia.org/wiki/Bilinear_interpolation.
 InverstmentMine, 2017. Lead Prices and Lead Price Charts. http://www.infomine.com/investment/metalprices/lead/.
 InverstmentMine, 2017. Iron Prices and News. http://www.infomine.com/investment/ironore/.
 G. S. Sidhu, P. S. Singh, and G. S. Mudahar, “Estimation of exposure buildup factor in iron using different method: A comparative study,” Journal of Nuclear Energy Science & Power Generation Technology, vol. 3, no. 2, 6 pages, 2014. View at: Google Scholar
Copyright
Copyright © 2017 Sarai Lekchaum and Kitsakorn Locharoenrat. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.