El objetivo de este trabajo es el desarrollo de técnicas para la simulación y optimización del diseño de circuitos osciladores y lazos de enganche en fase de microondas. La intención de estas técnicas es que puedan ser utilizadas por el diseñador para optimizar las prestaciones de este tipo de circuitos durante la etapa de diseño. Por este motivo, se ha intentado que en todo momento las técnicas puedan ser utilizadas en combinación con un programa comercial de simulación de circuitos de microondas.En el caso de los circuitos osciladores, inicialmente se han optimizado sus prestaciones cuando se utilizan como osciladores controlados por tensión. De esta forma, se han desarrollado una serie de técnicas que, en combinación con simulaciones en un programa comercial, permiten la linealización y extensión de la característica tensión-frecuencia. Mediante una técnica de control de estabilidad, se ha optimizado la respuesta dinámica del oscilador ante entradas variantes en el tiempo. En concreto, se ha aumentado la rapidez de respuesta eliminando transitorios lentos oscilantes que distorsionan la señal de salida deseada. Esta técnica se ha aplicado al caso particular de osciladores controlados por tensión utilizados para generar señales chirp, como puede ser en radares Frequency Modulated Continuos Wave (FMCW). Se ha analizado también el fenómeno del "injection-pulling", en el que una señal interferente desplaza la frecuencia de oscilación. Para ello, se ha desarrollado una formulación tipo transitorio de envolvente cuyos coeficientes pueden ser identificados mediante simulaciones de balance armónico en un simulador comercial. La técnica permite incrementar la robustez del circuito oscilador ante estas señales interferentes. Dados los problemas observados en el simulador comercial para simular la característica de ruido de fase en osciladores con estructuras acopladas, se ha desarrollado una técnica de simulación de ruido de fase que solventa estos problemas. La técnica obtiene la característica de ruido de fase a través de simulaciones de transitorio de envolvente en combinación con el uso de generadores auxiliares. Estas simulaciones pueden realizarse sin problemas usando un simulador comercial. Los resultados de todas las técnicas han sido corroborados mediante medidas en varios tipos de osciladores de microondas. Finalmente, se ha realizado un estudio preliminar para combinar el uso de series de Volterra con la técnica de transitorio de envolvente para la simulación de la respuesta transitoria de los osciladores.En el caso de los lazos de enganche en fase, se ha desarrollado un programa propio que realiza un análisis no lineal de lazos acoplados o "Coupled Phase-Locked Loops" (CPLL). Estos sistemas son utilizados en aplicaciones tales como en control de apuntamiento de antenas "phased-array". El programa, basado en una formulación tempo-frecuencial del sistema, permite la obtención de los rangos de operación del CPLL mediante una caracterización no lineal de los elementos que componen el lazo. Se delimitan los rangos de histéresis, y se analiza la variación de estos rangos en función de los parámetros del sistema. Se analiza la estabilidad de las soluciones estacionarias, teniendo en cuenta parámetros tales como el retardo del lazo. Mediante el control de la estabilidad y un análisis de tipo transitorio de envolvente, se optimiza la rapidez del sistema en el seguimiento de entradas moduladas. Finalmente, se analiza el ruido de fase, separando la perturbación en fase en distintas componentes. Esta separación permite clarificar el efecto del ruido en el control de apuntamiento de un array de antenas. Las predicciones de las técnicas han sido validadas mediante medidas en un sistema CPLL a 2 GHz. / The objective of this work is the development of techniques for the simulation and optimization of the design of microwave oscillator circuits and phase-locked loops. The intention of these techniques is that they can be used by the designer to optimize the features of these kinds of circuits during the design stage. For this reason, a lot of effort has been put along this thesis to ensure that the techniques can be used in combination with commercial microwave circuit simulators.In the case of the oscillator circuits, initially, their features have been optimized when used as voltage controlled oscillators (VCO). In this way, different techniques are proposed for the computer aided design of these circuits. The first technique allows setting the operation frequency band for specific values of the tuning voltage. The second technique imposes a linear frequency-voltage characteristic with the aid of an auxiliary generator. To follow this characteristic, the circuit is solved in terms of an ideal capacitance, synthesized, at a later stage, with the tuning varactor embedded in a linear network. In the third technique, the oscillator response to a sawtooth input, used to generate a chirp signal, is improved, eliminating spurious frequencies, not observable in steady state. To illustrate the techniques, a VCO operating in the C-band has been optimized and used to generate a chirp signal with low nonlinear frequency distortion. The injection-pulling phenomenon in oscillator circuits has been also analyzed. Injection pulling by interference signals is an undesired phenomenon in front-end oscillators, which causes a shift of the oscillation frequency and degrades the output spectrum. A semi-analytical formulation for the insightful analysis of injection-pulling phenomena in the presence of a modulated carrier or chirp signal is presented. The formulation enables an efficient analysis of interference problems difficult to simulate with harmonic balance or standard envelope transient. It allows the modification of the original design in order to reduce the injection pulling to the desired levels. The techniques have been applied to an oscillator at 6 GHz. Considering the problems found in commercial software to simulate the phase noise characteristic of coupled oscillator topologies, a numerical technique for the determination of the phase-noise spectrum of free-running oscillators is presented. The technique is based on envelope transient and can be applied to any commercial simulator on which this analysis method is available. The main advantage of the technique is that it allows simulating the near carrier phase noise spectrum, including possible resonances. The elements providing the oscillator phase-noise spectrum are obtained from envelope-transient simulations of low-computational cost. Comparisons are performed between the presented technique and other existing techniques, such as the carrier modulation approach. The technique has been successfully tested on the simulation of the near carrier phase noise spectrum of an oscillator circuit at 6.3 GHz. Finally, a preliminary study has been carried out to combine the use of Volterra series with the envelope transient technique for the simulation of oscillator transients.Regarding the phase-locked loops, in this thesis, harmonic-balance (HB) and envelope-transient formulations of coupled phase-locked loops (CPLLs) are presented. The CPLL has the added difficulty of its autonomous behavior since no reference oscillator is present. The new formulation takes into account the autonomy of the system, introducing a special set of state variables, which depend on the autonomous frequencies. The hysteresis phenomenon in CPLLs is analyzed in detail, efficiently obtaining the pull-in and hold-in ranges through HB. The pole analysis of the perturbed HB system enables an accurate prediction of instabilities and resonances. Due to the CPLL autonomy, there exists an inherent noise accumulation effect. This effect is taken into account, analyzing the perturbation in terms of accumulation and deviation components. The envelope formulation allows simulating the CPLL behavior in presence of modulation signals. The influence of the stability of the steady-state solution on the modulated signals is investigated. The simulation results have been successfully compared with the measurements in a manufactured CPLL system at 2 GHz.

