Les réseaux de distribution d’eau consistent en de grandes infrastructures réparties dans l’espace qui assurent la distribution d’eau potable en quantité et en qualité suffisantes. Les modèles mathématiques de ces systèmes sont caractérisés par un grand nombre de variables d’état et de paramètres dont la plupart sont incertains. Les temps de calcul peuvent s’avérer conséquents pour les réseaux de taille importante et la propagation d’incertitude par des méthodes de Monte Carlo. Par conséquent, les deux principaux objectifs de cette thèse sont l’étude des techniques de modélisation à ordre réduit par projection ainsi que la propagation spectrale des incertitudes des paramètres. La thèse donne tout d’abord un aperçu des méthodes mathématiques utilisées. Ensuite, les équations permanentes des réseaux hydrauliques sont présentées et une nouvelle méthode de calcul des sensibilités est dérivée sur la base de la méthode adjointe. Les objectifs spécifiques du développement de modèles d’ordre réduit sont l’application de méthodes basées sur la projection, le développement de stratégies d’échantillonnage adaptatives plus efficaces et l’utilisation de méthodes d’hyper-réduction pour l’évaluation rapide des termes résiduels non linéaires. Pour la propagation des incertitudes, des méthodes spectrales sont introduites dans le modèle hydraulique et un modèle hydraulique intrusif est formulé. Dans le but d’une analyse plus efficace des incertitudes des paramètres, la propagation spectrale est ensuite évaluée sur la base du modèle réduit. Les résultats montrent que les modèles d’ordre réduit basés sur des projections offrent un avantage considérable par rapport à l’effort de calcul. Bien que l’utilisation de l’échantillonnage adaptatif permette une utilisation plus efficace des états système pré-calculés, l’utilisation de méthodes d’hyper-réduction n’a pas permis d’améliorer la charge de calcul. La propagation des incertitudes des paramètres sur la base des méthodes spectrales est comparable aux simulations de Monte Carlo en termes de précision, tout en réduisant considérablement l’effort de calcul. / Water distribution systems are large, spatially distributed infrastructures that ensure the distribution of potable water of sufficient quantity and quality. Mathematical models of these systems are characterized by a large number of state variables and parameter. Two major challenges are given by the time constraints for the solution and the uncertain character of the model parameters. The main objectives of this thesis are thus the investigation of projection based reduced order modelling techniques for the time efficient solution of the hydraulic system as well as the spectral propagation of parameter uncertainties for the improved quantification of uncertainties. The thesis gives an overview of the mathematical methods that are being used. This is followed by the definition and discussion of the hydraulic network model, for which a new method for the derivation of the sensitivities is presented based on the adjoint method. The specific objectives for the development of reduced order models are the application of projection based methods, the development of more efficient adaptive sampling strategies and the use of hyper-reduction methods for the fast evaluation of non-linear residual terms. For the propagation of uncertainties spectral methods are introduced to the hydraulic model and an intrusive hydraulic model is formulated. With the objective of a more efficient analysis of the parameter uncertainties, the spectral propagation is then evaluated on the basis of the reduced model. The results show that projection based reduced order models give a considerable benefit with respect to the computational effort. While the use of adaptive sampling resulted in a more efficient use of pre-calculated system states, the use of hyper-reduction methods could not improve the computational burden and has to be explored further. The propagation of the parameter uncertainties on the basis of the spectral methods is shown to be comparable to Monte Carlo simulations in accuracy, while significantly reducing the computational effort.

