La motivation principale de cette thèse était de pouvoir simuler numériquement les procédés de mise en œuvre des propergols, comme le mélange et le remplissage de Booster. La viscosité de ce type de fluide implique durant la phase de remplissage l'apparition d'oscillations toroïdales caractéristiques du phénomène appelé Fluid Buckling. Ce phénomène est particulièrement difficile à représenter numériquement car l'interface entre le fluide et l'air présente une surface libre dont la complexité augmente avec le temps. C'est pourquoi, c'est cet exemple qui a été choisi pour valider l'efficacité d'une méthode numérique de représentation des surfaces libres. Nous avons pour cela choisit d'implémenter dans la librairie de calcul CimLib une méthode de type Level Set. Celle-ci a pour particularité d'effectuer l'étape de convection et l'étape de réinitialisation, indispensables à ce type de résolution, simultanément. De plus, la fonction distance a été tronquée aux alentours de l'interface à l'aide d'une fonction sinus, obtenant ainsi une fonction niveau sinusoïdale auto-déterminante. La résolution est stabilisée à l'aide d'une méthode de type SUPG. L'implémentation finale de la méthode est facilitée puisque le problème est réduit à une unique équation de convection d'une fonction Level Set Locale. Après une série de tests classiques de validation de méthode de surface libre, notre méthode a donc été validée en 2 et 3 Dimensions en simulant l'apparition du Fluid Buckling. La méthode a ensuite été implémentée dans les deux logiciels de simulation et validée sur des cas tests industriels prédéfinis de remplissage et de mélange. Les défauts obtenus numériquement correspondent à ceux attendus expérimentalement. La méthode est aujourd'hui utilisée pour simuler de nombreux autres procédés et a montré sa robustesse particulièrement dans des procédés multifluides. Afin d'optimiser la visualisation des défauts d'interfaces et de diminuer les temps de calcul, on peut envisager comme suite de ces travaux de recherche l'adaptation du maillage à proximité de l'interface et du pas de temps, le tout de façon automatique.

Ce travail concerne la simulation d'écoulements viscoélastiques compressibles appliquée à l'injection de polymères. La compressibilité est intégrée dans Rem3D en supposant que la densité du matériau suit une loi d'évolution du type loi de Tait. La conservation de la masse est écrite comme une équation en vitesse, pression et température, à travers des coefficients de compressibilité isotherme et de dilatation isobare. Le système obtenu est désigné "Stokes compressible" et sa résolution numérique est faite par la méthode des éléments finis mixtes. Le système obtenu est non-linéaire et non-symétrique. Le couplage thermique et l'extension à des problèmes avec surface libre sont aussi consid ér és. Le modèle viscoélastique choisi est le modèle Pom-Pom, issu de la dynamique moléculaire. L'extra-contrainte est fonction des propriétés microscopiques du matériau, comme l'orientation moléculaire et son étirement. L'élasticité est vue comme une perturbation dans le problème mécanique, et une méthode de stabilisation du type DEVSS est utilisée. L'orientation et l'étirement sont déterminés par la résolution de deux équations d'évolution via une méthode espace-temps Galerkin discontinu. Finalement, la thermoviscoélasticité est abordée brièvement. Dans le contexte de l'injection de polymères, REM3D couvre aujourd'hui toutes les phases du proc éd é. Néanmoins, la solidification et la transition liquide-solide sont approximées par un comportement du type liquide de très haute viscosité. L'introduction de la compressibilité permet de compenser le retrait du matériau par un apport supplémentaire de matière. D'un autre côté, la prise en compte d'un comportement viscoélastique détecte d'éventuelles anisotropies des propriétés de la pièce injectée. Les diverses comparaisons des résultats obtenus avec la littérature et l'expérience montre une bonne concordance, validant les modèles implémentés.