Lors de sa détente dans une turbine, la vapeur subit une chute d'enthalpie qui entraîne sa condensation spontanée sous forme d'un nuage de gouttelettes submicroniques. Ces gouttes vont se déposer sur les aubes aval et y former un film d'eau. Ce dernier est arraché sous l'effet de l'écoulement de vapeur environnant ce qui crée des gouttes de quelques dizaines de microns qui peuvent se redéposer sur les aubes aval. Ces phénomènes sont à l'origine de pertes, généralement regroupées sous le terme générique de "pertes par humidité", estimées grâce à la loi de Baumann. Le but de cette thèse est de mettre en place dans un code méridien des modèles permettant la prise en compte des phénomènes de condensation, de déposition et d'écoulement des films liquides afin de pouvoir estimer les pertes liées à chacun d'entre eux. Dans cette optique un modèle diphasique homogène permettant d'avoir accès à la fraction massique de liquide et au nombre de gouttes est implanté dans le code méridien. Ce modèle est validé sur un cas test expérimental de détente en tuyère et alimente le modèle de déposition. Les contributions de la diffusion, de la turbophorèse, de la thermophorèse, de la gravité et de l'inertie des gouttes à la déposition sont étudiées. Il apparaît nécessaire de prendre en compte la diffusion, la turbophorèse et la déposition inertielle sur les bords d'attaque des aubes. Un modèle permettant d'avoir accès à l'épaisseur et à la vitesse d'un film liquide soumis au cisaillement d'un écoulement environnant, à la force de frottement sur la paroi et aux effets de la rotation est mis en place. Ce modèle est validé par rapport à des résultats expérimentaux d'écoulement de film liquide sur une plaque plane dans des conditions proches de celles rencontrées en turbine à vapeur basse pression. Finalement, un calcul réalisé sur une géométrie réelle de turbine basse pression de 8 étages permet de démontrer l'applicabilité de la méthodologie mise en place sur un cas industriel. Les contributions des phénomènes étudiés aux pertes par humidité sont explicitées.

Les mousses cimentaires sont des matériaux solides poreux qui ont despropriétés isolantes intéressantes. Elles sont utilisées dans le domaine de laconstruction pour réduire la consommation énergétique des habitations. Cependant,lors de la prise cimentaire, certaines mousses se déstabilisent et s'effondrent. Laformulation de tels matériaux présente donc des difficultés techniques etscientifiques notamment du point de vue des interactions entre les nombreuxcomposés utilisés, organiques et minéraux. L’objectif de ce travail de thèse a doncété d’étudier spécifiquement les interactions entre la matière minérale et le tensioactifutilisé pour la mise en forme alvéolaire, afin de comprendre le rôle de l’agentmoussant dans les mécanismes de stabilisation des mousses. Pour ce faire, nousavons considéré un système réactif simplifié, mais proche de l’utilisation courante,composé d’un ciment classique, d'eau et d’un tensioactif dosable par spectroscopieUV-visible. Pour découpler les différents effets, nous avons étudié ce système à troiséchelles : la mousse cimentaire, le coulis non moussé et des films uniques. / Cement-based mineral foams are solid porous materials and presentvery interesting insulating properties. Such materials are used in the constructionsector in order to reduce the energy consumption of houses. However, theproduction of these materials can be difficult because of the kinetic competitionbetween foam stability and cement setting. Cement-based mineral foams areextremely complex materials because of the number of components and all theirinteractions. The objective of our study is to get insight into the interactions betweenmineral and organic compounds and to understand their role in the stabilizingmechanisms. Our strategy was to study a model system, sufficiently close the realone to be representative, and sufficiently simple to identify the contribution of thedifferent components. We used as a foaming agent a surfactant exhibiting acharacteristic peak in UV-visible spectroscopy in order to determine its function andposition in the complex system. We first worked on unfoamed slurries, then westudied thin liquid films and the foamed materials.

