Ces travaux de recherche s'intéressent à la caractérisation des transferts thermiques dans les milieux fluide-particules, et en particulier, les lits fluidisés au sein desquels un solide divisé est mis en suspension par un fluide. La grande diversité d'échelles spatiales et temporelles dans ces procédés nécessite d'étudier les interactions hydrodynamiques, thermiques et/ou chimiques entre les particules et le fluide à l'aide d'une approche multi-échelles. Une étude des transferts thermiques dans des lits fixes puis fluidisés, est réalisée à deux échelles : locale (Particle Resolved Simulation) et moyennée (Discrete Element Method-Computional Fluids Dynamics). L'étude PRS permet de caractériser les couplages locaux des transferts thermiques entre particules ainsi que la dynamique de ces transferts dans les configurations fluidisées. Une étude comparative entre les échelles met en évidence les limites du modèle DEM-CFD à capter les fluctuations des transferts thermiques observées dans les simulations PRS. Dans un dernier temps, les fermetures du modèle DEM-CFD sont améliorées de manière à réintroduire les fluctuations perdues par le changement d'échelles. / This work aims at characterizing heat transfer into fluid-solid flows, and more particularly fluidized beds, into which a solid phase is suspended by a flowing fluid. The wide range of spatial and temporal scales present in such processes encourage to study hydrodynamic, thermal and/or chemical interactions between the particles and the fluid through a multi-scale strategy. The analysis of thermal interactions was first carried out for fixed bed configurations and then, fluidized beds at two overlapping scales: local (PRS; Particle Resolved Simulation) and mesoscopic (DEMCFD; Discrete Element Method-Computional Fluids Dynamics). The PRS approach accounts for the local coupling of heat transfer between the particles and its dynamics into fluidized beds. A comparative study of the two scales indicated the limits of the DEM-CFD model to capture the heat transfer fluctuations observed into PRS. In a last step, the closure laws for DEM-CFD were improved to reintroduce the fluctuations lost at this scale.

Le chapitre 1 porte sur une équation des ondes monodimensionnelle soumise à une condition aux limites unilatérale. Sous des hypothèses de régularité appropriées sur les données initiales, une nouvelle preuve d’existence et d’unicité est proposée. La méthode de redistribution de masse qui repose sur une redistribution de la masse d’un corps de telle sorte qu’il n’y ait pas d’inertie au niveau du nœud de contact est introduite et sa convergence est prouvée. Une approximation de ce problème d’évolution combinant la méthode des éléments finis ainsi que la méthode de redistribution de masse est analysée dans le chapitre 2. Puis deux problèmes ainsi que leurs solutions analytiques respectives (l’une étant nouvelle) sont présentés et des discrétisations possibles en utilisant différentes méthodes d’intégration en temps sont décrites. Enfin, des simulations numériques de ces problèmes sont reportées. Dans le chapitre 3, la masse des nœuds de contact est redistribuée sur les autres nœuds et sa convergence ainsi qu’une estimation de l’erreur en temps sont établies. Ensuite, une solution analytique déjà introduite dans le chapitre 3 est comparée aux approximations obtenues pour plusieurs redistributions de masse possibles mettant ainsi en évidence que plus une redistribution de masse d’un corps se fait à proximité des nœuds de contact meilleures sont les solutions approchées obtenues. Les problèmes de contact élastodynamique en dimension d’espace deux et trois sont étudiés dans le chapitre 4. Comme pour les problèmes de contact monodimensionnels, une solution approchée combinant les éléments finis et la redistribution de masse est exposée. Quelques simulations numériques utilisant des méthodes d’intégration en temps mettent en évidence les propriétés de convergence de la méthode de redistribution de masse. / The chapter 1 focuses on a one–dimensional wave equation being subjected to a unilateral boundary condition. Under appropriate regularity assumptions on the initial data, a new proof of existence and uniqueness results is proposed. The mass redistribution method based on a redistribution of the body mass such that there is no inertia at the contact node is introduced and its convergence is proved. An approximation of this evolutionary problem combining the finite element method as well as the mass redistribution method is analyzed in chapter 2. Then two benchmark problems (one being new) with their analytical solutions are presented and some possible discretizations using different time–integration schemes are described. Finally, numerical experiments for these benchmark problems are reported. In chapter 3, the mass of the contact nodes is redistributed on the other nodes and its convergence as well as an error estimate in time are established. Then an analytical solution already introduced in chapter 3 is compared to approximate ones obtained for different choices of mass redistribution highlighting that more a mass redistribution of the body is done near the contact nodes better the approximate solutions are obtained. The two and three–dimensional elastodynamic contact problems are studied in chapter 4. As for the one–dimensional contact problems, an approximated solution combining the finite element and mass redistribution methods is exhibited. Some numerical experiments using time–integration methods highlighted the convergence properties of the mass redistribution method.

