Research Article  Open Access
Calculation of Radioactivity and Dose Rate of Activated Corrosion Products in WaterCooled Fusion Reactor
Abstract
In watercooled reactor, the dominant radioactive source term under normal operation is activated corrosion products (ACPs), which have an important impact on reactor inspection and maintenance. A threenode transport model of ACPs was introduced into the new version of ACPs source term code CATE in this paper, which makes CATE capable of theoretically simulating the variation and the distribution of ACPs in a watercooled reactor and suitable for more operating conditions. For code testing, MIT PWR coolant chemistry loop was simulated, and the calculation results from CATE are close to the experimental results from MIT, which means CATE is available and credible on ACPs analysis of watercooled reactor. Then ACPs in the blanket cooling loop of watercooled fusion reactor ITER under construction were analyzed using CATE and the results showed that the major contributors are the shortlife nuclides, especially Mn56. At last a point kernel integration code ARShield was coupled with CATE, and the dose rate around ITER blanket cooling loop was calculated. Results showed that after shutting down the reactor only for 8 days, the dose rate decreased nearly one order of magnitude, which was caused by the rapid decay of the shortlife ACPs.
1. Introduction
In the cooling loop of watercooled reactor, the oxidation and corrosion of the metal material by water are inevitable, and part of the corrosion products will be activated under neutron irradiation and become radioactive, which are called activated corrosion products (ACPs). Some ACPs continuously decay and emit harmful gammarays even after shutdown of the reactor for a long time. According to the surveillance data from French PWR plants, more than 90% of integrated dose under normal operation is due to ACPs [1], which means vital impact of ACPs on radiation protection. For the watercooled fusion reactor, it is also necessary to research the formation mechanism and transport process of ACPs and to predict the variation and distribution of radioactivity and dose rate of ACPs, which is important to radiation shielding design, inspection, and maintenance as well as accident analysis of the reactor.
Because of the importance of ACPs, the relevant research was started since 1960s in PWR plants, and many codes have been developed [2], for example, CORA from EPRI, PACTOLE from CEA, CRUDTRAN from KAERI, and MIGART from BSA, which were widely used in ACPs analysis of PWR plants and showed good effect. For ACPs analysis of watercooled fusion reactor, especially the International Thermonuclear Experimental Reactor (ITER) under construction, three codes were developed, which are PACTITER from CEA, TRACT from UKAEA, and CATE from NCEPU (North China Electric Power University) by the author. In the initial version of CATE, the model adopted is mainly based on the empirical coefficient method [3], which is simple but not universal for different operation condition. In the latest version of CATE, a threenode transport model was adopted, which is based on the theory that the main driving force for corrosion products transport is the temperature change of the coolant throughout the loop and the resulting change in metal ion solubility in the coolant [4], making CATE capable of theoretically simulating the variation and the distribution of ACPs and have a larger scope of application.
The basic theory and equations of the threenode transport model were described in the second section of this paper, and ACPs in MIT PWR coolant chemistry loop and ITER blanket cooling loop were calculated using the code CATE, respectively, in the third section and the fourth section. Then the dose rate around ITER blanket cooling loop caused by ACPs was calculated using the point kernel integration code ARShield in the fifth section. In the last section, a comprehensive comment was presented.
2. Description of the ThreeNode Transport Model
In the model, the cooling loop is divided into three nodes, which are the influx pipe, the coolant, and the outflux pipe. The corrosion products transport between the nodes is by means of dissolution and deposition. Nickel ferrite is the main corrosion products in the cooling loop, and it shows an inverse solubility curve with temperature. As the coolant temperature changes around the cooling loop, the saturation status of the corrosion products in the coolant also changes. In the region of outflux (e.g., the heatexchanger), as the coolant temperature decreases, the soluble species become undersaturated, so that there is a driving force for the corrosion products on the pipe surface to dissolve into the coolant to restore a saturation concentration of corrosion products in the coolant. On the other hand, in the region of influx (e.g., the core in PWR or the blanket in fusion reactor), where the coolant temperature increases, the soluble species exist in a supersaturated state, so that there is a driving force for the soluble species to deposit on the pipe surface. Through this process, the corrosion products are transported from the outflux pipe to the influx pipe.
In the region of influx, part of the corrosion products absorb neutrons and become ACPs. The ratio of ACPs to corrosion products on the pipe surface is larger than that in the coolant, so there is a net flow of ACPs from the pipe surface to the coolant by isotope exchange. In the region of outflux, the ratio of ACPs to corrosion products in the coolant is larger than that on the pipe surface, so there is a net flow of ACPs from the coolant to the pipe surface. Through this process, ACPs are transported from the influx pipe to the outflux pipe, which is described in Figure 1.
In brief, the threenode transport model is based on the theory that the main driving force for ACPs transport is the temperature change of the coolant throughout the loop and the resulting change in metal ion solubility in the coolant. So whether it is PWR temperature range (280~320°C) or ITER temperature range (140~180°C), the threenode transport model is all applicable.
In CATE, the corresponding equations are as follows:where is the mass of nuclide in corrosion products, kg, the subscript refers to the zone of influx pipe, refers to the zone of bulk coolant, and refers to the zone of outflux pipe; is the corrosion rate of the base metal, kg/(m^{2}·s); is the ratio of release rate to corrosion rate of the corrosion products, %; is the mass fraction of nuclide in the base metal, %; is the contact area between coolant and pipe, m^{2}; is the mass transfer coefficient of nuclide , m/s; is the coolant volume in the cooling loop, m^{3}; is the solubility of nuclide in the coolant close to the pipe surface, kg/m^{3}; is the collection efficiency of filter and resin in chemical and volume control system (CVCS), %; is the volume flow rate of coolant into CVCS, m^{3}/s; is the coolant volume in CVCS, m^{3}; is the atom number of nuclide in ACPs, atoms, and the superscript means nuclide is radioactive; is the decay branching ratio from nuclide to nuclide , %; is the decay rate constant of nuclide , s^{−1}; is the radiative capture cross section of converting nuclide to nuclide , cm^{2}; is the neutron flux, /(cm^{2}·s); is the molar mass of nuclide , kg/mol; is the Avogadro number; is the absorption cross section of nuclide , cm^{2}; is the ratio of coolant residence time in the zone of influx pipe to that in the whole cooling loop, %.
In the above model, some assumptions are made, including the following: (1) the concentration of nuclides in the coolant keeps the same along the cooling loop, which is because it takes the coolant only several seconds to circle around the cooling loop, which means the homogenization effect of nuclides in the coolant is dominant; (2) the influence of activation and decay on the mass of corrosion products is neglected because the quantity of ACPs is very limited compared to that of corrosion products; (3) the concentration of nuclides on the pipe surface is equal to its solubility in the adjacent coolant.
To solve the above differential equations, the fourthorder RungeKutta method is used, which is fast and can control the error well. Moreover, an algorithm of adaptive time step is adopted in CATE for determining the time step reasonably and the calculation speed can be improved.
3. Verification of CATE Code through MIT PWR Coolant Chemistry Loop (PCCL)
3.1. Description of MIT PCCL
MIT PCCL was chosen here to test the code CATE, which is a small scale loop operated under constant coolant chemistry in the MIT reactor to closely simulate the primary circuit of a typical PWR power plant. The main operation data of the loop is presented in Table 1 [5].

