L’analyse de sûreté des réacteurs nucléaires nécessite la modélisation fine des phénomènes y survenant et plus spécifiquement ceux permettant d’assurer l’intégrité des barrières de confinement. Les outils de modélisation et codes actuels favorisent une analyse fine du système réacteur par discipline dédiée, et couplée avec des modèles simplifiés. Néanmoins, le développement depuis plusieurs années d’une approche dite « Best Estimate », basée sur des calculs multiphysiques et multi-échelle, est en cours de réalisation. Cette approche permettra d’accéder au suivi et à l’analyse détaillée de problèmes complexes tels que l’étude des Réacteurs nucléaires en situation standard et accidentelle. Dans cette approche, les phénomènes physiques sont simulés aussi précisément que possible (selon la connaissance actuelle) par les modèles couplés. Par exemple, des codes disciplinaires existent et permettent la modélisation précise de la neutronique, de la thermohydraulique du cœur du réacteur ou de la thermohydraulique sur l'ensemble du système, de la thermomécanique du combustible ou des structures. Une approche « Best Estimate » consiste à coupler ces modèles afin de réaliser une modélisation globale et précise du système de réacteur nucléaire. Cette approche nécessite de bien définir les modèles qui sont utilisés afin de préciser exactement leurs limites, et donc préciser les incertitudes des résultats des modèles couplés afin de les assumer et de les optimiser.C’est dans ce contexte de travail que s’inscrit cette thèse. Elle consiste dans le développement d'un couplage multiphysique et multi-échelle « Best Estimate » afin d'obtenir une analyse précise des Réacteurs à Eau Légère en situations normale et accidentelle. Elle a consisté principalement en l’analyse des modèles et de leurs interactions et à la mise en œuvre d'un algorithme de couplage multiphysique entre une neutronique et une thermohydraulique exprimées à l'échelle du réacteur, ainsi qu’avec une thermomécanique fine à l'échelle élémentaire du crayon combustible. En outre, un travail spécifique a été effectué afin de préparer ou d'améliorer l’accés à l'information physique locale nécessaire à la mise en œuvre de modélisations couplées multi-échelles, à l'échelle du combustible. / The safety analysis of nuclear power plants requires a deep understanding of underlying key physical phenomena that determine the integrity of the physical containment barriers. At the present time, cutting edge models focus on a single aspect (discipline) of the physical system coupled with rough models of the other aspects needed to simulate the global system. But, safety analyses can be carried out based on Multiphysics and Multiscales modelling. This Best Effort approach would give a full and accurate (High Fidelity) comprehension of the reactor core under standard and accidental situations. In this approach, the physical phenomena are simulated as accurately as possible (according to present knowledge) by coupled models in the most efficient way. For example, codes exists that are accurate modellings of Neutronics, or modellings of thermal fluid mechanics inside the core, or modellings of thermal fluid mechanics over the whole system, or modellings of thermal mechanics of the fuel pin or over the whole device structure. A Best Estimate approach would couple these models in order to realize a global and accurate modelling of the Nuclear reactor. This approach requires to define well the models that are used in order to exactly specify their limits, and hence, specify uncertainties of the coupled model results in order to assume and optimize them.It is in this context that this PhD thesis work is being under taken. It consists in the development of a Multi-physics and multi-scale Best Estimate modelling in order to obtain an accurate analysis of Pressurized Water Reactor under standard and accidental operating situations. It mainly involves the understanding of each model and their interactions, followed by the implementation of multiphysics algorithms coupling Neutronics and Thermohydraulics at reactor scale to an accurate Thermomechanics at the elementary scale of the fuel pin. In addition, a work project has been carried out in order to prepare or improve the access to the local physical informations that are needed for the implementation of multiscale coupling scheme, at the elementary scale of the fuel pin.

