La croûte intra-océanique du sud-Kohistan comprend les complexes de Spat et Jijal, les amphibolites de Kamila et le complexe de Chilas. Le complexe de Jijal est une intrusion litée constituée de dunite, clinopyroxénite à olivine, websterite, webstérite à grenat et de gabbro à grenat. Mis en place à l'interface manteau croûte, il résulte de la cristallisation fractionnée sous haute pression et haute température (10-12 Kbars; 800- 1100°C) d'un liquide tholéiitique d'arc intra-océanique hautement magnésien. Les roches de Thak Gah (complexe de Chilas) ont cristallisé sous plus faible pression et température (5 - 6 kbars, 800 à 950°C) dans un réservoir magmatique plus superficiel. Les métaplutonites de Chilas et Kamila présentent également des signatures géochimiques analogues à des tholéiites d'arc insulaire. Au contraire, les amphibolites à grains fins du Groupe de Kamila apparaissent comme des métavolcanites de type N-MORB, elles représentent des reliques de croûte océanique et le substratum de l'arc du Kohistan. Plusieurs horizons enrichis en PGE caractérisent le complexe de Jijal. Des enrichissements en Ir,Ru, Ni corrélés au Cr dans les dunites, indiquent que l'irridium et le ruthénium précipitent précocément avec la chromite sous haute fugacité d'oxygène. Des enrichissements en Ir et Ru du même type caractérisent les dunites chromifères des complexes de Thak Gah et de Spat. Dans les webstérites à grenat de Jijal, des sulfures disséminés (chalcopyrite, pentlandite, pyrrhotite± pyrite) et des minéraux du Groupe du platine (témagamite, monchéite, mérenskyite et sperrylite) correspondent à des enrichissements en Ir, Ru, Pt, Pd, Au, Ni et Cu. Ils résultent de la formation d'un liquide sulfuré consécutif à la baisse de f02. Dans les gabbros à grenats, des enrichissements en Pt, Pd et Cu traduisent l'action de fluides hydrothermaux tardi-magmatiques. Ces minéralisations diffèrent de celles du Bushveld. La saturation en soufre induite par les processus de cristallisation sous haute pression apparaît le moteur essentiel des dépôts observés dans les webstérites à grenats. De telles minéralisations représentent un exemple type de dépôts de PGE associés aux arcs insulaires.

De la bactérie à l'homme, dispersées ou en tandem, les répétitions peuvent représenter jusqu'à 90 % de la séquence génomique. Malgré leur impact sur la plasticité et l'évolution des génomes eucaryotes, leurs mécanismes de formation sont encore très spéculatifs. L'insertion continue de nouvelles répétitions devrait conduire à une augmentation constante de la taille des génomes. Or, il ne semble pas que ce soit le cas. Y a t-il régulation de la taille des génomes? Le processus de régulation est-il le même dans l'euchromatine et l'hétérochromatine? Afin d'étudier la dynamique des répétitions, j'ai développé un ensemble de programmes informatiques pour la détection des duplications segmentaires (DS) et des répétitions en tandem (RT). A partir des caractéristiques des DS détectées chez Drosophila melanogaster, j'ai proposé un modèle de formation des DS, basé sur un modèle de recombinaison homologue non-allélique. J'ai également identifié les traces de l'implication des éléments transposables (ET) dans ce processus. Afin de caractériser la relation existante entre les répétitions et la structure de la chromatine, j'ai ensuite réalisé une analyse comparative de la dynamique des répétitions euchromatiques et hétérochromatiques. Pour ce travail, nous avons choisi comme modèle d'étude Arabidopsis thaliana. La construction d'arbres phylogénétiques des séquences répétées m'a permis de dater les répétitions. Nous suggérons ainsi une propagation par « vague » des ET. J'ai ensuite estimé les forces d'élimination des copies d'ET. Nos résultats suggèrent que dans l'euchromatine, la pression de sélection due aux gènes induit l'élimination des répétitions. Dans l'hétérochromatine, la faible densité en gènes permet de maintenir une forte densité en ET. Pourtant, les estimations du taux de perte en ADN, prédisent un turnover aussi rapide dans l'euchromatine que dans l'hétérochromatine. Afin de contre-balancer l'insertion des ET dans l'hétérochromatine, nous pouvons invoquer la recombinaison homologue non-allélique.

