• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.243-256
• /
• 2011

This paper considers the concept of analog computing based on a cellular neural network (CNN) paradigm to simulate nuclear reactor dynamics using a time-dependent second order form of the neutron transport equation. Instead of solving nuclear reactor dynamic equations numerically, which is time-consuming and suffers from such weaknesses as vulnerability to transient phenomena, accumulation of round-off errors and floating-point overflows, use is made of a new method based on a cellular neural network. The state-of-the-art shows the CNN as being an alternative solution to the conventional numerical computation method. Indeed CNN is an analog computing paradigm that performs ultra-fast calculations and provides accurate results. In this study use is made of the CNN model to simulate the space-time response of scalar flux distribution in steady state and transient conditions. The CNN model also is used to simulate step perturbation in the core. The accuracy and capability of the CNN model are examined in 2D Cartesian geometry for two fixed source problems, a mini-BWR assembly, and a TWIGL Seed/Blanket problem. We also use the CNN model concurrently for a typical small PWR assembly to simulate the effect of temperature feedback, poisons, and control rods on the scalar flux distribution.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1209-1217
• /
• 2019

Due to the limitation of the computers, the conventional homogenization method is based on many assumptions and approximations, and some tough problems such as energy spectrum and boundary condition are faced. To deal with those problems, the Monte Carlo global homogenization is adopted. The Reactor Monte Carlo code RMC is used to study the global homogenization method, and the one-step global homogenization method is proposed. The superimposed mesh geometry is also used to divide the physical models, leading to better geometric flexibility. A set of multigroup homogenization cross sections is online generated for each mesh under the real neutron energy spectrum and boundary condition, the cross sections are adjusted by the superhomogenization method, and no leakage correction is required. During the process of superhomogenization, the author-developed reactor core program NLSP3 is used for global calculation, so the global flux distribution and equivalent homogenization cross sections could be solved simultaneously. Meanwhile, the calculated homogenization cross section could accurately reconstruct the non-homogenization flux distribution and could also be used for fine calculation. This one-step global homogenization method was tested by a PWR assembly and a small reactor model, and the results show the validity.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.152-161
• /
• 2022

We analyzed mitigation strategies for steam generator tube rupture (SGTR) accidents using MARS code under both full-power and low-power and shutdown (LPSD) conditions. In general, there are two approaches to mitigating SGTR accidents: supplementing the reactor coolant inventory using safety injection systems and depressurizing the reactor coolant system (RCS) by cooling it down using the intact steam generator. These mitigation strategies were compared from the viewpoint of break flow from the ruptured steam generator tube, the core integrity, and the possibility of the main steam safety valves opening, which is associated with the potential release of radiation. The "cooldown strategy" is recommended for break flow control, whereas the "RCS make-up strategy" is better for RCS inventory control. Under full power, neither mitigation strategy made a significant difference except for on the break flow while, in LPSD modes, the RCS cooldown strategy resulted in lower break and discharge flows, and thus less radiation release. As a result, using the cooldown strategy for an SGTR under LPSD conditions is recommended. These results can be used as a fundamental guide for mitigation strategies for SGTR accidents according to the operational mode.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2714-2719
• /
• 2022

The present paper proposes some new computational methods and results in the framework of flow computation through congested domains seen as porous media, as it can be found in the core of a Pressurized Water Reactor (PWR). The flow is thus mostly governed by the distribution of pressure losses, both through the porous structures, such as fuel assemblies, and in the thin fluid layers between them. The purpose of the present paper is to consider the question of the interaction of a flow and a rod bundle from an analytical point of view gathering all the contributions through a set of equations as simple and representative as possible. It aims at demonstrating a sound understanding of the relevant phenomena governing the flow establishment in the geometry of interest instead of relying mainly on a posteriori observations obtained both experimentally and numerically. Comparison with two set of experimental results showed good agreement. The model proposed being analytical it appears easily implementable for studies needing an expression of fluid forces in a rod array as for fuel assembly bowing issue. It would be interesting to test the reliability of the model on other geometry with different P/R ratios.

• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4601-4619
• /
• 2022

The classical Probabilistic Safety Analysis (PSA) does not include any time dependence explicitly. However, the success criteria (SC) could evolve during the cycle for some initiating events. In that sense, there is a type of sequence in which this time-dependency is quite important, the family of Anticipated Transient without Scram (ATWS) sequences in Pressurized Water Reactors. Therefore, a new risk-informed approach is proposed in this paper, which makes it possible to obtain the time-dependent SC evolution of the safety functions affected by the Moderator Temperature Coefficient (MTC) value. Then, the evolution of the ATWS conditional core damage probability (CCDP) could be obtained using a PSA model. To quantify the CCDP, the average values of the time-dependent failure probabilities must be computed. Finally, the comparison between the CCDP obtained through the application of the classical PSA approach and the new one makes it possible to quantify the impact of time-dependence on the SC of the headers that this new risk-informed ATWS approach can provide.

• Nuclear Engineering and Technology
• /
• v.56 no.6
• /
• pp.2163-2173
• /
• 2024

