This work focusses on analysing the stability of the MSFR – a molten salt reactor with a fast neutron spectrum. The investigations are based on a model, which was published and studied by the Politecnico di Milano using a linear approach. Since linear methods can only provide stability information to a limited extent, this work continues the conducted investigations by applying nonlinear methods. In order to examine the specified reactor model, the system equations were implemented, adjusted and verified using MATLAB code. With the help of the computational tool MatCont, a so-called fixed-point solution was tracked and its stability monitored during the variation of selected control parameters. It was found that the considered fixed point does not change its stability state and remains stable. Coexisting fixed points or periodic solutions could not be detected. Therefore, the analysed MSFR model is considered to be a stable system, in which the solutions always tend towards a steady state.:1. Introduction 2. Molten Salt Reactor Technology 2.1. Introduction 2.2. Historical Development 2.3. Working Principle of Molten Salt Reactors 2.4. Molten Salt Coolants 2.5. Advantages and Drawbacks 2.6. Classification 2.7. Molten Salt Fast Reactor Design 3. Stability Characteristics of Dynamical Systems 3.1. Introduction 3.2. Dynamical Systems 3.3. Stability Concepts 3.3.1. Introduction 3.3.2. Lagrange Stability (Bounded Stability) 3.3.3. Lyapunov Stability 3.3.4. Poincaré Stability (Orbital Stability) 3.4. Fixed-Point Solutions 3.4.1. Stability Analysis of Fixed-Point Solutions 3.4.2. Bifurcations of Fixed-Point Solutions 3.5. Periodic Solutions 3.5.1. Stability Analysis of Periodic Solutions 3.5.2. Bifurcations of Periodic Solutions 4. Analysed Reactor System 4.1. Introduction 4.2. Specified Reactor Model 4.3. Implementation and Verification of the Linearised System of Equations 4.3.1. Linearised System of Delayed Differential Equations 4.3.2. Comparison with Reference Plots 4.3.3. Adaptation of Parameter Values 4.4. Implementation and Verification of the Nonlinear System of Equations 4.4.1. Nonlinear System of Delayed Differential Equations 4.4.2. Delayed Neutron Precursor Equation Adjustments 4.4.3. Salt Temperature Equation Adjustments 4.4.4. Nonlinear System of Ordinary Differential Equations 4.4.5. Verification of the Nonlinear System of Ordinary Differential Equations 5. Conducted Stability Analyses 5.1. Introduction 5.2. Nonlinear Stability Analysis 5.2.1. Implementation 5.2.2. Results 5.2.3. Interpretation 5.3. Linear Stability Analysis 5.3.1. Comparison Between the Linearised and Nonlinearised MSFR System of Equations 5.3.2. Stability Investigations Using a Linear Criterion 5.4. MatCont Reliability Test Using an MSBR Model 6. Conclusions and Recommendations for Future Studies / Im Fokus dieser Arbeit steht die Stabilitätsanalyse des MSFR – eines Flüssigsalzreaktors mit schnellem Neutronenspektrum. Als Grundlage wurde ein Modell verwendet, das am Politecnico di Milano erstellt und dort mittels linearer Methoden untersucht wurde. Da lineare Betrachtungen nur eingeschränkte Stabilitätsaussagen treffen können, erweitert diese Arbeit die Untersuchungen um die nichtlineare Stabilitätsanalyse. Zur Untersuchung des vorgegebenen Reaktormodells wurden die Systemgleichungen in MATLAB übertragen und verifiziert. Mithilfe der Rechensoftware MatCont wurde eine sogenannten Fixpunkt-Lösung des Modells unter der Variation ausgewählter Parameter verfolgt und deren Stabilität überprüft. Es hat sich gezeigt, dass der betrachtete Fixpunkt seinen Stabilitätszustand dabei nicht verändert und stabil bleibt. Koexistierende Fixpunkte oder periodische Lösungen konnten nicht nachgewiesen werden. Daher gilt das betrachtete MSFR-Modell als ein stabiles System, dessen Lösungen immer auf einen stationären Zustand zulaufen.:1. Introduction 2. Molten Salt Reactor Technology 2.1. Introduction 2.2. Historical Development 2.3. Working Principle of Molten Salt Reactors 2.4. Molten Salt Coolants 2.5. Advantages and Drawbacks 2.6. Classification 2.7. Molten Salt Fast Reactor Design 3. Stability Characteristics of Dynamical Systems 3.1. Introduction 3.2. Dynamical Systems 3.3. Stability Concepts 3.3.1. Introduction 3.3.2. Lagrange Stability (Bounded Stability) 3.3.3. Lyapunov Stability 3.3.4. Poincaré Stability (Orbital Stability) 3.4. Fixed-Point Solutions 3.4.1. Stability Analysis of Fixed-Point Solutions 3.4.2. Bifurcations of Fixed-Point Solutions 3.5. Periodic Solutions 3.5.1. Stability Analysis of Periodic Solutions 3.5.2. Bifurcations of Periodic Solutions 4. Analysed Reactor System 4.1. Introduction 4.2. Specified Reactor Model 4.3. Implementation and Verification of the Linearised System of Equations 4.3.1. Linearised System of Delayed Differential Equations 4.3.2. Comparison with Reference Plots 4.3.3. Adaptation of Parameter Values 4.4. Implementation and Verification of the Nonlinear System of Equations 4.4.1. Nonlinear System of Delayed Differential Equations 4.4.2. Delayed Neutron Precursor Equation Adjustments 4.4.3. Salt Temperature Equation Adjustments 4.4.4. Nonlinear System of Ordinary Differential Equations 4.4.5. Verification of the Nonlinear System of Ordinary Differential Equations 5. Conducted Stability Analyses 5.1. Introduction 5.2. Nonlinear Stability Analysis 5.2.1. Implementation 5.2.2. Results 5.2.3. Interpretation 5.3. Linear Stability Analysis 5.3.1. Comparison Between the Linearised and Nonlinearised MSFR System of Equations 5.3.2. Stability Investigations Using a Linear Criterion 5.4. MatCont Reliability Test Using an MSBR Model 6. Conclusions and Recommendations for Future Studies