Cette thèse expérimentale porte sur des impacts de gouttes en caléfaction, aussi appelée régime Leidenfrost. Dans ce cas, la goutte est isolée thermiquement et mécaniquement du substrat surchauffé par une fine couche de vapeur. On s'affranchit ainsi de la friction visqueuse. Les substrats présentent des textures micrométriques localisées. On retrouve un régime autosimilaire d'étalement aux temps courts. On caractérise des régimes de recouvrement d'un défaut ponctuel par la goutte. Ces régimes sont dictés par l'épaisseur de lamelle par rapport à celle du défaut. Les défauts génèrent des excroissances dans leur sillage dont la dynamique peut être approchée selon deux modèles inertio-capillaires ; l'un valable aux temps courts, l'autre aux temps plus longs. En présence de plusieurs défauts, on fragmente la lamelle depuis plusieurs sites selon le ratio entre largeur des défauts et épaisseur de la lamelle. On simule par un algorithme de pavage le motif final sur lequel le fluide se concentre à la fin de la fragmentation. Grâce à des rugosités plus complexes on peut canaliser l'étalement de la goutte. On exacerbe alors l'étalement maximal dans l'axe des canaux. On peut aussi inhiber l'étalement par des textures circulaires. Les rugosités affectent le temps avant rebond de la goutte. On exhibe une dépendance générale unique entre temps avant rebond et étalement maximal. La dépendance est valable lorsque les textures exacerbent l’étalement, aussi bien que lorsqu’elles l’inhibent / The presented work deals with drop impacts in Leidenfrost regime. In such a case, the drop is thermally and mechanically isolated from the overheated substrate by a thin vapor layer. Viscous friction can therefore be neglected. The substrates are shaped with localized micrometric textures. We rediscover a self-similar spreading regime at short times. We characterize covering regimes of a single defect by the drop. Those regimes are driven by the ratio between lamella thickness and defect thickness. Defects give rise to excrescences in their wake, whose dynamic can be approached by two inertial-capillary recoil models. One is valid at short times and the other at longer times. In the case of several defects, we break up the lamella from different sites according to the ratio between defect width and lamella thickness. We numerically predict, with a tessellation algorithm, the pattern on which the fluid is localized at the end of the fragmentation. Through more complex textures, we can channel the drop spreading. The spreading is increased in the directions of the channels. One can also inhibit the spreading with circular textures. The textures affect the time before drop rebound. We exhibit a general and unique dependency between time before rebound and maximal spreading. This dependency is valid when textures increase the spreading as well as when they inhibit it