L'objectif principal de cette étude consiste à réduire les temps de calcul des simulations de procédés de mise en forme continus sous le logiciel Forge3®. Ces procédés, tel que le laminage et le tréfilage, sont caractérisés par des pièces dont la longueur est très importante en comparaison des dimensions de la section ainsi que de la zone de contact. Une approche incrémentale générale implique des temps de calcul conséquents allant de quelques heures à plusieurs jours. En se concentrant sur le régime quasi permanant de ces procédés, une formulation stationnaire est développée pour accélérer leur simulation. Le domaine de calcul correspond initialement à une estimation de la forme de l'écoulement solution au voisinage des outils. Une étape de correction du domaine est ajoutée en plus du calcul stationnaire de l'écoulement. Comme les conditions aux limites sont modifiées, ces deux étapes sont répétées jusqu'à convergence.L'étude s'est concentrée principalement sur l'étape de correction du domaine correspondant à la résolution d'un problème de surface libre par la méthode des éléments finis. Le caractère purement convectif du problème ainsi que la prise en compte du contact nécessite l'utilisation de formulations faibles faisant apparaître un décalage amont (SUPG). Deux nouvelles formulations basées sur la méthode des moindres carrés sont développées avec succès (MC_supg et MC_lc). Pour appliquer la méthode à des géométries complexes, différentes méthodes de généralisations sont développées où un 2ème degré de liberté est ajouté aux nœuds de surface. La méthode la plus performante (CSL_dif) consiste à utiliser ce 2ème degré de liberté pour le calcul de surface libre uniquement sur les nœuds appartenant à une arête géométrique, pour les autres une régularisation du maillage dans la direction tangente y est effectuée. Des résultats excellents ont été observés sur un grand nombre de cas tests analytiques. Le contact est appliqué par une méthode de pénalisation aux nœuds. Afin de renforcer le couplage entre cette étape et celle du calcul de l'écoulement, un contact bilatéral glissant est attribué aux nœuds en compression alors que pour les autres nœuds un contact unilatéral est employé. Un algorithme spécifique est développé pour déterminer avec précision la zone de contact.Cette formulation itérative pour la recherche de l'état stationnaire a été appliquée avec succès sur un grand nombre de cas tests de mise en forme. Des accélérations comprises entre 10 et 60 ont été obtenues par rapport à Forge3®.

Les AMF à base fer se distinguent par la présence d'un couplage entre la transformation de phase et le glissement plastique. Pour caractériser leur comportement thermomécanique, des essais cycliques de chargement mécanique effectués à différentes températures et à différents niveaux de chargement, suivis par chauffage, ont été effectués. Ceux ci nous ont permis de distinguer les spécificités des AMF à base fer comme la non linéarité des interactions inter et intragranulaires, l'évolution des contraintes critiques avec la température, l'effet de la déformation plastique sur celle de transformation et l'activation de la transformation inverse uniquement par chauffage. En s'inspirant des résultats de ces essais, une loi de comportement, intégrant toutes ces spécificités, est proposée. Elle dérive d'une expression de l'énergie libre de Gibbs issue de considérations micromécaniques. Elle présente deux variables internes, la fraction volumique de martensite pour décrire la transformation de phase et le taux de plasticité pour le comportement plastique. Cette loi a été implémentée dans le code éléments finis Abaqus via la subroutine UMAT. Elle a été validée par comparaison des simulations numériques avec les résultats expérimentaux. Suite à cette validation, des applications en AMF à base fer ont été étudiées. Les résultats obtenus ont montré la capacité du modèle à être utilisé comme outil de dimensionnement de structures en AMF à base fer / The Fe-based shape memory alloys (SMAs) present a specific thermomechanical behaviour compared with classical SMAs. In this PhD thesis, experimental thermomechanical tests were performed in order to study such behavior. The applied loading is a tension followed by a significant heating. The loading cycle is repeated at different constant temperatures and maximum stresses. The experimental results show a coupling between two non linear inelastic mechanisms: phase transformation and plasticity. The reverse transformation activated only during heating, the effect of plastic strain on the transformation one, were also analysed. Taking into account all these specificities, a finite element numerical tool adapted to Fe-based SMA structural analysis is proposed. It is based on a developed constitutive model which describes the effect of phase transformation, plastic sliding and their interactions on the thermomechanical behavior. Two scalar internal variables were considered to describe phase transformation and plastic sliding effects. This model was derived from an assumed expression of the Gibbs free energy taking into account, in addition to mechanical and chemical quantities, the non linear interaction quantities related to inter- and intra-granular incompatibilities. The numerical tool derived from the implicit resolution of the non linear partial derivative constitutive equations was implemented in the Abaqus finite element code via the UMAT subroutine. After verification tests for homogeneous and heterogeneous thermo-mechanical loadings, two examples of Fe-based SMA applications were studied. They correspond to Fe-based SMA tightening systems: a fish plates for crane rails and a ring for tubes connection

