• Nuclear Engineering and Technology
• /
• v.45 no.6
• /
• pp.767-790
• /
• 2013

MAAP4 is a computer code that can simulate the response of a light water reactor power plant during severe accident sequences, including actions taken as part of accident management. The code quantitatively predicts the evolution of a severe accident starting from full power conditions given a set of system faults and initiating events through events such as core melt, reactor vessel failure, and containment failure. Furthermore, models are included in the code to represent the actions that could mitigate the accident by in-vessel cooling, external cooling of the reactor pressure vessel, or cooling the debris in containment. A key element tied to using a code like MAAP4 is an uncertainty analysis. The purpose of this paper is to present a MAAP4 based analysis to examine the sensitivity of a key parameter, in this case hydrogen production, to a set of model parameters that are related to a Level 2 PRA analysis. The Level 2 analysis examines those sequences that result in core melting and subsequent reactor pressure vessel failure and its impact on the containment. This paper identifies individual contributors and MAAP4 model parameters that statistically influence hydrogen production. Hydrogen generation was chosen because of its direct relationship to oxidation. With greater oxidation, more heat is added to the core region and relocation (core slump) should occur faster. This, in theory, would lead to shorter failure times and subsequent "hotter" debris pool on the containment floor.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.369-376
• /
• 2019

A complete solution for a soluble-boron-free (SBF) small modular reactor (SMR) is pursued with a new burnable absorber concept, namely centrally-shielded burnable absorber (CSBA). Neutronic flexibility of the CSBA design has been discussed with fuel assembly (FA) analyses. Major design parameters and goals of the SBF SMR are discussed in view of the reactor core design and three CSBA designs are introduced to achieve both a very low burnup reactivity swing (BRS) and minimal residual reactivity of the CSBA. It is demonstrated that the core achieves a long cycle length (~37 months) and high burnup (~30 GWd/tU), while the BRS is only about 1100 pcm and the radial power distribution is rather flat. This research also introduces a supplementary reactivity control mechanism using stainless steel as mechanical shim (MS) rod to obtain the criticality during normal operation. A further analysis is performed to investigate the local power peaking of the CSBA-loaded FA at MS-rodded condition. Moreover, a simple $B_4C$-based control rod arrangement is proposed to assure a sufficient shutdown margin even at the cold-zero-power condition. All calculations in this neutronic-thermal hydraulic coupled investigation of the 3D SBF SMR core are completed by a two-step Monte Carlo-diffusion hybrid methodology.

• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1813-1824
• /
• 2022

This paper shows the verification work of SARAX code system in the reactor core transient calculation based on the simplified EBR-II Benchmark. The SARAX code system is an analysis package developed by Xi'an Jiaotong University and aims at the advanced reactor R&D. In this work, a neutron-photon coupled power calculation model and a spatial-dependent reactivity feedback model were introduced. To verify the models used in SARAX, the EBR-II SHRT-45R test was simplified to an ULOF transient with an input flowrate change curve by fitting from reference. With the neutron-photon coupled power calculation model, SARAX gave close results in both power fraction and peak power prediction to the reference results. The location of the hottest assembly from SARAX and reference are the same and the relative power deviation of the hottest assembly is 2.6%. As for transient analysis, compared with experimental results and other calculated results, SARAX presents coincident results both in trend and absolute value. The minimum value of core net reactivity during the transient agreed well with the reported results, which ranged from -0.3$ to -0.35$. The results verify the models in SARAX, which are correct and able to simulate the in-core transient with reliable accuracy.

• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.980-992
• /
• 2024

Steam line break accident (SLB) in the nuclear reactor is one of the representative Non-LOCA accidents in which thermal-hydraulics and neutron kinetics are strongly coupled each other. Thus, the multi-scale and multi-physics approach is applied in this study in order to examine a realistic safety margin. An entire reactor coolant system is modelled by system scale node, whereas sub-channel scale resolution is applied for the region of interest such as the reactor core. Fuel performance code is extended to consider full core pin-wise fuel behaviour. The MARU platform is developed for easy integration of the codes to be coupled. An initial stage of the steam line break accident is simulated on the MARU platform. As cold coolant is injected from the cold leg into the reactor pressure vessel, the power increases due to the moderator feedback. Three-dimensional coolant and fuel behaviour are qualitatively visualized for easy comprehension. Moreover, quantitative investigation is added by focusing on the enhancement of safety margin by means of comparing the minimum departure from nucleate boiling ratio (MDNBR). Three factors contributing to the increase of the MDNBR are proposed: Various geometric parameters, realistic power distribution by neutron kinetics code, Radial coolant mixing including sub-channel physics model.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.641-654
• /
• 2014

A new computer code, named CORONA (Core Reliable Optimization and thermo-fluid Network Analysis), was developed for the core thermo-fluid analysis of a prismatic gas cooled reactor. The CORONA code is targeted for whole-core thermo-fluid analysis of a prismatic gas cooled reactor, with fast computation and reasonable accuracy. In order to achieve this target, the development of CORONA focused on (1) an efficient numerical method, (2) efficient grid generation, and (3) parallel computation. The key idea for the efficient numerical method of CORONA is to solve a three-dimensional solid heat conduction equation combined with one-dimensional fluid flow network equations. The typical difficulties in generating computational grids for a whole core analysis were overcome by using a basic unit cell concept. A fast calculation was finally achieved by a block-wise parallel computation method. The objective of the present paper is to summarize the motivation and strategy, numerical approaches, verification and validation, parallel computation, and perspective of the CORONA code.

• Nuclear Engineering and Technology
• /
• v.48 no.2
• /
• pp.318-329
• /
• 2016

