OECD/NEA ROSA/LSTF project tests are performed on the Large Scale Test Facility (LSTF). LSTF is a full-height, full-pressure and 1/48 volumetrically-scaled two-loop system which aims to simulate Japanese Tsuruga-2 Westinghouse-type 4-loop PWR. ROSA-V Test 6-1 simulates a pressure vessel (PV) upper-head small break loss-of-coolant accident (SBLOCA) with a break size equivalent to 1.9% of the volumetrically scaled cross-sectional area of the reference PWR cold leg.The main objective of present thesis is to build a TRACE calculation model for simulating thermal hydraulic behaviors in LSTF and PV upper-head SBLOCA, so as to assess different modeling options and parameters of TRACE code. The results show that TRACE code well reproduce the complex physical phenomena involved in this type of SBLOCA scenarios. Almost all the events in the experiment are well predicted by the model based on TRACE code. In addition, the sensitivity of different models and parameters are investigated. For example, the code slightly overestimates the break mass flow from upper head which affects the accuracy of the results significantly. The rising of core exit temperature (CET) is significantly influenced by the flow area of leakage between downcomer and hot leg. Besides, the effect of the break location, low pressure injection system (LPIS) and accumulator setup are also studied.

As a latest developed computational code, TRACE is expected to be useful and effective for analyzing the thermal-hydraulic behaviors in design, licensing and safety analysis of nuclear power plant. However, its validity and correctness have to be verified and qualified before its application into industry. Loss-of-coolant accident (LOCA) is a kind of transient thermal hydraulic event which has been emphasized a lot as a most important threat to the safety of the nuclear power plant. In the present study, based on FIX- II LOCA tests, simulation models for the tests of No. 3025, No. 3061 and No. 5052 were developed to validate the TRACE code (version 5.0 patch 2). The simulated transient thermal-hydraulic behaviors during the LOCA tests including the pressure in the primary system, the mass flow rate in certain key parts, and the temperature in the core are compared with experimental data. The simulation results show that TRACE model can well reproduce the transient thermal-hydraulic behaviors under different LOCA situations. In addition, sensitivity analysis are also performed to investigate the influence of particular models and parameters, including counter current flow limitation (CCFL) model, choked flow model , insulator in the steam dome, K-factor in the test section, and pump trip, on the results. The sensitivity analyses show that both the models and parameters have significant influence on the outcome of the model.

With the coming of Generation IV systems, there is a need for thermal-hydraulic codes to model such advanced reactors. Codes for High Temperature Gas-cooled Reactors (HTGRs) already exist, but often suffer from insufficient validation and little user experience. Therefore, some existing codes created for Light Water Reactors are updated with HTGR related features. In this study, the feasibility of providing the TRACE thermal-hydraulics code with those features is analysed. Two models were used, one of a pebble bed core, one of a prismatic block reactor. For this purpose the effective conductivity test of the High Temperature Test Unit was taken as a benchmark for the pebble bed core. For the prismatic block reactor a model of the GEMINI+ reactor was created. This would allow to simulate not only steady state, but also Depressurised Loss of Forced Cooling scenarios. For both models the effective conductivity is known to play a major role and, consequently, a model to incorporate such feature was developed and implemented within TRACE's control system module. Results show that TRACE has a good potential for HTGR simulation, but currently available models still provide unstable solutions. It is concluded that TRACE needs additional adjustments in order to be employed for HTGR safety analyses in the future. / Med fjärde generationens system på ingång, finns det ett behov av termisk-hydrauliska koder för att modellera sådana avancerade reaktorer. Koder för gaskylda högtemperaturreaktorer (HTGR) finns redan, men lider ofta av otillräcklig validering och liten användarupplevelse. Därför uppdateras vissa befintliga koder som skapats för lättvattenreaktorer med HTGR-relaterade funktioner. I denna studie analyseras möjligheten att tillhandahålla TRACE termisk-hydraulisk kod med dessa funktioner. Två modeller användes, den ena av en pebble-bed reaktor, den andra av en prismatisk blockreaktor. För detta ändamål togs det effektiva konduktivitetstestet för högtemperaturtestenheten som ett riktmärke för pebble-bedens härd. För den prismatiska blockreaktorn skapades en modell av GEMINI+-reaktorn. Detta skulle göra det möjligt att simulera inte bara steady state, utan även scenarier med trycklös förlust av forcerad kylning. För båda modellerna är den effektiva konduktiviteten känd för att spela en stor roll och följaktligen utvecklades och implementerades en modell för att införliva en sådan funktion inom TRACEs kontrollsystemmodul.. Resultaten visar att TRACE har en god potential för HTGR-simulering, men för närvarande tillgängliga modeller ger fortfarande instabila lösningar. Slutsatsen är att TRACE behöver ytterligare justeringar för att kunna användas för HTGR-säkerhetsanalyser i framtiden.

