Integrity Evaluation of a Reactor Pressure Vessel Based on a Sequential Abaqus-FRANC3D Simulation Method

+e safety-risk pressurized thermal shock (PTS) have on a reactor pressure vessel (RPV) is one of the most important studies for the lifetime ageing management of a reactor. Several studies have investigated PTS induced by postulated accidents and other anticipated transients. However, there is no study that analyzes the effect of PTS induced by one of the most frequent anticipated operational occurrences—inadvertent operation of the safety injection system. In this paper, a sequential Abaqus-FRANC3D simulation method is proposed to study the integrity status of an ageing pressurized water reactor subjected to PTS induced by inadvertent actuation of the safety injection system. A sequential thermal-mechanical coupling analysis is first performed using a three-dimensional reactor pressure vessel finite element model (3D-FEM). A linear elastic fracture mechanics submodel with a postulated semielliptical surface crack is then created from the 3D-FEM. Subsequently, the submodel is used to evaluate the stress intensity factors based on the M-integral approach coupled within the proposed simulation method. Finally, the stress intensity factors (SIFs) obtained using the proposedmethod are compared with the conventional extended finite elementmethod approach, and the result shows a good agreement. +e maximal thermomechanical stress concentration was observed at the inlet nozzleinner wall intersection. In addition, +e ASME fracture toughness of the reactor vessel’s steel compared with SIFs show that the PTS event and crack configuration analysed may not pose a risk to the integrity of the RPV.+is work serves as a critical reference for the ageing management and fatigue life prediction of reactor pressure vessels.