Reducing critical boron concentration in a commercial pressurized water reactor core offers many advantages in view of safety and economics. This paper presents a preliminary investigation of a reduced-boron pressurized water reactor core to achieve a clearly negative moderator temperature coefficient at hot zero power using the newly-proposed "Burnable absorber-Integrated Guide Thimble" (BigT) absorbers. The reference core is based on a commercial OPR1000 equilibrium configuration. The reduced-boron ORP1000 configuration was determined by simply replacing commercial gadolinia-based burnable absorbers with the optimized BigT-loaded design. The equilibrium cores in this study were directly searched via repetitive Monte Carlo depletion calculations until convergence. The results demonstrate that, with the same fuel management scheme as in the reference core, application of the BigT absorbers can effectively reduce the critical boron concentration at the beginning of cycle by about 65 ppm. More crucially, the analyses indicate promising potential of the reduced-boron OPR1000 core with the BigT absorbers, as its moderator temperature coefficient at the beginning of cycle is clearly more negative and all other vital neutronic parameters are within practical safety limits. All simulations were completed using the Monte Carlo Serpent code with the ENDF/B-VII.0 library.

• Nuclear Engineering and Technology
• /
• v.36 no.5
• /
• pp.420-429
• /
• 2004

The local power density should be estimated accurately to prevent fuel rod melting. The local power density at the hottest part of a hot fuel rod, which is described by the power peaking factor, is more important information than the local power density at any other position in a reactor core. Therefore, in this work, the power peaking factor, which is defined as the highest local power density to the average power density in a reactor core, is estimated by fuzzy neural networks using numerous measured signals of the reactor coolant system. The fuzzy neural networks are trained using a training data set and are verified with another test data set. They are then applied to the first fuel cycle of Yonggwang nuclear power plant unit 3. The estimation accuracy of the power peaking factor is 0.45% based on the relative $2_{\sigma}$ error by using the fuzzy neural networks without the in-core neutron flux sensors signals input. A value of 0.23% is obtained with the in-core neutron flux sensors signals, which is sufficiently accurate for use in local power density monitoring.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3990-4002
• /
• 2021

CSPACE (Core meltdown, Safety and Performance Analysis CodE for nuclear power plants) for a simulation of severe accident progression in a Pressurized Water Reactor (PWR) is developed by coupling of verified system thermal hydraulic code of SPACE (Safety and Performance Analysis CodE for nuclear power plants) and core damage progression code of COMPASS (Core Meltdown Progression Accident Simulation Software). SPACE is responsible for the description of fluid state in nuclear system nodes, while COMPASS is responsible for the prediction of thermal and mechanical responses of core fuels and reactor vessel heat structures. New heat transfer models to each phase of the fluid, flow blockage, corium behavior in the lower head are added to COMPASS. Then, an interface module for the data transfer between two codes was developed to enable coupling. An implicit coupling scheme of wall heat transfer was applied to prevent fluid temperature oscillation. To validate the performance of newly developed code CSPACE, we analyzed typical severe accident scenarios for OPR1000 (Optimized Power Reactor 1000), which were initiated from large break loss of coolant accident, small break loss of coolant accident, and station black out accident. The results including thermal hydraulic behavior of RCS, core damage progression, hydrogen generation, corium behavior in the lower head, reactor vessel failure were reasonable and consistent. We demonstrate that CSPACE provides a good platform for the prediction of severe accident progression by detailed review of analysis results and a qualitative comparison with the results of previous MELCOR analysis.

• Nuclear Engineering and Technology
• /
• v.45 no.5
• /
• pp.573-580
• /
• 2013

The combined effects of reactivity coefficients, along with other core nuclear characteristics, determine reactor core behavior in normal operation and accident conditions. The Power Coefficient of Reactivity (PCR) is an aggregate indicator representing the change in reactor core reactivity per unit change in reactor power. It is an integral quantity which captures the contributions of the fuel temperature, coolant void, and coolant temperature reactivity feedbacks. All nuclear reactor designs provide a balance between their inherent nuclear characteristics and the engineered reactivity control features, to ensure that changes in reactivity under all operating conditions are maintained within a safe range. The $CANDU^{(R)}$ reactor design takes advantage of its inherent nuclear characteristics, namely a small magnitude of reactivity coefficients, minimal excess reactivity, and very long prompt neutron lifetime, to mitigate the demand on the engineered systems for controlling reactivity and responding to accidents. In particular, CANDU reactors have always taken advantage of the small value of the PCR associated with their design characteristics, such that the overall design and safety characteristics of the reactor are not sensitive to the value of the PCR. For other reactor design concepts a PCR which is both large and negative is an important aspect in the design of their engineered systems for controlling reactivity. It will be demonstrated that during Loss of Regulation Control (LORC) and Large Break Loss of Coolant Accident (LBLOCA) events, the impact of variations in power coefficient, including a hypothesized larger than estimated PCR, has no safety-significance for CANDU reactor design. Since the CANDU 6 PCR is small, variations in the range of values for PCR on the performance or safety of the reactor are not significant.

• Nuclear Engineering and Technology
• /
• v.30 no.3
• /
• pp.212-221
• /
• 1998

In the present study, a method for the dynamic analysis of a reactor core is developed. Peak responses for the motions induced from earthquake are obtained for a core model. The dynamic responses such as fuel assembly shear force, bending moment, axial force and displacement, and spacer grid impact loads are investigated. Prediction of fuel assembly stress during an earthquake requires development of a fuel assembly stress analysis model capable of interfacing with the models and results discussed in the dynamic analysis of a reactor core. This analysis uses beam characteristics which describe the overall fuel assembly response. The stress analysis method and its application for the case of an increased seismic level are also presented.