The best estimate codes are used for licensing of Nuclear Power Plants (NPP), but with conservative assumptions. It is claimed that the uncertainties are covered by the conservatism of the calculation. Nowadays, it is possible to estimate certain parameters using non-conservative data with the complement of uncertainty evaluation, and these calculations can also be used for licensing. As NPPs are applying for power up-rates and life extension, new licensing calculations need to be performed. In this case, evaluation of the uncertainties could help improve the performance, while staying below the limit of the safety margins. Given the problem of unstable behavior of Boiling Water Reactors (BWR), which is known to occur at certain power and flow conditions, it could cause SCRAM and decrease the economic performance of the plant. Performing an uncertainty analysis for BWR stability would give better understating of the phenomenon and it would help to verify and validate (V&amp;V) the codes used to predict the NPP behavior. This thesis, reports a sensitivity/uncertainty study of numerical, neutronics, and thermal-hydraulics parameters on the prediction of the BWR stability within the framework of OECD Ringhals-1 (R1 stable reactor) and OECD Oskarshamn-2 (O2 unstable reactor) stability benchmarks. The time domain code TRACE/PARCS was used in the analyses. This thesis is divided in three parts: space-time convergence; uncertainty; sensitivity. A space-time convergence study was done for the numerical parameters (nodalization and time step). This was done by refining nodalization of all components and time step until obtaining space-time converged solution, i.e. further refinement doesn’t change the solution. When the space-time converged solutions were compared to the initial models, much better solution accuracy has been obtained for the stability measures (decay ratio and frequency), for both stable (R1) and unstable (O2) reactors with the space-time converged models. Further on, important neutronics and thermal-hydraulics parameters were identified and an uncertainty calculation was performed using the Propagation of Input Errors (PIE) methodology. This methodology, also known as the GRS method, has been used because it has been extensively tested and verified by the industry, and because it allows identifying the most influential parameters using the spearman rank correlation method. Using the uncertainty method’s results, an attempt has been done to identify the most influential parameters affecting the stability. A methodology using the spearman rank correlation coefficient has been implemented, which helps to identify the most influential parameters on the stability (decay ratio and frequency). Additional sensitivity calculations have been performed for better understanding of BWR stability and parameters that affect it. / <p>This work has been preformed thanks to the support of the Swedish Radiation Safety Authority (SSM) and EU project NURISP. QC 20121129</p>

Best Estimate codes are used for licensing, but with conservative assumptions. It is claimed that the uncertainties are covered by the conservatism of the calculation. As Nuclear Power Plants are applying for power up-rates and life extension, evaluation of the uncertainties could help improve the performance, while staying below the limit of the safety margins.   Given the problem of unstable behavior of Boiling Water Reactors (BWRs), which is known to occur during operation at certain power and flow conditions, it could cause SCRAM and decrease the economic performance of the plant. Performing an uncertainty analysis for BWR stability would give better understating of the phenomenon and it would help to verify and validate (V&amp;V) the codes used to predict the NPP behavior.   This thesis reports an uncertainty study of the impact of Thermal-Hydraulic, Neutronic, and Numerical parameters on the prediction of the stability of the BWR within the framework of OECD Ringhals-1 stability benchmark. The time domain code TRACE/PARCS was used in the analysis. This thesis is divided in two parts: Sensitivity study on Numerical Discretization Parameters (Nodalization, Time Step, etc.) and Uncertainty part.   A Sensitivity study was done for the Numerical Parameters (Nodalization and Time step). This was done by refining all possible components until obtaining Space-Time Converged Solution, i.e. further refinement doesn’t change the solution. When the space-time converged solution was compared to the initial discretization, a much better solution has been obtained for both the stability measures (Decay Ratio and Frequency) with the space-time converged model.   Further on, important Neutronic and Thermal-Hydraulic Parameters were identified and the uncertainty calculation was performed using the Propagation of Input Errors (PIE) methodology. This methodology, also known as the GRS method, has been used because it has been tested and extensively verified by the industry, and because it allows identifying the most influential parameters using the Spearman Rank Correlation. / QC 20101126