Introduction
e effect of ageing degradation mechanisms in nuclear power plants (NPPs) may result in a substantial loss of plant availability and costly part replacement [1]. In addition, the re-licensing regime of ageing NPPs is largely premised on the integrity assessment of critical components such as the reactor pressure vessel (RPV). e structural integrity of the RPV is a key safety priority in the operation of ageing pressurized water reactor (PWR) NPPs since it technically determines the feasible lifetime of the reactor [2]. Furthermore, under some NPP transient conditions, a small defect of the size of the nondestructive testing limit in an ageing RPV may rapidly grow leading to damage or failure. erefore, structural integrity analysis of vital components such as RPV subjected to potential transient loading from anticipated operational occurrences (AOOs) or postulated accidents (PAs) is essential to guarantee the safety of the whole NPP. In addition, integrity assessment results inform operators about the development of predictive maintenance and ageing management strategies that may prevent catastrophic failures of the RPV equipment.
A major risk to ageing RPVs is the exposure to pressurized thermal shock (PTS) induced by temperature and internal pressure loads from emergency cooling water triggered by AOOs or PA events [3].
e PTS initiating events that may arise during the operation of NPPs are mainly categorized based on their frequency of occurrence in the plant's lifetime. AOO transients, categorized as condition II events, must not escalate to condition III or IV postulated accidents under any NPP incident. Generally, AOO transients such as the inadvertent actuation of the safety injection system (SIS) or make-up systems occur in pressurized water reactors. is event involves the filling of the PWR's pressurizer and subsequent discharge of the water content through water safety or relief valves [4]. e filling of the pressurizer of PWRs may be initiated through an AOO heat-up or water mass addition event. In both events, the reactor coolant surges into the pressurizer, fills it, and departs the RCS through power-operated relief or pressurizer safety valves. e heat-up events usually cause the pressure and water levels of the pressurizer to rise, as the reactor coolant expands, until the reactor is automatically tripped by the sense of a pressurizer high-level signal. In a water mass addition event, such as an inadvertent operation of the emergency core cooling system (ECCS), the inadvertent opening of a pressurizer power-operated relief valve, or a chemical and volume control system (CVCS) malfunction, the tripped reactor will not cause the stoppage of the flow of water from the ECCS into the RCS. e pressurizer level will continue to increase until it reaches the water-solid pressurizer condition or a reactor operator terminates the safety injection flow, following the required plant emergency procedures, before reaching a more severe condition III event [4,5]. It is worth noting that a condition II small break loss of coolant accident (SBLOCA) that is initiated from an anticipated transient has a higher frequency of occurrence in an NPP than a condition III SBLOCA. Also, Mukin et al. averred from a thermal-hydraulic screening analysis that the inadvertent opening of a pressurizer power-operated relief valve, which is an AOO event, contributes significantly to PTS due to its high chance of occurrence [6]. erefore, the high frequency of AOO events occurring in the lifetime of a NPP, in comparison with the extensively studied postulated accidents, may have a critical effect on the physical characteristics of the RPV and other essential reactor equipment [4]. Several studies on the integrity assessment of RPV subjected to PTS loadings are mainly on the PTS initiated by postulated accidents [3]. Chen et al. have performed the structural integrity assessment of the RPV under the PTS loading initiated by a loss of coolant accident (LOCA) [7]. Also, Jhung et al. have evaluated the structural integrity of RPV based on pressurized thermal shock caused by SBLOCA [8]. Murtaza and Hyder have investigated the fracture mechanics analysis of the set-in nozzle of a RPV based on SBLOCA and Rancho-Seco transient events [9]. Zhu et al. have evaluated the fatigue failure of turbine bladed disks under cyclic loadings using a probabilistic modeling approach [10,11]. Similarly, He and Isozaki have performed an integrity assessment of the RPV of Qinshan NPP under SBLOCA, LBLOCA, and Rancho-Seco PTS transient events [12]. However, the evaluation of the structural integrity of an ageing RPV subjected to PTS from anticipated transients is scarcely studied, despite its relatively high frequency [4,13,14]. Furthermore, nuclear safety requirements obligate operators to investigate all potential AOO-and PA-induced PTS events. is forms part of the overall integrity assessment and justification of ageing NPP RPV safety margins in license extension applications [15].
Additionally, the repeated injection of emergency coolant in the SIS inadvertent operation induces thermal stresses in regions of the PWR RPV in contact with the cooling water [16]. ese thermal stresses coupled with high internal pressure can cause preexisting defects in the RPV material to increase, leading to damage, leakage, or the possibility of failure [17,18]. erefore, a complete safety assessment of ageing NPP must include the integrity analysis of RPV subjected to PTS loads from both anticipated transients and postulated accidents. is is critical for decisions on continuous operation, ageing management, shutdown, or service life extension for ageing NPPs.
Consequently, this study presents a sequential Abaqus-FRANC3D FE method for the evaluation of the structural integrity of an ageing double-circuit PWR reactor pressure vessel subjected to PTS loadings caused by the inadvertent operation of a safety injection system. A sequential thermal-mechanical coupling analysis is first performed to determine the highest stress location in a 3D-FE RPV model developed in Abaqus. A numerical submodeling technique is then applied to evaluate the stress intensity factors (K I ) of an assumed semielliptical surface crack set at the highest stress location in the RPV submodel. Subsequently, K I values at the deepest crack point are evaluated using the M-integral approach and extended finite element method (XFEM). Finally, K I is compared with ASME fracture toughness of the RPV steel (K IC ) to determine the safety risk the inadvertent operation of the SIS poses on the RPV. e novelties in this work are summarized as follows: (1) e work proposes a sequential Abaqus-FRANC3D simulation method to investigate the integrity risk posed by the most frequent anticipated operational occurrence (inadvertent operation of SIS) on an ageing RPV. is approach reduces the computational cost in the FE estimation of stress intensity factors in an integrity evaluation of complex geometries.
(2) e integrity parameter, K I , calculated using the M-integral approach coupled within the proposed method is in agreement with the conventional XFEM approach, confirming its accuracy. (3) is investigation serves as a preliminary framework for evaluating the effect of other anticipated operating transients on the integrity of an ageing RPV.

