KnE Energy & Physics | International Conference on Nuclear Energy Technologies and Sciences (2015) | pages: 163-173

, , and

Introduction

AP1000 is a pressurized water reactor (PWR) that can produce a nominal power of 1154 MWe (3415 MWth). The reactor is designed by Westing house Co. based on the perfomance of the proven PWR. AP1000 reactor can be operated for 18 months for each cycle and has a life time of about 60 years. Today, AP1000 reactors are still under construction in several countries, such as Bulgaria, England, China, and U.S.A [1,2].

In the frame of evaluation of the core safety system in National Nuclear Energy Agency (BATAN), a simulation calculation using computer code on safety design parameter of AP1000 core has been done. Some related researches that have been done are analysis on the criticality of AP1000 core [3], analysis on the coefficient reactivity of the 1000 MWe PWR [4], kinetic parameter calculation of AP1000 core [5,6], and analysis of using mixed oxide (MOX) fuel on AP1000 core design [7,8,9].

At the first operation cycle, AP1000 core uses 3 type enrichments of UO2 fuels, that are 2.35 w/o, 3.40 w/o and 4.50 w/o. To compensate the effective multiplication factor (keff) or excess reactivity in the beggining of cycle of the AP1000 core, borid acid which is dissolved in moderator, Pyrex absorber rod and Integrated Fuel Burnable Absorber (IFBA) are used. IFBA is a burnable absorber made frome ZrB2 which is integrated into UO2 fuel, whereas Pyrex is an absorber rod made from B2O3 which is inserted into the guide tube [10,11]. In the core, UO2 fuel, IFBA and Pyrex are arranged to form a configuration in such away that it can generate the neutronic parameters which meet the safety criteria, so the core configuration is safe and feasible to be operated. Figure 1 shows distribution of fuel and Pyrex layout on the quarter core of the AP1000 according to the inial design.

Figure 1

Fuel and Pyrex layout in a quarter core of the AP1000

Images/Fig1.jpg

Part of UO2 pellets are coated with absorbent material ZrB2 (called with IFBA fuel) and the remaining are without absorbent material. Fuel assembly is composed by a 17×17 fuel arangement with a certain ratio between the number of UO2 to IFBA. Based on the UO2 enrichment and number of IFBA, there will be 9 different types of fuel assemblies. Each fuel assembly with their respective positions are predetermined arranged to form the reactor core. Due to the possibility of human errors, misposition of fuel assemblies might happen when operators load thoes fuel assemblies into the core. Therefore, it is necessary to investigate the effects of those misposition to the AP1000 neutronic parameters.

In this research, analysis of changes in the neutronic parameters of the AP1000 core due to the mispositioning of two adjacent fuels was performed. Analyze is imposed to the core at hot zero power (HZP), beginning of cycle (BOC), and zero xenon condition. Neutronic parameters calculation was carried out by means of CITATION module of the Standard Reactor Analysis Code 2006 (SRAC2006) [12] computer code for ¼ core model in 2 dimension geometry. CITATION is a computer code, which uses diffusion theory in the multiplication factor calculation with finite different method to solve Boltzman equation. One of input data that is fuel macroscopic crossection was prepared using PIJ module. PIJ is a computer code based on transport theory with neutron collision probability method. The CITATION module of the SRAC2006 computer code has been validated for the criticality value of conventional PWR with good results [13]. From the analysis results, it is expected to know the safety characteristics of the AP1000 core’s design.

Methodology

Fuel Crossection Calculation

The homogenization macroscopic crossection of the fuel calculation performed by PIJ module of the SRAC2006 through condensation 107 to 10 group of neutron energy. Cross section data library used in the calculation was JENDL-3.3. Input data required for those calculation is the size, dimensions and material that make up the fuel. Cross section calculations carried out on the &¼ part of the fuel lattice model in the two-dimentional geometry. Fuel lattice of the AP1000 core composed by 17×17 of fuel pin cell arrangement with water as moderator. Model of the fuel pin cell as showed at Figure 2. is composed from fuel pellets, gap, cladding and moderators in the square form with size of 1.25984 cm × 1.25984 cm. The outer part of pin cell is a moderator that consist a mixing of H2O and boron solution with concentration of 1382 ppm. There are two sort of fuel pellet, standard UO2 fuel and IFBA. Fuel pellets standard are UO2 (2,35 w/o, 3,40 w/o dan 4,50 w/o) and IFBA are UO2 fuel coated with ZrB2. Dimension of pellet, gap of the He gas and cladding of the ZIRLO are r1 = 0,409575 cm, r2 = 0,426085 cm, and r3 = 0,483235 cm, respectively. Whereas IFBA fuel is same as standard fuel which ZrB2 coated with amount of B-10 content of 0,772 mg/cm. Composition of ZIRLO cladding are Zr 97,85%, Fe 0,15%, Sn 1,0%, dan Nb 1%.

