Development and validation of the code COUPLE3.0 for the coupled analysis of neutron transport and burnup in ADS

Development and validation of the code COUPLE3.0 for the coupled analysis of neutron transport and burnup in ADS 1 2 first-author 1 corresp Corresponding author. E-mail address: yangyongwei@impcas.ac.cn Corresponding author. E-mail address: yangyongwei@impcas.ac.cn yangyongwei@impcas.ac.cn 3 4 1 1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China Institute of Applied Physics and Computed Mathematics, Beijing 100088, China Institute of Applied Physics and Computed Mathematics, Beijing 100088, China CAEP Software Center for High Performance Numerical Simulation, Beijing 100088, China CAEP Software Center for High Performance Numerical Simulation, Beijing 100088, China The analysis of the fuel depletion behavior is critical for maintaining the safety of accelerator-driven subcritical systems (ADSs). The code COUPLE2.0 coupling 3-D neutron transport and point burnup calculations was developed by Tsinghua University. A Monte Carlo method is used for the neutron transport analysis and the burnup calculation is based on a deterministic method. The code can be used for the analysis of targets coupled with a reactor in ADSs. In response to additional ADS analysis requirements at the Institute of Modern Physics at the Chinese Academy of Sciences, the COUPLE3.0 version was developed to include the new functions of 1. a module for the calculation of proton irradiation for the analysis of cumulative behavior using the residual radionuclide operating history, 2. a fixed-flux radiation module for hazard assessment and analysis of the burnable poison, and 3. a module for multi-kernel parallel calculation, which improves the radionuclide replacement for the burnup analysis to balance the precision level and computational efficiency of the program. This paper introduces the validation of the COUPLE3.0 code using a fast reactor benchmark and ADS benchmark calculations. Moreover, the proton irradiation module was verified by a comparison with the analytic method of calculating the 210 Po accumulation results. The results demonstrate that COUPLE3.0 is suitable for the analysis of neutron transport and the burnup of nuclides for ADSs. COUPLE3.0 Neutron transport Burnup Accelerator-driven subcritical system 1. Introduction 1 An accelerator-driven subcritical system (ADS) is considered suitable for the transmutation of minor actinides (MA) and it also achieves energy production with good safety characteristics due to the sub-criticality [ 1 1 , 2 2 ] . Therefore, numerous research studies on ADSs have been conducted worldwide. Fuel depletion analysis is an important part of the evaluation of a subcritical reactor core. Domestic and foreign research institutions have developed several different programs [ 3 3 - 5 5 3-5 ] for the analysis of ADSs. Verified by benchmark calculations, these programs can be used to analyze a subcritical reactor of ADSs. During the continuous operation of ADSs, large amounts of spallation products are produced in the target zone by protons and neutrons. Most of them are toxic with high levels of radioactivity and the level of the external exposure dose rate is one of the important issues of radiation protection during operation and maintenance. For the lead-bismuth-eutectic (LBE) target of the ADS, the 210 Po is of high concern [ 6 6 ] ; it is produced by a (n,γ) reaction of 209 Bi together with the β-decay of 210 Bi. When tungsten is used as the spallation target material, we should be concerned about Hf because it has a large cross-section for neutron absorption [ 7 7 ] , which has serious impacts on the target-reactor coupling behavior. Therefore, it is important to consider the accumulation behavior of the spallation products in the target. The code COUPLE2.0 [ 8 8 - 10 10 8-10 ] , developed by the Institute of Nuclear and New Energy Technology (INET) at Tsinghua University, couples the neutron transport program, MCNPX, based on the Monte Carlo method with the nuclide depletion calculation code, ORIGEN2.1. It can be used for the analysis of a subcritical reactor in an ADS [ 11 11 ] . In response to additional ADS research requirements at the Institute of Modern Physics at the Chinese Academy of Sciences, the COUPLE3.0 version was developed in collaboration with the High Performance Computation Center at the China Academy of Engineering Physics. For the COUPLE3.0 code, the following functions were added: a). A proton radiation calculation module was developed that can be used to analyze the accumulation behavior of residual nuclides in the spallation target based on the running history. b). A fixed-flux neutron irradiation module was added for the analysis of burnable poisons. c). According to the depletion of some important nuclides and considering the contribution of the reaction of some minor actinides (n, absorption) to the overall neutron absorption in the reactor, the user can set an Importance Cut-off Value (ICV) to balance the precision level and the computation time. The transport process studies were achieved with a parallel message-passing interface (MPI) to improve the calculation efficiency. The addition of the post-processing function of the nuclide information is convenient for the extraction of the nuclide information. COUPLE3.0 consists of two modules, the neutron irradiation module and the proton irradiation module. In the first module, a constant fission power is used for the burnup analysis of the fuel zone; a constant flux model is also used for the burnup analysis for the burnable poison zones that do not produce fission power. The second module allows for determining the effects of the bare proton target and the target-reactor coupled system on the accumulation behavior of the residual nuclides. Using a thorium-based ADS benchmark, we verified the COUPLE3.0 code with regard to the burnup behavior in the ADS system that was induced by a fixed neutron source. A validation was conducted for the applicability of the code for the MA transmutation analysis using a fast reactor benchmark. In order to verify the accuracy of the proton irradiation module, the accumulation of 210 Po in the target was calculated using COUPLE3.0 and the results were compared with the estimated values from an analytical method; the results were in good agreement. 2. Neutron irradiation module 1 The neutron irradiation burnup calculation is one of the functions of the COUPLE3.0 code, which is based on the COUPLE2.0 code; in the new version, many aspects of the physical model were improved to enhance the accuracy of the expansion method for solving the equation of the burnup chains. Because of the existence of a spatial difference in the neutron flux distribution, the nuclear fuel depletion is different for each location in actual systems. Generally, an approximation can be obtained by dividing the fuel depletion area. The internal neutron flux in the fuel depletion area does not change significantly and there is no significant difference in the nucleon density in this region; therefore, the region can be simulated as a point. In a relatively short period of time, the fuel depletion is not very obvious, the change in the neutron spectrum is not very significant, and the change in a single-group cross-section is not obvious; therefore, the solution of the linear depletion equation can be obtained. In COUPLE3.0, MCNPX2.7.0 [ 12 12 ] is used for the transport calculation and the single-group section for the fuel burnup calculation based on the neutron flux and ORIGEN2.1 [ 13 13 ] is used for the burnup calculation to obtain the nuclide compositions at the end of the time step. The main part of COUPLE3.0 is responsible for delivering the necessary data for the neutron transport calculation by transferring the merged single-group cross-sections from MCNPX to ORIGEN2.1, calculating the power level of the depletion areas, and transferring them to ORIGEN2.1; at the end of the time step of the burnup calculation, ORIGEN2.1 transfers the material compositions to MCNPX to perform the neutron transport calculation in the next step. The two programs run alternately until the fuel depletion time reaches the end. The process of the neutron irradiation calculation is shown in Fig. 1. Fig. 1. COUPLE3.0 can run in a parallel MPI environment. The data used for MCNPX come from the ENDF/B-VII.0 library, which is currently suitable for fuel depletion calculations. However, it should be pointed out that, because of the large dataset of ENDF/B-VII.0 [ 14 14 ] , the calculations are computationally intensive. In COUPLE3.0, all nuclides are assumed to be able to change over time, which allows the program to deal with practical problems, such as fuel containing neutron absorption or the depletion of poisonous materials. In order to improve the calculation accuracy, a branching ratio is introduced for the excited states of the minor actinides. 3. Proton radiation calculation module 1 A proton irradiation function was added to COUPLE3.0 and can simulate the high-energy proton bombardment on the spallation target, the spallation reaction of the residual nuclides with depletion, and the accumulation of heavy metals. The spallation reaction generates neutrons and large amounts of residual nuclides, which accumulate and disappear by decay or other types of nuclear reactions such as absorption. The above process is complex and for an accurate simulation, the cross-section of the proton reaction and the residual nuclear yield of the different nuclides have to be known. Using these parameters, an equation of the proton irradiation burnup can be written similarly to the neutron irradiation equation; however, because some data are not available, it is impractical to develop an equation of the proton irradiation burnup. COUPLE3.0 coupled with MCNPX2.7.0 makes the simulation possible, not because of the proton reaction data of MCNPX2.7.0, but because COUPLE3.0 can calculate the residual nuclear yield after the different incident particles have bombarded the nuclides. The meaning of the yield of the residual nuclides is the same as that of the recharge rate of the F i in the burnup equation (1) used in ORIGEN2.1, and it can be used for the simulation of the burnup process of the proton irradiation. <math display="block" id="M1"><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mtext>d</mtext><msub><mi>X</mi><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><msub><mi>l</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>λ</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub></mrow></mstyle><mo>+</mo><mi>ϕ</mi><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><msub><mi>f</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>σ</mi><mi>k</mi></msub><msub><mi>X</mi><mi>k</mi></msub></mrow></mstyle><mo>−</mo><mo stretchy="false">(</mo><msub><mi>λ</mi><mi>i</mi></msub><mo>+</mo><mi>ϕ</mi><msub><mi>σ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>r</mi><mi>i</mi></msub><mo stretchy="false">)</mo><msub><mi>X</mi><mi>i</mi></msub><mo>+</mo><msub><mi>F</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mi>i</mi><mo>=</mo><mn>1...</mn><mi>N</mi></mtd></mtr></mtable></math> where X i is the nucleon density of the I th nuclides. N is the upper limit of the type of nuclide. L ij is the j th nuclide acting as the parent of the I th nucleus, which produces a fraction of the I th nuclide after decay. λ j is the decay constant of the j th species. ϕ is the neutron flux averaged over the energy zone and the burnup area. f ik is the share of the I th nuclide resulting from the nuclear reaction of the k th nuclide. σ k is the single-group reaction cross-section of the k th nuclide averaged over the neutron energy spectrum. r i is the removal ratio of the i th nuclide. F i is the recharge rate of the i th nuclide. When COUPLE3.0 is used in the calculation of proton irradiation, it is similar to the neutron irradiation calculation method and the concept of a transport step is also introduced. In the transport step, the accumulation of the spallation product is not obvious and it has little effect on the neutron yield, the neutron energy spectrum, and the residual nuclear recharge rate. Thus, the burnup calculation is performed with a constant recharge rate. In the next step, it is considered that the accumulation of the spallation products is considerable and that the influence on the transport cannot be ignored. At the same time, the materials in the burnup zone are renewed and the recharge rate is calculated for the burnup calculation. MCNPX2.7.0 is used for the transport calculation to obtain the distribution of the spallation neutron flux, the single-group cross-section, and the remaining nuclear recharge rate. ORIGEN2.1 is used for the burnup calculation in order to obtain the nuclide content of the burnup zones at the end of the transport step. The main portion of the code is responsible for transferring data between MCNPX2.7.0 and ORIGEN2.1 to achieve the coupling calculation of transport and burnup. It calculates the neutron flux and the recharge rate at a given proton irradiation level that are transferred to ORIGEN2.1 for the burnup calculation. Subsequently, COUPLE3.0 extracts the calculation results of ORIGEN2.1 at the end of the time step, selects the important nuclides, and updates the fuel compositions of the materials in the burnup zones for the transport calculation in the next step. It alternately and repeatedly calls the MCNPX2.7.0 and ORIGEN2.1 codes until the end of the proton irradiation cycle is reached. The flowchart for the proton irradiation calculation in COUPLE3.0 is shown in Fig. 2. Fig. 2.   The proton irradiation burnup calculation in COUPLE3.0 is divided into two modes. The first one considers the transmutation of the residual nuclides and the second one only considers the decay of the residual nuclides. The first method is closer to the actual physical process. Because the neutron flux in the spallation neutron target area is high, the transmutation reactions of the residual nuclides should not be neglected; in addition, the neutrons from the subcritical reactor core result in the transmutation of the residual nuclides.  It is necessary to point out that, due to the lack of a cross-section in high-energy neutron data, the reaction rates are underestimated. For the second kind of treatment, because of the neglect of the disappearance of the residual nuclides due to the nuclear reaction induced by the neutrons, the accumulation of the residual nuclides is overestimated. In view of the safety analysis, the effect of the specific nuclides is often overestimated. Therefore, this method is conservative. For a liquid or a flowing target, the residual nucleus will spread to other areas due to the flow and the actual cumulative amount of the residual nuclei in the target area is less than the amount in the stationary target. To solve this problem, COUPLE3.0 uses a rough approximation and considers that the migration rate of the residual nucleus is due to the flow of the fluid and that each spallation zone corresponds to a migration rate; this is represented by a coefficient between 0 and 1. At first, COUPLE3.0 calculates the recharge rate of the remaining nuclides without flow by calling MCNPX2.7.0; subsequently, the recharge rate is multiplied by the rate of fuel consumption and then transferred to ORIGEN2.1 for the burnup calculation. The migration rate needs to be manually set by the user. When the migration rate is 1, all of the nuclides are lost. When the migration rate is 0, there is no loss of nuclides. 4. Preliminary verification of COUPLE3.0 1 4.1 Fast reactor benchmark verification 2 In order to verify the applicability and reliability of COUPLE3.0 when used in the calculation of a reactor with a fast neutron spectrum, we performed the calculation and analysis of the benchmark of the fast spectrum plutonium-burning reactor proposed by the Organization for Economic Cooperation and Development Nuclear Energy Agency (OECD/NEA). Fig. 3 Fig. 3 shows the geometric model of the benchmark problem. In the calculation, the overall burnup time is set as 625 effective full power days (EFPD) and the thermal power is set as 1200 MW. Table 1 Table 1 provides the results from different research institutions and programs [ 15 15 ] . It can be seen that, in terms of the change in the amount of actinides or the change in reactivity between the Beginning Of the Circle (BOC) and Ending Of the Circle (EOC), the calculation results using COUPLE3.0 are in agreement with the results of all the research institutions. This indicates that the program is reliable for the burnup analysis of a fast neutron spectral device. 1 1 Institute or CODE 1 1 ANL (U.S.A) 1 1 CEA (FRANCE) 1 1 PNC (Jeff2) 1 1 PNC (Jeff3.2) 1 1 PSI (Swiss) 1 1 MOCUP 1 1 COUPLE3.0 8 1 Reactivity Changes between the BOC and EOC [%Δk/(k 1 k 2 )] 1 1   1 1 12.85 1 1 13.27 1 1 13.60 1 1 13.39 1 1 13.06 1 1 12.44 1 1 13.11 8 1 Nuclides Changes between the BOC and EOC [Δkg] 1 1 235 U 1 1 -5.6 1 1 -5.9 1 1 -5.8 1 1 -5.7 1 1 -5.5 1 1 -5.61 1 1 -5.78 1 1 238 U 1 1 -420 1 1 -411 1 1 -392 1 1 -390 1 1 -384 1 1 -394.5 1 1 -404.9 1 1 238 Pu 1 1 -50 1 1 -45 1 1 -50 1 1 -47 1 1 -43.4 1 1 -49.0 1 1 -49.4 1 1 239 Pu 1 1 -149 1 1 -174 1 1 -170 1 1 -131 1 1 -173 1 1 -160.9 1 1 -165.5 1 1 240 Pu 1 1 -38 1 1 -21 1 1 -32 1 1 -13 1 1 -34 1 1 -34.7 1 1 -29.4 1 1 241 Pu 1 1 -133 1 1 -139 1 1 -133 1 1 -122 1 1 -137 1 1 -136.5 1 1 -136.5 1 1 242 Pu 1 1 -29 1 1 -42 1 1 -31 1 1 -20 1 1 -25 1 1 -30.5 1 1 -29.5 1 1 241 Am 1 1 9.1 1 1 7.5 1 1 8.8 1 1 9.5 1 1 8.6 1 1 9.58 1 1 8.43 1 1 243 Am 1 1 31 1 1 44 1 1 33 1 1 34 1 1 33 1 1 31.2 1 1 31.7 1 1 242 Cm 1 1 3.7 1 1 5.2 1 1 4.8 1 1 4.5 1 1 4.1 1 1 4.63 1 1 4.54 1 1 244 Cm 1 1 4.1 1 1 7.4 1 1 5.3 1 1 5.3 1 1 5.4 1 1 4.48 1 1 4.95 4.2 ADS burnup benchmark test 2 In order to verify the use of COUPLE3.0 for ADS research, we analyzed the ADS benchmark problem of the International Atomic Energy Agency (IAEA) [ 16 16 ] . The benchmark relates to an ADS device, including a cylindrical subcritical 233 U-fueled reactor with 232 Th as a fertile material and natural lead as a reflector. The neutron source is placed in the center of the subcritical reactor to simulate the spallation neutron source. The benchmark problem is to test the burnup characteristics of the ADS under different subcritical conditions. In this study, we focus on the behavior changes of k eff and the intensity of the beam efficiency with an initial k eff value of 0.94. The results are compared with those of other institutions  [ 17 17 , 18 18 ] . Fig 4 Fig 4 shows the geometry of the benchmark problem. The core is divided into three zones: the first and second zones are for fuel with different enrichments of 233 U and the third zone is the fertile zone with 232 Th. The initial k eff value of 0.94 is obtained by adjusting the enrichment of 233 U. Fig. 5a Fig. 5a shows the changes in the k eff value during burnup. The comparison shows that the change in k eff calculated by COUPLE3.0 is close to the results of the Japanese institution. Fig. 5b Fig. 5b shows the changes in the beam intensity during the burnup in order to maintain the operating power at the level of 1500 MW th . The results show that the required beam intensity is largest for COUPLE3.0, and this is determined by the characteristics of the changes in k eff . A preliminary analysis indicates that this difference may be attributed to the use of the ENDF/B-VII.0 library in our calculations; the other institutions used the ENDF/B-VI library. Further studies should be conducted to determine the selection of a cross-section library.  In COUPLE3.0, the selection of the important nuclides is one approach for increasing the computational efficiency. Based on the IAEA-benchmark results, we performed several calculations to compare the computation times for different simulation results. These calculations were carried out using 64 CPUs with 1×10 6 source particles for each burnup step. In order to balance the precision level and the computation time, we set Importance Cut-off Values ( ICV) which represent the impacts on neutron absorption for all the nuclides. For each nuclide, its influence on the neutron transport is not considered in the MCNPX when its contribution fraction to the total value of neutron absorption in the specific zone is less than the ICV. This reduces the computation time and results in an acceptable precision level. Table 2 Table 2 shows that, when the ICV increases, the number of nuclides for the selected burnup zone decreases and the computation time decreases too. However, when the ICV is less than 1×10 -8 , there is no significant difference in the effective multiplication factors during the burnup ( Fig. 6 Fig. 6 ). When we set the cut-off value as 1×10 -8 , we can reduce the computation time by almost 40% when compared to an ICV of 1×10 -13 . 1 1 ICV 1 1 Number of nuclides in selected burnup zone 1 1 Computation time [mins] 1 1 1×10 -13 1 1 309 1 1 3701.7 1 1 1×10 -8 1 1 185 1 1 2269.8 1 1 1×10 -7 1 1 152 1 1 1900.4 1 1 1×10 -5 1 1 69 1 1 930.3 1 1 1×10 -4 1 1 22 1 1 164 4.3 Preliminary verification of proton irradiation 2 In order to verify the reliability of the newly added proton irradiation function in COUPLE3.0, we chose the target nuclide 210 Po as an example due to the lack of benchmark examples for the accumulation of radionuclides in the target area. Our analysis is based on the study of lead bismuth with a window target [ 19 19 ] , which is part of the project "Future advanced fission energy - ADS transmutation system" launched by the Chinese Academy of Sciences. The relevant parameters are as follows: the proton beam energy is 250 MeV and the proton beam intensity is 0.725 mA. Two methods were used for the analysis: method 1 is an analysis of the target nuclide accumulation behavior using the newly developed COUPLE3.0; method 2 is an analytical approach for the analysis of the 210 Po accumulation characteristics. The analysis and comparison are carried out for the bare target and the coupled condition of the target and the reactor, respectively; the k eff for the coupled condition is about 0.70. The analytical method is used to calculate the production of 210 Bi and 210 Po, which can be obtained by substituting MCNPX to study the evolution of the nuclides in the target.   <math display="block" id="M2"><mrow><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><msub><mi>N</mi><mrow><mtext>Bi</mtext></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mtext>p</mtext></msub><msub><mi>Y</mi><mrow><mtext>Bi</mtext></mrow></msub></mrow><mrow><msub><mi>λ</mi><mrow><mtext>Bi</mtext></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>λ</mi><mrow><mtext>Bi</mtext></mrow></msub><mo>⋅</mo><mi>t</mi></mrow></msup></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>N</mi><mrow><mtext>Po</mtext></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mtext>p</mtext></msub><msub><mi>Y</mi><mrow><mtext>Bi</mtext></mrow></msub><mo>+</mo><msub><mi>I</mi><mtext>p</mtext></msub><msub><mi>Y</mi><mrow><mtext>Po</mtext></mrow></msub></mrow><mrow><msub><mi>λ</mi><mrow><mtext>Po</mtext></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>λ</mi><mrow><mtext>Po</mtext></mrow></msub><mo>⋅</mo><mi>t</mi></mrow></msup></mrow><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>I</mi><mtext>p</mtext></msub><msub><mi>Y</mi><mrow><mtext>Bi</mtext></mrow></msub></mrow><mrow><msub><mi>λ</mi><mrow><mtext>Po</mtext></mrow></msub><mo>−</mo><msub><mi>λ</mi><mrow><mi mathvariant="normal">Bi</mi></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>λ</mi><mrow><mtext>Po</mtext></mrow></msub><mo>⋅</mo><mi>t</mi></mrow></msup><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>λ</mi><mrow><mtext>Bi</mtext></mrow></msub><mo>⋅</mo><mi>t</mi></mrow></msup></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mrow></mrow></math> where I p is the proton beam intensity, the unit is mA. <math id="M3"><mrow><msub><mi>Y</mi><mrow><mtext>Po</mtext></mrow></msub></mrow></math> and <math id="M4"><mrow><msub><mi>Y</mi><mrow><mtext>Bi</mtext></mrow></msub></mrow></math> are the yields of 210 Po and 210 Bi respectively, normalized to per mA. <math id="M5"><mrow><msub><mi>λ</mi><mrow><mtext>Po</mtext></mrow></msub></mrow></math> and <math id="M6"><mrow><msub><mi>λ</mi><mrow><mtext>Bi</mtext></mrow></msub></mrow></math> are the decay constants of 210 Po and 210Bi respectively, normalized to each day.  Fig. 7a Fig. 7a shows the accumulation of the 210 Po nuclides in the target area for the bare target operation, and Fig. 7b Fig. 7b shows the accumulation of the target 210 Po for the target reactor coupled condition. It can be seen that the results of the COUPLE3.0 simulation are in good agreement with the results obtained by the analytical method. This is because the cumulative production of the 210 Po nuclides is relatively low and the neutron nuclear reaction cross-section is low for the analysis situation; this means that the changes in the nuclides induced by the neutron can be neglected. In addition, the burnup of the core is not considered for the verification calculation. Therefore the core flux is constant.  The preliminary verification calculation shows that COUPLE3.0 is suitable for the simulation of proton irradiation. If the appropriate physical models of the nuclear reaction for medium- and high-energy protons and neutrons are chosen, we can accurately analyze the accumulation of the spallation products. At present, we use the CEM03 model [ 20 20 ] in MCNPX, and the calculated spallation products data fit well with the experimental data [ 13 13 ] .  4.4 Verification of burnable poison 2 The pressurized water reactor (PWR) fuel assembly of the Qinshan's phase II nuclear power plant, which contains 16 poisonous boron rods [ 21 21 ] , is used for the verification of the fixed-flux method for the burnable poison calculation function; the enrichment of the 235 U is 2.6% ( Fig. 8 Fig. 8 ). The main parameters are listed in Table 3 Table 3 . 1 1 Parameters 1 1 Value 1 1 Assembly type 1 1 17×17 1 1 Subtend distance (cm) 1 1 21.504 1 1 Cell pitch (cm) 1 1 1.26 1 1 Moderator temperature (K) 1 1 583.85 1 1 Fuel temperature (K) 1 1 924.15 1 1 Boron concentration (ppm) 1 1 0 1 1 Radius of fuel rods (cm) 1 1 0.4095/0.418/0.475 1 1 Radius of burnable poison (cm) 1 1 0.21375/0.23075/0.2413/0.4267/ 0.4369/0.48385/0.5625/0.6025 1 1 Composition of the burnable poison 1 1 0.74% 10 B, 55.16%O, 1.11%Al, 39.694%Si A comparison of the results of COUPLE3.0 and CASMO-3 ( Fig. 9 Fig. 9 ) indicates that the trend is similar for the two methods. In the beginning, the two curves fit well and at values greater than 15 GWd/tHM and the COUPLE3.0 calculation yields lower values. These differences may be caused by the differences in the nuclear database because we used ENDF/B-VII.0 in COUPLE3.0. 5. Conclusion 1 COUPLE3.0 can run the transport module in parallel in the MPI environment, which greatly improves the calculation accuracy and efficiency. Two types of neutron irradiation modes with constant power and a constant flux radiation model are included in the program and can be used for the calculation of the burnup of fuel and burnable poison.  The proton irradiation calculation module can be used to analyze the accumulation behavior of the residual nuclides in the spallation target based on the running history and the influence of the accumulated nuclides in the target area on the neutron behavior of the subcritical system. The improvements in COUPLE3.0 also include the selection of important radionuclides and the calculation of the branching ratios for the excited states of the radionuclides, which makes the calculation more accurate. The program has a user-friendly interface and is convenient for post-processing, extracting information, and determining activity information for the nuclides.  The calculation and analysis of a fast spectrum plutonium benchmark for incineration indicate that COUPLE3.0 is applicable and reliable for the analysis of the fuel burnup. The analysis of the IAEA-ADS burnup benchmark problem and the comparison with the results from other research institutes show that the k eff is lower than the average and the beam intensity is larger. We need to perform further analysis on the database selection and the calculation model. COUPLE3.0 cannot only be used for burnup analysis in ADS research, but also for the nuclide cumulative analysis of the target, the evaluation of the radiation characteristics of the spallation products, and their influences on neutronic characteristics, which can provide important reference data for the design and construction of ADS systems. References C.D. Bowman , E.D. Arthur , P.W. Lisowski , et al . Nuclear energy generation and waste transmutation using an accelerator-driven intense thermal neutron source . Nucl Instrum Meth A , 1992 , 320 : 336 - 367 . doi: 10.1016/0168-9002(92)90795-6 http://doi.org/10.1016/0168-9002(92)90795-6 Nuclear energy generation and waste transmutation using an accelerator-driven intense thermal neutron source H.A. Abderrahim , J. Galambos , Y. Gohar , et al . Accelerator and target technology for accelerator driven transmutation and energy production . DOE white paper on ADS , 2010 , 1 ( 1 ): 1 - 23 . http://192.107.175.200/~/media/hep/pdf/files/pdfs/ADS_White_Paper_final.pdf http://192.107.175.200/~/media/hep/pdf/files/pdfs/ADS_White_Paper_final.pdf Accelerator and target technology for accelerator driven transmutation and energy production K. Okumura , T. Mori , M. Nakagawa , et al . Validation of a continuous-energy Monte Carlo burn-up code MVP-BURN and its application to analysis of post irradiation experiment . J Nucl Sci Technol , 2000 , 37 ( 2 ): 128 - 138 . doi: 10.1080/18811248.2000.9714876 http://doi.org/10.1080/18811248.2000.9714876 Validation of a continuous-energy Monte Carlo burn-up code MVP-BURN and its application to analysis of post irradiation experiment J. Cetnar , W. Gudowski , J. Wallenius . MCB: A continuous energy Monte Carlo Burnup simulation code . Actinide and fission product partitioning and transmutation, 6th Information Exchange Meeting . Madrid, Spain , Dec.11-13,2000 . http://www.oecd-nea.org/pt/docs/iem/mol98/postersession/POSTSpaper12.pdf http://www.oecd-nea.org/pt/docs/iem/mol98/postersession/POSTSpaper12.pdf MCB: A continuous energy Monte Carlo Burnup simulation code S.C. Zhou , H.C. Wu , L.Z. Cao , et al . LAVENDER: a steady-state core analysis code for design studies of accelerator driven subcritical reactors . Nucl Eng Des , 2014 , 278 : 434 - 444 . doi: 10.1016/j.nucengdes.2014.07.027 http://doi.org/10.1016/j.nucengdes.2014.07.027 LAVENDER: a steady-state core analysis code for design studies of accelerator driven subcritical reactors L.F. Mao , T.Q. Dang , L Pan , et al . Preliminary analysis of polonium-210 contamination for China LEAd-based Research Reactor .  Prog Nucl Energ , 2014 , ( 70 ): 39 - 42 . doi: 10.1016/j.pnucene.2013.07.009 http://doi.org/10.1016/j.pnucene.2013.07.009 Preliminary analysis of polonium-210 contamination for China LEAd-based Research Reactor I. Remec , T. A. Gabriel , M. W. Wendel , et al . (2006) . Particle and radiation simulations for the proposed rare isotope accelerator facility .  Nucl Instrum Meth A: Accelerators, Spectrometers, Detectors and Associated Equipment ,  2006 , 562 ( 2 ): 896 - 899 . doi: 10.1016/j.nima.2006.02.103 http://doi.org/10.1016/j.nima.2006.02.103 Particle and radiation simulations for the proposed rare isotope accelerator facility Y.W. Wang , Y.W. Yang , P.F. Cui . Calculation and Analysis of Burnup and Optimum Core Design in Accelerator Driven Sub-critical System . Atom Energ Sci Technol , 2011 , 45 ( 6 ): 700 - 705 . doi: 10.7538/yzk.2011.45.09.1065 http://doi.org/10.7538/yzk.2011.45.09.1065 Calculation and Analysis of Burnup and Optimum Core Design in Accelerator Driven Sub-critical System H.Q. LI  and  Y. W.  YANG , The Development of Burn-up Program COUPLE by Combining MCNP and ORIGEN , Conference of Chinese Reactor Physics , Qinshan, China , Augest.8-12 (2004 ). The Development of Burn-up Program COUPLE by Combining MCNP and ORIGEN J. Zhao , Y.W. Yan , S.C. Xiao , et al . Burnup Analysis of Thorium-Uranium Based Molten Salt Blanket in a Fusion-Fission Hybrid Reactor . Fusion Sci Technol , 2013 , 64 ( 3 ): 521 - 524 . doi: 10.13182/FST13-A19145 http://doi.org/10.13182/FST13-A19145 Burnup Analysis of Thorium-Uranium Based Molten Salt Blanket in a Fusion-Fission Hybrid Reactor Y.G. Fu , J. Zhao , Y.W. Yang , et al . Physical scheme study of 10 MW accelerator driven sub-critical reactor . Atom Energ Sci Technol , 2013 , 47 ( suppl. ): 261 - 264 . doi: 10.7538/yzk.2013.47.SO.0261 http://doi.org/10.7538/yzk.2013.47.SO.0261 Physical scheme study of 10 MW accelerator driven sub-critical reactor D. Pelowitz , MCNPX User's Manual 2.7.0 . Los Alamos National Laboratory , Los Alamos, New Mexico  ( 2011 ). MCNPX User's Manual 2.7.0 A.G. Croff . ORIGEN 2.1—Isotope Generation and Depletion Code—Matrix Exponential Method. ORNL/NUREG/CSD-2 , 1980 . ORIGEN 2.1—Isotope Generation and Depletion Code—Matrix Exponential Method. ORNL/NUREG/CSD-2 M.B. Chadwick , P. Obložinský , M. Herman , et al . ENDF/B-VII. 0: next generation evaluated nuclear data library for nuclear science and technology . Nucl Data Sheet , 2006 , 107 ( 12 ): 2931 - 3060 . doi: 10.1016/j.nds.2006.11.001 http://doi.org/10.1016/j.nds.2006.11.001 ENDF/B-VII. 0: next generation evaluated nuclear data library for nuclear science and technology K. Wang , T. P. Lou , E. Greenspan , et al . Benchmarking and Validation of MOCUP , Proceedings of the International Topical Meeting on Advances in Reactor Physics and Mathematics and Computation into the Next Millennium , Pittsburgh, PA, U.S.A , 2000 . https://www.ipen.br/biblioteca/cd/physor/2000/physor/134.pdf https://www.ipen.br/biblioteca/cd/physor/2000/physor/134.pdf Benchmarking and Validation of MOCUP I. Slessarev , A. Tchistiakov . IAEA ADS-Benchmark Results and Analysis . Madrid: TCM-Meeting, Madrid, Spain , 1997 . IAEA ADS-Benchmark Results and Analysis X.F. Jiang , Z. Xie . Transport-burnup code systems and their applications for IAEA ADS benchmark . Ann Nucl Energy , 2004 , 31 ( 2 ): 213 - 225 . doi: 10.1016/S0306-4549(03)00205-6 http://doi.org/10.1016/S0306-4549(03)00205-6 Transport-burnup code systems and their applications for IAEA ADS benchmark L. Zhang , Y.W. YANG , Y.C. GAO . Preliminary physics study of the Lead-Bismuth Eutectic spallation target for China Initiative Accelerator-Driven System .  Nucl Sci Tech , 2016 , 27 ( 5 ): 120 . doi: 10.1007/s41365-016-0114-6 http://doi.org/10.1007/s41365-016-0114-6 Preliminary physics study of the Lead-Bismuth Eutectic spallation target for China Initiative Accelerator-Driven System S.G. Mashnik , A.J. Sierk , K.K. Gudima , et al . " CEM03 and LAQGSM03—new modeling tools for nuclear applications." J. Phys.: Conf. Ser , 2006 , 41 ( 1 ): 340 . doi: 10.1088/1742-6596/41/1/037 http://doi.org/10.1088/1742-6596/41/1/037 CEM03 and LAQGSM03—new modeling tools for nuclear applications." L. Chen , F. Ma , X.Y. Zhang , et al . Spallation yield of neutrons produced in thick lead target bombarded with 250MeV protons . Nucl Instrum Meth B , 2015 , 342 : 87 - 90 . doi: 10.1016/j.nimb.2014.09.020 http://doi.org/10.1016/j.nimb.2014.09.020 Spallation yield of neutrons produced in thick lead target bombarded with 250MeV protons Z.Y. Zhang , X.D. Li , Y.M. W , et al . Reactor Core Nuclear Design of Qinshan-II NPP . Nucl Power Eng , 1995 , 16 ( 1 ). http://en.cnki.com.cn/Article_en/CJFDTOTAL-HDLG199501000.htm http://en.cnki.com.cn/Article_en/CJFDTOTAL-HDLG199501000.htm Reactor Core Nuclear Design of Qinshan-II NPP 10.1007/s41365-018-0471-4 This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDA03030102). 2017-06-29 2017-11-22 2018-01-27 2018-07-24 2018-09-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0471-4 http://dx.doi.org/10.1007/s41365-018-0471-4 . 2018 Nuclear Science and Techniques 09 29