Les équations de Navier-Stokes linéarisées sont utilisées dans l’industrie aéronautique pour l’optimisation de forme aérodynamique, l’aéroélasticité et l’aéroacoustique. Deux axes ont été suivis pour accélérer et rendre plus robuste la résolution de ces équations. Le premier est l’amélioration de la méthode itérative de résolution de systèmes linéaires utilisée, et le deuxième la formulation du schéma numérique conduisant à ce système linéaire. Dans cette première partie, l’extension de l’algorithme GMRES avec déflation spectrale à des systèmes à plusieurs seconds membres a été testée sur des cas tests industriels. L’amélioration du préconditionnement de la méthode GMRES par l’utilisation d'une méthode de Schwarz additive avec préconditionneur ILU(k) a permis une accélération du temps de résolution allant jusqu’à un facteur dix, ainsi que la convergence de cas jusqu’alors impossibles à résoudre. La deuxième partie présente d’abord un travail sur la stabilisation SUPG du schéma élément fini utilisé. La forme proposée de la matrice de stabilisation, dite complète, a donné des résultats encourageants en non-linéaire qui ne se sont pas transposés en linéarisé. Une étude sur les conditions aux limites de Dirichlet clôt cette partie. Une méthode algébrique d’imposition de conditions non homogènes sur des variables non triviales du calcul, qui a permis l’application industrielle à l’aéroacoustique, y est détaillée. De plus, la preuve est apportée que le caractère transparent d’une condition de Dirichlet homogène sur toutes les variables s’explique par le schéma SUPG. / The linearized Navier-Stokes equations are solved at Dassault Aviation within numerical simulations for aerodynamic shape optimisation, flutter calculations and aeroacoustics. In order to improve the robustness and efficiency of the Navier-Stokes solver, this thesis followed two complementary paths. The first is work on the iterative methods used to solve linear systems, and the second is the improvement of the numerical scheme underlying these linear systems. In the first part, the extension to multiple right-hand sides of the GMRES algorithm with spectral deflation was tested on industrial test cases. The use of the ILU(k) preconditioner within an additive Schwarz method led to a reduction of the time needed to solve the systems with GMRES by a factor ten. It also enabled the convergence of some numerically very difficult cases which could not be solved by the software available before this thesis. The second part of the manuscript begins with work on the SUPG method used to stabilise the finite element scheme. A new way of computing the stabilisation matrix gave promising results on non-linear cases, which were however not observed for linear cases. A study on Dirichlet boundary conditions concludes this part. An algebraic method to impose non homogeneous Dirichlet boundary conditions on non-trivial variables is introduced. It enables the use, in an industrial context, of linearized Navier-Stokes for aeroacoustics. Moreover, the transparent behaviour of a homogeneous Dirichlet boundary conditions on all variables is proved to be due to the SUPG stabilisation.

