Ce travail entre dans le cadre d'un programme de recherche consistant à développer des outils d'optimisation de la prise de décision des régatiers face aux situations micro-météorologiques complexes rencontrées lors des compétitions olympiques ou de match-racing. Notre objectif est l'amélioration des connaissances de l'évolution du vent sur les plans d'eau de régates côtières. Un dispositif d'étude de la micro-météorologie en baie de Quiberon incluant un volet expérimental et un volet numérique, a été mis en place. Plusieurs situations météorologiques ont été rencontrées pendant la campagne de mesures Quiberon2006 (16-29 Juin 2006). L'analyse des données collectées montre, à partir d'une classification en quadrants, l'influence du vent synoptique sur les transitions et les développements des phénomènes de brise. Des simulations à haute résolution spatiale de certaines journées sélectionnées ont été menées avec le modèle atmosphérique ARPS pour appréhender les effets de sites liés au découpage de la côte et à la présence de la presqu'île. La représentation fine de la topographie et des modes d'occupation des sols s'appuie sur l'analyse de la base de données CORINE Land Cover. Les champs météorologiques d'échelle régionale, servant de forçage au modèle à haute-résolution, proviennent d'assimilations effectuées avec le modèle WRF. La capacité des modèles à reproduire les phénomènes micro-météorologiques côtiers en site complexe est évaluée par comparaison aux données expérimentales. L'approche de modélisation à haute résolution permet de caractériser les comportements du vent à l'échelle du plan d'eau, caractéristiques difficilement quantifiables expérimentalement.

Les matériaux composites à matrice céramique CMC, sont généralement formés d'au moins deux matériaux ayant une forte capacité d'adhésion. Ces matériaux sont principalement composés de renforts fibreux assurant la tenue mécanique de la structure et d'une matrice qui permet sa cohésion. Utilisées principalement dans le domaine de l'aéronautique, elles sont reconnues pour leur bonne tenue mécanique, leur réfractarité élevée tout en conservant une densité faible. Par contre, l'inconvénient majeur associé à ces matériaux est l'apparition de fissures qui sont dues soit au procédé de fabrication soit aux sollicitations mécaniques externes.Dans ce travail, une attention toute particulière est consacrée aux composites à matrice auto-cicatrisante dont la principale propriété est l'aptitude à "réparer" les effets de la fissure par formation d'un verre visqueux.Ces verres visqueux se forment au sein de la fissure grâce à l'oxydation des éléments qui constituent la matrice. Selon la température, différents verres peuvent être formés.Leur fonction est de reboucher les fissures de taille micrométrique de façon à ralentir la diffusion de l'oxygène en direction des fibres et éviter leur rupture par oxydation.Cependant, pour des systèmes en rotation rapide tels que les turbines basse pression des moteurs d'avion (pièce étudiée actuellement), on peut s'interroger quant à la mobilisation du verre visqueux cicatrisant dans un système complexe géométriquement et inhomogène du point de vue de la nature des surfaces. Pour approcher le comportement du verre cicatrisant dans un système modèle mais néanmoins réaliste, une approche de modélisation numérique a été entreprise. L'outil numérique utilisé pour cette étude est le code de calcul Thétis développé à l'I2M. Celui-ci est adapté à ce type de simulation puisqu'il permet la modélisation d'écoulements diphasiques incluant des phénomènes physiques complexes tels que le mouillage. Ainsi, l'objectif de ce travail est-il de déterminer les limites d'utilisation de ce type de matériau en fonction des conditions auxquels il est soumis en évaluant la mobilité du verre cicatrisant dans la fissure.