No description available.

The effects of appendages on flow characteristics and scalar mixing in gap-connected twin-subchannel geometries has been assessed. The assessment considers a symmetric, rectangular compound channel geometry connected by a single rectangular gap using computational fluid dynamics (CFD). Detailed numerical models (geometry and turbulence), characterizing the full test section from a reference experimental study, are generated and validated against measurements. Time varying details of the gap induced periodic structures and appendage induced vortices are captured through calculations in an unsteady Reynolds averaged Navier-Stokes (RANS) framework coupled with the Spalart-Allmaras (SA) turbulence model closing the RANS equations. Companion simulations are performed at each of two Reynolds numbers (2690 and 7500), one with and one without a gap-centered appendage. The appendage size modelled is representative of CANDU endplates. The appendage effects on flow characteristics and mixing are isolated through comparison of the associated simulations. In the absence of appendages, fluid exchange between subchannels is dominated by quasi-periodic flow pulsations through the gap formed due to flow instability in the near gap region. Without a gap-centered appendage, the magnitude, frequency, and structure length of the gap flow pulsations are well predicted by the model at both Reynolds numbers. The total tracer transfer between subchannels is reasonably well predicted for Re = 2690 (within approximately 17% of the experimental value). The model fails to capture the measured increase in scalar transfer through the gap with increased Reynolds number, underpredicting scalar mixing by 55% at Re = 7500. An argument is presented that the use of an isotropic turbulence model in the channel (SA), which precludes the development of channel secondary flows, is the source of the discrepancy between modelled and measured mixing at Re = 7500. Appendages, such as those introduced by end plates or bearing pads in CANDU fuel bundles, augment the exchange process between subchannels. With an appendage representative of a CANDU fuel bundle endplate introduced into the gap region, crossflow velocity and frequency are predicted to increase immediately downstream of the appendage due to flow diversion and vortex shedding. The higher local frequency is shown to be consistent with the vortex shedding frequency calculated for a stationary rectangular cylinder at the gap conditions. Further downstream, gap induced instabilities begin to re-establish as the dominant contributor to crossflow pulsations although they are not fully recovered by the test section exit. Mixing is augmented more by the appendage with increasing Reynolds number for the range examined. / Thesis / Master of Applied Science (MASc) / The fuel bundle and pressure tube assembly in the core of a CANDU reactor forms an intricate web of subchannels of varying geometries with interconnecting gaps. Heat generated within the fuel bundles is removed by coolant flowing through the pressure tube and within the bundle subchannels. Although flow is nominally axial along the length of the rod bundles, coolant is free to move between subchannels through the gaps by a variety of mechanisms. Detailed fluid flow in these rod bundle geometries is a complex 3D phenomenon, strongly affected by fluid turbulence and flow instabilities associated with the subchannel geometry. This flow is investigated in the current study and extended to include the effect of appendages, which hold the fuel rods in place, to determine their impact on mixing along the length of the bundle. Particular applications of the results of this study are in the areas of nuclear reactor performance and safety. The extent of coolant exchange between subchannels affects the local subchannel flow and temperature and, as a result, local cooling at the fuel element surfaces. Fuel element cooling is a principal component of reactor analysis under both normal operating conditions and postulated accident scenarios.