Avec le développement des matériaux fonctionnels (multi-matériaux, multicouches,…), la caractérisation du comportement mécanique par des moyens macroscopiques conventionnels est devenue de plus en plus difficile. Ces méthodes conventionnelles sont donc substituées progressivement par des moyens de caractérisation multi-échelles. Parmi ces moyens, la nanoindentation, qui peut résoudre certains défis de la micro-caractérisation tels que la présence de phases indissociables, les systèmes multicouches, les revêtements ultra-minces, etc. Cet outil est devenu une technique de haute précision capable de solliciter des volumes de matière très faibles et fournir des informations riches pour la caractérisation des matériaux. Cependant, cet outil est utilisé majoritairement pour identifier les propriétés élastiques et qualitativement certains paramètres tels que la dureté, la ductilité et les contraintes internes.Ce travail de thèse s’intéresse à la caractérisation du comportement élastoplastique par nanoindentation à deux échelles : l’échelle macroscopique et l’échelle du cristal.Le premier défi de ce travail est expérimental. Il s’agit de générer des surfaces avec des propriétés représentatives de la microstructure étudiée. Ce défi est d’autant plus relevé que le matériau utilisé comme modèle est l’acier 316L très ductile et dont la surface est sensible au moindre changement. Un protocole expérimentale a été mis en place, à l’issu de ce travail, et les erreurs et dispersions de la réponse en nanoindentation introduites par les différentes étapes de génération de surface ont été quantifiés.Une base de données étendue a été mise en place, par la suite. Différentes géométries d’indent ont été appliquées à plusieurs profondeurs. Cette base de données va alimenter des stratégies d’identification inverse basée sur un couplage entre des algorithmes d’optimisation et une modélisation éléments finis de l’essai. Deux types d’algorithme ont été appliqués : Levenberg-Marquardt et l’algorithme génétique. Ce dernier est très consommateur en temps de calcul. Différents modèles EF axisymétrique et 3D ont été utilisés. Ces modèles ont été soigneusement optimisés par rapport au temps de calcul.Plusieurs stratégies d’identification ont été employées en se basant sur différentes données expérimentales issues de l’essai de nanoindentation telles que la courbe de charge-décharge, la forme de l’empreinte résiduelle et l’association de plusieurs géométries d’indent. Plusieurs modèles d’écrouissage isotrope ont été identifiés. À l’échelle macroscopique, les modèles d’écrouissage isotrope classiques ont été déterminés. À l’échelle du grain, la loi cristalline de Méric et Cailletaud a été identifiée. Les résultats obtenus ont été confrontés, à l’échelle macroscopique, à des identifications réalisées sur le même matériau à partir des essais de traction et de compression et ont montré que l’association de multiples géométries d’indentation permet de reproduire le comportement volumique du 316L avec une précision acceptable. Pour le comportement du cristal, des essais de compression de micropilliers ont été utilisé pour se procurer des données de référence à cette échelle. La comparaison montre beaucoup de dispersion dans les deux cas. En effet, certains phénomènes liés à la densité de dislocation très variables d’un grain à l’autre sont responsables de cette dispersion. Cette densité de dislocation n’est pas prise en compte, en tant que variable, dans le modèle cristallin utilisé. L’utilisation d’un modèle plus physique intégrant la densité de dislocation et son évolution permet d’améliorer ces résultats. Enfin, une nouvelle méthode d’identification a été proposée. Cette méthode est basée sur l’estimation et l’introduction de la géométrie réelle de l’indent dans le modèle EF utilisé pour l’identification. La méthode a été validée dans le cas de la pointe Berkovich et elle montre des résultats très prometteurs. / With the development of functional materials (multi-materials, multilayers, ...), the mechanical behavior characterization by conventional macroscopic methods has become progressively difficult. These conventional methods are therefore gradually substituted by multiscale characterization processes. Among these methods, the nanoindentation, this can solve certain challenges of micro-characterization such as the presence of indissociable phases, multilayer systems, ultra-thin coatings, etc. This tool has become a high-precision technique capable of testing very small volumes of matter and providing rich information for material characterization. However, this tool is used mainly to identify the elastic properties and, qualitatively, some parameters such as hardness, ductility and internal stresses.This thesis work focuses on the characterization of elastoplastic behavior by nanoindentation at two scales: the macroscopic scale and the crystal scale.The first challenge of this work is experimental. It involves generating surfaces with properties representative of the studied microstructure. This challenge is important because the material used as a model is 316L steel which is very ductile and whose surface is sensitive to small perturbations. An experimental protocol was implemented at the end of this work, and the errors and dispersions of the nanoindentation response introduced by the different surface generation steps were quantified. Then, a wide database was implemented with different indenter geometries and several depths. This database will feed inverse identification strategies based on a coupling between optimization algorithms and finite element modeling of this test. Two types of algorithm have been applied: Levenberg-Marquardt and genetic algorithms. The latter is very consumer in computing time. Different axisymmetric and 3D FE models have been used. These models have been carefully optimized with respect to computation time.Several identification strategies were employed based on various experimental databases from the nanoindentation test such as the loading-unloading curve, the residual imprint shape and the association of several indent geometries. Some models of isotropic hardening have been identified. On the macroscopic scale, classical isotropic hardening models have been determined. At the grain scale, the crystal plasticity constitutive model of Méric and Cailletaud has been identified. The results obtained were compared on the macroscopic scale with identifications carried out on the same material from the tensile and compression tests. The comparison showed that the combination of multiple indentation geometries makes it possible to reproduce the volume behavior of the 316L with acceptable accuracy. For crystal behavior, micropillar compression tests were used to obtain reference data at this scale. The comparison shows a lot of dispersion in both cases. Indeed, some phenomena related to the density of dislocation very variable from one grain to another are responsible of this dispersion. This dislocation density is not taken into account, as a variable, in the used crystal constitutive model. The use of a more physical law integrating the dislocation density and its evolution makes it possible to improve these results. Finally, a new identification method has been proposed. This method is based on estimating and introducing the real indent geometry in the FE model used for identification. The method has been validated in the case of Berkovich tip and shows very promising results.