Les instabilités fluides peuvent dégrader le confinement du plasma au sein des tokamaks. Étant données les échelles spatio-temporelles, on choisit les modèles fluides obtenus à partir de la dérivation des modèles cinétiques. On dérive plusieurs modèles hiérarchiques de la MagnétoHydroDynamique (MHD) et en particulier les modèles de la MHD réduite du Current Hole et de l'équilibre de Grad-Shafranov. Une des difficulté de l'ensemble de ces modèles est de respecter l'équation modélisant l'absence de monopôles magnétiques. Pour assurer cette condition en tout point du domaine, le champ magnétique est réécrit avec un potentiel vecteur. L'utilisation de potentiels fait apparaître des équations faisant intervenir des dérivées d'ordre supérieurs. La stratégie numérique développée est l'utilisation de la méthode des éléments finis avec des éléments C1 de Bell. Sur un maillage non structuré, ces éléments ont l'intérêt de présenter une base réduite définir exclusivement avec des variables aux noeuds du maillage. Les modèles de MHD réduite du Current Hole et de Grad-Shafranov ont été résolus avec ces éléments. La résolution du cas test de Grad-Shafranov avec les conditions de bords exactes a permis d'obtenir l'ordre optimale de 5. La résolution du système du Current Hole avec ces éléments, validée par l'obtention du paramètre η1/3, a permis l'observation de développement d'instabilités en dents de scies. / Fluid instabilities can degrade plasma confinement in tokamaks. Given the spatial and temporal scales, we choose the fluid models obtained from the derivation of kinetic models. We derived several hierarchical models of MagnetoHydroDynamic (MHD) and in particular models of reduced MHD like the Current Hole and the Grad-Shafranov equilibrium. One of the difficulty of all these models is to respect the absence of magnetic monopoles equation. To ensure this condition at any point, the magnetic field is rewritten with a vector potential. The use of vector portential implies that higher order derivatives appear in the equation. The numerical strategy is developed using the finite element method with C1 Bell's elements. On a unstructured mesh, these have the advantage to present a reduced basis with degrees of freedom defined exclusively on the nodes of the mesh. The reduced MHD models of the Current Hole and Grad-Shafranov have thus been resolved with these elements. The resolution of a Grad-Shafranov test case with exact boundary conditions yields the optimal order of 5. The resolution of the Current Hole system with thesse elements has been validated by obtaining physical parameter η1/3 and allowed the observation of the development of sawtooth instabilities.

