Les défauts d'aspect des ouvrants automobiles sont des perturbations géométriques de la surface extérieure dont les dimensions caractéristiques ont une profondeur de l'ordre de quelques dizaines de micromètres et une étendue de quelques dizaines de millimètres. Ils entrainent des irrégularités de réflexion de la lumière sur la surface de la carrosserie et influent fortement sur la qualité perçue du véhicule. Le terme " ouvrants automobiles " regroupe les portes, le capot moteur et le coffre à bagages, pièces propices à l'apparition des défauts d'aspect à cause de leurs caractéristiques géométriques tourmentées. Cette étude s'appuie sur des essais expérimentaux pour caractériser les défauts d'aspect et sur la simulation par éléments finis pour les prédire numériquement. Les matériaux étudiés sont des aciers d'épaisseur comprise entre 0,6 et 0,8 mm. Leur carac- térisation mécanique est réalisée à partir d'essais de traction uniaxiale et de cisaillement simple (monotone et Bauschinger), avec pour objectif d'identifier les paramètres matériau des lois de comportement élastoplastiques utilisées dans la simulation numérique. Un dispositif de laboratoire a été conçu pour reproduire à petite échelle le défaut d'aspect apparaissant lors de l'opération de tombage d'un coin de porte. Sa modularité permet de faire varier différents paramètres comme la hauteur et le rayon de tombage. La géométrie des éprou- vettes tombées est mesurée avec une machine à mesurer tridimensionnelle. Une nouvelle méthode de détection et de caractérisation expérimentale est proposée à partir d'un calcul de courbure de sections. Cette méthode, moins subjective que la méthode traditionnelle par déméritage, permet de caractériser les défauts d'aspect au travers d'une cartographie des profondeurs sur des sections de la surface étudiée par des plans sécants. On montre que la simulation numérique reproduit l'allure du défaut d'aspect de coin de porte à condition d'utiliser un maillage fin, mais minimise sa profondeur, quel que soit le type d'élément fini utilisé. La méthode a été validée sur un coin de porte industriel présentant un défaut d'aspect lors de l'opération de tombage. Le défaut est observé qualitativement et mesuré par digitalisation optique. On montre que l'allure du défaut est reproduite par la simulation numérique sous réserve d'un maillage raffiné et de la prise en compte de toute l'histoire de déformation. Enfin, les défauts d'aspect apparaissant autour d'une poignée de porte ont aussi été caractérisés par digitalisation optique. La simulation numérique de ce défaut récurrent sur les automobiles, sous les mêmes conditions que pour le coin de porte, permet de prédire la majorité des défauts d'aspect observés.

.

