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.

Les différentes structures composant les moteurs aéronautiques requièrent des analyses dynamiques afin de prédire leur durée de vie. Pour des raisons d'allègement, les roues aubagées fixes de turbomachines, appelées redresseurs, sont conçus comme des ensembles de secteurs comportant plusieurs aubes. Cette architecture rompt la symétrie cyclique empêchant l'application des méthodes numériques l'exploitant. De plus, les dispersions géométriques et matériaux génèrent un désaccordage involontaire impliquant des zones de forte densité modale, dans lesquelles est observée une amplification de la réponse vibratoire, accrue par le caractère monobloc, et donc peu amorti, des secteurs. Une méthodologie statistique de prédiction du niveau vibratoire d'un secteur de redresseur désaccordé aléatoirement est développée ici. La modélisation des incertitudes est basée sur une approche paramétrique de la théorie probabiliste : des paramètres matériaux aléatoires suivant une loi uniforme sont associés à différentes parties du secteur. Une expansion de Karhunen-Loève permet de réduire le champ stochastique à un petit nombre de variables aléatoires et donc de diminuer les temps de calcul. Les modes stochastiques sont ensuite projetés sur ces espaces aléatoires par le biais de deux méthodes d'interpolation non-intrusives. La première est basée sur une projection sur une base du chaos polynomial tandis que la deuxième est une méthode de régression non-paramétrique (méthode MARS). Afin d'appliquer les deux méthodes de calcul à un modèle industriel, une méthode de double synthèse modale est appliquée permettant de diviser le temps de calcul des modes par un facteur d'environ 300. La sous-structuration adoptée s'adapte à la méthode de modélisation des incertitudes et s'avère robuste vis-à-vis du désaccordage. De plus, les deux méthodes permettent d'obtenir des résultats prédictifs en termes de moments statistiques tout en réduisant les temps de calculs. Enfin, la méthodologie est validée expérimentalement puisque l'enveloppe vibratoire numérique encadre la réponse fréquentielle expérimentale au niveau de la zone des modes d'intérêt. Une stratégie de positionnement des jauges de déformation est proposée à partir d'une distribution statistique des déplacements maximaux à mi-hauteur de veine sur une plage fréquentielle donnée. / Aircraft engine components necessitate extensive dynamical analyses in order to obtain life cycle prediction. In order to lighten the structure, turbomachinery stator bladed disks, called stator vanes, are designed as a set of multiple blades clusters. This architecture implies a loss of cyclic symmetry condition and prevents the use of numerical methods using it. Moreover, geometric dispersions and materials defaults generate an involuntary mistuning involving high modal density areas, in which is observed an amplification of the vibratory response, enhanced by the monobloc character - and hence low damped - of stator vanes. A statistical methodology for predicting the vibratory level of a randomly mistuned industrial stator vanes is developed here. Uncertainties modelization is based on a parametric approach of the probability theory : material random parameters following a uniform distribution are associated with different cluster's parts. A Karhunen-Loeve expansion reduces the stochastic field to a small number of random variables and therefore reduces the computation time. Stochastic modes are then projected on these random spaces through two non-intrusive methods of interpolation. The first is based on a projection on a polynomial chaos basis while the second is non-parametric regression method (MARS method). In order to implement both numerical methods to an industrial model, a double modal synthesis method is applied to divide the calculation time of modes by a factor around 300. The sub-structuring way adopted fits the uncertainties modelization method and is robust towards mistuning. Moreover, both methods yield predictive results in terms of statistical moments while reducing computation time. Finally, the methodology is experimentally validated because the numerical vibratory envelope frames the experimental frequency response at the area of the modes of interest. A positioning strategy of strain gauges is proposed based on a statistical distribution of the maximum displacements in vein halfway over a given frequency range.

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