La fissuration à chaud est un défaut majeur rencontré en solidification des alliages d'aluminium. Elle est liée à l'incapacité du liquide de s'écouler dans les zones où des porosités sont présentes, ne permettant pas de les refermer avant qu'elles gagnent en volume. Pour comprendre la fissuration à chaud, il est crucial de développer nos connaissances du comportement mécanique de la zone pâteuse. Pour cela, il est très utile d'effectuer des expériences de microtomographie aux rayons X et des simulations mécaniques sur des volumes élémentaires représentatifs. Dans cette thèse, nous proposons de coupler les deux approches en initialisant une simulation par éléments finis grâce à des données de microtomographie issues d'un test de traction isotherme d'un alliage d'aluminium-cuivre à l'état pâteux. Cette approche originale nous donne directement accès à la réalité expérimentale et permet des comparaisons des évolutions numérique et expérimentale de l'éprouvette. Nous expliquons dans un premier temps comment obtenir la représentation numérique à l'aide de l'algorithme des marching cubes et de la méthode d'immersion de volume. Nous présentons ensuite notre modèle numérique qui s'appuie sur une résolution monolithique des équations de Stokes. Une fois le champ de vitesse obtenu dans l'ensemble des phases solide, liquide et gazeuse, nous utilisons une méthode level set dans un formalisme eulérien afin de faire évoluer la morphologie de notre échantillon numérique. Malgré la simplicité du modèle, les résultats expérimentaux et numériques montrent un accord raisonnable en ce qui concerne la propagation de l'air à l'intérieur de l'échantillon. / Hot tearing is a major defect arising during solidification of aluminium alloys. This defect is associated with the inability of liquid to feed areas where voids have started to appear, not allowing to heal small defects before they grow bigger. To understand hot tearing, it is mandatory to develop a good knowledge of the semi-solid mechanical behaviour. It is thus very useful to carry out X-ray microtomographies experiments and mechanical simulations on representative elementary volumes. In this work, we couple the both approaches by initialising a finite element simulation with the help of microtomography data obtained during an isothermal tensile testing of an aluminium-copper alloy in the mushy state. This innovative approach gives a direct access to the experimental reality and allows comparisons of numerical and experimental evolutions of the sample. We explain in a first time how to get the numerical representation thanks to a marching cubes algorithm and the immersed volume method. Then, we present our numerical model for which we solve the Stokes equations in a monolithic way. Once the velocity computed in all the solid, liquid and gaseous phases, we use a level set method in a Eulerian formalism to obtain the morphological evolution of our numerical sample. Despite the model simplicity, numerical and experimental results show a reasonable agreement concerning the air propagation inside the sample.