Figure 2

Model of the fuel pin cell

Images/Fig2.jpg

Figure 3 shows a quarter model geometry of the fuel lattice. Whole of the fuel lattice on the AP1000 core have size 21.50364 cm × 21.50364 cm that contains 264 fuel, 25 guide tubes and with moderator H2O. Guide tube is provided for absorber Pyrex, or Ag-In-Cd control rod, or moderator. Based on the amount of the fuel IFBA, then fuel lattice are divided into fuel lattice without IFBA, fuel lattice with 28 IFBA, 44 IFBA, 72 IFBA, 112 IFBA and 88 IFBA.

Figure 3

A Quarter geometry model of the fuel lattice

Images/Fig3.jpg

Neutronic Parameters Calculation of The AP100 Core

Neutronic parameters has been calculated for a quarter 2 dimensional geometry of AP1000 core. Input needed to core calculation using module CITATION of SRAC2006 are nuclide density of non fuel material, dimension size of the core, and homogenization macroscopic crossection of the fuel. Output from the core calculation using CITATION are k-eff, power distribution of the fuel, prompt neutron life time and generation time, as well as delayed neutron fraction.

As shown in Figure 4, two-dimentional model geometry of the AP1000 core is built from fuel lattice of 3 type of fuel enrichment that are 2.35 w/o, 3.40 w/o dan 4.50 w/o which relates to region number 1, 2 and 3 respectively. Then, the core is surrounded by reactor baffle (no.4), coolant (no.5 and 7), reactor barrel (no.6), reactor vessel (no.8) and vacum (no.9). Size and material composition of the core structure are reactor baffle from SS304, reactor barrel from SS304 with din/dout 339,725 cm/349,885 cm, reactor vessel from SS304 with din/dout 398,8 cm/420.1 cm.

Figure 4

Radial core structure of the two-dimentional model AP1000 core

Images/Fig4.jpg

Mispositioning of Two Adjacent Fuels Assumption

In this study, analysis of the neutronic parameters change due to mispositioning of two adjacent fuel has been done. Such as illustrated in Figure 5., it is assumed that mispositioning have occur of two adjacentt fuel in the AP1000 core (Case 1 to 6). While the absorber rod Pyrex position was fixed, do not change from the design data. Case 1 to 6 are all located in column 8 start from the center to the edge of the core. Case 1 to 5, the mispositioning accour between UO2 3,40 w/o 24 Pyrex 88 IFBA fuel with UO2 2,3 w/o fuel. The mispositioning fuel in the core were G-8↔F-8, F-8↔E-8, E-8↔D-8, D-8↔C-8 and C-8↔B-8. Whereas, mispotition fuel in the case 6 was between UO2 4,45 w/o 12 Pyrex 88 IFBA and UO2 2,35 w/o fuel.

Figure 5

Misposition fuel located in column 8 of the quarter AP1000 core for case 1 to 6

Images/Fig5.jpg

Results and Discussions

Change in the Effective Multiplication factor

Calculation results of the keff value changes due to the misposition fuel of AP1000 core are shown in Table 1. The standard core that has same configuration of the fuel with design data showed that the value of keff is 1,001910. Whereas the value of keff for the first operating cycle of the AP1000 core at Hot Full Power, BOC, Zero Xenon condition in the reference data is 1,00. Therefore, differences among calculation result with reference data is 0,001910 or 0,19%. So that, it can be said that the keff value of AP1000 core calculation result using module CITATION of the SRAC2006 computer code is about the same with refferences value.

Table 1

Change of the keff value due to the mispositioning fuel in the AP1000 core

No. Case of AP1000 Core keff Change of the keff value
∆k=kn- k0
1. Standard Core (k0) 1.001910
2. Case 1 (k1) 1.004707 0.002797 0.28
3. Case 2 (k2) 1.004228 0.002318 0.23
4. Case 3 (k3) 1.003691 0.001791 0.18
5. Case 4 (k4) 1.003488 0.001578 0.16
6. Case 5 (k5) 1.004056 0.002146 0.21
7. Case 6 (k6) 1.001718 −0.000192 −0.02