This thesis is concerned with nonlinear analyses of BWR stability behaviour, contributing to a deeper understanding in this field. Despite negative feedback-coefficients of a BWR, there are operational points (OP) at which oscillatory instabilities occur. So far, a comprehensive and an in-depth understanding of the nonlinear BWR stability behaviour are missing, even though the impact of the significant physical parameters is well known. In particular, this concerns parameter regions in which linear stability indicators, like the asymptotic decay ratio, lose their meaning. Nonlinear stability analyses are usually carried out using integral (system) codes, describing the dynamical system by a system of nonlinear partial differential equations (PDE). One aspect of nonlinear BWR stability analyses is to get an overview about different types of nonlinear stability behaviour and to examine the conditions of their occurrence. For these studies the application of system codes alone is inappropriate. Hence, in the context of this thesis, a novel approach to nonlinear BWR stability analyses, called RAM-ROM method, is developed. In the framework of this approach, system codes and reduced order models (ROM) are used as complementary tools to examine the stability characteristics of fixed points and periodic solutions of the system of nonlinear differential equations, describing the stability behaviour of a BWR loop. The main advantage of a ROM, which is a system of ordinary differential equations (ODE), is the possible coupling with specific methods of the nonlinear dynamics. This method reveals nonlinear phenomena in certain regions of system parameters without the need for solving the system of ROM equations. The stability properties of limit cycles generated in Hopf bifurcation points and the conditions of their occurrence are of particular interest. Finally, the nonlinear phenomena predicted by the ROM will be analysed in more details by the system code. Hence, the thesis is not focused on rendering more precisely linear stability indicators like DR. The objective of the ROM development is to develop a model as simple as possible from the mathematical and numerical point of view, while preserving the physics of the BWR stability behaviour. The ODEs of the ROM are deduced from the PDEs describing the dynamics of a BWR. The system of ODEs includes all spatial effects in an approximated (spatial averaged) manner, e.g. the space-time dependent neutron flux is expanded in terms of a complete set of orthogonal spatial neutron flux modes. In order to simulate the stability characteristics of the in-phase and out-of-phase oscillation mode, it is only necessary to take into account the fundamental mode and the first azimuthal mode. The ROM, originally developed at PSI in collaboration with the University of Illinois (PSI-Illinois-ROM), was upgraded in significant points: • Development and implementation of a new calculation methodology for the mode feedback reactivity coefficients (void and fuel temperature reactivity) • Development and implementation of a recirculation loop model; analysis and discussion of its impact on the in-phase and out-of-phase oscillation mode • Development of a novel physically justified approach for the calculation of the ROM input data • Discussion of the necessity of consideration of the effect of subcooled boiling in an approximate manner With the upgraded ROM, nonlinear BWR stability analyses are performed for three OPs (one for NPP Leibstadt (cycle7), one for NPP Ringhals (cycle14) and one for NPP Brunsbüttel (cycle16) for which measuring data of stability tests are available. In this thesis, the novel approach to nonlinear BWR stability analyses is extensively presented for NPP Leibstadt. In particular, the nonlinear analysis is carried out for an operational point (OP), in which an out-of-phase power oscillation has been observed in the scope of a stability test at the beginning of cycle 7 (KKLc7_rec4). The ROM predicts a saddle-node bifurcation of cycles, occurring in the linear stable region, close to the KKLc7_rec4-OP. This result allows a new interpretation of the stability behaviour around the KKLc7_rec4-OP. The results of this thesis confirm that the RAM-ROM methodology is qualified for nonlinear BWR stability analyses. / Die vorliegende Dissertation leistet einen Beitrag zum tieferen Verständnis des nichtlinearen Stabilitätsverhaltens von Siedewasserreaktoren (SWR). Trotz der Tatsache, dass in diesem technischen System nur negative innere Rückkopplungskoeffizienten auftreten, können in bestimmten Arbeitspunkten oszillatorische Instabilitäten auftreten. Obwohl relativ gute Kenntnisse über die signifikanten physikalischen Einflussgrößen vorliegen, fehlt bisher ein umfassendes Verständnis des SWR-Stabilitätsverhaltens. Das betrifft insbesondere die Bereiche der Systemparameter, in denen lineare Stabilitätsindikatoren, wie zum Beispiel das asymptotische Decay Ratio (DR), ihren Sinn verlieren. Die nichtlineare Stabilitätsanalyse wird im Allgemeinen mit Systemcodes (nichtlineare partielle Differentialgleichungen, PDG) durchgeführt. Jedoch kann mit Systemcodes kein oder nur ein sehr lückenhafter Überblick über die Typen von nichtlinearen Phänomenen, die in bestimmten System-Parameterbereichen auftreten, erhalten werden. Deshalb wurde im Rahmen der vorliegenden Arbeit eine neuartige Methode (RAM-ROM Methode) zur nichtlinearen SWR-Stabilitätsanalyse erprobt, bei der integrale Systemcodes und sog. vereinfachte SWR-Modelle (ROM) als sich gegenseitig ergänzende Methoden eingesetzt werden, um die Stabilitätseigenschaften von Fixpunkten und periodischen Lösungen (Grenzzyklen) des nichtlinearen Differentialgleichungssystems, welches das Stabilitätsverhalten des SWR beschreibt, zu bestimmen. Das ROM, in denen das dynamische System durch gewöhnliche Differentialgleichungen (GDG) beschrieben wird, kann relativ einfach mit leistungsfähigen Methoden aus der nichtlinearen Dynamik, wie zum Beispiel die semianalytische Bifurkationsanalyse, gekoppelt werden. Mit solchen Verfahren kann, ohne das DG-System explizit lösen zu müssen, ein Überblick über mögliche Typen von stabilen und instabilen oszillatorischen Verhalten des SWR erhalten werden. Insbesondere sind die Stabilitätseigenschaften von Grenzzyklen, die in Hopf-Bifurkationspunkten entstehen, und die Bedingungen, unter denen sie auftreten, von Interesse. Mit dem Systemcode (RAMONA5) werden dann die mit dem ROM vorhergesagten Phänomene in den entsprechenden Parameterbereichen detaillierter untersucht (Validierung des ROM). Die Methodik dient daher nicht der Verfeinerung der Berechnung linearer Stabilitätsindikatoren (wie das DR). Das ROM-Gleichungssystem entsteht aus den PDGs des Systemcodes durch geeignete (nichttriviale) räumliche Mittelung der PDG. Es wird davon ausgegangen, dass die Reduzierung der räumlichen Komplexität die Stabilitätseigenschaften des SWR nicht signifikant verfälschen, da durch geeignete Mittlungsverfahren, räumliche Effekte näherungsweise in den GDGs berücksichtig werden. Beispielsweise wird die raum- und zeitabhängige Neutronenflussdichte nach räumlichen Moden entwickelt, wobei für eine Simulation der Stabilitätseigenschaften der In-phase- und Out-of-Phase-Leistungsoszillationen nur der Fundamentalmode und der erste azimuthale Mode berücksichtigt werden muss. Das ROM, welches ursprünglich am Paul Scherrer Institut (PSI, Schweiz) in Zusammenarbeit mit der Universität Illinois (USA) entwickelt wurde, ist in zwei wesentlichen Punkten erweitert und verbessert worden: • Entwicklung und Implementierung einer neuen Methode zur Berechnung der Rückkopplungsreaktivitäten • Entwicklung und Implementierung eines Modells zur Beschreibung der Rezirkulationsschleife (insbesondere wurde der Einfluss der Rezirkulationsschleife auf den In-Phase-Oszillationszustand und auf den Out-of-Phase-Oszillationszustand untersucht) • Entwicklung einer physikalisch begründeten Methode zur Berechnung der ROM-Inputdaten • Abschätzung des Einflusses des unterkühlten Siedens im Rahmen der ROM-Näherungen Mit dem erweiterten ROM wurden nichtlineare Stabilitätsanalysen für drei Arbeitspunkte (KKW Leibstadt (Zyklus 7) KKW Ringhals (Zyklus 14) und KKW Brunsbüttel (Zyklus 16)), für die Messdaten vorliegen, durchgeführt. In der Dissertationsschrift wird die RAM-ROM Methode ausführlich am Beispiel eines Arbeitspunktes (OP) des KKW Leibstadt (KKLc7_rec4-OP), in dem eine aufklingende regionale Leistungsoszillation bei einem Stabilitätstest gemessen worden ist, demonstriert. Das ROM sagt die Existenz eines Umkehrpunktes (saddle-node bifurcation of cycles, fold-bifurcation) voraus, der sich im linear stabilen Gebiet nahe der Stabilitätsgrenze befindet. Mit diesem ROM-Ergebnis ist eine neue Interpretation der Stabilitätseigenschaften des KKLc7_rec4-OP möglich. Die Resultate der in der Dissertation durchgeführten RAM-ROM Analyse bestätigen, dass das weiterentwickelte ROM für die Analyse des Stabilitätsverhaltens realer Leistungsreaktoren qualifiziert wurde.

This thesis is concerned with nonlinear analyses of BWR stability behaviour, contributing to a deeper understanding in this field. Despite negative feedback-coefficients of a BWR, there are operational points (OP) at which oscillatory instabilities occur. So far, a comprehensive and an in-depth understanding of the nonlinear BWR stability behaviour are missing, even though the impact of the significant physical parameters is well known. In particular, this concerns parameter regions in which linear stability indicators, like the asymptotic decay ratio, lose their meaning. Nonlinear stability analyses are usually carried out using integral (system) codes, describing the dynamical system by a system of nonlinear partial differential equations (PDE). One aspect of nonlinear BWR stability analyses is to get an overview about different types of nonlinear stability behaviour and to examine the conditions of their occurrence. For these studies the application of system codes alone is inappropriate. Hence, in the context of this thesis, a novel approach to nonlinear BWR stability analyses, called RAM-ROM method, is developed. In the framework of this approach, system codes and reduced order models (ROM) are used as complementary tools to examine the stability characteristics of fixed points and periodic solutions of the system of nonlinear differential equations, describing the stability behaviour of a BWR loop. The main advantage of a ROM, which is a system of ordinary differential equations (ODE), is the possible coupling with specific methods of the nonlinear dynamics. This method reveals nonlinear phenomena in certain regions of system parameters without the need for solving the system of ROM equations. The stability properties of limit cycles generated in Hopf bifurcation points and the conditions of their occurrence are of particular interest. Finally, the nonlinear phenomena predicted by the ROM will be analysed in more details by the system code. Hence, the thesis is not focused on rendering more precisely linear stability indicators like DR. The objective of the ROM development is to develop a model as simple as possible from the mathematical and numerical point of view, while preserving the physics of the BWR stability behaviour. The ODEs of the ROM are deduced from the PDEs describing the dynamics of a BWR. The system of ODEs includes all spatial effects in an approximated (spatial averaged) manner, e.g. the space-time dependent neutron flux is expanded in terms of a complete set of orthogonal spatial neutron flux modes. In order to simulate the stability characteristics of the in-phase and out-of-phase oscillation mode, it is only necessary to take into account the fundamental mode and the first azimuthal mode. The ROM, originally developed at PSI in collaboration with the University of Illinois (PSI-Illinois-ROM), was upgraded in significant points: • Development and implementation of a new calculation methodology for the mode feedback reactivity coefficients (void and fuel temperature reactivity) • Development and implementation of a recirculation loop model; analysis and discussion of its impact on the in-phase and out-of-phase oscillation mode • Development of a novel physically justified approach for the calculation of the ROM input data • Discussion of the necessity of consideration of the effect of subcooled boiling in an approximate manner With the upgraded ROM, nonlinear BWR stability analyses are performed for three OPs (one for NPP Leibstadt (cycle7), one for NPP Ringhals (cycle14) and one for NPP Brunsbüttel (cycle16) for which measuring data of stability tests are available. In this thesis, the novel approach to nonlinear BWR stability analyses is extensively presented for NPP Leibstadt. In particular, the nonlinear analysis is carried out for an operational point (OP), in which an out-of-phase power oscillation has been observed in the scope of a stability test at the beginning of cycle 7 (KKLc7_rec4). The ROM predicts a saddle-node bifurcation of cycles, occurring in the linear stable region, close to the KKLc7_rec4-OP. This result allows a new interpretation of the stability behaviour around the KKLc7_rec4-OP. The results of this thesis confirm that the RAM-ROM methodology is qualified for nonlinear BWR stability analyses. / Die vorliegende Dissertation leistet einen Beitrag zum tieferen Verständnis des nichtlinearen Stabilitätsverhaltens von Siedewasserreaktoren (SWR). Trotz der Tatsache, dass in diesem technischen System nur negative innere Rückkopplungskoeffizienten auftreten, können in bestimmten Arbeitspunkten oszillatorische Instabilitäten auftreten. Obwohl relativ gute Kenntnisse über die signifikanten physikalischen Einflussgrößen vorliegen, fehlt bisher ein umfassendes Verständnis des SWR-Stabilitätsverhaltens. Das betrifft insbesondere die Bereiche der Systemparameter, in denen lineare Stabilitätsindikatoren, wie zum Beispiel das asymptotische Decay Ratio (DR), ihren Sinn verlieren. Die nichtlineare Stabilitätsanalyse wird im Allgemeinen mit Systemcodes (nichtlineare partielle Differentialgleichungen, PDG) durchgeführt. Jedoch kann mit Systemcodes kein oder nur ein sehr lückenhafter Überblick über die Typen von nichtlinearen Phänomenen, die in bestimmten System-Parameterbereichen auftreten, erhalten werden. Deshalb wurde im Rahmen der vorliegenden Arbeit eine neuartige Methode (RAM-ROM Methode) zur nichtlinearen SWR-Stabilitätsanalyse erprobt, bei der integrale Systemcodes und sog. vereinfachte SWR-Modelle (ROM) als sich gegenseitig ergänzende Methoden eingesetzt werden, um die Stabilitätseigenschaften von Fixpunkten und periodischen Lösungen (Grenzzyklen) des nichtlinearen Differentialgleichungssystems, welches das Stabilitätsverhalten des SWR beschreibt, zu bestimmen. Das ROM, in denen das dynamische System durch gewöhnliche Differentialgleichungen (GDG) beschrieben wird, kann relativ einfach mit leistungsfähigen Methoden aus der nichtlinearen Dynamik, wie zum Beispiel die semianalytische Bifurkationsanalyse, gekoppelt werden. Mit solchen Verfahren kann, ohne das DG-System explizit lösen zu müssen, ein Überblick über mögliche Typen von stabilen und instabilen oszillatorischen Verhalten des SWR erhalten werden. Insbesondere sind die Stabilitätseigenschaften von Grenzzyklen, die in Hopf-Bifurkationspunkten entstehen, und die Bedingungen, unter denen sie auftreten, von Interesse. Mit dem Systemcode (RAMONA5) werden dann die mit dem ROM vorhergesagten Phänomene in den entsprechenden Parameterbereichen detaillierter untersucht (Validierung des ROM). Die Methodik dient daher nicht der Verfeinerung der Berechnung linearer Stabilitätsindikatoren (wie das DR). Das ROM-Gleichungssystem entsteht aus den PDGs des Systemcodes durch geeignete (nichttriviale) räumliche Mittelung der PDG. Es wird davon ausgegangen, dass die Reduzierung der räumlichen Komplexität die Stabilitätseigenschaften des SWR nicht signifikant verfälschen, da durch geeignete Mittlungsverfahren, räumliche Effekte näherungsweise in den GDGs berücksichtig werden. Beispielsweise wird die raum- und zeitabhängige Neutronenflussdichte nach räumlichen Moden entwickelt, wobei für eine Simulation der Stabilitätseigenschaften der In-phase- und Out-of-Phase-Leistungsoszillationen nur der Fundamentalmode und der erste azimuthale Mode berücksichtigt werden muss. Das ROM, welches ursprünglich am Paul Scherrer Institut (PSI, Schweiz) in Zusammenarbeit mit der Universität Illinois (USA) entwickelt wurde, ist in zwei wesentlichen Punkten erweitert und verbessert worden: • Entwicklung und Implementierung einer neuen Methode zur Berechnung der Rückkopplungsreaktivitäten • Entwicklung und Implementierung eines Modells zur Beschreibung der Rezirkulationsschleife (insbesondere wurde der Einfluss der Rezirkulationsschleife auf den In-Phase-Oszillationszustand und auf den Out-of-Phase-Oszillationszustand untersucht) • Entwicklung einer physikalisch begründeten Methode zur Berechnung der ROM-Inputdaten • Abschätzung des Einflusses des unterkühlten Siedens im Rahmen der ROM-Näherungen Mit dem erweiterten ROM wurden nichtlineare Stabilitätsanalysen für drei Arbeitspunkte (KKW Leibstadt (Zyklus 7) KKW Ringhals (Zyklus 14) und KKW Brunsbüttel (Zyklus 16)), für die Messdaten vorliegen, durchgeführt. In der Dissertationsschrift wird die RAM-ROM Methode ausführlich am Beispiel eines Arbeitspunktes (OP) des KKW Leibstadt (KKLc7_rec4-OP), in dem eine aufklingende regionale Leistungsoszillation bei einem Stabilitätstest gemessen worden ist, demonstriert. Das ROM sagt die Existenz eines Umkehrpunktes (saddle-node bifurcation of cycles, fold-bifurcation) voraus, der sich im linear stabilen Gebiet nahe der Stabilitätsgrenze befindet. Mit diesem ROM-Ergebnis ist eine neue Interpretation der Stabilitätseigenschaften des KKLc7_rec4-OP möglich. Die Resultate der in der Dissertation durchgeführten RAM-ROM Analyse bestätigen, dass das weiterentwickelte ROM für die Analyse des Stabilitätsverhaltens realer Leistungsreaktoren qualifiziert wurde.