Ce travail de thèse a pour oblectif de simuler numériquemetn les impacts de recharges sédimentaires sur la morphodynamique du tronçon du Vieux-Rhin situé entre les barrages de Kembs et Breisach. Dans le cadre du projet InterReg de redynamisation du Vieux-Rhin, ces simulations sont primordiales pour évaluer l'efficacité des scénarios de recharge sédimentaire ainsi que lesrisques associés , tels que la fixation des dépôts ou encore la déstabilisation du lit pavé. a l'échelle du vieux-Rhin entioer (50 km) les simulations numériques 1D des années 1950, 1990 et 2009 ont permis d'appréhender la dynamique fluviale du lit en tenant compte de l'impact de la végétation. Les modélisations des processus granulométriques en jeu dans les recharges sédimentaires ont été préalablement testés et améliorés sur des expériences de laboratoires avec le modèles numériques 1D RubarBE. Dans ces modélisations, le tri granulométrique longitudinal, les évolutions morphologiques et le phénomène de capture des sédiments les plus fins ont été reproduits. Les modélisations bi-dimensionnelle à l'échelle des sites d'injections ont été réalisées avec le modèle numérique Rubar20TS. Les évolutions morphologiques des recharges ont été simulées et les zones favorisant les dépôts sédimentaires ont été identifiées. Les conclusions suivantes ont été tirées des simulations d'un banc sédimantaire et des scénarios de redynamisation : des flancs de banc pentus facilitent la mobilisation sédimentaire ; le choix du site et de la largeur du banc semblent très imporants afi de créer des conditions suffisament dynamiques en bordure de banc ; l'injection d'un banc de galets apparaît plus efficace qu'une érosion latérale de berge. Un scénario de recharge sédimentaire optimal pourrait consister en des injections de bancs latéraux étroits d'un volume d''environ 20 000 m3 répétées dans le temps et/ou dans l'espace.

Cette thèse s'intéresse à la turbulence incompressible homogène et anisotrope, et plus particulièrement à l'effet d'une rotation solide et d'un champ magnétique uniforme et stationnaire. En plus des simulations numériques directes (DNS), nous présentons également un modèle synthétique de turbulence (Kinematic Simulation ou KS) construit à partir d'une superposition de modes de Fourier, et au sein duquel la dynamique linéaire (basée sur la Rapid Distorsion Theory) peut être incluse. Dans un premier temps, l'effet de la rotation solide est étudié avec un effort particulier porté sur les corrélations en deux-temps. Une comparaison entre DNS et KS est proposée dans le cas isotrope comme dans le cas en rotation. Dans le contexte aéroacoustique, on montre dans quelle mesure la rotation modifie l'émission acoustique d'une turbulence homogène. L'effet d'un champ magnétique est ensuite considéré et l'anisotropie de l'écoulement est étudiée en fonction du nombre de Reynolds magnétique. Il est montré dans quelle mesure la turbulence magnétohydrodynamique quasi-statique est similaire à la turbulence "deux-dimensions trois-composantes" du fait de la dissipation Joule anisotrope et comment l'écoulement restaure son isotropie lorsque le nombre de Reynolds magnétique croît. Enfin, l'effet couplé d'une rotation et d'un champ magnétique est considéré. Les propriétés énergétiques ainsi que l'anisotropie sont étudiés et une étude paramétrique en fonction du nombre d'Elsasser est proposée.

Le phénomène de démouillage est caractérisé par la croissance d'un cristal sans contact avec la paroi du creuset due à l'existence d'un ménisque liquide au niveau de l'interface solide-liquide. Ceci crée un espace de quelques dizaines de micromètres entre le cristal et le creuset dans lequel il est élaboré et une des conséquences immédiates de ce phénomène est la nette amélioration de la qualité du cristal. Dans le cas du démouillage, la forme du ménisque est décrite par l'équation de Young-Laplace à partir de laquelle on a pu établir le système des équations différentielles non-linéaires qui décrivent l'évolution du rayon du cristal et de la hauteur du ménisque en fonction de diverses perturbations. L'intérêt réside dans la possibilité d'étudier numériquement, à partir de ces équations, les solutions stationnaires, leur stabilité, statique et dynamique, et d'effectuer une étude de la sensibilité des solutions vis-à-vis des paramètres du procédé et des conditions initiales.