En automobile, les suspensions semi-actives à base de fluides magnétorhéologiques (MR) proposent un compromis très intéressant entre des performances élevées et une Consommation électrique faible. Par le contrôle et l'ajustement de leurs paramètres d'amortissement, elles permettent de réaliser une optimisation en temps réel du comportement sur une large bande de fréquences de sollicitations. Cette technologie utilise les fluides MR dont la viscosité augmente sous l'action d'un champ magnétique, en quelques ms. Ce mémoire présente la conception et le dimensionnement d'un amortisseur MR ainsi que le banc de test associé. La caractérisation dynamique et statique de l'amortisseur a été effectuée. L'évaluation expérimentale de la force d'amortissement en fonction des grandeurs dynamiques et électriques a été mesurée et comparée au modèle théorique utilisé dans la phase de conception. Le comportement en régime statique de l'amortisseur a été modélisé par un modèle de Bingham et le comportement dynamique par un modèle de Bouc-Wen. La dernière partie de l'étude concerne la mise en place et la validation de plusieurs lois de contrôle pour les suspensions semi-actives, notamment de nouvelles lois basées sur l'inversion du modèle REM. Une comparaison des performances en termes de confort et de consommation électrique de chaque loi est présentée. Ces nouvelles lois de contrôle ont montré de très bonnes performances tout en réduisant la consommation électrique par rapport aux autres lois déjà existantes.

Le principal défi auquel sont confrontés aujourd'hui les équipementiers aéronautiques est d'augmenter l'utilisation des systèmes électriques à bord de l'avion. De nos jours, le remplacement des systèmes hydrauliques par des actionneurs électriques conduit à placer les systèmes électriques dans un environnement hostile, par exemple dans la nacelle du moteur. L'équipement est soumis à des contraintes sévères telles que des températures élevées et basses, un cycle thermique étendu, une humidité élevée et une basse pression. En conséquence, des efforts doivent être faits pour réduire le poids et le volume du convertisseur de puissance sans perdre ses performances. Pour atteindre cet objectif, la conception de modules de puissance doit permettre un haut niveau d'intégration, d'efficacité et de fiabilité. On s’intéresse en particulier aux dommages causés par la fatigue qui ont une influence significative sur les performances électriques de ces modules. Les tests de performance liés à la fatigue restent des efforts coûteux pour l'équipement aéronautique. Un nombre fini de tests destructifs, par vieillissement accéléré, peut être effectué pour un nombre assez faible de configurations. Le but de ces tests est d'étudier les modes de défaillance apparaissant lors du vieillissement accéléré. Par conséquent, des simulations numériques ont été envisagées, facilement évolutives et utilisables pour un grand nombre de configurations, mais nécessitant des données d'essais expérimentaux. Dans ce manuscrit, quelques modes de défaillances sont étudiés. On propose une méthode numérique intégrant les contraintes principales dans les équipements, à savoir la simulation électrique, thermique et mécanique. Ces trois problèmes physiques ont des temps caractéristiques différents et sont fortement couplés avec un comportement non trivial. Pour optimiser l'utilisation des ressources et avoir une représentation pertinente du problème, un procédé couplé électrique 1D / thermique 3D / mécanique 3D a été implémenté sur un bus de cosimulation. Différents pas de temps, différents niveaux d'abstraction et différentes compétences sont utilisés pour fournir un modèle multiphysique de modules de puissance. / Today’s main challenge for aeronautical equipment manufacturers is to respond to the more electrical aircraft regulations. Moreover, there are many applications in aircraft area where high temperature technologies are needed. Nowadays, the replacement of hydraulic systems for electric ones leads to place the power inverters in a harsh environment, for example in the engine nacelle. The equipment is under high constraints such as high and low temperatures, wide temperature cycling, high humidity and low pressure. Combined to these environmental constraints, the new aircraft system is submitted to weight and operating cost reduction. As a consequence, efforts shall be done to reduce weight and volume of the power converter without losing its performance. To reach such a goal, the design of the converter must enable a high level of integration, efficiency and reliability. In particular, fatigue damage has a significant influence on such modules electrical power performance. And fatigue-related performance testing remains a costly endeavor for aeronautical equipment. A finite number of destructive tests can be carried out in specific facilities for a fairly low number of configurations. The purpose of these destructive tests is to investigate the failure modes appearing regarding this accelerated ageing. Therefore numerical simulations have been envisaged since non-destructive, easily evolving and usable for a high number of configurations, though needing data from experimental assays. In this study, we propose a method dealing with the main constraints for such equipment, i.e. electrical, thermal and mechanical simulation. Those three physical problems have different characteristic time and are strongly coupled with a non-trivial behavior. To optimize the resources usage and have a relevant representation of the problem, a 1D electrical / 3D thermal / 3D mechanical coupled method has been implemented over a co-simulation bus. Different time steps, different abstraction levels and different skills are used to provide predictions of the multiphysical fatigue behavior of power modules.