The base metal of influx pipe is Zircaloy4 and the corrosion rate of this material is nearly 0. The base metal of outflux pipe is mainly Inconel and the corrosion rate of it is .
3.2. Results of ACPs Radioactivity in MIT PCCL
The surveillance data from MIT PCCL for 42 days of full power operation was published, and the radioactivity of ACPs in the loop is as shown in Table 2.

Calculation results from CATE are close to the experimental results from MIT, which means CATE is available and credible on ACPs analysis of watercooled nuclear reactor. The calculation results of mass of corrosion products in the influx pipe, the coolant, and the outflux pipe are, respectively, kg, kg, and kg for CATE. Using these data, we can calculate the specific radioactivity of ACPs in the influx pipe, the coolant, and the outflux pipe, which is, respectively, GBq/kg, GBq/kg, and GBq/kg for CATE. The values are degressive, which is consistent with the transport direction of ACPs in the cooling loop and is rational.
4. Calculation of ACPs Radioactivity in ITER Blanket Cooling Loop Using CATE Code
4.1. Description of ITER Blanket Cooling Loop
The International Thermonuclear Experimental Reactor (ITER) is under construction now. In its design plan, the primary cooling loops are most watercooled, such as the blanket, the divertor, the NB injector, and the vacuum vessel. So ITER can be treated as a representative of watercooled fusion reactor. China is developing its own fusion reactor CFETR (China Fusion Engineering Test Reactor), whose design plan of heat transfer and radiation shielding partly refers to ITER, so ACPs analysis of ITER will benefit the work of source term control for CFETR.
The planned operation data of ITER blanket cooling loop is presented in Table 3 [6].