Le réchauffement climatique est lié aux augmentations des concentrations de gaz à effet de serre dans l'atmosphère terrestre et en particulier aux émissions anthropiques de CO₂. La séquestration géologique a la capacité et la longévité potentielles pour diminuer de façon significative les émissions anthropiques de CO₂. Cette séquestration à grande profondeur induit des risques de fuite des réservoirs géologiques. Parmi les scénarios de fuite envisagés, celui d'une fuite diffuse est le plus inquiétant puisque sans surveillance, cette fuite pourrait perdurer et entrainer des séquelles sur l'environnement ainsi que des risques pour les populations. Des outils et protocoles de surveillance doivent donc être mis au point pour la surveillance en proche surface. Ce travail de thèse s'inscrit dans le cadre de cette problématique. Il a pour objectif la caractérisation, la quantification et la modélisation du transport et des interactions du CO₂ dans une zone non saturée carbonatée. Ce travail a suivi une approche expérimentale sur un site pilote naturel à Saint-Emilion (Gironde, France), avec la réalisation de fuites diffuses en ZNS carbonatée. Cette étude aborde plusieurs thématiques: la description et l'instrumentation du site pilote naturel ; la caractérisation physico-chimique de l'hétérogénéité du réservoir carbonaté ; l'étude du fonctionnement naturel de la ZNS carbonatée et en particulier la mise en place d'une ligne de base des concentrations en CO₂ ; la caractérisation de l'extension des panaches de gaz suite à des expériences de fuite diffuse dans la ZNS carbonatée et l'étude par simulation numérique des interactions gaz-eau-roche lors d'une fuite diffuse de CO₂ dans une ZNS carbonatée. Les résultats de ces travaux montrent l'importance de la caractérisation de l'hétérogénéité du réservoir carbonaté ainsi que des techniques d'échantillonnage et d'analyse des différentes phases en présence. L'établissement de la ligne de base a une importance particulière pour permettre de distinguer les variations naturelles de celles induites par une fuite diffuse de CO₂ dans la ZNS carbonatée. Les modes de transport du CO₂ vont évoluer en fonction des paramètres physico-chimiques. Ce transport se fait par advection et/ou par diffusion. L'utilisation de gaz inertes au niveau du site de séquestration géologique est très importante puisque la détection de ces traceurs permettrait de prédire les arrivées de panaches de CO₂ en proche surface. Par ailleurs, les interactions chimiques doivent être prises en compte dans les modèles de transport afin de pouvoir définir les facteurs de retard et l'impact d'une fuite diffuse de CO₂ sur une ZNS carbonatée. / Global warming is related to atmospheric greenhouse gas concentration increase and especially anthropogenic CO₂ emissions. Geologic sequestration has the potential capacity and the longevity to significantly diminish anthropogenic CO₂ emissions. This sequestration in deep geological formation induces leakage risks from the geological reservoir. Several leakage scenarios have been imagined. Since it could continue for a long period, inducing environmental issues and risks for human, the scenario of a diffusive leakage is the most worrying. Thus, monitoring tools and protocols are needed to set up a near-surface monitoring plan. The present thesis deals with this problematic. The aims are the characterisation, the quantification and the modelling of transport and interactions of CO₂ in a carbonate unsaturated zone. This was achieved following an experimental approach on a natural pilot site in Saint-Emilion (Gironde, France), where diffusive gas leakage experiments were set up in a carbonate unsaturated zone. Different aspects were investigated during the study: natural pilot site description and instrumentation; the physical and chemical characterisation of carbonate reservoir heterogeneity; the natural functioning of the carbonate unsaturated zone and especially the set-up of a CO₂ concentrations baseline; the characterisation of gas plume extension following induced diffusive leakage in the carbonate unsaturated zone and the study of gas-water-rock interactions during a CO₂ diffusive leakage in a carbonate unsaturated zone through numerical simulations. The results show the importance of the carbonate reservoir heterogeneity characterisation as well as the sampling and analysing methods for the different phases. The baseline set-up is of main interest since it allows discrimination between the induced and the natural CO₂ concentrations variations. The transfer of CO₂ in a carbonate unsaturated zone is varying in function of physical and chemical properties. This transfer is done by diffusion and/or advection. Because the detection of the noble gases allows the prediction of CO₂ plume arrival, the use of tracers in the sequestration site is of main importance. The chemical interactions have to be taken under account in transport models in order to predict delay factors and the impact of a CO₂ leakage in a carbonate unsaturated zone.