La compréhension de l'état turbulent d'un plasma de fusion, responsable du faible temps de confinement observé, constitue un enjeu fondamental vers la production d'énergie par cette voie. Pour les machines les plus performantes, les tokamaks, les conductivités thermiques ionique et électronique mesurées sont du même ordre de grandeur. Les sources potentielles de la turbulence sont les forts gradients de température, de densité,... présents au coeur d'un plasma de tokamak. Si les pertes de chaleur par le canal ionique sont relativement bien comprises, l'origine du fort transport de chaleur électronique est quant à elle largement inconnue. En plus des fluctuations de vitesses électrostatiques, il existe des fluctuations de vitesses magnétiques, auxquelles des particules rapides sont particulièrement sensibles. Expérimentalement, le temps de confinement peut s'exprimer en fonction de paramètres non adimensionnels. Ces lois d'échelle sont encore trop imprécises, néanmoins de fortes dépendances en fonction du rapport de la pression cinétique à la pression magnétique, β et du rayon de Larmor normalisé, ρ* sont prédites.<br /><br />La thèse proposée ici cherche à déterminer la pertinence d'un modèle fluide non linéaire, électromagnétique, tridimensionnel, basé sur une instabilité particulière pour décrire les pertes de chaleur par le canal électronique et de déterminer les dépendances du transport turbulent associé en fonction de paramètres adimensionnels, dont β et ρ*. L'instabilité choisie est une instabilité d'échange générée par le gradient de température électronique (Electron Temperature Gradient (ETG) driven turbulence en anglais). Ce modèle non linéaire est construit à partir des équations de Braginskii. Le code de simulation développé est global au sens qu'un flux de chaleur entrant est imposé, laissant les gradients libres d'évoluer.<br /><br />A partir des simulations non linéaires, nous avons pu mettre en évidence trois caractéristiques principales pour le modèle ETG fluide: le transport de chaleur turbulente est essentiellement électrostatique; les fluctuations de potentiel et de pression forment des structures radialement allongées; le niveau de transport observé est beaucoup plus faible que celui mesuré expérimentalement.<br /><br />L'étude de la dépendance du transport de chaleur en fonction du rapport de la pression cinétique à la pression magnétique a montré un faible impact de ce paramètre mettant ainsi en défaut la loi empirique d'Ohkawa. En revanche, il a été montré sans ambiguïté le rôle important du rayon de Larmor électronique normalisé dans le tranport de chaleur: le temps de confinement est inversement proportionnel à ce paramètre. Enfin, la faible dépendance du transport de chaleur turbulent en fonction du cisaillement magnétique et de l'inverse du rapport d'aspect a été mise en évidence.<br /><br />Bien que le niveau de transport observé dans les simulations soit plus faible que celui mesuré expérimentalement, nous avons tenté une confrontation directe avec un choc de Tore Supra. Ce tokamak est particulièrement bien désigné pour étudier les pertes de chaleur électronique. En conservant la plupart des paramètres d'un choc bien référencé de Tore Supra, la simulation non linéaire obtenue donne un seuil en gradient de température proche de la valeur expérimentale. Le niveau de transport observé est plus faible d'un facteur cinquante environ que le transport mesuré. Un paramètre important qui n'a pu être conservé est le rayon de Larmor normalisé.<br /><br />La limitation en ρ* devra être franchie afin de confirmer ces résultats. Enfin une rigoureuse confrontation avec des simulations girocinétiques permettra de disqualifier ou non l'instabilité ETG pour rendre compte des pertes de chaleur observées.<br /><br />Mots-clés: fusion thermonucléaire, tokamak, plasma, turbulence ETG, simulations numériques

Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark « Ecoulement diphasique » proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé / The two-phase flow in porous media is a complex phenomenon and which relate to many industrial problems. EDF works on the feasibility and the safety of a storage in deep geologic layer of nuclear waste. In this domain the simulation of the two-phase flow in porous media is particularly important in at least three domains : first of all during the phase of ventilation of the galleries of the storage which could desaturate the rock and so modify its properties, but also during the phase of resaturation of the materials and finally during the arrival of the water on the metal parts contained in the storage which will then involve phenomena of corrosion and a hydrogen release. In this context, EDF wishes to obtain robust numerical methods without restrictive condition on the mesh. This work is dedicated at first to the development of the finite volume scheme SUSHI (Scheme Using Stabilization and Hybrid Interfaces) in the code of mechanics of EDF, Code_Aster in order to simulate the two-phase flow in porous media. This scheme was developed in 2D and in 3D. At the same time a new formulation which allows to simulate in a uniform way the flows in saturated and unsaturated porous media for miscible and immiscible problems is proposed. Various studies simulating difficulties related to the problems of the storage of nuclear waste in deep geological layers were study. We can quote the study of a bi-material which advances the capillary rebalancing of a material by an other one possessing properties and initial very heterogeneous conditions in saturation. We will also quote the study of the injection of hydrogen in an porous media initially saturated in pure water which is proposed by the benchmark "two-phase Flow " proposed by the GNR MOMAS. This study had for objective to bring to light the good treatment of the appearance of a phase in a saturated porous media and thus the relevance of our new formulation to study with a way unified a problem of saturated flow and a problem of unsaturated flow

