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

2021 ◽  
Vol 2021 ◽  
pp. 1-12
Author(s):  
M. Annor-Nyarko ◽  
Hong Xia

The 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 proposed method are compared with the conventional extended finite element method approach, and the result shows a good agreement. The maximal thermomechanical stress concentration was observed at the inlet nozzle-inner wall intersection. In addition, The 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. This work serves as a critical reference for the ageing management and fatigue life prediction of reactor pressure vessels.

2000 ◽  
Vol 199 (1-2) ◽  
pp. 101-111 ◽  
Author(s):  
S.N. Choi ◽  
K.S. Jang ◽  
J.S. Kim ◽  
J.B. Choi ◽  
Y.J. Kim

Author(s):  
Minghao Qin ◽  
Shu Tang ◽  
Francis Ku ◽  
Daniel V. Sommerville ◽  
Hal Gustin

T-stress is used as an indicator of the condition of crack tip constraint. In current fracture mechanics engineering applications in the U.S. nuclear industry, T-stress generally has been ignored during the calculation of applied stress intensity factors (SIF). Consideration of this crack tip constraint component could affect the evaluation of material fracture behavior, under either plane strain or plane stress or plane strain and plane stress combination. When the T-stress shows that the condition of crack front constraint is not plane strain, incorporation of T-stress may allow reduction of unnecessary conservatisms in such calculations. Under this condition, the allowable stress intensity factor is modified by increasing it above the KIc value, and it potentially increases the predicted allowable flaw sizes. In this paper, T-stress has been calculated using 3-D finite element analyses (FEA) with a typical semi-elliptical crack in a reactor pressure vessel (RPV) nozzle blend radius. Both thermal and internal pressure load cases are considered. To verify this finite element analysis approach, this method is applied to comparable literature models. The FEA results are consistent with closed-form solutions for T-stress calculation.


Author(s):  
Diego Fernando Mora Méndez ◽  
Markus Niffenegger ◽  
Guian Qian ◽  
Michal Jaros ◽  
Bojan Niceno

To perform the integrity assessment of a reactor pressure vessel (RPV) related to Pressurized Thermal Shock (PTS), we model the RPV using the 3D finite element method (FEM). Accurate prediction of temperature and stress fields is determined by using 2-Phase computational fluid dynamics (CFD) simulation in combination with an appropriate finite element discretization of the RPV wall. The cladding and the ferritic low alloy steel are considered as two separated layers, which can be intersected by superficial cracks. The calculation of the stress intensity factor (SIF) in mode I is based on the linear fracture mechanics theory (LEFM) and hypothetical cracks are located in different locations to consider the most critical cases. In the present study, the sub-modeling technique is implemented to refine the mesh required by the fracture analysis in the region of interest. Three types of cracks are considered: axial, circumferential and inclined. The performed integrity assessment uses the master curve approach. The stress intensity factor in the deepest point of a surface crack was compared with the material’s fracture toughness. In previous studies the integrity of the RPV subjected to medium and small break Loss-of-Coolant Accident (MBLOCA and SBLOCA, respectively) has been evaluated, therefore the concern in this contribution is the large break of Loss-Of-Coolant Accident (LBLOCA). The combination of 3D FEM with CFD simulations allows us to study the influence of the dynamic cooling plume on the stress intensity in more detail than with the conventional one dimensional method.


Author(s):  
Etienne de Rocquigny ◽  
Yoan Chevalier ◽  
Silvia Turato ◽  
Eric Meister

The structural integrity assessment of a nuclear Reactor Pressure Vessel (RPV) during accidental conditions such as loss-of-coolant accident (LOCA) is a major safety concern. Besides conventional deterministic calculations to justify as a nuclear operator the RPV integrity, Electricite´ de France (EDF) carries out probabilistic analyses. Probabilistic analyses become most interesting when some key variables, albeit conventionally taken at conservative values, can be modelled more accurately through statistical variability. In the context of low failure probabilities, this requires however a specific coupling effort between a specific probabilistic analysis method (e.g. Form-Sorm method) and the thermo-mechanical model to be reasonable in computing time. In this paper, the variability of a key variable — the mid-transient cooling temperature, tied to a climate-dependent tank — has been modelled, in some flaw configurations (axial sub-clad) for a French vessel. In a first step, a simplified analytical approach was carried out to assess its sensitivity upon the thermo-mechanical phenomena; hence, a direct coupling had to be implemented to allow a probabilistic calculation on the finite-element mechanical model, taking also into account a failure event properly defined through minimisation of the instantaneous failure margin during the transient. Comparison with the previous (indirectly-coupled) studies and the simplified analytical approach is drawn, demonstrating the interest of this new modelling effort to understand and order the sensitivity of the probability of crack initiation to the key variables. While being noticeable in the cases studied, sensitivity to the safety injection temperature variability proves to be less than the choice of the toughness model. Finally, regularity of the thermo-mechanical model is evidenced by the coupling exercise, suggesting that a modified response-surface based method could replace direct coupling for further investigation.


