Simulation tests for temperature mixing in a core bottom model of the HTR-Module

Interatom and Siemens are developing a helium-cooled Modular High Temperature Reactor. Under nominal operating conditions temperature differences of up to 120°C will occur in the 700°C hot helium flow leaving the core. In addition, cold gas leakages into the hot gas header can produce even higher temperature differences in the coolant flow. At the outlet of the reactor only a very low temperature difference of maximum ±15°C is allowed in order to avoid damages at the heat exchanging components due to alternating thermal loads. Since it is not possible to calculate the complex flow behaviour, experimental investigations of the temperature mixing in the core bottom had to be carried out in order to guarantee the necessary reduction of temperature differences in the helium. The presented air simulation tests in a 1:2.9 scaled plexiglass model of the core bottom showed an extremely high mixing rate of the hot gas header and the hot gas duct of the reactor. The temperature mixing of the simulated coolant flow as well as the leakage flows was larger than 95%. Transfered to reactor conditions this means a temperature difference of only ±3°C for the main flow at a quite reasonable pressure drop. For the cold gas leakages temperature differences in the hot gas up to 400°C proved to be permissible. The results of the simulation experiments in the Aerodynamic Test Facility of Interatom permitted to design a shorter bottom reflector of the core.     

Numerical investigations of thermal mixing performance of a hot gas mixing structure in high-temperature gas-cooled reactor

A numerical simulation study was performed to clarify the thermal mixing characteristics of coolant in the core bottom structure of the high-temperature gas-cooled reactor(HTR). The flow field and temperature field in the hot gas chamber and the hot gas duct of the HTR were obtained based on the commercial computational fluid dynamics(CFD) program. The numerical simulation results showed that the helium flow with different temperatures in the hot gas mixing chamber and the hot gas duct mixed intensively, and the mixing rate of the temperature in the outlet of the hot gas duct reached 98 %. The results indicated many large-scale swirling flow structures and strong turbulence in the hot gas mixing chamber and the entrance of the hot gas duct, which were responsible for the excellent thermal mixing of the hot gas chamber and the hot gas duct. The calculated results showed that the temperature mixing rate of the hot gas chamber decreased only marginally with increasing Reynolds number.     

CFD investigations of helical wire-wrap fuel pin bundle and its comparison with straight wire bundle

Computational Fluid Dynamics (CFD) investigations of a fast reactor fuel pin bundle wrapped with helical and straight spacer wires have been carried out and the advantages of using helical spacer wire have been assessed. The flow and temperature distributions in the fuel pin bundle are obtained by solving the statistically averaged 3-Dimensional conservation equations of mass, momentum and energy along with high Reynolds number k-ε turbulence model using a customized CFD code CFDEXPERT. It is seen that due to the helical wire-wrap spacer, the coolant sodium not only flows in axial direction in the fuel pin bundle but also in a transverse direction. This transverse flow enhances mixing of coolant among the sub channels and due to this, the friction factor and heat transfer coefficient of the coolant increase. Estimation of friction factor, Nusselt number, sodium temperature uniformity at the outlet of the bundle and clad hot spot factor which are measures of the extent of coolant mixing and non-homogeneity in heat transfer coefficient around fuel pin are paid critical attention. It is seen that the friction factor and Nusselt number are higher (by 25% and 15% respectively) for the helical wire wrap pin bundle compared to straight wire bundle. It is seen that for 217 fuel pin bundle the maximum clad temperature is 750 K for straight wire case and the same for helical wire is 720 K due to the presence of transverse flow. The maximum temperature occurs at the location of the gap between pin and wire. The ΔT between the bulk sodium in the central sub-channel and peripheral sub-channel is 30 K for straight wire and the same for helical wire is 18 K due to the presence of secondary transverse flow which makes the outlet temperature more uniform. The hotspot factor and the hot channel factors predicted by CFD simulation are 10% lower than that used in conventional safety analysis indicating the conservatism in the safety analysis.     

高温气冷实验堆堆芯出口热气联箱混合性能对比试验

１０ＭＷ高温气冷实验堆堆芯出口冷却剂温度径向分布很不均匀，若不使之均匀化，将造成蒸汽发生器件部上过大的热应力，设置在堆底反射层中的堆芯出口热气联箱的作用之一是使冷却剂氦气在其中得到充分的热混合。     

Uncontrolled withdrawal of a control rod without scram