Depuis quinze ans, de multiples instruments de mesure ont été envoyés vers Mars et explorent aujourd'hui sa surface ou y sont en orbite. Les données, renvoyées par ces sondes ont fourni une grande quantité d'informations quant aux conditions régnant aujourd'hui sur Mars, mais aussi sur la géologie des terrains et sur la minéralogie associée. Les minéraux découverts sur Mars se sont formés dans des conditions tout à fait différentes de celles qui y règnent aujourd'hui.Une approche associant modélisations numériques et expériences est proposée pour permettre d'identifier l'influence de divers processus comme le rayonnement UVC, la température et la pression partielle de CO2 ou de SO2. Des expériences en boîte à gants sont réalisées à des pressions partielles de CO2 de 10-6 bar, 10-3 bar et 10-2 bar, en présence ou non de rayonnement UVC. Conjointement à ces expériences, des expériences en capsules en or sont menées pour des températures de 20°C et 50°C avec des pressions partielles de CO2 d'environ 10-3 bar, 10-2 bar et 1 bar ainsi qu'avec une pression partielle de SO2 d'environ 1 bar. Les résultats ont montré qu'il était possible de reproduire expérimentalement et par les codes de calculs, la minéralogie supposée résultant de l'altération d'un verre basaltique de composition proche de roches martiennes. L'observation de carbonates de calcium dans les expériences réalisées à des pressions partielles de CO2 supérieures à 10-3 bar est en accord avec la prédiction de ces mêmes minéraux par des simulations numériques réalisées avec le logiciel DISSOL. Néanmoins, l'abondance de ces carbonates dans les deux types de modélisations, ainsi que leur composition à dominance calcique, contrastent avec les observations des rares minéraux carbonatés magnésiens à la surface de Mars. Ce contraste pose les questions concernant la composition de l'atmosphère, la composition initiale du verre et la préservation des carbonates après leur précipitation. Bien que des argiles n'aient pas été observées dans les expériences à cause d'un taux de dissolution du verre trop faible, des montmorillonites ferrifères et magnésiennes sont prédites dans les simulations numériques quelle que soit la pression partielle de CO2. De la beidellite sodique est également prédite dans le cas où les conditions atmosphériques martiennes seraient réductrices quelle que soit la pression partielle de CO2. Contrairement aux observations de la surface de Mars, aucune nontronite ne précipite. Les sulfates sont absents des paragenèses minérales obtenues lors des expériences avec une atmosphère à CO2. Dans les modélisations, aucun minéral sulfaté n'est néoformé pour les pressions partielles de CO2 imposées (> 10-6 bar) et la teneur en soufre du verre expérimental de synthèse. Toutefois, la précipitation de gypse peut être observée dans les simulations si la teneur en soufre est 100 fois supérieure à celle du verre utilisé pour ce travail. L'absence de sulfates magnésiens, de même que de carbonates magnésiens suggère que la composition chimique choisie pour le verre soumis à l'altération ne correspond pas totalement celle de la roche qui a pu donner les différents assemblages minéralogiques observés à la surface de Mars. En revanche, des sulfates de calcium et, des sulfates de fer et d'aluminium sont observés pour des expériences en capsules avec une atmosphère à SO2. Le rôle des atmosphères est essentiel quant à l'abondance des minéraux carbonatés et des minéraux sulfatés. Les expériences montrent que les radiations UVC accélèrent la dissolution du verre et favorise l'oxydation partielle du soufre du verre en pyrite / Process of deterioration on the primitive surface of Mars: experimental and digital simulations and mineralogical implications

