Predicting Thermal Mixing and Fatigue Inside Control Rod Guide Tubes

Author(s):  
Eric Lillberg

The cracked control rods shafts found in two Swedish NPPs were subjected to thermal fatigue due to mixing of cold purge flow with hot bypass water in the upper part of the top tube on which the control rod guide tubes rests. The interaction between the jets formed at the bypass water inlets is the main source of oscillation resulting in low frequency downward motion of hot bypass water into the cold purge flow. This ultimately causes thermal fatigue in the control rod shaft in the region below the four lower bypass water inlets. The transient analyses shown in this report were done to further investigate this oscillating phenomenon and compare to experimental measurements of water temperatures inside the control rod guide tube. The simulated results show good agreement with experimental data regarding all important variables for the estimation of thermal fatigue such as peak-to-peak temperature range, frequency of oscillation and duration of the temperature peaks. The results presented in this report show that CFD using LES methodology and the open source toolbox OpenFOAM is a viable tool for predicting complex turbulent mixing flows and thermal loads.

Author(s):  
Hernan Tinoco ◽  
Hans Lindqvist ◽  
Ylva Odemark ◽  
Carl-Maikel Ho¨gstro¨m ◽  
Kristian Angele

Two broken control rods and a large number of rods with cracks were found at the inspection carried out during the refueling outage of the twin reactors Oskarshamn 3 and Forsmark 3 in the fall of 2008. As a part of an extensive damage investigation, time dependent CFD simulations of the flow and the heat transfer in the annular region formed by the guide tube and control rod stem were carried out, [1]. The simulations together with metallurgical and structural analyses indicated that the cracks were initiated by thermal fatigue. The knowledge assembled at this stage was sufficient to permit the restart of both reactors at the end of year 2008 conditioned to that further studies to be carried out for clarifying all remaining matters. Additionally, all control rods were inserted 14% to protect the welding region of the stem. Unfortunately, this measure led to new cracks a few months later. This matter will be explained in the second part of this work, [2]. As a part of the accomplished complementary work, new CFD models were developed in conformity with the guidelines of references [3] and [4]. The new results establish the simulation requirements needed to accomplish accurate conjugate heat transfer predictions. Those requirements are much more rigorous than the ones needed for flow simulations without heat transfer. In the present case, URANS simulations, which are less resource consuming than LES simulations, seem to rather accurately describe the mixing process occurring inside the control rod guide tube. Structure mechanics analyses based on the CFD simulations show that the cracks are initiated by thermal fatigue and that their propagation and growth are probably enhanced by mechanical vibrations.


Author(s):  
Hidemasa Yamano ◽  
Yoshiharu Tobita

This paper describes experimental analyses using SIMMER-III/IV, which are two/three-dimensional multi-component multi-phase Eulerian fluid-dynamics codes, for the purpose of the code validation. Two topics of key phenomena in core disruptive accidents were presented in this paper: duct-wall failure and fuel discharge/relocation behavior. To analyze the duct-wall failure behavior, the SCARABEE BE+3 in-pile experiments were selected. The SIMMER-III calculation was in good agreement with the overall event progression; which was characterized by coolant boiling, clad melting, fuel failure, molten pool formation, duct-wall failure, etc.; observed in the experiment. The CAMEL C6 experiment investigated the fuel discharge and relocation behavior through a simulated control rod guide tube, which is important in evaluating the neutronic reactivity. SIMMER-IV well simulated fuel-coolant interaction, sodium voiding, fuel relocation behavior observed in the experiment. These experimental analyses indicated the validity of the SIMMER-III/IV computer code for the duct wall failure and fuel discharge/relocation behavior.


Author(s):  
Kristian Angele ◽  
Mathias Cehlin ◽  
Carl-Maikel Ho¨gstro¨m ◽  
Ylva Odemark ◽  
Mats Henriksson ◽  
...  

