Seismic analysis of nuclear power plants by using three-dimensional finite element models: a review

ABSTRACT

Seismic design of nuclear power plants (NPPs) is important for ensuring their integrity during earthquakes. Seismic analysis has been conducted using lumped mass beam models (LMBMs) for the design of plants in Japan, whereas three-dimensional (3D) finite element models (FEMs) have been used for novel plants outside Japan. The purposes of this study are to organize issues related to the development and application of 3D FEMs for seismic analysis of Japanese NPPs and to indicate future study directions. To organize these issues, the authors systematically investigated: (1) international guides and standards related to seismic analysis and (2) 3D FEMs of novel NPPs outside Japan. By considering other studies on the issues, the authors suggest directions for future studies. Resolving the issues will contribute to application of 3D FEMs for seismic analysis in the design of Japanese NPPs.     

Seismic analysis of the APR 1400 reactor vessel internals using the model reduction method

To secure reliability of the seismic design of the reactor vessel internals (RVIs) through the finite element analysis, it is important to develop the accurate analysis model which can represent the geometric complexity of the RVIs. However, the seismic analysis requires too large computation cost to solve the complex equations; thus, it needs to reduce the overall size of the analysis model. Here, we apply a model reduction method based on the fixed-interface component mode synthesis (CMS) method to practical RVIs to solve complex numerical problems efficiently. To verify the model reduction method, several cases of the RVIs with different conditions are analyzed for the static and dynamic problems. Finally, the seismic analysis was performed with the suggested reduced model with the CMS method. The time history analysis is performed to extract important seismic responses at the specified locations, and the stress analysis is also performed to identify that the RVIs satisfy the seismic design. In the last part of the paper, an example of the design modification is suggested to reduce the stress intensity at the support locations.     

核电厂大型组合结构的有限元抗震分析方法研究

在现代核电站抗震设计中,有限元法是各类相关设备抗震分析与评价的重要数值仿真工具。对于形状复杂、部件众多的大型组合结构,采用整体三维建模的有限元模型通常需要很大的存储和计算规模,超出现有的计算条件。因此需要首先研究组合结构各个部件的动力学特性,从而建立合理的三维简化力学模型,并以该模型为基础进行有限元数值仿真。本文以某地车-吊车组合结构为例,给出此类大型组合结构的抗震分析方法,并将等效静力法与反应谱法相结合,对该结构进行分析,最后根据相关法规对各子结构进行评价,以确保总体组合结构在极限安全地震条件下能够保持结构完整性。     

3D finite element analysis of a nuclear fuel rod with gap elements between the pellet and the cladding

Nuclear fuel rods which comprises an important component of a nuclear power plant are composed of nuclear fuel and cladding. Simulating the nuclear fuel rod using a computer program is the universal method to verify its safety. The computer program used for this is called the fuel performance code. The main objective of this study is to simulate the nuclear fuel rod behavior considering the gap conductance using three-dimensional gap elements. Gap elements are used because, unlike other methods, this approach does not require special methods or other variables such as the Lagrange multiplier. In this work, a nuclear fuel rod has been simulated and the results are compared with the experimental results.     

Mass flow rate sensitivity and uncertainty analysis in natural circulation boiling water reactor core from Monte Carlo simulations

Our aim was to evaluate the sensitivity and uncertainty of mass flow rate in the core on the performance of natural circulation boiling water reactor (NCBWR). This analysis was carried out through Monte Carlo simulations of sizes up to 40,000, and the size, i.e., repetition of 25,000 was considered as valid for routine applications. A simplified boiling water reactor (SBWR) was used as an application example of Monte Carlo method. The numerical code to simulate the SBWR performance considers a one-dimensional thermo-hydraulics model along with non-equilibrium thermodynamics and non-homogeneous flow approximation, one-dimensional fuel rod heat transfer. The neutron processes were simulated with a point reactor kinetics model with six groups of delayed neutrons. The sensitivity was evaluated in terms of 99% confidence intervals of the mean to understand the range of mean values that may represent the entire statistical population of performance variables. The regression analysis with mass flow rate as the predictor variable showed statistically valid linear correlations for both neutron flux and fuel temperature and quadratic relationship for the void fraction. No statistically valid correlation was observed for the total heat flux as a function of the mass flow rate although heat flux at individual nodes was positively correlated with this variable. These correlations are useful for the study, analysis and design of any NCBWR. The uncertainties were propagated as follows: for 10% change in the mass flow rate in the core, the responses for neutron power, total heat flux, average fuel temperature and average void fraction changed by 8.74%, 7.77%, 2.74% and 0.58%, respectively.     

Analysis of the MOX critical experiments BASALA by pin-by-pin core analysis method using three-dimensional direct response matrix