Les sels fondus ont récemment été proposés comme milieux caloporteurs à haute température. Dans l'industrie nucléaire, le concept de réacteur à sels fondus (MSR en anglais) est le seul concept de quatrième génération qui propose l'utilisation d'un sel fondu liquide comme combustible nucléaire. Cette innovation présente des aspects positifs pour la conception et la sûreté nucléaire, mais impose de nouveaux défis. Le réacteur rapide à sels fondus (MSFR en anglais) est un concept qui est actuellement étudié dans le projet européen H2020 SAMOFAR, incluant le développement et la validation expérimentale (dans la plateforme expérimentale SWATH) de modèles plus performants pour les sels fondus : tel est l'objectif de ce travail de thèse. En outre les modèles développés peuvent s'appliquer à d'autres MSRs et à d'autres applications énergétiques utilisant des sels fondus comme milieux caloporteurs.La thèse suivante est divisée en trois parties :Premièrement, le développement de modèles pour décrire de façon réaliste certains des phénomènes thermiques microscopiques et macroscopiques associés à l’utilisation de sels liquides fondus comme milieux caloporteurs. Cette partie comprend l’utilisation et le développement de nouveaux modèles neutroniques pour étudier la production d'énergie nucléaire, ainsi que la modélisation des phénomènes turbulents dans les sels fondus, l’étude de l’interaction du rayonnement thermique et la turbulence dans les sels fondus. Enfin, cette partie traite également du développement d’une approche multi-échelle pour l'étude précise de la solidification/fusion dans les sels.Deuxièmement, la conception et la mise en œuvre d’expériences dédiées à la validation de ces modèles. Deux expériences clés ont été conçues au cours de cette thèse et ont été implémentées dans la plate-forme SWATH. L'objectif de ces expériences est d'étudier le comportement de différents modèles de turbulence et de tester les modèles de solidification développés dans les sels fondus.Troisièmement, les modèles développés ont été couplés dans une plateforme multi-physique pour l'étude précise du transitoire drainant du MSFR. / Molten salts have been recently proposed as high-temperature heat carrier media for energy applications. In the nuclear industry, the Molten Salt Reactors (MSRs) are the only fourth generation concept proposing the usage of a liquid nuclear fuel. This innovative aspect allows proposing improved safety and design features, but it leads to novel challenges. In particular, the Molten Salt Fast Reactor (MSFR) is a MSR concept that is currently being studied in the H2020 European project SAMOFAR. Among the project activities, there are the development of more performant molten salts models and their experimental validation through the SWATH platform. This is the objective of the present thesis. However, the models developed are appropriate for other MSRs and other energy applications using molten salts as heat carrier media.The following thesis is divided into three parts.The first part is dedicated to the development of models for describing realistically some of the microscopic and macroscopic thermal phenomena associated with the usage of liquid molten salts as heat carrier media. This part includes the development and implementation of neutronic models to study nuclear power production in the MSFR, the study of turbulence and turbulence-radiation interaction in molten salt flows and the development of a multiscale approach to model the solidification/melting phenomena in salts.The second part is devoted to the design and implementation of dedicated experiments for validating these models. Two key experiments are addressed: an experiment to study the behavior of different turbulence models after a boundary layer detachment and one to test the multiscale solidification models developed for molten salts.The third part is committed to the coupling of the models developed into a multiphysics platform for the precise study of the draining transient of the MSFR.

The system code TRACE is the latest development of the U.S. Nuclear Regulatory Commission (US NRC). TRACE, developed for the analysis of operational conditions, transients and accidents of light water reactors (LWR), is a best-estimate code with two fluid, six equation models for mass, energy, and momentum conservation, and related closure models. Since TRACE is mainly applied to LWR specific issues, the validation process related to innovative nuclear systems (liquid metal cooled systems, systems operated with supercritical water, etc.) is very limited, almost not existing. In this work, essential contribution to the validation of TRACE related to lead and lead alloy cooled systems as well as systems operated with supercritical water is provided in a consistent and corporate way. In a first step, model discrepancies of the TRACE source code were removed. This inconsistencies caused the wrong prediction of the thermo physical properties of supercritical water and lead bismuth eutectic, and hence the incorrect prediction of heat transfer relevant characteristic numbers like Reynolds or Prandtl number. In addition to the correction of the models to predict these quantities, models describing the thermo physical properties of lead and Diphyl THT (synthetic heat transfer medium) were implemented. Several experiments and numerical benchmarks were used to validate the modified TRACE version. These experiments, mainly focused on wall-to-fluid heat transfer, revealed that not only the thermo physical properties are afflicted with inconsistencies but also the heat transfer models. The models for the heat transfer to liquid metals were enhanced in a way that the code can now distinguish between pipe and bundle flow by using the right correlation. The heat transfer to supercritical water was not existing in TRACE up to now. Completely new routines were implemented to overcome that issue. The comparison of the calculations to the experiments showed, on one hand, the necessity of these changes and, on the other hand, the success of the new implemented routines and functions. The predictions using the modified TRACE version were close to the experimental data. After validating the modified TRACE version, two design studies related to the Generation IV International Forum (GIF) were investigated. In the first one, a core of a lead-cooled fast reactor (LFR) was analyzed. To include the interaction between the thermal hydraulic and the neutron kinetic due to temperature and density changes, the TRACE code was coupled to the program system ERANOS2.1. The results gained with that coupled system are in accordance with theory and helped to identify sub-assemblies with the highest loads concerning fuel and cladding temperature. The second design which was investigated was the High Performance Light Water Reactor (HPLWR). Since the design of the HPLWR is not finalized, optimization of vital parameters (power, mass flow rate, etc.) are still ongoing. Since most of the parameters are affecting each other, an uncertainty and sensitivity analysis was performed. The uncertainty analysis showed the upper and lower boundaries of selected parameters, which are of importance from the safety point of view (e.g., fuel and cladding temperature, moderator temperature). The sensitivity study identified the most relevant parameters and their influence on the whole system.