From this table it can also be known that except case 6, mispositioning of the fuel in tne AP100 core cause a small increament in the value of keff. As it is known that in the case 1 to 5 mispositioning of fuel occured between UO2 3.40 w/o 88 IFBA fuel with UO2 2.35 w/o fuel with each position are F-8↔E-8, E-8↔D-8, D-8↔C-8 and C-8↔B-8, repectively. It causes the location of 2 absorber Pyrex and IFBA are separated with different grid position. While in the initial design, 2 absorber Pyrex and IFBA located on the same grid position. Therefore, it would result in a decrease in thermal neutron absorption. Furthermore it will increase the value of thermal neutron utilization factor whichby, then, rising keff value of the core. The biggest keff value changes due to misposition fuel in the core is 0.002797 or approximately 0.28 % , that is for the case 1 (G - 8↔F - 8). This is caused by the UO2 fuel with 2.35 w/o enrichment without absorbent material IFBA moved its position to the centre of the core. So due to the influence of the interaction of neutrons produced by the two fuel UO2 2.35 w/o in the centre of core resulting the increase of the neutron thermal utilization. While the case 6, mispositioning of the fuel occurs in the outer position of the core. That are mispositioning among UO2 2.35w/o fuel (A-8) with UO2 4.5w/o 88 IFBA fuel (B-8). In this case 6, the keff value decreases about −0.000192 (−0.02 %), which is from 1.001718 to 1.001910. This is caused by the fuel with lowest enrichment (UO2 2.35 w/o) in outer position of the core (A-8) which also contained 12 absorber rods Pyrex. Thereby decreasing neutron thermal utilization factor and it will also decreasing of the keff value.

Power Factor Distribution

Figure 6 shows calculated power factor distribution at first cycle operation of the AP1000 at HFP, BOC, Zero Xenon conditions compared to the reference value. It shows that the maximum power factor from SRAC 2006 CITATION module calculation is 1.246 at position D-8. Whereas from the reference, its value is 1.279 and at position H-8. The difference of theis maximum value is 0.033 or 2.58%. At position D-8, the reference value is 1.254 or 0.037 (2.92%) difference. The maximum differences of the power factor values are occured at positions B-9 and G-14 which are amount to 0.088 (9.63%) and 0.086 (9.43%). The differences of the calculated and reference values can be resulted from the differences in the computer codes module and input data used by the reference.

Figure 6

Comparison calculation result and refference value of power factor distribution of the ¼ AP1000 core at hot full power, BOC, zero xenon condition

Images/Fig6.jpg

Nevertheless, based on the power distribution, it can be said that the both results are in a good agreement. Fuel lattice with high power factor (pi > 1,2) occurs at UO2 , 2.35 w/o enrichment without absorber IFBA or Pyrex. On the other hand, fuel lattice with small power factor (0.6) occurs for UO2 , 4.45 w/o enrichment. It is caused by its positions are in the outer part of the core that has lower neutron fluxes compared to that in the inner part. It can be concluded that the value of power factor depends besides on the enrichment also on the in-core position and the existence of absorber in the fuel as well.

Calculation results of power factor distribution for ¼ core of the AP1000 due to fuel mispositioning case 1 to 6 is presented in Figure 7. From the table it is found that maximum power factors for case 1 to 5 occurs in the same position that ic at H-12. Those values are 1.274, 1.281, 1.299, 1.382 and 1.377 consecutively. Thhese maximum power factor values still lays below the predetermined safety limit of 1.6. The figure shows also that mispositioning of fuel causes slightly changes in power factor distribution compared to initially designed.

Figure 7

Power Factor Distribution in the ¼ AP1000 Core for Case 1 to 6 in Term of The Misspositioning Fuel at Hot Full Power, BOC, Zero Xenon Condition

Images/Fig7.jpg

Prompt Neutron Parameter

Calculation results of prompt neutron lifetime and prompt neutron generation time of AP1000 core due to fuel mispositioning are illustated in Table 2. The value of the standard core (20.7006 μs) is slightly higher than the reference value (19.8 μs) or about 4.5% difference. For all cases of fuel mispositioning (1 – 6), the maximum changes of prompt neutron lifetime and prompt neutron generation time are 2.35% and 2.06% that is for Case 1. In other words, mispositioning of fuel does not give any significant changes for prompt neutron lifetime and prompt neutron generation time.

Table 2

Change of prompt neutron lifetime (𝓁) and prompt neutron generation time ()core AP1000 due to fuel mispositioning

prompt neutron lifetime,
ℓ (µs) prompt neutron generation time,
ʌ (µs)
AP1000 Core Calculation Difference (%) Calculation Difference (%)
Standard core 20.7006 20.6611
Case 1 21.1862 0.4856 2.35% 21.0870 0.4259 2.06%
Case 2 20.9238 0.2272 1.10% 20.8357 0.1746 0.85%
Case 3 21.0300 0.3294 1.59% 20.9527 0.2916 1.41%
Case 4 20.7698 0.0693 0.33% 20.6976 0.0365 0.18%
Case 5 20.9406 0.2400 1.16% 20.8560 0.1949 0.94%
Case 6 20.6070 −0.0936 −0.45% 20.5716 −0.0895 0.43%

Delayed Neutron Parameter

