La détection des défauts non francs est un passage obligé dans la gestion de la maintenance des câbles (Wire Health Monitoring) et permet d’anticiper l’apparition de défauts francs engendrés par l'aggravation de défauts non francs. Ces travaux de thèse visent à étudier la problématique de défauts non francs dans les câbles électriques. L'étude proposée consiste à étudier le problème direct : modélisation électromagnétique et compréhension des phénomènes physiques liés à la présence des défauts non francs, et impact des défauts sur leurs signatures obtenues par réflectométrie. Il est proposé dans ce travail de caractériser des défauts non francs représentatifs de situations réelles, ce qui peut servir ultérieurement dans la résolution du problème inverse : déterminer la localisation et la sévérité des défauts à partir du réflectogramme en vue du diagnostic. L'approche proposée se fait en deux temps. Dans un premier temps, une caractérisation électromagnétique d’une zone localisée d’un câble présentant un défaut non franc est réalisée grâce à une modélisation électromagnétique prenant en compte la géométrie tridimensionnelle du défaut. Deux méthodes numériques ont été exploitées : une méthode de différences finies dans le domaine temporel (FDTD) et une méthode de Galerkin discontinu. Les résultats des simulations permettent d'exprimer les perturbations engendrées par le défaut et de déduire leur influence sur la propagation des signaux dans le câble. Des validations expérimentales menées dans le domaine temporel et le domaine fréquentiel permettent de confronter les mesures aux prédictions numériques. Dans un deuxième temps, des modèles électriques de défauts sont exprimés sous forme de paramètres localisés R, L, C, et G qui peuvent être utilisés dans un modèle de ligne de transmission. Une telle analyse des défauts non francs permet de relier une faible variation locale des caractéristiques physiques et électriques de la ligne à une variation des signaux de réflectométrie et des paramètres électriques de la ligne. L’approche permet de fournir des informations utiles pour caractériser des défauts et peut ainsi contribuer à améliorer les performances des systèmes de détection / The soft fault detection feature is certainly a very important aspect of wire health monitoring and an important process required in electrical wiring system operation. It has a great influence on the security and quality of supply. In transmission line networks, this feature is needed to provide a timely identification of the faulted line thus anticipating the appearance of severe faults that are initially caused by soft fault degradation. This work focuses on soft fault problems in electrical fault diagnosis and their weak impact on coaxial transmission lines. The objective of this work is to carry out a soft fault forward model: electromagnetic modeling and investigating the behavior of the line after soft damages and then to analyze its effects on their signatures obtained by reflectometry. It is proposed in this work to characterize the representative soft damages in real situations, which could be used later in solving the inverse problem: determining the position and severity of defects from the reflectometry response for the diagnosis. The proposed approach is based on the following steps: first, an electromagnetic characterization of a faulty region of a cable is carried out by electromagnetic modeling, by taking into account the three-dimensional geometry of defect. For this purpose, two numerical methods have been used: Finite Difference Time Domain (FDTD) and a Discontinuous Galerkin. The simulation results allow to study the disruption initiated by the fault and to infer their influence on the signal propagation along the cable. The experimental validation provided in frequency- and time-domain allows to confront experimental measurements with simulation predictions. In a second step, electrical fault models are expressed in terms of lumped parameters R, L, C, and G, which can be used in a transmission line model. Such analysis of soft faults allows to relate low local variation of the physical and electrical characteristics of the line to a reflectometry signals variation and electrical parameters changes. The approach can provide useful information to characterize defects and can thus contribute to improve the performance of detection systems