Mycoplasma agalactiae est responsable de l'agalactie contagieuse, maladie des petits ruminants difficilement contrôlée et figurant sur la liste de l’OIE. Afin d’évaluer la diversité génétique de ce pathogène, 101 isolats ont été comparés par trois techniques (VNTR, RFLP, répertoire vpma). Les résultats révèlent une grande homogénéité génétique dont la souche type PG2 est représentative. Quelques isolats font exception telle la souche 5632 que nous avons séquencée et analysée ici. La comparaison des génomes et des protéomes entre 5632 et PG2 indiquent que la plasticité de ces génomes est liée à d’importants échanges d'ADN et à la présence de nombreux éléments génétiques mobiles (10% du génome). Ces analyses révèlent également une forte dynamique au sein de répertoires de gènes codant des protéines de surfaces. Pour les mycoplasmes, bactéries minimales dépourvues de paroi, ces évènements ont certainement joués un rôle dans leur survie et leur adaptation à des hôtes complexes. / Mycoplasma agalactiae is responsible of contagious agalactia, a disease of small ruminants that is still difficult to control and is listed by the OIE. In order to evaluate the genetic diversity of this pathogen, 101 isolates were compared using three techniques (VNTR, RFLP, vpma repertoire). Results revealed a high genetic homogeneity with the PG2 type strain as representative. Some isolates however diverged such as the 5632 which was sequenced and analysed here. Whole comparative genomic and proteomic analyses of the 5632 and PG2 strains indicate that their genomic plasticity resides in important genes flux and in the presence of several mobile genetic elements (10% of the genome). These analyses also revealed that specific loci encoding repertoire of surface proteins are highly dynamic. For these minimal bacteria that lack a cell-wall, these events have most likely played a major role in their survival and adaptation to complex hosts.