An analysis of the MOX critical experiments BASALA was performed to verify the pin-by-pin core analysis method using a three-dimensional direct response matrix. The BASALA experiments simulate full MOX BWR cores, and they were carried out in the EOLE critical facility of the French Atomic Energy Commission (CEA) by the Nuclear Power Engineering Corporation (NUPEC) in collaboration with CEA. The BASALA experimental cores are very heterogeneous because their size is much smaller than that of commercial power plants. The main features of the pin-by-pin core analysis method using the three-dimensional direct response matrix are that the response matrix can reflect the intra-assembly heterogeneous effect, the diffusion approximation is not involved, and the fuel rod fission rate can be directly evaluated. The maximum difference of the critical k-effective values among all nine cores analyzed was about 0.4% Δk. The root mean square differences between the calculated and measured radial fuel rod fission rate distributions in the test assembly of all cores were within 1.8% and nearly comparable to measurement error. The calculated results of the reactivity worth agreed with the measured results within 9%. These good agreements mean that the pin-by-pin core analysis method using the three-dimensional direct response matrix accurately reflects the effects of the intra- and inter-assembly heterogeneities in heterogeneous systems like the BASALA experimental cores.     

Development of numerical analysis method of flow-acoustic resonance in stub pipes of safety relief valves

The boiling water reactors (BWRs) have steam dryer in the upper part of the pressure vessel to remove moisture from the steam. The steam dryer in the Quad Cities Unit 2 nuclear power plant was damaged by high-cycle fatigue due to acoustic-induced vibration during extended power uprate operation. The principal source of the acoustic-induced vibration was flow-acoustic resonance at the stub pipes of the safety relief valves (SRVs) in the main steam lines (MSLs). The acoustic wave generated at the SRV stub pipes propagates throughout the MSLs and eventually reaches and damages the steam dryer. Therefore, for power uprate operation of the BWRs, it has been required to predict the flow-acoustic resonance at the SRV stub pipes. The purpose of this article was to propose a numerical analysis method for evaluating the flow-acoustic resonance in the SRV stub pipes. The proposed method is based on the finite difference lattice Boltzmann method (FDLBM). So far, the FDLBM has been applied to flow-acoustic simulations of laminar flows around simple geometries at low Reynolds number. In order to apply the FDLBM to the flow-acoustic resonance simulations of turbulent flows around complicated geometries at the high Reynolds number, we developed computationally efficient model by introducing new function into the governing equation. The proposed method was compared with the conventional FDLBM in the cavity-driven flow simulation. The proposed method was validated by comparisons with the experimental data in the 1/10-scale test of BWR-5 under atmosphere condition. The following three results were obtained; the first is that the proposed method can reduce the computing time by 30% compared with the conventional FDLBM; the second is that the proposed method successfully simulated the flow-acoustic resonance in the SRV stub pipes of the BWR-5, and the pressure fluctuations of the simulation results agreed well with those of the experimental data; and the third is the mechanism of the flow-acoustic resonance in the SRV stub pipes. Acoustic waves causing the flow-acoustic resonance in the SRV stub pipes are generated by the unsteady vortices in the SRV stub pipes.     

Developing architecture for upgrading I&C systems of an operating nuclear power plant using a quality attribute-driven design method

This paper presents the architecture for upgrading the instrumentation and control (I&C) systems of a Korean standard nuclear power plant (KSNP) as an operating nuclear power plant. This paper uses the analysis results of KSNP's I&C systems performed in a previous study. This paper proposes a Preparation–Decision–Design–Assessment (PDDA) process that focuses on quality oriented development, as a cyclical process to develop the architecture. The PDDA was motivated from the practice of architecture-based development used in software engineering fields. In the preparation step of the PDDA, the architecture of digital-based I&C systems was setup for an architectural goal. Single failure criterion and determinism were setup for architectural drivers. In the decision step, defense-in-depth, diversity, redundancy, and independence were determined as architectural tactics to satisfy the single failure criterion, and sequential execution was determined as a tactic to satisfy the determinism. After determining the tactics, the primitive digital-based I&C architecture was determined. In the design step, 17 systems were selected from the KSNP's I&C systems for the upgrade and functionally grouped based on the primitive architecture. The overall architecture was developed to show the deployment of the systems. The detailed architecture of the safety systems was developed by applying a 2-out-of-3 voting logic, and the detailed architecture of the non-safety systems was developed by hot-standby redundancy. While developing the detailed architecture, three ways of signal transmission were determined with proper rationales: hardwire, datalink, and network. In the assessment step, the required network performance, considering the worst-case of data transmission was calculated: the datalink was required by 120 kbps, the safety network by 5 Mbps, and the non-safety network by 60 Mbps. The architecture covered 17 systems out of 22 KSNP's I&C systems. The architecture is implementable with the equipment developed in South Korea. The architecture can be used as a model to upgrade the existing I&C systems in a planned, large-scale, and one-shot manner. A more detailed architecture down to software level will be developed in the future.     

An overview of instrumentation and control systems of a Korea standard nuclear power plant: A signal interface standpoint

This paper presents an overview of instrumentation and control (I&C) systems of a pressurized water reactor (PWR) type nuclear power plant (NPP) in Korea. Yonggwang unit 3, which was constructed as a basis model for a Korea standard nuclear power plant (KSNP), is selected as an example for the presentation. This overview is derived from analyzing the I&C systems based on a top-down approach. The I&C systems consist of 30 systems. The 183 I&C cabinets are also analyzed and mapped to the systems. The overview is focused on an interface between the systems and the cabinets. This information will be used to understand the implementation of the I&C systems and to group the systems for an upgrade.