Le domaine de la simulation numérique évolue vers des simulations de phénomènes physiques toujours plus complexes. Cela se traduit typiquement par le couplage de plusieurs codes de simulation, où chaque code va gérer une physique (simulations multi-physiques) ou une échelle particulière (simulations multi-échelles). Dans ce cadre, l'analyse des résultats des simulations est un point clé, que ce soit en phase de développement pour valider les codes ou détecter des erreurs, ou en phase de production pour confronter les résultats à la réalité expérimentale. Dans tous les cas, le pilotage de simulations peut aider durant ce processus d'analyse des résultats. L'objectif de cette thèse est de concevoir et de réaliser une plate-forme logicielle permettant de piloter de telles simulations. Plus précisément, il s'agit à partir d'un client de pilotage distant d'accéder ou de modifier les données de la simulation de manière cohérente, afin par exemple de visualiser "en-ligne" les résultats intermédiaires. Pour ce faire, nous avons proposé un modèle de pilotage permettant de représenter des simulations couplées et d'interagir avec elles efficacement et de manière cohérente. Ces travaux ont été validés sur une simulation multi-échelles en physique des matériaux. / The numerical simulations evolve more and more to simulations of complex physical phenomena through multi-scale or multi-physics codes. For these kind of simulations data analysis is a main issue for many reasons, as detecting bugs during the development phase or to understand the dynamic of the physical phenomena simulated during the production phase. The computational steering is a technique well suited to do all this kind of data analysis. The goal of this thesis is to design and develop a computational steering framework that take into account the complexity of coupled simulations. So, through a computational steering client we want to interact coherently with data generated in coupled simulations. This afford for example to visualize on-line the intermediate results of simulations. In order to make this possible we will introduce an abstract model that enables to represent coupled simulations and to know when we can interact coherently with them. These works have been validated on a legacy multi-scale simulation of material physics.

Ce mémoire s'intéresse aux propriétés de transfert, de résistance et de plasticité de roches argileuses. Deux enjeux industriels directs sont la détermination du potentiel de production de gaz naturel des roches mères et de la capacité de rétention des roches de couvertures utilisées pour les stockages de gaz ou de déchets en milieu souterrain. Dans un premier temps, la perméabilité est décrite comme résultant de l'homogénéisation de l'écoulement d'un fluide visqueux dans le réseau poreux. Un cadre variationnel analogue à celui de Hashin et Shtrikman en élasticité est établi et mis à profit pour proposer une méthode numérique basée sur la Transformée de Fourier Rapide (FFT) pour le calcul efficace de bornes sur la perméabilité à partir de représentations tridimensionnelles voxelisées de milieux poreux. Par ailleurs, des techniques d'homogénéisation analytiques basées sur la définition de cellules perméables équivalentes sont proposées pour fournir des briques élémentaires pour la modélisation micro-mécanique de la perméabilité.Le phénomène de glissement du gaz aux parois est considéré pour rendre compte quantitativement d'effets du type Klinkenberg sur la perméabilité au gaz. L'effet d'une saturation partielle du réseau poreux sur la perméabilité et la pression de percée sont étudiés dans le cas de morphologies particulières des pores. En parallèle, une campagne expérimentale est conduite sur différentes roches argileuses issues de forages pour caractériser leur porosité et leur perméabilité sous diverses conditions de confinement mécanique et desaturation partielle en eau. Dans un second temps, les capacités de résistance sont étudiées sous hypothèse de rupture ductile. Une méthode numérique efficace basée sur la FFT est proposée pour encadrer par l'extérieur le domaine de résistance de milieux hétérogènes à géométrie complexe. D'autre part, une modélisation analytique à trois échelles de la résistance d'un matériau granulaire renforcé par inclusions rigides avec interfaces imparfaites est présentée. Enfin, une relecture élasto-plastique de ce modèle de résistance a abouti à un modèle purement micro-mécanique, dont l'interprétation macroscopique est analogue au modèle Cam-clay, avec écrouissage par changement de porosité et ligne d'états critiques / This thesis is focused on the strength, plasticity and transport properties of mudrocks. Two industrial applications are the shale gas production and the underground gas or waste storage. In a first part, the permeability is described as resulting from the homogenisation of the flow of a Newtonian fluid within the pore space. A Hashin-Shtrikman like framework is derived for the permeability upscaling and used to propose a FFT-based numerical method for the efficient computation of bounds on the permeability, directly compatible with a voxelised representation of the pore space. As an alternative, analytic homogenisation techniques based on the definition of equivalent permeable cells are developed to provide building blocks for the micro-mechanical modelling of permeability. The gas slip at pore walls is accounted for to model Klinkenberg effects for gas permeability. Partial water saturation is also considered to model relative gas permeability and gas breakthrough pressure. In the mean time, a thorough experimental investigation of the evolution of porosity and permeability with confining pressure and partial water saturation has been carried out on several types of mudrocks. In a second part, the ductile strength properties is studied. An efficient FFT-based numerical method is proposed to compute the homogenised strength domain of heterogeneous media with complex micro-structures. Next, a three-scale analytic model of the strength of a granular media reinforced by rigid inclusions with imperfect interfaces is presented. In a third part, this strength model is re-interpreted in plasticity to propose a purely micro-mechanical model, whose macroscopic interpretation is similar to the Cam-clay model, including hardening or softening due to an evolving porosity and a critical state line