The base metal in influx pipe and outflux pipe is the same kind of stainless steel, SS316, and the corrosion rate of it is as follows:
The above equation and the element composition of SS316 are quoted from [7]. The relevant nuclear reaction data are quoted from the European Activation File EAF2007 [8, 9], and the cross section of activation reaction is collapsed from 172 groups into 1 group using the corresponding neutron spectrum, which is calculated with MCNP code [10] and based on an ITER blanket module [11]. The neutronics model and neutron spectrum of the chosen blanket module can be seen in Figures 2 and 3.
4.2. Results of ACPs Radioactivity in ITER Blanket Cooling Loop
After 1.2 years of full power operation, the mass of corrosion products and the radioactivity of ACPs in the loop are as shown in Table 4.

From Table 4, we can see that mass of corrosion products in the coolant is much lower than that on the pipe surface, which is due to the limitation of solubility of corrosion products in the coolant. It should be noticed that the corrosion products on the pipe surface include oxides and deposits. The oxides are compact and fixed, while the deposits are loose and mobilizable and their mass is usually one order of magnitude lower than that of oxides. So, although the total mass of corrosion products on the pipe surface is tens of kg, the actual mass of deposits is only in the magnitude of kg, which is consistent with the experience in safety analysis.
We also can see that ACPs radioactivity on the pipe surface is much higher than that in the coolant, which means the pipe surface is the main radioactive contamination region and should be decontaminated regularly through water chemistry method. The specific radioactivity of ACPs in the influx pipe, the coolant, and the outflux pipe is calculated as GBq/kg, GBq/kg, and GBq/kg. The values are degressive, which is consistent with the transport direction of ACPs in the cooling loop and is rational.
From Table 5, we can see that the shortlife nuclides (V52, V53, Cr55, Mn56, Co58m, Co60m, and Ni57) are the major contributors to radioactivity in all the three regions, especially Mn56, which alone contributes 63.33% of the radioactivity in the region of influx pipe and nearly 80% of the radioactivity in the regions of coolant and outflux pipe. That is very different from PWR, in which the longlife nuclides Co58 and Co60 are the main ACPs. And when the reactor shuts down for a period of time (about several days), ACPs radioactivity will decrease obviously due to the rapid decay of the shortlife nuclides, and then the right time for workers to do maintenance comes.