No description available.

Τα κοκκώδη υλικά είναι αναπόσπαστο κομμάτι του κόσμου μέσα στον οποίο ζει ο άνθρωπος, και συνεπώς, για την καλύτερη κατανόηση του κόσμου αυτού, επιβάλλεται η μελέτη τους. Αυτός είναι και ο σκοπός της παρούσας διδακτορικής διατριβής. Επικεντρωνόμαστε σε διάδρομο μεταφοράς ο οποίος αποτελεί αντιπροσωπευτικό μοντέλο για πληθώρα εφαρμογών τόσο στην βιομηχανία όσο και στο φυσικό περιβάλλον. Αποτελεί επίσης χαρακτηριστικό παράδειγμα της οικογένειας ανοικτών πολυσωματιδιακών συστημάτων, η οποία βρίσκεται στην καρδιά της σύγχρονης επιστήμης της Πολυπλοκότητας. Αρχικά εισάγουμε το μοντέλο ροής στο οποίο το κοκκώδες υλικό αντιμετωπίζεται ως ένα ειδικό ρευστό (συνεχές μέσο) με εσωτερική απώλεια ενέργειας. Μελετάμε τη δυναμική ισορροπία που επικρατεί στο σύστημα υπό σταθερές συνθήκες, καθώς και την κατάρρευση της ομαλής ροής μέσω του σχηματισμού συσσωματώματος. Ειδική μνεία γίνεται στα πρόδρομα φαινόμενα της συσσωμάτωσης, τα οποία ερμηνεύουμε μέσω μίας αντίστροφης διακλάδωσης διπλασιασμού περιόδου. Διερευνώντας την εξάρτηση μεταξύ της μορφής της ροϊκής συνάρτησης και του τρόπου με τον οποίο το σύστημα μεταβαίνει σε καθεστώς συσσωμάτωσης αποκαλύπτουμε τόσο ποιοτικές όσο και ποσοτικές διαφορές σε σχέση με τον παραπάνω τύπο διακλάδωσης. Μια σημαντική παραλλαγή του συστήματος μεταφοράς προκύπτει εφαρμόζοντας ανατροφοδότηση του πρώτου δοχείου με το συνολικό υλικό που εκρέει από το τελευταίο. Η μαθηματική επεξεργασία αποδεικνύει ότι σε αυτήν την περίπτωση η δημιουργία συσσωματώματος συντελείται μέσω μιας διακλάδωσης Hopf αντί για διακλάδωσης διπλασιασμού περιόδου. Επιστρέφοντας στο αρχικό μας σύστημα, μελετάμε και το συνεχές όριο, θεωρώντας το διάδρομο μεταφοράς να έχει «άπειρο» μήκος. Η δυναμική ισορροπία, που ισοδυναμεί με το ισοζύγιο της μάζας ανάμεσα σε διαδοχικά δοχεία του διακριτού συστήματος, τώρα παίρνει τη μορφή μιας μη γραμμικής μερικής διαφορικής εξίσωσης δεύτερης τάξης με μη σταθερούς συντελεστές. Η προσεκτική μελέτη της εξίσωσης και των συντελεστών της, σε συνδυασμό πάντα με τις συνοριακές συνθήκες στην είσοδο και έξοδο του διαδρόμου, μας επιτρέπει όχι μόνο να αναπαραγάγουμε τα προηγούμενα αποτελέσματα υπό το πρίσμα του συνεχούς ορίου αλλά και να τα ερμηνεύσουμε βάσει φυσικών διεργασιών όπως είναι η μεταφορά (drift) και η διάχυση (diffusion). Ειδικότερα, η συσσωμάτωση συμβαίνει σε καθεστώς αρνητικής διάχυσης (antidiffusion). Κλείνουμε την διατριβή προτείνοντας γενικεύσεις των συστημάτων που ερευνήσαμε. Επεκτείνουμε το διάδρομο μεταφοράς σε πλέγματα δύο διαστάσεων και μελετάμε άλλα μοντέλα που σχετίζονται με ροές διακριτών σωματιδίων όπως είναι η κυκλοφορία οχημάτων στους αυτοκινητοδρόμους. / Granular materials are ubiquitous in nature and in our daily lives, and understanding their behavior is therefore of crucial importance. The present thesis wants to contribute to this. We focus on a conveyor belt, which is not only a representative model for numerous applications both in industry and the natural environment, but also a prime example of an open multi-particle system prone to spontaneous pattern formation. This places our study right in the center of the modern science of complexity. Initially we introduce the flux model, in which the granular material is treated as a special fluid (a continuous medium) with internal energy losses. We examine the dynamic equilibrium that exists in the system under steady state conditions and also the breakdown of this equilibrium when the inflow rate exceeds a certain critical threshold value, resulting in the formation of a cluster and the obstruction of the conveyor belt. We focus especially on the pre-clustering phenomena and find that these can be described mathematically by a reverse period doubling bifurcation. Investigating the relation between the precise form of the flux function and the way in which the transition to the clustered state takes place, we reveal that the above scenario via a reverse period doubling bifurcation is not universal. Also other bifurcation types are possible. An important variation of our transport system is obtained by applying a feedback mechanism: All the particles that flow out from the last compartment are inserted into the first, making the system closed with respect to matter (mass conservation). The mathematical analysis proves that in this case the cluster formation occurs via a Hopf bifurcation instead of a period doubling. Returning to our original system, we study its continuum limit by considering a conveyor belt of ‘infinite’ length. The dynamics of the system is now described by a second-order nonlinear partial differential equation with non-constant coefficients. A careful analysis of this PDE and its coefficients, in combination with the special boundary conditions at the entrance and exit of the system, allows us not only to reproduce the results of the discrete system in the setting of differential equations but also to interpret these results in terms of physical processes such as drift and diffusion. In particular, the clustering occurs when the diffusion coefficient becomes negative, which gives antidiffusion. We close this thesis by discussing several generalizations of the system investigated. Among other things we expand the one-dimensional conveyor belt to a two-dimensional lattice. We further propose to use a similar flux model for the study of other, non-granular instances of discrete particle flows, such as vehicles on a highway.

