L’intérêt scientifique porté depuis plusieurs décennies au développement des plateformes carbonatées tropicales s’explique par la complexité des facteurs de contrôle de la sédimentation, qui montre une influence de la tectonique, du climat, de l’eustatisme, de l’hydrodynamisme etc… Les plateformes carbonatées ne sont pas seulement des environnements enclins à accumuler des sédiments, mais sont aussi de prolifiques « usines » à carbonates qui produisent davantage que ce qu’elles peuvent stocker. Les sédiments en excès sont transférés vers le bassin profond. Ainsi, les transitions plateforme – bassin illustrent les interactions entre production carbonatée in-situ, transfert, et accumulation sédimentaire issue de la re-sédimentation. Dans le sud-ouest de l’Albanie, la série carbonatée du Crétacé Supérieur est composée de dépôts sédimentaire de plateforme et de bassin. Ces successions sont à présent intégrées dans le système chevauchant de la chaîne Dinarides - Hellénides, dont la mise en place est contemporaine d’une phase de déformation orogénique Plio-Holocène. La paléogéographie de la région péri-Adriatique au Crétacé Supérieur révèle une juxtaposition de méga-plateformes et de bassins adjacents dont le remplissage est en partie conditionné par la remobilisation gravitaire des sédiments de la bordure de plateforme. Les investigations sont menées sur cinq zones d’étude. Elles intègrent des descriptions faciologiques macro- et microscopiques pour chacune des huit successions étudiées. Les unités définies font l’objet d’une cartographie à l’échelle de la zone d’étude qui précise l’architecture stratigraphique du système. La datation des dépôts repose sur une étude biostratigraphique, elle est renforcée par les données chronostratigraphiques des isotopes du strontium. Sur la plateforme, les descriptions faciologiques sont appuyées par des méthodes complémentaires de pétrographie (cathodoluminescence, épifluorescence, microscopie électronique à balayage), de pétrophysique (porosimétrie à injection mercure) et de géochimie (isotopes stables du carbone et de l’oxygène) visant à révéler la dynamique de sédimentation et les propriétés réservoir de la succession.Les conditions de sédimentation sont précisées, i) sur la plateforme, dix faciès spécifiques sont identifiés, représentatifs de contidions de dépôt supra-, inter-, et subtidales. Ils s’intègrent dans des motifs d’empilement distinctifs (small-scale sequences) qui attestent le caractère cyclique de la sédimentation, contrôlé par des fluctuations haute-fréquence et basse amplitude du niveau marin relatif; ii) dans le bassin, quatorze faciès sont classés suivant le mécanisme de transport sédimentaire dominant. La répartition spatiale des dépôts permet d’identifier la provenance et les sources préférentielles du matériel calci-clastique. L’évolution tectono-sédimentaire plateforme – bassin au Crétacé Supérieur est déterminée par deux séquences: i) du Cénomanien au Turonien, la sédimentation de plateforme est caractérisée par l’aggradation plus de 700 mètres de faciès intertidaux organisés en séquences (small-scale sequences) émersives ou sub-émersives. Aucun transfert significatif n’est attesté dans le bassin adjacent, témoignant d’une relative stabilité du système évoluant dans un contexte subsident ; ii) l’intervalle Coniacien – Santonien marque l’installation d’une plateforme à rudistes favorable au transfert sédimentaire. Ce dernier se traduit par une progradation marquée des dépôts gravitaires dans le bassin pendant le Campanien. Cette dynamique de transfert est accentuée par la mise en place de slumps résultants du démantèlement tectonique de la bordure de plateforme au Campanien Supérieur et Maastrichtien. / The evolution of tropical carbonate platforms depends on complex interacting factors influencing the sedimentation, such as tectonism, climate, eustacy, hydrodynamism etc… Due to this complexity, it has been of scientific interest for decades. Carbonate platforms are not only prone to accumulate sediments, but also represent prolific carbonates “factories” producing more than they can store on their tops. Excess sediments are shed basinward. Platform-to-basin transitions exemplify interactions between in-situ carbonate production, transfer and sedimentary accumulations resulting from re-sedimentation. In south-west Albania, the Upper Cretaceous carbonate series are made up of platform and basinal deposits. Sedimentary successions are presently integrated in the Dinarides-Hellenides fold-and thrust belt that originate from a Plio-Holocene phase of the Alpine Orogeny. The paleogeographic setting during the Late Cretaceous reveals a juxtaposition of mega-platforms and adjacent basins, partly filled with sediments derived from the shelf edge. Our scientific investigations focused on five study areas. Macro- and microfacies descriptions are provided for eight platform and basin successions. Sedimentary units are defined and mapped on each study area, supporting the stratigraphic reconstruction of the system. Dating is based on biostratigraphy and supported by strontium-isotope data. Regarding platform deposits, facies descriptions are seconded by complementary methods of petrography (cathodoluminescence, epifluorescence, scanning electron microscopy), petrophysics (mercury intrusion porosimetry), and geochemistry (stable carbon and oxygen isotopes), aiming to precise the sedimentation dynamics and reservoir properties of the succession. The sedimentation conditions are specified: i) within the platform, ten specific facies are identified, precising the depositional setting that comprise supra-, inter- and subtidal environments. These facies are integrated in distinctive stacking patterns (small-scale-sequences) pointing to a cyclic sedimentary dynamic controlled by high-frequency and low-amplitude sea level changes; ii) the basinal deposits comprise fourteen facies classified according to the dominant grain-support mechanism. Spatial distribution of the deposits allowed identifying the provenance and preferential sources of calciclatic sediments. The tectono-sedimentary evolution of the platform-to-basin system during the Late Cretaceous can be subdivided into two sequences: i) from the Cenomanian to the Turonian, the platform sedimentation is characterized by a substantial aggradation (700 meters) of intertidal small-scale sequences. No significant transfer was evidenced in the adjacent Ionian Basin, pointing to stable conditions in a subsiding context; ii) the Coniacian-Santonian interval evidences the establishment of a rudist platform massively shedding sediments basinward during the Campanian. A clear progradation of gravity-flow deposits is attested during this period. During the Upper Campanian and Maastrichtian, this transfer is accentuated by the setting of tectonically-triggered slumps resulting from the dismantling of the platform edge.