This study introduces a transient analysis module developed for RAST-V and validates it using the Kalinin-3 benchmark problem. For the benchmark analysis, RAST-V standalone and STREAM/RAST-V calculations were performed. STREAM supplies the few-group constants and RAST-V conducts a 3D core simulation utilizing few-group cross-sectional data. To improve accuracy, the main solver was developed based on the advanced semi-analytic nodal method. To evaluate the computational capability of the transient analysis module in RAST-V, Kalinin-3 benchmark is employed. Kalinin-3 represents a coolant transient benchmark that offers experimental data during the deactivation of the Main Circulation Pumps. Consequently, the transient calculations reflected the changes in the reactor flow rate. Benchmark comprising steady-state and transient calculations. During the steady state, the STREAM/RAST-V combination demonstrated a 30 ppm root mean square difference from 0 to 128.50 EFPD. For the transient calculations, STREAM/RAST-V showed power differences within ±7 % over a range of 0-300 s. Axial offset differences were within ±3 %, and the RMS difference in radial power ranged within 2.596 % at both 0 and 300 s. Overall, this study effectively demonstrated the newly developed transient solver in RAST-V and validated it using the Kalinin-3 benchmark problem.

• Nuclear Engineering and Technology
• /
• v.28 no.2
• /
• pp.137-149
• /
• 1996

The present study is to analyze an integral test, BETHSY test 6.9c, which represent loss of RURS accident under mid-loop operation. Both the pressurizer manway and the steam generator outlet plenum manway are opened as vent paths in order to prevent the system from pressurization by removing the steam generated in the core. The main purposes are to gain insights into the physical phenomena and identify sensitive parameters. Assessment of capability of CATHARE2 prediction can be established the effective recovery procedures using the code in an actual plant. Most of important physical phenomena in the experiment could be predicted by the CATHARE2 code. The peak pressure in the upper plenum is predicted higher than experimental value by 7 kPa since the differential pressure between the pressurizer and the surge line is overestimated. The timing of core uncovery is delayed by 500 seconds mainly due to discrepancy in the core void distribution. It is demonstrated that openings of the pressurizer manwey and the steam generator manway can prevent the core uncovery using only gravity feed injection. Although some disagreements are found in the detailed phenomena, the code prediction is considered reasonable for the overall system behaviors.

• Nuclear Engineering and Technology
• /
• v.18 no.3
• /
• pp.209-217
• /
• 1986

The inlet coolant temperature reduction technique has been used to extend the load follow operation further in the end-of-cycle-life(EOL). In order to simulate the technique and calculate the nuclear characteristics of a PWR core according to the load follow operation, the three dimensional computing system has been established. The analysis was performed in both MINB and SPINR modes of typical 12-3-6-3 locad follow operation for the EOL of KNU-1 plant. Moreover, the capability of return-to-power has been also tested for those two modes with the system analysis by the RETRAN-02 code. The results show that it has no difficulty to extend the load follow operation further in the EOL by applying the inlet coolant temprature reduction, and also the spinning reserve capacity(SRC) increases by 13% in MINB mode and 14% in SPINR mode Bore that used by control rods only, for 14$^{\circ}$ F drop in the inlet temperature.

• Nuclear Engineering and Technology
• /
• v.23 no.3
• /
• pp.352-362
• /
• 1991

Recent popular trends in pressurized water reactor(PWR) fuel management are to extend the cycle length and to employ the low-leakage core designs for the optimal utilization of the uranium resources. In control strategy incorporated with the fuel management, turnable absorbers are required to control the power peaking and to ensure a negative moderator temperature coefficient during reactor operation. In this study, the nuclear characteristics and the optimal allocation of gadolinium-poisoned rods within the fuel assembly are considered using KWU SAV 79 A Code Package. First, analyses are carried out to compare the nuclear characteristics of the fuel assemblies contain-ing WABA(Wet Annular Burnable Absorber) and Gadolinium burnable absorbers respectively. The analyses show that the gadolinium-bearing fuel assembly has peculiar depletion characteristics ensuing from the very large thermal neutron absorption cross section. Peculiar characteristics of gadolinium provide basis for the optimal allocation of Gd rods in fuel assembly. Second, the methodology of an optimal allocation of gadolinium-poisoned rods within the fuel assembly is developed and applied to some nuclear designs.

• Nuclear Engineering and Technology
• /
• v.42 no.5
• /
• pp.474-499
• /
• 2010

This paper presents the mostortant developments implemented in the APOLLO2 spectral code since its last general presentation at the 1999 M&C conference in Madrid. APOLLO2 has been provided with new capabilities in the domain of cross section self-shielding, including mixture effects and transfer matrix self-shielding, new or improved flux solvers (CPM for RZ geometry, heterogeneous cells for short MOC and the linear-surface scheme for long MOC), improved acceleration techniques ($DP_1$), that are also applied to thermal and external iterations, and a number of sophisticated modules and tools to help user calculations. The method of characteristics, which took over the collision probability method as the main flux solver of the code, allows for whole core two-dimensional heterogeneous calculations. A flux reconstruction technique leads to fast albeit accurate solutions used for industrial applications. The APOLLO2 code has been integrated (APOLLO2-A) within the $ARCADIA^{(R)}$ reactor code system of AREVA as cross section generator for PWR and BWR fuel assemblies. APOLLO2 is also extensively used by Electricite de France within its reactor calculation chain. A number of numerical examples are presented to illustrate APOLLO2 accuracy by comparison to Monte Carlo reference calculations. Results of the validation program are compared to the measured values on power plants and critical experiments.