System Analysis and Reference Transient.
e SIS is a critical safety feature in a nuclear power plant that injects cold water into the RCS during loss of coolant accidents (LOCA). A typical SIS of a double-circuit PWR has four high-pressure safety injection (HPSI) pumps. e inlets of two of the pumps are connected to the refueling water storage tank (RWST) via power operation isolation and check valves. Also, all inlets of the two pumps are connected to the outlets of the shutdown margin pumps via power operation isolation valves. e outlets of the pumps are connected to a common pipe header with four discharge pipelines connected to the hot legs or cold legs of the RCS. Under the condition of the safety injection, the HPSI pumps enable cold water from the RWST to be pumped into the reactor core [19]. Events leading to the inadvertent injection of colder water directly into the downcomer result in the increase of the core power and pressurizer level. Also, the frequent occurrence of inadvertent operation of SIS during the lifetime of PWRs either by a control system malfunction or by an operator error may cause thermal stresses in the inlet nozzle and inner surface of the RPV wall [4].

RPV Material Properties.
e RPV material is made of ferritic low alloy SA508 Class 3 steel. e material is characterized as homogenous, isotropic, elastic, and temperature-dependent. e RPV thermomechanical properties (thermal conductivity λ, mean thermal expansion coefficient α, elastic modulus E, and specific heat capacity C) are listed in Table 1 [12]. e mean linear thermal expansion coefficients used for the thermomechanical stress calculations were converted to the stress-free reference temperature of 289°C following the conversion approach described in [20]. In addition, the yield strength of vessel material at room temperature, its density, and its Poisson's ratio are 450 MPa, 7600 kgm −3 , and 0.3, respectively, as reported in [12,21]. Figure 1 represents the stress-strain characteristics of the RPV material obtained from experimental tensile tests performed at room and high temperatures [9]. Also, the Ramberg-Osgood power law expressed in equation (1) depicts the nonlinear behavior of the vessel material: where ε, E, σ 0 , n, and α, respectively, denote the total strain, Young's modulus, yield strength, hardening exponent, and material constant.

ASME RPV Material Fracture
Toughness. e high-energy neutron from the reactor core causes a decrease in the RPV material's fracture toughness and an upward shift in the nil-ductility transition temperature. In this study, the fracture behavior of the referenced RPV material (K IC ) is estimated using the main steel parameters in Table 1 and ASME empirical formula expressed as where T and RT NDT are material temperature and nilductility transition temperature, respectively [22][23][24]. e temperature, RT NDT , is given as follows: where where RT NDT(i) , M, ΔRT NDT , f, and CF are the reference temperature of un-irradiated RPV material, average safety margin, mean transition temperature shift, neutron irradiation fluence (10 19°n /cm 2 ), and chemistry factor of RPV material, respectively. Also, n, A i , and f i are the radiation monitoring points, measured values for RT NDT(i) , and fast neutron fluence per irradiation monitoring point, respectively. σ Δ and σ I represent the standard deviation for ΔRT NDT and RT NDT(i) , fs is the neutron fluence at the inner surface of the vessel, and a (mm) is the crack depth. It is conservatively assumed that the inner surface and inlet nozzle of the RPV will be subjected to a neutron fluence of 4 × 10 19 n/cm 2 at the end of life of the referenced doublecircuit PWR. Also, the main random variables of the RPV steel used in this study, shown in Table 2, were derived using the prediction function [25]. From equations (2)-(6) and assumed random variables in Table 2, RT NDT at 1/4 vessel wall thickness is calculated as In addition, the fracture toughness expression for the inlet nozzle-vessel wall intersection steel, based on the ASME code IWB-3613, is given as e upper bound limiting criteria for the fracture toughness curve (equations (9) and (10)) were estimated, based on the API 579-1/ASME FFS-1 standard, as K ICv � 220 MPa �� m √ and K ICn � 155.6 MPa �� m √ , respectively. Figure 2 show the fracture toughness behaviour of the RPV steel material.