Das Vorhaben bestand in der Validierung des gekoppelten neutronenkinetisch-thermohydraulischen Programmkomplexes ATHLET/DYN3D für Siedewasserreaktoren durch Teilnahme an dem OECD/NRC Benchmark zum Turbinenschnellschluss. Das von der OECD und der amerikanischen NRC definierte Benchmark basiert auf einem Experiment mit Schließens des Turbinenschnellschlussventils, das 1977 im Rahmen einer Serie von 3 Experimenten im Kernkraftwerk Peach Bottom 2 durchgeführt wurde. Im Experiment erzeugte das Schließen des Ventils eine Druckwelle, die sich unter Abschwächung bis in den Reaktorkern ausbreitete. Die durch den Druckanstieg bewirkte Kondensation von Dampf im Reaktorkern führte zu einem positiven Reaktivitätseintrag. Der folgende Anstieg der Reaktorleistung wurde durch die Rückkopplung und das Einfahren der Regelstäbe begrenzt. Im Rahmen des Benchmarks konnten die Rechenprogramme durch Vergleiche mit den Messergebnissen und den Ergebnissen der anderen Teilnehmer an dem Benchmark validiert werden. Das Benchmark wurde in 3 Phasen oder Exercises eingeteilt. Die Phase I diente der Überprüfung des thermohydraulischen Modells für das System bei vorgegebener Leistungsfreisetzung im Kern. In der Phase II wurden 3-dimensionale Berechnungen des Reaktorkerns für vorgegebene thermohydraulische Randbedingungen durchgeführt. Die gekoppelten Rechnungen für das ausgewählte Experiment und für 4 extreme Szenarien erfolgten in der Phase III. Im Rahmen des Projekts nahm FZR an Phase II und Phase III des Benchmarks teil. Die Rechnungen für Phase II erfolgten mit dem Kernmodell DYN3D unter Berücksichtigung der Heterogenitätsfaktoren und mit 764 thermohydraulischen Kanälen (1 Kanal/Brennelement). Der ATHLET-Eingabedatensatz für die Reaktoranlage wurde von der Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) übernommen und für die Rechnungen zu Phase III, die mit der parallelen Kopplung von ATHLET mit DYN3D erfolgten, geringfügig modifiziert. Für räumlich gemittelte Parameter wurde eine gute Übereinstimmung mit den Messergebnissen und den Resultaten anderer Codes erzielt. Der Einfluss der Modellunterschiede wurde mit Hilfe von Variantenrechnungen zu Phase II untersucht. So können Unterschiede in der Leistungs- und Voidverteilung in einzelnen Brennelementen auf die unterschiedliche neutronenkinetische und thermohydraulische Modellierung des Reaktorkerns zurückgeführt werden. Vergleiche zwischen ATHLET/DYN3D (parallele Kopplung) und ATHLET/QUABOX-CUBBOX (interne Kopplung) zeigen für räumlich gemittelte Parameter nur geringe Unterschiede. Abweichungen in den lokalen Parametern können im wesentlichen mit der unterschiedlichen Modellierung des Reaktorkerns erklärt werden (geringere Anzahl von modellierten Kühlkanälen, keine Berücksichtigung der Heterogenitätsfaktoren und ein anderes Siedemodell in der Rechnung mit ATHLET/QUABOX-CUBBOX). Die Rechnungen für die extremen Szenarien von Phase III zeigen die Anwendbarkeit des gekoppelten Programms ATHLET/DYN3D für die Bedingungen bei Störfällen, die weit über das Experiment hinausgehen.