Lors de sa détente dans une turbine, la vapeur subit une chute d’enthalpie qui entraîne sa condensation spontanée sous forme d’un nuage de gouttelettes submicroniques. Ces gouttes vont se déposer sur les aubes aval et y former un film d’eau. Ce dernier est arraché sous l’effet de l’écoulement de vapeur environnant ce qui crée des gouttes de quelques dizaines de microns qui peuvent se redéposer sur les aubes aval. Ces phénomènes sont à l’origine de pertes, généralement regroupées sous le terme générique de ”pertes par humidité”, estimées grâce à la loi de Baumann. Le but de cette thèse est de mettre en place dans un code méridien des modèles permettant la prise en compte des phénomènes de condensation, de déposition et d’écoulement des films liquides afin de pouvoir estimer les pertes liées à chacun d’entre eux. Dans cette optique un modèle diphasique homogène permettant d’avoir accès à la fraction massique de liquide et au nombre de gouttes est implanté dans le code méridien. Ce modèle est validé sur un cas test expérimental de détente en tuyère et alimente le modèle de déposition. Les contributions de la diffusion, de la turbophorèse, de la thermophorèse, de la gravité et de l’inertie des gouttes à la déposition sont étudiées. Il apparaît nécessaire de prendre en compte la diffusion, la turbophorèse et la déposition inertielle sur les bords d’attaque des aubes. Un modèle permettant d’avoir accès à l’épaisseur et à la vitesse d’un film liquide soumis au cisaillement d’un écoulement environnant, à la force de frottement sur la paroi et aux effets de la rotation est mis en place. Ce modèle est validé par rapport à des résultats expérimentaux d’écoulement de film liquide sur une plaque plane dans des conditions proches de celles rencontrées en turbine à vapeur basse pression. Finalement, un calcul réalisé sur une géométrie réelle de turbine basse pression de 8 étages permet de démontrer l’applicabilité de la méthodologie mise en place sur un cas industriel. Les contributions des phénomènes étudiés aux pertes par humidité sont explicitées. / During the expansion in a turbine, the enthalpy of the steam fall. This fall leads the steam to cross the saturation line and brings about its spontaneous condensation and the appearance of fog droplets. The deposition of these droplets on downstream blades lies at the root of the creation of a liquid film. This film is torn off by steam flow and creates coarse water. These big droplets can impact downstream blades. Each of these phenomena induces some losses which are generally grouped in the ”wetness losses” estimated thanks to Baumann’s rule. The aim of this work is to develop, in a throughflow code, some models which allow to take into account the condensation, deposition and liquid film flow and the losses linked to each of these phenomena. An homogeneous two phases flow model is implemented in the throughflow code. Two transport equations on the mass fraction of liquid and on the droplets number are added to Euler’s equations written for gas phase. This model is validated on an experimental test case of expansion in nozzle and feeds the model of deposition. The influences of diffusion, turbophoresis, thermophoresis, gravity and of inertia of droplets on the deposition are studied. It appears to be necessary to take into account diffusion, turbophoresisand deposition due to inertia on the leading edges of blades. A model which allows to evaluate the thickness and the velocity of a liquid film submitted to aerodynamic shear, friction on blades and rotational effects has been developed. This model is validated on experimental results of a liquid film flow on a flat plate in some conditions representatives of those encountered in low pressure steam turbines. Finally a calculation realised on a real geometry of a 8-stages low pressure steam turbine has demonstrated that the methodology developed during this work can be used on an industrial test case. The contributions of the phenomena studied to wetness losses are evaluated.

Les caractéristiques des transferts de masse et de chaleur par convection forcée lors de l'évaporation du film liquide dans le canal en présence d'un cylindre carré poreux parcouru par un écoulement transversal sont étudiées numériquement. L'objectif principal de la présente étude est d'évaluer l'effet de l'introduction d'un cylindre carré poreux sur le transfert de chaleur et de masse. Plus précisément, cette étude examine l'influence de paramètres tels que l'humidité relative de l'air ambiant, la température de l'air à l'entrée, le flux de chaleur imposé, la position du cylindre, le taux de blocage et le nombre de Reynolds. Une comparaison entre les deux configurations, avec et sans cylindre carré poreux a été réalisée pour mettre en évidence l'effet de son ajout. Pour modéliser ce phénomène nous avons résolu l'équation classique de convection forcée et le modèle de Darcy-Brinkman-Forchheimer dans les milieux poreux. Nous avons trouvé que l'insertion d'un cylindre carré poreux dans le canal pourrait perturber l'écoulement et améliorer de manière significative les taux de transfert de masse et de chaleur au niveau des parois du canal. Les transferts de chaleur et de masse deviennent plus importants avec la diminution du nombre de Darcy et lorsque l'obstacle poreux est placé au milieu du canal. Par contre, l'augmentation de l’humidité relative de l’air ambiant et de la température d'entrée réduit le transfert de masse. A Da = 10-6, l'écoulement ne pénètre plus à travers le cylindre poreux, la structure de l'écoulement est similaire à celui d'un cylindre carré plein. Enfin, nous proposons des lois de corrélations qui permettent de prédire les valeurs des nombres de Sherwood et de Nusselt en fonction des nombres de Reynolds, de Biot et du taux de blocage. / The characteristics of mass and heat transfers by forced convection during liquid film evaporation in the channel with a built in porous square cylinder in a cross flow are investigated numerically. The main objective of the present study is to evaluate the effect of introducing a porous square cylinder on the heat and mass transfer. Specifically, this study examines the influence of parameters such as the relative humidity of the ambient air, the air inlet temperature, the imposed heat flux, the variation of the cylinder position, blockage ratio and Reynolds number on the performance of the evaporation at the channel wall. A comparison between the two configurations, with and without, porous square cylinder has been performed to highlight the effect of its addition. To achieve this, we solved the classic equation of forced convection and the Darcy-Brinkman-Forchheimer model in the porous media. We find that the insertion of a porous square cylinder in the channel could make the flow more disturbed and significantly improve mass and heat transfer rates at the channel walls. The heat and mass transfer enhancements is greater with a decrease of the Darcy number and for γ=1 when the porous obstacle is placed in the middle of the channel. It is also greater with a decrease of the temperature and relative humidity of the air at the inlet. At Da = 10-6, the flow does not penetrate through the porous cylinder, the flow pattern is similar to that of a solid square cylinder. Finally, we propose correlations that allow us to define the Sherwood and Nusselt numbers based on the Reynolds, Biot numbers and the blockage ratio.