Author(s):  
Kiminobu Hojo ◽  
Naoki Ogawa ◽  
Yoichi Iwamoto ◽  
Kazutoshi Ohoto ◽  
Seiji Asada ◽  
...  

A reactor pressure vessel (RPV) head of PWR has penetration holes for the CRDM nozzles, which are connected with the vessel head by J-shaped welds. It is well-known that there is high residual stress field in vicinity of the J-shaped weld and this has potentiality of PWSCC degradation. For assuring stress integrity of welding part of the penetration nozzle of the RPV, it is necessary to evaluate precise residual stress and stress intensity factor based on the stress field. To calculate stress intensity factor K, the most acceptable procedure is numerical analysis, but the penetration nozzle is very complex structure and such a direct procedure takes a lot of time. This paper describes applicability of simplified K calculation method from handbooks by comparing with K values from finite element analysis, especially mentioning crack modeling. According to the verified K values in this paper, fatigue crack extension analysis and brittle fracture evaluation by operation load were performed for initial crack due to PWSCC and finally structural integrity of the penetration nozzle of RPV head was confirmed.


Author(s):  
B. Richard Bass ◽  
Paul T. Williams ◽  
Terry L. Dickson ◽  
Hilda B. Klasky

This paper describes further results from an ongoing study of a simplified engineering model that is intended to account for effects of clad residual stresses on the propensity for initiation of preexisting inner-surface flaws in a commercial nuclear reactor pressure vessel (RPV). The deposition of stainless steel cladding during fabrication of an RPV generates residual stresses in the cladding and the heat affected zone of the under-lying base metal. In addition to residual stress, thermal strains are generated by the differential thermal expansion (DTE) of the cladding and base material due to temperature changes during normal operation. A simplified model used in the ORNL-developed FAVOR probabilistic fracture mechanics (PFM) code accounts for the clad residual stress by incorporating a stress-free temperature (SFT) approach. At the stress-free temperature (Ts-free), the model assumes there is no thermal strain, i.e., the thermal expansion stresses and clad residual stresses offset each other. For normal cool-down transients applied to the RPV, interactions of the latter stresses generate additional crack driving forces on shallow, internal surface-breaking flaws near the clad/base metal interface; those flaws tend to dominate the RPV failure probability computed by FAVOR. In a previous report from this study (PVP2015-45086), finite element analysis was used to compare the stresses and stress-intensity factors (SIF) during a cool-down transient for two cases: (1) the existing SFT model of FAVOR, and (2) directly applied RPV clad residual stress (CRS) distribution obtained from empirical (hole-drilling) measurements made at room temperature on an RPV that was never put into service. However, those analyses were limited in scope and focused on a single flaw orientation. In this updated study, effects of CRS on the SIF histories computed for both circumferential and axial flaw orientations subjected to a cool-down transient were determined from an extended set of finite element analyses. Specifically, comparisons were made between results from applying CRS experimental data to ABAQUS two-dimensional, inner-surface flaw models and those generated by the FAVOR SFT model. It is demonstrated that the FAVOR-recommended SFT value of 488 °F produces conservatively high values of SIF relative to the use of CRS profiles in the ABAQUS models. For the vessel and flaw geometry and transient under study, the circumferential flaw (360° continuous) required a decrease of SFT down to 390 °F to match the CRS SIF histories. For the infinite axial flaw model, a decrease down to 300 °F matched the CRS SIF histories. Future plans are described to develop more general conclusions regarding the FAVOR model.


Sign in / Sign up

Export Citation Format

Share Document