Diese Arbeit stützt sich auf die theoretische und experimentelle Untersuchung der Thermodiffusion im Magnetfeld. Bei magnetischen Flüssigkeiten als kolloidalen Suspensionen versteht man unter der Thermodiffusion einen durch einen Temperaturgradienten angestoßenen unidirektionalen Partikeltransport, der zur Separation des Fluids führt. Beschrieben wird die Thermodiffusion theoretisch über das Konzentrationsprofil der Partikel in Abhängigkeit von Zeit und Ort in einer Fluidschicht. Die Experimente detektieren die Separation des Fluids über die Konzentrationsdifferenz zwischen zwei Fluidkammern. Die Bestimmung des Soret-Koeffizienten erfolgt über einen Datenfit zwischen experimentellen und theoretischen Daten. Für das kerosinbasierte Ferrofluid EMG905 wurden zwei Effekte festgestellt. Bei kleinen Magnetfeldstärken wandern die Partikel zum kalten Rand der Schicht (ST>0), bei steigenden Feldstärken kehrt sich diese Richtung um (ST<0). Die Ergebnisse der Untersuchungen zur Thermodiffusion gehen dann in eine lineare Stabilitätsanalyse einer Ferrofluidschicht bei anliegendem Temperaturgradienten und Magnetfeld ein. Dabei wird festgestellt, dass die kritische Rayleigh-Zahl als charakteristische Größe zum Einsetzen von Konvektion von dem Soret-Koeffizienten abhängt. Ist letzterer positiv, wird das Einsetzen von Konvektion begünstigt, ist er wiederum negativ, so kann Konvektion vollständig unterdrückt werden.