Dans un système de machines rotatives comme un moteur à turbine, les paliers d'éléments roulants à grande vitesse jouent un rôle important dans le support de l'arbre ou du rotor rotatif, et ont besoin d'une lubrification pour assurer leur fonction. Sauf qu'une petite quantité d'huile est nécessaire pour former le film lubrifiant élastohydrodynamique dans la zone de contact, la plus grande partie du lubrifiant reste en suspension dans l'air, formant un mélange huile/air. Ce phénomène entraîne des pertes hydrauliques parasitaires excessives lorsque les éléments roulants se translatent et tournent dans l'environnement fluide, ce qui peut constituer une partie relativement importante de la perte de puissance totale du roulement, appelée traînée d'enroulement et pertes de barattage. Pour une vitesse de rotation jusqu'à 3× 106 Ndm, la contribution de la traînée/dérive au total peut atteindre 50%. Cependant, jusqu'à présent, il existe peu d'approches utilisées directement pour l'estimation des pertes par traînage, qui ne pouvait fournir qu'une approximation plutôt grossière. Dans cette thèse, la méthode CFD est utilisée pour analyser d'abord l'écoulement autour d'un cylindre circulaire de longueur finie avec deux extrémités libres dans un espace ouvert. Ensuite, le modèle est remplacé par plusieurs cylindres circulaires en ligne pris en sandwich par deux parois plates, ce qui représente une approche simplifiée. Le fluide est ici considéré comme incompressible, représentant un fluide monophasé équivalent pour l'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec des propriétés de fluide spécifiées. Les résultats indiquent que l'écoulement autour de l'élément de rouleau de longueur finie est perturbé par ses deux extrémités libres, les anneaux environnants, la cage et d'autres éléments roulants. Il est proposé une relation entre le coefficient de traînée et le nombre de Reynolds approprié pour un cylindre circulaire dans les roulements à rouleaux (1<L/D<6), ainsi qu'une formulation pour la prévision des pertes de barattage. L'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec lubrification sous la course est également étudié dans ce travail. Le volume couplé de niveau de fluide (CLS-VOF) est utilisé pour démontrer la distribution du lubrifiant le long de la circonférence du palier. L'effet de divers facteurs est étudié, par ex. la vitesse d'injection d'huile, le diamètre de la buse, les propriétés de l'huile et l'angle d'injection de l'huile. La vitesse de rotation de tous les composants du palier est étudiée en particulier pour quantifier leur influence sur la fraction du volume d'huile à l'intérieur de la cavité du palier. Les résultats démontrent que non seulement la vitesse de rotation relative de l'anneau interne, mais la vitesse de la cage elle-même pourrait changer la distribution d'huile. / In a rotating machinery system like turbine engine, high speed rolling element bearings play an important role in supporting the rotating shaft or rotor, and need lubrication to insure their function. Except a small quantity of oil is needed to form the elastohydrodynamic lubricant film in the contact zone, most of lubricant remains in suspension in air, forming an oil/air mixture. This phenomenon leads to excessive parasitic hydraulic losses when rolling elements translate and rotate into the fluid environment, which may constitute a relatively large portion of the bearing's total power loss, named windage drag and churning losses. For high speed applications, i.e. for rotational speed up to 3× 10^6 Ndm, the contribution of drag/windage loss to the total one may reach up to 50%. However, so far there are few approaches used directly for drag and churning losses estimation, which could only provide a rather gross approximation. In this thesis, the Computational Fluid Dynamics (CFD) method is employed to analyze first the flow around one finite-length circular cylinder with two free ends in an open space. Then the model is changed to several in-line circular cylinders sandwiched by two flat walls, which represents a simplified approach. The fluid here is regarded as incompressible, representing an equivalent one-phase fluid for the oil/air two-phase flow inside the bearing cavity with specified fluid properties. The results indicate that the flow around the finite length roller element is perturbed by its two free ends, the surrounding rings, the cage and other rolling elements. A relationship between the drag coefficient and the Reynolds number suitable for circular cylinder in roller bearings (1<L/D<6) is proposed, as well as a formulation for churning losses prediction. The oil/air two phase flow inside the bearing cavity with under-race lubrication is also studied in this work. The coupled level-set volume of fluid (CLS-VOF) method is employed to demonstrate the lubricant distribution along the bearing circumference. The effect of various factors is studied, e.g. the oil injection velocity, the nozzle diameter, the oil properties, and the oil injecting angle. Rotational speed of all the bearing components are studied particularly to quantify their influence to the oil volume fraction inside the bearing cavity. The results demonstrate that not only the inner-ring relative rotational speed, but the cage speed itself could change the oil distribution. The results can be used for the precise lubrication design to optimate the oil distribution inside the bearing.