L'objectif de ce travail de thèse est de dégager des méthodes et des outils permettant de mieux comprendre le rôle joué par les phénomènes de transport (cellulaire, hydraulique et chimique) dans la mise au point de stratégies thérapeutiques de réparation osseuse. Pour cela, nous avons choisi d'associer deux approches : la réalisation d'études expérimentales et la mise au point de modèles numériques. Nous avons ainsi pu, lors d'une première étude présentée dans le chapitre 2 de ce document, relier la perméabilité intrinsèque d'un milieu poreux, paramètre déterminant dans l'étude du transport de fluide en son sein, à la structure géométrique de ses pores. Nous avons également mis en évidence l'importance des interactions électrochimiques lors de la progression d'une solution ionique (telle que les fluides physiologiques) à travers le tissu osseux, en raison de la structure poreuse et de la composition chimique (présence de fibres de collagènes chargées par exemple) de ce dernier. Ces outils ont ensuite permis d'analyser, en première approche, les résultats expérimentaux obtenus lors de la réalisation de tests de perméabilité sur des échantillons de périoste fémoral ovin, dans le but d'identifier les phénomènes physico-chimiques à l'origine du comportement particulier de cette membrane (chapitre 5). Nous nous sommes par ailleurs intéressés au développement d'implants osseux associant un substrat minéral biocompatible et des cellules souches mésenchymateuses, afin de favoriser une reconstruction tissulaire en volume des lésions de grande taille. Nous avons ainsi pu mettre en place, dans le chapitre 3, un dispositif expérimental permettant de réaliser de manière reproductible un test d'ensemencement cellulaire et d'évaluer le nombre, la répartition et le taux de viabilité des cellules greffées sur le biomatériau utilisé. A partir des résultats expérimentaux issus des tests d'ensemencement cellulaire, nous avons ensuite développé un modèle numérique dans le chapitre 4, pour dégager un ensemble de critères à respecter dans l'élaboration d'un substitut osseux qui favoriserait un développement tissulaire homogène contrôlé lors des premières étapes de la culture in vitro de ce type d'implants. Ce modèle constitue une première étape dans la détermination d'un cahier des charges géométrique de tels substrats / This study aims to set up methods and tools to improve our understanding of the role played by transport phenomena (transport of cells, fluid and chemical species) in the development of new therapeutic protocols for bone reconstruction, using a double approach: experimental studies and numerical simulations. Hence, in the second chapter of this document, we have been able to link the intrinsic permeability of a porous medium – a key parameter regarding fluid transport through porous media – to the geometric structure of its pores. We have also highlighted the influence of electrochemical interactions on the flow of an ionic solution (such as physiologic fluids) through cortical bone, due to its porous structure and its chemical composition (presence of electrically charged fibers). These tools have then enabled us to analyze, at first glance, the experimental results of permeability tests conducted on ovin femoral periosteum, to identify the chemical-physical phenomena responsible for the specific behavior of this membrane (chapter 5). We also focused on the development of large bone implants coupling a mineral substitute and mesenchymal stem cells to enhance a volumic reconstruction of critical-sized bone defects. We have therefore designed, in chapter 3, a custom experimental set up that allows one to perform a reproducible cell seeding test on a porous scaffold and quantify the number of seeded cells as well as their viability rate. The experimental results provided by these tests have then initiated the numerical model exposed in chapter 4, that aims to highlight criteria to meet regarding the design of new bone substitutes that would enhance a homogeneous volumic tissue growth during the first stages of the extit [in vitro} development of coupled implants