Les vibrations des systèmes de forage pétrolier sont à l'origine de nombreux dysfonctionnements (ruptures des tiges par une fatigue accélérée, réduction des performances, endommagement des outils de mesures, endommagement des parois du puits, etc.). Face à la complexité des puits forés aujourd'hui, la maîtrise des vibrations des systèmes de forage est plus que jamais un enjeu majeur dans la réussite économique d'un projet pétrolier. Durant l'opération de forage, les tiges en rotation entrent en interaction avec les parois du puits (tubage et/ou formation) et encaissent dans certains cas des vibrations sévères. On distingue généralement trois modes de vibrations suivant le plan de leur occurrence : axiales, latérales et de torsion. Nous ne nous intéressons dans ce mémoire qu'aux vibrations latérales des tiges de forage. Si les vibrations latérales ont fait l'objet de nombreuses études dans le passé, il reste néanmoins des axes d'amélioration possible, tant sur la compréhension des phénomènes (contact garniture-puits par exemple) que sur la recherche de méthodes numériques permettant de réduire les temps de calcul. Dans le cadre de cette thèse, un modèle a été développé pour étudier les vibrations latérales des garnitures de forage dans des forages à trajectoires complexes. Ce modèle permet de prédire les vibrations latérales des tiges pour des paramètres opératoires donnés (vitesse de rotation, poids sur l'outil de forage). Un modèle numérique en éléments finis a été développé pour résoudre les équations du mouvement et analyser ainsi la sensibilité des vibrations aux paramètres opératoires du forage en particulier la vitesse de rotation et l'effort axial dans les tiges. Le modèle permet en outre d'analyser la réponse dynamique d'une garniture en cours du forage (conception). En outre, cette étude a permis de mieux élucider le phénomène sévère de précession des tiges (whirling), très nuisible à l'intégrité mécanique des systèmes de forage. Un nouveau banc d'essais a été mis au point par le Centre de Géosciences de Mines ParisTech pour reproduire les vibrations latérales, mieux comprendre le phénomène du whirling et valider les résultats numériques du modèle. Par ailleurs, des mesures dynamiques en surfaces et en fond de puits au cours de forages réels ont été analysées afin de mettre en évidence les vibrations latérales les plus sévères et en particulier le whirling. Ces données de terrain ont permis de comparer les fréquences propres du système mesurées et celles fournies par le modèle numérique. / Drillstring vibrations are commonly observed during oil & gas well drilling operations. Vibrations are a major cause of drilling tools dysfunction (drillstring breaking because of fatigue, reduced drilling efficiency, measurement-while-drilling tools failure, damaging of drill bits, etc.). Because of the increasing complexity of oil & gas wells drilled nowadays, operators need to mitigate efficiently the drillstring vibrations in order to successfully achieve the drilling process. During the drilling operation, rotating drillstrings are in interaction with the well borehole (casing and/or rock) which may lead to severe vibrations. Different vibrations modes occur simultaneously while drilling, we identify mainly three modes: axial, torsional and lateral. This work deals only with lateral vibrations. Literature survey papers show numerous experimental and numerical studies carried out on drillstring dynamics. The developed models don't take into account sufficiently the complex drillstring-borehole interactions or the efficient numerical methods needed to reduce the computation time. A new drillstring dynamics model has been developed within this thesis in order to compute the lateral vibrations of drillstrings in a complex well trajectory. Given the operating parameters (rotary speed, weight on bit) the model predicts the dynamics response of the drillstring in terms of lateral vibrations. A finite element model has been implemented to solve for the equations of motion of the dynamics model and study the dependence of the lateral vibrations on some operating parameters of the system, mainly the rotary speed and the axial load on the drillstrings. The finite element model can be used to compute and enhance the dynamic response of a given drillstring configuration for design issues. Besides, the model can be used to understand some dynamic phenomena encountered while drilling (post-analysis). Moreover, this study was useful to better understand the “whirling” phenomenon which is very harmful for the drilling system components. A new lateral vibrations simulator was built at Mines ParisTech in order to understand the whirling phenomenon and validate the numerical results provided by the dynamics model. Surface and downhole fieldmeasurements have been analyzed in order to understand the occurrence of whirling. The eigenfrequencies evaluated from the field data have been found very close to those provided by the dynamics model.

Dans cette thèse, nous avons étudié la simulation numérique des phénomènes d'interaction fluide-structure entre un fluide compressible et une structure déformable. En particulier, nous nous sommes intéressés au couplage par une approche partitionnée entre une méthode de Volumes Finis pour résoudre les équations de la mécanique des fluides compressibles et une méthode d'Eléments discrets pour le solide, capable de prendre en compte la fissuration. La revue des méthodes existantes de domaines fictifs ainsi que des algorithmes partitionnés couramment utilisés pour le couplage conduit à choisir une méthode de frontières immergées conservative et un schéma de couplage explicite. Il est établi que la méthode d'Eléments Discrets utilisée permet de retrouver le comportement macroscopique du matériau et que le schéma symplectique employé assure la préservation de l'énergie du solide. Puis nous avons développé un algorithme de couplage explicite entre un fluide compressible non-visqueux et un solide indéformable. Nous avons montré des propriétés de conservation exacte de masse, de quantité de mouvement et d'énergie du système ainsi que de consistance du schéma de couplage. Cet algorithme a été étendu au couplage avec un solide déformable, sous la forme d'un schéma semi-implicite. Cette méthode a été appliquée à l'étude de problèmes d'écoulements non-visqueux autour de structures mobiles : les comparaisons avec des résultats numériques et expérimentaux existants démontrent la très bonne précision de notre méthode / This work aims at the numerical simulation of compressible fluid/deformable structure interactions. In particular, we have developed a partitioned coupling algorithm between a Finite Volume method for the compressible fluid and a Discrete Element method capable of taking into account fractures in the solid. A survey of existing fictitious domain methods and partitioned algorithms has led to choose an Embedded Boundary method and an explicit coupling scheme. We first showed that the Discrete Element method used for the solid yielded the correct macroscopic behaviour and that the symplectic time-integration scheme ensured the preservation of energy. We then developed an explicit coupling algorithm between a compressible inviscid fluid and an undeformable solid. Mass, momentum and energy conservation and consistency properties were proved for the coupling scheme. The algorithm was then extended to the coupling with a deformable solid, in the form of a semi-implicit scheme. Finally, we applied this method to unsteady inviscid flows around moving structures: comparisons with existing numerical and experimental results demonstrate the excellent accuracy of our method

