<p>Severe accidents in a Light Water Reactor (LWR) have been a subject of the research for the last three decades. The research in this area aims to further understanding of the inherent physical phenomena and reduce the uncertainties surrounding their quantification, with the ultimate goal of developing models that can be applied to safety analysis of nuclear reactors. The research is also focusing on evaluation of the proposed accident management schemes for mitigating the consequences of such accidents.</p><p>During a hypothetical severe accident, whatever the scenario, there is likelihood that the core material will be relocated and accumulated in the lower plenum in the form of a debris bed or a melt pool. Physical phenomena involved in a severe accident progression are complex. The interactions of core debris or melt with the reactor structures depend very much on the debris bed or melt pool thermal hydraulics. That is why predictions of heat transfer during melt pool formation in the reactor lower head are important for the safety assessment.</p><p>The main purpose of the present study is to advance a method for describing turbulent natural convection heat transfer of a melt pool, and to develop a computational platform for cost-effective, sufficiently-accurate numerical simulations and analyses of Core Melt-Structure-Water Interactions in the LWR lower head during a postulated severe core-melting accident.</p><p>Given the insights gained from Computational Fluid Dynamics (CFD) simulations, a physics-based model and computationally-efficient tools are developed for multi-dimensional simulations of transient thermal-hydraulic phenomena in the lower plenum of a Boiling Water Reactor (BWR) during the late phase of an in-vessel core melt progression. A model is developed for the core debris bed heat up and formation of a melt pool in the lower head of the reactor vessel, and implemented in a commercial CFD code. To describe the natural convection heat transfer inside the volumetrically decay-heated melt pool, we advanced the Effective Convectivity Conductivity Model (ECCM), which was previously developed and implemented in the MVITA code. In the present study, natural convection heat transfer is accounted for by only the Effective Convectivity Model (ECM). The heat transport and interactions are represented through an energy-conservation formulation. The ECM then enables simulations of heat transfer of a high Rayleigh melt pool in 3D large dimension geometry.</p><p>In order to describe the phase-change heat transfer associated with core debris, a temperature-based enthalpy formulation is employed in the ECM (the phase-change ECM or so called the PECM). The PECM is capable to represent possible convection heat transfer in a mushy zone. The simple approach of the PECM method allows implementing different models of the fluid velocity in a mushy zone for a non-eutectic mixture. The developed models are validated by a dual approach, i.e., against the existing experimental data and the CFD simulation results.</p><p>The ECM and PECM methods are applied to predict thermal loads to the vessel wall and Control Rod Guide Tubes (CRGTs) during core debris heat up and melting in the BWR lower plenum. Applying the ECM and PECM to simulations of reactor-scale melt pool heat transfer, the results of the ECM and PECM calculations show an apparent effectiveness of the developed methods that enables simulations of long term accident transients. It is also found that during severe accident progression, the cooling by water flowing inside the CRGTs plays a very important role in reducing the thermal load on the reactor vessel wall. The results of the CFD, ECM and PECM simulations suggest a potential of the CRGT cooling as an effective mitigative measure during a severe accident progression.</p>