L'OBJECTIF DE CE TRAVAIL DE THESE EST DE METTRE EN EVIDENCE LE ROLE DU TRANSFERT THERMIQUE DANS LA PAROI CHAUFFANTE LORS DE L'EBULLITION D'UN LIQUIDE STAGNANT A LA SATURATION. CET EFFET, HABITUELLEMENT NEGLIGE DANS LES MODELES D'EBULLITION, PEUT EXPLIQUER CERTAINES CARACTERISTIQUES DU CYCLE D'EBULLITION ET DE L'ACTIVITE DES SITES DE NUCLEATION. IL A ETE MONTRE QUE LE CYCLE D'EBULLITION PEUT ETRE DECRIT PAR DEUX PHASES: (1) LORS DE LA CROISSANCE DE LA BULLE, LA TEMPERATURE DE LA PAROI DIMINUE GLOBALEMENT PAR EVAPORATION DE LA MICROCOUCHE A LA BASE DE LA BULLE; (2) LA SURCHAUFFE INITIALE EST RETABLIE PRINCIPALEMENT PAR CONDUCTION RADIALE DE LA CHALEUR DANS LA PAROI. ON PEUT ALORS DONNER UNE EXPLICATION DE L'EVOLUTION DE LA TEMPERATURE DE LA PAROI MIS EN EVIDENCE EXPERIMENTALEMENT A L'AIDE DE CRISTAUX LIQUIDES THERMOCHROMIQUES ET DE L'INFLUENCE DE L'ANGLE DE CONTACT SUR LES ECHANGES THERMIQUES.

Les turbines à gaz doivent satisfaire des normes d'émission polluantes toujours en baisse. La formation de polluants est directement liée à la qualité du mélange d'air et de carburant en amont du front de flamme. Ainsi, leur réduction implique une meilleure prédiction de la formation du spray et de son interaction avec l'écoulement gazeux. La Simulation aux Grandes Échelles (SGE) semble un outil numérique approprié pour étudier ces mécanismes. Le but de cette thèse est de développer des modèles phénoménologiques décrivant la phase liquide notamment le film et son atomisation en bout de lèvre d'injecteur, dans un contexte SGE. Ces modèles sont validés ou calibrés sur l'expérience académique réalisée par l'Institut für Thermische Strömungsmaschinen (ITS) de l'université technologique de Karlsruhe (KIT), et appliqués dans une configuration réelle de moteur d'hélicoptère. Dans un premier temps, le film liquide mince est décrit par une approche Lagrangienne. Les particules de film représentent un volume élémentaire de liquide adhérent à la paroi. L'équation du mouvement est donnée par l'intégration des équations de Saint-Venant sur l'épaisseur du film. La dynamique du film est donnée par le gradient de pression longitudinal, le cisaillement interfacial du gaz et la gravité. Dans un second temps, l'atomisation du film est caractérisée par la distribution de taille de gouttes du spray généré. Celle-ci est décrite par une distribution de Rosin-Rammler dont les coefficients sont paramétrés par la vitesse du gaz, la tension superficielle du liquide et l'épaisseur de la lèvre de l'injecteur. Les constantes de ce modèle, baptisé PAMELA, sont calibrées sur l'expérience du KIT-ITS. La simulation de l'expérience KIT-ITS permet de valider le modèle de film, de vérifier la robustesse du modèle PAMELA, et de comparer qualitativement l'angle du spray généré. L'application de ces modèles dans une configuration réelle partiellement instrumentée permet de valider PAMELA sans modification de ses constantes, et d'étudier leur impact sur la structure de flamme, comparé aux méthodes traditionnelles d'injection liquide.