Le comportement non-linéaire des sols observé lors des mouvements sismiques forts est maintenant bien établi et le déploiement des puits accélérométriques a permis des analyses détaillés de la propagation des ondes ainsi qu'une évaluation quantitative des paramètres physiques tels que la vitesse de cisaillement et de compression des ondes et les facteurs d'amortissements en fonction de la déformation. En dépit du nombre grandissant d'´etudes sur ce phénomène, sa connaissance est encore récente et les recherches sur les données de puits accélérométriques restent une étape importante vers la compréhension du comportement complexe in-situ des sédiments soumis à des mouvements sismiques forts. L'objectif de ces travaux est triple. Premièrement, un code d'inversion par algorithme génétique est développé afin d'inverser des données de puits accélérométriques via la théorie des matrices de propagation de Thomson- Haskell. Cette technique nous permet dans un premier temps de valider la structure en une dimension (1D) (e.g., vitesse des ondes de cisaillement, facteurs d' amortissements) d'un puits accélérométrique dans le domaine linéaire et dans un second temps de mettre en évidence de manière quantitative le comportement non-linéaire des sédiments lors du séisme de Fukuoka, 2005, Japon. Deuxièmement, les résultats de l'inversion sont utilisés pour tester des lois de comportement simples et avancées en utilisant la Méthode des Eléments Finis. Les résultats montrent clairement que l'hypothèse bilinéaire de la loi de comportement simple produit des séries temporelles non réalistes en vitesse et en accélération. L'utilisation d'une loi de comportement avancée mène à de meilleurs résultats, cependant, le nombre de paramètres ajustables pour obtenir des résultats consistants avec l'observation est un obstacle inévitable. Troisièmement, afin d'étendre l'étude des effets de site à des dimensions supérieures, des codes 2D et 3D de la Méthode en Eléments Spectraux sont développés et validés en comparant leurs résultats dans le domaine linéaire avec ceux obtenus théoriquement ou via d'autres méthodes numériques.

Les propriétés physiques des milieux granulaires trouvent leur origine à l'échelle locale des contacts entre les grains. Cette étude est consacrée à l'analyse de l'influence de ces contacts sur le comportement global du matériau à différentes échelles. Nous développons dans la première partie de ce mémoire, une approche continue qui s'affranchit des limitations de la théorie de Janssen et permet de calculer les contraintes dans un matériau granulaire ensilé. Cette approche permet également de représenter aussi bien qualitativement que quantitativement l'effet d'écrantage dans les silos. Dans une seconde partie, une modélisation du phénomène d'écrantage est effectuée à partir des simulations numériques discrètes réalisées avec le code MULTICOR. Nous calculons les contraintes moyennes s'exerçant au niveau des parois lors du remplissage d'un silo par un milieu granulaire polydisperse. Une bonne concordance est observée entre les résultats des simulations numériques discrètes et ceux de l'approche continue développée dans la première partie. Une autre conséquence de l'existence de contacts privilégiés entre les grains est l'anisotropie caractéristique des matériaux granulaires. Dans la dernière partie, nous développons une approche micromécanique qui permet de modéliser cette anisotropie par l'intermédiaire d'un tenseur de texture d'ordre quatre. Nous proposons une hypothèse cinématique générale qui est incorporée dans un schéma d'homogénéisation des milieux granulaires anisotropes et comparée, dans le cas isotrope, à des simulations numériques discrètes existant dans la littérature. Les tendances qualitatives et quantitatives des résultats sont tout à fait satisfaisantes.

Le comportement mécanique de composites chanvre/époxy à renfort tissu a été étudié au moyen d'essais de traction aux différentes échelles des constituants. La variabilité des fibres naturelles a été prise en compte grâce à des analyses statistiques qui ont été intégrées dans un modèle analytique permettant d'obtenir une estimation des caractéristiques mécaniques des composites à partir de ses constituants élémentaires. Les mécanismes d'endommagement des composites ont été analysés par des observations microscopiques couplées à la technique d'émission acoustique et au suivi de la perte de rigidité. Les résultats ont montré une progression plus rapide de l'endommagement dans les composites chanvre/époxy [0/90]7 et [-45/+45]7 que dans leurs équivalents en composite verre/époxy pour un même taux volumique de fibres. De plus, les endommagements de l'interface fibre/matrice sont les plus nombreux et les plus précoces dans les composites chanvre/époxy. Les propriétés interfaciales et le transfert de charge à l'interface fibre/matrice des composites chanvre/époxy ont été étudiés grâce à des essais spécifiques instrumentés par photoélasticimétrie et corrélation d'images numériques. Les résultats de l'étude expérimentale ont montré une contrainte critique de cisaillement interfacial plus faible dans les composites chanvre/époxy que dans les composites verre/époxy. Les simulations numériques des tests de fragmentation ont fourni une estimation correcte de cette grandeur sur éprouvette monofilamentaire chanvre/époxy. Des tests d'optimisation ont également été réalisés en effectuant des traitements thermiques et chimiques des renforts.

Il a été récemment montré que le rapport B/Ca des foraminifères planctoniques dépend du pH des eaux de surface, tandis que celui des foraminifères benthiques dépend de la saturation en ions carbonates des eaux de fond. Ce nouveau traceur pourrait permettre de mieux contraindre les variations du cycle océanique du carbone, en particulier lors des changements glaciaire-interglaciaires (G/IG) du Pléistocène. Les objectifs de cette thèse étaient donc (1) de développer l'analyse du B/Ca des foraminifères par spectrométrie de masse (ICP-QMS) au LSCE ; (2) d'identifier certains biais potentiels de ce traceur pour en apprécier les limites; et (3) d'appliquer le B/Ca à des reconstructions paléocéanographiques. Les contaminations et la facilité avec laquelle le bore s'adsorbe (effet mémoire) rendent difficile son analyse dans les coquilles de foraminifères où le rapport B/Ca est faible, typiquement compris autour de ~ 40-120 µmol/mol pour les foraminifères planctoniques et ~ 100-250 µmol/mol pour les foraminifères benthiques. Les techniques de nettoyage des foraminifères et d'analyse par ICP-QMS que j'ai employées et améliorées permettent au final de mesurer le B/Ca des foraminifères avec une précision de l'ordre de 3.5%.A partir de l'étude de sédiments de surface prélevés le long d'un profil bathymétrique dans l'Atlantique (Sierra Leone Rise), j'ai pu montrer que la dissolution des tests de l'espèce planctonique Globigerinoides sacculifer se traduit par une perte préférentielle du bore. La dissolution croissante des coquilles de G. sacculifer le long du profil se traduit ainsi par une diminution de ~ 0.1 unités dans les estimations du pH des eaux de surface à partir du B/Ca, un ordre de grandeur similaire aux variations G/IG. Pour s 'affranchir de cet effet de la dissolution, je propose une procédure de correction basée sur le B/Ca des foraminifères benthiques (ΔCO₃²⁻ des eaux de fond). Cependant, au-delà de l'effet de la dissolution, la reconstruction du paléo-pH des eaux de surface se heurte surtout aux contradictions des études récentes concernant la calibration du coefficient apparent de partition (KD) du bore dans la calcite des foraminifères planctoniques. Pour s'affranchir de la co-variance de la température, de la salinité et du CO₃²⁻ dans les eaux de surface Atlantique, qui peut masquer le(s) facteur(s) important(s) pour la calibration du KD, j'ai réalisé l'analyse du B/Ca de l'espèce Globigerinoides ruber sur des sommets de carottes de l'Atlantique et de l'Indo-Pacifique. Mes données suggèrent qu'il n'y a pas de relation robuste entre le KD et la température ou le CO₃²⁻ des eaux de surface. Pour souligner les incertitudes sur les reconstructions de paléo-pH depuis le dernier stade glaciaire, j'ai effectué une étude de sensibilité basée sur les données B/Ca de l'espèce G. ruber prélevée dans une carotte du Canal du Mozambique et je propose plusieurs pistes pour tenter d'améliorer, à terme, la compréhension du proxy B/Ca dans les foraminifères planctoniques.Si l'utilisation du B/Ca des foraminifères planctoniques est complexe et demande encore un gros travail de compréhension et de calibration, celle du B/Ca des foraminifères benthiques est plus simple, le B/Ca étant directement lié à la saturation en ions carbonates des eaux de fond. L'analyse du B/Ca de Cibicidoides wuellestorfi à travers deux transitions G/IG (MIS16/15 et MIS12/11) au site ODP 849 (Pacifique équatorial) m'a permis de montrer que les changements de la chimie des eaux de fond estimés à partir de traceurs sédimentaires de la dissolution (e.g. teneur en CaCO3, poids et fragmentation des foraminifères,..) ont des amplitudes erronées et présentent des biais temporels. Les profils de B/Ca et du d13C des foraminifères benthiques permettent de mettre en évidence i) le lien entre dégazage océanique de CO2 et préservation des carbonates et ii) la dynamique de la compensation des carbonates sédimentaires marins.