In the present work the thermal-hydraulics of reactivity-induced transients in low enriched uranium (LEU) core of a typical material test research reactor (MTR) are analyzed using the previous program developed by Khater et al. The analysis was done for uncontrolled withdrawal of a control rod with scram-disabled conditions. Initiating reactivity events with and without the influence of reactivity efficiency curve (“S” curve) were considered. The results of the proposed transients are analyzed and compared with each other. In transient without the “S” curve influence, a high primary peak power of 406.18 MW is attained and a clad melt down takes place after 1.85 s. In the transient with the “S” curve influence, a high super prompt-critical situation is produced (1.762$ at 0.895 s) with a very high primary peak power of 801.05 MW at 0.912 s. Also, a fast clad melt down is resulted in the hot channel at 1.088 s and a stable film boiling is established. This study indicates that, compared to the application of linear reactivity curve, the application of the reactivity efficiency curve results in the prediction of higher peaks in power and temperatures (fuel, clad and coolant) with a fast clad melt down.     

Analysis, verification, and benchmarking of the transient thermal hydraulic ITERTHA code for the design of ITER divertor

An analytical study for the International Thermonuclear Experimental Reactor Thermal Hydraulic Analysis code (ITERTHA) is carried out for a copper divertor with a 5 mm tungsten tile. The influence of the incident heat flux, swirl-tape insertion in cooling channels as well as the coolant flow velocity on the divertor thermal response is analyzed and discussed. The ITERTHA code results are verified by the commercial finite element code, COSMOS. The heat transfer coefficients at the nodes located on the cooling channel-wall are determined outside COSMOS code by the same methodology used in ITERTHA. A good agreement is achieved under different incident heat fluxes. The ITERTHA code is also benchmarked against the thermal-hydraulic calculation of the outer divertor of the Fusion Ignition Research Experiment, FIRE for an incident heat flux of 20 MW/m² and coolant flow velocity of 10 m/s in a cooling channel of 8 mm diameter with swirl-tape inserts of 2 ratio and 1.5 mm thickness. The results show excellent agreement for both steady and transient states and prove the successful implementation of both the hydraulic and heated diameters of the swirl-tape channels in the used heat transfer correlations.     

Simulating test for thermal mixing in the hot gas chamber of the HTR-10

The helium coolant at the outlet of the pebble bed core of the 10 MW High Temperature Gas-cooled Reactor-Test Module exhibits a severe radial temperature deviation. In order to avoid damages at the downstream components due to alternating thermal loads such as the steam generator, a hot gas chamber is especially designed to solve the problem. Thermal mixing performance of the coolant in the hot gas chamber is experimentally investigated on a 1:1.5 scale model by air. The experimental result shows that within the Reynolds number range of 1.4×10⁵–5.8×10⁵, the hot gas chamber with a radial mixer reaches excellent thermal mixing of the coolant of about 94%. The flow resistance coefficient for the hot gas chamber is also presented.     

Boron dilution transients during natural circulation flow in PWR—Experiments and CFD simulations

Coolant mixing in the cold leg, downcomer and the lower plenum of pressurized water reactors is an important phenomenon mitigating the reactivity insertion into the core. Therefore, mixing of the de-borated slugs with the ambient coolant in the reactor pressure vessel was investigated at the four loops 1:5 scaled Rossendorf coolant mixing model (ROCOM) mixing test facility. In particular thermal hydraulics analyses have shown, that weakly borated condensate can accumulate in the pump loop seal of those loops, which do not receive a safety injection. After refilling of the primary circuit, natural circulation in the stagnant loops can re-establish simultaneously and the de-borated slugs are shifted towards the reactor pressure vessel (RPV).In the ROCOM experiments, the length of the flow ramp and the initial density difference between the slugs and the ambient coolant was varied. From the test matrix experiments with 0 resp. 2% density difference between the de-borated slugs and the ambient coolant were used to validate the CFD software ANSYS CFX. To model the effects of turbulence on the mean flow a higher order Reynolds stress turbulence model was employed and a mesh consisting of 6.4 million hybrid elements was utilized. Only the experiments and CFD calculations with modeled density differences show stratification in the downcomer. Depending on the degree of density differences the less dense slugs flow around the core barrel at the top of the downcomer. At the opposite side, the lower borated coolant is entrained by the colder safety injection water and transported to the core. The validation proves that ANSYS CFX is able to simulate appropriately the flow field and mixing effects of coolant with different densities.     