Introduction: Le but de l’étude était d’examiner l’effet des matériaux à empreintes sur la précision et la fiabilité des modèles d’études numériques. Méthodes: Vingt-cinq paires de modèles en plâtre ont été choisies au hasard parmi les dossiers de la clinique d’orthodontie de l’Université de Montréal. Une empreinte en alginate (Kromopan 100), une empreinte en substitut d’alginate (Alginot), et une empreinte en PVS (Aquasil) ont été prises de chaque arcade pour tous les patients. Les empreintes ont été envoyées chez Orthobyte pour la coulée des modèles en plâtre et la numérisation des modèles numériques. Les analyses de Bolton 6 et 12, leurs mesures constituantes, le surplomb vertical (overbite), le surplomb horizontal (overjet) et la longueur d’arcade ont été utilisés pour comparaisons. Résultats : La corrélation entre mesures répétées était de bonne à excellente pour les modèles en plâtre et pour les modèles numériques. La tendance voulait que les mesures répétées sur les modèles en plâtre furent plus fiables. Il existait des différences statistiquement significatives pour l’analyse de Bolton 12, pour la longueur d’arcade mandibulaire, et pour le chevauchement mandibulaire, ce pour tous les matériaux à empreintes. La tendance observée fut que les mesures sur les modèles en plâtre étaient plus petites pour l’analyse de Bolton 12 mais plus grandes pour la longueur d’arcade et pour le chevauchement mandibulaire. Malgré les différences statistiquement significatives trouvées, ces différences n’avaient aucune signification clinique. Conclusions : La précision et la fiabilité du logiciel pour l’analyse complète des modèles numériques sont cliniquement acceptables quand on les compare avec les résultats de l’analyse traditionnelle sur modèles en plâtre. / Introduction: The purpose of this study was to examine the effects of different impression materials on the accuracy and reliability of digital models. Methods: Models from 25 patients selected at random from the files of the Department of Orthodontics at the University of Montreal were used in this study. One alginate (Kromopan 100 alginate, Italy), 1 alginate alternative (Alginot, Kerr Dentistry, Orange, CA), and 1 PVS (Aquasil, Dentsply Caulk, Milford, DE) impression was taken of both arches of each patient and sent to Orthobyte (Othobyte Digital Technology inc., Calgary, AB) for fabrication of a plaster model and scanning for production of a digital model. The Bolton 6 and 12 analyses and their constituent measurements, overbite, overjet, and arch length were used for the comparison. Results: The repeatability of measurements using both the plaster and digital methods was good to excellent, with the plaster measurements tending to be more reliable. There were statistically significant differences in the Bolton 12 and mandibular arch length and spacing measurements for all impression materials, with the plaster models tending to give a smaller measurement for the Bolton 12 and a higher measurement for the mandibular arch length and spacing. Although statistically significant differences in some measurements were found for the reliability and validity of digital models, none was clinically significant. Conclusions: The accuracy and reliability of the software for comprehensive cast analysis is clinically acceptable and reproducible when compared with traditional plaster study model analysis.

La dynamique instable des jets tournants est étudiée, en tenant compte des variations de masse volumique au sein de l'écoulement. Un code de simulation numérique directe permettant de résoudre les équations de Navier-Stokes à masse volumique variable a été développé, en utilisant une méthode originale et efficace pour résoudre le champs de pression. Analytiquement, deux modes instables bidimensionnels ont été mis en évidence, et sont identifiés comme des modes de Couette-Taylor et de Rayleigh-Taylor, ainsi qu'un troisième mode tridimensionnel, du à un couplage de vitesse. La dynamique instable de cet écoulement résulte d'une compétition entre ces trois modes, et les simulations numériques montrent que ces modes perdurent non linéairement. Ensuite, le comportement spatio-temporel de cette instabilité est étudiée par simulation numérique directe, et il a été montré qu'il existe une transition vers des modes absolument instables, sous l'effet du rapport de densité s ainsi que du taux de rotation q. Cette dynamique est également étudiée expérimentalement au travers de plusieurs méthodes de mesures, et la présence de mode globaux auto-entretenus est mise en évidence qui sont en bon accord avec les résultats numériques. Finalement, le phénomène de l'éclatement tourbillonnaire est étudié, et montre le rôle prépondérant de la viscosité réelle. En effet, l'éclatement tourbillonnaire est un mécanisme permettant de soulager le système de l'intensification de la vorticité, au travers de la viscosité, alors qu'il n'apparaît pas en traitant les équations d'Euler tronquées. / The unstable dynamics of a swirling jet flow is studied, including density variations within the flow. A direct numerical simulation method was developed to solve variable density Navier-Stokes equations, using an accurate and efficient pressure solver. Analitically, two unstable bi-dimensionnal modes are highlighted, and are identified as Couette-Taylor and Rayleigh-Taylor modes. A three-dimensionnal mode is also highlighted, wich is created by the shear. Numerical simulations show that those modes are nonlinearly persistant. Then, the spatio-temporal instability behaviour is studied numerically, and show that the instability undergoes to a convective/absolute transition with density ratio s and rotation rate q. This dynamic is also studied experiementally through different methods, and Global selfsustained modes are highlighted wich are in ggod agreement with numerical results. Finally, the vortex breakdown phenomenon is studied, and show the crucial role of real viscosity. Indeed, the vorticity intensification is relaxed through the viscosity effect, while it is not treating the truncated Euler Equations.