Les batteries lithium-ion s'affichent comme de bons candidats pour assurer le stockage réversible de l'énergie électrique sous forme électrochimique. Toutefois, elles sont à l'origine d'un certain nombre d'incidents aux conséquences plus ou moins dramatiques. Ces incidents sont souvent liés au phénomène d'emballement thermique. La sécurité des batteries Li-ion représente par conséquent un enjeu technique et sociétal très important. C'est dans ce contexte que vient s'inscrire ce travail de thèse dans le cadre d'une collaboration entre IFPEN, l'INERIS et le LISE. Une double approche de modélisation et expérimentation a été retenue. Un modèle 3D du comportement thermique a été développé à l'échelle de la cellule, couplant les phénomènes thermiques et chimiques, et prenant en compte le vieillissement par croissance de la SEI sur l'électrode négative. Le modèle a été calibré pour la chimie LFP/C sur deux technologies A123s (2,3 Ah) et LifeBatt (15 Ah), puis validé expérimentalement. Le modèle permet d'identifier les paramètres critiques d'emballement de cellules, il permet également de discuter l'effet du vieillissement sur l'emballement thermique. Grâce à l'expérimentation, les connaissances en termes d'amorçage et de déroulement d'un emballement thermique d'une batterie Li-ion, ont pu être enrichies, en particulier pour les cellules commerciales LFP/C cylindriques A123s, LifeBatt, et pour les cellules NMC/C prismatiques en sachet souple PurePower (30 Ah). Cette étude ouvre de nouvelles possibilités pour améliorer la prédiction des différents événements qui ont lieu lors de l'emballement thermique des batteries Li-ion, à différentes échelles. / Li-ion secondary batteries are currently the preferred solution to store energy since a decade for stationary applications or electrical traction. However, because of their safety issues, Li-ion batteries are still considered as a critical part. Thermal runaway has been identified as a major concern with Li-ion battery safety. In this context, IFPEN, INERIS and LISE launched a collaboration to promote a PhD thesis so called « understanding and modeling of thermal runaway events pertaining to new and aged Li-ion batteries ». To achieve this goal, a double approach with modeling and experimental investigation is used. A 3D thermal runaway model is developed at cell level, coupling thermal and chemical phenomena, and taking into account the growth of the SEI layer as main ageing mechanism on negative electrode. Advanced knowledge of cells thermal behavior in over-heated conditions is obtained particularly for commercial LFP / C cylindrical cells: A123s (2,3Ah), LifeBatt (15Ah), and NMC / C pouch cells: PurePower (30 Ah). The model was calibrated for LFP / C cells, and then it was validated with thermal abuse tests on A123s and LifeBatt cells. This model is helpful to study the influence of cell geometry, external conditions, and even ageing on the thermal runaway initiation and propagation. This study opens up new possibilities for improving the prediction of various events taking place during Li-ion batteries thermal runaway, at various scales for further practical applications for safety management of LIBs.