Dans cette thèse on s'intéresse aux jets en écoulement transverse dans une configuration générique de celle du refroidissement par effusion de chambres de combustion aéronautiques. L'amélioration des modèles de paroi avec transfert de masse passe par une meilleure connaissance de l'interaction entre les jets et l’écoulement principal. Nous avons donc réalisé la simulation numérique directe d'un jet issu d'un perçage incliné avec ou sans giration, pour des écoulements isothermes, turbulents et à bas nombre de Mach, dans un contexte compressible. Pour cela nous avons travaillé avec la bibliothèque AeroSol d'éléments finis continus et discontinus sur maillage hybride. En particulier nous nous sommes intéressés à la stabilité des flux numériques pour le compressible instationnaire associés à la méthode de Galerkin discontinue lorsque le nombre de Mach tend vers zéro. Nous avons pu mettre en évidence des comportements instables lors de l'utilisation de discrétisation temporelle explicite que nous avons corrigés en proposant un nouveau flux. Dans un deuxième temps, nous avons effectué les développements nécessaires à la réalisation des calculs. Nous nous sommes en particulier intéressés à la génération d'un champ de vitesse turbulent synthétique par la méthode SEM (Synthetic Eddy Method) que nous avons implantée dans AeroSol et validée. Grâce aux outils de post-traitement développés, nous avons conduit l'analyse de nos résultats. Dans le cas sans giration, les comparaisons avec les résultats expérimentaux et les résultats de simulations RANS que nous avons obtenus en parallèle sur la configuration du banc d'essai MAVERIC sont encourageants. La structure moyenne d'ensemble du jet est notamment correctement reproduite. En ce qui concerne la cas avec giration, le comportement attendu de déflexion successive du jet dans les deux plans caractéristiques (plan d'injection et plan de l'écoulement transverse) est bien reproduit et illustre tout le potentiel prévisionnel de la librairie de calcul que nous avons contribué à développer. / In this work we are interested in jet in crossflow in a generic configuration to the one used in effusion cooling for combustion chambers. Improved wall models with mass transfer requires a better knowledge of the interaction between the jets and the main flow. We therefore carried out the direct numerical simulation of a jet issuing from an inclined hole with or without gyration, for isothermal turbulent flow at low Mach number, in a compressible context. To achieved this, we worked with the continuous and discontinuous finite element library : AeroSol on hybrid grid. In particular we studied the stability of numerical flux for the unsteady compressible flow associated with discontinuous Galerkin method when the Mach number tends to zero. We were able to demonstrate unstable behavior when using explicit time discretization and we corrected them by providing a new flux. In a second time, we have performed the necessary development to achieve the calculations. We have been especially interested in the generation of a synthetic turbulent velocity field using the SEM method (Synthetic Eddy Method) that we have implemented in aerosol and validate. Thanks to the developed post-processing tools, we have conducted an analysis of our results. In the case without gyration, comparisons with experimental results and the results of RANS simulations we obtained on the Maveric test-bench configuration are encouraging. The mean flow of the jet is correctly reproduced. In the case with gyration, the expected behavior of successive deflection of the jet in both planes (injection plane and transverse plane of the flow) is reproduced and shows all the potential of the AeroSol library we helped to develop.