Proper orthogonal decomposition of the flow in a tight lattice rod-bundle

Axial coolant flow inside a tightly packed rod bundle presents a complex behavior; experimental analysis had clearly shown that when reducing the pitch-to-diameter ratio (P/D) the turbulence field in rod bundles deviates significantly from that in a circular tube. Moreover for extremely tight configurations the existence of large-scale periodic “flow oscillations” has been shown, which is responsible for the high inter-sub-channel heat and momentum exchange. A complete understanding of these oscillations has still to be achieved; the evidence shown to this point suggests that the oscillations are connected to interactions between coherent structures in adjacent sub-channels. Moreover, the coherent structures show a truly three-dimensional pattern (i.e., structures in different gaps tend to interact) that has not been fully investigated up to this point.A fully transient simulation of turbulence has been performed for an infinite tight triangular lattice. In this case it has been performed with Large Eddy Simulation (LES) at Re = 6400 and P/D = 1.05. A database of snapshots of the flow field has then been collected and proper orthogonal decomposition (POD), a powerful statistical technique, has been applied in order to obtain the most energetic modes of turbulence. The results obtained highlight the presence of several travelling waves propagating in the streamwise direction. The spatial modulation of the travelling waves offers a phenomenological explanation for observed de-coherence effects between the velocity signal in different gaps. It also provides additional insight into the three dimensional structure of the flow oscillations.     

三环路压水堆压力容器上腔室交混矩阵数值研究

使用STAR-CCM+软件对三环路压水堆压力容器上腔室流场进行了大规模、精细化三维数值模拟，并采用组分跟踪方法分别对157个燃料组件出口冷却剂流动进行计算，构造了一个具有3×157个元素的“上腔室交混矩阵”，用该矩阵即可定量、精确地描述冷却剂从堆芯流出后，经上腔室内交混并再分配到各热管道的复杂流动过程。研究发现堆芯流出的冷却剂在压力容器上腔室内的交混是并不充分的，径向上不同位置燃料组件流出的冷却剂会在上腔室同热管道的接口区域存在明显的对应关系，而燃料组件径向功率分布的差异必然导致热管道中冷却剂热分层现象的产生。      

From discovery to recognition of periodic large scale vortices in rod bundles as source of natural mixing between subchannels—A review

The mixing of cooling fluid in rod bundles from one subchannel to another through the gaps between the rods reduces the temperature differences in the coolant as well as along the perimeter of the rods. The phenomenon of natural mixing was first intensively investigated in the 1960s and remains a topic of research up to the present time. The paper describes the main stations on the way to understand the nature of the flow in rod bundles and generally in compound channels with the focus on work performed at Research Center Karlsruhe (FZK).¹Earlier, it was noticed that the mixing rates where higher than could be accounted for by turbulent diffusion alone. For more than 20 years attempts were made to prove experimentally and by code application that secondary flows could account for the measured mixing rates, although the measured secondary flow velocities were much too low. Measurements of the turbulence structure by hot wire anemometry confirmed the existence of cyclic flow pulsations, which had been postulated earlier on the basis of thermocouple measurements. More sophisticated hot wire measurements revealed the nature of these pulsations as periodic, coupled to gap width and Reynolds number. Finally, the extension of the investigation to other compound channel types and flow visualization revealed the true nature of the mixing process as a vortex train moving along the gap between rods or in the narrow part of a compound channel. These findings have been confirmed by LES calculations. Based on these results CFD codes with improved turbulence models calculated successfully the flow in rod bundles including the macroscopic oscillations.     

Safety analysis of a supercritical-pressure water-cooled fast reactor under supercritical pressure

The supercritical-pressure water-cooled fast reactor (SWFR) is a fast spectrum supercritical water-cooled reactor (SCWR) studied by the University of Tokyo. The SWFR is designed as a two-pass core with an outlet temperature 500 °C. The SWFR has fuel channels cooled by downward flow, higher power density, and smaller coolant density reactivity feedback compared with Super LWR. This paper describes the safety analyses of abnormal events for the SWFR. SPRAT-F code is used for the safety analysis at supercritical pressure considering the downward flow cooled seed fuel channel. This code is based on a 1-D node junction model with point kinetics and decay heat calculations. Flow redistribution among parallel paths is calculated by pressure-loss balance and momentum conservation. The initiating events are selected from those of LWRs. For the safety analysis, nine abnormal transients and four accidents are selected with considering types of abnormality. By the numerical analyses, it was found that the loss of flow events can be mitigated by the “water source” effect of the downward flow blanket channels in the abnormal transients and accidents. All the abnormal events satisfy the criteria with margin.     