L'objectif principal de ce travail de thèse concerne l'étude de l' élimination de la DCO et de la nitrification dans une nouveau lit bactérien Multi-Section ( MSB ) . Après une caractérisation de l’hydrodynamique et du transfert d’oxygène de ce lit bactérien, les expériences biologiques menées sous des conditions opératoires contrastées (fortes et faibles charges organiques eteaux usées contenant ou pas des matières particulairs) ont été menées. En parallèle, des simulations avec le logiciel Biowin® ont été réalisées. Les principaux résultats sont résumés en suivant :- La rétention de liquide statique est majoritaire par rapport à la rétention dynamique que ce soit en présence ou en absence de biofilm. Le biofilm joue le rôle d’une "éponge" permettant un maintien de l’humidité du lit même à faible débit. Les expériences de DTS ont montré que le biofilm accroit le temps de séjour du liquide et conduit à une diminution de l’épaisseur du film liquide permettant ainsi de promouvoir le transfert de l'oxygène.- Le réacteur MSB montre une élimination efficace de la DCO (> 95 % ) et de la nitrification ( > 60 % de l’azote entrant), mais une accumulation de DCO particulaire a lieu dans le filtre ce qui conduira à un colmatage à terme. La nitrification cohabite avecl’élimination de la DCO même dans la première section et pour une charge organique élevée ce qui implique une bonne capacité d’oxygénation du MSB par l’aération naturelle.- Un modèle dynamique de MSB a été utilisé implémenté sur le simulateur - BioWin , afin d'obtenir la répartition des biomasses au sein du réacteur et d'évaluer le processus limitant dans chaque section. Le modèle partiellement calibré peut aider à estimer les besoins minimum d'oxygène pour la nitrification et peut rendre compte de la compétition entre la croissance hétérotrophe et la nitrification. / The main objective of this PhD work focused on the study of the COD removal and nitrification in a new designed Multi-Section Bioreactor (MSB). Hydrodynamic characterization of the reactor, biological experiments under contrasted conditions and simulations by Biowin® software were carried out:- Firstly, it was found that static liquid retention is the predominant part both without and with the presence of biofilm. Biofilm acts like a "sponge". RTD experiments showed that biofilm can promote liquid residence time, decrease the liquid film andpromote the oxygen transfer consequently.- Secondly, the MSB operated at contrasted organic loading rate (OLRs) and nitrogen loading rate (NLRs) showed that COD can be effectively removed (removal efficiency > 95%) and nitrification (> 60% of the N removal) occurred in this biofilter.Nitrification is efficient even in the first section implying no drastic oxygen limitation though only natural aeration is occurring.- Thirdly, a TF dynamic model has been used from a simulator - BioWin, in order to get more insights on the biomass distribution in the pilot and to assess the limiting process in each section of the bioreactor. Calibration of the model can help us to estimate theminimum oxygen requirement for nitrification for each zone inside the pilot and it can well represent the competition between heterotrophic growth and nitrification.