Dans ce mémoire, nous étudions le problème centre-foyer sur un système polynomial. Nous développons ainsi deux mécanismes permettant de conclure qu’un point singulier monodromique dans ce système non-linéaire polynomial est un centre. Le premier mécanisme est la méthode de Darboux. Cette méthode utilise des courbes algébriques invariantes dans la construction d’une intégrale première. La deuxième méthode analyse la réversibilité algébrique ou analytique du système. Un système possédant une singularité monodromique et étant algébriquement ou analytiquement réversible à ce point sera nécessairement un centre. Comme application, dans le dernier chapitre, nous considérons le modèle de Gauss généralisé avec récolte de proies. / In this thesis, we study the center-focus problem in a polynomial system. We describe two mechanisms to conclude that a monodromic singular point in this polynomial system is a center. The first one is the method of Darboux. In this method, one uses invariant algebraic curves to build a first integral. The second method is the algebraic (and analytic) reversibility. A monodromic singularity, which is algebraically or analytically reversible at the singular point, is necessarily a center. As an application, in the last chapter, we consider the generalized Gause model with prey harvesting and a generalized Holling response function of type III.

Ce mémoire concerne la modélisation mathématique de l’érythropoïèse, à savoir le processus de production des érythrocytes (ou globules rouges) et sa régulation par l’érythropoïétine, une hormone de contrôle. Nous proposons une extension d’un modèle d’érythropoïèse tenant compte du vieillissement des cellules matures. D’abord, nous considérons un modèle structuré en maturité avec condition limite mouvante, dont la dynamique est capturée par des équations d’advection. Biologiquement, la condition limite mouvante signifie que la durée de vie maximale varie afin qu’il y ait toujours un flux constant de cellules éliminées. Par la suite, des hypothèses sur la biologie sont introduites pour simplifier ce modèle et le ramener à un système de trois équations différentielles à retard pour la population totale, la concentration d’hormones ainsi que la durée de vie maximale. Un système alternatif composé de deux équations avec deux retards constants est obtenu en supposant que la durée de vie maximale soit fixe. Enfin, un nouveau modèle est introduit, lequel comporte un taux de mortalité augmentant exponentiellement en fonction du niveau de maturité des érythrocytes. Une analyse de stabilité linéaire permet de détecter des bifurcations de Hopf simple et double émergeant des variations du gain dans la boucle de feedback et de paramètres associés à la fonction de survie. Des simulations numériques suggèrent aussi une perte de stabilité causée par des interactions entre deux modes linéaires et l’existence d’un tore de dimension deux dans l’espace de phase autour de la solution stationnaire. / This thesis addresses erythropoiesis mathematical modeling, which is the process of erythrocytes production and its regulation by erythropeitin. We propose an erythropoiesis model extension which includes aging of mature cells. First, we consider an age-structured model with moving boundary condition, whose dynamics are represented by advection equations. Biologically, the moving boundary condition means that the maximal lifespan varies to account for a constant degraded cells flux. Then, hypotheses are introduced to simplify and transform the model into a system of three delay differential equations for the total population, the hormone concentration and the maximal lifespan. An alternative model composed of two equations with two constant delays is obtained by supposing that the maximal lifespan is constant. Finally, a new model is introduced, which includes an exponential death rate depending on erythrocytes maturity level. A linear stability analysis allows to detect simple and double Hopf bifurcations emerging from variations of the gain in the feedback loop and from parameters associated to the survival function. Numerical simulations also suggest a loss of stability caused by interactions between two linear modes and the existence of a two dimensional torus in the phase space close to the stationary solution.