Le travail de cette thèse a porté sur le développement de techniques numériques pour l'homogénéisation de matériaux présentant à une petite échelle des hétérogénéités aléatoires. Sous certaines hypothèses, la théorie mathématique de l'homogénéisation stochastique permet d'expliciter les propriétés effectives de tels matériaux. Néanmoins, en pratique, la détermination de ces propriétés demeure difficile. En effet, celle-ci requiert la résolution d'équations aux dérivées partielles stochastiques posées sur l'espace tout entier. Dans cette thèse, cette difficulté est abordée de deux manières différentes. Les méthodes classiques d'approximation conduisent à approcher les propriétés effectives par des quantités aléatoires. Réduire la variance de ces quantités est l'objectif des travaux de la Partie I. On montre ainsi comment adapter au cadre de l'homogénéisation stochastique une technique de réduction de variance déjà éprouvée dans d'autres domaines. Les travaux de la Partie II s'intéressent à des cas pour lesquels le matériau d'intérêt est considéré comme une petite perturbation aléatoire d'un matériau de référence. On montre alors numériquement et théoriquement que cette simplification de la modélisation permet effectivement une réduction très importante du coût calcul / In this thesis we investigate numerical methods for the homogenization of materials the structures of which, at fine scales, are characterized by random heterogenities. Under appropriate hypotheses, the effective properties of such materials are given by closed formulas. However, in practice the computation of these properties is a difficult task because it involves solving partial differential equations with stochastic coefficients that are additionally posed on the whole space. In this work, we address this difficulty in two different ways. The standard discretization techniques lead to random approximate effective properties. In Part I, we aim at reducing their variance, using a well-known variance reduction technique that has already been used successfully in other domains. The works of Part II focus on the case when the material can be seen as a small random perturbation of a periodic material. We then show both numerically and theoretically that, in this case, computing the effective properties is much less costly than in the general case