La relation scalaire entre contrainte et déformation obtenue par le cisaillement simple dans les rhéomètres classiques n'est pas assez riche pour décrire les écoulements complexes. Pour obtenir plus d'information, nous superposons deux écoulements orthogonaux en utilisant la géométrie plan-plan. Le fluide, sous forme cylindrique, peut alors être cisaillé par la rotation mais aussi écrasé par le rapprochement (ou étiré par l'éloignement) des disques. Nous détaillons les calculs théoriques permettant de déterminer les liens entre contraintes et taux de cisaillement et les efforts et vitesses macroscopiques associés. Ensuite, nous décrivons précisément le dispositif expérimental mis en place pour imposer toutes sortes d'écoulements combinant des cisaillements stationnaires ou oscillants, en rotation ou en écrasement. Puis nous présentons les résultats de la comparaison entre l'écoulement de rotation et l'écoulement d'écrasement. Nous présentons enfin les expériences de superposition des deux écoulements. Nous créons des écoulements complexes et divers afin, entre autres, de mesurer et comprendre la loi d'écoulement 3D et le critère d'écoulement 3D des fluides à seuil / The scalar relation between stress and strain obtained by simple shear within classical rheometers is not rich enough to describe complex flows. To obtain more information, we shall superimpose two orthogonal flows and for that, we shall use the parallel-plate geometry. The volume of fluid has a cylindrical shape and can be sheared by the rotation of a plate or squeezed by the move closer (or stretched by the move away) of the plates. We detail theoretical considerations that allow us to determine the link between stresses and shear rates and associated macroscopic efforts and velocities. Then we precisely describe the experimental setup developed to impose all kind of flows : steady or oscillating flows, rotation or squeeze. After that, we present the results of comparison between rotational and squeeze flow. Finally, we present the experiments with both superimposed flows. We create complex flows to measure and understand, among others, 3D behavior and 3D flow criteria of yield stress fluids

Les écoulements diphasiques dans les centrales de type réacteur à eau pressurisée appartiennent à des régimes très variés allant du faible nombre de Mach jusqu'aux ondes de chocs. Calculer des solutions approchées précises de ces écoulements peut s'avérer délicat dans certains régimes. On s'intéresse dans cette thèse à la conception et à l'étude de méthodes numériques robustes et stables à grand pas de temps, capables de calculer des solutions approchées précises quel que soit le régime d'écoulement, y compris sur maillage grossier. Une stratégie pour construire de tels schémas consiste à : utiliser un schéma semi implicite basé sur un splitting d’opérateurs pour séparer la résolution approchée des phénomènes rapides de celles des phénomènes lents ; corriger les flux numériques afin d’améliorer la précision du schéma dans certains régimes. Deux approches sont utilisées pour analyser la capacité du schéma numérique à gérer plusieurs régimes d'écoulement. L’approche des schémas asymptotic preserving est utilisée pour traiter le système de la dynamique des gaz avec termes sources raides. On utilise ensuite la notion de schéma tout-régime pour le système de la dynamique des gaz et les systèmes diphasiques homogénéisés HRM et HEM à bas nombre de Mach. Des propriétés garantissant la stabilité et la robustesse des schémas ont été obtenues, et en particulier des inégalités d'entropie discrètes. L'implémentation de ces méthodes a permis de mener des expériences numériques en 1D et 2D sur maillage non structuré qui confirment le gain en précision et en temps de calcul des schémas asymptotic preserving et tout-régime ainsi construits par rapport à des schémas numériques classiques. / Two-phase flows in Pressurized Water Reactors belong to a wide range of Mach number flows. Computing accurate approximate solutions of those flows may be challenging from a numerical point of view as classical finite volume methods are too diffusive in the low Mach regime. In this thesis, we are interested in designing and studying some robust numerical schemes that are stable for large time steps and accurate even on coarse meshes for a wide range of flow regimes. An important feature is the strategy to construct those schemes. We use a mixed implicit-explicit strategy based on an operator splitting to solve fast and slow phenomena separately. Then, we introduce a modification of a Suliciu type relaxation scheme to improve the accuracy of the numerical scheme in some regime of interest. Two approaches have been used to assess the ability of our numerical schemes to deal with a wide range of flow regimes. The first approach, based on the asymptotic preserving property, has been used for the gas dynamics equations with stiff source terms. The second approach, based on the all-regime property, has been used for the gas dynamics equations and the homogeneous two-phase flows models HRM and HEM in the low Mach regime. We obtained some robustness and stability properties for our numerical schemes. In particular, some discrete entropy inequalities are shown. Numerical evidences, in 1D and in 2D on unstructured meshes, assess the gain in term of accuracy and CPU time of those asymptotic preserving and all-regime numerical schemes in comparison with classical finite volume methods.

Dans ce projet de recherche, on s'intéresse au développement et à l'évaluation de nouvelles méthodes numériques pour les écoulements peu profonds. De nouvelles techniques de discrétisation spatiales et temporelles des équations sont proposées. La première partie de la thèse est dédiée au développement d'une méthode des volumes finis explicite d'ordre élevé et d'une famille de schémas semi-implicites qui sont efficaces pour la modélisation des processus lents et rapides dans les écoulements océaniques et atmosphériques. La deuxième partie du projet de recherche concerne la construction d'un schéma numérique efficace sans solveur de Riemann pour les écoulements peu profonds avec une topographie variable sur un maillage non structuré. Dans cette partie de la thèse, une nouvelle approche est proposée pour l'analyse de stabilité des schémas numériques non structurés pour les équations en eaux peu profondes. Dans la troisième partie de la thèse, deux schémas de volumes finis sont développés pour les lois de conservation sur des surfaces courbes qui ont un large potentiel d'être appliqués aux écoulements peu profonds sur la sphère. Dans ces cas, les schémas numériques sont développés en adoptant la démarche suivie par Stanley Osher. Cette démarche consiste à utiliser des systèmes hyperboliques simples qui génèrent des phénomènes d'ondes complexes et des solutions qui ont différentes structures. Ces solutions sont très efficaces pour tester les méthodes numériques. Dans notre cas, nous avons utilisé les équations de Burgers qui ont joué un rôle très important dans le développement des schémas numériques à capture de chocs en mécanique des fluides. / This research project focuses on the development and evaluation of numerical methods for shallow flows by proposing new spatial and temporal discretization techniques. First, a new high-order explicit finite volume method and a class of semi-implicit schemes are introduced which are effective for modelling fast and slow waves in oceanic and atmospheric flows. In the second part of the research project, a central-upwind scheme is proposed for shallow water flows on variable topography using unstructured grids. In this part of the project, a new approach is proposed for the stability analysis of unstructured numerical schemes for shallow water equations. In the third part of the thesis, two finite volume methods are developed for the conservation laws on curved geometries which are potentially applicable to shallow flows on a sphere. For such cases, numerical schemes are developed by using the approach followed by Stanley Osher. This approach employs simple hyperbolic systems which generate complex wave phenomena, and solutions that are effective for assessing numerical methods. In our case, Burgers’ equations are used since they have played an important role in the development of shock-capturing schemes in fluid mechanics.

Actuellement, les codes de calcul de composants thermohydrauliques de réacteurs nucléaires du CEA et d'EDF utilisent des modèles physiques diphasiques de mélange à 3 ou 4 équations. Or, il existe un fort besoin industriel pour des modèles physiques plus sophistiqués tels que le modèle diphasique à 6 équations voire des modèles multichamps. Par ailleurs, le code système CATHARE du CEA, un des codes systèmes les plus utilisés aujourd'hui sur le plan international, utilise un modèle physique diphasique à 6 équations et un schéma numérique semi-implicite de type ICE à maillages décalés. Le schéma de CATHARE est connu pour sa robustesse dans une large gamme de configurations d'écoulement. En s'inspirant de l'expérience de CATHARE, on propose de mettre en œuvre un schéma volumes finis colocalisés de type « pressure based ». Le but est d'obtenir un décentrement des flux qui assure la robustesse du schéma tout en gardant une bonne précision. De plus, le fait de pouvoir utiliser des maillages colocalisés (structurés ou non-structurés) permet de traiter de différentes configurations complexes de cœur de réacteurs et de réaliser des calculs fin d'inter assemblage. Le schéma doit conserver exactement la masse et l'énergie et la solution numérique doit converger lorsque l'on raffine le maillage. Le schéma doit être capable de traiter des cas d'apparition et de disparition des phases, par exemple le cas de la colonne bouillante où il y a changement de phase dû au transfert de chaleur, des cas de tuyères avec changement de phase dû à un élargissement ou un rétrécissement brusque, ou de séparation de phase par gravité. En outre le schéma doit être capable de traiter des configurations de calcul à faible nombre de Mach, par exemple le cas du renoyage d'un cœur de réacteur. L'objet de la thèse consistera à développer une méthode Volumes Finis co-localisés (dans l'esprit de Ghidaglia et al.) et la direction d'investigation s'inspirera des travaux de Jeong et al. qui a conduit au code CUPID. Références. Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605. / Currently, there is a strong industrial need for sophisticated physical models such as the two-phase model with 6 equations or multi-field models for the thermo-hydraulic calculation of nuclear reactor components. In addition, the system code CATHARE of CEA, one of the most utilized system codes at the international level, employs a two-phase 6 equations and a semi-implicit numerical scheme of ICE type on staggered grids. The code CATHARE is known for its robustness in a wide range of flow configurations. Drawing on the experience of CATHARE, we propose to implement a collocated finite volume pressure based scheme. The aim is to obtain a decentering of flux which ensures the robustness of the scheme while keeping good accuracy. In addition, being able to use collocated grids (structured or unstructured) can handle different complex configurations. The scheme should conserve exactly the mass and energy and the numerical solution needs to converge when the mesh is refined. The scheme should be able to handle cases of phase appearance and disappearance, for example in the case of boiling column where phase change is due to heat transfer, the case of nozzles with phase change due to a widening or abrupt narrowing, or phase separation by gravity. In addition the scheme should be capable of calculation configurations at low Mach number. The purpose is to develop a co-located Finite Volume method (in the spirit of Ghidaglia et al.) And direction of investigation is build on work of Jeong et al. which led to CUPID code. References: Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605.

Dans le contexte des écoulements micro- et nano-fluidiques, ce travail porte sur l'étude des interactions à l'interface entre des flux de gaz rares et des surfaces métalliques via une approche de modélisation multi-physique et multi-échelle. Elle se concentre tout particulièrement sur l'interaction entre l'argon et une surface d'or. Pour ce faire la modélisation a été effectuée en deux étapes, une première partie utilisant la mécanique quantique à l'échelle atomique et une deuxième partie de dynamique moléculaire à l'échelle nanométrique. La première partie est consacrée à l'obtention de potentiels d'interaction entre un atome d'argon et les atomes d'or de la surface par des méthodes de calculs théoriques basés sur la DFT comportant des effets à longues distances. Deux approches, donnant des résultats comparables, ont été utilisées : la première est liée à la description périodique de la surface d'or par un modèle basé sur la description des électrons par des ondes planes alors que la seconde permet de récupérer séparément les parties répulsives et attractives de l'interaction d'un atome d'argon avec un petit cluster d'or. Ces potentiels d'interactions ont été décomposés en potentiels de paires Ar-Au utilisables par des simulations de dynamique moléculaire. Ces simulations ont consisté en la projection d'atomes d'argon sur des surfaces d'or ‘parfaites' dites lisses ou des surfaces rugueuses plus représentatives de la technologie actuelle. L'analyse statistique des vitesses réfléchies permet de déterminer le coefficient d'accommodation tangentiel de l'argon sur des surfaces d'or. Ce coefficient est la traduction du phénomène de glissement qui peut ainsi être modélisé dans une description plus macroscopique de l'écoulement d'un gaz dans une micro-conduite. L'approche multi-physique utilisée dans ce travail a permis la détermination numérique de coefficients d'accommodations tangentiels très précis et comparables à l'expérience pour le couple argon-or, et doit pouvoir être appliquée à d'autres couples / In the context of micro- and nano-flows, this work concentrates on the study of interactions at the interface of noble gas and metal surfaces by a multi-physics and multiscale model. Particularly, the interaction of an argon atom with a gold surface is the focus of the study. The work has been made in two steps: the first one occurred at the atomic scale in which Quantum Mechanics is employed and the second one at the nanoscale with the use of Molecular Dynamics.The first part of the work was devoted to the determination of interaction potentials between an argon atom and gold atoms from the surface by DFT calculation methods comporting long range effects. Two approaches, leading similar results, have been used: the first one is linked to a periodic description of the gold surface where electrons are defined by plane waves, the second one gives independently repulsive and attractive parts of the interaction of an argon atom with a small gold cluster. Those interaction potentials are then decomposed in pair potentials suitable for Molecular Dynamics simulations. These last ones consisted in multiple times projecting argon atoms on smooth or rough gold surfaces (which are more representative of the roughness of actual technologies). The statistical analysis of the reflected velocities yielded the tangential momentum accommodation (TMAC) coefficient of argon on gold surfaces. This coefficient is the transcription of slip phenomena which occur at the interface, and it can then be used in nano-flow simulations. The multi-physics approach of the thesis gives accurate TMAC values which are comparable to experiments. The accounted method could then be applied to other noble gas metal surface couples

Le lien entre la force fluide s'appliquant sur un obstacle et les structures cohérentes est évident ; il n'existe pourtant pas de relation quantitative entre chaque structure de l'écoulement et la force qu'elle produit. Nous étudions cette problématique selon deux grands axes : d'abord la définition d'une densité volumique d'effort, puis une définition des structures cohérentes adaptée à cette densité de force. Nous montrons que pour caractériser les structures cohérentes dans un domaine englobant un obstacle, la formulation utilisée ne doit pas faire apparaître de terme de pression et ne doit pas comporter d'information sur les contours du domaine. Nous étudions plus précisément deux formulations de la force et montrons que seule la formulation diffusive faisant intervenir une impulsion du laplacien de vorticité est adaptée à un calcul volumique des efforts. En appliquant cette formulation à un écoulement 2D autour d'un cylindre dans le régime de von Kármán, nous pouvons ensuite étudier les contributions à la force de deux parties de l'écoulement définies à partir du laplacien de vorticité : les couches de vorticité et la zone de recirculation, dans un domaine restreint autour du cylindre. Les couches de vorticité captent ainsi les effets visqueux, tandis que la zone de recirculation capte les effets de l'enroulement des couches limites en tourbillons. Il est intéressant de noter que la zone de recirculation n'englobe pas les tourbillons développés du sillage, l'influence de ces tourbillons étant en fait captée par la densité de force dans une petite zone en aval du cylindre. L'application 3D de la formulation diffusive nécessite une très bonne définition de la vorticité dans la zone proche paroi, et nous montrons que cette approche ne peut pas encore être facilement mise en oeuvre sur des géométries complexes avec les moyens numériques actuels.

Cette thèse comporte trois parties. <br />Dans la première, nous présentons une méthode de simulation d'écoulements fluide/particules. Nous montrons que la pénalisation du tenseur des contraintes, associée à une discrétisation en temps par la méthode des caractéristiques, conduit à une formulation variationnelle de type Stokes généralisée. Des tests numériques sont effectués sous FreeFem++ afin d'étudier la convergence. Nous en présentons également trois exemples d'utilisation. <br />Dans la seconde partie nous proposons un modèle permettant de prendre en compte les forces de lubrification dans les simulations directes d'écoulements fluide/particules. Nous présentons d'abord un modèle de contact visqueux dans le cas particule/plan, obtenu comme limite, à viscosité nulle, du modèle de lubrification. Nous décrivons ensuite un algorithme reposant sur une étape de projection des vitesses, à chaque instant, sur un espace dit de vitesses admissibles. On montre alors la convergence du schéma et on généralise l'algorithme au cas multi-particules. Nous en présentons également un exemple de programmation orientée objet. <br />Dans la dernière partie, nous considérons un système discret de sphères (boulier en 1D) qui interagissent à travers la force de lubrification. Le modèle microscopique repose sur le développement de cette force à courte distance. Nous proposons une équation constitutive marcoscopique, de type Newtonien, reposant sur une viscosité linéique proportionnelle à l'inverse de la fraction locale de fluide. Nous établissons la convergence du modèle microscopique vers le modèle macroscopique proposé.

L'objectif de c e mémoire est d'apporter une contribution à la modélisation et la simulation numérique de la convection thermosolutale de mélanges binaires de gaz parfaits contenus dans des cavités. Un modèle a été élaboré en se basant sur l'approximation de faible compressibilité. Le premier chapitre précise la démarche suivie dans la modélisation et une formulation originale en est déduite afin de traiter les différents types de conditions aux limites et de conditions de références hydrostatiques analysés dans le mémoire. Les variations de masse volumique sont déduites de la loi des gaz parfaits et la pression thermodynamique est calculée à partir de la conservation de la masse totale. La méthode numérique repose sur la méthode des volumes finis mise en uvre sur des maillages décalés. Le couplage vitesse-pression est traité par un nouvel algorithme dont l'efficacité est discutée en détail. La démarche numérique est validée via des comparaisons avec des solutions de références, en régime stationnaire comme en régime transitoire pour des écoulements transitionnels. Dans la seconde partie du mémoire, on considère d'abord la convection thermosolutale dans une cavité rectangulaire verticale dans le cas où les écoulements sont induits par des gradients horizontaux de température et de concentration. On discute en particulier les limites de l'approximation d'extrême dilution. La condensation de vapeur d'eau et l'évaporation d'un film d'eau liquide sur les parois d'une cavité sont ensuite étudiées en régime transitoire. Ces changements de phase surfaciques sont associés à la convection naturelle dans une cavité dont les températures des quatre parois varient au cours du temps / The aim of this dissertation is at modeling and numerically simulating thermosolutal convection within cavities filled by binary gas mixtures of ideal gases. A new problem formulation, based on the weakly compressible approximation, has been derived in order to account for the changes in density as well as thermodynamic pressure. The ideal gas law and global conservation of mass are invoked for the predictions of density field and thermodynamic pressure. The first part of the manuscript is devoted to the mathematical derivation of the governing equations and boundary conditions, numerical procedure implemented and, checks of the accuracy of the results through deeply comparisons with updated benchmark solutions. The emphasis is put on the efficiency of algorithm used for solving the pressure-velocity coupling for unsteady, transitional flow regimes. Thermosolutal convection without phase changes at the cavity surfaces is first considered in the second part of the manuscript. We investigated the case of vertical cavities with horizontal gradients of concentration and temperature. The results clearly show that the dilute approximation fails to be valid as soon as the maximum concentration difference exceeds a critical value function of the flow parameters. Surface condensation of water vapor or evaporation of liquid water film in vertical cavities is then considered. The specific thermal boundary conditions of uniform but time-varying temperatures of the four walls are considered. Periodic variations of the flow, temperature and relative humidity fields are discussed in detail. The evolutions of thicknesses of the water film over the four walls are also predicted

Les bâtiments résidentiels et industriels munis d'un réseau de ventilation constituent des installations complexes, susceptibles d'être le siège de transferts de masse et d'énergie variés, selon les situations de fonctionnement. Afin d'étudier ces transferts de masse, une méthodologie permettant d'établir des expérimentations à échelle réduite pour l'étude des écoulements isothermes, en régime permanent ou transitoire, a été développée. Cette méthodologie a été validée numériquement et expérimentalement sur des configurations simples, puis appliquée à deux configurations de référence, représentatives de celles rencontrées dans le domaine nucléaire.L'influence du vent sur les transferts de masse au sein de ces configurations, en situation de fonctionnement normale, dégradée (arrêt de la ventilation) ou accidentelle (surpression interne), a été étudiée dans la soufflerie climatique Jules Verne du CSTB. Les effets du vent, couplés ou non à une surpression interne, peuvent alors entraîner une perte partielle ou globale du confinement des polluants au sein des installations. De plus, la turbulence du vent peut induire des inversions instantanées des débits de fuite, qui ne sont pas identifiées en régime permanent. Par ailleurs, l'analyse de sollicitations transitoires montre la faible influence de l’inertie des branches sur les écoulements transitoires, pour des grandeurs caractéristiques d'une installation réelle. Enfin, des essais de traçage gazeux ont été réalisés afin d’étudier la dispersion d’un polluant au sein d’une configuration de référence soumise aux effets couplés du vent, de la ventilation mécanique et d'une surpression interne.La robustesse du code à zones SYLVIA, utilisé notamment pour appuyer les évaluations de sûreté des installations nucléaires, a été analysée à partir de ces résultats expérimentaux. La prise en compte des phénomènes physiques observés expérimentalement a été validée, en régimes permanent et transitoire. Toutefois, quelques limitations ont été identifiées pour l'étude de la dispersion d'un scalaire passif, du fait des hypothèses utilisées dans le code SYLVIA, comme dans tout code à zones (concentration homogène dans les locaux, propagation instantanée dans les branches et dans les locaux). / Residential and industrial buildings equipped with a ventilation system are complex facilities, where various heat and mass transfers could occur according to the operating conditions. In order to study these mass transfers, a methodology has been developed so as to carry out reduced-scale experiments for the study of isothermal flows, in steady or transient state. This methodology has been numerically and experimentally validated on simple configurations, and then applied to two standard configurations, representing nuclear facilities.The wind influence on mass transfers inside these configurations, in normal, damaged (stopping ventilation) or accidental (internal overpressure) situations, has been studied in the Jules Verne climatic wind tunnel of the CSTB. The wind effects, coupled or not with an internal overpressure, can lead to a partial or a total loss of the pollutant's containment inside buildings. Moreover, the wind turbulence can bring about instantaneous reversal leakage flowrates, which cannot be identified in steady state. In addition, the study of transient phenomena has highlighted the low influence of the branch inertia on transient flows, for typical values of real facilities. Finally, tracer tests have been carried out in order to study the pollutant dispersion inside a standard configuration subjected to wind, mechanical ventilation and internal overpressure effects.The reliability of the zonal code SYLVIA, used notably to support safety assessment in nuclear buildings, has been analyzed from these experimental results. The modelling of the physical phenomena experimentally observed has been validated, in steady and transient states. However, limitations have been identified for the study of pollutant dispersion, due to hypothesis used in SYLVIA code, as in all zonal codes (homogenous concentration inside rooms, instantaneous propagation inside branches and rooms).