Au cours de cette thèse, nous avons développé un outil numérique, basé sur une formulation intégrale en éléments de frontière, qui permet une analyse classique du comportement d'une fissure 3D soumise à des sollicitations mécaniques complexes. Cet outil industriel est destiné à être intégré dans un code de calcul à usage industriel. Dans le but d'appréhender l'impact de la microstructure sur le comportement en fissuration fragile, nous nous sommes intéressés aux modèles de comportement non local. Nous avons commencé par adopter le modèle de comportement élastique non local de Eringen, qui permet de décrire plus finement le comportement élastique au voisinage de la fissure en prenant en compte les interactions à longue distance au sein du matériau. Cette modélisation du comportement conduit, contrairement à l'approche classique, à un un champ de contrainte fini sur le front de la fissure et localement maximal en avant du front. Ces résultats montrent qu'il est possible de prévoir la stabilité et la direction de propagation de la fissure à l'aide d'un critère plus simple et plus naturel, basé sur les variations du champ de contrainte au voisinage du front de la fissure. La stratégie numérique adoptée permet de traiter indifféremment des cas de fissure en traction, compression, cisaillement ou sollicitation mixte. L'intérêt de l'approche non-locale étant clairement démontré, nous avons considéré la version améliorée du modèle de Eringen telle que proposée par Polizzotto. Cette modélisation est la plus appropriée pour les milieux finis et requiert une mise en oeuvre numérique particulière. Les bases d'une méthodologie numérique, initiée par R. Kouitat ont été établies. Cette méthode est fondée sur un couplage des éléments de frontière avec une méthode de collocation par points d'équations aux dérivées partielles. Les premiers résultats obtenus dans ce cadre sont très encourageants et montrent qu'il sera effectivement possible de traiter le phénomène irréversible de fissuration de la même façon que les problèmes de plasticité / In this thesis, we have developed a numerical tool, based on a classical boundary elements method, which allows a conventional analysis of a stationary crack in a 3D specimen under complex mechanical loading. In order to assess the impact of the microstructure on the brittle fracture, we were interested in non local models of behavior. First, we have adopted the non local elastic model due to Eringen. This refined constitutive equation allows to account for long range interactions in the description of the elastic behavior in the vicinity of the crack front. Unlike the traditional approach, this type of model leads to a finite stress field at the crack front. Moreover, the stress is locally maximal ahead of the front. These interesting results indicate that it is possible to predict the stability and direction of crack propagation in a simple and more naturel way by using stress based criteria. The implemented numerical strategy can handle cases of crack in tension or compression, under shear stress or mixed loadings. Having clearly highlighted the interest of non local models, we have adopted the improved version of Eringen elastic model as proposed by Polizzotto. This elastic model is applicable to finite domains and requires a specific numerical treatment. The basis of such a numerical strategy initiated by R. Kouitat has been established. The method couples the conventional boundary element method with local point interpolation of a strong form differential equation. Promising results are obtained and show that with such modeling of material behavior, it is possible to describe the irreversible process of fracturing in a similar way as plasticity