Dans ce projet, nous étudions la reconstruction de milieux terrestres souterrains.L’imagerie sismique est traitée avec un problème de minimisation itérative àgrande échelle, et nous utilisons la méthode de l’inversion des formes d’ondes(Full Waveform Inversion, FWI method). La reconstruction est basée sur desmesures d’ondes sismiques, car ces ondes sont caractérisées par le milieu danslequel elles se propagent. Tout d’abord, nous présentons les méthodesnumériques qui sont nécessaires pour prendre en compte l’hétérogénéité etl’anisotropie de la Terre. Ici, nous travaillons avec les solutions harmoniques deséquations des ondes, donc dans le domaine fréquentiel. Nous détaillons leséquations et l’approche numérique mises en place pour résoudre le problèmed’onde.Le problème inverse est établi afin de reconstruire les propriétés du milieu. Ils’agit d’un problème non-linéaire et mal posé, pour lequel nous disposons de peude données. Cependant, nous pouvons montrer une stabilité de type Lipschitzpour le problème inverse associé avec l’équation de Helmholtz, en considérantdes modèles représentés par des constantes par morceaux. Nous explicitons laborne inférieure et supérieure pour la constante de stabilité, qui nous permetd’obtenir une caractérisation de la stabilité en fonction de la fréquence et del’échelle. Nous revoyons ensuite le problème de minimisation associé à lareconstruction en sismique. La méthode de Newton apparaît comme naturelle,mais peut être difficilement accessible, dû au coup de calcul de la Hessienne.Nous présentons une comparaison des méthodes pour proposer un compromisentre temps de calcul et précision. Nous étudions la convergence de l’algorithme,en fonction de la géométrie du sous-sol, la fréquence et la paramétrisation. Celanous permet en particulier de quantifier la progression en fréquence, en estimantla taille du rayon de convergence de l’espace des solutions admissibles.A partir de l’étude de la stabilité et de la convergence, l’algorithme deminimisation itérative est conduit en faisant progresser la fréquence et l’échellesimultanément. Nous présentons des exemples en deux et trois dimensions, etillustrons l’incorporation d’atténuation et la considération de milieux anisotropes.Finalement, nous étudions le cas de reconstruction avec accès aux données deCauchy, motivé par les dual sensors développés en sismique. Cela nous permetde définir une nouvelle fonction coût, qui permet de prometteuses perspectivesavec un besoin minimal quant aux informations sur l’acquisition. / In this project, we investigate the recovery of subsurface Earth parameters. Weconsider the seismic imaging as a large scale iterative minimization problem, anddeploy the Full Waveform Inversion (FWI) method, for which several aspects mustbe treated. The reconstruction is based on the wave equations because thecharacteristics of the measurements indicate the nature of the medium in whichthe waves propagate. First, the natural heterogeneity and anisotropy of the Earthrequire numerical methods that are adapted and efficient to solve the wavepropagation problem. In this study, we have decided to work with the harmonicformulation, i.e., in the frequency domain. Therefore, we detail the mathematicalequations involved and the numerical discretization used to solve the waveequations in large scale situations.The inverse problem is then established in order to frame the seismic imaging. Itis a nonlinear and ill-posed inverse problem by nature, due to the limitedavailable data, and the complexity of the subsurface characterization. However,we obtain a conditional Lipschitz-type stability in the case of piecewise constantmodel representation. We derive the lower and upper bound for the underlyingstability constant, which allows us to quantify the stability with frequency andscale. It is of great use for the underlying optimization algorithm involved to solvethe seismic problem. We review the foundations of iterative optimizationtechniques and provide the different methods that we have used in this project.The Newton method, due to the numerical cost of inverting the Hessian, may notalways be accessible. We propose some comparisons to identify the benefits ofusing the Hessian, in order to study what would be an appropriate procedureregarding the accuracy and time. We study the convergence of the iterativeminimization method, depending on different aspects such as the geometry ofthe subsurface, the frequency, and the parametrization. In particular, we quantifythe frequency progression, from the point of view of optimization, by showinghow the size of the basin of attraction evolves with frequency. Following the convergence and stability analysis of the problem, the iterativeminimization algorithm is conducted via a multi-level scheme where frequencyand scale progress simultaneously. We perform a collection of experiments,including acoustic and elastic media, in two and three dimensions. Theperspectives of attenuation and anisotropic reconstructions are also introduced.Finally, we study the case of Cauchy data, motivated by the dual sensors devicesthat are developed in the geophysical industry. We derive a novel cost function,which arises from the stability analysis of the problem. It allows elegantperspectives where no prior information on the acquisition set is required.