Afin d'estimer les forces statiques équivalentes du vent, qui produisent les réponses quasi-statiques et dynamiques extrêmes dans les structures soumises au champ de pression instationnaire induit par les effets du vent, une nouvelle méthode probabiliste est proposée. Cette méthode permet de calculer les forces statiques équivalentes du vent pour les structures avec des écoulements aérodynamiques complexes telles que les toitures de stade, pour lesquelles le champ de pression n'est pas gaussien et pour lesquelles la réponse dynamique de la structure ne peut être simplement décrite en utilisant uniquement les premiers modes élastiques (mais nécessitent une bonne représentation des réponses quasi-statiques). Généralement, les mesures en soufflerie du champ de pression instationnaire appliqué à une structure dont la géométrie est complexe ne suffisent pas pour construire une estimation statistiquement convergée des valeurs extrêmes des réponses dynamiques de la structure. Une telle convergence est nécessaire pour l'estimation des forces statiques équivalentes afin de reproduire les réponses dynamiques extrêmes induites par les effets du vent en tenant compte de la non-gaussianité du champ de pression aléatoire instationnaire. Dans ce travail, (1) un générateur de réalisation du champ de pression instationnaire non gaussien est construit en utilisant les réalisations qui sont mesurées dans la soufflerie à couche limite turbulente; ce générateur basé sur une représentation en chaos polynomiaux permet de construire un grand nombre de réalisations indépendantes afin d'obtenir la convergence des statistiques des valeurs extrêmes des réponses dynamiques, (2) un modèle d'ordre réduit avec des termes d'accélération quasi-statique est construit et permet d'accélérer la convergence des réponses dynamiques de la structure en n'utilisant qu'un petit nombre de modes élastiques, (3) une nouvelle méthode probabiliste est proposée pour estimer les forces statiques équivalentes induites par les effets du vent sur des structures complexes décrites par des modèles éléments finis, en préservant le caractère non gaussien et sans introduire le concept d'enveloppes des réponses. L'approche proposée est validée expérimentalement avec une application relativement simple et elle est ensuite appliquée à une structure de toiture de stade pour laquelle des mesures expérimentales de pressions instationnaires ont été effectuées dans la soufflerie à couche limite turbulente / In order to estimate the equivalent static wind loads, which produce the extreme quasi-static and dynamical responses of structures submitted to random unsteady pressure field induced by the wind effects, a new probabilistic method is proposed. This method allows for computing the equivalent static wind loads for structures with complex aerodynamic flows such as stadium roofs, for which the pressure field is non-Gaussian, and for which the dynamical response of the structure cannot simply be described by using only the first elastic modes (but require a good representation of the quasi-static responses). Usually, the wind tunnel measurements of the unsteady pressure field applied to a structure with complex geometry are not sufficient for constructing a statistically converged estimation of the extreme values of the dynamical responses. Such a convergence is necessary for the estimation of the equivalent static loads in order to reproduce the extreme dynamical responses induced by the wind effects taking into account the non-Gaussianity of the random unsteady pressure field. In this work, (1) a generator of realizations of the non-Gaussian unsteady pressure field is constructed by using the realizations that are measured in the boundary layer wind tunnel; this generator based on a polynomial chaos representation allows for generating a large number of independent realizations in order to obtain the convergence of the extreme value statistics of the dynamical responses, (2) a reduced-order model with quasi-static acceleration terms is constructed, which allows for accelerating the convergence of the structural dynamical responses by using only a small number of elastic modes of the structure, (3) a novel probabilistic method is proposed for estimating the equivalent static wind loads induced by the wind effects on complex structures that are described by finite element models, preserving the non-Gaussian property and without introducing the concept of responses envelopes. The proposed approach is experimentally validated with a relatively simple application and is then applied to a stadium roof structure for which experimental measurements of unsteady pressures have been performed in boundary layer wind tunnel