Severe accidents in a Light Water Reactor (LWR) have been a subject of intense research for the last three decades. The research in this area aims to reach understanding of the inherent physical phenomena and reduce the uncertainties in their quantification, with the ultimate goal of developing models that can be applied to safety analysis of nuclear reactors, and to evaluation of the proposed accident management schemes for mitigating the consequences of severe accidents.  In a hypothetical severe accident there is likelihood that the core materials will be relocated to the lower plenum and form a decay-heated debris bed (debris cake) or a melt pool. Interactions of core debris or melt with the reactor structures depend to a large extent on the debris bed or melt pool thermal hydraulics. In case of inadequate cooling, the excessive heat would drive the structures' overheating and ablation, and hence govern the vessel failure mode and timing. In turn, threats to containment integrity associated with potential ex-vessel steam explosions and ex-vessel debris uncoolability depend on the composition, superheat, and amount of molten corium available for discharge upon the vessel failure. That is why predictions of transient melt pool heat transfer in the reactor lower head, subsequent vessel failure modes and melt characteristics upon the discharge are of paramount importance for plant safety assessment.  The main purpose of the present study is to develop a method for reliable prediction of melt pool thermal hydraulics, namely to establish a computational platform for cost-effective, sufficiently-accurate numerical simulations and analyses of core Melt-Structure-Water Interactions in the LWR lower head during a postulated severe core-melting accident. To achieve the goal, an approach to efficient use of Computational Fluid Dynamics (CFD) has been proposed to guide and support the development of models suitable for accident analysis.   The CFD method, on the one hand, is indispensable for scrutinizing flow physics, on the other hand, the validated CFD method can be used to generate necessary data for validation of the accident analysis models. Given the insights gained from the CFD study, physics-based models and computationally-efficient tools are developed for multi-dimensional simulations of transient thermal-hydraulic phenomena in the lower plenum of a LWR during the late phase of an in-vessel core melt progression. To describe natural convection heat transfer in an internally heated volume, and molten metal layer heated from below and cooled from the top (and side) walls, the Effective Convectivity Models (ECM) are developed and implemented in a commercial CFD code. The ECM uses directional heat transfer characteristic velocities to transport the heat to cooled boundaries. The heat transport and interactions are represented through an energy-conservation formulation. The ECM then enables 3D heat transfer simulations of a homogeneous (and stratified) melt pool formed in the LWR lower head. In order to describe phase-change heat transfer associated with core debris or binary mixture (e.g. in a molten metal layer), a temperature-based enthalpy formulation is employed in the Phase-change ECM (so called the PECM). The PECM is capable to represent natural convection heat transfer in a mushy zone. Simple formulation of the PECM method allows implementing different models of mushy zone heat transfer for non-eutectic mixtures. For a non-eutectic binary mixture, compositional convection associated with concentration gradients can be taken into account. The developed models are validated against both existing experimental data and the CFD-generated data. ECM and PECM simulations show a superior computational efficiency compared to the CFD simulation method. The ECM and PECM methods are applied to predict thermal loads imposed on the vessel wall and Control Rod Guide Tubes (CRGTs) during core debris heatup and melting in a Boiling Water Reactor (BWR) lower plenum. It is found that during the accident progression, the CRGT cooling plays a very important role in reducing the thermal loads on the reactor vessel wall. Results of the ECM and PECM simulations suggest a high potential of the CRGT cooling to be an effective measure for severe accident management in BWRs. / <p>QC 20100812</p>

Fundamental understanding of the behaviour of nuclear fuel has been of great importance. Enhancing this knowledge not only by means of experimental observations, but also via multi-scale modelling is of current interest. The overall goal of this thesis is to understand the impact of atomic interactions on the nuclear fuel material properties. Two major topics are tackled in this thesis. The first topic deals with non-stoichiometry in uranium dioxide (UO2) to be addressed by empirical potential (EP) studies. The second fundamental question to be answered is the effect of the atomic fraction of americium (Am), neptunium (Np) containing uranium (U) and plutonium (Pu) mixed oxide (MOX) on the material properties.UO2 has been the reference fuel for the current fleet of nuclear reactors (Gen-II and Gen-III); it is also considered today by the Gen-IV International Forum for the first cores of the future generation of nuclear reactors on the roadmap towards minor actinide (MA) based fuel technology. The physical properties of UO2 highly depend on material stoichiometry. In particular, oxidation towards hyper stoichiometric UO2 – UO2+x – might be encountered at various stages of the nuclear fuel cycle if oxidative conditions are met; the impact of physical property changes upon stoichiometry should therefore be properly assessed to ensure safe and reliable operations. These physical properties are intimately linked to the arrangement of atomic defects in the crystalline structure. The first paper evaluates the evolution of defect concentration with environment parameters – oxygen partial pressure and temperature by means of a point defect model, with reaction energies being derived from EP based atomic scale simulations. Ultimately, results from the point defect model are discussed, and compared to experimental measurements of stoichiometry dependence on oxygen partial pressure and temperature. Such investigations will allow for future discussions about the solubility of different fission products and dopants in the UO2 matrix at EP level.While the first paper answers the central question regarding the dominating defects in non-stoichiometry in UO2, the focus of the second paper was on the EP prediction of the material properties, notably the lattice parameter of Am, Np containing U and Pu MOX as a function of atomic fractions.The configurational space of a complex U1-y-y’-y’’PuyAmy’Npy’’O2 system, was assessed via Metropolis-Monte Carlo techniques. From the predicted configuration, the relaxed lattice parameter of Am, Np bearing MOX fuel was investigated and compared with available literature data. As a result, a linear behaviour of the lattice parameter as a function of Am, Np content was observed, as expected for an ideal solid solution. These results will allow to support and increase current knowledge on Gen-IV fuel properties, such as melting temperature, for which preliminary results are presented in this thesis, and possibly thermal conductivity in the future. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