L'objectif de ce travail est d'approfondir et de mieux comprendre les mécanismes responsables de la formation et de la croissance de la structure des grains dans le silicium multicristallin pour des applications photovoltaïques. Lors de la solidification du silicium multicristallin, la sélection des grains, le contrôle de la distribution de leur taille et leur direction de croissance sont des paramètres importants pour obtenir un matériau de bonne qualité et homogène. Ces paramètres influencent directement le rendement de conversion des cellules photovoltaïques, au travers de la capture et de la recombinaison des porteurs de charges et des interactions avec les impuretés. La structure de grains dans le silicium photovoltaïque évolue au cours de la solidification : des grains vont disparaître, d'autres vont apparaître, d'autres vont grossir pour donner au final une structure composée de gros grains, de petits grains dénommés 'grits', de joints de grains, et de macles. Il est donc important de comprendre les relations entre les différents paramètres du procédé industriel et leur influence sur les phénomènes physico-chimiques qui se produisent lors de la croissance afin de pouvoir influer sur la structure de grains dans le silicium, et de prévoir ses propriétés. Dans une première étape, nous avons établi un modèle de développement des grains basé sur le type de croissance (facettée, rugueuse ou mixte), la cinétique de ces divers types de croissances, le phénomène de maclage et la sélection des grains, dont nous montrons qu'ils sont, avec la germination initiale, à l'origine de la taille et de la structure des grains. Ensuite, nous proposons une approche de modélisation numérique de l'évolution de la structure des grains au cours de la solidification. Cette méthode se base sur l'analyse dynamique bidimensionnelle du joint de grains au niveau de la ligne triple grain-grain-liquide (rugueuse, facettée) tout en prenant en compte les phénomènes produits à l'échelle macroscopique (le champ de température local) et microscopique (la cinétique des grains). Le modèle résulte du couplage thermique et des mécanismes cinétiques de croissance. Nous avons donc développé un modèle numérique de croissance des grains en 2 dimensions et nous l'avons introduit dans le code 2D-MiMSiS qui se déroule en 2 étapes : Premièrement, le calcul en régime transitoire de la solidification macroscopique d'un lingot de silicium nous permet d'obtenir le champ thermique dans le lingot et la position précise de l'interface solide-liquide à différents instants ainsi que sa vitesse, son orientation (sa forme) et les gradients de température dans le liquide et le solide. Deuxièmement, la modélisation de la croissance est basée sur la description géométrique des joints de grains qui dépend de la cinétique des grains qui les bordent. Elle suit des critères dépendants de la morphologie (rugueuse ou facettée) de l'interface. Elle s'appuie sur le réseau d'isothermes du calcul thermique sans l'influencer dans un premier temps. Un des objectifs de ce modèle est de faire varier différents paramètres du procédé et d'en mesurer l'impact sur la structure cristalline finale. Des résultats de calculs 2D sont présentés et discutés par rapport à l'expérience.