Les conceptions récentes de turbomachines ont tendance à évoluer vers des liaisons entre étages de plus en plus souples et des niveaux d'amortissement faibles, donnant lieu à des configurations où les modes sont susceptibles de présenter des niveaux de couplages inter-étages forts. En général, les ensembles disques aubes multi-étagés n'ont aucune propriété de symétrie cyclique d'ensemble et l'analyse doit porter sur un modèle de la structure complète donnant lieu à des calculs très coûteux. Pour palier ce problème, une méthode récente appelée symétrie cyclique multi-étages peut être utilisée pour réduire le coût des calculs des rotors composés de plusieurs étages, même lorsque les étages ont un nombre différent de secteurs. Cette approche profite de la symétrie cyclique inhérente à chaque étage et utilise une hypothèse spécifique qui aboutit à des sous-problèmes découplés pour chaque ordre de Fourier spatial. La méthodologie proposée vise à étudier l'effet des incertitudes sur le comportement dynamique des rotors en utilisant l'approche de symétrie cyclique multi-étages et l'expansion en Chaos Polynomial. Les incertitudes peuvent découler de l'usure des aubes, des changements de température ou des tolérances de fabrication. En première approche, seules les incertitudes provenant de l'usure uniforme de l'ensemble des aubes sont étudiées. Celles-ci peuvent être modélisées en considérant une variation globale des propriétés du matériau de l'ensemble des aubes d'un étage particulier. L'approche de symétrie cyclique multi-étages peut alors être utilisée car l'hypothèse de secteurs identiques est respectée. La positivité des matrices aléatoires concernées est assurée par l'utilisation d'une loi gamma très adaptée à la physique du problème impliquant le choix des polynômes de Laguerre comme base pour le chaos polynomial. Dans un premier temps des exemples numériques représentatifs de différents types de turbomachines sont introduits dans le but d'évaluer la robustesse de la méthode de symétrie cyclique multi-étages. Ensuite, les résultats de l'analyse modale aléatoire et de la réponse aléatoire obtenus par le chaos polynomial sont validés par comparaison avec des simulations de Monte-Carlo. En plus des résultats classiquement rencontrés pour les fréquences et réponses forcées, les incertitudes considérées mettent en évidence des variations sur les déformées modales qui évoluent entre différentes familles de modes dans les zones de forte densité modale. Ces variations entraînent des modifications sensibles sur la dynamique globale de la structure analysée et doivent être considérées dans le cadre des conceptions robustes. / Recent designs in turbomachinery tend to have more flexible inter-stage rims and to be more lightly damped, resulting in configurations where modes might not be confined to only one stage. In general, multi-stage rotors have no particular axial symmetry property and the computationally costly analysis of the whole structure becomes mandatory. However, a multi-stage cyclic symmetry approach can be used for reducing the cost of modeling rotors composed of several stages even when the stages have different numbers of sectors. This approach takes advantage of the inherent cyclic symmetry of each stage and uses a specific assumption that results in decoupled subproblems for each spatial Fourier harmonic. The methodology proposed in this work allows including uncertainties in the analysis of multi-stage rotors using the multi-stage cyclic symmetry approach and the Polynomial Chaos Expansion. Uncertainties in rotors may arise from in-use wear of blades, temperature changes or manufacturing tolerances. As a first approach, only uncertainties arising from uniform in-use wear of the set of blades are included. These may be modeled by considering a global variation of the material properties of the set of blades of a particular stage. The multi-stage cyclic symmetry approach can then be used since the underlying assumption of identical sectors is respected. The positiveness of the random matrices involved is reached by using gamma-distributed random variables which imply the use of Laguerre's polynomials as basis for the polynomial chaos. Numerical examples representative of various types of turbomachinery are introduced in order to assess the robustness of the method of multi-stage cyclic symmetry. Uncertainties results for the free and forced response analyses obtained by the polynomial chaos are validated by comparison with Monte Carlo simulations. The considered uncertainties induce variations on the mode shapes that evolve between different families of modes in areas of high modal density. These variations result in significant changes in the global dynamics of the structure and must be considered in the context of robust designs.