La contribution des processus de fusion de galaxies lors de l'assemblage cosmologique de masses est encore mal comprise. Dans ce contexte, l'étude du support dynamique des galaxies est un moyen permettant de contraindre les différents scénarios d'évolution.Dans une première partie, je présente ma contribution à l'analyse cinématique des galaxies de l'échantillon MASSIV. Je définis une technique d'analyse qui vise à quantifier la fraction d'objets susceptibles d'être le produit d'une fusion.J'introduis dans une seconde partie une méthode de construction de conditions initiales de galaxies idéalisées.Dans une troisième partie, je présente la définition de l'échantillon de simulations MIRAGE, reproduisant les propriétés physiques des galaxies jeunes. Les simulations montrent en particulier l'absence de flambées de formation stellaire durant les fusions ce qui suggère un mécanisme de saturation.Une analyse des propriétés des complexes granulaires issus d'instabilités gravitationnelles est réalisée dans la quatrième partie, soulignant la compatibilité du scénario granules âgés avec les observations en terme de taux d'éjection de gaz et d'âge de populations stellaires. Je mets en évidence que le mécanisme de fusion de granules permet de contrôler la croissance de la masse du bulbe, d'aplatir le profil central de densité du halo de matière sombre ainsi que d'éjecter de grandes quantités de gaz.Dans une cinquième partie, je complète ce travail en réalisant une étude comparative des données cinématiques MASSIV à plus de 4000 pseudo-observations construites à partir des simulations de l'échantillon MIRAGE afin de déterminer la capacité à détecter des signatures de fusion. / The contribution of the fusion processes to the cosmological mass assembly is still poorly understood. In this context, the study of the dynamical support of galaxies is a way to constrain different evolution scenarios.In the first part, I introduce my contribution to the kinematical analysis of galaxies MASSIV. I define an analysis that aims at quantifying the fraction of objects that are likely to be the product of a recent merger.I introduce in the second part a method to build initial conditions of idealized galaxies.In the third part, I present the definition of the simulation sample MIRAGE, designed to reproduce the physical properties of young galaxies. The simulations show in particular the absence of a star formation bursts in mergers of fragmented and turbulent disks, suggesting a saturation mechanism.An analysis of the properties of clumps that are generated by gravitational instabilities is performed in the fourth part of this thesis, highlighting the compatibility of long-lived clumps scenarios with observations in terms of outflowing material and stellar population ages. I also highlight that the gas rich clump fusion mechanism is able to control the bulge mass growth, to erode the central profile of the dark matter halo and to drive massive gas outflows in the disk plane.Finally, in the fifth part, I complete this work by carrying out a comparative study of the MASSIV kinematical data to a set of more than 4000 pseudo-observations at z=1.7 built from simulations of the MIRAGE sample to determine the ability to detect galaxy merger signatures.

Cette étude s'inscrit dans le domaine des moteurs à doublets d'hélices contrarotatives, aussi appelés open rotors. Elle a porté sur le développement des méthodologies à coût modéré pour la simulation aérodynamique des open rotors. Elle avait pour objectif, d'un côté, la mise en place et la validation de ces méthodologies rapides ; et d'un autre côté, l'approfondissement de la compréhension de l'origine des efforts dans le plan hélice, aussi appelés les efforts 1P. Pour le premier des objectifs, le code HOST-MINT, basé sur la méthode de la ligne portante a été adapté et amélioré pour la simulation de l'aérodynamique instationnaire des hélices et des open rotors. Des validations ont été réalisées avec succès par comparaison avec des données expérimentales et des simulations CFD plus avancées. Cette étude a ainsi ouvert des nombreuses perspectives d'application de ce type de méthodologies rapides dans la conception aérodynamique des futurs open rotors. En outre, cette méthode s'est révélée aussi adaptée pour d'autres domaines autour de l'aérodynamique, comme par exemple, pour les problèmes d'aéroélasticité ou pour les prédictions aéroacoustiques préliminaires. / This study is devoted to the development of moderate-cost methodologies for the aerodynamic simulation of open rotors. The main goals are, on one side, to develop and validate these rapid methodologies, and, on the other side, to better understand the mechanisms behind propeller in-plane loads, also called the 1P loads. To reach the first goal, the HOST-MINT code, based on the lifting-line theory, has been adapted and improved for the unsteady simulation of propellers and open rotors. The code has been assessed by comparison with experimental data and more complex and precise CFD simulations. Finally, the first developments and tests of a Lagrangian/Eulerian coupling strategy between HOST-MINT and the elsA CFD code have been performed. These studies enable to endeavor a number of applications of this type of rapid methodologies in the aerodynamic design of future open rotors. Moreover, these methodologies may be adapted for other domains linked to aerodynamics, such as aeroelastic problems or preliminary aeroacoustic predictions.