Thermal and stability considerations for a supercritical water-cooled fast reactor with downward-flow channels during power-raising phase of plant startup

This paper describes study on the procedure of raising the reactor thermal power and the reactor coolant flow rate during the power-raising phase of plant startup for the supercritical water-cooled fast reactor (SWFR), which is selected as one of the Generation IV reactor concepts. Since part of the seed fuel assemblies and all the blanket fuel assemblies of the SWFR are cooled by downward flow, the feedwater from the reactor vessel inlet nozzle to the mixing plenum located below the core is distributed among these fuel assemblies and the downcomer. The flow rate distribution as the function of both the reactor thermal power and the feedwater flow rate, which are the design parameters for the power-raising phase, is obtained by the thermal hydraulic calculations. Based on the flow rate distribution, thermal analyses and thermal-hydraulic stability analyses are carried out in order to obtain the available region of the reactor thermal power and the feedwater flow rate for the power-raising phase. The criteria for the “available” region are the maximum cladding surface temperature (MCST) and the decay ratio of thermal-hydraulic stability in three “hot” channels; two seed assemblies with upward/downward flow and a blanket assembly. The effects of various heat transfer correlations and axial power distributions are also studied.     

Establishment of the design requirements for a flow blockage detection system through a LES analysis of the temperature fluctuation in the upper plenum

When an assembly in a core is partially blocked, the temperature in the upper plenum fluctuates at an early stage. Therefore, the temperature fluctuation in the upper plenum can provide the information about a local blockage of an assembly. For developing the detection algorithm for the partial blockage, we analyzed the temperature fluctuation in the upper plenum due to the partial blockage in an assembly. The LES turbulence model in the CFX code was used for analyzing the temperature fluctuation in the upper plenum because the LES is suitable for analyzing the time dependent turbulence variables. After analyzing the temperature fluctuations in the upper plenum, we established basic design requirements for the flow blockage detection system through a FFT analysis and some statistical analysis. We concluded that response time of a measuring device was less than 13 m s and that it should cover a high temperature range of 1000 K. In addition, the resolution of the thermocouple was less than 2 K and its location should be within 25 cm from the exit of each assembly.     

Undesirable flow behavior in a proposed validation data set

The next generation nuclear plant (NGNP), whose development is supported by the U.S. Department of Energy, will be a very high temperature reactor (VHTR). The VHTR is a single-phase helium-cooled reactor that will provide helium at up to 1000 °C. The prospect of a coolant at these temperatures circulating in the reactor vessel demands that careful analysis be performed to ensure that excessively hot spots are not created and that sufficient mixing of the coolant is obtained. Computational fluid dynamics (CFD) coupled with heat transfer will be used to perform the desired analyses. However, primarily because of the imperfect nature of modeling turbulent flow, any CFD calculations used to perform nuclear reactor safety analysis must be validated against experimental data. Experimental data have been taken in a scaled section of the lower plenum of a prismatic VHTR at the matched index of refraction (MIR) facility at the Idaho National Laboratory. These data were taken with the intent that they be examined for use as validation data. A series of investigations have been conducted to assess the MIR data. Issues that have already been examined include the extent of the required computational domain, the outlet boundary condition, the inlet data and the effect of the turbulence model. One of the jets that flow into the model impacts on a wedge, which represents a portion of a hexagonal graphite block that lines the inner wall of the lower plenum. The nature of the flow below this particular jet is such that a randomly varying recirculation zone is created. This recirculation zone is seen to change in size, causing a relatively long-time scale of motion or disturbance of the flow in the model. It is concluded that such a feature is undesirable in a validation data set, firstly because of its apparent random nature and, secondly, because to obtain an appropriate long-time average would be impractical because of the compute time required. It is predicted computationally that by eliminating the first of the four inlet jets into the scaled model, the resulting recirculation zone is rendered stable.     

Design and analysis of a thermal core for a high performance light water reactor