A large number of control rod cracks were detected during the refuelling outage of the twin reactors Oskarshamn 3 and Forsmark 3 in the fall of 2008. The extensive damage investigation finally lead to the restart of both reactors at the end of 2008 under the condition that further studies would be conducted in order to clarify all remaining matters. Also, all control rods were inserted 14% in order to locate the welding region of the control rod stem away from the thermal mixing region of the flow. Unfortunately, this measure led to new cracks a few months later due to a combination of surface finish of the new stems and the changed flow conditions after the partial insertion of the control rods. The experimental evidence reported here shows an increase in the extension of the mixing region and in the intensity of the thermal fluctuations. As a part of the complementary work associated with the restart of the reactors, and to verify the CFD simulations, experimental work of the flow in the annular region formed by the guide tube and control rod stem was carried out. Two full-scale setups were developed, one in a Plexiglass model at atmospheric conditions (in order to be able to visualize the mixing process) and one in a steel model to allow for a higher temperature difference and heating of the control rod guide tube. The experimental results corroborate the general information obtained through CFD simulations, namely that the mixing region between the cold crud-removal flow and warm by-pass flow is perturbed by flow structures coming from above. The process is characterized by low frequent, high amplitude temperature fluctuations. The process is basically hydrodynamic, caused by the downward transport of flow structures originated at the upper bypass inlets. The damping thermal effects through buoyancy is of secondary importance, as also the scaling analysis shows, however a slight damping of the temperature fluctuations can be seen due to natural convection due to a pre-heating of the cold crud-removal flow. The comparison between numerical and experimental results shows a rather good agreement, indicating that experiments with plant conditions are not necessary since, through the existing scaling laws and CFD-calculations, the obtained results may be extrapolated to plant conditions. The problem of conjugate heat transfer has not yet been addressed experimentally since complex and difficult measurements of the heat transfer have to be carried out. This type of measurements constitutes one of the main challenges to be dealt with in the future work.


2017 ◽  
Vol 2017 ◽  
pp. 1-8 ◽  
Author(s):  
Tagor Malem Sembiring ◽  
Surian Pinem ◽  
Peng Hong Liem

The in-house coupled neutronic and thermal-hydraulic (N/T-H) code of BATAN (National Nuclear Energy Agency of Indonesia), NODAL3, based on the few-group neutron diffusion equation in 3-dimensional geometry using the polynomial nodal method, has been verified with static and transient PWR benchmark cases. This paper reports the verification of NODAL3 code in the NEA-NSC PWR uncontrolled control rods withdrawal at zero power benchmark. The objective of this paper is to determine the accuracy of NODAL3 code in solving the continuously slow and fast reactivity insertions due to single and group of control rod bank withdrawn while the power and temperature increment are limited by the Doppler coefficient. The benchmark is chosen since many organizations participated using various methods and approximations, so the calculation results of NODAL3 can be compared to other codes’ results. The calculated parameters are performed for the steady-state, transient core averaged, and transient hot pellet results. The influence of radial and axial nodes number was investigated for all cases. The results of NODAL3 code are in very good agreement with the reference solutions if the radial and axial nodes number is 2 × 2 and 2 × 18 (total axial layers), respectively.


2021 ◽  
Vol 247 ◽  
pp. 07017
Author(s):  
Andreas Pautz ◽  
Winfried Zwermann

Cold-startup and hot-standby reactivity accident tests conducted at the SPERT III E-core research reactor are analysed with the coupled neutron-kinetic/thermal-hydraulic code system DYN3D-ATHLET. Homogenised 2-group cross sections for DYN3D are thereby generated with the Monte Carlo neutron transport code Serpent 2 in combination with the ENDF/B-VII.1 cross section library. Results in terms of maximum power, energy release, and reactivity compensation are in good agreement with the experimental values. The time-dependent contributions to the reactivity feedback are investigated for both a cold-startup test and a hot-standby test. These findings prove the suitability of the combined application of the simulation codes to predict the reactor dynamic behaviour in the event of prompt-critical and super-prompt critical transients even for small reactor cores. Furthermore, static core characteristics of the SPERT III E-core reactor at cold-startup condition are analysed with using a static DYN3D model, a detailed Serpent reference model, and a simplified Serpent model consistent with the DYN3D model. The critical control rod position and the excess reactivities of both the control rods and the transient rod obtained with the Serpent reference model are consistent with the experimental values. For the same parameters, the DYN3D model is in good agreement with the Serpent simplified model.