La faible exploitation de l'irradiation inter-rangée limite la production des modules photovoltaïques (PV). Le projet "Aleph" explore l'intérêt d'ajouter des réflecteurs plans entre les rangées pour augmenter la production, et dégage des règles claires permettant l'optimisation géométrique de l'ensemble. Ce travail présente une modélisation multiphysique du système, des simulations numériques de son comportement, et la comparaison avec des données expérimentales. Deux technologies de module PV sont considérées : silicium amorphe (a-Si) et silicium polycristallin (p-Si). Les mesures montrent des gains énergétiques importants grâce aux réflecteurs. Les gains sont plus importants pour les modules a-Si que p-Si. La modélisation associe un modèle optique de lancers de rayons par méthode Monté-Carlo sous EDStaR, un modèle photoélectrique sous SPICE, et un modèle thermique empirique. Le modèle complet est calibré avec des données expérimentales en utilisant un algorithme évolutif. Une fois calibré, le modèle démontre une bonne performance en simulant la puissance générée par les modules en fonction des données atmosphériques et radiatives. / The poor utilisation of the inter-row irradiation limits the production of photovoltaic (PV) modules. The "Aleph" project explores the potential of adding inter-row planar reflectors to increase the system yield, and defines clear rules for optimal settings of such systems in a given location and under a given climate. This work presents a multiphysics model of the system, numerical simulations of its behaviour, and the comparison with experimental data. Two PV module technologies are tested: amorphous silicon (a-Si) and polycrystalline silicon (p-Si). The experimental data show significant gains in produced energy brought by the reflectors. The gains are higher for a-Si modules compared to p-Si. The modelling work combines a Monte-Carlo ray-tracing optical model (EDStaR), a photo-electric model (SPICE), and an empirical thermal model. The complete model is calibrated with measurements using an evolutionary algorithm. Once calibrated, the model demonstrates good performance in predicting the module power output as a function of atmospheric and irradiance data.

Le Haut-Parleur Digital Matriciel (HPDM) est un moyen de transduction électroacoustique qui reçoit comme entrée un signal numérique et qui effectue la conversion vers l'analogique directement dans l'air. Il est constitué de plusieurs éléments rayonnants disposés au sein d'une matrice. Ces éléments seront désignés par le terme «speaklet» lorsqu'ils sont de tailles réduites. Le rayonnement acoustique du HPDM est en effet très sensible à la taille de la matrice, ce qui le rend tout particulièrement adapté à la technologie MEMS. Cette thèse porte sur l'étude de HPDM MEMS piézoélectriques. Après une introduction qui débute par certaines généralités jusqu'à se focaliser sur le sujet, la thèse aborde les modélisations multiphysiques des HPDM étudiés, le dimensionnement des speaklets puis les tests expérimentaux. Des modèles analytiques ainsi que des simulations numériques et par éléments finis sont mis en place et permettent de prédire le comportement mécanique des speaklets présentés, les pressions rayonnées par les HPDM et les puissances électriques consommées. Les speaklets sont ensuite dimensionnés à partir de l'empilement technologique afin de maximiser le niveau de pression qu'ils génèrent. Des tests expérimentaux permettent alors de valider la majorité des modèles ou au contraire de revenir sur certains d'entre eux pour les optimiser ou montrer leurs limitations. Les résultats ont en effet montré l'importance de la prise en compte des contraintes résiduelles, qui provoquent une déformée initiale des speaklets et modifient leurs fréquences propres, rendant alors l'utilisation de grands rayons inefficaces. En accord avec les modèles, les speaklets possèdent un comportement dynamique linéaire, ce qui permet de les caractériser à l'aide de fonctions de transfert. La théorie et les enregistrements sonores montrent alors qu'un HPDM composé de tels speaklets permet dans le meilleur des cas de produire une pression identique à celle générée par la même matrice pilotée en analogique. Dans notre cas, des taux de distorsions supérieurs ont été obtenus lors des reconstructions digitales, à cause des réponses non uniformes des speaklets, dues à des résistances d'accès différentes. Le HPDM présenté possède cependant d'autres avantages, le plus important étant la très faible consommation électrique qu'il est théoriquement possible d'atteindre en utilisant les méthodes de charges et de décharges adiabatiques. Le HPDM piézoélectrique MEMS apparait donc comme étant une technologie prometteuse. L'optimisation de notre premier prototype à l'aide des outils développés doit en effet conduire à un HPDM capable de générer une pression équivalente à celle obtenue en mode analogique, mais avec un rendement électroacoustique beaucoup plus important. Les futurs travaux devront ensuite se concentrer sur la conception de speaklets non-linéaires et sur la forme du pulse de pression qu'ils génèrent, afin de gagner en niveau sonore. / The Digital Loudspeaker Array (DLA) is an electroacoustic transducer which receives as input a digital signal and performs the analog conversion directly into the air. It consists of a plurality of radiating elements arranged in a matrix. These elements will be designated by the term “speaklet” when they are reduced in size. The acoustic radiation of a DLA is indeed very sensitive to the size of the matrix due to differences in path length, which makes it especially suitable for MEMS technology. This thesis is on the study of a piezoelectric MEMS DLA. After an introduction that is increasingly focused on the subject, the thesis addresses the multiphysics modeling of the DLA, dimensioning of the speaklets and experimental tests. Analytical formulas, numerical simulations and finite element models are developed and used to predict the mechanical behavior of the presented speaklets, the pressure radiated by the DLA and the electrical power consumption. The speaklet are then dimensioned from the technological stack (set in advance) in order to maximize the pressure level. Experimental tests involving the use of an anechoic chamber, an optical interferometer, a vibrometer and an impedancemeter validate most of the models. Otherwise, these tests are usefull for improving some of them or for showing their limitations. The results have shown the importance of the residual stresses, which cause an initial deformation of the speaklets and modify their resonance frequencies, thus rendering ineffective the use of large radii. In accordance with the models, the static deflection of the speaklets is nonlinear but their dynamic behavior is linear. This enables characterizations using transfer functions. Theory and sound recordings show that a DLA made of such speaklets can produce in the best case the same pressure to that generated by the same matrix driven in an analog way. In our case, more distortions were obtained in digital reconstructions because of non-uniform responses of the speaklets, due to different access resistances. However, the presented DLA has other advantages, the most important being the very low power consumption it is theoretically possible to achieve using the adiabatic charge principle. The piezoelectric MEMS DLA thus appears as a promising technology. The optimization of our first prototype using the developed tools should indeed lead to a DLA able to generate an equivalent presure to that obtained with analog control, but with a far greater electroacoustic efficiency. Future work should then focus on the design of nonlinear speaklets and on the shaping of the pulse of pressure they generate, in order to increase the total pressure level.