Le rôle des tourbillons anticycloniques dans l'évolution des disques protoplanétaires et, en particulier, dans les mécanismes de formation des planétésimaux, est au coeur des défis actuels de l'astrophysique moderne. C'est pourquoi une étude approfondie de leur structure et de leur dynamique est primordiale.Grâce à un outil numérique spécifiquement développé pour l'étude des disques, nous avons revisité l'Instabilité en Ondes de Rossby dans le régime non linéaire, et découvert l'existence d'une cascade des modes de perturbation qui permet de mieux comprendre la formation des tourbillons par cette instabilité.Leur structure à été décrite par un modèle gaussien innovant, remarquablement en accord avec les résultats numériques. Grâce à un échantillon de près de 300 tourbillons, nous avons borné le domaine des dimensions radiales, azimutales et de la vorticité. Deux familles de tourbillons possibles ont été distinguées : (i) les tourbillons incompressibles, stables et quasi-stationnaires; (ii) les tourbillons compressibles, très mobiles et associés à l'émission d'ondes de densité. Leur persistance sur plus de 1000 rotations confirme l'observabilité de tous ces tourbillons. Enfin, nous avons caractérisé leur migration vers l'étoile en fonction de leur géométrie, du gradient de pression et de l'échelle de hauteur du disque. Pour la première fois, une expression analytique permet d'estimer le taux de migration en fonction de ces paramètres; l'échelle de temps pour tomber sur étoile peut aller de 10^6 à 100 rotations. Suivant un modèle de viscosité alpha, la perte de moment cinétique pourrait être suffisante pour maintenir un taux d'accrétion significatif dans la zone morte. / The role of anticyclonic vortices in the protoplanetary disk evolution and in how do planetesimals form are among the most important chalenges of the modern astrophysics. That is why an exhaustive study of the structure and the evolution of these vortices is necessary.Thanks to a numerical code specificly designed for the study of these disks, we have revisited the Rossby Wave Instability in the nonlinear regime, and have discovered that a cascade of the perturbation modes can explain the formation of the vortices created by this instability.We have described the structure of these Rossby vortices with a new gaussian vortex model, which accurately fits the numerical results. A sample of 300 different vortices led us to define the bondaries of the radial and azimuthal extent as well as the vorticity of the vortices. We have distinguished two main families : (i) the incompressible family, which is stable and quasi stationnary ; (ii) the compressible family, moving and exciting density waves. We found them surviving more than 1000 orbits, a clear confirmation of their observability.Finaly, we have caracterized the inward migration of the vortices as a fonction of their shape, their vorticity, but also of the pressure gradient and the scale height of the disk. For the first time, we exhibit a equation relating the migration rate to these parameters. The time scale of the migration ranges from 10^6 to just 100 rotations of the disk. Extremely steep pressure gradients are needed to reverse the migration to an outward regime. Following the alpha viscosity approch, the loss of kinetic momentum due to this migration would be sufficient to sustain the accretion in the dead zone.

La dégradation géométrique de la voie ballastée sous circulation commerciale nécessite des opérations de maintenance fréquentes et onéreuses. La caractérisation du comportement des pro-cédés de maintenance comme le bourrage, la stabilisation dynamique, est nécessaire pour proposer des améliorations en terme de méthode, paramétrage pour augmenter la pérennité des travaux. La simulation numérique d'une portion de voie soumise à un bourrage ou une stabilisation dynamique permet de comprendre les phénomènes physiques mis en jeu dans le ballast. Toutefois, la complexité numérique de ce problème concernant l'étude de systèmes à très grand nombre de grains et en temps de sollicitation long, demande donc une attention particulière pour une résolution à moindre coût. L'objectif de cette thèse est de développer un outil de calcul numérique performant qui permet de réaliser des calculs dédiés à ce grand problème granulaire moins consommateur en temps. La méthodologie utilisée ici se base sur l'approche Non Smooth Contact Dynamic s(NSCD) avec une discrétisation par Éléments Discrets (DEM). Dans ce cadre, une méthode de dé-composition de domaine (DDM) alliée à une parallélisation adaptée en environnement à mémoire partagée utilisant OpenMP sont appliquées pour améliorer l'efficacité de la simulation numérique. / The track deterioration rate is strongly influenced by the ballast behaviour under commercial traffic. In order to restore the initial track geometry, different maintenance processes are performed, like tamping, dynamic stabilisation. A better understanding of the ballast behaviour under these operations on a portion of railway track is a key to optimize the process, to limit degradationand to propose some concept for a better homogeneous compaction. The numerical simulation isdeveloped here to investigate the mechanical behaviour of ballast. However, the main difficultiesof this research action concerns the size of the granular system simulation increasing both in termof number of grains and of process duration. The purpose of this thesis is to develop an efficient numerical tool allows to realize faster computations devoted to large-scale granular samples. In this framework, the Non-Smooth Contact Dynamics (NSCD) of three-dimensional Discrete ElementMethod (DEM) simulations, improved by Domain Decomposition Method (DDM) and processedwith the Shared Memory parallel technique (using OpenMP) has been applied to study the ballast media mechanics.