L’approximation dite « anélastique » permet de filtrer les ondes acoustiques grâce à un développement asymptotique deséquations de Navier-Stokes, réduisant ainsi le pas en temps moyen, lors de la simulation numérique du développement d’instabilités hydrodynamiques. Ainsi, les équations anélastiques sont établies pour un mélange de deux fluides pour l’instabilité de Rayleigh-Taylor. La stabilité linéaire de l’écoulement est étudiée pour la première fois pour des fluides parfaits, par la méthode des modes normaux, dans le cadre de l’approximation anélastique. Le problème de Stokes issu des équations de Navier-Stokes sans les termes non linéaires (une partie de la poussée d’Archiméde est prise en compte) est défini ; l’éllipticité est démontrée, l’étude des modes propres et l’invariance liée à la pression sont détaillés. La méthode d’Uzawa est étendue à l’anélastique en mettant en évidence le découplage des vitesses en 3D, le cas particulier k = 0 et les modes parasites de pression. Le passage au multidomaine a permis d’établir les conditions de raccord (raccord Co de la pression sans condition aux limites physiques). Les algorithmes et l’implantation dans le code AMENOPHIS sont validés par les comparaisons de l’opérateur d’Uzawa développé en Fortran et à l’aide de Mathematica. De plus des résultats numériques ont été comparés à une expérience avec des fluides incompressibles. Finalement, une étude des solutions numériques obtenues avec les options anélastique et compressible a été menée. L’étude de l’influence de la stratification initiale des deux fluides sur le développement de l’instabilité de Rayleigh-Taylor est amorcée. / The « anelastic » approximation allows us to filter the acoustic waves thanks to an asymptotic development of the Navier-Stokes equations, so increasing the averaged time step, during the numerical simulation of hydrodynamic instabilitiesdevelopment. So, the anelastic equations for a two fluid mixture in case of Rayleigh-Taylor instability are established.The linear stability of Rayleigh-Taylor flow is studied, for the first time, for perfect fluids in the anelastic approximation.We define the Stokes problem resulting from Navier-Stokes equations without the non linear terms (a part of the buoyancyis considered) ; the ellipticity is demonstrated, the eigenmodes and the invariance related to the pressure are detailed.The Uzawa’s method is extended to the anelastic approximation and shows the decoupling speeds in 3D, the particular casek = 0 and the spurius modes of pressure. Passing to multidomain allowed to establish the transmission conditions.The algorithms and the implementation in the existing program are validated by comparing the Uzawa’s operator inFortran and Mathematica langages, to an experiment with incompressible fluids and results from anelastic and compressiblenumerical simulations. The study of the influence of the initial stratification of both fluids on the development of the Rayleigh-Taylor instability is initiated.

No description available.