Cette thèse s'inscrit dans le cadre général de l'étude de la dynamique de populations. Elle porte sur la modélisation et l'analyse mathématique de l'hématopoïèse, le processus de production et de régulation des cellules sanguines. La population de cellules est perçue comme un milieu continu avec une structuration en âge et en espace. Nous avons commencé par analyser des modèles d'équations différentielles et aux différences à retard discret et distribué. Ces modèles à retard permettent de mettre en évidence des comportements particuliers tels que l'existence de solutions périodiques. Ensuite, nous avons pris en compte l'aspect spatial et la diffusion des cellules dans ces modèles, tout en sachant que la structuration en espace, dans le cas de l'hématopoïèse, a été très peu abordée par le passé. Un nouveau modèle a été obtenu du point de vue mathématique. Une étude d'existence d'ondes progressives est effectuée lorsque le domaine est non borné et lorsque le domaine est borné une étude de stabilité des états stationnaires ainsi que de l'existence d'une bifurcation de Hopf est réalisée / This thesis focuses on the study of population dynamics. It is devoted to the mathematical analysis and modeling of hematopoiesis, which is the process leading to the production and regulation of blood cells. The cell's population is seen as a continuous medium structured in age and space. We analyzed models of differential-difference system with discrete- and distributed -delay. These models can exhibit specific behaviors such as the existence of periodic solutions. Then we consider a space structuration and the diffusion of cells in such models, knowing that the space structure has not been widely studied in the case of hematopoiesis. A new model is obtained from the mathematical point of view. We studied the existence of traveling waves when the domain is unbounded. When the domain is bounded, the stability of stationary solutions and the existence of a Hopf bifurcation are obtained