Les chaussées subissent des sollicitations liées au trafic et au climat conduisant à leur dégradation, par fissuration notamment. Il est nécessaire dans le contexte actuel de pouvoir modéliser le comportement de ces structures multicouches endommagées afin de prévoir leur durée de vie résiduelle ou dimensionner des solutions de renforcement. L’objectif de la thèse est ainsi de proposer un outil de calcul dédié à l'analyse 3D des chaussées fissurées ou comportant des discontinuités. L’approche retenue repose sur la modélisation simplifiée d’une chaussée par un empilement de plaques du Modèle Multi-particulaire des Matériaux Multicouches à 5n équations d’équilibre (M4-5n). Un outil de calcul rapide de référence de chaussées 2D fissurées et une méthode de résolution générale du M4-5n par Eléments Finis mixtes sont développés. Le point de départ de la méthode de résolution est l’écriture, pour le M4-5n, du principe variationnel basé sur le théorème de l'énergie complémentaire où la condition de contraintes statiquement admissibles est assurée à partir de multiplicateurs de Lagrange. La discrétisation des efforts généralisés utilise des espaces d’interpolation permettant le bon conditionnement du système d’équations algébriques à résoudre et garantissant la stabilité de la solution. La méthode est implémentée dans FreeFem++. Elle ramène le problème 3D initial à une modélisation EF 2D et conduit à des valeurs finies des efforts généralisés au niveau des fissures ou décollement d’interface. L’outil de calcul final ainsi développé est validé et appliqué à l’étude de la réponse d’une structure fissurée,représentative d’une chaussée testée en vraie grandeur sur le site de l’IFSTTAR. / Pavements are multilayer structures which undergo cracking distress due to traffic and climatic factors. It is important nowadays to be able to model the mechanical response of such damaged pavements in order to assess their residual lifetime or to design reinforcement solutions. In this context, the present thesis aims at developing a numerical tool dedicated to the analysis of pavements incorporating cracks or discontinuities. In the developed approach, the pavement structure is modeled as a stacking of “plate” elements of typeM4-5n (Multi-Particle Models of Multilayer Materials) which considers 5n equilibrium equations. A reference quick 2D calculation tool for cracked pavements and a general solving of M4-5n by the mixed Finite Element (FE) method was developed. The starting point for this method is the derivation for M4-5n of the variational principle based on the complementary energy theorem whose condition of statically admissible stress is taken into account using Lagrange multipliers. Discretization of the generalized stresses involves interpolation spaces, proposed to avoid ill-conditioned system of algebraic equations after discretization and to insure stability of the solution. The developed method is implemented in a FreeFem++ script. In this method, the initial 3D problem can be handled through FE simulations in 2D and finite values of the generalized stresses are obtained at crack and interlayer debonding locations. The developed numerical tool was validated and applied to the study of the mechanical response of a structure with cracks representative of a pavement tested underfull-scale conditions during an accelerated fatigue test performed at IFSTTAR.