Cette thèse se concentre sur l'étude numérique de l'interaction laser ultrabref avec les diélectriques transparents. En particulier, le phénomène d'auto-organisation des nanoréseaux dans la silice est discuté et un modèle multiphysique est proposé pour expliquer le mécanisme de leur formation. Les nanoréseaux en volume sont des nanostructures périodiques de périodicité sub-longueur d'onde, qui consistent en un matériau moins dense et sont générés par une irradiation laser multi-impulsionnelle femtoseconde dans certains verres, cristaux et semiconducteurs. Leur origine physique ainsi que les conditions d'irradiation laser pour leur formation et leur effacement sont investiguées dans ce travail théorique. Pour simuler la propagation nonlinéaire dans les verres, les équations de Maxwell sont couplées avec l'équation d'évolution de la densité électronique. Il est démontré que les nanoplasmas périodiques 3D sont formés pendant l'interaction laser ultrabref avec les inhomogénéités de la silice fondue. Les nanopores induits par laser sont supposés jouer le rôle de centres inhomogènes de diffusion. La périodicité sub-longueur d'onde et l'orientation des nanoplasmas dépendante de la polarisation, révélées dans cette thèse, font d'eux un excellent candidat pour expliquer la formation des nanoréseaux en volume. En plus, il est demontré que les nano-ripples sur la surface de silice fondue et les nanoréseaux en volume ont des mécanismes de formation similaires. Pour justifier la présence de nanopores dans la silice fondue irradiée par laser, les processus de décomposition du verre sont étudiés. Premièrement, les profils de température sont calculés sur la base d'un modèle électron-ion. Ensuite, à partir des températures calculées, des critères de cavitation et de nucléation dans le verre ainsi que des équations hydrodynamiques de Rayleigh-Plesset, les conditions pour la formation des nanopores et la survie des nanoréseaux en volume sont élucidées. Pour établir les dépendances des paramètres du laser de formation et d'effacement des nanoréseaux en volume, l'approche multiphysique est développée comprenant la propagation du laser ultrabref dans le verre, les processus d'excitation/relaxation électroniques et le modèle à deux températures. Les résultats numériques fournissent les paramètres du laser en fonction de l'énergie de l'impulsion, sa durée et le taux de répétition pour induire des nanoréseaux en volume, en bon accord avec les expériences nombreuses et indépendantes de la littérature. Le travail réalisé a non seulement permis de déterminer les mécanismes de formation des nanostructures périodiques mais améliore également notre connaissance du contrôle optimal des paramètres du laser sur la réponse ultrarapide d matériau, en ouvrant des nouvelles opportunités de traitement des diélectriques par laser ultrabref / This thesis is focused on the numerical modeling of ultrashort laser interaction with transparent dielectrics. More particularly, the phenomenon of self-organized volume nanogratings in fused silica bulk is discussed and a multiphysical model is proposed to explain the mechanism of their formation. Volume nanogratings are sub-wavelength periodic nanostructures, consisting of less dense material, which are commonly induced by multipulse femtosecond laser irradiation in some glasses, crystals and indirect semiconductors. Their physical origin as well as the laser irradiation conditions for theirformation and erasure are investigated in this theoretical work. To model the nonlinear propagation inside glass, Maxwell's equations are coupled with rate equation. It is shown that three-dimensional periodic nanoplasmas are formed during ultrashort laser interaction with fused silica inhomogeneities. Laser-induced nanopores are proposed to play the role of inhomogeneous scattering centers. Subwavelength periodicity and polarization dependent orientation of the nanoplasmas, revealed in this thesis, make them a strong candidate for explaining volume nanogratings formation. Additionally, it is demonstrated that the nanoripples on fused silica surface and volume nanogratings have similar formation mechanisms. To justify the presence of nanopores in laser-irradiated fused silica bulk, glass decomposition processes are investigated. Firstly, the temperature profiles are found by incorporating the electron-ion temperature model. Then, based on the calculated temperatures, criteria for cavitation and nucleation in glass and also hydrodynamic Rayleigh-Plesset equation, the conditions for nanopores formation and for volume nanogratings survival are elucidated. To define the laser parameter dependencies on the volume nanogratings formation/erasure, a selfconsistent multiphysical approach is developed including ultrafast laser propagation in glass, multiple rate equation to take into account excitation/relaxation processes and two-temperature model. The numerical results provide a laser parameter window as a function of laser pulse energy, laser pulse duration and repetition rate for volume nanogratings consistent with numerous independent experiments. The performed work not only provides new insights into the formation mechanisms of periodic nanostructures but also improves our knowledge of the optimal laser parameter control over ultrafast material response, opening new opportunities in ultrashort laser processing of dielectrics