The High Performance Light Water Reactor is a Generation IV light water reactor concept, operated at a supercritical pressure of 25 MPa with a core outlet temperature of 500 °C. A thermal core design for this reactor has been worked out by a consortium of Euratom member states within the 6th European Framework Program. Aiming at peak cladding temperatures of less than 630 °C, including uncertainties and allowances for operation, the coolant is heated up in three steps with intermediate coolant mixing to eliminate hot streaks. Different from conventional reactors, the radial power profile is intended to be non-uniform, with the highest power in the first heat-up step in the core center and the lowest power in the second superheater step to result in the same peak cladding temperatures in each region. The concept has been studied with neutronic, thermal-hydraulic and structural analyses to assess its feasibility. Coupled neutronic/thermal-hydraulic analyses are defining the initial distribution of enrichment, control rod positions and the use of burnable poisons. Sub-channel analyses predict the coolant mixing inside assemblies, and a porous media approach simulates the flow of moderator water between assembly boxes. Finally, structural analyses of the assembly boxes are needed to minimize deformations during operation. Even though the core design cannot yet considered to be final, this state of the art review shall summarize the progress achieved so far and outline the remaining challenges.     

Numerical simulation of three-dimensional thermal-hydraulic behavior for HTTR (High Temperature Engineering Test Reactor)

Safety demonstration tests using the HTTR (High Temperature Engineering Test Reactor) are now in progress in order to verify the inherent safety features and to improve safety designing and analysis technologies for future HTGR (high temperature gas-cooled reactor). Coolant flow reduction test is one of the safety demonstration tests for the purpose of demonstration of inherent HTGR safety features in the case that coolant flow is reduced by tripping of helium gas circulators. If reactor core element temperature and core internal structure temperature during abnormal events are estimated by numerical simulation with high-accuracy, developed numerical simulation method can be applied to future HTGR designing efficiently. In the present research, three-dimensional in- and ex-vessel thermal-hydraulic calculations for the HTTR are performed with a commercially available thermal-hydraulic analysis code “STAR-CD^®” with finite volume method. The calculations are performed for normal operation and coolant flow reduction tests of the HTTR. Then calculated temperatures are compared with measured ones obtained in normal operation and coolant flow reduction test. As the result, calculated temperatures are good agreement with measured ones in normal operation and coolant flow reduction test.     

Influence of injection temperature and flow rate on mixing and stratification in small passively cooling enclosures

Thermal mixing and stratification phenomena may occur during the loss of a coolant accident or main steam line break accident in the containment of a Passive Containment Cooling System, or in the suppression pools in BWR. However, the present study pays insufficient attention to the thermal stratification phenomena in the containment of small modular reactors(SMR). In this paper, an investigation on the mixing and thermal stratification phenomena caused by the plumes or buoyant jets in SMR containments was carried out. The experiments were both conducted under non-adiabatic and adiabatic conditions for a steel containment. In each condition, two key parameters, inlet temperature, and flow rate were tested by controlling variables to identify their influence on the thermal stratification phenomenon. The visualization experiments illustrated the jet mixing and stratification development. The experiment results were compared with the numerical computation and they reached a good agreement.     

In situ determination of radon concentration and total gamma radiation in Kastel Gomilica, Croatia

To determine current radiation background of the environment at the “Giricic” location in Kastel Gomilica, Croatia, in situ measurement of radon concentration (²²²Rn and ²²⁰Rn) in an open atmosphere on a ground level and at the height of 1.5 m has been made as well as total gamma radiation at the height of 1 m in an energy range of 15 keV to 2 MeV. The researched location was divided in three specific parts: (i) regulated area with the bottom ash and flying ash in the basis (“old” depot), (ii) unregulated area with waste materials, including bottom ash and flying ash, in the basis (“new” depot), (iii) uncontaminated area with no waste materials deposited on. Average radon concentration on a ground level was 213 Bq/m³ for the “old” depot, 214 Bq/m³ for the “new” depot and 59 Bq/m³ for the uncontaminated area and at the height of 1.5 m 20 Bq/m³ for the “old” depot, 34 Bq/m³ for the “new” depot and 26 Bq/m³ for the uncontaminated area. Average total gamma radiation values in selected energy range were 109.92 cps (counts per second) for the “old” depot, 357.76 cps for the “new” depot and 65.97 cps for the uncontaminated area. For selected radionuclides (²¹⁴Pb, ¹³⁷Cs, ²²⁸Ac, ^234mPa, ⁴⁰K and ²¹⁴Bi) average gamma radiation values at characteristic energies have been determined as well.