HTLV-I is the first discovered human retrovirus and a causative agent of both adult T-cell leukemia (ATL) and HTLV-I-associated myelopathy (or tropical spastic paraparesis) (HAM/TSP). Previous models have been successful in providing insight into the progression of HTLV-I infection. The relative simplicity of HTLV as well as its similarities to HIV and other diseases allow HTLV-I research to have diverse applications. The development of HAM/TSP is precipitated by a CTL immune response. Previous models for CTL response to HTLV-I infection have had relatively simple behaviours. A novel sigmoidal CTL response function results in complex behaviours previously unobserved. We establish the existence of bistability between solutions corresponding to carrier and endemic states. In addition, both super- and sub-critical Hopf bifurcations as well as the resulting stable and unstable periodic solutions are observed. Analytical and numerical results are discussed, as well as the biological consequences of the aforementioned behaviours. / Applied Mathematics

This thesis considers the problem of making dynamic models for industrial processes by combining physical modelling with experimental data. The focus is on distributed parameter systems, that is, systems for which the model structure involves partial differential equations (PDE). Distributed parameter systems are important in many applications, e.g., in chemical process systems and in intracellular biochemical processes, and involve for instance all forms of transport and transfer phenomena. For such systems, the postulated model structure usually requires a finite dimensional approximation to enable identification and validation using experimental data. The finite dimensional approximation involves translating the PDE model into a set of ordinary differential equations, and is termed model reduction. The objective of the thesis is two-fold. First, general PDE model reduction methods which are efficient in terms of model order for a given level of accuracy are studied. The focus here is on a class of methods called moving mesh methods, in which the discretization mesh is considered a dynamic degree of freedom that can be used for reducing the model reduction error. These methods are potentially highly efficient for model reduction of PDEs, but often suffer from stability and robustness problems. In this thesis it is shown that moving mesh methods can be cast as standard feedback control problems. Existing moving mesh methods are analyzed based on tools and results available from control theory, and plausible explanations to the robustness problems and parametric sensitivity experienced with these methods are provided. Possible remedies to these problems are also proposed. A novel moving finite element method, Orthogonal Collocation on Moving Finite Elements (OCMFE), is proposed based on a simple estimate of the model reduction error combined with a low order linear feedback controller. The method is demonstrated to be robust, and hence puts only small demands on the user. In the second part of the thesis, the integration of PDE model reduction methods with grey-box modelling tools available for finite dimensional models is considered. First, it is shown that the standard approach based on performing model reduction using some ad hoc discretization method and model order, prior to calibrating and validating the reduced model, has a number of potential pitfalls and can easily lead to falsely validated PDE models. To overcome these problems, a systematic approach based on separating model reduction errors from discrepancies between postulated model structures and measurement data is proposed. The proposed approach is successfully demonstrated on a challenging chromatography process, used for separation in biochemical production, for which it is shown that data collected at the boundaries of the process can be used to clearly distinguish between two model structures commonly used for this process. / QC 20101020

Electrical and optical characterizations of heterostructures and thin films based on group III-V compound semiconductors are presented. Optical properties of GaMnN thin films grown by Metalorganic Chemical Vapor Deposition (MOCVD) on GaN/Sapphire templates were investigated using IR reflection spectroscopy. Experimental reflection spectra were fitted using a non - linear fitting algorithm, and the high frequency dielectric constant (ε∞), optical phonon frequencies of E1(TO) and E1(LO), and their oscillator strengths (S) and broadening constants (Γ) were obtained for GaMnN thin films with different Mn fraction. The high frequency dielectric constant (ε∞) of InN thin films grown by the high pressure chemical vapor deposition (HPCVD) method was also investigated by IR reflection spectroscopy and the average was found to vary between 7.0 - 8.6. The mobility of free carriers in InN thin films was calculated using the damping constant of the plasma oscillator. The terahertz detection capability of n-type GaAs/AlGaAs Heterojunction Interfacial Workfunction Internal Photoemission (HEIWIP) structures was demonstrated. A threshold frequency of 3.2 THz (93 µm) with a peak responsivity of 6.5 A/W at 7.1 THz was obtained using a 0.7 µm thick 1E18 cm−3 n - type doped GaAs emitter layer and a 1 µm thick undoped Al(0.04)Ga(0.96)As barrier layer. Using n - type doped GaAs emitter layers, the possibility of obtaining small workfunctions (∆) required for terahertz detectors has been successfully demonstrated. In addition, the possibility of using GaN (GaMnN) and InN materials for terahertz detection was investigated and a possible GaN base terahertz detector design is presented. The non - linear behavior of the Inter Pulse Time Intervals (IPTI) of neuron - like electric pulses triggered externally in a GaAs/InGaAs Multi Quantum Well (MQW) structure at low temperature (~10 K) was investigated. It was found that a grouping behavior of IPTIs exists at slow triggering pulse rates. Furthermore, the calculated correlation dimension reveals that the dimensionality of the system is higher than the average dimension found in most of the natural systems. Finally, an investigation of terahertz radiation efect on biological system is reported.