L'Accident Vasculaire Cérébral (AVC) - pathologie résultant d'une perturbation de l'apport sanguin dans le cerveau - est un problème de santé publique majeur, représentant la troisième cause de mortalité dans les pays industrialisés. Afin d'améliorer la prise en charge des patients atteints d'un AVC, il est important de posséder des méthodes efficaces pour l'identification des patients éligibles aux différentes thérapies et pour l'évaluation du rapport bénéfice/risque associé à ces thérapies. Dans ce contexte, l'Imagerie par Résonance Magnétique (IRM) dynamique de perfusion par contraste de susceptibilité, une modalité d'imagerie utile pour apprécier l'état de la perfusion cérébrale, peut aider à identifier les tissus à risque de s'infarcir. Cependant, l'intégralité de la chaîne de traitement, de l'acquisition à l'analyse et l'interprétation de l'IRM de perfusion demeure complexe et plusieurs limitations restent encore à surmonter. Durant ces travaux de thèse, nous contribuons à l'amélioration de la chaîne de traitement de l'IRM de perfusion, avec comme objectif final, l'obtention d'une meilleure prédiction de l'évolution de la lésion ischémique dans l'AVC. Dans une première partie, nous travaillons principalement sur l'étape de déconvolution des signaux temporels, une des étapes clefs à l'amélioration de l'IRM de perfusion. Cette étape consiste en la résolution d'un problème inverse mal-posé, et permet le calcul de paramètres hémodynamiques qui sont des biomarqueurs importants pour la classification de l'état final des tissus dans l'AVC. Afin de comparer de façon objective les performances des différents algorithmes de déconvolution existants et d'en valider des nouveaux, il est nécessaire d'avoir accès à une information sur la vérité terrain après déconvolution. Dans ce but, nous avons développé un simulateur numérique pour l'IRM de perfusion, avec une vérité terrain générée automatiquement. Ce simulateur est utilisé pour démontrer la faisabilité d'une automatisation du réglage des paramètres de régularisation, et établir la robustesse d'un algorithme de déconvolution avec régularisation spatio-temporelle d'introduction récente. Nous proposons également un nouvel algorithme de déconvolution globalement convergent. Enfin, la première partie de ces travaux se termine avec une discussion sur une autre étape de la chaîne de traitement en IRM de perfusion, à savoir, la normalisation des cartes de paramètres hémodynamiques extraites des images déconvoluées / Stroke – a neurological deficit resulting from blood supply perturbations in the brain – is a major public health issue, representing the third cause of death in industrialized countries. There is a need to improve the identification of patients eligible to the different therapies, as well as the evaluation of the benefit-risk ratio for the patients. In this context, perfusion Dynamic Susceptibility Contrast (DSC)-MRI, a prominent imaging modality for the assessment of cerebral perfusion, can help to identify the tissues at risk of infarction from the benign oligaemia. However, the entire pipeline from the acquisition to the analysis and interpretation of a DSC-MRI remains complex and some limitations are still to be overcome. During this PhD work, we contribute to improving the DSC-MRI processing pipeline with the ultimate objective of ameliorating the prediction of the ischemic lesion evolution in stroke. In a first part, we primarily work on the step of temporal signal deconvolution, one of the steps key to the improvement of DSC-MRI. This step consists in the resolution of an inverse ill-posed problem and allows the computation of hemodynamic parameters which are important biomarkers for tissue fate classification in stroke. In order to compare objectively the performances of existing deconvolution algorithms and to validate new ones, it is necessary to have access to information on the ground truth after deconvolution. To this end, we developed a numerical simulator of DSC MRI with automatically generated ground truth. This simulator is used to demonstrate the feasability of a full automation of regularization parameters tuning and to establish the robustness of a recent deconvolution algorithm with spatio-temporal regularization. We then propose a new globally convergent deconvolution algorithm. Then, this first part ends with a discussion on another processing step in the DSC-MRI pipeline, the normalisation of the hemodynamic parameters maps extracted from the deconvolved images. In a second part, we work on the prediction of the evolution of the tissue state from longitudinal MRI data. We first demonstrate the interest of modeling longitudinal MRI studies in stroke as a communication channel where information theory provides useful tools to identify the hemodynamic parameters maps carrying the highest predictive information, determine the spatial observation scales providing the optimal predictivity for tissue classification as well as estimate the impact of noise in prediction studies. We then demonstrate the interest of injecting shape descriptors of the ischemic lesion in acute stage in a linear regression model for the prediction of the final infarct volume. We finally propose a classifier of tissue fate based on local binary pattern for the encoding of the spatio-temporal evolution of the perfusion MRI signals