The search for Accident tolerant fuels (ATF) which is the first part of this thesis is currently one of the most actual topics in the field of nuclear fuels. These fuels must be first successfully tested in operational and also accident conditions for their possible inclusion in commercial use. Following part of the thesis specifically focuses on the boiling crisis in nuclear reactors which can damage the nuclear fuel cladding. Therefore, it is necessary to know the critical heat flux value and the departure from nuclear boiling ratio. Calculations which determine critical heal flux value are placed in the practical part of the thesis. Calculations are compared with the data obtained during experiments. The ALTHAMC12 and the other correlations which are based on the previous measurements are used for the computational analysis.

Severe accidents in a Light Water Reactor (LWR) have been a subject of the research for the last three decades. The research in this area aims to further understanding of the inherent physical phenomena and reduce the uncertainties surrounding their quantification, with the ultimate goal of developing models that can be applied to safety analysis of nuclear reactors. The research is also focusing on evaluation of the proposed accident management schemes for mitigating the consequences of such accidents. During a hypothetical severe accident, whatever the scenario, there is likelihood that the core material will be relocated and accumulated in the lower plenum in the form of a debris bed or a melt pool. Physical phenomena involved in a severe accident progression are complex. The interactions of core debris or melt with the reactor structures depend very much on the debris bed or melt pool thermal hydraulics. That is why predictions of heat transfer during melt pool formation in the reactor lower head are important for the safety assessment. The main purpose of the present study is to advance a method for describing turbulent natural convection heat transfer of a melt pool, and to develop a computational platform for cost-effective, sufficiently-accurate numerical simulations and analyses of Core Melt-Structure-Water Interactions in the LWR lower head during a postulated severe core-melting accident. Given the insights gained from Computational Fluid Dynamics (CFD) simulations, a physics-based model and computationally-efficient tools are developed for multi-dimensional simulations of transient thermal-hydraulic phenomena in the lower plenum of a Boiling Water Reactor (BWR) during the late phase of an in-vessel core melt progression. A model is developed for the core debris bed heat up and formation of a melt pool in the lower head of the reactor vessel, and implemented in a commercial CFD code. To describe the natural convection heat transfer inside the volumetrically decay-heated melt pool, we advanced the Effective Convectivity Conductivity Model (ECCM), which was previously developed and implemented in the MVITA code. In the present study, natural convection heat transfer is accounted for by only the Effective Convectivity Model (ECM). The heat transport and interactions are represented through an energy-conservation formulation. The ECM then enables simulations of heat transfer of a high Rayleigh melt pool in 3D large dimension geometry. In order to describe the phase-change heat transfer associated with core debris, a temperature-based enthalpy formulation is employed in the ECM (the phase-change ECM or so called the PECM). The PECM is capable to represent possible convection heat transfer in a mushy zone. The simple approach of the PECM method allows implementing different models of the fluid velocity in a mushy zone for a non-eutectic mixture. The developed models are validated by a dual approach, i.e., against the existing experimental data and the CFD simulation results. The ECM and PECM methods are applied to predict thermal loads to the vessel wall and Control Rod Guide Tubes (CRGTs) during core debris heat up and melting in the BWR lower plenum. Applying the ECM and PECM to simulations of reactor-scale melt pool heat transfer, the results of the ECM and PECM calculations show an apparent effectiveness of the developed methods that enables simulations of long term accident transients. It is also found that during severe accident progression, the cooling by water flowing inside the CRGTs plays a very important role in reducing the thermal load on the reactor vessel wall. The results of the CFD, ECM and PECM simulations suggest a potential of the CRGT cooling as an effective mitigative measure during a severe accident progression. / QC 20101119