La présente étude s'inscrit dans la thématique de la pollution particulaire de l'air intérieur et plus particulièrement dans le développement et la mise au point d'outils de prédiction du comportement de l'air chargé en particules à l'intérieur des bâtiments. La première partie présente un état de l'art des connaissances de la pollution particulaire. Les différentes particules rencontrées à l'intérieur des bâtiments sont recensées et classées selon leurs origines et leurs tailles. Le problème de la prédiction du mouvement des particules est introduit par l'étude de la diversité des forces élémentaires qui s'exercent sur les particules. La seconde partie est consacrée au développement d'un modèle de prédiction des concentrations en particules dans l'air à l'intérieur d'une enceinte monozone. Ce modèle considère l'air pollué comme une entité unique, et traite les échanges avec l'extérieur et les parois du bâtiment dans leurs globalités. Une expérimentation sur maquette a été mise au point afin d'étudier le dépôt des particules sur les revêtements du bâtiment. Cette étude a permis à la fois de quantifier ce dépôt et de mettre au point un principe de classification des revêtements vis-à-vis du dépôt particulaire. La dernière partie traite de la modélisation numérique du comportement des particules au sein d'écoulements confinés tels que ceux rencontrés dans le bâtiment. L'extension à la troisième dimension de l'espace du code de calcul développé au LEPTAB, a été validée par l'étude de la dispersion des particules au sein d'un écoulement turbulent de grille. Deux lois caractérisant l'interaction entre une particule et un revêtement ont pu être testées et comparées par la simulation de l'expérimentation sur maquette. Les perspectives à la fois sur le développement et la validation du modèle global ainsi que sur l'amélioration du code numérique de calcul de la dispersion des particules sont présentées et argumentées en conclusion.

La dynamique des flammes de prémélange est étudiée par deux approches numériques différentes. La première résout les équations compressibles de Navier-Stokes avec une chimie simplifiée (DNS). Afin de réduire les coûts de calcul, nous analysons et développons un schéma numérique à grille décalée. Le traitement des ondes acoustiques aux sorties est connu pour rendre les flammes cylindriques légèrement carrées. Ces déformations non-physiques sont expliquées en mettant en évidence la modélisation insuffisamment précise de l'accélération du fluide lorsque l'écoulement est oblique à la sortie. Une étude paramétrique et statistique de flammes turbulentes est menée en 2D et une simulation parallèle 3D est réalisée dans un domaine de (3cm)3. En considérant la flamme infiniment mince, l'approche EEM diminue considérablement les coûts de calcul. Les mêmes simulations sont réalisées et comparées aux résultats de DNS pour tester la capacité du modèle EEM à fournir des résultats quantitatifs.

Dans le cadre du projet INTERMARC (INTERaction dans les Moteurs à Allumage Commandé), la tâche du CORIA a consisté à produire une base de données à l'échelle RANS (provenant de données DNS) afin de tester, valider et modifier le modèle d'interaction développée par IFPen. Ce modèle vise l'ajout d'une composante d'interaction, phénomène non pris en compte par les lois de paroi actuelles.Ce projet repose sur l'interaction forte entre les différents protagonistes présents. Le CORIA et le CETHIL ont travaillé ensemble à la réalisation d'une base de données pour tester les modèles initiaux proposés par IFPen, puis en fonction des résultats obtenus, à itérer avec IFPen pour modifier et améliorer les modèles. Ces tests ont inclus des simulations 2D laminaires, 2D turbulentes, et 3D turbulentes.

Cette étude est consacrée à l'analyse des écoulements réactifs supersoniques cisailléset, plus particulièrement, des couches de mélange compressibles pouvant se développerdans les moteurs ramjet et scramjet. Des méthodes numériques appropriées ont été implémentéeset vérifiées pour aboutir au développement d'un code de calcul numériquemassivement parallèle, appelé CREAMS (compressible reactive multi-species solver). Cedernier a été spécialement conçu pour conduire des simulations numériques haute précision(simulations numériques directes ou DNS) de ce type d'écoulements. Une attentionparticulière a été portée à la description des termes de transport moléculaire et des termessources chimiques de façon à considérer la description physique la plus fidèle possible desmélanges des gaz réactifs à haute vitesse, au sein desquelles les temps caractéristiqueschimiques et de mélange aux petites échelles sont susceptibles d'être du même ordre degrandeur. Les simulations des couches de mélange bidimensionnelles et tridimensionnelles,inertes et réactives, confirment l'importance des effets associés à la compressibilité et autaux de dégagement de chaleur. Les résultats ainsi obtenus diffèrent en certains points deceux issus d'autres simulations qui introduisaient certaines hypothèses simplificatrices :développement temporel, emploi d'une chimie globale ou encore lois de transport simplifiées.En revanche, ils reproduisent certains tendances déjà observées dans un certainnombre d'études expérimentales conduites dans des conditions similaires.