Benchmark calculations for the validation of the coupled neutron kinetics/thermohydraulic code complex DYN3D-ATHLET are described. Two benchmark problems concerning hypothetical accident scenarios with leaks in the steam system for a VVER-440 type reactor and the TMI-1 PWR have been solved. The first benchmark task has been defined by FZR in the frame of the international association "Atomic Energy Research" (AER), the second exercise has been organised under the auspices of the OECD. While in the first benchmark the break of the main steam collector in the sub-critical hot zero power state of the reactor was considered, the break of one of the two main steam lines at full reactor power was assumed in the OECD benchmark. Therefore, in this exercise the mixing of the coolant from the intact and the defect loops had to be considered, while in the AER benchmark the steam collector break causes a homogeneous overcooling of the primary circuit. In the AER benchmark, each participant had to use its own macroscopic cross section libraries. In the OECD benchmark, the cross sections were given in the benchmark definition. The main task of both benchmark problems was to analyse the re-criticality of the scrammed reactor due to the overcooling. For both benchmark problems, a good agreement of the DYN3D-ATHLET solution with the results of other codes was achieved. Differences in the time of re-criticality and the height of the power peak between various solutions of the AER benchmark can be explained by the use of different cross section data. Significant differences in the thermohydraulic parameters (coolant temperature, pressure) occurred only at the late stage of the transient during the emergency injection of highly borated water. In the OECD benchmark, a broader scattering of the thermohydraulic results can be observed, while a good agreement between the various 3D reactor core calculations with given thermohydraulic boundary conditions was achieved. Reasons for the differences in the thermohydraulics were assumed in the difficult modelling of the vertical once-through steam generator with steam superheating. Sensitivity analyses which considered the influence of the nodalisation and the impact of the coolant mixing model were performed for the DYN3D-ATHLET solution of the OECD benchmark. The solution of the benchmarks essentially contributed to the qualification of the code complex DYN3D-ATHLET as an advanced tool for the accident analysis for both VVER type reactors and Western PWRs.

Der Reaktordynamikcodes DYN3D wurde in der neu entwickelten Mehrgruppen-Version DYN3D-MG für die Anwendung auf wassergekühlte Reaktoren alternativ zu industriellen DWR und SWR ertüch-tigt. Es wurde die Anwendbarkeit für den graphitmoderierten Druckröhrenreaktor EGP-6 (KKW Bilibi-no), eine Konzeptstudie eines fortgeschrittenen Siedewasserreaktors mit schnellem Neutronenspekt-rum (RMWR) und das Reaktorkonzept RUTA-70 zur Wärmeversorgung nachgewiesen. Beim RUTA-Reaktor geht es vor allem um die Modellierung des Naturumlaufs des Kühlmittels bei niedrigen Sys-temdrücken. Zur Validierung wurden Experimente zu flashing-induzierten Naturumlaufinstabilitäten an der Versuchsanlage CIRCUS der TU Delft mit RELAP5 nachgerechnet. Für die Anwendung von DYN3D auf die alternativen Reaktorkonzepte wurden Modellerweiterungen und Anpassungen vorgenommen, u.a. Modifikationen in den Wärmeleitungs- und -übergangsmodellen. Vergleichsrechnungen mit dem stationären russischen Feingitter-Diffusionscode ACADEM ergänzen die Verifikationsdatenbasis von DYN3D-MG. Zur Validierung wurden zwei reak-tordynamische Experimente am Reaktor EGP-6 nachgerechnet. Für Reaktoren EGP-6, RMWR und RUTA wurden verschiedene Transienten mit Ausfahren von Re-gelstäben mit und ohne Reaktorschnellabschaltung gerechnet. Weiterhin wurden Analysen für den ATWS-Störfall \"Abschalten aller Hauptkühlmittelpumpen bei Vollleistung\" für den RUTA-Reaktor mit den gekoppelten Programmkomplexen DYN3D/ATHLET und DYN3D/RELAP5 durchgeführt. Der Reaktor geht in einen sicheren Zustand mit reduzierter Leistung bei Naturumlauf des Kühlmittels über. Die Ergebnisse von Analysen zum unkontrollierten Ausfahren einer Regelgruppe für den RMWR lassen dagegen eine belastbare Schlussfolgerung bezüglich der Beherrschbarkeit des Aus-fahrens einer Regelgruppe nicht zu. Abschließend wurde der Nutzen der Programmertüchtigung von DYN3D für die Anwendung auf GenIV -Konzepte und LWR mit hohem Konversionsfaktor bewertet.