De nos jours, la plupart des analyses CEM d’équipements et systèmes électroniques sont basées sur des approches quasi-déterministes dans lesquelles les paramètres internes et externes des modèles sont supposés parfaitement connus et où les incertitudes les affectant sont prises en compte sur les réponses par le biais de marges de sécurité importantes. Or, l’inconvénient de telles approches est qu’elles sont non seulement trop conservatives, mais en outre totalement inadaptées à certaines situations, notamment lorsque l’objectif de l’étude impose de prendre en compte le caractère aléatoire de ces paramètres via des modélisations stochastiques appropriées de type variables, processus ou champs aléatoires. Cette approche probabiliste a fait l’objet ces dernières années d’un certain nombre de recherches en CEM, tant au plan national qu’au plan international. Le travail présenté dans cette thèse est une contribution à ces recherches et a un double objectif : (1) développer et mettre en œuvre une méthodologie probabiliste et ses outils numériques d’accompagnement pour l’évaluation de la fiabilité et l’analyse sensibilité des équipements et systèmes électroniques en se limitant à des modélisations stochastiques par variables aléatoires ; (2) étendre cette étude au cas des modélisations stochastiques par processus et champs aléatoires dans le cadre d’une analyse prospective basée sur la résolution de l’équation aux dérivées partielles des télégraphistes à coefficients aléatoires.L’approche probabiliste mentionnée au point (1) consiste à évaluer la probabilité de défaillance d’un équipement ou d’un système électronique vis-à-vis d’un critère de défaillance donné et à déterminer l’importance relative de chacun des paramètres aléatoires en présence. Les différentes méthodes retenues à cette fin sont des adaptations à la CEM de méthodes développées dans le domaine de la mécanique aléatoire pour les études de propagation d’incertitudes. Pour le calcul des probabilités de défaillance, deux grandes catégories de méthodes sont proposées : celles basées sur une approximation de la fonction d’état-limite relative au critère de défaillance et les méthodes de Monte-Carlo basées sur la simulation numérique des variables aléatoires du modèle et l’estimation statistique des probabilités cibles. Pour l’analyse de sensibilité, une approche locale et une approche globale sont retenues. Ces différentes méthodes sont d’abord testées sur des applications académiques afin de mettre en lumière leur intérêt dans le domaine de la CEM. Elles sont ensuite appliquées à des problèmes de lignes de transmission et d’antennes plus représentatifs de la réalité.Dans l’analyse prospective, des méthodes de résolution avancées sont proposées, basées sur des techniques spectrales requérant les développements en chaos polynomiaux et de Karhunen-Loève des processus et champs aléatoires présents dans les modèles. Ces méthodes ont fait l’objet de tests numériques encourageant, mais qui ne sont pas présentés dans le rapport de thèse, faute de temps pour leur analyse complète. / Nowadays, most EMC analyzes of electronic or electrical devices are based on deterministic approaches for which the internal and external models’ parameters are supposed to be known and the uncertainties on models’ parameters are taken into account on the outputs by defining very large security margins. But, the disadvantage of such approaches is their conservative character and their limitation when dealing with the parameters’ uncertainties using appropriate stochastic modeling (via random variables, processes or fields) is required in agreement with the goal of the study. In the recent years, this probabilistic approach has been the subject of several researches in the EMC community. The work presented here is a contribution to these researches and has a dual purpose : (1) develop a probabilistic methodology and implement the associated numerical tools for the reliability and sensitivity analyzes of the electronic devices and systems, assuming stochastic modeling via random variables; (2) extend this study to stochastic modeling using random processes and random fields through a prospective analysis based on the resolution of the telegrapher equations (partial derivative equations) with random coefficients. The first mentioned probabilistic approach consists in computing the failure probability of an electronic device or system according to a given criteria and in determining the relative importance of each considered random parameter. The methods chosen for this purpose are adaptations to the EMC framework of methods developed in the structural mechanics community for uncertainty propagation studies. The failure probabilities computation is performed using two type of methods: the ones based on an approximation of the limit state function associated to the failure criteria, and the Monte Carlo methods based on the simulation of the model’s random variables and the statistical estimation of the target failure probabilities. In the case of the sensitivity analysis, a local approach and a global approach are retained. All these methods are firstly applied to academic EMC problems in order to illustrate their interest in the EMC field. Next, they are applied to transmission lines problems and antennas problems closer to reality. In the prospective analysis, more advanced resolution methods are proposed. They are based on spectral approaches requiring the polynomial chaos expansions and the Karhunen-Loève expansions of random processes and random fields considered in the models. Although the first numerical tests of these methods have been hopeful, they are not presented here because of lack of time for a complete analysis.