La montée en puissance des systèmes de production de l'énergie électrique, à partir de centrales éoliennes ou photovoltaïques, a fait apparaître un besoin aigu en moyens de stockage de cette énergie. Les technologies d'accumulateurs électrochimiques sont à même de répondre à cette problématique, en particulier les batteries sodium chlorure de nickel, une technologie d'accumulateur à haute température. Mais afin de gérer au mieux les performances de ces batteries, il est nécessaire d'avoir une connaissance approfondie du comportement électrothermique des cellules unitaires les composant. Ce travail de thèse présente en détail le développement d'un modèle multiphysique 2D d'une cellule Na-MCl2 commerciale, à même de simuler son fonctionnement en décharge à courant constant. Ce modèle a été construit sur une étude approfondie des mécanismes électrochimiques à l'œuvre dans ce type de technologie, et des effets thermiques associés. Il repose également sur la synthèse et la critique des modèles proposés précédemment dans la littérature. Le modèle inédit proposé dans ce manuscrit permet de prendre en compte deux réactions électrochimiques simultanées, ainsi que le comportement thermique de la cellule. Le modèle a été validé en confrontant les résultats de simulations aux performances électriques réelles d'une cellule commerciale Na-MCl2 (incluant une caractérisation des phénomènes électrochimiques par voltammetrie cyclique). Ce modèle permet de simuler avec succès des décharges à courant constant à différents régimes, et à partir de différentes températures. Il est également capable de prévoir les effets électrothermiques qui seraient consécutifs à un changement de design de la cellule (design structurel et de constitution). / The growing production capacity of wind turbines or photo-voltaic plants has revealed an acute need for electric energy storage systems. Electrochemical accumulator technologies are able to address this issue, in particular sodium nickel chloride batteries, a high temperature battery system. But to enhance the energy management of these batteries, it is necessary to have a deep knowledge of the electro-thermal behavior of its unit cell component. This thesis presents the development of a comprehensive multiphysics 2D model of a commercial Na-MCl2 cell, able to simulate its operation in constant current discharge. This model was built on a thorough study of electrochemical mechanisms at work with this type of technology, and its associated thermal effects. It is also based on the synthesis of existing models, and their critique. The new model proposed in this manuscript can take into account two simultaneous electrochemical reactions, as well as the thermal behavior of the cell. The model was validated by comparison between simulation results and the actual electrical performances of a commercial Na-MCl2 cell (including the characterization of electrochemical phenomena by cyclic voltammetry). This model successfully simulates constant current discharges at different rates, and from different temperatures. It is also capable of predicting the electro-thermal effects consecutive to design changes of the cell (structural and relevant to the active material).

