2.1. AP1000 Reactor Core Description
The AP1000 pressurized water reactor has a thermal power of 3400 MW. The core consists of 157 fuel assemblies arranged in a 17 × 17 elongated fuel assemblies arrangement [
21]. The AP1000 reactor core is loaded with three fuel assemblies at different enrichment levels staggered in a tessellated pattern, as shown in
Figure 1.
The numbers indicate the amounts of burnable absorber rods of each type in a given fuel assembly: IFBA (I) and Pyrex (P). Each assembly has 289 positions, of which 264 positions are occupied by fuel elements and 24 positions are occupied by guide tubes to provide positions for the core functional assemblies. The guide tubes, together with the centrally located neutron measurement tubes and eight positioning grids along the height direction, form the skeleton of the fuel assembly to provide support.
Figure 2 shows the RMC model of the fuel assembly without the core functional assembly and with control rods. In addition,
Appendix A shows the 9 fuel assemblies in the AP1000.Details of the A1000 fuel assembly are shown in
Table 1 [
22]. The center of the single fuel rod is a low-enrichment-sintered UO2 ceramic core block, externally filled with helium and encapsulated in a zirconium–niobium alloy cladding.
The burnable poison rods are primarily used for the first loading of the core, insertion into the fuel assembly, and occupying the fuel rods. The AP1000 uses separated burnable poison rods, known as Wet Annular Burnable Absorbent (WABA) [
23]. It is used to suppress the initial excess reactivity of the core at the start of the cycle and to reduce the RTC.
Figure 1 shows the number of the Integral Fuel Burnable Absorber (IFBA) and the Pyrex Burnable Absorber.
The control rods are used to regulate the reactivity of the reactor during operation or for fast shutdown in emergencies. According to the purpose of control rods, control rods are categorized into regulating rods and shutdown rods. There are four sets of shutdown rods (SD) [
24,
25], each with eight bundles of control rod assemblies for rapid shutdown. There is only one set of Axial Offset Rods (AO) [
24,
25], consisting of nine bundles of control rod assemblies for axial power distribution control. There are six sets of Shim Rods (M) [
24,
25], which are used to compensate for daily reactivity changes.
Figure 3 shows the distribution of control rods in the 1/4 core. There are 69 control rods in the AP1000 core. Reactor startup, shutdown, and fuel consumption compensation all require regulation of boron concentrations. Material and geometry data describing the AP1000 core are taken from the existing literature [
26].
2.2. Calculation Methods
This paper investigates the effects of fuel temperature, boron concentration, control rod position, and moderator density variations on AP1000 reactor reactivity. The paper calculates the temperature reactivity coefficients and control rod worth for different influences.
We also studied the variation of effective delayed neutron fraction (βeff) and effective neutron production time (Λeff) under different conditions. The effects of the core parameter change on the core neutron flux and relative power density distributions are studied comparatively to analyze the state of the AP1000 reactor core under accident conditions.
The study uses the Monte Carlo particle transport program, RMC [
27], to create a detailed three-dimensional model of the AP1000 reactor core based on the AP1000 reactor data from Westinghouse [
26]. RMC is the physics computational core of the multi-physics, multi-size coupled nuclear energy system numerical analysis platform. RMC uses the ENDF/B-VII.1 database for complex 3D models. In the calculations, the RMC critical card is set to 20,000 particles per iteration and 500 iterations in total, and the initial 100 inactive generations of neutron records are skipped before counting. Due to the incomplete cross-section data of nuclides at different temperatures in the RMC3.5 software, the NJOY2016 software was used to generate the nuclear cross-section data required for the temperature of the study.
Figure 4 shows the detailed 3D model of the AP1000 core created by the RMC3.5 software based on AP1000 reactor data provided by Westinghouse [
26].
After the core model was completed, the reliability of the AP1000 core model first needed to be verified. The K
eff results calculated by the RMC3.5 software are compared with the values reported by Westinghouse [
26]. In the second step, the variation of reactivity was investigated computationally over a wide range of temperatures, limiting values of boron concentration, different densities of moderator, and different positions of the control rods. As shown in
Table 2, twenty study conditions are set up according to the three states of AP1000 (CZP, HZP, and HFP). The study conditions include normal and core overheating, boron dilution, and other accidents. In this part of the study, the reactivity of the AP1000 core is evaluated by calculating K
eff (including the Doppler Effect, temperature effect, moderator density effect, boron worth, and control rod worth). The above core configurations are symbolized as C1–20.
The main calculations are as follows.
The Doppler effect calculation conditions are from normal operating conditions to core overheating (from 565 K to 732.5 K, 900 K, 1500 K, and 2000 K). Two different sets of calculations were conducted. One is based on soluble boron concentrations of 0 ppm (C1\C2, C1\C3, C1\C4, and C1\C5), and the other is based on soluble boron concentrations of 800 ppm (C9\C10, C9\C11, C9\C12, and C9\C13).
The effect of moderator density is based on the presence or absence of soluble boron, and keeping the fuel temperature constant was explored and calculated (C3\C6, C5\C7, C13, C15, and C11\C14).
The research on soluble boron worth and boron concentration coefficient mainly focuses on the impact of changes in soluble boron concentration on reactivity at different temperatures (C9\C1, C10\C2, C11\C3, C12\C4, and C13\C5).
The study of the control rod worth is determined by the reactivity difference with and without control rod insertion in the core (C9, C10, C17\C11, C13\C19, and C15\C20). For C9 and C10, additional comparisons were made with the control rods fully inserted under the same conditions.
In addition, the effects of temperature variation, coolant density variation, boron concentration reduction, and the position of the control rods in the duct on the neutron flux and relative power density distributions are also investigated.
The neutron flux and relative power density distributions of the fuel assemblies at each operating condition are calculated by using the notation card Type in RMC. The variation of neutron flux and relative power density peak is also studied.