L'utilisation croissante de véhicules au cours des dernières décennies a causé une augmentation dans les émissions de gaz d'échappement nocifs provoquant des problèmes de santé et d'environnement. Cette problématique a conduit les gouvernements à mettre en place des limites d'émissions plus strictes. Afin de respecter ces limites, il est nécessaire de développer des systèmes de post-traitement plus performants. Parmi plusieurs solutions possibles, l'analyse et l'optimisation de l'impact de l'écoulement sur les réactions chimiques est une approche importante pour résoudre ce problème complexe. Au point de vue expérimental, il est très difficile de mesurer le champ de vitesse détaillé et la distribution de la température et de la concentration. Les simulations numériques peuvent fournir des informations supplémentaires pour comprendre l'interaction entre la distribution de l'écoulement et l'efficacité des réactions. La simulation numérique de ces systèmes est très coûteuse due aux phénomènes physiques et chimiques complexes ayant lieu simultanément dans tout le système. Afin de diminuer ce coût, on peut développer des modèles physiques et chimiques simplifiés ainsi que des techniques numériques spécifiques pour simuler le système tout en réduisant le temps de calcul. Dans cette thèse, nous développons une approche générale pour modéliser et simuler le système de post-traitement des gaz d'échappement automobile contenant les aspects physiques et chimiques. L'approche présente combine des modèles simplifiés et des techniques numériques de multi-résolutions afin de capturer correctement les caractéristiques de l'écoulement dans le système. Alors que la méthode de raffinement de maillage adaptatif (AMR) est optimisée afin de minimiser l'effort de calcul dans les régions divergente et convergente, un modèle de sous-maille est développé pour décrire l'écoulement dans les micro-canaux du substrat catalytique placé entre la région divergente et la région convergente. La performance du modèle est validée par rapport aux résultats expérimentaux obtenus par Benjamin et al. (2002). Le couplage effectif de méthode AMR et le modèle de sous-maille permet de capturer les caractéristiques de l'écoulement dans le système avec un temps de calcul réduit de manière significative. L'impact de l'écoulement pulsé et de l'écoulement en expansion sur le taux de réaction dans un canal de monolithe est étudié numériquement. La méthode AMR est montrée pour capturer les couches mince de diffusion près de la paroi à l'interface solide-liquide. Sur la base des résultats numériques, nous proposons un modèle simplifié de transport capturant les effets de l'écoulement, la diffusion et la réaction catalytiques à la paroi. Le modèle simplifié de transport peut être directement appliqué dans la forme de modèle de sous-maille pour la description complète de tous les processus physiques et chimiques ayant lieu à l'intérieur du système. Grâce aux approches de modélisation numériques physiques et chimiques développées, il est possible de simuler le système en trois dimensions avec un temps de calcul raisonnable capturant encore la physique principale du problème. / The increasing use of passenger vehicles over the past decades has caused an increase in harmful exhaust gas emissions which give rise to environmental and health problems. This problematic has led governments to establish very stringent emission limits. The emission restrictions require more performing after-treatment systems. Among many other solutions, analysing and optimising the flow impact on the conversion efficiency is an important step towards the solution of the complex engineering problem. Detailed velocity, temperature and concentration distributions are very difficult to measure experimentally. Numerical simulations can provide additional information to understand the interaction of flow distribution and conversion efficiency. The simulation of these systems is computationally very expensive due to complex physical and chemical phenomena occurring simultaneously throughout the system. To overcome this cost, one can resort to some simplified physical and chemical models together with specific numerical techniques to simulate the system with a reduced computational time. In this dissertation, we develop a general approach to model and simulate the automotive catalytic converter system including all the physical and chemical processes. The present approach combines sub-grid models and numerical multi-resolution techniques in order to correctly capture the flow characteristics inside an automotive catalytic converter. While Adaptive Mesh Refinement (AMR) techniques are optimized in order to minimise the computational effort in the divergent and convergent regions, a sub-grid model is developed to describe the flow inside the catalytic substrate placed between the convergent and divergent regions. The performance of the sub-grid model is validated against the experimental results obtained by Benjamin et al. (2002). The effective coupling of AMR techniques and the sub-grid model allows to capture the flow features with significantly reduced computational time. The impact of pulsating and expansion flow on the conversion efficiency within a single monolith channel is investigated numerically. AMR techniques are shown to capture the small boundary layers near the wall at the solid-fluid interface. Based on the numerical results, we propose a simplified transport model that captures the effects of flow, diffusion and catalytic wall reactions. The simplified transport model can be directly applied as a sub-grid model for the complete description of all the physical and chemical processes taking place inside the system. The developed physical, chemical and numerical modelling approaches make the three-dimensional simulations possible with a reasonable computational time still capturing the main physics of the problem.

