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.

This thesis consists of two main parts. First of all, the condition number of rectangular Vandermonde matrices with nodes on the complex unit circle is studied. The first time quantitative bounds for the extreme singular values are proven in the multivariate setting and when nodes of the Vandermonde matrix form clusters. In the second part, an optimized presentation of the deterministic stability analysis of the subspace method ESPRIT is given and results from the first part are applied.

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.

Η παρούσα διατριβή επανεξετάζει την οριακή ανάλυση ευστάθειας των λίθινων αψίδων. Η οριακή ανάλυση ευστάθειας χρησιμοποιείται σήμερα ως το βασικό εργαλείο αποτίμησης της ευστάθειας τόξων και θολωτών κατασκευών από τοιχοποιία, όπως ακριβώς συνέβαινε και τους τελευταίους τέσσερις αιώνες. Παρά την τόσο μακρόχρονη ιστορία της μεθόδου, δεν έχουν πλήρως διασαφηνιστεί στην επιστημονική κοινότητα θεμελιώδης έννοιες και δεν έχουν σαφώς απαντηθεί ερωτήματα όπως: Ποιες είναι οι φυσικά πραγματοποιήσιμες γραμμές ώθησης και ποιες όχι; Ποια είναι η επίδραση της στερεοτομίας ενός τόξου στην οριακή του ευστάθεια; Ποιος είναι ο ρόλος της αλυσοειδούς καμπύλης και κατά πόσο αυτή είναι μία φυσικά αποδεκτή γραμμή ώθησης; Τι σχέση υπάρχει ανάμεσα στην κλίση της συνισταμένης θλιπτικής δύναμης και στην κλίση της γραμμής ώθησης στο σημείο εφαρμογής της; Η παρούσα διατριβή αναζητά απαντήσεις στα ερωτήματα αυτά, και έχει ως στόχο τη βαθύτερη κατανόηση της οριακής ανάλυσης ευστάθειας των τόξων, με παράλληλη ανάδειξη νέων υπολογιστικών διαδικασιών. Η δομή της παρουσιάζεται συνοπτικά παρακάτω. Στο πρώτο κεφάλαιο γίνεται ιστορική ανάλυση της μεθόδου μέσα από παρουσίαση και σχολιασμό των εργασιών με τη σημαντικότερη συμβολή, από τα μέσα του 17ου αιώνα μέχρι σήμερα. Στο δεύτερο κεφάλαιο επανεξετάζεται ένα από τα πιο κλασικά προβλήματα της μηχανικής: ποιο είναι το ελάχιστο επιτρεπτό πάχος ενός ημικυκλικού τόξου υπό τη δράση του ιδίου βάρους του για να είναι ευσταθές. Παράλληλα απαντώνται τα ερωτήματα που τέθηκαν παραπάνω αναπτύσσοντας νέες κλειστές μαθηματικές εκφράσεις των γραμμών ώθησης μέσω γεωμετρικής προσέγγισης, αλλά και μέσω του λογισμού των μεταβολών. Στο τρίτο κεφάλαιο χρησιμοποιείται παρόμοια διαδικασία για την ανάλυση της γενικής περίπτωσης των ελλειπτικών τόξων, οποιουδήποτε γεωμετρικού λόγου ύψος προς βάση, καθώς δεν είναι διαθέσιμα αναλυτικά αποτελέσματα στη διεθνή βιβλιογραφία, όπως συμβαίνει για τα κυκλικά τόξα. Στο τέταρτο κεφάλαιο εξετάζεται η οριακή ευστάθεια κυκλικών τόξων οποιασδήποτε γωνίας εναγκαλισμού, υπό την ταυτόχρονη δράση του ιδίου βάρους τους και σταθερής οριζόντιας εδαφικής επιτάχυνσης, ενώ υπολογίζεται με ακρίβεια η μορφή που θα έχει ο επικείμενος μηχανισμός κατάρρευσης μαζί με το οριακό πάχος, συναρτήσει της σεισμικής φόρτισης. Τα αποτελέσματα της μαθηματικής ανάλυσης (Κεφ. 2-4) επιβεβαιώνουν την ακρίβεια του λογισμικού που αναπτύχθηκε για τις ανάγκες της διατριβής, καθώς και τα αποτελέσματα που προκύπτουν από εμπορικό λογισμικό της μεθόδου των διακριτών στοιχείων. Στο πέμπτο κεφάλαιο γίνεται εφαρμογή και σύγκριση των πιο αντιπροσωπευτικών υπολογιστικών μεθόδων που απαντώνται σήμερα στη βιβλιογραφία για την αποτίμηση της ευστάθειας και φέρουσας ικανότητας της υπόγειας Θολωτής Διόδου του Σταδίου της Αρχαίας Νεμέας, ενώ η οριακή ανάλυση ευστάθειας αναδεικνύεται ως ένα μοναδικό εργαλείο για την κατανόηση της αλληλεπίδρασης της κατασκευής με το περιβάλλον έδαφος. Επιπλέων των συμπερασμάτων στο τέλος κάθε κεφαλαίου (Κεφ. 2 έως 5), στο έκτο κεφάλαιο παρουσιάζονται τα πιο σημαντικά συμπεράσματα και η συνεισφορά της παρούσας διατριβής. / This doctoral thesis revisits the limit equilibrium analysis of masonry arches. Limit equilibrium analysis is used today as the main analysis method for the assessment of the stability of masonry arches and vaulted structures, and is the outcome of important contributions that happened during the last four centuries. Although this method has a long history and a rich literature, there are still fundamental concepts that have not been thoroughly clarified, such as: What are the physically admissible thrust lines of an arch? How the stereotomy of an arch affects its limit stability? What is the role of the catenary curve (the alysoid)? Is the catenary curve a physically admissible thrust line? What is the relation between the direction of the thrust force and the slope of the thrust line at the point of application of the force? This thesis investigates these questions and aims to a better understanding of the limit equilibrium analysis of masonry arches, and at the same time, to present innovative methodologies and new analysis tools. Chapter 1 presents the work of other authors that have contributed the most to the stability analysis of masonry arches and vaulted structures over the last centuries. Chapter 2 revisits one of the most classical problems of Mechanics—what is the minimum thickness of a semicircular masonry arch subjected to its own weight. At the same time, the analysis presented in this chapter answers to the aforementioned questions through the development of closed-form expressions of the thrust line and the application of calculus of variation. Chapter 3 is focused on the limit equilibrium state of elliptical masonry arches, using the same approaches that were used in Chapter 2. This analysis was motivated from the fact that numerical results have been available in literature only for circular and not for elliptical masonry arches. Chapter 4 computes the location of the imminent hinges and the minimum thickness of circular masonry arches, for every given embrace angle, which can just sustain their own weight, together with a given level of horizontal ground acceleration. The numerical results presented in Chapters 2 to 4 confirm the accuracy of the in-house software that was developed for the needs of this thesis and the results obtained with a representative, commercially available software of the distinct element method. Chapter 5 present a comprehensive structural analysis of the Tunnel-Entrance to the Stadium of Ancient Nemea which ranges from the thrust line limit analysis and the discrete element method, to a 3-dimensional finite-element analysis. Limit equilibrium analysis emerges as a unique analysis method for the assessment of the stability of the structure and its interaction with the surrounding soil. While at the end of every chapter (Chapters 2 to 5) are presented detailed comments and conclusions, Chapter 6 is focused on outlining the most important conclusions and the main contribution of this thesis.