By exploiting the concept of emergent network properties and the hierarchical nature of networks, we have constructed several levels of models facilitating the investigations of issues in the area of respiratory neural control. The first of such models is an intracellular second messenger pathway model, which has been shown to be an important contributor to intracellular calcium metabolism and mediate responses to neuromodulators such as serotonin. At the next level, we have constructed new single neuron models of respiratory-related neurons (e.g. the pre-Btzinger complex neuron and the Hypoglossal motoneuron), where the electrical activities of the neurons are linked to intracellular mechanisms responsible for chemical homeostasis. Beyond the level of individual neurons, we have constructed models of neuron populations where the effects of different component neurons, varying strengths and types of inter-neuron couplings, as well as network topology are investigated. Our results from these simulation studies at different structural levels are in line with experiment observations. The small-world topology, as observed in previous anatomical studies, has been shown here to support rhythm generation along with a variety of other network-level phenomena. The interactions between different inter-neuron coupling types simultaneously manifesting at time-scales orders of magnitude apart suggest possible explanations for variations in the outputs measured from the XII rootlet in experiments. In addition, we have demonstrated the significance of pacemakers, along with the importance of considering neuromodulations and second-messenger pathways in an attempt to understand important physiological functions such as breathing activities.

Materials when subjected to deformation exhibit unstable plastic flow beyond the elastic limit. In certain range of temperature and strain rates many solid state solutions, both interstitial as well as substitutional, exhibit the phenomenon of serrated yielding which also goes by the name, the Portevin - Le Chatelier (PLC) effect. The origin of this plastic instability is due to the interaction of dislocations with solute atoms. The objective of the thesis is to provide a dynamical systems approach to the study of this plastic flow instability. The thesis work discusses, within the framework of a model, the connection between microscopic dislocation mechanisms and macroscopic mechanical response of the specimen as stress drops in stress-strain curves. An extension of the model to the associated deformation bands is also considered. The emphasis is on the dynamical aspects of the instability. The methods of nonlinear dynamics like geometrical slow manifold and Poincare map formalism are applied for the first time to study the PLC effect. However, the approach and techniques transcend this particular application as the techniques are equally well applicable for many other physical systems as well, in particular, systems involving multiple time scales. The material covered should be of interest to investigators in the materials science, in particular, those, involved in the dislocation patterning and self organization of dislocations. Many theoretical models for the PLC effect exist in literature. Although the physical phenomenon is inherently dynamic, the conventional theoretical models do not involve any dynamical aspect. A dynamical model for this effect, due to Ananthakrishna, Sahoo and Valsakumar provides an explanation in terms of the dynamic interactions between different dislocation species and evolution of densities of these dislocation species. This model is known to reproduce several of the experimental results. It is within the perspective of this model and its extensions we analyze the PLC effect. The macroscopic manifestation of the PLC effect is the repeated load drops or serration in stress-strain curves (beyond the yield point). Each of the load drop is associated with the formation of a spatial dislocation band and its subsequent propagation. From the perspective of a dynamical system, the changeover from the stress-strain curve with single yield drop to repeated yield drops (the PLC effect) corresponds to a Hopf bifurcation wherein equilibrium state changes over to a periodic steady state. These repeated load drops correspond to auto oscillations of the applied stress (in the absence of any periodic driving force). In particular, as implied by the slow loading and sudden load drops, these oscillations are classified as relaxation oscillations. Relaxation oscillations are a result of disparate time scales of dynamics of the participating modes. Within the context of the model, this refers to very different time scales of evolution of densities of mobile (fast), immobile (slow) dislocations and those with a cloud of solute atoms (not too slow). The focus of attention in the thesis work is on these auto relaxation oscillations. There are several methodologies in nonlinear dynamical systems to study the oscillatory behavior of multidimensional systems with multiple time scales. An effective way is to study the reduced dynamical system in an appropriate space without sacrificing the required dynamical information. To this end, we discuss two techniques which compliment each other. 1.Slow manifold approach: This method utilizes the presence of multiple time scales dynamics. Advantage is that the information on the nature of evolution of the periodic orbit is retained. The limitation is that the transition from one stable state to another as parameter is varied cannot be dealt with. 2.Poincare maps:This approach utilizes the recurrent behavior of the period orbit. This is a convenient methodology to study the nature of stability of periodic orbits. However, in this, the information about the nature of evolution is lost. Both the above techniques provide good description in the presence of high dissipation or larger separation of time scales of the participating modes. For slow manifold analysis, this leads to exact slow manifold structure while in the case of Poincare maps, it leads to simpler, lower dimensional attractors. Specific issues that are dealt with using these approaches and others in this thesis are the following. To start with, we first provide a comprehensive overview of the dynamical behavior as envisaged by the model system in physically relevant two parameter space. The existence of relaxation oscillations bounded by back-to-back Hopf bifurcation is a good representation of the fact that the PLC effect manifests only in a window of strain rates. Within this boundary of Hopf bifurcations relaxation oscillations destabilize to give rise to new states of order, including the chaotic states. The changes in the nature of these oscillations with control parameters is projected through the bifurcation diagrams and analyzed using techniques like Floquet multipliers, Lyapunovs exponents etc. After the identification of the relevant parameter space for the monoperiodic relaxation oscillations, we focus our attention on the time scales involved in these relaxation oscillations and its connection to the time scales apparent in serrations of the stress-strain curve of the PLC effect. This characteristic feature of the PLC effect, the stick-slip nature of stress-strain curves, is believed to result from the negative strain rate dependence of the flow stress. The latter is assumed to arise from a competition of the relevant time scales involved in the phenomenon. However, in the previous works, the identification and the role of the time scales in the dynamical phenomenon is not clear. The motivation of this part of the work is to identify the time scales involved in the stress drops of the time series and their origin. Since the dynamics involves distinct time scales, in the long time limit, the evolution is controlled only by the slow modes. Hence, the adiabatic elimination or quasi-steady state approximation of the fast modes leads to an invariant manifold, the slow manifold which is useful for the analysis of time scales. The geometry of the slow manifold which is atypical with two connected pieces is shown to be at the root of the relaxation oscillations. The analysis of the slow manifold structure helps to understand the time scales of the dynamics operating in different regions of the slow manifold. The analysis also helps us to provide a proper dynamical interpretation for the negative branch of the strain rate sensitivity of the flow stress. The slow-fast dynamical nature manifests itself through multiperiodic oscillations also, in the form of mixed mode oscillations (MMOs), which are oscillations with both large amplitude excursions as well as small amplitude loops. In MMOs, the small amplitude oscillatory loops are confined to one part of the slow manifold (around the fixed point) and the large amplitude excursions arise as jumps from one piece of the slow manifold to the other. More generally, MMOs are a characteristic feature of a family of dynamical systems which also exhibit alternate periodic-chaotic sequences in bifurcation portraits. Usually, the origin of these features is explained in terms of either the approach to a homoclinic bifurcation duo to a saddle fixed point (Shilnikov scenario) or a saddle orbit (Gavrilov-Shilnikov scenario). However, the dynamical model exhibits features from both the above scenarios. The emphasis of this study is on explaining the origin of the incomplete approach to a global bifurcation in the dynamical model. Apart from attempting to understand the complex bifurcation sequences, an additional motivation for this study is the apparent lack of systematic investigation into the incomplete approach to global bifurcation exhibited by a variety of physical systems. The method of the analysis is general and applicable to the family of MMO systems. In the model, using the structure of the bifurcation sequences, and the equilibrium fixed point, a local analysis shows that the approach to homoclinicity is asymptotic at best, and is a result of the ‘softening' of eigenvalues of the saddle equilibrium point. This softening, in turn, is a consequence of back-to-back Hopf bifurcation which reflects the constraint of the physical phenomenon, namely, the occurrence of the multiple stress drops only in an interval of the strain rates. The characteristic features, namely, MMOs, alternate periodic-chaotic sequences, and incomplete approach to homoclinicity are related to each other and arise as a consequence of the atypical slow manifold structure. The slow manifold structure analysis assumes that the evolution of the system is constrained within the neighborhood of the slow manifold which also implies that the dynamical system involves high dissipation. Hence, the dimension of the effective dynamics in the long time limit is reduced. The analysis reveals information regarding the structure of the periodic orbit for a given set of parameter values but does not provide any information regarding the nature of stability of the periodic orbits. However, any insight into the mechanism of the instability of the periodic orbits in the model may lead to a better understanding of the underlying physical phenomenon. Poincare maps and equivalent discrete dynamical systems provide a convenient means to obtain such an insight on the nature of the periodic solutions of the dynamical system. This methodology compliments the invariant slow manifold analysis, since in Poincare maps, the nature of the stability information is preserved at the expense of the structure of the periodic orbit. However, these two methodologies are not exclusive to each other, since the slow manifold structure as well as Poincare maps may be constructed using a common factor, namely, extremal values of the fast variable of the dynamical system. The methodologies adopted for the analysis assumes large dissipation arising out of the multiple time scale behavior such that the next maximal amplitude (NMA) maps can be modeled by one dimensional discrete dynamical systems. The dynamical portrait of the model shows differing nature of dynamics and consequently Poincare maps with different geometrical shapes in the {m,c) plane. Within the framework of one dimensional maps, these shapes can be schematically reconstructed using minimal information regarding the principal periodic orbit embedded in higher dimension and its nature of stability. This suggests that one dimensional maps might be sufficient to represent the higher dimensional dynamical system. For most of the parameter space, the NMA maps of the dynamical model possess characteristic features of a locally smooth maximum and asymptotically long tail. These features have been observed in many other physical systems, both experimental and model systems. Hence, this analysis is focused on a broader issue of Poincare maps in a family of dynamical systems with multiple time scale dynamics and mixed mode oscillations. Here, the dynamical model has been used as a representative dynamical system for this family. The scope of the study is to understand the dynamical features of the MMO systems within the framework of one dimensional systems. Specifically, by using some general constraints on the one dimensional map, we first analyze the basic mechanism that is responsible for the reversal of periodic sequences of RLk type which corresponds to the dominant periodic states of the MMO systems. This in turn allows us to understand the period adding sequences as well. The analysis also helps to demonstrate that the width of the periodic states contained within the chaotic regions bounded by two successive periodic states of the form RLk is smaller than that for RLk .To this end, we first construct a model map which mimics the dominant bifurcation sequences of MMO systems. This map is utilized to verify the analytical results for the parameter width of the periodic windows. This analysis also throws light on the origin of the ordered structure of the isolas of RLk periodic orbits, in MMO systems, which was shown to be the result of a back-to-back Hopf bifurcation. The results indicate the ubiquity in the qualitative dynamical features of physical systems from widely differing origin, exhibiting alternate periodic-chaotic sequences. Although the model for the PLC effect is successful in describing the features of the phenomenon, a shortcoming of the dynamical model has been the absence of the spatial aspect. A dominant process in the PLC effect is the movement of dislocations (mainly through cross glide) which is essentially nonlocal. This feature has been incorporated into the dynamical model through a 'diffusive' term for the mobile dislocations. Preliminary results indicate that various types of band propagation, as seen in experiments, are recovered. It is known that the solute atmosphere aggregation occurs primarily during the waiting time of the mobile dislocations after its arrest. As another extension, the present model has been revised to incorporate these aging effects also. An outline of the thesis is as follows. Focus of this thesis work is on the dynamical aspects of the PLC effect. The phenomenology and few techniques in nonlinear dynamics are introduced in Chapters 1 and 2. Chapter 3 provides a comprehensive tour of dynamical behavior of the model in physically relevant two-parameter space. The rest of the work is presented in three parts (six chapters). In the first part of the thesis, the structure of the relaxation oscillations in the phase space is analyzed using the topology of the slow manifold. A connection between the slow manifold structure and the negative strain rate sensitivity of the flow stress is attempted using this analysis (Chapter 4). As a natural extension, the approach is utilized for the analysis of multiperiodic relaxation oscillations also. The emphasis is on the connection between the dynamical behavior of the model and incomplete approach to a global bifurcation (Chapter 5). In the second part of the thesis, the stability properties of periodic orbits are analyzed in detail using the Poincare map formalism, complimenting the study on the structure of periodic orbits using slow manifold. The structure and gross features of the Poincare map are reproduced utilizing only minimum information regarding the principal periodic orbit in the multidimensional space (Chapter 6). Within the framework of one dimensional systems, we analyze the mechanisms responsible for the structure of bifurcation portraits of MMO systems (Chapter 7). Third and the last part, of work focuses on modeling the spatial aspect of the PLC effect and refinement of the dynamical model (Chapters). The last chapter, Chapter9, is devoted for discussion of the results and scope for future work.

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.