5. Calculation of Dose Rate Caused by ACPs Using ARShield Code
The ARShield code developed by NCEPU of China is applied to dose rate calculation, which is a new version of the point kernel integration code QADCG developed by Los Alamos National Laboratory. ARShield breaks some restrictions of QADCG, such as complicated modeling, complicated source setting, 3D fine mesh results statistics, and largescale computing efficiency, and is proved to be reliable and efficient on dose rate calculation.
The density of each radionuclide at chosen regions calculated by CATE is introduced into ARShield and then converted to dose rate using point kernel integration method, which is as follows:where is the point at which gamma dose rate is to be calculated; is the location of source in volume ; is the volume of source region; is the total attenuation coefficient at energy ; is the distance between source point and point at which gamma intensity is to be calculated; is the fluxtodose conversion factor; is the dose buildup factor.
The geometry of hot leg pipe in ITER blanket cooling loop is adopted to represent the region of outflux pipe, which has an internal diameter of 0.527 m, thickness of 0.02 m, and length of 33 m. The dose rate around the hot leg pipe calculated by ARShield is shown in Figure 4.
We can see that after shutting down the reactor for 8 days, the dose rate is much lower than that of normal operation for 1.2 years, which is because that contribution from the shortlife ACPs decreases almost to zero after shutdown for 8 days. The typical values of dose rate are as follows: during normal operation, the dose rate at the outer surface of the pipe is 8.80 mSv/h, and the dose rate at a distance 1 m away from the outer surface of the pipe is 2.52 mSv/h; after shutdown for 8 days, the values drop to 1.05 mSv/h (contact) and 0.31 mSv/h (1 m away). When these values are compared to the annual permissible worker dose rate of 20 mSv/year, recommended by the ICRP [12], the contact dose rate value after shutdown for 8 days would allow approximately 19 h exposure per year.
6. Conclusions
In this paper, a threenode transport model was introduced into the ACPs source term code CATE, making CATE capable of theoretically simulating the variation and the distribution of ACPs in a watercooled reactor and suitable for more operating conditions. MIT PCCL was chosen to test the new version of CATE, and the calculation results from CATE are close to the experimental results from MIT, which means CATE is available and credible on ACPs analysis of watercooled reactor. Then the radioactivity and composition of ACPs in ITER blanket cooling loop were analyzed using CATE and the results showed that the major contributors are the shortlife nuclides (V52, V53, Cr55, Mn56, Co58m, Co60m, and Ni57) for ITER, especially Mn56. That is very different from PWR, in which the longlife nuclides Co58 and Co60 are the main ACPs. At last, the dose rate around ITER blanket cooling loop caused by ACPs was calculated through coupling the code CATE with a point kernel integration code ARShield. The results showed that after shutting down the reactor only for 8 days, the dose rate can decrease nearly one order of magnitude compared to that of normal operation, which is caused by the rapid decay of the shortlife ACPs. In the future, CATE will be applied for ACPs analysis of Chinese watercooled fusion reactor CFETR as soon as its design plan is completed.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was sponsored jointly by “National Special Project for Magnetic Confined Nuclear Fusion Energy” with Grant no. 2014GB119000 and the “Fundamental Research Funds for the Central Universities” with Grant no. 2014QN26 in China. The authors would like to express their gratitude for the support.
References
 A. Rocher, J. L. Bretelle, and M. Berger, “Impact of main radiological pollutants on contamination risks (ALARA) optimization of physico chemical environment and retention technics during operation and shutdown,” in Proceedings of the European Workshop on Occupational Exposure Management at NPPs (ISOE '04), Session 2, EDF, Lyon, France, March 2004. View at: Google Scholar
 M. Rafique, N. M. Mirza, S. M. Mirza, and M. J. Iqbal, “Review of computer codes for modeling corrosion product transport and activity buildup in light water reactors,” Nukleonika, vol. 55, no. 3, pp. 263–269, 2010. View at: Google Scholar
 L. Li, J. Zhang, W. Song, Y. Fu, X. Xu, and Y. Chen, “CATE: a code for activated corrosion products evaluation of watercooled fusion reactor,” Fusion Engineering and Design, vol. 100, pp. 340–344, 2015. View at: Publisher Site  Google Scholar
 IAEA, “Modelling of transport of radioactive substances in the primary circuit of watercooled reactors,” IAEATECDOC1672, 2012. View at: Google Scholar
 C. B. Lee, Modeling of Corrosion Product Transport in PWR Primary Coolant, Nuclear Engineering Department, MIT, 1990.
 S. Nisan, L. D. Pace, and J.C. Robin, “Evaluation of the activated corrosion products for the ITER heat transfer systems,” in Fusion Technology, pp. 1763–1766, Elsevier, 1996. View at: Google Scholar
 P. J. Karditsas, “Activation product transport using TRACT: ORE estimation of an ITER cooling loop,” Fusion Engineering and Design, vol. 45, no. 2, pp. 169–185, 1999. View at: Publisher Site  Google Scholar
 R. A. Forrest, J. Kopecky, and J.Ch. Sublet, “The European Activation File: EAF2007 neutroninduced cross section library,” UKAEA FUS 535, 2007. View at: Google Scholar
 R. A. Forrest, “The European Activation File: EAF2007 decay data library,” UKAEA FUS 537, 2007. View at: Google Scholar
 J. F. Briesmeister, MCNP: a general Monte Carlo Nparticle transport code version 4C, Los Alamos National Laboratory, LA13709M, 2000.
 H. Iida, V. Khripunov, and L. Petrizzi, “Nuclear analysis report,” ITER Document NAG201010617FDR, 2001. View at: Google Scholar
 ICRP, “Recommendations of the international commission on radiological protection, publication 60,” Annals of the ICRP 21, Pergamon Press, Oxford, UK, 1991. View at: Google Scholar
Copyright
Copyright © 2016 Jingyu Zhang et al. 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.