Les éléments transposables (ETs) sont des éléments génétiques mobiles capables de se déplacer et de se multiplier au sein d’un génome. Identifiés dans la plupart des espèces vivantes incluant les bactéries, mais longtemps considérés comme de l’ADN poubelle, aujourd’hui les ETs sont indéniablement des acteurs majeurs impliqués dans l’évolution des gènes, des génomes et des organismes. Si à l’échelle des individus les ETs peuvent avoir des effets délétères pouvant entrainer des maladies, à plus grande échelle ils sont de puissants agents évolutifs impliqués dans la plasticité génomique. Ces « parasites » peuvent également être sources de nouveaux matériels génétiques comme des promoteurs ou même de nouveaux gènes avec de nouvelles fonctions pour l’hôte. Les objectifs majeurs de mon travail de thèse ont été de déterminer les différentes familles d’ETs présentes dans les génomes de poissons, la part que chacune d’entre elles occupe dans ces génomes et enfin de comprendre l’histoire évolutive des familles d’ETs dans les génomes de poissons en comparaison avec les autres génomes de vertébrés. Cette comparaison à grande échelle permettra de comprendre les différentes stratégies évolutives des ETs. D’autre part, j’ai étudié deux gènes de vertébrés, Gin-1 et Gin-2 dérivés d’ETs, dans le but de comprendre leurs origines et évolution au sein des vertébrés ainsi que d’émettre des hypothèses quant à leur fonction moléculaire potentielle encore inconnue. Pour cela, des analyses in silico ont permis de mieux comprendre les origines de ces gènes. Gin-1, présent chez les amniotes, et Gin-2, absent uniquement des mammifères placentaires, dérivent tous deux de transposons GIN. / Transposable elements (TEs) are mobile genetic elements - able to move and to multiply within genomes - identified in almost all living organisms including bacteria. Considered as junk DNA for long, nowadays they are undeniably major players of gene, genome and host evolution. TEs can be deleterious causing diseases but these “parasites” can also be source of new genetic materials as promoters or even new genes bringing new functions for hosts. The objectives of my thesis was to determine the presence or not of the different TE families in vertebrate genomes, as well as their respective content to understand their evolutionary history. I performed a large-Scale comparative analysis to highlight the various evolutionary strategies of TEs. I showed that TE content is highly variable in vertebrate genomes, the smallest and the largest being found in fish, and may contribute to their genome sizes especially in fish. These superfamilies underwent differential waves of activity in vertebrate species highlighting TE dynamics. On another hand, I focused on the study of a vertebrate-Specific TE-Derived gene, named Gin-2, to understand its origin, evolution, and its potential function in vertebrates. In silico analyses showed that Gin-2 is a very ancient gene (500 My, only absent from placentals) derived from GIN transposons. Further analyses present a particular expression in brain and gonads during adulthood, while a strong expression during gastrulation suggests a potential role of Gin-2 in zebrafish development. All together, the different analyses contribute to a better view of TE evolution and their evolutionary impacts in vertebrate genomes.

