Dans le cadre général de la simulation du procédé de frittage en phase solide, ce travail propose une approche numérique, à l'échelle des particules, de la consolidation d'un compact pulvérulent céramique. Le frittage est un procédé mettant en jeu plusieurs chemins de diffusion activés thermiquement. Parmi ces chemins de diffusion, les plus importants sont : la diffusion surfacique, la diffusion aux joints des grains et la diffusion volumique. La mise en place de cette physique dans un cadre de calcul intensif doit permettre de mieux comprendre ces mécanismes de diffusion ainsi que leur influence sur l'évolution de la microstructure. Le but de ce travail consiste à développer un modèle ainsi qu'une stratégie numérique capable d'intégrer les différents mécanismes de diffusion dans un cadre de calcul intensif. Le flux de matière est calculé en fonction du Laplacien de la courbure dans les cas de la diffusion surfacique, tandis que pour la diffusion volumique ce flux est proportionnel au gradient de la pression hydrostatique. Le modèle physique est tout d'abord présenté dans le cadre de la mécanique des milieux continus. Ensuite, la stratégie numérique développée pour la simulation du frittage d'un empilement granulaire est détaillée. Cette stratégie est basée sur une discrétisation du problème par des éléments finis stabilisés couplée avec une méthode Level-set pour décrire la surface libre des particules. Cette stratégie nous permet de faire des simulations avec un "grand" nombre de particules. Plusieurs simulations en 3D, menées dans un cadre de calcul parallèle, montrent l'évolution qui a lieu sur un empilement granulaire réaliste. / Within the general context of solid-state sintering process, this work presents a numerical modeling approach, at the particle scale, of ceramic particle packing consolidation. Typically, the sintering process triggers several mass transport paths that are thermally activated. Among those diffusion paths, the most important ones are: surface diffusion, grain boundary diffusion and volume diffusion. Including this physics into a high-performance computing framework would permit to gain precious insights about the driving mechanisms. The aim of the present work is to develop a model and a numerical strategy able to integrate the different diffusion mechanisms into continuum mechanics framework. In the cases of surface diffusion and volume diffusion, the mass flux is calculated as a function of the surface curvature Laplacian and the hydrostatic pressure gradient, respectively. The physical model describing these two transport mechanisms is first presented within the framework of continuum mechanics. Then the numerical strategy developed for the simulation of the sintering of many particles is detailed. This strategy is based on a discretization of the problem by using a finite element approach coupled with a Level-Set method used to describe the particles free surface. This versatile strategy allows us to perform simulations involving a relatively large number of particles. Furthermore, a mesh adaptation technique allows the particles surface description to be improved, while the number of mesh elements is kept reasonable. Several 3D simulations, performed in a parallel computing framework, show the changes occurring in the structure of 3D granular stacks.

Ce travail de recherche traite du couplage entre un liquide incompressible, irrotationnel et son contenant : une structure élastique. Cette interaction fluide-structure est traitée dans le cadre des petites déformations autour d'un état d'équilibre.Dans un premier temps, on présente une méthode d'introduction des sources dissipatives visqueuses dans le liquide à partir des équations du système couplé conservatif en s'appuyant sur une approche de type fluide potentiel généralement utilisée pour traiter les problèmes de couplage fluide-structure linéarisés non amortis. Un modèle d'amortissement diagonal est alors choisi pour le liquide et les effets dissipatifs de celui-ci sont pris en compte en calculant les coefficients d'amortissement modaux. Seuls les effets dissipatifs liées à la viscosité du liquide sont alors pris en compte. Le système couplé dissipatif obtenu possède une matrice d'amortissement non symétrique. Une résolution de ce système à amortissement non classique est alors présentée et les expressions des réponses fréquentielle et temporelle linéarisées sont données pour différents types d'excitations.Dans un deuxième temps, le liquide est supposé non visqueux et les forces de tension surfacique sont prises en compte. Cette configuration concerne principalement les satellites où le système couplé est en situation de microgravité. Une formulation du problème conservatif permettant de prendre en compte l'incompressibilité du fluide, la condition de continuité à l'interface fluide structure, les effets de capillarité du fluide ainsi que les effets éventuels de précontraintes statiques est alors établie. On se propose pour cela d'utiliser une méthode énergétique via le Principe de Moindre Action. La démarche est alors décomposée en deux étapes : une étude statique afin de déterminer la position de référence, puis une étude dynamique linéarisée autour de cette position d'équilibre. Cette formulation forme notamment une base pour l'introduction des sources dissipatives liées aux effets de capillarité via la méthode précédemment introduite.

Ce travail de recherche traite du couplage entre un liquide incompressible, irrotationnel et son contenant : une structure élastique. Cette interaction fluide-structure est traitée dans le cadre des petites déformations autour d'un état d'équilibre.Dans un premier temps, on présente une méthode d'introduction des sources dissipatives visqueuses dans le liquide à partir des équations du système couplé conservatif en s'appuyant sur une approche de type fluide potentiel généralement utilisée pour traiter les problèmes de couplage fluide-structure linéarisés non amortis. Un modèle d'amortissement diagonal est alors choisi pour le liquide et les effets dissipatifs de celui-ci sont pris en compte en calculant les coefficients d'amortissement modaux. Seuls les effets dissipatifs liées à la viscosité du liquide sont alors pris en compte. Le système couplé dissipatif obtenu possède une matrice d'amortissement non symétrique. Une résolution de ce système à amortissement non classique est alors présentée et les expressions des réponses fréquentielle et temporelle linéarisées sont données pour différents types d'excitations.Dans un deuxième temps, le liquide est supposé non visqueux et les forces de tension surfacique sont prises en compte. Cette configuration concerne principalement les satellites où le système couplé est en situation de microgravité. Une formulation du problème conservatif permettant de prendre en compte l'incompressibilité du fluide, la condition de continuité à l'interface fluide structure, les effets de capillarité du fluide ainsi que les effets éventuels de précontraintes statiques est alors établie. On se propose pour cela d'utiliser une méthode énergétique via le Principe de Moindre Action. La démarche est alors décomposée en deux étapes : une étude statique afin de déterminer la position de référence, puis une étude dynamique linéarisée autour de cette position d'équilibre. Cette formulation forme notamment une base pour l'introduction des sources dissipatives liées aux effets de capillarité via la méthode précédemment introduite. / This study deals with the coupling between an incompressible, irrotational fluid and an elastic container in the context of small amplitude vibrations.Firstly, we present a method to introduce the viscous dissipative sources in the liquid directly from the equations of the conservative coupled problem using a fluid potential approach generally used to treat linear undamped problems. A diagonal damping model is chosen for the liquid and its dissipative effects are taken into account through modal damping coefficients. Only the viscous effects are considered here. The coupled system obtained has a non symmetric damping matrix. This system with non classical damping is solved and expressions of the frequency and linearized time responses are given for different load examples.Secondly, the liquid is supposed to be inviscid and surface tension forces are considered. This configuration is related to satellite applications where the coupled system is in microgravity conditions. A unified formulation of the conservative problem taking into account the fluid incompressibility, the contact condition at the fluid structure interface, capillarity and prestress effects is given. Thus, we propose to use an energy method via the Least Action Principle. The reasoning is then divided into two parts: a static study to determine the reference state and a linearized dynamic study around this equilibrium state. This formulation is a good framework to introduce the dissipative sources associated with the capillary effects by using the method previously introduced.