Les écoulements diphasiques en fracture se produisent dans de nombreuses applications industrielles et environnementales. Pour modéliser ces écoulements, il est important de connaître le rôle des caractéristiques géométriques des surfaces constituant les fractures.L’objectif de ce travail de thèse est d’étudier expérimentalement les écoulements diphasiques et la dispersion de traceurs dans des fractures réelles dont on connaît la cartographie des surfaces et d’en proposer une modélisation adéquate.L’écoulement diphasique inertiel a été modélisé à l’aide de la généralisation de la loi de Darcy couplée à l’approche de la fonction F. Le modèle proposé permet de prédire les perméabilités relatives et les saturations en fonction des débits des deux fluides avec deux paramètres d’ajustement.Les profils de concentration obtenus lors des expériences de dispersion ont été interprétés avec différents modèles afin d'évaluer leur capacité à décrire la dispersion dans les fractures et donc l’hétérogénéité de celles-ci. Les résultats montrent que le modèle stratifié, avec un seul paramètre estimé, donne des résultats satisfaisants et permet de quantifier l'hétérogénéité.Enfin, avec cinq paramètres estimés, à partir de la courbe de la concentration de traceur à la sortie des fractures, le modèle CTRW permet de décrire avec une bonne précision les profils de concentration de toutes les abscisses / The two-phase flows through fractures occur in many industrial and environmental processes.Modelling these flows needs understanding the role of geometric characteristics of the fractures surfaces. The objective of this thesis is to study experimentally the two-phase flow and the dispersion of tracers in fractures which the aperture map is known and propose a suitable model.The inertial two-phase flow model is based on the generalization of the single-phase full cubiclaw that accounts for non-Darcian effects by using the F function approach. On this approach,the presence of a second fluid is taken into account through a multiplier function introduced into the superficial velocity of each fluid. This model let us predict the relative permeability and the saturation as a function of the two fluid flows with two adjustment parameters.The concentration profiles obtained during dispersion experiments were interpreted with different models to study their ability to describe the dispersion in fractures and hence the heterogeneity of them. The results show that the classical convection-dispersion is not appropriate due to the heterogeneity of the fractures. On the other hand, the stratified model,with only one estimated parameter, gives satisfactory results and allows quantifying heterogeneity. Finally, with five estimated parameters, from the concentration curve at the outlet of the fractures, the CTRW model can describe the concentration profiles of all abscissas with good accuracy

Ce travail s’inscrit dans le cadre de nos recherches sur les écoulements à bulles. C’est l’étude expérimentale du mélange d’un traceur très peu diffusif (fluoresceine dans l’eau) dans l’écoulement instationnaire engendré par le passage d’un écoulement uniforme au travers d’un réseau d’obstacles sphériques (2% de fraction volumique) répartis aléatoirement dans l’espace. Cet écoulement reproduit correctement les caractéristiques de l’agitation dans un essaim de bulles en ascension. La vitesse du fluide est mesurée par PIV de manière assez classique. Le traceur est injecté en amont du réseau et l’´evolution de sa concentration est mesurée par PLIF. L’utilisation de la PLIF pour mesurer des champs de grande extension (15 cm) et avec une grande dynamique d’intensité lumineuse fluorescée constitue une contribution importante de ce travail. Les résultats ainsi obtenus montrent que, à petit nombre de Reynolds, le mélange est régi par les forts gradients de vitesse qui existent dans le voisinage des sphères. A grand nombre de Reynolds, il est maîtrisé par la turbulence qui se développe alors. L’analyse des résultats comporte deux parties principales : (1) une analyse statistique des profils de concentration aboutissant à la détermination d’un coefficient de diffusion effectif et (2) une description de la distribution spectrale des fluctuations de vitesse et de concentration. / This research is a part of our research about bubbly flows. Experiments are performed about mixing of a high Schmidt scalar component (fluorescein in water) by the agitation generated by the flow through a random array of fixed spheres (at high Re and with a volume fraction of solid equal to 2%). This flow mimics for a great part the agitation in the liquid phase of a bubble swarm rising in a liquid otherwise at rest. The velocity of the liquid is estimated from PIV measurements. The scalar is injected through a point source in the array and the evolution of its concentration is estimated by PLIF method. An important part of this research is the measurement of large fields of concentration (15 cm) with a good precision by PLIF. The results show that, at moderate Reynolds number (100), mixing is controled by the steep velocity gradients located near the spheres, while, at large Reynolds number, it is controled by the turbulence that develops. The analysis of the results is composed of two parts : (1) the statistical analysis of the spatial distribution of scalar concentration, and the determination of an effective diffusion coefficient, (2) a spectral analysis of the velocity and concentration fields.