Reference Transient Event and Sequential Analysis
Implementation.
e AOO PTS event initiated by an inadvertent SIS operation selected for this study was investigated by Wang et al. [19]. is type of event frequently occurred in the lifetime of PWR NPPs as reported in IAEA operational experience feedback reports [14,26]. e unexpected action of the SIS caused the injection of cold water Science and Technology of Nuclear Installations into the RCS leading to a rise in the loop pressure. e maximum injection flow rate was approximately 13 kg/s. e cooling water temperature and internal pressure distributions obtained from the AOO event are shown in Figure 3. Also, the critical phase during the transient event was observed at the pressure value of 14.7 MPa at 642 s. e initial temperatures of the inner and outer surface of the RPV model are conservatively set to 289°C and 20°C, respectively. In addition, the heat transfer coefficient between the air and the outer surface vessel is assumed to be 20 W/m 2°C . e structural integrity evaluation was done following the sequential simulation process shown in Figure 4.

Heat Transfer Coefficient Estimation.
To evaluate the heat transferred from the injected cold water to the inner parts of the RPV model, an approximate heat transfer coefficient (h) was estimated using the following equation [27,28]: where N u , h, D, and k is the Nusselt number, heat transfer coefficient, hydraulic diameter (the diameter of the hot leg, 0.699 m, of the referenced PWR is assumed in the present FE model), and thermal conductivity of water, respectively. e Nusselt number is expressed as However, the Grashof and Prandtl numbers are determined from [27] Gr � D 3 ρ 2 gβΔT where β, C p , g, ΔT, ρ, and μ are the isobaric cubic expansion coefficient, specific heat capacity, gravity, temperature change, density, and viscosity, respectively [29,30]. e heat transfer coefficient was then derived as follows: C � 0.590, n � 0.3, and k � 0.6096 W/m K are assumed in computing the heat transfer coefficient. Also, the physical properties of the coolant water; β, C p , g, ρ, and μ are estimated at a PWR operational pressure of 14.7 MPa and ∆T � 289°C using STEAMEST software tool developed by [29].

Reactor Pressure Vessel Model.
e core beltline and connecting nozzle area are considered as part of the key regions of the RPV due to the exposure to high neutron irradiation that leads to an increase in nil-ductility transition temperature and embrittlement [31]. erefore, the threedimensional (3D) design of one-half of a double-circuit PWR RPV, without cladding, capturing a section of the beltline and inlet nozzle areas shown in Figure 5 was adopted in this study. e inner diameter of the RPV, vessel wall thickness, and nozzle's inner diameter are 3.374 m, 0.170 m, and 0.70 m, respectively.

ermal and ermomechanical Models.
e RPV FE models used for the sequential thermal-mechanical coupling analysis were preprocessed and meshed using the finely discretized hexahedral element type in HyperMesh software [32]. e thermal and thermomechanical FE models were meshed using element types DC3D20 and C3D20R, respectively, as shown in Figure 6(a). e bottom section of the FE model is constrained to limit the displacement of the rigid body, but the upper section is left in free state. Symmetricity boundary conditions are also applied. e heat transfer coefficient between the ambient air and the outer surface of the RPV was taken to be 20 W/m 2°C and outer vessel temperature as 20°C. e temperature and pressure histories obtained from the reference anticipated transient event, in Figure 5, and a heat expansion coefficient of 8450 w/m 2°C conservatively estimated from equation (15) were assumed as loads and input for the sequential coupling analysis. e loads and boundary conditions depicted in the RPV global model shown in Figure 6(b) were used to simulate the thermomechanical stresses in Abaqus software.