Este trabalho tem a finalidade de estudar a influência da parede de concreto na análise de estabilidade de contenções atirantadas bem como discutir sobre segurança nestas análises. Para isto foram elaborados modelos em estado plano de deformação por meio do método dos elementos finitos, MEF, para análise. A parede de concreto foi modelada com variações de rigidez e modelos reológico, com o fim de se entender sua influência no fator de segurança. Por fim foi realizado um breve estudo sobre a utilização dos métodos estatísticos na análise de estabilidade de contenções. / This work has the purpose of study the influence of the concrete wall in the stability analysis of tieback retaining walls and to discuss these safety analysis. Models were developed using plane strain state via the finite element method, FEM, for analysis. The concrete wall was modeled with variations of stiffness and rheological models, in order to bore its influence on the safety factor. Finally a brief study was conducted on the use of statistical methods in stability analysis of retaining walls.

La présente thèse est consacrée à l’étude des peignes optiques de Kerr dans les résonateurs àmodes de galerie, au sein desquels la lumière peut être excitée par pompage externe. L’effet Kerrexistant dans ces résonateurs engendre des modes latéraux équidistants (dans le domaine spectral)de part et d’autre du mode excité, c’est à dire un peigne de fréquence. Cette thèse est diviséeen trois chapitres. Le premier est dédié à l’introduction de la génération de ces peignes et leurapplications. Le deuxième chapitre présente l’analyse de l’équation de Lugiato-Lefever, décrivantde manière analytique le système, et conduit à la construction de deux diagrammes de bifurcationpour les dispersions normale et anomale. Ils sont tracés en fonction des deux seuls paramètresexpérimentalement contrôlables une fois le résonateur fabriqué : la puissance du laser et sondécalage de fréquence. Ces diagrammes indiquent les plages de paramètres pour lesquels une,deux, ou trois solutions existent ainsi que leur stabilité. Les simulations numériques renseignentle type exact de solution associée à chaque aire (notamment les solitons brillants et sombres, lesbreathers, les peignes optiques de Kerr de premier et deuxième ordre, et un régime chaotique) ; cesdiagrammes indiquent donc les paramètres du laser à choisir afin de générer la solution souhaitée.Le troisième chapitre est dédié aux peignes de Kerr optique secondaires, lignes additionnelles dansle domaine spectral générées entre les lignes du peigne principal. Ils apparaissent en dispersionanormale, lorsque la quantité de photon pompe excède un seuil dit de second ordre, qui a étédéterminé numériquement. / This thesis is dedicated to the study of the Kerr optical frequency combs in whispering gallery moderesonators, where the light can be excited by the extern pump. Due to the Kerr effect existing in theseresonators, the quasi-equidistant lines in the spectral domain are generated around the excited mode,that is the frequency comb. This thesis is devided in three chapters. The first one is dedicated to theintroduction of the Kerr comb generation and their applications.The second one presents the analysisof the Lugiato-Lefever equation used for the analytical study of the system, leading to the constructionof two bifurcation diagrams for the normal and anomalous dispersions. They are plotted for twoparameters, which can be controlled during experiments once the resonator has been fabricated,which are the pump power of the laser and its frequency detuning. These diagrams show the areas ofthe parameters for which one, two, or three solutions exist and their stability. The additional numericalsimulations show the exact type of the solution in each area (such as the bright and dark solitons,the breathers, the primary and secondary Kerr combs and chaotical regimes), finally these diagramsshow the parameters of the laser needed to be choosen for the generation of the desired solution.The third chapter is dedicated to the secondary Kerr combs, which are the additional lines generatedbetween the lines of the primary comb. They appear in the anomalous dispersion regime, when thequantity of the pump photons crosses the second-order threshold, which has been found numerically.