En France, un tiers des personnes tuées sur la route le sont lors d'un accident sur un obstaclefixe. Dans 90% des cas, ces accidents surviennent après une perte de contrôle du véhicule.Les dispositifs de retenue de véhicule ont pour but de maintenir les véhicules en perdition surla chaussée en limitant la sévérité de l'impact.Ces dispositifs doivent subir des essais de chocs normatifs afin de pouvoir être installés sur lebord des routes européennes et d'évaluer leurs performances en termes de sévérité et dedéflexion.Les tolérances existantes sur les paramètres d'essai (véhicule, masse du véhicule, vitesse,angle et point d'impact ...) et les incertitudes sur les caractéristiques mécaniques desmatériaux constituant le dispositif ont un effet sur les performances de ce dispositifs etdoivent être prises en compte lors des calculsLes dispositifs mixtes (acier-bois) présentent une difficulté supplémentaire en raison del'hétérogénéité du matériau et de sa sensibilité aux variables d'environnement telles que latempérature et l'humidité.Afin de prendre en compte cette variabilité et d'évaluer son impact sur les performances d'undispositif, des essais dynamiques sur des échantillons de structure ont été réalisés et modélisésnumériquement.Enfin, un modèle complet d'un dispositif de retenue de véhicule a été effectué et corrélé surun essai de choc réel à l'aide d'une méthode prenant en compte la variation de paramètresphysiques liés à l'apparition des modes de ruine de la structure. Une fois corrélé, le modèle aété utilisé afin d'évaluer l'incidence de la modification des caractéristiques mécaniques dubois liée aux variations des conditions environnementales.

Une fistule artérioveineuse (FAV) est un accès vasculaire permanent créé par voie chirurgicale en connectant une veine et une artère chez le patient en hémodialyse. Cet accès vasculaire permet de mettre en place une circulation extracorporelle partielle afin de remplacer les fonctions exocrines des reins. En France, environ 36000 patients sont atteint d'insuffisance rénale chronique en phase terminale, stade de la maladie le plus grave qui nécessite la mise en place d'un traitement de suppléance des reins : l'hémodialyse. La création et présence de la FAV modifient significativement l'hémodynamique dans les vaisseaux sanguins, au niveau local et systémique ainsi qu'à court et à plus long terme. Ces modifications de l'hémodynamiques peuvent induire différents pathologies vasculaires, comme la formation d'anévrysmes et de sténoses. L'objectif de cette étude est de mieux comprendre le comportement mécanique et l'hémodynamique dans les vaisseaux de la FAV. Nous avons étudié numériquement les interactions fluide-structure (IFS) au sein d'une FAV patient-spécifique, dont la géométrie a été reconstruite à partir d'images médicales acquises lors d'un précédent doctorat. Cette FAV a été créée chez le patient en connectant la veine céphalique du patient à l'artère radiale et présente une sténose artérielle réduisant de 80% la lumière du vaisseau. Nous avons imposé le profil de vitesse mesuré sur le patient comme conditions aux limites en entrée et un modèle de Windkessel au niveau des sorties artérielle et veineuse. Nous avons considéré des propriétés mécaniques différentes pour l'artère et la veine et pris en compte le comportement non-Newtonien du sang. Les simulations IFS permettent de calculer l'évolution temporelle des contraintes hémodynamiques et des contraintes internes à la paroi des vaisseaux. Nous nous sommes demandées aussi si des simulations non couplées des équations fluides et solides permettaient d'obtenir des résultats suffisamment précis tout en réduisant significativement le temps de calcul, afin d'envisager son utilisation par les chirurgiens. Dans la deuxième partie de l'étude, nous nous sommes intéressés à l'effet de la présence d'une sténose artérielle sur l'hémodynamique et en particulier à ses traitements endovasculaires. Nous avons dans un premier temps simulé numériquement le traitement de la sténose par angioplastie. En clinique, les sténoses résiduelles après angioplastie sont considérées comme acceptables si elles obstruent moins de 30% de la lumière du vaisseau. Nous avons donc gonflé le ballonnet pour angioplastie avec différentes pressions de manière à obtenir des degrés de sténoses résiduelles compris entre 0 et 30%. Une autre possibilité pour traiter la sténose est de placer un stent après l'angioplastie. Nous avons donc dans un deuxième temps simulé ce traitement numériquement et résolu le problème d'IFS dans la fistule après la pose du stent. Dans ces simulations, la présence du stent a été prise en compte en imposant les propriétés mécaniques équivalentes du vaisseau après la pose du stent à une portion de l'artère. Dans la dernière partie de l'étude nous avons mis en place un dispositif de mesure par PIV (Particle Image Velocimetry). Un moule rigide et transparent de la géométrie a été obtenu par prototypage rapide. Les résultats expérimentaux ont été validés par comparaison avec les résultats des simulations numériques.

La réalisation par des procédés robotisés de pièces composites à destination, notamment, de l’aéronautique vise à en augmenter la productivité. Cependant le Placement de Fibres Robotisé (PFR) est encore en phase de maturation et requiert de nombreux développements, en particulier dans le cas des composites à matrice thermoplastique ou à fibre sèche. Dans cette thèse, nous proposons différents outils qui permettent in fine de déterminer à l’avance la meilleure puissance de chauffe pour la mise en œuvre de ces composites. La difficulté réside dans le fait que cette puissance dépend de nombreux paramètres, non seulement au niveau de la matière utilisée (densité, chaleur spécifique. . .), mais encore au niveau du procédé lui-même (vitesse de déplacement, nombre et orientation des plis. . .). Nous avons construit un modèle réduit multi-paramétrique en utilisant la technique appelée Proper Generalised Decomposition (PGD). Les résultats ont été confrontés à ceux obtenus par des méthodes plus conventionnelles et aussi à des données expérimentales. / Realising composite parts mainly dedicated to aeronautics sector by robotised processes tends to increase productivity. However, Automated Tape Placement (ATP) is still in ripening stage and requires various developments, particularly in case of thermoplastic composites or dry fibres. In this manuscript, we propose different tools allowing to determine in advance the best heating power to drape composite fibres. Difficulty arises from dependence on many parameters, from material (density, specific heat) or from process itself (velocity, number and orientation of plies). Therefore we construct a multi-parametric reduced order model using the so called Proper Generalised Decomposition (PGD). Results have been faced to those obtained from more conventional methods but also to experimental data.

Les travaux présentés dans cette thèse concernent l'étude de l'évolution des régimes d'écoulements diphasiques lors de l'évaporation progressive dans un canal horizontal. Le but est de mieux comprendre l'écoulement à l'intérieur d'un tube récepteur d'une centrale solaire à concentration à génération directe de vapeur. Cette technologie, présentée comme une amélioration des systèmes actuels pouvant permettre une réduction des coûts, consiste en la production de vapeur directement sous l'effet du rayonnement solaire concentré. La prévision de l'écoulement liquide-vapeur alors généré dans le tube est encore de nos jours difficile, c'est pourquoi le recours à la simulation numérique est intéressant. Pour cela un modèle a été développé permettant la simulation de ces écoulements, depuis le début de la création de la vapeur jusqu'à l'existence de larges poches. Basé sur le modèle diphasique VOF du code Fluent, par l'ajout de fonctions personnalisées et d'une phase dispersée supplémentaire, il permet de modéliser différents phénomènes liés au processus d'évaporation : création en paroi, transport, recondensation et création de larges structures. Ce développement a été mis en oeuvre pour simuler des écoulements en évaporation, permettant de reproduire l'évolution des régimes d'écoulement. La validation est faite grâce à une étude expérimentale de la littérature, en comparant les régimes d'écoulements obtenus pour différents débits de liquide et sous l'effet de différents flux de chaleur. Enfin, le modèle a été appliqué à la simulation de la génération de vapeur dans le tube récepteur d'une centrale solaire, mettant en évidence l'apparition et l'évolution des différents régimes d'écoulement. Au vu du peu d'installations expérimentales trouvées dans la littérature sur le sujet, et afin de valider au mieux les fonctions développées, une installation expérimentale a été conçue et dimensionnée. / This PhD thesis is about the study of two-phase flow patterns evolution during progressive evaporation in horizontal tubes. The goal is to better understand the flow regimes inside a receiver tube of a concentrated solar power plant with direct steam generation. This technological evolution allows vapor production directly inside the solar field, which can lead to coast reductions. A two-phase liquid-vapor flow occurs inside the tube, which is currently still difficult to predict. Numerical simulation is an interesting way to investigate these complex phenomena. A model has been developed in order to simulate the flow patterns, from first vapor generation to large vapor slugs. It is based on Fluent software's two-phase VOF model, to which are added user-defined functions and a new dispersed phase. Different phenomena linked to the evaporating process are taken into account: vapor creation at the wall, its transport, recondensation and large structures creation. The model is used to simulate evaporating flows, and retrieves well two-phase flow patterns evolution. Validation is made using experimental data from the literature, by comparing flow regimes obtained for different flow rates and heat fluxes. Finally numerical simulation of direct steam generation inside a concentrated solar plant receiver is conducted, clearly showing apparition and evolution of two-phase flow patterns. Because few experimental data where found in the literature concerning evaporating two-phase flows visualization, a new experimental apparatus has been conceived and sized in order to better validate our numerical results.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

La variation spatiale du mouvement sismique a des effets significatifs sur la réponse dynamique des structures de génie civil de grandes dimensions. Dans la pratique courante, l'excitation du mouvement sismique le long de la fondation de la structure est considérée uniforme, approche cependant inadéquate pour les structures de large portance au sol localisées à proximité des failles ou sur des sites présentant une structure du sous-sol latéralement hétérogène. Cette thèse se propose donc de comprendre les facteurs clefs contrôlant localement la variabilité spatiale du mouvement sismique, avec en ligne de mire la mise en place de recommendations en vue d'incorporer ces effets dans l'estimation de l'aléa sismique et le dimensionnement des structures. La première partie de cette thèse s'intéresse à la composante intra-évènement de l'écart-type de la distribution du mouvement sismique en champ proche à l'aide de simulations numériques du mouvement sismique pour des sources étendues présentnat une cinématique de rupture réaliste. Les résultats suggèrent que la variabilité intra-évènement dépend significativement du type de rupture, cette variabilité augmentant avec la distance pour les ruptures bilatérales et diminuant pour les ruptures unilatérales. La seconde partie traite de la caractérisation de la composition du champ d'onde dans la vallée de Koutavos-Argostoli, qui est une vallée de petite dimension et d'épaisseur sédimentaire faible, située sur l'île – sismiquement active - de Céphalonie en Grèce. Les champs d'onde générés par 46 séismes, ayant des magnitudes variant entre 2 et 5 et des distances épicentrales jusqu'à 200 km, ont été analysés à partir de l'enregistrement par deux réseaux denses de capteurs sismologiques. L'algorithme de traitement d'antenne MUSIQUE est utilisé pour extraire la vitesse, l'azimut, le type et la polarisation des ondes dominantes se propageant à travers le réseau. Les résultats montrent clairement d'importantes diffractions d'ondes de surface aux bords de vallée au-delà de la fréquence de résonance de la vallée. Tandis que les ondes de Love dominent clairement le champ d'ondes proche de la fréquence de résonance, les ondes de Rayleigh dominent à plus haute fréquence dans des gammes de fréquences étroites. Par ailleurs, un excellent accord est observé entre les champs d'onde de surface diffractés localement et les caractéristiques d'amplification du site. La “cohérence décalée” de la partie la plus énergétique du signal a été quantifiée pour chaque paire de stations du réseau. En général, la cohérence calculée sur les composantes horizontales diminue avec la distance entre stations et la fréquence. La cohérence sur la composante verticale indique des valeurs relativement fortes à haute fréquence. Les valeurs de cohérence apparaissent très faiblement corrélées à la magnitude, l'azimut et la distance épicentrale du séisme, mais sont au contraire liées aux caractéristiques géométriques de la vallée. La coherence est systématiquement plus forte pour les couples de capteurs orientés selon la direction parallèle à l'axe de la vallée, et moins forte pour des couples de capteurs orientés dans la direction perpendiculaires. Cette observation est en accord avec les résultats du traitement d'antenne: la prédominance d'ondes de surface diffractées sur les bords de la vallée conduit à des mouvements en phase le long de la direction parrallèle à l'axe de la vallée. Les résultats de cette thèse apporte des elements de compréhension sur la variabilité spatiale du movement sismique et ouvrent de nombreuses perspectives d'application. / Spatial variation of earthquake ground motion over short distances significantly affects the dynamic response of engineered structures with large dimensions. In current practices, the ground motion excitation across the foundation of a structure is assumed to be spatially uniform, which becomes inadequate for spatially extended structures in the near-fault region or on sites with lateral inhomogeneity. This PhD thesis seeks to understand the key parameters that locally control the ground motion spatial variability with the intent of putting forth practical propositions for incorporating such effects in seismic design and hazard assessment. The first part of the thesis addresses the within-event component of the standard deviation of ground-motion distribution in near source region by means of numerical simulation of ground motions for extended sources with realistic rupture kinematics. The results suggest that the within-event variability significantly depends on the rupture type, depicting an increase with distance for bilateral ruptures and a decrease for unilateral ruptures. The second part deals with the characterization of seismic wave field at the Koutavos-Argostoli site, a small-size, shallow, alluvial valley located in the seismically active Cephalonia Island in Western Greece. The seismic wave field was investigated from the recordings of a dense seismological array for a set of 46 earthquakes, with magnitude 2 to 5 and epicentral distance up to 200 km. The MUSIQUE array analysis algorithm was used to extract the phase velocity, back-azimuth, type and polarization of the dominant waves crossing the array. The results clearly indicate dominant scattering of seismic surface waves, mainly from the valley-edge directions, beyond the fundamental frequency of the valley. While Love surface waves clearly dominate the wave field close to the resonance frequency, Rayleigh waves strongly dominate only in relatively narrow frequency bands at higher frequency. Besides, an excellent consistency is observed between the dominance of the identified surface wave type in the wave field and the site amplification. The "lagged coherency" of the most energetic part of the ground motion has been quantified for each station-pair within the array. In general, spatial coherency estimated from the horizontal components exhibit decays with frequency and interstation distance. Estimates from the vertical component exhibit rather larger values at some higher frequencies. Although coherency does not show any consistent trend indicating dependence on the magnitude, back-azimuth or site-to-source distance of the event, it seems to be primarily controlled by the site geometry. Larger coherency is systematically observed when the station pair is oriented parallel to the valley axis, while lower values are observed in the perpendicular direction. This observation proves to be consistent with the MUSIQUE analysis results: the predominance of scattered surface waves propagating across the valley implies an in-phase motion along valley-parallel direction and out-of-phase motion along valley-perpendicular direction. The findings of the present research are expected to contribute in enhancing our understanding of spatial variability of ground motion and improving the coherency models used in engineering. This work also opens up new insights and many questions in need of further investigation.

Cette thèse s’inscrit dans le cadre du projet ANR « Uncertain flow optimization » (UFO) consacré au développement et à l’application de méthodes efficaces de quantification d’incertitudes pour l’analyse et l’optimisation d’écoulements. Dans ce contexte, ces méthodes sont appliquées à des gazosiphons de grande échelle utilisés comme pompe. Plus particulièrement, on s’intéresse à de tels gazosiphons choisis pour constituer l’organe central d’un système géothermique innovant,  compatible avec un environnement urbain. On souhaite en quantifier le potentiel énergétique par voie numérique avec la recherche d’un compromis entre justesse des résultats et efficacité optimale. La simulation de l’écoulement diphasique produit dans le gazosiphon est fondée sur un modèle quasi-1D à flux de dérive et s’appuie sur une démarche de résolution implicite. Les résultats sont validés sur les études expérimentales les plus pertinents de la littérature, dont aucune toutefois n’atteint les longueurs requises de l'ordre du kilomètre. Le code de simulation du gazosiphon fait ensuite l’objet d’une démarche de prise en compte d’incertitudes physiques et de modélisation, précédée par une analyse de deux méthodes de quantification d’incertitude : une méthode non-intrusive de type chaos polynomial, et une méthode plus récente dite semi-intrusive qui fut développée en amont du projet UFO. Cet outil est intégré dans une modélisation simplifiée du système géothermique urbain dans son ensemble impliquant les composants en surface, notamment le compresseur d'air. Il en résulte une optimisation énergétique robuste préliminaire de deux variantes du système géothermique urbain proposé, respectivement de récupération de chaleur et de production d’électricité. / This PhD thesis is part of the ANR project « Uncertain Flow Optimization » (UFO). The project is devoted to the development and application of efficient uncertainty quantification methods for flow analysis and optimization. In this framework, these methods are applied to the study of a large-scale airlift pump. The airlift pump is selected to be part of an innovative geothermal system, which can be exploited within an urban environment. We wish to quantify and optimize the energy potential of this new system with numerical tools. They provide both good accuracy and efficiency properties. The airlift two-phase flow simulation is based on a quasi one-dimensional drift flux model, which is implicitly solved. The solver is validated by comparison with relevant experimental airlift studies from the literature. However, these studies remain below the kilometric-targeted pipe length. Thanks to the analysis of two uncertainty quantification methods, a non-intrusive approach relying on polynomial chaos expansion and a new semi-intrusive method developed ahead of the UFO project, we perform airlift pump simulations taking into account physical and modelling uncertainties. This numerical tool is inserted into a simplified model of the complete urban geothermal system that involves surface devices, such as an air compressor. Finally, a robust preliminary optimization process is performed for two versions of the proposed geothermal urban system. They are designed respectively for heat recovery and electricity production.

La production pariétale de bulles de gaz et son impact sur la dynamique de la phase liquide en canal vertical est étudiée numériquement et expérimentalement. Dans un premier temps, un modèle de mélange 2D stationnaire est utilisé pour décrire l’évolution moyenne des panaches de gaz. Grâce à cette approche, un modèle de couche limite a pu être développé et a permis l’identification des nombres adimensionnels pertinents (analogues aux nombres de Rayleigh et de Prandtl pour la thermique) afin de caractériser les écoulements à bulles dispersées. Dans un second temps, un modèle Eulérien-Lagrangien 3D instationnaire, prenant en compte le couplage quadrilatéral (interactions bullesliquide et bullesbulles) est résolu par Simulation Numérique Directe (DNS) et permet ainsi une description plus fine de l’écoulement à l’échelle de la phase dispersée. Enfin, ces approches numériques sont complétées par des mesures de Spectroscopie d’Impédance Electrochimique (SIE) lors de la production de dihydrogène et de dioxygène par électrolyse alcaline. Les modèles d’écoulement proposés ici montrent globalement un très bon accord avec les résultats expérimentaux tirés de la littérature. Les approches homogènes et DNS présentent toutefois quelques disparités sur l’évaluation du taux de vide dans certaines conditions. Parallèlement, les mesures et simulations de SIE ont montré être clairement affectées par les évolutions du panache de bulles, les spectres d’impédance ont notamment mis en évidence une contribution basse fréquence fortement dépendante de la nature de la phase dispersée (taille de bulle et lois de dispersion). Les trois approches (modèle homogène, DNS et SIE) menées conjointement sont donc fortement complémentaires. Elles permettent non seulement une meilleure compréhension de la physique de l’écoulement diphasique, mais offrent aussi une capacité d’analyse de la pertinence des modèles existants tout en ouvrant la voie à leurs futures améliorations / The wall production of gas bubbles and its impact on the liquid dynamics in a vertical channel is studied by means of numerical simulations and experimentation. First, a 2D stationary mixture model is used to describe the averaged plumes evolutions. Through this approach, a boundary layer model has been developed and identified dimensionless numbers (Raleigh-like and Prandtl-like) characteristic of bubbly flows. Secondly, a 3D non-stationary four-way coupled (with bubblesliquid and bubblesbubbles interactions) Eulerian-Lagrangian model is solved by Direct Numerical Simulation (DNS) and allows a finer description of the two-phase flows at bubble-scale. Finally, the numerical methods are completed by Electrochemical Impedance Spectroscopy (EIS) measurements during hydrogen and oxygen production by alkaline electrolysis.The two-phase flow models are in good agreement with experimental results from literature. There are still some disparities between the homogeneous model and the DNS about the void fraction calculation under certain conditions, though. In the meantime, both EIS measurements and simulations were clearly affected by bubbles plume evolutions, the impedance spectra highlighted a low frequency contribution highly sensible to the nature of the dispersed phase (bubble size a dispersion laws). The 3 approaches (homogeneous model, DNS and EIS) used collectively are strongly complementary. They allow not only a better comprehension of the physics of the two-phase flow, but also serve the analysis of existing models while leading the way for further improvements