Fracture Mechanics Submodel.
e majority of failures associated with the RPV have been traced to shallow cracks. Axial surface cracks are also known to be more critical than circumferential surface cracks due to the maximum stress they experience in a vessel configuration [33]. Furthermore, probabilistic fracture mechanics studies reported in many literature studies show that shallow cracks contribute more to the probability of crack initiation in RPV than deep cracks [10]. is is because thermomechanical loading and neutron irradiation are more severe at the surface of the RPV [16]. In this study, an axially oriented semielliptical surface crack is assumed. Also, the crack geometry applied as shown in Figure 7(a) had a depth (a) of one-fourth the RPV wall thickness with an aspect ratio of (2c/a) � 6 [34]. ese Science and Technology of Nuclear Installations postulated crack size and aspect ratio are mostly observed during nondestructive testing or expected during the lifetime of the NPP [35]. e node-based submodels shown in Figures 8 and 9 which comprise the highest stress region obtained from the thermomechanical analysis ( Figure 10) were used for the K I calculation. e integrity parameter, K I , was evaluated based on the linear elastic fracture mechanics (LEFM) approach. e submodels incorporated all the boundary conditions of the RPV global model, shown in Figure 6(b), by interpolating the temporal displacements obtained from the global model as mechanical boundary conditions on the boundaries of the submodel during the simulation process. e stress intensity factor (K 1 ) was estimated as a function of the thermomechanical stress distribution, the assumed surface crack geometry, and the orientation using the principle of superposition. e K I values were estimated at the deepest crack point using both M-integral and XFEM approaches implemented in Abaqus-FRANC3D co-simulation and Abaqus software, respectively. e surface crack geometry shown in Figure 7(b) was used for the XFEM K I modeling process. e FE mesh surrounding the crack geometry in the submodel (Figures 8 and 9) was finely meshed so as to obtain accurate K I and also to alleviate oscillations associated with it using the XFEM method. e mathematical computation and   Highest stress location S, Mises (Avg: 75%) +6.400e+08 +5.868e+08 +5.336e+08 +4.804e+08 +4.272e+08 +3.739e+08 +3.207e+08 +2.675e+08 +2.143e+08 +1.079e+08 +1.611e+08 +5.464e+07 +1.420e+06 Science and Technology of Nuclear Installations application of the M-integral energy method in computing K I are comprehensively discussed in [36]. Also, details on the application of the XFEM method can be found in [37,38].

Sequential ermomechanical Coupling
Analysis. e sequential thermomechanical analysis using the FE models described in Section 3.4 was performed in Abaqus code following the multistep procedure shown in Figure 4. e thermal loads were determined separately and then coupled with the pressure loads to compute the thermomechanical stresses. is sequential coupling analysis was to determine among others the location of the highest stress from the numerical simulation results. Figure 10 shows the von Mises stress observed at 1600 s during the AOO event. e inlet nozzle and inner RPV wall intersection (depicted in Figure 10) experienced the maximal stress concentration. e high stress was due to the complex temperature field in that joint area. erefore, a crack-like defect was most likely to be initiated by a PTS loading at that location. e heat flux, temperature, and hoop and axial stress profiles through the thickness of inlet nozzle-inner RPV wall intersection (path), starting from the highest stress location, at different transient time periods are shown in Figures 11-14. Figures 15 and 16 also show the evolution of the temperature, hoop, and axial stresses extracted from the maximal stress element of interest in the FE model for the duration of the AOO event.
e heat flux decrease, observed through the measured thickness, was due to the steady decrease of the temperature gradient and thermal conductivity. Also, the temperature gradient observed in the course of transient time decreased steadily from the inlet nozzle-inner wall RPV intersection to the outer part of the RPV model. e temperature drop was because of the sudden cooling and RPV's large thermal capacity. Furthermore, the temperature varied in the range of 60°C to 250°C through the measured thickness. is variation was due to the relatively high velocity of the cooling water injected by the actuation of HPSI pumps. Also, temperature readings in this range during this AOO transient are rated as very high in the nuclear power field. e hoop and axial stress distributions also decreased from the same maximal location to compressive stress on the outer part of the RPV wall. e hoop stress ( Figure 16) exhibited a steady change throughout the transient time. Similarly, the axial stress increased steeply before remaining relatively stable from 100 s till the end of the AOO event. is phenomenon was largely due to the combined action of the thermal and internal pressure variations during the transient event.