2021 ◽  
Vol 11 (3) ◽  
pp. 1243
Author(s):  
Hongseok Jeong ◽  
Jeung-Hoon Lee ◽  
Yong-Hyun Kim ◽  
Hanshin Seol

The dominant underwater noise source of a ship is known to be propeller cavitation. Recently, attempts have been made to quantify the source strength using on-board pressure sensors near the propeller, as this has advantages over conventional noise measurement. In this study, a beamforming method was used to estimate the source strength of a cavitating propeller. The method was validated against a model-scale measurement in a cavitation tunnel, which showed good agreement between the measured and estimated source levels. The method was also applied to a full-scale measurement, in which the source level was measured using an external hydrophone array. The estimated source level using the hull pressure sensors showed good agreement with the measured one above 400 Hz, which shows potential for noise monitoring using on-board sensors. A parametric study was carried out to check the practicality of the method. From the results, it was shown that a sufficient recording time is required to obtain a consistent level at high frequencies. Changing the frequency resolution had little effect on the result, as long as enough data were provided for the one-third octave band conversion. The number of sensors affected the mid- to low-frequency data.


1988 ◽  
Vol 15 (3) ◽  
pp. 149-156 ◽  
Author(s):  
R. A. Cavina ◽  
N. E. Waters

The angular stiffness of a multiple looped span, subject to rotation about a vertical axis (torsion) and also to rotation about a horizontal or radial axis (mesio-distal tilt), have been derived using the complementary (strain) energy method. Experimental measurements on enlarged models were in good agreement with the values calculated from the theoretical relationships obtained. The variations in angular stiffness resulting from changes in the loop height, width, and position of clinical sized loops are discussed.


1992 ◽  
Vol 06 (06) ◽  
pp. 309-321 ◽  
Author(s):  
W.Y. CHING ◽  
MING-ZHU HUANG ◽  
YONG-NIAN XU ◽  
FANQI GAN

The electronic structure and optical properties of crystalline C 60 and their pressure dependence have been studied by first-principles local density calculations. It is shown that fcc C 60 has a low dielectric constant and an optical spectrum rich in structures. The spectrum shows five disconnected absorption bands in the 1.4 to 7.0 eV region with sharp structures in each band that can be attributed to critical point transitions. This is a manifestation of the localized molecular structure coupled with long range crystalline order unique to the C 60 crystal. At a sufficient high pressure, the structures in the optical spectrum start to merge due to the merging of the bands. These results are in good agreement with some recent experimental measurements.


Polymers ◽  
2021 ◽  
Vol 13 (11) ◽  
pp. 1703
Author(s):  
Michael Coja ◽  
Leif Kari

A waveguide model for a pre-compressed cylindrical natural rubber vibration isolator is developed within a wide frequency range—20 to 2000 Hz—and for a wide pre-compression domain—from vanishing to the maximum in service, that is 20%. The problems of simultaneously modeling the pre-compression and frequency dependence are solved by applying a transformation of the pre-compressed isolator into a globally equivalent linearized, homogeneous, and isotropic form, thereby reducing the original, mathematically arduous, and complex problem into a vastly simpler assignment while using a straightforward waveguide approach to satisfy the boundary conditions by mode-matching. A fractional standard linear solid is applied as the visco-elastic natural rubber model while using a Mittag–Leffler function as the stress relaxation function. The dynamic stiffness is found to depend strongly on the frequency and pre-compression. The former is resulting in resonance phenomena such as peaks and troughs, while the latter exhibits a low-frequency magnitude stiffness increase in addition to peak and trough shifts with increased pre-compressions. Good agreement with nonlinear finite element results is obtained for the considered frequency and pre-compression range in contrast to the results of standard waveguide approaches.


2001 ◽  
Vol 1 (2) ◽  
pp. 65-72 ◽  
Author(s):  
Ömer Akgiray ◽  
Ahmet M. Saatçı

A new approach to model media expansion during filter backwash is presented. The proposed approach is based on the assumption that the Ergun equation remains valid after fluidization. Mathematical formulas are derived for predicting expanded porosity for a given backwash velocity or backwash velocity for a given expanded porosity. These formulas can be easily used by the engineer. Values predicted using the proposed approach are in good agreement with experimental measurements.


Sign in / Sign up

Export Citation Format

Share Document