Le problème de la dimension stochastique élevée est récurrent dans les analyses probabilistes des structures. Il correspond à l’augmentation exponentielle du nombre d’évaluations du modèle mécanique lorsque le nombre de paramètres incertains est élevé. Afin de pallier cette difficulté, nous avons proposé dans cette thèse, une approche à deux étapes. La première consiste à déterminer la dimension stochastique efficace, en se basant sur une hiérarchisation des paramètres incertains en utilisant les méthodes de criblage. Une fois les paramètres prépondérants sur la variabilité de la réponse du modèle identifiés, ils sont modélisés par des variables aléatoires et le reste des paramètres est fixé à leurs valeurs moyennes respectives, dans le calcul stochastique proprement dit. Cette tâche fut la deuxième étape de l’approche proposée, dans laquelle la méthode de décomposition de la dimension est utilisée pour caractériser l’aléa de la réponse du modèle, par l’estimation des moments statistiques et la construction de la densité de probabilité. Cette approche permet d’économiser jusqu’à 90% du temps de calcul demandé par les méthodes de calcul stochastique classiques. Elle est ensuite utilisée dans l’évaluation de l’intégrité d’une toiture à ossature bois d’une habitation individuelle installée sur un site d’aléa sismique fort. Dans ce contexte, l’analyse du comportement de la structure est basée sur un modèle éléments finis, dans lequel les assemblages en bois sont modélisés par une loi anisotrope avec hystérésis et l’action sismique est représentée par huit accélérogrammes naturels fournis par le BRGM. Ces accélérogrammes permettent de représenter différents types de sols selon en se référant à la classification de l’Eurocode 8. La défaillance de la toiture est définie par l’atteinte de l’endommagement, enregistré dans les assemblages situés sur les éléments de contreventement et les éléments d’anti-flambement, d’un niveau critique fixé à l’aide des résultats des essais. Des analyses déterministes du modèle éléments finis ont montré que la toiture résiste à l’aléa sismique de la ville du Moule en Guadeloupe. Les analyses probabilistes ont montré que parmi les 134 variables aléatoires représentant l’aléa dans le comportement non linéaire des assemblages, 15 seulement contribuent effectivement à la variabilité de la réponse mécanique ce qui a permis de réduire la dimension stochastique dans le calcul des moments statistiques. En s’appuyant sur les estimations de la moyenne et de l’écart-type on a montré que la variabilité de l’endommagement dans les assemblages situés dans les éléments de contreventement est plus importante que celle de l’endommagement sur les assemblages situés sur les éléments d’anti-flambement. De plus, elle est plus significative pour les signaux les plus nocifs sur la structure. / The problem of the curse of dimensionality is frequently encountered in practical applications. It can be defined as the significant increase of the number of mechanical model calls with the number of uncertain parameters. To overcome this difficulty, a two-steps stochastic approach has been developed in this work. The first step of this approach consists in calculating the stochastic effective dimension by the means of Morris screening method. Once the most significant uncertain parameters on the variability of the mechanical responses are identified, they are modeled as random variables and the remaining parameters are fixed to their respective mean values. This allows us to reduce significantly the stochastic dimension of the problem in the second step of the approach where the decomposition method is used to estimate the statistical characteristics of the mechanical responses. The efficiency and the accuracy of this approach are evaluated through an academic problem dealing with the assessment of the integrity of a three-span five-story frame structure subjected to horizontal loads. We have demonstrate that we can reduce about 90% of the computation time required by the classical stochastic methods. Then, the proposed approach is used to the analysis of the integrity of timber roofs under seismic loading. The behaviour of this structure is described through a finite element model where the timber joints are modeled by anisotropic hysteresis law, and the seismic action is represented by eight real earthquake ground motion records. These accelerograms provided by the French institution involved in geosciences BRGM allow us to take into account different soil types according to the classification provided by the europeen design code dealing with seismic events Eurocode 8. The failure of timber roofs is reached when the damage levels in the timbers joints localized on the buckling and bracing members reach the critical value. It is shown, through a deterministic analysis, that the structure resists the seismic hazard representing the city of Le Moule in Guadeloupe. The stochastic analysis has shown that, among the 134 random variables representing the uncertainty in the nonlinear behaviour of the timber joints, only 15 have a significant effect on the variability of the structural response, which allow us to reduce the stochastic dimension in the computation of the statistical moments. According to the estimates of the mean and the standard deviation, we have shown that the variability of bracing members damage is greater than the variability of buckling members damage. Moreover, the variability of the bracing members damage is more significant for the earthquake ground motion records having the lowest collapse PGA.