L'application du procédé SLM est limitée par la difficulté à contrôler le procédé. Son application aux céramiques est particulièrement difficile en raison de leur faible absorption au laser et de leur faible résistance au choc thermique. La maîtrise de ce procédé nécessite une compréhension complète du transfert de chaleur, de la dynamique des fluides et de la mécanique des solides. Dans ce travail, nous proposons un modèle numérique pour la simulation du procédé SLM appliqué aux céramiques. Le modèle est développé à l'échelle du cordon et avec l'hypothèse d'un lit de poudre continu. Il est basé sur la méthode level set et l'homogénéisation multiphasique, avec laquelle nous sommes capables de suivre l'évolution de l'interface gaz/matière et les transformations de phase. La simulation dévelopée permet d'étudier l'influence des propriétés du matériau et des paramètres du procédé sur la température, la forme du bain liquide, la dynamique des fluides et la mécanique des solides. En dehors de la puissance du laser et de la vitesse de balayage, l'absorption du matériau est également importante pour la thermique et la forme du bain liquide. Avec la dynamique des fluides, la forme convexe du cordon est obtenue sous tension de surface. Les gouttelettes liquides se forment lors de la fusion de la poudre et créent une instabilité du bain. Ceci entraîne une irrégularité du cordon après solidification. L'effet Marangoni, provoqué par le gradient surfacique de la tension de surface, est étudié. Son influence sur la répartition de la température, la forme du bain liquide et la régularité du cordon est évoquée. Cet effet peut lisser la surface du cordon avec ∂γ/∂T négatif. En augmentant la vitesse de balayage, la surface du cordon devient plus irrégulière. L'effet de « balling » est reproduit avec une vitesse de balayage élevée. Cela peut être utile pour trouver le régime donnant une forme de cordon régulière étant données la puissance et la vitesse du laser. Le défaut de fissuration est délétère dans la fabrication additive. L'utilisation d'un laser auxiliaire peut aider à éviter ce défaut en diminuant la contrainte de traction maximale. Le mode de fonctionnement de ce laser auxiliaire reste un sujet intéressant à étudier et quelques pistes ont été données par les simulations présentées. Le modèle est validé par la comparaison de la forme du bain liquide avec des expériences dans différentes conditions de procédé. Les simulations peuvent également révéler la tendance de variation de la surface du cordon dans certains cas. Par la simulation de la déposition de cordons multiples, l'influence de taux de recouvrement sur la surface d'une couche, la température et l'évolution de contrainte est soulignée. / The application of SLM process is limited by the difficulty of process control. Its application to ceramics is especially challengeable due to their weak absorption to laser and weak resistance to thermal shock. The mastery of this process requires a full understanding of heat transfer, fluid dynamics in melt pool and solid mechanics. In this work, we propose a numerical model for the simulation of SLM process applied to ceramics. The model is developed at the track scale and with the assumption of continuous powder bed. It is based on level set method and multiphase homogenization, with which we are able to follow the evolution of gas/material interface and phase transformation. Simulations are performed to study the influence of material properties and process parameters on temperature, melt pool shape, fluid dynamics and solid mechanics. Apart from the laser power and scanning speed, material absorption is also found to be important to the thermal behavior and the melt pool shape. With the fluid dynamics, convex shape of track cross section is achieved under surface tension. Besides that, liquid droplets collapsing formed by the melting of powder create melt pool instability when falling, thus leading to track irregularity after solidification. The Marangoni effect, caused by surface tension gradient at gas/material interface, is investigated. Its influence on temperature distribution, melt pool shape and track regularity is recognized. One interesting finding is the smoothing effect of track surface with negative ∂γ/∂T. When combine surface tension with scanning speed, track surface becomes more irregular with the increase of scanning speed. The well-known balling effect is reproduced with high scanning speed. This can be helpful to find the regime for regular track shape with given laser power and scanning speed. Cracking defect is deleterious in additive manufacturing. The use of an auxiliary laser can help to avoid this defect by decreasing the maximum tensile stress. The process mode of this auxiliary laser remains an interesting subject to be studied and some guidelines have been given by the presented simulations. The model is validated by the comparison of melt pool shape with experiments under different process conditions. Simulations can also reveal the tendency of track surface variation for certain cases. By the application to multi-track deposition, the influence of hatch distance on layer surface, temperature and stress evolution is emphasized.