Les simulations numériques à l'échelle du pore sont fréquemment utilisées pour étudier le comportement des écoulements multiphasiques largement rencontrées dans phénomènes naturels et applications industrielles. Dans ce travail, la morphologie de structures liquides et l'action capillaire sont examinées à l'échelle des pores par la méthode de Boltzmann sur réseau (LBM) à plusieurs composants selon le modèle de Shan-Chen. Les résultats numériques obtenus sont en bon accord avec les solutions théoriques. Les simulations numériques sont étendues à microstructures complexes au-delà du régime pendulaire.La LBM a été utilisée pour modéliser l'écoulement multiphasique à travers un milieu poreux idéalisé dans des conditions de drainage primaire quasi-statique. Les simulations LBM ont fourni une excellente description du déplacement de l'interface fluide-fluide à travers les grains. Pendant le drainage, les simulations LBM sont capables de reproduire la déconnexion d'une phase dans le milieu granulaire sous la forme de ponts pendulaires ou structures liquides complexes. Malheureusement, le temps de calcul nécessaire pour ce type de simulations est assez élevé. Afin d’optimiser les ressources de calcul, nous présentons un modèle 2D (modèle Throat-Network) basé sur des solutions analytiques pour décrire l'écoulement biphasique à travers un ensemble de disques dans un temps de calcul très réduit, donc le modèle 2D est susceptible de remplacer les simulations LBM lorsque les ressources de calcul sont limitées. L'approche souligne l'importance de simuler le problème a l'échelle de la gorge du pore pour obtenir les relations volume - pression capillaire locales. Le modèle Throat-Network est un point de départ pour le modèle hybride proposé pour résoudre les problèmes en 3D. Le modèle hybride combine l’efficacité de l’approche réseau de pores et la précision du LBM à l’échelle des pores. Le modèle hybride est basé sur la décomposition de l’échantillon en petits sous-domaines, dans lesquels des simulations LBM sont effectuées pour déterminer les propriétés hydrostatiques principales (pression capillaire d'entrée, courbe de drainage primaire et morphologie du liquide pour chaque gorge du pore). Malgré la réduction significative des temps de calcul obtenus avec le modèle hybride, le temps n’est pas négligeable et les modélisations numériques d'échantillons de grandes tailles ne sont pas réalistes. Les approximations données par les méthodes Incircle et MS-P, qui prédisent les propriétés hydrostatiques, sont comparées à celles de LBM et du modèle hybride. / Numerical simulations at the pore scale are a way to study the behavior of multiphase flows encountered in many natural processes and industrial applications. In this work, liquid morphology and capillary action are examined at the pore-scale by means of the multicomponent Shan-Chen lattice Boltzmann method (LBM). The accuracy of the numerical model is first contrasted with theoretical solutions. The numerical results are extended to complex microstructures beyond the pendular regime.The LBM has been employed to simulate multiphase flow through idealized granular porous media under quasi-static primary drainage conditions. LBM simulations provide an excellent description of the fluid-fluid interface displacement through the grains. Additionally, the receding phase trapped in the granular media in form of pendular bridges or liquid clusters is well captured. Unfortunately, such simulations require a significant computation time. A 2D model (Throat-Network model) based on analytical solutions is proposed to mimic the multiphase flow with very reduced computation cost, therefore, suitable to replace LBM simulations when the computation resources are limited. The approach emphasizes the importance of simulating at the throat scale rather than the pore body scale in order to obtain the local capillary pressure - liquid content relationships. The Throat-Network model is a starting point for the a hybrid model proposed to solve 3D problems. The hybrid model combines the efficiency of the pore-network approach and the accuracy of the LBM at the pore scale to optimize the computational resources. The hybrid model is based on the decomposition of the granular assembly into small subsets, in which LBM simulations are performed to determine the main hydrostatic properties (entry capillary pressure, capillary pressure - liquid content relationship and liquid morphology for each pore throat). Despite the reduction of computation time, it is still not negligible and not affordable for large granular packings. Approximations by the Incircle and the MS-P method, which predict hydrostatic properties, are contrasted with the results provided by LBM and the hybrid model. Relatively accurate predictions are given by the approximations.

Une synthèse bibliographique a permis de faire un bilan sur les mécanismes d'échanges entre phases et de transport d'un composé organo-volatil présent en zone non saturée, de mettre en équations son transfert et de définir les dispositifs expérimentaux d'étude ainsi que les essais pertinents.<br /><br />Le polluant considéré est le trichloréthylène, le milieu poreux est un sable moyen et les expériences ont été réalisées en colonne.<br /><br />Après la mise au point d'une méthode d'analyse du TCE aqueux, gazeux et sorbé, les mécanismes prédominants des transferts ont été quantifiés isolement essentiellement par le biais d'essais de diffusion en sable sec, à saturation résiduelle en eau et saturé. Puis, ces mécanismes ont été couplés au sein d'une expérience de transfert de TCE en sable avec un profil de teneur en eau réel (hydrostatique). Chaque type d'essai a été dimensionné, si possible dupliqué et interprété grâce au code multiphasique Comsol dont l'équation de l'écoulement a été modifiée afin de considérer la convection gravitaire.<br /><br />En contrôlant rigoureusement les facteurs externes influents ainsi que les conditions aux limites, il a été possible de réaliser des expériences reproductibles et interprétables de transfert de TCE, un composé volatil et réactif en milieu poreux très perméable. Une bonne reproductibilité des résultats expérimentaux par la simulation a été obtenue moyennant de légères modifications des paramètres de base : rapport perméabilité sur viscosité, tortuosité(Millington, 1959) et paramètre de calage de la courbe de conductivité aéraulique (Thomson et al., 1997). Ce travail a permis une compréhension fine des transferts gazeux en ZNS, en particulier en frange capillaire.

Un bassin sédimentaire est un milieu poreux de grande dimension (plusieurs dizaines de kilomètres de long et de large pour une profondeur d'environ cinq kilomètres) qui évolue au cours du temps par les effets de compaction et de sédimentation. Au cours de cette évolution, des hydrocarbures vont se former et s'écouler dans le bassin. On établit alors un modèle permettant de simuler cette évolution de bassin ainsi que la création, la migration et le piégeage des hydrocarbures dans des roches appelées roches magasins. Ces phénomènes se déroulant sur des centaines de millions d'années, on s'est attaché à étudier principalement une discrétisation temporelle de ces équations. On a ainsi mis en avant un raffinement local du pas de temps dont le principe est de recalculer la solution sur une zone jugée "mauvaise". A l'extérieure de cette zone, la solution est admissible. La difficulté vient de la détermination de la zone qui doit - être suffisamment "grande" pour avoir une bonne qualité de la solution, mais suffisamment "petite" pour obtenir un gain calcul. Les estimateurs a posteriori permettent de contourner cette difficulté. On a donc entrepris une étude théorique de ces estimateurs a posteriori dans le cas des équations linéaires elliptique et parabolique. Des simulations numériques montrent l'efficacité de ces estimateurs dans des cas académiques.

La simulation numérique du comportement et l'analyse de stabilité des ouvrages de géotechnique renforcés par inclusions linéaires "rigides" ou "souples", c'est-à-dire capables ou non de reprendre des efforts de flexion et cisaillement, reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions employées, ainsi qu'à leurs dimensions caractéristiques très petites vis à vis de celles de l'ouvrage dans son ensemble. Une approche alternative, fondée sur une modélisation dite "multiphasique", a été récemment développée, qui permet de s'affranchir des difficultés précédentes et d'aboutir ainsi à des méthodes de dimensionnement accessibles à l'ingénieur. On s'intéresse dans le cadre de ce travail aux développements théorique et numérique de cette approche en prenant en compte non seulement les effets de flexion et de cisaillement des inclusions, mais également les différents types d'interaction entre ces dernières et le sol environnant. La première partie du mémoire est consacrée à la construction des approches du calcul à la rupture et à la mise en œuvre numérique élasto-plastique du modèle multiphasique des ouvrages renforcés par inclusions "souples", qui est alors appliqué à l'analyse de la stabilité d'ouvrage en "terre armée". L'extension de l'approche du calcul à la rupture et du code éléments finis en élastoplasticité, développés dans la première partie, à la prise en compte des effets de flexion et de cisaillement fait l'objet de la seconde partie du mémoire. La version la plus complète du modèle multiphasique élaborée est appliquée à la simulation du comportement d'une fondation renforcée par inclusions "rigides", soumise à un séisme induisant un chargement latéral. Partant de l'idée d'optimiser le schéma de renforcement de l'ouvrage traité à la deuxième partie, en inclinant les inclusions, un modèle "triphasique" est formulé et développé. Ce modèle vise à décrire le comportement macroscopique d'un sol renforcé par deux familles d'inclusions ayant des orientations différentes.

L'objectif du travail de thèse a été la mise en place de deux modélisations multi-physiques fondées sur des fondements théoriques communs mais appliquées à deux domaines de la recherche scientifique très différents: i) l'étude du comportement du béton au jeune âge pour la prévention de la fissuration précoce; ii) l'analyse des phénomènes physiques, chimiques et biologiques qui gouvernent la croissance et l'évolution de la tumeur cancéreuse. Le développement d'un outil numérique pour la modélisation du béton au jeune âge est très important pour la conception de structures durables. Le modèle développé pendant la thèse doctorale a été implanté sur le code aux éléments finis Cast3M, puis validé expérimentalement. Il permet de multiples applications: étude des sollicitations et des phénomènes de fissuration au jeune âge, gradients thermiques et hydriques, prédiction du retrait endogène et de dessiccation, étude de l'inhibition de l'hydratation causée par le séchage, prédiction du fluage et de la redistribution des contraintes associées, étude des réparations. Les équations qui gouvernent le comportement thermo-hydro-chemo-mécanique du béton au jeune âge ont plusieurs analogies formelles avec celles qui sont typiquement à la base de la modélisation de la croissance des tumeurs cancéreuses. L'élargissement de l'analyse numérique dans le domaine médical est d'un grand intérêt social en complément de l'intérêt scientifique. Les équations utilisées pour le béton ont été réadaptées, et le modèle mathématique obtenu a été implanté dans Cast3M. Les premiers résultats du modèle ont été satisfaisants et qualitativement très proches des données expérimentales de la littérature dans ce domaine.

La technique de renforcement des sols par inclusions linéaires rigides, c'est-à-dire capables de reprendre des efforts de flexion et cisaillement, connaît depuis quelques décennies un essor important bien que la simulation et le calcul de tels ouvrages reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions mis en jeu dans de tels ouvrages. Une modélisation multi phasique a récemment été développée permettant d'aborder le dimensionnement de tels ouvrages de géotechnique avec des temps de calcul considérablement réduits par rapport à ceux d'une modélisation numérique classique, où les inclusions de renforcement doivent être discrétisées individuellement. On s'intéresse dans le cadre de ce travail à la mise en oeuvre de ce modèle dans le cadre du comportement élastoplastique des différents constituants, lorsque les effets de flexion et de cisaillement des inclusions sont pris en compte. Après avoir développé un certain nombre de solutions analytiques de référence, un outil numérique basé sur la méthode des éléments finis est mis au point, puis appliqué à la simulation du comportement de différentes structures renforcées, telles que les radiers de fondation sur groupes de pieux, ou la stabilisation de pentes instables par l'incorporation de telles inclusions.

Les Champs Phlégréens, situés dans la métropole napolitaine (Italie du sud), forment l’une des plus grandes structures volcaniques au monde. Depuis 1950, ce complexe volcanique manifeste un regain d’activité, qui s’est amplifié au cours de la dernière décennie. Cette accélération s’exprime au travers d’une intensification de la sismicité, de la déformation du sol ainsi qu’une extension de la zone de dégazage. L’ensemble des récentes études s’accorde à dire que le système s’achemine actuellement vers un point critique, sans toutefois pouvoir préciser quand et où pourrait avoir lieu une éventuelle éruption. Cette difficulté à prédire l’état réel du système est principalement associée à la présence d’un système hydrothermal relativement développé. Aux Champs Phlégréens, il est en effet difficile de déconvoluer les signaux provenant du forçage magmatique de ceux résultant de la réponse hydrothermale. L’objectif de cette thèse est donc d’améliorer les connaissances actuelles du système hydrothermal superficiel du volcan de la Solfatara, lieu où se concentre actuellement la reprise d’activité. Pour cela, une approche multidisciplinaire a été menée en deux phases : l’imagerie géophysique du volcan puis la modélisation de son système hydrothermal.La tomographie haute-résolution de résistivité électrique 3-D du cratère a permis de reconnaître les principales formations géologiques et leurs connexions avec les structures et écoulements hydrothermaux. L’interprétation du modèle de résistivité électrique a été réalisée grâce à un ensemble de mesures superficielles complémentaires : flux de CO2, température, potentiel spontané, capacité d’échange cationique et pH du sol. Deux panaches à dominante liquide ont été identifiés : la mare de boue de la Fangaia et la fumerole de Pisciarelli. À la Fangaia, une étude conjointe des modèles de résistivité électrique et de vitesses du sous-sol (obtenues par l’INGV) établit la présence de forts gradients, à la frontière entre panache hydrothermal et zone de dégazage diffus. Au niveau du principal secteur fumerolien, le modèle de résistivité électrique et la localisation des sources acoustiques révèlent clairement l’anatomie d’une zone fumerolienne. Deux conduits séparés, saturés en gaz, alimentent les fumeroles de Bocca Grande et de Bocca Nuova, depuis un même réservoir de gaz situé à ~50 mètres de profondeur. L’intense dégazage diffus produit à proximité de ces fumeroles occasionne la condensation de vapeur. Le modèle de résistivité électrique met en évidence la circulation souterraine de cet important volume d’eau, canalisée à l’intérieur d’une zone fracturée.En utilisant l’ensemble de ces informations structurelles, un modèle thermodynamique des écoulements multiphasiques de la principale zone fumerolienne a été réalisé. Ce modèle reproduit fidèlement les observables des fumeroles : température, flux et rapport CO2/H2O. Il valide l’imagerie géophysique et confirme l’interaction entre la circulation d’eau de condensation et l’un des conduits fumeroliens. Ainsi, cette simulation explique, pour la première fois par un effet d’interaction superficiel, les différentes signatures géochimiques des deux fumeroles : Bocca Nuova et Bocca Grande. L’approche multidisciplinaire, employée dans cette thèse, constitue une nouvelle étape vers une meilleure connaissance des interactions hydrothermales. Celles-ci doivent être prise en compte dans l’objectif de réaliser des modélisations dynamiques précises permettant d’appréhender in fine l’état réel du système volcanique. / The Campi Flegrei caldera is located in the metropolitan area of Naples (Italy), and it is one of the largest volcanic systems on Earth. Since 1950, this volcanic complex shows significant unrest, which accelerated over the last decade with a rise in the seismic activity, ground deformation, and the extent of the degassing area. Recent studies indicate that the volcanic system is potentially moving toward a critical state, although their authors remain unable to point out when and where a possible eruption could take place. The difficulty of predicting the real volcanic state is here mainly related to the hydrothermal system. Indeed, at the Campi Flegrei, it is difficult to separate the magmatic input signal from the hydrothermal response. Hence, the aim of this thesis is to improve our knowledge on the shallow hydrothermal system of the Solfatara volcano, where most of the renewal activity takes place. A multidisciplinary approach has been performed in two steps: first a geophysical imagery of the volcano and second the modeling of its hydrothermal system.The 3-D electrical resistivity tomography of the crater allows to recognize the main geological units, and their connection with hydrothermal fluid flow features. The interpretation of the resistivity model has been realized thanks to numerous soil complementary measurements: CO2 flux, temperature, self-potential, Cation Exchange Capacity and pH. We identify two liquid-dominated plumes: the Fangaia mud pool and the Pisciarelli fumarole. In the Fangaia area, the comparison between electrical resistivity and velocity models reveals strong gradients related to a sharp transition at the border between the hydrothermal plume and the high diffuse degassing region. Combining electrical resistivity model with hydrothermal tremor sources localization reveal the anatomy of the main fumarolic area. Two separated conduits, gas-saturated, feed the two fumaroles Bocca Grande and Bocca Nuova. These conduits originate from the same gas reservoir located 60 m below the surface. The intense degassing activity, produced in the vicinity of fumaroles, creates large amounts of vapor condensation. The resistivity model reveals this condensate circulation, within a fractured area.All these results are incorporated into a multiphase flow model of the main fumarolic area. The simulation accurately reproduces the fumaroles observables: temperature, flux and CO2/H2O ratio. The model validates the geophysical imagery and confirms the interaction between Bocca Nuova fumarolic conduit and the condensate flow. Hence, this simulation explains for the first time the distinct geochemical signature of the two fumaroles due to a shallow water-interaction. The multidisciplinary approach performed in this thesis constitutes a new step toward a better understanding of hydrothermal interactions. Those phenomena have to be taken into account in order to perform dynamic modelling, and thus apprehend the real state of the volcanic system.

Depuis le début du millénaire, la simulation numérique directe est apparue comme un outil précieux capable d'étudier l’atomisation d’une goutte isolée par un écoulement à grande vitesse. L’atomisation peut être divisée en deux phases distinctes : l'éclatement se produit d'abord sous la forme d'aplatissement de la goutte, formant également des filaments, puis il se poursuit via l'obtention d'une multitude de gouttes de tailles réduites ce qui complète le processus d’atomisation. Les principaux objectifs pour le présent travail étaient donc d’établir un modèle et une méthode numérique capables d’étudier au mieux ces phénomènes. L'atomisation d’une goutte isolée est présentée et est accompagnée d’une comparaison avec l’expérience qui confirme les capacités du modèle et de la méthode à simuler numériquement les différents processus physiques mis en jeu. Des informations essentielles quant aux mécanismes d’atomisation, non exploitables avec l’expérience, sont décrites et l’objectif d’obtenir des gouttes de tailles réduites est atteint. / Only at the beginning of the millennium, direct numerical simulation has emerged as a valuable tool capable of studying the atomization of an isolated droplet by a high-speed flow. The atomization can be divided into two distinct phases: the aerobreakup occurs first in the form of flattening of the droplet, also forming filaments, and then it continues via the obtaining of a multitude of reduced sizes droplets what completes the process of atomization. The main objectives of this work were therefore to establish a model and a numerical method able to study these phenomena as well as possible. The atomization of an isolated droplet is presented and is accompanied by a comparison with the experiment which confirms the capacities of the model and the method to numerically simulate the different physical processes involved. Essential information on atomization mechanisms, which cannot be exploited with experiments, is described and the objective of obtaining droplets of reduced sizes is achieved.

Dans les travaux de recherche développés durant cette thèse, une première approche desécoulements magnétohydrodynamiques (MHD) multiphasiques est proposée. Cette approche seconcentre sur les phénomènes liés à l’interaction entre une dynamique interfaciale et un écoulementde coeur MHD. Le couplage induit entre rhéologie de surface etMHDde sous-phase est particulièrementillustré par le développement d’un viscosimètre annulaire surfacique, dédié à l’étude desmétaux liquides progressivement oxydés. En premier lieu sont introduits les éléments théoriquespropres à laMHDet à la rhéologie de surface. La modélisation de leur couplage fait intervenir deuxparamètres interfaciaux : les viscosités surfaciques dilatationnelle et de cisaillement. L’influence respectivede ces deux paramètres sur l’écoulement MHD de sous-phase est étudiée analytiquementet numériquement dans le cas d’un écoulement (stratifié) annulaire MHD permanent. Dans la configurationretenue, un champ magnétique uniforme vertical est imposé, perpendiculairement à lasurface liquide graduellement oxydée. Le rôle décisif des contraintes visqueuses interfaciales concernantl’(in)activation des couches de Hartmann est démontré, conduisant à des topologies atypiquesd’écoulement MHD. Le viscosimètre annulaire MHD est ensuite proposé en tant que méthodeexpérimentale originale, permettant la mesure sélective des viscosités surfaciques de fluides électroconducteurs.Les premières campagnes expérimentales aboutissent à une estimation de la viscositéde cisaillement interfacial d’un alliage métallique (GaInSn) liquide à température ambiante. Finalement,une ouverture sur l’écoulement MHD à proximité d’une inclusion gazeuse sphérique rigideest discutée en annexe de ce projet, en lien avec des conditions mécaniques variables à l’interfaceliquide/gaz pilotées par la rhéologie de surface, constituant une première approche vers la descriptiondes écoulements MHD dispersés. / In this thesis work, a first approach towards the description of magnetohydrodynamic(MHD) multiphase flows is proposed, based on the investigation of the role of surface rheologyin permanent regime. The study of the coupling between bulk MHD and surface rheology is particularlymotivated by the development of an annular surface viscometer devoted to liquid metalstopped with an oxidation layer. First, theoretical foundations of MHD and of surface rheologyare separately introduced. The modelling brings out a strong coupling between bulk and surfacevelocities, the latter being dependent on two interfacial parameters: the surface shear and dilatationalviscosities. Their respective influence is analytically and numerically investigated in the caseof a permanent annular (stratified) MHD end-driven flow. In the considered geometry, a vertical(uniform) magnetic field is imposed, perpendicular to the gradually oxidising liquid surface. Thecontribution of planar surface viscous stresses to the possible electrical activation of Hartmannlayers is demonstrated, leading to a variety of atypical MHD flow patterns. The annular MHDviscometer is then developed as a first device able to perform selective measurement of the surfaceviscosities of electroconductive fluids. First experiments lead to an estimated value of the surfaceshear viscosity for a liquid alloy (GaInSn) at room temperature. Some hints are finally given to investigatethe study case of aMHDflow past a rigid fluid sphere, with varying interfacial conditionsgoverned by surface rheology, as an additionalwork constituting a first step towards the descriptionof dispersed MHD flows.

La présente étude est dédiée à la compréhension des mécanismes de malaxage intervenant lors du soudage de matériaux dissimilaires par des sources de haute énergie et en particulier sur deux couples de matériaux présentant des problèmes métallurgiques différents : • cuivre - inox (lacune de miscibilité, différence de propriétés thermophysiques),• TA6V- inox (oxydation, formation de phases intermétalliques fragilisant la soudure).Pour le premier couple de matériaux, le soudage par laser Nd:YAG continu et par faisceau d'électrons a été utilisé. L'étude des évolutions de la morphologie des soudures, de la composition et de la microstructure des zones fondues ainsi que des propriétés mécaniques a permis de proposer des hypothèses sur les mécanismes de formation du mélange hétérogène à solubilité limitée. Afin de quantifier les phénomènes physiques intervenant en soudage continu de matériaux dissimilaires, la modélisation numérique a été mise en œuvre en utilisant le logiciel FEM "Comsol Multiphysics". Une série des modèles simulant les champs de températures, les mouvements convectifs et le malaxage (diffusion, méthode level set, méthode des champs de phases) a été créée. Dans le cas du laser, la formulation pseudo-stationnaire du transfert de chaleur basée sur la géométrie du capillaire simplifiée et la convection a été couplée avec les problèmes 2D de diffusion et de malaxage des matériaux dans différents plans horizontaux. En soudage par faisceau d'électrons, la morphologie de la microstructure a nécessité une formulation temporelle. Le modèle multiphysique final en couplage complet (solution multiphysique simultanée) reproduit le processus de formation d'une structure périodique de solidification lors du soudage par faisceau d'électrons et permet d'expliquer l'aspect des structures alternées entre matériaux immiscibles ou présentant de grandes différences de propriétés thermophysiques.Le deuxième couple de matériaux présente des problèmes métallurgiques majeurs liés à la formation des phases intermétalliques rendant l'assemblage direct par fusion impossible. La composition locale devient donc l'aspect-clef de la formation d’une soudure correcte : l'introduction d’un troisième matériau (cuivre) ayant une meilleure compatibilité avec le titane est nécessaire. Pour pouvoir déterminer les fenêtres optimales des conditions opératoires, les modèles numériques, créés précédemment, ont été adaptés pour quatre procédés de l’assemblage : faisceau d'électrons, soudage lasers Nd:YAG continu et pulsé, brasage par laser avec apport de fil. L'analyse élémentaire des microstructures dans les soudures résistantes mécaniquement a permis de développer le scénario de la solidification d'une zone fondue et de comprendre l'influence de la composition aux interfaces sur la résistance mécanique des assemblages.Les modèles numériques multiphysiques créés au cours de cette étude permettent l'accès rapide à la grande quantité d'information sur le comportement de la zone fondue en fonction des paramètres de soudage en se basant sur le nombre des données de départ relativement limité et sur quelques hypothèses simplificatrices. L'approche multiphysique à la modélisation de soudage permet de reproduire la forme de la zone fondue, visualiser les écoulements du liquide et cartographier la distribution de certains éléments avec une bonne corrélation avec les résultats expérimentaux. L'ensemble des modèles permet de déterminer les conditions opératoires répondant aux critères fixes en fonction de la métallurgie d'un couple hétérogène. / The present study is dedicated to the comprehension of the mechanism of materials mixing during dissimilar welding by high power beam sources. We have been interested in joining of two couples of metallic materials which present different metallurgical problems: • copper- stainless steel (miscibility gap, important difference in physical properties);• TA6V- stainless steel (oxidation on air, formation of intermetallic phases which made the joint brittle).For the first couple of materials, continuous laser Nd:YAG welding and electron beam welding have been applied. The experimental study of morphology evolution, composition, microstructure and mechanical properties has allowed establishing the hypotheses on formation of heterogeneous mixture between the materials having limited solubility. To quantify the physical phenomena of continuous dissimilar welding, the numerical modeling has been carried out by means of FEM software package "Comsol Multiphysics". A number of models reproducing temperature field, convection movements and mixing (diffusion, level set method, phase field method) between the materials has been created. In case of continuous laser welding, the pseudo-stationary formulation of heat transfer based on simplified key-hole geometry and convection has been coupled with two-dimensional problems of diffusion and mixing in horizontal planes. The electron beam welding presenting the nonlinear development of the weld has needed employing of temporary formulation. Final model including complete coupling (simultaneous multiphysical solving) reproduces the process of development of periodic solidification structure during electron beam welding and allows explaining the mechanism of formation of altered structures between immiscible materials which have important difference in thermophysical properties.The second couple of materials presents weldability problems due to formation of brittle intermetallic phases making direct joining by fusion impossible. The local elementary composition becomes the key-aspect of successful joining: the introduction of the third material (pure copper) having better compatibility with titanium is necessary. To determine the ranges of optimal operational conditions, numerical models created previously have been adapted to the case of four joining techniques: electron beam and laser Nd:YAG (continuous and pulsed) welding and laser brazing with filler wire. Elementary analysis of microstructures of resistant welds has allowed developing the solidification scenario and understanding the influence of local composition of heterogeneous interfaces on tensile properties of the joints. The multiphysical models created during this study allow rapid access to high quantity of data on behavior of melted zone in function of welding parameters basing on relatively limited input data and several simplification hypotheses. The multiphysical approach to welding modeling allows recreating the shape of melted zone, to visualization the convection movements and providing the cartography of several elements in good correspondence with experimental results. A set of models allows determination of operational parameters respecting fixed criterions determined by metallurgy of dissimilar couple.