Tesis por compendio / [ES] El flujo anular es uno de los regímenes de flujo bifásico más importantes y se caracteriza por que una fracción de líquido muy pequeña conocida como película de líquido que viaja cerca de la pared y un núcleo gaseoso. El flujo anular se puede observar durante la operación de plantas nucleares y en diferentes escenarios transitorios, aunque también en muchas otras aplicaciones industriales. La película de líquido es determinante en muchas de ellas ya que posee una alta capacidad de transferencia de masa, momento y energía. Parte de estas propiedades se deben a que la película presenta un comportamiento interfacial no linear con desarrollo de ondas interfaciales. Además, en determinadas instalaciones donde se la película de líquido actúa como refrigerante, es esencial conocer su comportamiento tanto por motivos de optimización como por razones de seguridad. Para estudiar los fundamentos del comportamiento de la película de líquido se han llevado a cabo una serie de experimentos en una instalación diseñada para generar flujo anular aire-agua en tubería circular vertical. En esta instalación se ha medido la evolución temporal del espesor de la película de líquido bajo diferentes condiciones y subrégimenes, como flujo en caída libre o flujo en cocorriente ascendente y descendente. El sistema de medida empleado se ha diseñado y construido para esta aplicación y consiste en sondas de conductancia de 3 electrodos rasantes a la pared y dispuestas en diferentes partes de la sección de test. Tanto el sistema electrónico como el dispositivo de calibración se diseñaron específicamente para trabajar con estas son-das de conductancia. La instalación cuenta con dos diámetros diferentes para poder comparar también el efecto del diámetro de la tubería así como aumentar el rango de medidas disponibles en bases de datos. Una de las características más particulares de la película de líquido son sus ondas interfaciales. Las principales ondas que se pueden diferenciar son las disturbance waves, ondas coherentes de gran calibre; y las ripple waves, ondas de pequeño tamaño, no coherentes que se generan constantemente antes de desaparecer al ser absorbidas por otras ondas. Las variables principales de la película de líquido que se han analizado en la instalación experimental son el espesor medio, la altura y frecuencia de las disturbance waves, la altura de las ripple waves y la altura de líquido no perturbado. Se han lleva-do a cabo diferentes estudios experimentales con objeto de añadir un valor adicional a las medidas. Para flujo anular descendente se ha estudiado el desarrollo de la película a través de diferentes zonas de medida y se han comparado las secciones de test de diferente diámetro. Además, múltiples correlaciones se han propuesto y los resultados se han comparado con estudios similares de otros autores. Para el análisis del flujo anular ascendente, se ha añadido un estudio del efecto de la tensión superficial en las variables de la película de líquido mediante la adición de pequeñas cantidades de 1-butanol. Es objeto de esta tesis también la modelación del flujo anular mediante análisis numérico. Los códigos de fluidodinámica computacional (CFD) son herramientas computacionales que permiten analizar el comportamiento de los fluidos. Han experimentado una fuerte evolución a lo largo de los últimos años gracias a los avances tecnológicos y los resultados que se obtienen de su correcta utilización son muy prometedores. No obstante, el flujo multifásico sigue siendo difícil de modelar y es necesario contrastar las predicciones de los códigos CFD con medidas experimentales. Por lo tanto, la fenomenología de flujo anular desarrollado se ha estudiado también mediante el código ANSYS CFX. Existe un importante vacío de conocimiento en la cuantificación de la incertidumbre (UQ) de dichos códigos CFD. Esta tesis también presenta los fundamentos del método UQ Polynomial Chaos Expansion (PCE) aplicado a dos casos prácticos. / [CA] El flux anular és un dels règims de flux bifàsic més importants i es caracteritza perquè la fracció de líquid és molt xicoteta coneguda com a pel·lícula de líquid. El flux anul·lar es pot observar durant l'operació de plantes nuclears i en diferents escenaris transitoris, encara que també en moltes altres aplicacions industrials. La pel·lícula de líquid és determinant en moltes d'elles ja que posseeix una alta capacitat de transferència de massa, moment i energia. Part d'aquestes propietats es deuen al fet que la pel·lícula presenta un comportament interfacial no linear amb desenvolupament d'ones interfacials. A més, en determinades instal·lacions on li la pel·lícula de líquid actua com a refrigerant, és essencial conéixer el seu comportament tant per motius d'optimització com per raons de seguretat. Per a estudiar els fonaments del comportament de la pel·lícula de líquid s'han dut a terme una sèrie d'experiments en una instal·lació dissenyada per a generar flux anular aïre-aigua en canonada circular vertical. En aquesta instal·lació s'ha mesurat l'evolució temporal de la grossària de la pel·lícula de líquid sota diferents condicions i subrégimenes, com a flux en caiguda lliure o flux en cocorriente ascendent i descendent. El sistema de mesura emprat s'ha dissenyat i construït per a aquesta aplicació i consisteix en sondes de conductància de 3 elèctrodes i disposades en diferents parts de la secció de test. Tant el sistema electrònic com el dispositiu de calibratge es van dissenyar específicament per a treballar amb aquestes sondes de conductància. La instal·lació compta amb dos diàmetres diferents per a poder comparar també l'efecte del diàmetre de la canonada així com augmentar el rang de mesures disponibles en bases de dades. Una de les característiques més particulars de la pel·lícula de líquid són les seues ones interfacials. Les principals ones que es poden diferenciar són les disturbance waves, ones coherents de gran calibre; i les ripple waves, ones de xicoteta grandària, no coherents que es generen constantment abans de desaparéixer en ser absorbides per altres ones. Les variables principals de la pel·lícula de líquid que s'han analitzat en la instal·lació experimental són la grossària mitjana, l'altura i freqüència de les disturbance waves, l'altura de les ripple waves i l'altura de líquid no pertorbat. S'han dut a terme diferents estudis experimentals a fi d'afegir un valor addicional a les mesures. Per a flux anul·lar descendent s'ha estudiat el desenvolupament de la pel·lícula a través de diferents zones de mesura i s'han comparat els diferents diàmetres. A més, múltiples correlacions s'han proposat i els resultats s'han comparat amb estudis similars d'altres autors. Per a l'anàlisi del flux anul·lar ascendent, s'ha afegit un estudi de l'efecte de la tensió superficial en les variables de la pel·lícula de líquid mitjançant l'addició de xicotetes quantitats de 1-butanol. És objecte d'aquesta tesi també el modelatge del flux anul·lar mitjançant anàlisi numèrica. Els codis de fluidodinámica computacional (CFD) són eines computacionals que permeten analitzar el comportament dels fluids. Han experimentat una forta evolució al llarg dels últims anys gràcies als avanços tecnològics i els resultats que s'obtenen de la seua correcta utilització són molt prometedors. No obstant això, el flux multifásico continua sent difícil de modelar i és necessari contrastar les prediccions dels codis CFD amb mesures experimentals. Per tant, la fenomenologia de flux anular desenvolupat s'ha estudiat també mitjançant el codi ANSYS CFX. Existeix un important buit de conei-xement en la quantificació de la incertesa d'aquests codis CFD. En aquesta tesi es mostren els fonaments del Polynomial Chaos Expansion (PCE) com a mètode per a calcular la incertesa dels resultats de simulació mitjançant propagació. El PCE per quadratura de Gauss-Hermite s'ha aplicat a les simulacions de dos experiments. / [EN] Annular flow is one of the most important two-phase flow regimes and is characterized by a very small liquid fraction known as a liquid film travelling close to the wall and a gas core. Annular flow can be observed during the operation of nuclear plants, in different transient scenarios, and many other industrial applications. The liquid film is decisive in many of them as it has a high mass, momentum and energy transfer capacity. Many of these properties are due to the film exhibiting nonlinear interfacial behavior with the generation of interfacial waves. In addition, in certain facilities where the liquid film acts as a coolant, it is essential to know its behavior both for optimization and safety reasons. In order to study the fundamentals of the liquid film, a series of experiments have been carried out in a facility designed to generate air-water annual flow in a vertical circular pipe. In this facility, the time evolution of the liquid film thickness has been measured under different conditions and sub-regimes, such as free-fall flow or upward and downward cocurrent flow. The measurement system used has been designed and built for this application and consists of 3-electrode conductance probes mounted flush to the wall and arranged at different distances from the entrance of the test section. Both the electronics and the calibration device were specifically designed to work with these conductance probes. The facility has two different diameters to compare the effect of the pipe diameter and increase the range of measurements available in databases. One of the main characteristics of the liquid film is its interfacial waves. The two primary types of waves that can be distinguished are the disturbance waves, which are large coherent waves, and the ripple waves, small, non-coherent waves that are constantly generated before disappearing when absorbed by other waves. The main variables of the liquid film analyzed in the experimental setup are the mean film thickness, the height and frequency of the disturbance waves, the height of the ripple waves and the height of the unperturbed liquid. Different experimental studies have been carried out to add additional value to the measurements. For downward annular flow, the development of the film through different measuring zones has been studied, and the different test section diameters have been compared. In addition, multiple correlations have been proposed, and the results have been compared with similar studies by other authors. To analyze the upward annular flow, a study of the effect of surface tension on the liquid film variables by adding small amounts of 1-butanol has been added. The modelling of annular flow by numerical analysis is also the subject of this thesis. Computational Fluid Dynamics (CFD) codes are computational tools that allow the analysis of fluid behavior. They have undergone a strong evolution over the last few years thanks to technological advances, and the results obtained from their correct use are very promising. However, multiphase flow remains challenging to model, and it is necessary to contrast the predictions of CFD codes with experimental measurements. Therefore, the developed annular flow phenomenology has also been studied using the ANSYS CFX code. There is a significant knowledge gap in the uncertainty quantification of CFD codes. Some methodologies are available, although many are in the early stages or have not been explored by researchers. All applications of CFD codes in nuclear safety require extensive knowledge of the uncertainty of the predictions, so developing these methodologies is crucial. This thesis shows the fundamentals of Polynomial Chaos Expansion (PCE) as a method to calculate the uncertainty of simulation results by propagation. The PCE by Gauss-Hermite quadrature has also been applied to the simulations of two experiments: the experimental setup of this thesis, and an international benchmark. / I would like to acknowledge the support provided by the Ministerio de Economía, Industria y Competitividad and the Agencia Nacional de Investigación under the FPI grant BES-2017-080031, which provided funding for my research. / Rivera Durán, Y. (2023). Experimental and Modelling Study of Interfacial Phenomena in Annular Flow with Uncertainty Quantification [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/194606 / Compendio