Le présent travail apporte une contribution à la compréhension des mécanismes des transferts combinés de chaleur et de masse en convection naturelle lors de l’évaporation d’un film liquide d’eau d’épaisseur négligeable dans un canal vertical ondulé. L’écoulement est laminaire et bidimensionnel. Les équations régissant le phénomène sont résolues à l’aide d’une méthode aux volumes finis et le traitement du couplage vitesse-pression est réalisé par la méthode de projection. Les influences de la densité de flux de chaleur, de la température ainsi que l’humidité de l’air à l’entrée et la forme de la paroi du canal sur les transferts sont étudiées. Les résultats sont présentés sous la forme de ligne de courant, d’isothermes et d’iso-concentrations.Les simulations numériques effectuées ont permis l’étude détaillée de la structure de l’écoulement ainsi que des champs thermiques et massiques. Nous représentons également, les nombres de Nusselt et de Sherwood. / The present work is a contribution to the understanding of the mechanisms of combined heat and mass transfers in natural convection during the evaporation of a liquid film with negligible thickness in a wavy vertical channel. The flow is laminar and two-dimensional. The equations governing the phenomenon are resolved using the finite volumes method and the treatment of the coupling between velocity and pressure is carried out by the projection method. The influences of the heat flux density, the temperature and the humidity of the inlet air and the shape of the channel wall on the transfers are studied. The results are presented in the form of cstreamlines, isotherms and iso-concentrations.The numerical simulations carried out have allowed the detailed study of the flow structure as well as the thermal and mass fields. We also represent the Nusselt and Sherwood numbers.