Table 3 shows the comparison of calculation and reference value of the delayed neutron fraction (2eff) for all 6 cases of fuel mispositioning at AP 1000 core. The design or reference value is 7.5×10−3, so deviation of calculation result is 0.5759 or 7.68%. Then, it can be said that calculation of kinetics parameter using SRAC2006-CITATION give a good result.

it is also shown from the table that delayed neutron fraction (2eff) of AP1000 core a very small change either up or down with maximum change of 5.02×10−6 or 0.073% for cases of fuel mispositioning.

Table 3

Change of delayed neutron fraction, βeff core AP1000 due to fuel mispositioning

Delayed neutron fraction, βeff
AP1000 Core Calculation Difference (%)
Standard core 6.92410×10−3 - -
Case 1 6.92046×10−3 −3.64×10−6 −0.050%
Case 2 6.92592×10−3 1.82×10−6 0.026%
Case 3 6.92028×10−3 −3.82×10−6 −0.055%
Case 4 6.92912×10−3 5.02×10−6 0.073%
Case 5 6.92391×10−3 −1.90×10−7 0.003%
Case 6 6.92342×10−3 −6.80×10−7 0.010%

Conclusion

Neutronic parameter analysis has been done due to the mispositioning of the fuel in the first operating cycle of the AP1000 at Hot Full Power, BOC, Zero Xenon condition. From the analysis, it can be concluded that misposition of the fuel in the AP1000 core insignificantly change the neutronic parameters, so it will not reduce the performance of the core.

Acknowledgements

Authors would like to say thank you to Prof. Drs. Surian Pinem, M.Si. for his priceless advice and also to the Head of the Nuclear Reactor Physic And Technology Division as well as the Director of the Center for Nuclear Reactor Technology and Safety for supporting this research.

References

[1] 

Balendra Sutharshan Meena Mutyala Ronald P Vijuk Alok MishraAsian Nuclear Prospects 2010-The AP1000TM Reactor: Passive Safety and Modular Design-Journal of Energy Procedia002010000000.

[2] 

W. E. Cummins M. M. Corletti T. L. SchulzWestinghouse AP1000 Advanced Passive PlantProceedings of ICAPP ’03, Cordoba, Spain, May47, 2003, Paper 3235

[3] 

Jati Susilo RokhmadiAnalysis of Criticallity Parameter of AP1000 Reactor Core Using SRAC-CITATIONProsidings of Scientific Meeting and Presentation on Basic Research and of Nuclear Science and TechnologyYogyakarta, 24–25September201394102

[4] 

T. M. Sembiring S. PinemEvaluation of Moderator Temperature Coefficient of Reactivity for The 1000 MW PWR Nuclear PlantProsidings of National Seminar 17th Technology And Safety of NPP and Nuclear FacilityYogyakarta, October12011164174

[5] 

Jati Susilo Lily SuparlinaAnalysis of Kinetic Parameter As a Function of Burn up of the AP1000 CoreProsidings of National Seminar 19th Technology And Safety of NPP and Nuclear FacilityYogyakartaSeptember24-2520131323

[6] 

Tagor Malem Sembiring Surian PinemEvaluasi of Delayed Neutron Parameter of AP1000 ReactorProsidings of National Seminar 18th Technology And Safety of NPP and Nuclear FacilitySeptember2920126775

[7] 

Tukiran SurbaktiOptimization of MOX Fuel for Core Design of AP1000 ReactorProsidings of National Seminar 19th Technology And Safety of NPP and Nuclear FacilityYogyakarta, September24–25, 20132334

[8] 

RokhmadiAnalisys of Control Rod Reactivity on The PWR AP1000 Design Fueled MOXProsidings of National Seminar 19th Technology And Safety of NPP and Nuclear FacilityYogyakarta, September24–25, 20133445

[9] 

Tukiran SurbaktiEvaluation of Kinetic Parameter for The Safety of The AP1000 Core Used MOX FuelProsidings of National Seminar 20th Technology And Safety of NPP and Nuclear FacilityPontianakJune192014859868

[10] 

[11] 

T. L. SchulzWestinghouse AP1000 advanced passive plantJournal of Nuclear Engineering and Design236200615471557.

[12] 

Keisuke Okumura Teruhiko Kugo Kunio Kaneko Keichiro TsuchihashiSRAC2006; A Comprehensive Neutronics Calculation Code SystemJAERI-Data/Code 2007-004Japan Atomic Energy AgencyJanuary2007

[13] 

Jati Susilo ”Verification of MVP-II and SRAC06 Code for the Case VERA Benchmark Reactor CoreJournal TDM162June2014ISSN 1411-240X, No. Acreditation: 402/AU2/P2MI-LIPI/04/20127588

FULL TEXT

Statistics

  • Downloads 28
  • Views 183

Navigation

Refbacks



ISSN: 2413-5453