Les procédés de soudage conduisent à des modifications de microstructure et induisent des contraintes et distorsions résiduelles jouant un rôle important sur la tenue mécanique des assemblages. Cette étude a pour objectif la modélisation des phénomènes thermomécaniques liés à l'opération de soudage du TA6V, ce qui implique d'identifier et de modéliser les sollicitations thermiques imposées par le soudage, le comportement métallurgique et thermomécanique du matériau. L'étude thermique et métallurgique a permis de montrer les principales évolutions de phases de l'alliage, durant le soudage laser. Ainsi, en première approximation, nous pouvons considérer que 3 phases distinctes sont à prendre en compte dans le calcul thermomécanique (la phase a, la phase ß et la phase a'). Afin de caractériser les phases ou mélange de phase en présence, des essais thermomécaniques sont réalisés pour diverses températures et conditions de déformation. Concernant le choix du modèle numérique, on adopte pour chacune des phases un modèle à variables internes de type Lemaître et Chaboche avec écrouissage cinématique non linéaire, un critère de plasticité de type Von Mises et une loi d'écoulement viscoplastique de type Norton, couplé à une loi de transition d'échelle. L'identification de chaque phase se fait par méthode inverse, à partir des essais thermomécaniques, grâce au module d'optimisation de ZSeT/ZMaT. La loi de localisation choisie est basée sur un modèle polycristallin en ß qui permet le passage du comportement microscopique de chaque phase au comportement macroscopique de l'ensemble. La validité de ce modèle décrit précédemment passera par la définition d'un essai expérimental représentatif des différents phénomènes rencontrés lors de l'opération de soudage.

Il existe deux stratégies principales pour fabriquer des engrenages, liées à l'ordre d'enchaînement de l'usinage de finition et du traitement thermique. Dans le cas qui nous occupe ici chez Renault, la majorité des pièces produites adopte un processus où le traitement thermique termine la gamme de fabrication, ce qui n'autorise aucune modification de la géométrie de la denture en sortie de trempe. Il est donc primordial d'analyser les déformations qui apparaissent lors de cette étape. Indépendamment de ces stratégies de fabrication, la vie d'une pièce peut être divisée en deux phases : la mise au point et la vie série. Chaque période génère des problématiques séparées. Pour la première, il est utile de connaître au préalable la déformation de la pièce lors du traitement thermique, voire des autres opérations de la gamme. Un état de l'art des modélisations phénoménologiques actuelles est donc dressé, ainsi qu'une liste des données nécessaires au calcul. L'applicabilité de la simulation numérique à ce problème est alors étudiée. Une bonne corrélation qualitative a été observée entre simulation et données expérimentales. Cependant, étant donné la complexité des donnée nécessaires et la détermination des conditions aux limites, une autre possibilité a également été développée, structurée autour d'une base de données recensant les déformations des engrenages de chaque site. Cette méthodologie a été implémentée dans une application qui est maintenant utilisée sur tous les sites de mécanique Renault. Une fois cette mise au point terminée, la pièce entre dans la deuxième phase, la vie série. Celle-ci voit apparaître alors un second type de problème : la dérive du processus nominal. Cette dérive peut être occasionnée par exemple par une matière légèrement modifiée, des conditions d'usinage optimisées ou une variation du traitement thermique qui échappe aux contrôles classiques mis en place. Nous nous attachons ici principalement aux dérives qui provoquent une variation dimensionnelle de la pièce mais le principe peut être appliqué à d'autres variations. Retrouver le plus rapidement possible l'origine d'une telle dérive et la corriger alors que la production est arrêtée est impératif. Nous proposons donc ici les fondements d'une méthode permettant de reconnaître l'origine d'une dérive identifiée à partir de la seule mesure de la déformation de la pièce. Une base des dérives possibles est créée en réalisant une décomposition modale (POD) des simulations numériques représentant les variations du process. En projetant la mesure réelle sur cette base, nous sommes ainsi capables d'identifier l'origine de la dérive et donc de réduire les délais d'analyse du problème.