Les travaux présentés dans cette thèse concernent la simulation numérique du démarrage et de la saturation d'un moteur thermoacoustique à ondes stationnaires couplé à une charge résistive placée à une de ses extrémités. Le modèle utilisé est hybride : l'écoulement dans la cellule est décrit dans le cadre bidimensionnel et non linéaire ; il est couplé à un modèle monodimensionnel d'acoustique linéaire du résonateur, par raccordement dans la limite asymptotique de faible nombre de Mach. Le système sélectionne naturellement les modes acoustiques instables. L'analyse du signal temporel de pression issu de la simulation numérique permet de calculer les taux de croissance des modes dominants. On peut ainsi déterminer la température critique de l'échangeur chaud permettant au moteur de démarrer ainsi que la fréquence du mode associé. Les variations de la température critique et de la fréquence d'oscillation sont caractérisées pour toutes les valeurs possibles de la charge résistive. L'influence de paramètres physiques tels que la pression moyenne, ou de divers paramètres géométriques de la cellule active est également examinée. Les résultats sont confrontés avec succès à la théorie linéaire et avec des résultats expérimentaux de la littérature. Le modèle hybride a permis, dans certains cas, de mener les simulations jusqu'à l'obtention du régime périodique, ce qui représente plusieurs dizaines de milliers de périodes acoustiques. Enfin, deux simulations en régime périodique sont détaillées pour analyser la dynamique de l'écoulement (formation de tourbillons) au voisinage des extrémités des plaques du stack et des échangeurs, pour des ondes de faible ou de forte intensité. / The work presented in this thesis concerns the numerical simulation of the starting phase and saturation regime of a thermoacoustic engine equipped with a resistive load. A hybrid model is used: the flow in the active cell, described by two-dimensional nonlinear equations, is coupled to a one-dimensional linear acoustics model in the resonator, using matched asymptotic expansions in the low Mach number limit. Unstable acoustic modes develop spontaneously in the system. The computed acoustic pressure signal in the active cell is analyzed in order to extract the growth rate and frequency of the dominant modes. Therefore the critical hot exchanger temperature and frequency of the associated mode allowing the engine to start can be determined. Those critical parameters are characterized for all possible values of the resistive load. The effects of physical parameters such as mean pressure or of geometrical parameters of the active cell are also investigated. Results are found in agreement with linear theory and with experimental results from the literature. In some instances, the hybrid model enables to carry the simulations up to the periodic regime, which represents tens of thousands of acoustic periods. Finally, two simulations of the periodic regime are detailed in order to analyze flow dynamics (vortex formation) in the vicinity of the stack plate/heat exchanger extremities, for small and large drive ratio.