La macroségrégation est un défaut connu dans les procédés de coulées industrielles. La genèse de ce défaut est la conséquence de l'interaction complexe entre la microségrégation ou la distribution des espèces chimiques à l'échelle de la microstructure et les mouvements des phases liquide et solides. Les hétérogéneités de concentration en solutés à l'échelle de la pièce peuvent être rédhibitoires vis-à-vis de la qualité du produit.Dans ce travail, on propose un modèle numérique pour simuler et prédire la formation des macroségrégations en coeur des pièces d'alliages multi-constitués, induites par des variations thermiques et solutales dans la phase liquide. Dans un premier temps, on considère que le métal solidifie à volume constant. Dans ce contexte, la convection thermosolutale est étudiée ainsi que son influence sur la formation des canaux ségrégées à différentes échelles de modélisation. Dans un deuxième temps, le modèle vise à prédire les macroségrégations en présence de changement de volume du métal, dont la cause principale est le retrait à la solidification, pouvant être à l'origine du phénomène de ségrégation inverse. La surface entre le métal et le gaz environnant au cours du retrait évolue pendant le retrait en fonction du chemin de solidification qui varie avec la macroségrégation. Cette évolution d'interface est suivi par la méthode level set. Des prédictions de concentration moyenne, couplées aux bases de données thermodynamiques pour mieux prédire les chemins de solidification des alliages multi-constitués, sont analysées et comparées avec des résultats expérimentaux. Finalement, des calculs de solidification en microgravité sont présentées, simulant un essai expérimental dans le contexte du projet CCEMLCC lancé par l'Agence Spatiale Européenne. Les résultats en fin de solidification montrent un accord acceptable quant à la forme et l'élongation des échantillons solidifiés. Ces calculs sont faits avec des approximations binaire, ternaire et quaternaire d'une même nuance d'acier utilisée dans les essais en microgravité. / Macrosegregation is key defect in industrial casting processes. During solidification, solute redistribution at the scale of microstructure, also known as microsegregation, take place with complex interactions, in order to form one or more solid phases. These interactions between microsegregation and movements of liquid and solid phases may lead to macrosegregations. These solute heterogeneities spanning on a larger scale, may result in a bad casting quality. In this thesis, we propose a numerical model to simulate and predict macrosegregations occurring in the centre of multicomponent alloys, caused by thermal and solutal variations in the liquid phase. First, we assume that the metallic alloy solidifies with a constant volume. In this context, we study the influence of thermosolutal convection on the formation of channel segregations, at different modelling scales. The second part of this modelling work consider solidification while the metallic alloy's volume is decreasing, mainly due to overall density variation, also known as solidification shrinkage, possibly leading to the so-called inverse segregation phenomenon, appearing on the alloy's skin. In the context of solidification shrinkage, the shape of the metal's boundary with surrounding gases varies according to a constantly changing solidification path due to macrosegregation. The level set method is therefore used to track its evolution with time. Composition predictions, coupled with thermodynamic database mappings for more accurate multicomponent solidification paths, are analysed and compared to existing experimental setups. Finally, simulations of a reduced-gravity solidification cases are performed, mocking an experimental benchmark from the CCEMLCC project launched by the European Space Agency. The results after complete solidification show acceptable agreement for the final shape, compared to experimental results. These computations were performed with binary, ternary and quaternary approximations of the same steel grade which was used in reduced-gravity experiments.