La fabrication additive est l’un des secteurs industriels les plus en développent ces dernières années. L’une de ces technologies de fabrication les plus prometteuses est la fusion laser sélective (SLS), et relève d’un intérêt croissant aussi bien industriel qu’académique. Néanmoins, beaucoup de phénomène mis en jeu par ce procédé demeure non encore bien compris, entravant ainsi son développement pour la production de pièces de bonne qualité pour des applications industrielles. L’objectif de cette thèse est de développer un cadre de simulation numérique permettant la simulation du procédé SLS pour des poudres de polymère afin de comprendre les multiples et complexes phénomènes physiques qui se produise lors du frittage laser et d’étudier l’influence des paramètres du procédé sur la qualité du produit final. Contrairement aux approches classiques de modélisation numérique, basées sur la définition de matériaux homogène équivalents pour la résolution des équations de bilan, nous proposons une simulation globale du procédé du frittage laser de poudres, en utilisant la méthode des Eléments Discrets (DEM). Cela consiste en un couplage entre quatre sous-modèles : transferts radiatif dans le milieu granulaire semi-transparent, conduction thermique dans les milieux discrets, coalescence puis densification. Le modèle de transferts par rayonnement concerne l’interaction du faisceau laser avec le lit de poudre. Plusieurs phénomènes sont ainsi pris en compte, notamment la réflexion, la transmission, l’absorption et la réfraction. De plus, une méthode de Monte-Carlo couplée à la méthode du Lancer de rayons est développée afin d’étudier l’influence de la réfraction sur la distribution de l’énergie du laser dans le lit de poudre. Le modèle de conduction dans des milieux discrets décrit la diffusion thermique inter-particules. Finalement, le modèle de frittage décrit les cinétiques de coalescence et de diffusion de l’air dans le polymère et densification du milieu. Cela permet de décrire les cinétiques de fusion des grains, dont l’énergie de surface et la diffusons de l’air sont les deux moteurs principaux. Le couplage entre les différents modèles nous a permis de proposer un modèle numérique global, validé grâce à des comparaisons à des résultats de simulations théoriques et expérimentales, trouvés dans la littérature. Une analyse paramétrique est alors proposée pour la validation du modèle et l’étude du procédé. L’influence de différents paramètres aussi bien du procédé que du matériau sur le champ de température, la densité relative du matériau sa structure, etc , est ainsi investiguée. Les résultats montrent une bonne précision dans la modélisation des différents phénomènes complexes inhérents à ce procédé, et ce travail constitue un potentiel réel pour la modélisation et l’optimisation des procédés de fabrication additive par matériaux granulaires. / Many industrial and academic interests concerning the additive manufacturing processes are developed in the last decades. As one of the most promising technique of additive manufacturing, the Selective Laser Sintering (SLS) has been valued by both industry and academic. However, it remains that several phenomena are still not well understood in order to properly model the process and propose quality improvement of parts made. The goal of this Ph.D. project is to develop a framework of numerical simulation in order to model the SLS process in polymer powder bed, meanwhile understanding multiple physical phenomena occurring during the process and studying the influence of process parameters on the quality of final product. In contrast to traditional approach, based on the equivalent homogeneous material in numerical modeling of partial differential equations derived from conservation laws, we propose a global model to simulate powder-based additive manufacturing by using the Discrete Element method (DEM). It consists in a coupling between four different physical models: radiative heat transfer, discrete heat conduction, sintering and granular dynamics models. Firstly, the submodel of radiative heat transfer concerns the interaction between the laser beam and powder bed. Several phenomena are considered, including the reflection, transmission, absorption and scattering. Besides, a modified Monte Carlo ray-tracing method is developed in order to study the influence of scattering on the distribution of the deposited laser energy inside the powder bed Furthermore, the submodel of discrete heat conduction describes the inter-particles heat diffusion. Moreover, the sintering submodel concerns the phenomena of coalescence and air diffusion. It describes the melting kinetics of grains, driven by surface tension and the release of entrapped gases inside powder bed. Finally, the granular dynamics submodel concerns the motions and contacts between particles when depositing a new layer of powders. The coupling between these submodels leads to propose a global numerical framework, validated by comparing the results to both simulated and experimental ones from literatures. A parametric study is then proposed for model validation and process analysis. The Influence of different material and process parameters on the evolution of temperature, relative density and materials structure and characteristics are investigated. The results exhibit accurate modeling of the complex phenomena occurring during the SLS process, and the work constitute a great potential in modeling and optimization of additive processes.