Parmi les options en cours d’investigation, le stockage souterrain de l'hydrogène dans les formations sédimentaires comme les grès pourrait offrir un potentiel unique pour stocker de grandes quantités d'énergie. L'évaluation des modalités de stockage souterrain de l'hydrogène nécessite donc à la fois une connaissance précise des transformations minéralogiques dues à la présence de l'hydrogène et l’acquisition de données sur le comportement hydrodynamique des fluides. Par conséquent, cette étude se composera de trois parties : 1- Etude des interactions géochimiques de l’hydrogène dans des formations sédimentaires gréseuses : Les produits expérimentaux portent la marque d'une réaction très limitée entre les minéraux du grès et l'hydrogène. Si les résultats expérimentaux sont combinés aux résultats numériques, l’étude démontre que l'hydrogène, une fois injecté, peut être considéré comme relativement inerte. De façon globale, nos résultats renforcent la faisabilité du confinement de l'hydrogène dans des réservoirs géologiques comme les grès. 2- Etude de la migration de l'hydrogène dans les grès : détermination de la perméabilité relative et de la pression capillaire du système hydrogène-eau : Afin de fournir des données quantitatives pour le développement du stockage souterrain de l'hydrogène, la pression capillaire et la perméabilité relative ont été mesurées pour le système hydrogène-eau en deux conditions potentielles. Les résultats indiquent que les données obtenues sont applicables à l’ensemble des conditions de stockage de l'hydrogène. 3- Modélisation numérique d’un site de stockage géologique d’hydrogène : La simulation numérique a été effectuée pour caractériser l'évolution dynamique d’un site de stockage d'hydrogène pur. Une fluctuation saisonnière du fonctionnement du réservoir et l'effet des fuites d'hydrogène dus aux réactions ont été pris en compte. / Underground hydrogen storage has been introduced as storage solution for renewable energy systems as it offers a unique potential to store large amounts of energy, especially in sedimentary formations such as sandstones. However, evaluating the underground hydrogen storage requires a precise knowledge of the hydrodynamic behavior of the fluids and of mineralogical transformations due to the presence of hydrogen that may affect the storage properties. Therefore, this study is consists in three parts: 1- Study of geochemical reactivity of hydrogen in sandstone sedimentary formations: The experimental products bear the mark of only very limited reaction between sandstone minerals and hydrogen. Taken together with the numerical results, this study demonstrates that hydrogen, once injected, can be considered as relatively inert. Overall, our results support the feasibility of hydrogen confinement in geological reservoirs such as sandstones. 2- Study of the migration of hydrogen in sandstone: determination of relative permeability and capillary pressure of hydrogen-water system: To provide quantitative data for the development of underground hydrogen storage, capillary pressures and relative permeabilities of hydrogen-water system have been measured at two potential conditions. The interpretation of the results would suggest that the obtained data are applicable for the entire range of hydrogen storage conditions. Interfacial tensions and contact angles for the hydrogen-water system have been also derived. 3- Numerical simulation of a geological hydrogen storage site: The numerical simulation was performed to characterize the evolution of pure hydrogen storage, by considering the seasonal fluctuation of renewable energy and the effect of hydrogen loses due to the biotic reactions.