Depuis l'invention du soudage par frottement malaxage (FSW) en 1991 par The Welding Institute (TWI), l'industrie aéronautique envisage de remplacer l'assemblage par rivetage des tôles par ce procédé. L'objectif est d'alléger les structures, plus particulièrement les panneaux des structures d'aéronefs qui sont souvent raidis avec des lisses et des profilés assemblés à la peau extérieure. La qualification du procédé FSW dans l'aéronautique requiert de produire des joints de très haute qualité avec une longue durée de vie des outils. Ceci nécessite un réglage fin des paramètres procédé tels que la géométrie de l'outil (dimensions de l'épaulement et du pion, filetage) et les vitesses d'avance et de rotation. Le but de ce travail est d'accompagner la conception du procédé par la simulation numérique. Il comporte trois parties: (1) le développement de la formulation arbitrairement lagrangienne ou eulérienne (ALE) dans la version parallèle du logiciel Forge® pour l'adapter au FSW, (2) la modélisation du procédé FSW et (3) la calibration du modèle EF et la validation des résultats numériques suite à la réalisation d'essais FSW en configuration par recouvrement.Ce travail est basé sur la formulation ALE parallèle développée dans Forge® afin de modéliser les différents défauts caractéristiques du FSW (bavures et trous). L'algorithme de transport des champs continus, indispensable au caractère eulérien, est affiné afin de considérer les particularités de l'écoulement du procédé (nœud du plan d'entrée, en contact avec l'outil ou dans les bavures et sortant du domaine d'étude). Un nouveau schéma d'intégration temporelle basée sur les coordonnées cylindriques plus adapté à ce procédé fortement rotationnel est introduit. Enfin, un filetage sur l'outil est modélisé par l'introduction d'une nouvelle loi de frottement à l'interface outil / matière.Cette formulation ALE générale et parallèle a montré sa robustesse pour modéliser le procédé FSW. Les résultats thermomécaniques sont en accord avec de précédents résultats validés par l'expérience. De plus l'outil numérique obtenu a montré ses capacités pour la modélisation des bavures ou des pertes de contact derrière le pion.Finalement la configuration de soudage par recouvrement industrielle a été étudiée. Pour cela des essais expérimentaux ont été réalisés pour divers paramètres du procédé et divers types de joints de recouvrement. Puis les mesures de forces et du couple dans l'outil et de la température dans l'outil et la tôle ont permis de calibrer les paramètres du modèle (frottement, loi de comportement, coefficients thermiques) et de valider les résultats EF. / Since the invention of Friction Stir Welding (FSW) in 1991 by The Welding Institute (TWI), aeronautics industry has been investigating the possibility to use this process instead of riveting with the objective to lighten its structures and more particularly the aircrafts structure panels. Aircraft panels are often straightened with stringers and profiles which are joined to the outer skin. The qualification of FSW process in aeronautics requires producing very high quality joints with the longevity of tools, which requires fine tuning of process parameters such as the geometry of the welding tool (shoulder and pin dimensions, threads on pin and shoulder) and the advancing and rotating speeds. The aim of this work is to support the design of the process by numerical simulation. It has three parts: (1) developing an efficient and accurate Arbitrary Lagrangian or Eulerian (ALE) formulation within the parallel version of Forge® software, (2) modelling the FSW process and (3) calibrating the F.E. model and validating simulation results thanks to FSW experiments on lap joints.This work is based on the parallel ALE formulation developed in Forge® to model the different possible defects taking place in FSW (flashes and worm holes). The transport algorithm of continuous fields, required by the Eulerian frame, is enhanced to take the special characteristics of the FSW's flow into account (nodes located in input plan or flashes or in contact with the tool). A new time integration scheme based on cylindrical coordinates, which are best suited for this process, is introduced. Finally, the pin and shoulder threads are modelled by introducing a new friction law at the tool / material interface.This general and parallel ALE formulation is robust enough to model the FSW process. Thermomechanical results obtained are in agreement with previous results validated by experiences. And the numerical tool demonstrated its ability to model flashes formation and losses of contact behind the pin. Finally industrial welding lap joints configuration was studied. Experimental tests were conducted with several process parameters and type of lap joint. And measure of torque and forces in tool, and temperature in tool and sheet metal allowed to calibrate model parameters (friction, behavior law, thermal coefficients) and to validate FE results.

La première partie de ce travail a porté sur la mesure du gradient de vitesse et la direction de l'écoulement local de liquide dans un empilement de sphères. Cette étude a été menée grâce à des micro-électrodes tri-segmentées implantées à la paroi de sphères supports. Ces expériences originales ont en particulier permis de caractériser les recirculations dans le pore suivant les conditions opératoires. La seconde partie a porté sur la simulation numérique de ce même écoulement local. Un écoulement monophasique dans une géométrie complexe représentant exactement le milieu poreux a été simulé. Ces simulations ont été confrontées aux expériences; l'accord est satisfaisant, surtout à faible vitesse. L'inadéquation des modèles de turbulence usuels a été montrée dans notre cas. Enfin, des mesures par Vélocimétrie Laser Doppler ont montré l'intérêt de cette technique, bien que des améliorations soient nécessaires pour obtenir des résultats exploitables sur l'ensemble du milieu / The first part of this work deals with local measurement of flow direction and velocity gradients in a structured packing of spheres. Tri-segmented micro-electrodes flush mounted at the wall of supporting spheres are used for the measurements. Among other results, these original experiments allowed to charatcerize the flow recirculation within the pore.In the second part, the local liquid flow in the packing was numerically computed. More specifically, the local flow of liquid in a complex geometry was modelled and simulated. Numerical results are in agreement with experiments, especially for low flow rates. Usual turbulence models were shown to be inappropriate for the present case. Finally, some Laser Doppler Velocimetry measurements were carried out and showed the interesting potential of this technique, but improvements are yet to be made to get velocities in the bulk of the bed