No description available.

Dans ce mémoire, nous étudions le problème centre-foyer sur un système polynomial. Nous développons ainsi deux mécanismes permettant de conclure qu’un point singulier monodromique dans ce système non-linéaire polynomial est un centre. Le premier mécanisme est la méthode de Darboux. Cette méthode utilise des courbes algébriques invariantes dans la construction d’une intégrale première. La deuxième méthode analyse la réversibilité algébrique ou analytique du système. Un système possédant une singularité monodromique et étant algébriquement ou analytiquement réversible à ce point sera nécessairement un centre. Comme application, dans le dernier chapitre, nous considérons le modèle de Gauss généralisé avec récolte de proies. / In this thesis, we study the center-focus problem in a polynomial system. We describe two mechanisms to conclude that a monodromic singular point in this polynomial system is a center. The first one is the method of Darboux. In this method, one uses invariant algebraic curves to build a first integral. The second method is the algebraic (and analytic) reversibility. A monodromic singularity, which is algebraically or analytically reversible at the singular point, is necessarily a center. As an application, in the last chapter, we consider the generalized Gause model with prey harvesting and a generalized Holling response function of type III.

Dans cette thèse, le comportement d'un modèle mathématique permettant de transcrire la dynamique neuronale est étudié : le système de FitzHugh-Nagumo. En particulier, nous nous intéressons au caractère aléatoire d'ouverture et de fermeture des canaux ioniques d'un neurone qui reçoit ou non un stimulus. Ce caractère aléatoire de la dynamique neuronale est considéré, dans notre modèle, comme un bruit. Dans un premier temps, le comportement du modèle de FitzHugh-Nagumo a été caractérisé au voisinage de la bifurcation d'Andronov-Hopf qui traduit la transition entre l'état d'activation et l'état de repos du neurone. Classiquement, un neurone positionné à l'état de repos ne produit aucun potentiel d'action. Cependant, il a été montré un phénomène pour lequel une quantité appropriée de bruit permet la production de potentiels d'action des plus réguliers : la résonance cohérente. Le deuxième effet observé lors de simulations numériques permet au neurone d'améliorer la détection et l'encodage d'un signal subliminal : il s'agit de la résonance stochastique. De plus, cette thèse s'inscrit dans un contexte électronique puisqu'en plus de simuler numériquement le système de FitzHugh-Nagumo, les résultats de simulations ont également été confirmés en réalisant un circuit électronique. En effet, nous avons reproduit la dynamique non linéaire du système de FitzHugh-Nagumo à l'aide de ce circuit électronique. Cela a permis de mettre en évidence expérimentalement les deux phénomènes de résonance cohérente et de résonance stochastique pour lesquelles le bruit peut avoir une influence constructive sur le comportement de notre circuit électronique. / We study the nonlinear FitzHugh-Nagumo model witch describes the dynamics of excitable neural element. It is well known that this system exhibits three different possible responses. Indeed, the system can be mono-stable, oscillatory or bistable. In the oscillatory regime, the system periodically responds by generating action potential. By contrast, in the mono-stable state the system response remains constant after a transient. Under certain conditions, the system can undergo a bifurcation between the stable and the oscillatory regime via the so called Andronov-Hopf bifurcation. In this Phd thesis, we consider the FitzHugh-Nagumo model in the stable state, that is set near the Andronov-Hopf bifurcation. Moreover, we take into account the contribution of noise witch can induces two phenomena coherence resonance and stochastic resonance. First, without external driving, we show the effect of coherence resonance since a critical noise level enhances the regularity of the system response. Another numerical investigation reports how noise can allow to detect a subthreshold deterministic signal applied to the system. In this case, an appropriate amount of noise maximizes the signal to noise ratio reveling the stochastic resonance signature. Besides this numerical studies, we have also built a non linear circuit simulating the FitzHugh-Nagumo model under the presence of noise. This circuit has allowed to confirm experimentally the numerical observation of stochastic resonance and coherence resonance. Therefor, this electronic circuit contributes a framework for further experimental investigation in the field of neural sciences to better understand the role of noise in neural encoding.

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.