Fracture Mechanics Analysis.
e integrity parameter, K I , was evaluated using a linear elastic fracture mechanics 3D-FE submodel with loadings obtained from the previous sequential thermomechanical coupling analysis (Section4.1). e mesh size and type of each submodel (Figures 8 and 9) were sensitively determined to ensure computational efficiency and accuracy of the stress intensity factors, K I . Also, the K I simulation followed the multistep coupling procedure shown in Figure 4. A semielliptical surface crack with a depth of a � 42 mm and an aspect ratio of (2c/a) � 6 was assumed for the analysis . e surface crack was inserted at the highest stress location in the RPV submodels obtained from the preceding thermomechanical analysis. e FE submodels shown in Figures 8 and 9 capture the maximal stress concentration area from the RPV global model ( Figure 10) and the postulated crack orientation. In addition, the thermomechanical loads and boundary conditions applied to the FE submodels were imposed from the 3D FE model shown in Figure 6(b). e K I values at the deepest tip of the crack were then computed using the M-integral approach coupled within the proposed sequential Abaqus-FRANC3D simulation method. To verify the K I values obtained from the proposed method, the conventional XFEM method was also applied to compute K I . To accurately estimate the ASME vessel material's fracture toughness, K IC , the random parameters of RPV's steel shown in Table 2, were used to calculate RT NDT (equation (8)). Figure 4 represents the comparison of the ASME curves generated based on RT NDT values inputted in equations (9) and (10). Figure 17 shows the temperature evolution at the deepest point of the surface crack during the AOO transient. It was observed that the crack tip temperature decreased steeply for the first 200 s before remaining relatively stable through the transient time of 1600 s. e slow rate of temperature change was due to the higher heat capacity of the RPV material compared with the overcooling water. e comparison of K I values along the postulated crack front and crack tip temperature calculated using M-integral and XFEM methods is presented in Figure 18, and the absolute error between the K 1 estimation of both methods is shown in Figure 19. e results presented in Figures 18 and  19 show that both methods are in good agreement irrespective of the varied assumptions and numerical methods. Also, the small K 1 values observed at the deepest crack tip was due to low compressive stresses. From the comparison of K I and K IC shown in Figure 20, the maximum K I of 43.01 MPam 0.5 was less than the upper bound fracture toughness of the vessel, K ICn . In summary, the semielliptical surface crack configuration investigation results show that although the pressurized thermal shock induced by inadvertent operation of the safety injection system may not pose an immediate safety risk to the structural integrity of the RPV, it may accelerate the damage of a highly embrittled RPV with time.

Conclusion
One of the most frequently anticipated operational transients in a pressurized water reactor is the inadvertent operation of the safety injection systems. is transient event can induce significant pressurized thermal shocks in the reactor pressure vessel; hence, it is vital to estimate the structural mechanical changes associated with such cyclic loadings. In this paper, a sequential Abaqus-FRANC3D simulation approach was followed in the determination of the potential structural integrity risk an inadvertent operation of SIS may pose to an ageing PWR RPV.
e key conclusions from the study are summarized as follows: (1) e applied sequential Abaqus-FRANC3D simulation method efficiently reduced the computational cost in the evaluation of the AOO PTS-induced stress intensity factors. is was due to the simplified LEFM submodeling technique and FE mesh type adopted in the simulation process.
(2) e integrity parameter, K I , calculated with M-integral and XFEM methods at the deepest crack tip was in good agreement with absolute errors of less than 5. erefore, the K I value obtained by the proposed sequential Abaqus-FRANC3D approach has a high accuracy and meets the precision criterion in structural engineering calculation. induced PTS may not pose an immediate integrity risk to the RPV, but it can cumulatively accelerate ageing mechanisms such as fatigue in the reactor pressure vessel.
is study is useful for both operators and regulators, for decision-making in reactor safety margin, ageing management, and service life extension. is work also serves as a foundation for our future study on the effect of AOO PTS loads on RPV cladding material, crack growth analysis, and fatigue life prediction.

Data Availability
Data are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.