Le travail actuel traite le problème de l'écoulement gaz-liquide compositionnel pour la représentation d'un puits dans des simulations de réservoir. L'objectif est de développer les rapports analytiques qui pourraient lier la pression de puits, la saturation et les concentrations de composant à leurs valeurs moyennes dans chaque zone de l'influence de puits. Nous avons montre que des N-2 équations décrivant le transport des concentrations de phase peuvent être transformées en équations ordinaires (différentiation en ce qui concerne la pression) indépendantes du temps et de l'espace examinant le long des lignes de courant. Ces équations transformées représentent des relations thermodynamiques additionnelles réduisant le degré de liberté thermodynamique. En raison de ceci le variance thermodynamique du modèle compositionnel limite s'avère égal à 1 pour tout nombre de composants chimiques. Cette transformation assurent se découplage total du modèle compositionnel limite dans le nouveau modèle thermodynamique et le modèle hydrodynamique, qui peut être resoved inedpendently d'un un autre. Le modèle thermodynamique décompose est totalement indépendant sur l'hydrodynamique, et décrit le comportement d'équilibre d'un système gaz-liquide ouvert. Ce modèle contient les équations d'équilibre et la EOS classiques, aussi bien que les N-2 nouvelles équations appelées la "Delta-loi", qui déterminent la variation de composition d'un système ouvert dans lequel la masse de chaque composant n'est pas conservée. Le modèle hydrodynamique décompose a été utiliser pour développer les solutions asymptotiques des problèmes d'écoulement de gaz-condensat. Le problème a été montré perturbé singulièrement avec la formation d'une couche limite à voisinage du puits. Dans cette couche la propriété basique de contraste des mobilities de gaz et de liquide est perturbée. On développe une technique spéciale qui permet de construire des expansions asymptotiques sous forme de deux diverses séries: le primer est valide loin du puits (l'expansion extérieure), alors que le second dans valide à voisinage du puits (la couche limite ou l'expansion intérieure). En appliquant la méthode asymptotique suggérée, nous avons développé les solutions asymptotiques pour le problème de l'écoulement multicompositionnel de gaz-condensat àu puits dans un domaine borné à un débit variable. En plusieurs cas la solution peut être obtenue sous la forme analytique, alors qu'en cas général de l'écoulement la méthode mène à une solution de semi-analytical, présentée comme problème initial pour une équation. Cette solution, même étant présenté en forme non-analytique, est beaucoup plus simple que le modèle compositionnel original, car l'équation pour la saturation ne dépend pas de la pression locale, mais dépend de la pression de bord seulement. Dans le dernier chapitre nous avons prolongé cette approche au cas quand la pression capillaire n'est pas négligée. Nous avons supposé cependant que les forces capillaires sont inférieures à la différence de pression entre le puits et la bord de réservoir, dus à ce que nous avons appliqué la méthode de perturbation pour petit nombre capillaire inverse. On obtient les solutions asymptotiques améliorées qui tiennent compte de l'effet capillaire. Simulations numériques montrées que ces effets sont maximaux àu voisinage du puits. Le cas d'une exploitation à long terme du réservoir. Tout d'abord, la simulation traditionnelle du comportement de réservoir peut être effectuée avec l'ECLIPSE en ajoutant la méthode de Peaceman de représentation bonne, qui est une relation analytique pour la pression de puits par l'intermédiaire du débit de production. Cette relation inclut une saturation condensat qui peut être évaluée comme saturation moyenne de réservoir. Une telle simulation fournit un bon résultat pour la pression de puits (ou le débit de production), et un bon résultat pour la saturation de bord, mais des données faibles pour la saturation de puits. Cette valeur peut être calculée en utilisant les solutions asymptotiques suggérées dans le présent projet. Le cas d'un puits de production à court terme (un essai de puits). Il est suffisant de simuler le comportement de réservoir dans le domaine de l'influence de puits, en supposant que l'état de frontière demeure invariable (et connu a priori). Dans ce cas-ci les solutions asymptotiques suggérées dans le travail de presnet peuvent être directement employées pour simuler le problème (sans employer l'ECLIPSE). Le problème de l'écoulement de gaz-condensat à une fracture. Nous avons construit un champ plutôt arbitraire avec des lignes de courant orientées à la fracture, en supposant que la fracture joue le rôle d'une décharge, et les lignes de courant sont stationnaire. Pour une ligne de courant arbitraire nous avons reformulé le modèle d'écoulement de gaz-condensat dans des coordonnées cartésiennes. Pour ce problème nous avons développé les expansions asymptotiques / The present work deals with the problem of the compositional gas-liquid flow for the well representation in reservoir simulations. The objective is to develop analytical relationships which would be able to link the wellbore pressure, saturation and component concentrations to their mean values within each zone of the well influence. It is shown that N-2 equations describing the transport of phase concentrations can be transformed into the space- and time-independent ordinary differential equations (differentiation with respect to pressure) when examined along flow streamlines. These transformed equations represent additional thermodynamic relations reducing the thermodynamic degree of freedom. Due to this the thermodynamic variance of the limit compositional model is shown to be equal to 1 for any number of chemical components. This transformation ensure a total splitting of the limit compositional model into the new thermodynamic model and a hydrodynamic model, which may be resoved inedpendently of one another. The split thermodynamic model is totally independent on the hydrodynamic one, and describes the equilibrium behaviour of an open gas-liquid system. This model contains the classic equilibrium equations and EOS, as well as N-2 new differential equations called the "delta-law" which determine the composition variation in an open system, in which the mass of each component is not conserved. The split hydrodynamic model consists of two equations for pressure and saturation. The split hydrodynamic model was used to develop asymptotic solutions of gas-condensate flow problems. The problem was shown to be singularly perturbed with formation of a boundary layer in the vicinity of the well. In this layer the basic contrast property of gas and liquid mobilities is perturbed. A special technique is developed which enables to construct asymptotic expansions in the form of two various series, one of them is valid far from the well (the exterior expansion), while the second one in valid in the vicinity of the well (the boundary-layer or interior expansion). By applying the suggested asymptotic method, we have developed the asymptotic solutions for the problem of multi-component gas-condensate flow to a well in a bounded domain at a variable flow rate. In several cases the solution may be obtained in the analytical form, while in general case of flow the method leads to a semi-analytical solution presented as an initial problem for a differential equation. This solution, even being presented in non-analytical form, is much simpler than the original compositional model, as the equation for saturation does not depend on the local pressure, but on the boundary pressure only. In the last chapter we extended this approach to the case when the capillary pressure is not neglected. We assumed however that the capillary forces are lower than the pressure difference between the wellbore and reservoir boundary, due to which we applied the perturbation method over the small inverse capillary number. The improved asymptotic solutions are obtained which take into account the capillary effect. Numerical simulations shown that these effects are maximal in the vicinity of the well. For the practice, the obtained asymptotic solutions may be used in the following way to resolve the problem of gas-condensate well representation. The case of a long-term exploitation of the reservoir}. First of all, the traditional simulation of the reservoir behaviour can be performed with ECLIPSE by adding the Peaceman method of well representation, which is an analytical relation for the wellbore pressure via the production rate. This relation includes a condensate saturation which can be evaluated as a mean reservoir saturation. Such a simulation provides a good result for the wellbore pressure (or the production rate), and a good result for the boundary saturation, but poor data for the wellbore saturation. This value can be calculated next by using the asymptotic solutions suggested in the presented project. The case of a short-term well production (a well test). It is sufficient to simulate the reservoir behaviour in the domain of the well influence, by assuming that the boundary state remains invariable (and known a priori). In this case the asymptotic solutions suggested in the presnet work can be directly used to simulate the problem (without using ECLIPSE)

L’objectif est de construire le modèle homogénéisé d’un écoulement diphasique en milieu poreux et fracturé, en mettant en évidence le phénomène de mélange dynamique (mixing) entre les phases, provoqué par l’hétérogénéité du milieu. L’attention est concentrée sur l’influence de la capillarité. L’homogénéisation à double échelle a été appliquée. Le mixing se manifeste sous forme de la dispersion hydrodynamique et de l’advection renormalisée. Le tenseur de dispersion, déterminé à travers le problème cellulaire, est une fonction non linéaire de la saturation, vitesse d’écoulement, rapport de viscosité et du nombre capillaire. Pour les milieux fracturés, une méthode streamline configurations a été avancée pour le cas diphasique. Elle permet d’obtenir la dispersion et la perméabilité effective sous forme analytique pour des réseaux de fracture périodiques, ou semi-analytique pour des réseaux aléatoires. La simulation d’un déplacement diphasique à la base du nouveau modèle a été réalisée / The objective of the thesis is to develop the homogenized model of a two-phase flow through a porous and fractured medium by highlighting the dynamic mixing between the phases, caused by the medium heterogeneity. Attention is focused on the influence of the capillarity. The two-scale homogenization is applied. The mixing is manifested in form of the hydrodynamic dispersion and renormalized advection. The dispersion tensor, determined by the cell problem, is a nonlinear function of saturation, flow velocity, viscosity ratio and capillary number. For a fractured medium the method of streamline configurations was advanced for a two- phase case. This method enables to obtain the dispersion tensor and the effective permeability in analytical form for periodic fractured networks or in semi-analytical form for random networks. The simulation of two- phase displacement based on the new model is performed

Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements instationnaires de deux fluides visqueux non miscibles, séparés par une interface mobile. Plus particulièrement des écoulements sans choc constitués d'une phase gazeuse et d'une phase liquide sont considérés. Pour modéliser de tels écoulements, une approche dans laquelle le gaz est décrit par les équations de Navier-Stokes compressible et le liquide par les équations de Navier-Stokes incompressible est proposée. C'est le couplage de ces deux modèles qui constitue l'originalité et l'enjeu principal de de cette thèse. Pour traiter cette difficulté majeure, une méthode globale (i.e. la même dans chaque phase) et simple à mettre en oeuvre est élaborée. L'utilisation des équations de Navier-Stokes formulées de façon unifiée pour les inconnues primitives (pression, vitesse et température) constitue le point de départ pour la construction de notre méthode qui repose sur les composants suivants: une méthode d'éléments finis stabilisés pour la discrétisation spatiale des équations de Navier-Stokes; une approche Level Set pour représenter précisément l'interface dont l'équation de transport a été résolue par une méthode de type Galerkin Discontinu; et des grandeurs moyennes pour traiter les discontinuités à l'interface. Le bon comportement de notre approche est illustré sur différents tests mono et bi-dimensionnels. / In this work, we are interested in the numerical simulation of instationnary viscous flows of two immiscible fluids, separated by a mobile interface. In particular, flows without shock composed of a gas phase and a liquid phase are considered. In order to modelize such flows, an approach in which the gaz is described by compressible Navier-Stokes equations and the liquid by incompressible Navier-Stokes équations is proposed. The coupling between these two models is the originality and the stake of this thesis. To treat this important difficulty, a global (i.e. the same for each phase) and simple method is elaborated. In our procedure we propose, using the Navier-Stokes equations formulated in set of primitives unknowns (pressure, velocity and temperature), to elaborate a strategy that relies on the follow components: the stabilized finite element method to discretize spatially the Navier-Stokes equations; the Level Set method for tracking the interface precisely with a discontinuous Galerkin method to solve the associated transport equation; and some averaged quantities to treat the discontinuities at the interface. The good behaviour of this approach is performed on both one and two spatial dimensions.

Certaines applications industrielles nécessitent des niveaux d’étanchéité exceptionnels pour permettre la réalisation d’un vide poussé ou pour répondre à des enjeux de sécurité radiologique par exemple. Ces niveaux de haute étanchéité statique sur des assemblages démontables sont obtenus à l’aide de joints entièrement métalliques. La fuite résultante de l’assemblage n’est due qu’à la persistance d’un champ des ouvertures à l’interface entre le joint d’étanchéité et la bride d’assemblage, conséquence d’un contact imparfait entre les deux surfaces rugueuses. Le champ des ouvertures à l’interface de contact est assimilable à une fracture rugueuse hétérogène, de nature multi-échelle, et peut en principe être obtenu par un calcul de déformations mécaniques préalable. Dans ce travail, on s’intéressera plus particulièrement à l’écoulement gazeux raréfié dans le régime glissant au sein de cette fracture. Pour les régimes modérément raréfiés,l’écoulement est modélisé par l’équation de Reynolds faiblement compressible avec correction de glissement de premier ordre aux parois que l’on développe. On effectue ensuite un changement d’échelle par la méthode de la prise de moyenne volumique, permettant d’établir un modèle macroscopique d’écoulement à l’échelle d’un élément représentatif, où le débit massique est relié au gradient de pression via le tenseur de transmissivité. Celui-ci, caractéristique de l’élément représentatif de fracture, est calculé par résolution d’un problème de fermeture et est dépendant de la microstructure ainsi que du libre parcours moyen représentatif sur l’élément. Pour remonter à l’écoulement dans l’ensemble de la fracture, hétérogène à cette échelle, celle-ci est subdiviséeen pavés sur chacun desquels est calculé un tenseur de transmissivité local par la méthode sus-citée. Ensuite, l’écoulement dans ce champ de tenseurs est résolu par une méthode des éléments finis de frontière, donnant la transmissivité apparente glissante du joint dans son ensemble. Cette approche à deux échelles, vue comme outil d’aide à la conception, permet une réduction de la complexité de calcul par rapport à une simulation directe, rendant possible une analyse plus efficace du comportement d’un système d’étanchéité. Pour valider l’utilisation du modèle de glissement d’un point de vue macroscopique et s’affranchir des incertitudes sur le calcul de déformation mécanique, des puces nanofluidiques de type réseau hétérogène de canaux droits sont fabriquées par photolithographie par niveaux de gris. Des essais expérimentaux de mesure de fuite sont réalisés sur ces géométries modèles, représentant des joints idéalisés. Ces essais sont effectués en appliquant une forte différence de pression d’hélium par utilisation d’un spectromètre de masse mesurant la fuite, produisant une condition de vide en sortie de puce.Selon les puces, les régimes de raréfaction atteints vont alors du régime glissant au régime moléculaire. Le débit de fuite mesuré est alors supérieur à celui prédit par le modèle de premier ordre, l’écart restant inférieur à un ordre de grandeur quel que soit le régime

Le dimensionnement des digues paravalanches s’appuie sur la connaissance des processus physiques liés au ressaut, qui se forme lorsqu’un écoulement fin et rapide rencontre un obstacle suffisamment haut pour ralentir et épaissir l’écoulement incident. La hauteur du ressaut est aujourd’hui calculée à partir d’équations strictement valides pour des écoulements de matériaux non frictionnels et non compressibles tels que l’eau, sur fond plat et lisse. Or, les avalanches de neige dense sont des écoulements compressibles qui ne peuvent avoir lieu qu’en pente, et au sein desquels se produit de la dissipation d’énergie par friction et collisions entre les grains. Il est donc essentiel de mieux connaître le comportement des ressauts dans les écoulements granulaires en pente. Pour cela, la thèse s’appuie sur plusieurs approches. Les ressauts granulaires stationnaires sont d’abord étudiés de manière purement théorique, à l’aide des équations de conservation de la masse et de la quantité de mouvement moyennées dans l’épaisseur, afin de trouver une relation générale pour prédire la hauteur des ressauts quelques soient les conditions d’entrée. Nous simulons ensuite numériquement un grand nombre de ressauts granulaires en faisant varier plusieurs paramètres (la pente du plan incliné, le débit, le diamètre des grains, la friction entre les grains) à l’aide de la méthode aux éléments discrets. Cette méthode permet d’accéder à la structure interne des ressauts, et notamment à la mesure des champs de vitesse, de fraction volumique, ou encore de la dissipation d’énergie. Les simulations sont réalisées en deux dimensions. Enfin, un dispositif de mesure innovant, qui utilise la radiographie à rayons X dynamique, a été adapté à une expérience de laboratoire existante pour créer des ressauts granulaires stationnaires. Cette technique de mesure permet, en particulier, de mesurer la distribution spatiale moyennée dans la largeur de l’écoulement de la fraction volumique avant, à l’intérieur et après le ressaut granulaire. La comparaison du nouveau cadre théorique proposé avec les résultats expérimentaux et numériques nous permet de mettre en évidence une grande diversité des types de ressauts granulaires en fonction des conditions initiales. Pour chaque type de ressaut, les lacunes du cadre théorique classique, qui ne tient pas compte des forces mises en jeu dans le ressaut ni de la compressibilité, sont clairement établies. / The design of avalanche protection dams relies on the understanding and modelling of physical processes related to the formation of jumps that form when a thin and fast flow meets an obstacle high enough to slow down and thicken the incoming flow. The jump height is nowadays calculated through equations that are strictly valid for non-frictional incompressible flows on a horizontal and smooth bottom. However, dense-snow avalanches are compressible granular flows taking place on a slope, and inside which energy is dissipated through enduring frictional contacts and collisions between grains. It is then essential to decipher the behaviour of jumps formed during granular flows down inclines. To this extent, the thesis relies on several approaches. Standing granular jumps are first studied in a purely theoretical way, with the help of depth-averaged mass and momentum conservation equations, in order to find a relation to predict the height of the jumps regardless of the input conditions. A great number of granular jumps are then simulated by varying several parameters (the slope angle of the incline, the discharge, the grain diameter, the grain-grain friction) thanks to the discrete element method. This method allows us to access to the internal structure of the jumps, and in particular to the spatial distributions of velocity, volume fraction and energy dissipation. Those simulations are done in two dimensions. Finally, an innovative measurement technique using dynamic X-ray radiography was adapted to an existing small-scale laboratory device to produce standing granular jumps. This technique allows in particular the measurement of the width-averaged spatial distribution of volume fraction before, inside and after the granular jumps. The comparison between the new theoretical framework proposed and both the experimental and numerical data, allows us to evidence a rich variety of granular jump patterns as a function of the input conditions. For each type of jump pattern, the shortcomings of the classical theoretical framework, which does not account for the forces at stake within the jump volume nor the compressibility, are well established.

Ce manuscrit présente une étude couplée, numérique et théorique, portant sur les écoulements tournants transitionnels et turbulents. L'accent y est mis sur la formation de structures macroscopiques cohérentes au sein de l'écoulement, générées par des procédés rendus fortement tri-dimensionnels par la présence des couches limites sur les disques et le long des parois cylindriques extérieure (carter) et/ou intérieure (moyeu). La complexité de ces écoulements pose de véritables difficultés en recherche fondamentale mais les résultats de ces travaux ont aussi une importance non négligeable pour les machines industrielles tournantes, depuis les disque-durs jusqu'aux turbopompes spatiales, la conception de ces dernières étant la motivation première pour ces travaux de thèse. Ce travail peut être divisé en deux sous-parties. Dans un premier temps, les cavités industrielles sont modélisées par de simples cavités rotor/stator lisses pour y étudier la dynamique de l'écoulement. Comme les campagnes expérimentales sur les machines industrielles ont révélé de dangereux phénomènes instationnaires en leur sein, l'accent est mis sur l'obtention et l'étude des fluctuations de pression dans les écoulements modèles. Ensuite, les SGE de trois configurations de turbine industrielle réelle sont réalisées pour étudier les fluctuations de pression in situ et appliquer les diagnostiques éprouvés sur les géométries modèles. / This thesis provides a numerical and theoretical investigation of transitional and turbulent enclosed rotating flows, with a focus on the formation of macroscopic coherent flow structures. The underlying processes are strongly threedimensional due to the presence of boundary layers on the discs and on the walls of the outer (resp. inner) cylindrical shroud (resp. shaft). The complexity of these flows poses a great challenge in fundamental research however the present work is also of importance for industrial rotating machinery, from hard-drives to space engines turbopumps - the design issues of the latter being behind the motivation for this thesis. The present work consists of two major investigations. First, industrial cavities are modeled by smooth rotor/stator cavities and therein the dominant flow dynamics is investigated. For the experimental campaigns on industrial machinery revealed dangerous unsteady phenomena within the cavities, the emphasis is put on the reproduction and monitoring of unsteady pressure fluctuations within the smooth cavities. Then, the LES of three configurations of real industrial turbines are conducted to study in situ the pressure fluctuations and apply the diagnostics already vetted on academic problems.

Les écoulements diphasiques liquide-gaz et liquide-vapeur sont présents dans de nombreuses problématiques industrielles. De fait, ils sont rencontrés dans des configurations diverses, tant en termes de dimensions que d’orientation et présentent des caractéristiques variées. Pour autant, il peut être distingué dans la littérature deux cas limites de modélisation de l’écoulement diphasique liées à son confinement : l’échelle micro et l’échelle macro. Cependant, la caractérisation de la transition entre ces deux échelles reste un enjeu majeur de compréhension de l’écoulement diphasique. De plus, les comportements thermiques en ébullition convective sont impactés par le phénomène de nucléation, fortement influencé par la géométrie de l’application et les conditions de saturations. De ce fait, l’objectif principal de la thèse est de contribuer à la compréhension des écoulements diphasiques adiabatiques et non adiabatiques à une échelle intermédiaire entre les échelles micro et macro, où les comportements tant dynamiques que thermiques sont mal définis. Pour ce faire, dans le travail présent, les effets du confinement et de l’inclinaison de l’écoulement diphasique sur les régimes d’écoulement, les chutes de pression et les échanges thermiques sont étudiés et comparés avec les observations et modèles de la littérature. Une section test inclinable de 1,6 mm de diamètre intérieur a été installée sur le banc d’essais de Charnay (2014), permettant l’étude de l’ébullition convective du R-245fa. Celle-ci permet de visualiser l’écoulement à l’intérieur de l’évaporateur ainsi que d’acquérir simultanément les chutes de pression et le champ de température de paroi extérieure. Pour ce faire, un tube en saphir chauffé par un dépôt d’ITO a été employé. Le champ de température a été obtenu par une procédure d’acquisition et de post-traitement infrarouge. Au préalable, une campagne d’essais a été menée en conditions adiabatiques sur la section de Charnay (2014). La section test présente, quant à elle, a permis deux campagnes d’essais en conditions adiabatiques et en ébullition convective. Des études menées en conditions adiabatiques est apparu un fort effet de l’inclinaison et du confinement sur les régimes d’écoulement ainsi que sur les chutes de pression. L’effet de l’inclinaison s’apparente à ce qui a pu être observé en macro-canal dans la littérature. De plus, l’inclinaison n’affecte que peu les échanges thermiques au regard des incertitudes inhérentes au dispositif infrarouge. Enfin, il est observé un effet de l’inclinaison sur le régime d’assèchement. Dans cette configuration, l’évolution temporelle du champ de température de paroi extérieure apparaît synchrone avec l’évolution dynamique de l’écoulement diphasique. / Liquid-gas and liquid-vapour two-phase flows are encountered in a large range of industrial applications. They are observed in multiple configurations, in terms of dimension and orientation and then have various characteristics. In the literature, it can be distinguished two limit cases of the two-phase flow linked to its confinement: micro and macro-scale. Nonetheless, characterisation of the micro-to-macro scale transition is still a main issue for the comprehension of two-phase flow. In addition, thermal behaviours of flow boiling are affected by the nucleation phenomenon, which is strongly influenced by the geometry of the application and the saturation conditions. The main objective of the thesis is to contribute to the comprehension of adiabatic and diabatic two-phase flows in mini-scale, where the general behaviours are not well definite. In the present study, the effects of confinement and orientation of two-phase flows on flow patterns, pressure drops and heat exchanges are studied and compared with observations and models of the literature. An 1.6 mm inner diameter inclinable test section was installed on the test bench of Charnay (2014), which enabled to study R-245fa flow boiling. This test section enables to visualize the flow directly in the evaporator and the simultaneous acquisition of the pressure drops and the outside wall temperature field. A sapphire tube, heated by a transparent ITO coating, is unemployed to insure the transparency. The temperature field is obtained by an IR image processing. Previously, tests were led on the section of Charnay (2014) in adiabatic conditions. The present test section was used for two series of tests, led in both adiabatic and diabatic conditions. It appeared a strong effect of confinement and orientation on the flow patterns and pressure drops. The effect of the orientation is likely the same that one observed in macro-scale. In addition, the orientation slightly affects heat exchanges in comparison with the uncertainties of the IR dispositive. Finally, it is observed an effect of orientation on dryout flow pattern. In this configuration, the temperature field evolution with time is synchronized with the dynamic evolution of the two-phase flow.

Le tissu adipeux est traditionnellement décrit comme étant constitué de lobules : des entités de formes ovoïdales composées de cellules et de vaisseaux et faiblement connectées entre elles.Récemment, il a été montré qu’un potentiel métabolique spécifique (le browning) co-localise avec cette organisation en lobules au sein d’un même tissu. Dans ce travail de thèse, nous nous intéressons à décrire plus précisément l’organisation structurelle et fonctionnelle du tissu adipeux selon plusieurs aspects. Dans un premier temps, on s’attache à segmenter les lobules du tissu adipeux en utilisant une méthode de traitement d’image originale. Nous mettons en évidence une organisation 3D complexe et suivant plusieurs échelles. En particulier, il semble que le potentiel de browning soit également lié à une organisation structurelle particulière en clusters de lobules. Dans un second temps, à partir d’imagerie 3D, nous reconstruisons le réseau vasculaire entier du tissu adipeux et réalisons une simulation d’écoulements sanguins micro-vasculaires. Plusieurs hétérogénéités structurelles et fonctionnelles sont alors mises en valeurs à l’aide d’une analyse en communautés qui composent le tissu adipeux (par algorithme de clustering). Ces résultats confirment l’existence d’une zone centrale fortement vascularisée et qui se démarque également comme étant le lieu d’une perfusion sanguine d’intensité différente. Dans une dernière partie, nous abordons la question de transferts thermiques entre vaisseaux sanguins suivant des géométries simples mais pertinentes. Nous réalisons une étude systématique des paramètres adimensionnels clés du problème et mettons en évidence un invariant des échanges de chaleur : un optimum à faible nombre de Péclet (convection de même ordre que la diffusion). Nous introduisons également une méthode de calibration de paramètres effectifs dans le contexte des modèles homogénéisés de température à travers des tissus vascularisés.