[ES] La tecnología nuclear para el uso civil genera más preocupación por la seguridad que muchas otras tecnologías que se usan a diario. La Autoridad Nuclear define las bases de cómo debe realizarse la operación segura de una Central Nuclear. De acuerdo a las directrices establecidas por la Autoridad Nuclear, una Central Nuclear debe analizar una envolvente de escenarios hipotéticos y comprobar de manera determinista que los criterios de aceptación para dicho evento se cumplen. El Análisis Determinista de Seguridad utiliza herramientas de simulación que aplican la física conocida sobre el comportamiento de la Central Nuclear para evaluar la evolución de una variable de seguridad y asegurar que los límites no se sobrepasan. El desarrollo de la tecnología informática, de los métodos matemáticos y de la física que envuelve el comportamiento de una Central Nuclear han proporcionado herra-mientas de simulación potentes que son capaces de predecir el comportamiento de las variables de seguridad con una importante precisión. Esto permite analizar escenarios de manera más realista evitando asumir condiciones conservadoras que hasta la fecha compensaban la falta de conocimiento modelado en las herramientas de simulación. Las herramientas conocidas como De Mejor Estimación son capaces de analizar even-tos transitorios en diferentes escalas. Además, emplean modelos analíticos de las dife-rentes físicas más detallados, así como correlaciones experimentales más realistas y actuales. Un paso adelante en el Análisis Determinista de Seguridad pretende combinar las diferentes herramientas de Mejor Estimación que se emplean para analizar las dis-tintas físicas de una Central Nuclear, considerando incluso la interacción entre ellas y el análisis progresivo a diferentes escalas, llegando a analizar fenómenos más locales si es necesario. Para este fin, esta tesis presenta una metodología de análisis multi-físico y multi-escala que emplea diferentes códigos de simulación analizando el escenario propuesto a dife-rentes escalas, es decir, desde un nivel de planta que incluye los distintos componentes, hasta el volumen de control que supone el refrigerante pasando entre las varillas de combustible. Esta metodología permite un flujo de información que va desde el análi-sis a mayor escala hasta el de menor escala. El desarrollo de esta metodología ha sido validado con datos de planta para poder evaluar el alcance de esta metodología y pro-porcionar nuevas líneas de trabajo futuro. Además, se han añadido los resultados de los distintos procesos de validación y verificación que han surgido a lo largo de este trabajo. / [CA] La tecnologia nuclear per a l'ús civil genera més preocupació per la seguretat que moltes altres tecnologies d'ús quotidià. L'Autoritat Nuclear defineix les bases de com ha de realitzar-se l'operació segura d'una Central Nuclear. D'acord amb les directrius establertes per l'Autoritat Nuclear, una Central Nuclear ha d'analitzar una envoltant d'escenaris hipotètics I comprovar de manera determinista que els criteris d'acceptació per a l'esdeveniment seleccionat es compleixen. L'Anàlisi Determinista de Seguretat utilitza eines de simulació que apliquen la física coneguda sobre el comportament de la Central Nuclear per avaluar l'evolució d'una variable de seguretat i assegurar que els límits no es traspassen. El desenvolupament de la tecnologia informàtica, els mètodes matemàtics i de la física que envolta el comportament d'una Central Nuclear han proporcionat eines de simulació potents amb capacitat de predir el comportament de les variables de seguretat amb una precisió significativa. Això permet analitzar escenaris de manera realista evitant assumir condicions conservadores que fins al moment compensaven la mancança de coneixement. Les eines de simulació conegudes com De Millor Estimació son capaces d'analitzar esdeveniment transitoris a diferent escales. A més, utilitzen models analítics per a les diferents físiques amb més detall així com correlacions experimentals més actualitzades i realistes. Un pas més endavant en l'Anàlisi Determinista de Seguretat pretén combinar les diferents eines de Millor Estimació que se utilitzen per analitzar les distintes físiques d'una Central Nuclear, considerant inclús la interacció entre ells i l'anàlisi progressiu a diferents escales, amb la finalitat de poder analitzar fenòmens locals. Per a aquest fi, esta tesi presenta una metodologia d'anàlisi multi-física i multi-escala que utilitza diferents codis de simulació analitzant l'escenari proposat a diferents escales, és a dir, des d'un nivell de planta que inclou els distints components, fins al volum de control que suposa el refrigerant passant entre les varetes de combustible. Esta metodologia permet un flux de informació que va des de l'anàlisi d'una escala major a una menor. El desenvolupament d'aquesta metodologia ha sigut validada i verificada amb dades de planta i els resultats han sigut analitzats a fi d'avaluar la capacitat de la metodologia i les possibles línies de treball futur. A més s'han afegit els principals resultats de verificació i validació que han sorgit en les distintes etapes d'aquest treball. / [EN] The nuclear technology for civil use has generated more concerns for the safety than several other technologies applied to the daily life. The Nuclear Regulators define the basis of how the Safety Operation of Nuclear Power Plants is to be done. According to these guidelines, a Nuclear Power Plant must analyze an envelope of hypothetical events and deterministically define if the acceptance criteria for these events is met. The Deterministic Safety Analysis uses simulation tools that apply the physics known in the behavior of the Nuclear Power Plant to evaluate the evolution of a safety varia-ble and assure that the safety limits will not be exceeded. The development of the computer science, the numerical methods and the physics involved in the behavior of a Nuclear Power Plant have yield powerful simulation tools that are capable to predict the evolution of safety variables which significant accuracy. This allows to consider more realistic simulation scenarios instead of con-servative approaches in order to compensate the lack of knowledge in the applied prediction methods. The so called Best Estimate simulation tools are capable to analyze the transient events in different scales. Furthermore, they account more detailed analytical models and experimental correlations. A step forward in the Deterministic Safety Analysis intends to combine the Best Estimate simulation tools of the different physics considering the interaction among them and analyzing the different scales, considering more local approaches if necessary. For this purpose, this thesis work presents a multi-scale and multi-physics methodology that uses different physics codes and has the aim of modeling postulated scenarios in different scales, i.e. from system models representing the components of the plants to the subchannel models that analyze the behavior of the coolant between the fuel rods. This methodology allows a flow of information where the output of one scale is used as input in a more detailed scale to predict a more local analysis of parameters, such as the Critical Power Ratio, which are of great importance for the estimation of safety margins. The development of this methodology has been validated against plant data with the aim of evaluating the scope of this methodology and in order to provide future lines of development. In addition, different results of the validation and verifi-cation yielded in the development of the parts of this methodology are presented. / Hidalga García-Bermejo, P. (2020). Development and validation of a multi-scale and multi-physics methodology for the safety analysis of fast transients in Light Water Reactors [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/160135

Considering the late in-vessel phase of an unlikely core melt down scenario in a light water reactor (LWR) with the formation of a corium pool in the lower head of the reactor pressure vessel (RPV) the possible failure modes of the RPV and the time to failure have to be investigated to assess the possible loadings on the containment. In this work, an integral model was developed to describe the processes in the lower plenum of the RPV. Two principal model domains have to be distinguished: The tem-perature field within the melt and RPV is calculated with a thermodynamic model, while a mechanical model is used for the structural analysis of the vessel wall. In the introducing chapters a description is given of the considered accident scenario and the relevant analytical, experimental, and numerical investigations are discussed which were performed worldwide during the last three decades. Following, the occurring physical phenomena are analysed and the scaling differences are evaluated between the FOREVER-experiments and a prototypical scenario. The thermodynamic and the mechanical model can be coupled recursively to take into account the mutual influence. This approach not only allows to consider the temperature dependence of the material parameters and the thermally induced stress in the mechanical model, it also takes into account the response of the temperature field itself upon the changing vessel geometry. New approaches are applied in this work for the simulation of creep and damage. Using a creep data base, the application of single creep laws could be avoided which is especially advantageous if large temperature, stress, and strain ranges have to be covered. Based on experimental investigations, the creep data base has been developed for an RPV-steel and has been validated against creep tests with different scalings and geometries. It can be stated, that the coupled model is able to exactly describe and predict the vessel deformation in the scaled integral FOREVER-tests. There are uncertainties concerning the time to failure which are related to inexactly known material parameters and boundary conditions. The main results of this work can be summarised as follows: Due to the thermody-namic behaviour of the large melt pool with internal heat sources, the upper third of the lower RPV head is exposed to the highest thermo-mechanical loads. This region is called hot focus. Contrary to that, the pole part of the lower head has a higher strength and therefore relocates almost vertically downwards under the combined thermal, weight and internal pressure load of the RPV. On the one hand, it will be possible by external flooding to retain the corium within the RPV even at increased pressures and even in reactors with high power (as e.g. KONVOI). On the other hand, there is no chance for melt retention in the considered scenario if neither internal nor external flooding of the RPV can be achieved. Two patents have been derived from the gained insights. Both are related to passively working devices for accident mitigation: The first one is a support of the RPV lower head pole part. It reduces the maximum mechanical load in the highly stressed area of the hot focus. In this way, it can prevent failure or at least extend the time to failure of the vessel. The second device implements a passive accident mitigation measure by making use of the downward movement of the lower head. Through this, a valve or a flap can be opened to flood the reactor pit with water from a storage reservoir located at a higher position in the reactor building. With regard to future plant designs it can be stated - differing from former presumptions - that an In-Vessel-Retention (IVR) of a molten core is possible within the reactor pressure vessel even for reactors with higher power. / Für das unwahrscheinliche Szenario eines Kernschmelzunfalls in einem Leichtwasserreaktor mit Bildung eines Schmelzesees in der Bodenkalotte des Reaktordruckbehälters (RDB) ist es notwendig, mögliche Versagensformen des RDB sowie Versagenszeiträume zu ermitteln, um die daraus resultierende mögliche Belastung des Sicherheitsbehälters bestimmen zu können. In dieser Arbeit wird ein integrales Modell entwickelt, das die Vorgänge im unteren Plenum beschreibt. Dabei sind zwei prinzipielle Modellbereiche zu unterscheiden: Das Temperaturfeld in der Schmelze und im RDB wird mit einem thermodynamischen Modell berechnet, während für die Strukturanalyse des RDB ein mechanisches Modell verwendet wird. Zunächst werden das betrachtete Unfallszenario dargestellt und die bisher in den letzten drei Dekaden weltweit durchgeführten wesentlichen analytischen, experimentellen und numerischen Untersuchungen diskutiert. Anschließend werden die auftretenden physikalischen Vorgänge analysiert. Gleichzeitig werden Skalierungsunterschiede zwischen den in dieser Arbeit betrachteten Experimenten der FOREVER-Reihe und einem prototypischen Szenario herausgearbeitet. Das thermodynamische und das mechanische Modell können rekursiv gekoppelt werden, wodurch die wechselseitige Beeinflussung berücksichtigt werden kann. Insbesondere werden damit neben der Temperaturabhängigkeit der Materialparameter und den thermisch induzierten Spannungen im mechanischen Modell auch die Rückwirkungen der Behälterverformung auf das Temperaturfeld selber erfasst. Für die Kriech- und Schädigungssimulation werden in dieser Arbeit neue Verfahren angewendet. Durch die Entwicklung und den Einsatz einer Kriechdatenbasis konnte die bei sehr unterschiedlichen Temperaturen, Spannungen und Dehnungen ungeeignete Verwendung einzelner Kriechgesetze umgangen werden. Aufbauend auf experimentellen Untersuchungen wurde eine Kriechdatenbasis für einen RDB-Stahl entwickelt und an Hand von Kriechversuchen verschiedener Geometrie und Dimension validiert. Als Ergebnis lässt sich festhalten, dass das gekoppelte Modell prinzipiell in der Lage ist, die Behälterdeformation im Falle der skalierten FOREVER-Experimente exakt zu beschreiben bzw. vorherzusagen. Unsicherheiten bezüglich der Versagenszeit resultieren aus nicht exakt bekannten Materialparametern und Randbedingungen. Die wesentlichen Ergebnisse dieser Arbeit lassen sich wie folgt zusammenfassen: Aufgrund des thermodynamischen Verhaltens eines großen Schmelzesees mit inneren Wärmequellen erfolgt die höchste thermomechanische Belastung des RDB im oberen Drittel der Bodenkalotte. Dieser Bereich wird als heißer Fokus bezeichnet. Der untere Bereich der Kalotte weist hingegen eine höhere Festigkeit auf und verlagert sich deswegen bei entsprechender Belastung des RDB im wesentlichen senkrecht nach unten. Bei einer externen Flutung besteht auch bei hohen Innendrücken für einen Reaktor großer Leistung (KONVOI) die Möglichkeit, die Schmelze im RDB zurückzuhalten. Ohne interne oder externe Flutung besteht für das betrachtete Szenario keine Aussicht für eine Schmelzerückhaltung im RDB. Aus den gewonnenen Erkenntnissen wurden zwei Patente abgeleitet. Dabei handelt es sich um passiv wirkende Einrichtungen zur Schadensbegrenzung: Die erste reduziert durch Abstützen des unteren Kalottenzentrums die Maximalspannungen im hochbeanspruchten Bereich des heißen Fokus und kann damit ein Versagen verhindern oder zumindest verzögern. Die zweite Einrichtung ermöglicht die passive Auslösung einer Flutung, indem die Abwärtsbewegung der Kalotte zur Steuerung genutzt wird. Hierdurch kann beispielsweise ein Ventil geöffnet werden, um Wasser aus im Gebäude höher angeordneten Reservoirs in die Reaktorgrube zu leiten. Abweichend von bisherigen Annahmen kann im Hinblick auf die Entwicklung zukünftiger Baulinien festgehalten werden, dass eine Kernschmelzerückhaltung im Reaktordruckbehälter auch für Reaktoren größerer Leistung möglich ist.

Considering the late in-vessel phase of an unlikely core melt down scenario in a light water reactor (LWR) with the formation of a corium pool in the lower head of the reactor pressure vessel (RPV) the possible failure modes of the RPV and the time to failure have to be investigated to assess the possible loadings on the containment. In this work, an integral model was developed to describe the processes in the lower plenum of the RPV. Two principal model domains have to be distinguished: The tem-perature field within the melt and RPV is calculated with a thermodynamic model, while a mechanical model is used for the structural analysis of the vessel wall. In the introducing chapters a description is given of the considered accident scenario and the relevant analytical, experimental, and numerical investigations are discussed which were performed worldwide during the last three decades. Following, the occurring physical phenomena are analysed and the scaling differences are evaluated between the FOREVER-experiments and a prototypical scenario. The thermodynamic and the mechanical model can be coupled recursively to take into account the mutual influence. This approach not only allows to consider the temperature dependence of the material parameters and the thermally induced stress in the mechanical model, it also takes into account the response of the temperature field itself upon the changing vessel geometry. New approaches are applied in this work for the simulation of creep and damage. Using a creep data base, the application of single creep laws could be avoided which is especially advantageous if large temperature, stress, and strain ranges have to be covered. Based on experimental investigations, the creep data base has been developed for an RPV-steel and has been validated against creep tests with different scalings and geometries. It can be stated, that the coupled model is able to exactly describe and predict the vessel deformation in the scaled integral FOREVER-tests. There are uncertainties concerning the time to failure which are related to inexactly known material parameters and boundary conditions. The main results of this work can be summarised as follows: Due to the thermody-namic behaviour of the large melt pool with internal heat sources, the upper third of the lower RPV head is exposed to the highest thermo-mechanical loads. This region is called hot focus. Contrary to that, the pole part of the lower head has a higher strength and therefore relocates almost vertically downwards under the combined thermal, weight and internal pressure load of the RPV. On the one hand, it will be possible by external flooding to retain the corium within the RPV even at increased pressures and even in reactors with high power (as e.g. KONVOI). On the other hand, there is no chance for melt retention in the considered scenario if neither internal nor external flooding of the RPV can be achieved. Two patents have been derived from the gained insights. Both are related to passively working devices for accident mitigation: The first one is a support of the RPV lower head pole part. It reduces the maximum mechanical load in the highly stressed area of the hot focus. In this way, it can prevent failure or at least extend the time to failure of the vessel. The second device implements a passive accident mitigation measure by making use of the downward movement of the lower head. Through this, a valve or a flap can be opened to flood the reactor pit with water from a storage reservoir located at a higher position in the reactor building. With regard to future plant designs it can be stated - differing from former presumptions - that an In-Vessel-Retention (IVR) of a molten core is possible within the reactor pressure vessel even for reactors with higher power. / Für das unwahrscheinliche Szenario eines Kernschmelzunfalls in einem Leichtwasserreaktor mit Bildung eines Schmelzesees in der Bodenkalotte des Reaktordruckbehälters (RDB) ist es notwendig, mögliche Versagensformen des RDB sowie Versagenszeiträume zu ermitteln, um die daraus resultierende mögliche Belastung des Sicherheitsbehälters bestimmen zu können. In dieser Arbeit wird ein integrales Modell entwickelt, das die Vorgänge im unteren Plenum beschreibt. Dabei sind zwei prinzipielle Modellbereiche zu unterscheiden: Das Temperaturfeld in der Schmelze und im RDB wird mit einem thermodynamischen Modell berechnet, während für die Strukturanalyse des RDB ein mechanisches Modell verwendet wird. Zunächst werden das betrachtete Unfallszenario dargestellt und die bisher in den letzten drei Dekaden weltweit durchgeführten wesentlichen analytischen, experimentellen und numerischen Untersuchungen diskutiert. Anschließend werden die auftretenden physikalischen Vorgänge analysiert. Gleichzeitig werden Skalierungsunterschiede zwischen den in dieser Arbeit betrachteten Experimenten der FOREVER-Reihe und einem prototypischen Szenario herausgearbeitet. Das thermodynamische und das mechanische Modell können rekursiv gekoppelt werden, wodurch die wechselseitige Beeinflussung berücksichtigt werden kann. Insbesondere werden damit neben der Temperaturabhängigkeit der Materialparameter und den thermisch induzierten Spannungen im mechanischen Modell auch die Rückwirkungen der Behälterverformung auf das Temperaturfeld selber erfasst. Für die Kriech- und Schädigungssimulation werden in dieser Arbeit neue Verfahren angewendet. Durch die Entwicklung und den Einsatz einer Kriechdatenbasis konnte die bei sehr unterschiedlichen Temperaturen, Spannungen und Dehnungen ungeeignete Verwendung einzelner Kriechgesetze umgangen werden. Aufbauend auf experimentellen Untersuchungen wurde eine Kriechdatenbasis für einen RDB-Stahl entwickelt und an Hand von Kriechversuchen verschiedener Geometrie und Dimension validiert. Als Ergebnis lässt sich festhalten, dass das gekoppelte Modell prinzipiell in der Lage ist, die Behälterdeformation im Falle der skalierten FOREVER-Experimente exakt zu beschreiben bzw. vorherzusagen. Unsicherheiten bezüglich der Versagenszeit resultieren aus nicht exakt bekannten Materialparametern und Randbedingungen. Die wesentlichen Ergebnisse dieser Arbeit lassen sich wie folgt zusammenfassen: Aufgrund des thermodynamischen Verhaltens eines großen Schmelzesees mit inneren Wärmequellen erfolgt die höchste thermomechanische Belastung des RDB im oberen Drittel der Bodenkalotte. Dieser Bereich wird als heißer Fokus bezeichnet. Der untere Bereich der Kalotte weist hingegen eine höhere Festigkeit auf und verlagert sich deswegen bei entsprechender Belastung des RDB im wesentlichen senkrecht nach unten. Bei einer externen Flutung besteht auch bei hohen Innendrücken für einen Reaktor großer Leistung (KONVOI) die Möglichkeit, die Schmelze im RDB zurückzuhalten. Ohne interne oder externe Flutung besteht für das betrachtete Szenario keine Aussicht für eine Schmelzerückhaltung im RDB. Aus den gewonnenen Erkenntnissen wurden zwei Patente abgeleitet. Dabei handelt es sich um passiv wirkende Einrichtungen zur Schadensbegrenzung: Die erste reduziert durch Abstützen des unteren Kalottenzentrums die Maximalspannungen im hochbeanspruchten Bereich des heißen Fokus und kann damit ein Versagen verhindern oder zumindest verzögern. Die zweite Einrichtung ermöglicht die passive Auslösung einer Flutung, indem die Abwärtsbewegung der Kalotte zur Steuerung genutzt wird. Hierdurch kann beispielsweise ein Ventil geöffnet werden, um Wasser aus im Gebäude höher angeordneten Reservoirs in die Reaktorgrube zu leiten. Abweichend von bisherigen Annahmen kann im Hinblick auf die Entwicklung zukünftiger Baulinien festgehalten werden, dass eine Kernschmelzerückhaltung im Reaktordruckbehälter auch für Reaktoren größerer Leistung möglich ist.

The nuclear industry has long relied upon bounding parametric analyses in predicting the safety margins of reactor designs undergoing design-basis accidents. These methods have been known to return highly-conservative results, limiting the operating conditions of the reactor. The Best-Estimate Plus Uncertainty (BEPU) method using a modernized version of the Code-Scaling, Applicability, and Uncertainty (CSAU) methodology has been applied to more accurately predict the safety margins of the Oregon State University Advanced Plant Experiment (APEX) facility experiencing a Loss-of-Feedwater Accident (LOFA). The statistical advantages of the Bayesian paradigm of probability was utilized to incorporate prior knowledge when determining the analysis required to justify the safety margins. RELAP5 Mod 3.3 was used to accurately predict the thermal-hydraulics of a primary Feed-and-Bleed response to the accident using assumptions to accompany the lumped-parameter calculation approach. A novel coupling of thermal-hydraulic and statistical software was accomplished using the Symbolic Nuclear Analysis Package (SNAP). Uncertainty in Peak Cladding Temperature (PCT) was calculated at the 95/95 probability/confidence levels under a series of four separate sensitivity studies. / Graduation date: 2013