Dans le cas des sols et plus généralement des milieux granulaires, qui sont des matériaux<br />non-associés, des ruptures diffuses existent pour des états de contrainte strictement inclus à<br />l'intérieur de la condition limite de plasticité. Nous proposons d'analyser ces ruptures comme<br />un phénomène de bifurcation : la rupture diffuse est un mode de déformation qui correspond<br />à une branche bifurquée avec perte d'unicité constitutive au point de bifurcation. Les points<br />de bifurcation sont détectés à l'aide du signe du travail du second ordre, soit la forme locale<br />du critère de stabilité de Hill. Les analyses présentées portent principalement sur des<br />simulations directes par la Méthode des Eléments Discrets. Pour des assemblages<br />granulaires numériques de différentes densités, un domaine de bifurcation est mis en<br />évidence à l'intérieur du critère de Mohr-Coulomb. Des conditions de sollicitation<br />conduisant à la rupture du matériau à partir d'un point de bifurcation sont données et des cas<br />de rupture diffuse sont modélisés. On parvient ainsi à reproduire et prévoir des ruptures non<br />décrites dans le cadre de l'élasto-plasticité classique. Finalement les origines microscopiques<br />(à l'échelle des contacts) du travail du second ordre sont analysées.

Ce travail est motivé par l'effort international actuel de créer expérimentalement une dynamo fluide auto-entretenue. L'effet dynamo, dont l'existence a été prevu par Larmor au début du XXème siècle, est considéré comme étant responsable de la production du champ matnétique terrestre et d'autres objets célestes par l'intermédiaire de l'écoulement d'un fluide conducteur. Afin d'étudier numériquement l'écoulement de von Kármán, qui modélise la configuration d'une expérience dynamo mise en place à Cadarache, nous avons développé une approche numérique originale permettant la résolution des équations magnétohydrodynamiques dans une géométrie dylindrique en formulation potentielle. La décomposition en potentiels poloïdal et toroïdal a été utilisée pour garantir la nature solénoïdale des champs de vitesse et magnétique. Nous utilisons la technique de la matrice d'influence pour satisfaire aux conditions aux limites et aux conditions de continuité du champ magnétique à la paroi du cylindre. La grande puissance de calcul, résultant de la parallélisation MPI du code, a presmis de l'appliquer deux problèmes concernant la turbulence dans la géométrie cylindrique : la turbulence axisymétrique et une bifurcation entre états turbulents.

No description available.

No description available.