Nous nous sommes intéressés dans la présente thèse à l’étude de l’écoulement gravitaire de nappesliquides non guidées qui s’écoulent verticalement dans l’air ambiant. Après une synthèse bibliographique,nous avons réalisé une double étude, expérimentale et numérique, en considérant des liquidesde viscosités différentes (allant de 1 à 50 fois celle de l’eau) et de tension superficielle proche de cellede l’eau. Le nombre de Reynolds Rel du liquide a varié de quelques unités à quelques milliers alorsque le nombre de Weber du liquide allait approximativement de 0.1 à 10. Le dispositif expérimentalque nous avons mis en place nous a permis de créer et d’étudier les formes géométriques des nappesliquides. Nous avons employé une méthode expérimentale originale pour mesurer le champ d’épaisseurdes nappes. Nous avons mené, parallèlement aux expériences, des simulations numériques 2D et 3Dinstationnaires et diphasiques (VOF), utilisant le calcul parallèle. Nous avons trouvé que les nappesexpérimentales et numériques deviennent plus courtes (verticalement) et moins épaisses lorsque le débitdiminue. Expérimentalement, lorsque le débit du liquide devient suffisamment faible, des filamentsliquides commencent à apparaître à coté d’une nappe moins large qu’auparavant. Dans le cas desnappes d’eau, cette transition de régime d’écoulement a été précédée de l’apparition systématique detrous dans la partie inférieure des nappes. Pour les autres liquides newtoniens plus visqueux (solutionsaqueuses de glycérine), l’apparition des filaments liquides a été précédée d’une déstabilisation des bourreletsqui délimitent la partie plane des nappes. Nous avons étendu par la suite l’étude expérimentaleà celle de fluides au comportement rhéologique plus complexe en utilisant un liquide non newtonienrhéofluidifiant. Nous avons montré pour ce fluide rhéofluidifiant que le débit associé à la transition versle régime des filaments diminue en comparaison avec celui associé à un liquide newtonien de viscositésimilaire. Cela nous a conduit à suggérer que la présence des propriétés rhéofluidifiantes des nappesliquides peut représenter une solution pour les applications de "coating" pour lesquelles on cherche àproduire des nappes stables et sans percement dans des configurations d’écoulement de faibles débits. / In this thesis, unguided plane liquid sheets flowing vertically by gravity in an ambient air atmosphereare studied experimentally and numerically. First of all a litterature survey clearly identified themain issues regarding the dynamics and modelling of such flows. Subsequently, different liquids exhibitinga wide range of viscosity (1 to 50 times that of water) and a surface tension close to that of waterwere selected. The liquid flow regimes were characterized by a Reynolds number Rel ranging from afew units to a few thousand while the Weber number Wel was varied between 0.1 to 10. A dedicatedexperimental system was designed and operated to study the relevant sheet features (geometry, thickness)by means of an original optical method. In parallel, finite volume based 2D and 3D simulationsof the flows were undertaken. All rely on the volume of fluid method (VOF) combined with adaptivemeshing. The experimental and numerical sheets became shorter (vertically) and thinner as the massflow rate decreased. Experimentally, when the mass flow rate of the liquid becomes sufficiently low,liquid threads begin to appear next to a narrower sheet than before. In the case of water, this flowregime transition was preceded by the systematic appearance of holes in the lower part of the sheets.For the other more viscous Newtonian liquids (mixtures of water and glycerin), the appearance of theliquid threads was preceded by a destabilization of the rims which delimited the flat part of the sheets.The experimental study was then extended to fluids featuring more complex rheological behavior e.g.by the use of a non-Newtonian shear-thinning fluid. For such a fluid, it was shown that the criticalmass flow rate associated with the transition towards the threads regime was lower than its Newtoniancounterpart of similar viscosity. It is suggested that the presence of shear-thinning properties in liquidsheets may represent a solution for "coating" applications for which stable and non-pierced curtainsin flow configurations of low mass flow rates are targeted.

Dans l'industrie aéronautique, les alliages de titane sont utilisés pour leur excellent comportement mécanique associé à une faible masse volumique. Ils sont largement employés sous forme de tôles dont la mise en forme peut se faire par le biais de trois procédés : à température ambiante par opération d'emboutissage, à très hautes températures (T≈900°C) par formage superplastique (SPF) et à des températures intermédiaires (T=730°C, 880°C) par formage à chaud (HF). Le projet repose sur le développement du procédé d'emboutissage à chaud d'une tôle d'alliage de titane Ti-6Al-4V en conditions isothermes à des températures inférieures à 700°C. Par conséquent, la détermination des paramètres procédés et matériaux constitue une étape importante pour la mise en œuvre de simulations numériques et contribue à la réussite des opérations d'emboutissage de pièces industrielles. Ces paramètres procédés sont liés à la vitesse du poinçon, aux efforts de serre-flan et au frottement induit entre le flan et l'outillage. Leur analyse a permis de déterminer deux niveaux de températures (400°C et 500°C) offrant une chute drastique du coût énergétique, en comparaison des procédés HF ou SPF, tout en conservant des niveaux d'allongement suffisants. Les paramètres matériaux influençant le comportement de l'alliage sont analysés et quantifiés. Ils peuvent être influencés par plusieurs mécanismes : élasticité, viscosité, anisotropie (Hill48, Barlat91) et nature de l’écrouissage (isotrope, cinématique). Dans cette étude, un modèle de comportement élasto-viscoplastique anisotrope, capable de considérer les trajets de chargement subis par la tôle lors de sa mise en forme, a été formulé pour les deux niveaux de température. L’implantation du modèle de comportement a été réalisée dans le code de calcul éléments finis Abaqus/Standard 6.14® interfacé avec le logiciel ZMAT®. Elle a permis d’une part des simulations d’emboutissage de profil Omega pour lesquelles des comparaisons avec les expériences ont été réalisées et d’autre part, des calculs sur une pièce de forme complexe. / In the aerospace industry, titanium alloys are used for their excellent mechanical behavior associated with low density. They are widely available in sheet form and the final shape can be obtained through three processes: at room temperature by stamping operation, at very high temperatures (T≈900°C) by superplastic forming (SPF) and at intermediate temperature (T=730°C, 880°C) by hot forming (HF). The project is based on the development of the hot stamping process of Ti-6Al-4V titanium alloy sheet under isothermal conditions at temperatures below than 700°C. Therefore, the determination of the process and material parameters constitutes an important stage for implementing the numerical simulation while contributing to the success of the stamping operation at the scale of an industrial part. The process parameters are related to the punch speed, the blank holder forces and the friction induced between the sheet and the tool. Their analysis allowed to determine two temperature levels (400°C et 500°C) leading a drastic drop in energy cost, compared to HF or SPF processes, while maintaining enough elongation levels. The material parameters influencing the behavior of the alloy are analyzed and quantified. They can be influenced by several mechanisms: elasticity, viscosity, anisotropy (Hill48, Barlat91) and nature of hardening (isotropic, kinematic). In this study, an anisotropic elasto-viscoplastic behavior model, able to consider the loading path undergone by sheet during forming, has been formulated for both temperature levels. The implementation of the behavior model is achieved in Abaqus/Standard 6.14® Finite Element code with the material library plugin ZMAT®. It enables, on the one hand, stamping numerical simulations of a simple shape Omega profile for which experimental comparisons were done, on the other hand, calculations on an industrial part with a complex shape.