Continued increase in system load leading to a reduction in operating margins, as well as the tendency to move towards a deregulated grid with renewable energy sources has increased the vulnerability of the grid to blackouts. Advanced intelligent techniques are therefore required to design new monitoring schemes that enable smart grid operation in a secure and robust manner. As the grid is highly interconnected, monitoring of transmission and distribution systems is increasingly relying on digital communication. Conventional security assessment techniques are slow, hampering real-time decision making. Hence, there is a need to develop fast and accurate security monitoring techniques. Intelligent techniques that are capable of processing large amounts of captured data are finding increasing scope as essential enablers for the smart grid. The research work presented in this thesis has evolved from the need for enhanced monitoring in transmission and distribution grids. The potential of intelligent techniques for enhanced system monitoring has been demonstrated for disturbed scenarios in an integrated power system. In transmission grids, one of the challenging problems is network partitioning, also known as network area-decomposition. In this thesis, an approach based on relative electrical distance (RED) has been devised to construct zonal dynamic equivalents such that the dynamic characteristics of the original system are retained in the equivalent system within the desired accuracy. Identification of coherent generators is another key aspect in power system dynamics. In this thesis, a support vector clustering-based coherency identification technique is proposed for large interconnected multi-machine power systems. The clustering technique is based on coherency measure which is formulated using the generator rotor measurements. These rotor measurements can be obtained with the help of Phasor Measurement Units (PMUs). In distribution grids, accurate and fast fault identification of faults is a key challenge. Hence, an automated fault diagnosis technique based on multi class support vector machines (SVMs) has been developed in this thesis. The proposed fault location scheme is capable of accurately identify the fault type, location of faulted line section and the fault impedance in the distributed generation (DG) systems. The proposed approach is based on the three phase voltage and current measurements available at all the sources i.e. substation and at the connection points of DGs. An approach for voltage instability monitoring in 3-phase distribution systems has also been proposed in this thesis. The conventional single phase L-index measure has been extended to a 3-phase system to incorporate information pertaining to unbalance in the distribution system. All the approaches proposed in this thesis have been validated using standard IEEE test systems and also on practical Indian systems.

The modern day power system is witnessing a tremendous change. There has been a rapid rise in the distributed generation, along with this the deregulation has resulted in a more complex system. The power demand is on a rise, the generation and trans-mission infrastructure hasn't yet adapted to this growing demand. The economic and operational constraints have forced the system to be operated close to its design limits, making the system vulnerable to disturbances and possible grid failure. This makes the study of voltage stability of the system important more than ever. Generally, voltage stability studies are carried on a single phase equivalent system assuming that the system is perfectly balanced. However, the three phase power system is not always in balanced state. There are a number of untransposed lines, single phase and double phase lines. This thesis deals with three phase voltage stability analysis, in particular the voltage stability index known as L-Index. The equivalent single phase analysis for voltage stability fails to work in case of any unbalance in the system or in presence of asymmetrical contingency. Moreover, as the system operators are giving importance to synchrophasor measurements, PMUs are being installed throughout the system. Hence, the three phase voltages can be obtained, making three phase analysis easier. To study the effect of unbalanced system on voltage stability a three phase L-Index based on traditional L-Index has been proposed. The proposed index takes into consideration the unbalance resulting due to untransposed transmission lines and unbalanced loads in the system. This index can handle any unbalance in the system and is much more realistic. To obtain bus voltages during unbalanced operation of the system a three phase decoupled Newton Raphson load ow was used. Reactive power distribution in a system can be altered using generators voltage set-ting, transformers OLTC settings and SVC settings. All these settings are usually in balanced mode i.e. all the phases have the same setting. Based on this reactive power optimization using LP technique on an equivalent single phase system is proposed. This method takes into account generator voltage settings, OLTC settings of transformers and SVC settings. The optimal settings so obtained are applied to corresponding three phase system. The effectiveness of the optimal settings during unbalanced scenario is studied. This method ensures better voltage pro les and decrease in power loss. Case studies of the proposed methods are carried on 12 bus and 24 bus EHV systems of southern Indian grid and a modified IEEE 30 bus system. Both balanced and unbalanced systems are studied and the results are compared.

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.