Ces travaux de thèse s'incrivent dans le cadre du traitement de gaz acides et captage CO2 dans les colonnes à garnissages structurés. Les gaz à traiter réagissent avec un liquide s'écoulant à contre-courant sur des plaques métalliques dont la compléxité géométrique permet d'accroître l'aire d'échange, et donc l'efficacité du procédé. Dans un contexte de modélisation multi-échelles des contacteurs à garnissages structurés, les écoulements gaz-liquide à la plus petite échelle géométrique des plaques de garnissages (de l'ordre de l'épaisseur du film liquide) sont étudiés, pour améliorer la compréhension et la modélisation des écoulements diphasiques et phénomènes de mouillage dans les garnissages. L'objectif final est de développer une méthodologie CFD pour reproduire des écoulements diphasiques 3D sur des géométries complexes telles que les plaques de garnissages. Pour ce faire, il est nécessaire de progresser en méthodes numériques et de proposer des méthodes expérimentales pour observer des écoulements de film liquide sur des géométries complexes. Ces travaux comprennent une partie numérique et une partie expérimentale. Un écoulement sur une plaque de garnissage structuré peut présenter des zones sèches, et donc des lignes de contact (dynamiques), ce qui présente un défi en simulation numérique à cause des différentes échelles de l'écoulement. La méthodologie employée ici en simulation numérique consiste à résoudre l'écoulement jusqu'à une échelle intermédiaire en modélisant les effets des plus petites échelles. Le code de calcul Two-Phase Level-Set a été utilisé et modifié dans ce but. Différentes méthodes level-set ont d'abord été testées de manière à identifier une méthode satisfaisante quant à la réduction des erreurs de conservation de masse, un problème rencontré en level-set. Il est ici montré que certaines combinaisons de schémas de discrétisation spatiale et temporelle permettent de réduire considérablement ces erreurs de conservation de masse. Après avoir réalisé de nombreux tests de validation, une nouvelle méthode numérique est proposée pour simuler les grandes échelles d'écoulements diphasiques 3D avec ligne de contact dynamique en level-set, dans des conditions réalistes. La méthode est ici validée pour des écoulements axisymétriques de gouttes simulés en 3D, en régime visqueux et en régime inertiel, et pour des écoulements de gouttes sur plan incliné. Les résultats sont en très bon accord avec d'autres travaux numériques et expérimentaux. Afin de faciliter l'utilisation de cette méthodologie pour des applications industrielles, un modèle sous-maille similaire a été implémenté dans un code VOF commercial; les résultats sont aussi en très bon accord avec d'autres travaux. En plus de ces développements numériques, une campagne expérimentale est mise en oeuvre pour observer des écoulements de film liquide sur une plaque de garnissage structuré. Les méthodes expérimentales employées sont d'abord testées et validées pour des écoulements de film plat ou ondulé sur plan incliné, et ensuite utilisées pour observer des écoulements de film sur des plaques de garnissage. L'épaisseur de film liquide est mesurée aux creux et aux crêtes des picots des plaques de garnissages, pour différents débits, par imagerie confocale chromatique. Des lois de puissance de l'épaisseur de film en fonction du Reynolds sont proposées; celles-ci sont très différentes suivant la position des relevés de mesure, aux creux ou aux crêtes des picots. La vitesse à l'interface de l'écoulement gaz-liquide est aussi mesurée, par PIV et PTV, en utilisant des particules hydrophobes. Les résultats montrent que le liquide a tendance à dévier du creux des canaux (corrugations), et la norme de la vitesse semblent présenter des extremums correspondant aux creux et crêtes des picots. [...] / The work described in this thesis is motivated by the use of structured packing columns in acid gas treatment and post-combustion CO2 capture. In a counter-current mode, flue gases react with the liquid that flows down over metal sheets, the geometrical complexity of which allows increasing the specific interfacial area, and thereby the overall efficiency of the process. In the context of multiscale modeling of structured-packing contacting devices, the focus in this work is on the gas-liquid flows at the smallest geometrical scale of packing sheets, of the order of the liquid film thickness, aiming to improve understanding and modeling of two-phase flows and wetting phenomena in structured packings. The ultimate objective is to build up a CFD methodology to reproduce 3D two-phase flows over complex surfaces such as structured packing sheets. For this purpose, progress is necessary both in pertinent computational methods and in the adaptation of experimental methods for observing liquid film flows over complex surfaces. This thesis therefore consists of computational and experimental parts. Flows over structured packing sheets may exhibit dry zones, and hence (moving) contact lines, the numerical simulation of which presents a computational challenge due to the disparity in length scales involved. Here, the methodology for large-scale numerical simulations of flows with moving contact lines consists in resolving the flow down to an intermediate scale and modeling effects of smaller ones. The parallelized freeware Two-Phase Level-Set has been extended for this purpose. First though, because some level-set methods have been reproached to yield mass conservation issues, an assessment is made of the mass conservation properties of a range of level-set methods. It is demonstrated that the combined use of some spatial and temporal discretization schemes allows to drastically reduce mass conservation errors in level-set methods. Having thus implemented a level-set method with satisfactory performance at such tests (and others), a novel numerical method is proposed to perform 3D large-scale simulations of flows with moving contact lines in level-set, under realistic conditions. Validation tests of axisymmetric droplet spreading in a viscous, and in an inertial regime, simulated in 3D, and sliding drops are shown to be in excellent agreement with prior experimental and numerical work. The results show that complex contact-line dynamics observed in prior experimental studies on sliding droplets can be simulated using the present large-scale methodology. To facilitate dissemination of this work in industrial applications, a similar subgrid model has been implemented in a commercial volume-of-fluid code; results of validation tests are shown to be in excellent agreement with other work. These computational developments are accompanied by an experimental campaign to observe liquid film flows over structured packing sheets. All experimental methods used herein are tested and validated for flat and wavy films down an inclined plane before being used for observing liquid film flows over packing sheets. The film thickness is measured at local troughs and crests of small-scale corrugations of the structured packing sheet, for different flow rates, by Chromatic Confocal Imaging. Power laws of the Reynolds number for the mean liquid film thickness are suggested, with significant differences for measurements at crests compared to that at troughs. Interface velocity measurements are also performed by PIV and PTV using hydrophobic particles. Results reveal that the liquid tends to deviate from troughs of large-scale corrugations, and seems to exhibit local extrema of the velocity magnitude corresponding to troughs and crests of small-scale corrugations. [...]