Beaucoup de procédés utilisés dans l'industrie sidérurgique (coulée de lingots,coulée continue, …) peuvent générer des défauts : macro-ségrégation, mauvaises propriétés de la microstructure, défauts surfaciques. Ces problèmes peuvent être résolus par un contrôle de la température et de l’écoulement d'acier liquide. Le brassage électromagnétique (EMS) est une technique largement utilisée pour contrôler l’écoulement d'acier liquide par l’imposition d'un champ électromagnétique. Cette technique est complexe car elle couple plusieurs types de problèmes physiques:écoulement multiphasique, solidification,transfert de chaleur et induction électromagnétique à basse fréquence.En outre, l’approche expérimentale est difficile de par la dimension,l'environnement et le coût des procédés considérés. Pour ces raisons, des simulations numériques efficaces sont nécessaires pour comprendre les applications EMS et améliorer les procédés évoqués. L'objectif de cette thèse est de développer une méthodologie numérique robuste,efficace et précise pour la simulation multi-physique de l'EMS, en particulier pour le brassage dans le moule dans le cadre de la coulée continue d'acier. Cette méthodologie a été mise en oeuvre dans le code commercial THERCAST® pour être utilisé dans le cadre d’applications industrielles / Many of the processes used in thesteelmaking industry (e.g. ingot casting,continuous casting, …) can lead todefects: macro-segregation, poormicrostructure properties, surfacedefects. These issues can be solved bycontrolling the temperature and the flowof molten steel. Electromagnetic stirring(EMS) is a widely used technique to steerthe flow of liquid steel by thesuperimposition of an electro-magneticfield. This application is complex becauseit couples several physical problems:multi-phase flow, solidification, heattransfer and low frequency electromagneticinduction. In addition,experimental work is difficult because ofthe size, environment and cost of theconsidered processes. For thesereasons, efficient and effective numericalsimulations are needed to understandEMS applications and improve theaforementioned processes.The objective of this thesis is to developa robust, efficient and accurate numericalprocedure for the multi-physicssimulation of EMS, especially for in-moldstirring in the framework of continuouscasting of steel. This procedure has beenimplemented in the commercial codeTHERCAST® in order to be used forindustrial applications.