La démarche expérimentale menée dans cette thèse CIFRE a eu comme objectif l’étude de la faisabilité de la séparation électrostatique d’une large gamme de mélanges de matériaux granulaires issus de déchets d’équipements électriques et électroniques (conducteur – isolant, conducteur – conducteur et isolant – isolant). Trois mécanismes physiques de charge électriques ont été mis en oeuvre : l’induction électrostatique, la décharge couronne et l’effet triboélectrique. Un dispositif de charge de laboratoire et deux séparateurs électrostatiques industriels ont été conçus et réalisés. Le dispositif de charge à cylindre(s) tournant(s) permet la charge triboélectrique de matériaux plastiques. Le prototype de séparateur industriel a d’abord été réalisé pour le traitement de mélanges granulaires conducteur – isolant puis a été adapté à la séparation isolant – isolant. Les trois mécanismes de charge peuvent y être utilisés. Ce prototype a permis la fabrication du premier séparateur électrostatique de série de la société CITF dédié à la séparation de déchets de câbles électriques. Des mesures de charge et de masse ainsi que des simulations et visualisations des trajectoires des particules par caméra rapide ont facilité l’analyse des phénomènes. La méthode des plans d’expérience a permis d’identifier les facteurs influents sur les processus de séparation et de définir leur point de fonctionnement optimal. L’influence d’une décharge à barrière diélectrique sur la charge triboélectrique et la séparation électrostatique de mélanges plastiques a été étudiée. Les travaux réalisés dans cette thèse ont confirmé la séparabilité électrostatique de mélanges de matériaux granulaires. / The experimental approach used in this thesis was aimed at the study of the feasibility of the electrostatic separation of a wide range of granular mixtures from waste electric and electronic equipment (conductive – insulator, conductive – conductive, insulator – insulator). Three physical mechanisms of electrical charging have been employed: electrostatic induction, corona discharge and triboelectric effect. A tribocharging device and two industrial electrostatic separators have been designed and engineered. The rotating cylinder tribocharging device allows the triboelectric charging of plastics. The industrial separator prototype was firstly built for processing conducting – insulator granular mixtures and then adapted for insulator – insulator separation. This prototype, which enables the use of all three electrical charging mechanisms, served as model for the CITF company first industrial electrostatic separator for the recycling of metals and plastics from electric cable wastes. Charge and mass measurement, as well as modelling and visualization of particle trajectories by a high-speed camera facilitated the investigations. Design of experiments method was used to define the influential factors and find the optimum operating conditions of the separation processes. Dielectric barrier discharge exposure of granular plastics may enhance the efficiency of the triboelectric charging of these materials and hence facilitate their electrostatic separation. These researches have confirmed the feasibility of the electrostatic separation of a wide variety of granular mixtures.