Le procédé d'hydroformage repose sur la déformation plastique d'une forme première (tôle, tube) sous l'action combinée d'un fluide sous pression et d'un outil rigide. L'objet de ces travaux est de contribuer à la modélisation de ce procédé. Cette contribution s'articule autour de trois problématiques essentielles: la caractérisation rhéologique, la prédiction des défauts de striction et le développement d'algorithmes adaptés à la simulation de l'hydroformage. - L'identification expérimentale des courbes d'écrouissage en grandes déformations est nécessaire afin d'éviter toute extrapolation intempestive au cours d'une simulation numérique. A cet effet, un essai de gonflage circulaire, analysé selon un modèle homogène, a été mis en place. Des résultats très satisfaisants ont pu ainsi être obtenus grâce notamment à la mesure locale du rayon de courbure de la tôle. - La prédiction de la striction au cours de l'opération d'hydroformage est essentielle pour la maîtrise de ce procédé. Dans cette optique, le Critère de Force Maximum Modifié a fait l'objet d'une étude détaillée. Il a été démontré que la forme locale de la surface de charge constitue le facteur le plus influant sur les prédictions du modèle. Le problème de la dépendance de la CLF au trajet a été également abordé, nous avons pu mettre l'accent sur l'utilité d'une représentation des CLF dans l'espace des déformations équivalentes afin de s'affranchir de cette dépendance. Enfin, nous avons proposé un modèle intégrant l'effet de la sensibilité à la vitesse de déformation. La confrontation avec les résultats expérimentaux de tôles en acier a montré une amélioration des prédictions mais le modèle reste relativement conservatif. - Le troisième volet a concerné la modélisation numérique. D'une part, un schéma d'intégration implicite reposant sur une actualisation itérative du maillage et tenant compte des dérivées du domaine a pu être testé. Les résultats ont montré la faible influence des dérivées du domaine. En revanche, l'actualisation itérative semble améliorer la conservation du volume. D'autre part, il est bien établi que certains matériaux peuvent présenter un potentiel de déformation au-delà du pic de pression. Pour une simulation correcte de ce potentiel de déformation, nous avons développé plusieurs variantes de la méthode de continuation. Reposant sur une technique de prédiction-correction, ces méthodes ont permis de dépasser le pic qui reste un point infranchissables pour la méthode conventionnelle de contrôle de pression.

De nombreux matériaux du génie civil, comme le béton ou les enrobés bitumineux, sont des pâtes granulaires, c'est-à-dire des suspensions dont la concentration volumique en inclusions solides est très élevée. Ce travail porte sur la modélisation et la simulation numérique discrète du comportement et du compactage de pâtes granulaires 2D non browniennes et non colloïdales constituées de grains élastiques rugueux de forme circulaire, baignant dans un liquide newtonien de très grande viscosité. Après avoir introduit un modèle de pâte granulaire dans lequel le seul effet du liquide est d'induire des interactions de lubrification entre paires de grains proches, nous présentons une méthode numérique quasi-statique permettant de simuler le comportement de ces systèmes lorsqu'ils sont soumis à un état de contrainte donné. L'étude de l'état d'équilibre adopté par une pâte granulaire lubrifiée lorsqu'on lui applique une pression isotrope montre que, dans la limite où la rugosité des grains est nulle, un tel système se comporte comme un assemblage sec de grains non frottants. C'est la raison pour laquelle nous étudions ensuite le comportement de ces assemblages en 3D lorsqu'on les soumet à des chargements statiques anisotropes ainsi qu'à des essais de cisaillement simple dans la limite quasi-statique. Nous montrons en particulier que les coefficients de frottement interne statique et dynamique sont non nuls et coïncident, que ces assemblages ne présentent aucune dilatance et qu'ils satisfont un critère de rupture de Lade-Duncan. Le comportement d'une pâte granulaire lubrifiée 2D faite de grains de rugosité finie soumise à un cisaillement maintenu est ensuite analysé et nous comparons nos résultats avec ceux obtenus en cisaillant des assemblages granulaires secs 2D. Finalement, nous étudions de manière détaillée les paramètres influençant le compactage d'une pâte granulaire lubrifiée et mettons en évidence la réduction de la mobilisation du frottement intergranulaire accompagnant la densification du matériau.

Nous considérons la simulation numérique des écoulements de fluides viscoélatiques. Développant une approximation en temps par la méthode des directions alternées, nous proposons un algorithme entièrement nouveau permettant de découpler le calcul des contraintes de celui des vitesses. D'ordre deux en temps, cette méthode permet de plus le calcul rapide des solutions stationnaires. L'éléments à divergence nulle de Thomas-Raviart est utilisé pour les vitesses, et celui de Lesaint-Raviart pour les contraintes. La méthode est appliquée au problème de l'écoulement de fluides du type Oldroyd dans une contraction brusque (problème de la marche).

L'amélioration des moteurs à allumage commandé représente un défi de première importance pour les ingénieurs afin de produire plus d'énergie, de consommer moins de matière première et de réduire les émissions polluantes. Les nouvelles technologies apparues ces dernières années amènent les moteurs de plus en plus proches de leurs limites de fonctionnement, favorisant ainsi des phénomènes néfastes qui doivent être contrôlés. Parmi ces phénomènes, les variations cycle-à-cycle (VCC) doivent être minimisées pour garder une performance optimale et éviter une dégradation rapide du moteur. La Simulation aux Grandes Echelles (SGE) est un outil prometteur afin de prédire numériquement les niveaux de variabilités obtenues lors du design d'un moteur (limitant ainsi les coûteuses campagnes de mesures expérimentales). Ce manuscrit s'est attaché à développer une méthodologie numérique pour la prédiction des variabilités cycliques, à simuler un nombre suffisant de cycles pour pouvoir estimer les niveaux de VCC et à valider les résultats obtenus par rapport aux résultats expérimentaux. La SGE semble capter les points de fonctionnements stable et instable étudiés. Les sources qui provoquent ces VCC ont aussi été analysées et une modification du fonctionnement du moteur a été proposée afin de réduire les VCC.

Ce travail porte sur l'étude par simulation en champ de phase du vieillissement isotherme de trois alliages intermétalliques. Nous avons montré que les alliages Ni-Al "inverses", dans lesquels les précipités sont désordonnés et coexistent avec la matrice ordonnée, reproduisent le chemin cinétique des alliages Ni-Al "directs". Une étude statistique de l'alliage Al-Sc a montré que la théorie LSW ne reproduit pas la fonction de distribution de taille de particules sur les stades de coalescence avancés. Nos simulations montrent que la compétition entre les forces motrices chimique, élastique et l'énergie d'interface est à l'origine de la croissance en forme de "papillon" des précipités Al3Sc isolés. Enfin, nous avons adapté le modèle de champ de phase pour l'alliage ternaire Al-Zr-Sc afin de caractériser l'évolution de la composition des précipités Al3(Zr,Sc) durant la coalescence. Nos simulations reproduisent la structure coeur-coquille de ces précipités observée expérimentalement.

Cette thèse porte sur l'utilisation de champs électriques pour faciliter l'élimination de l'eau coproduite avec le pétrole brut, sous la forme d'émulsions stables, lors des étapes d'extraction et de dessalement. Ce procédé, connu sous le nom d'électrocoalescence, s'appuie sur la capacité qu'ont les forces électrostatiques à promouvoir l'attraction et la coalescence de gouttelettes d'eau proches afin d'en augmenter la taille et ainsi d'en accélérer la sédimentation par gravité. Bien que les premières observations expérimentales datent déjà d'un siècle, de nombreuses zones d'ombres subsistent, notamment en ce qui concerne l'optimisation de l'efficacité des électrocoalesceurs de dernière génération. Dans ce contexte, une démarche, combinant simulation numérique multiphysique et expérimentation, a été mise en place pour étudier les phénomènes de mouvement, de déformation et d'instabilité d'interfaces eau-huile induit par la présence d'un champ électrique. La contribution la plus marquante concerne la modélisation et l'analyse de l'effet des forces électrostatiques sur le mécanisme d'amincissement du film d'huile séparant les gouttes. Les résultats numériques mettent en évidence la singularité du problème et l'inadaptabilité des modèles théoriques de lubrification classiquement adoptées pour représenter la coalescence dans les écoulements diphasiques. Une nouvelle expression asymptotique pour le calcul du temps de drainage entre les gouttelettes de l'émulsion est proposée et utilisée pour déduire un critère prédisant la probabilité d'électrocoalescence lors d'une collision dans un écoulement cisaillé. En parallèle, un dispositif sophistiqué, permettant de reproduire expérimentalement le phénomène et d'améliorer la représentativité du critère, a été construit. Enfin, en réponse à un point bloquant décelé lors de la phase de conception de ce dernier, une technique innovante d'injection à la demande de gouttes conductrices non chargées dans un liquide visqueux isolant, utilisant des impulsions électrostatiques, a été développé.

L'objet de cette thèse a été d'étudier la simulation numérique directe d'écoulements diphasiques de fluides non miscibles, incompressibles et isothermes avec la technique du raffinement adaptatif local de maillage. La résolution des équations de Navier-Stokes incompressibles est faite par le biais d'une méthode de projection explicite. La capture de l'interface est réalisée explicitement par la méthode Level Set. L'utilisation de la méthode Ghost Fluid pour le traitement des conditions de saut permet la résolution couplée des deux fluides. Le maillage adaptatif a été implémenté avec les librairies parallèles PARAMESH. Celles-ci gèrent la création d'un maillage construit par blocs qui peuvent être bissectés de façon récursive afin d'obtenir la résolution désirée. Les blocs ont tous le même nombre de cellules, ce qui permet une répartition efficace de la charge de travail en parallle. Un ensemble d'outils nécessaires à une correcte résolution des équations ont été développés, à partir d'un robuste solveur elliptique BiCG-stab préconditionné par méthode multigrille. Le code a été vérifié sur des cas tests académiques capable de maintenir la précision sur les grilles plus fines. Les performances du code en termes de réduction de charge de travail et d'efficacité de la parallelisation ont été également illustrées. Le code a été enfin testé sur la désintégration assistée d'une nappe liquide bidimensionnelle cisaillée par des courants gazeux à haute vitesse. Le code s'est avéré capable de retrouver certains phénomènes physiques comme l'oscillation longitudinale de la nappe, ainsi que de permettre une simulation multi échelles, ce qui permet de plus se rapprocher des conditions des injecteurs réels.

Au cours des procédés de forgeage et d'usinage, le matériau subit des déformations rapides et importantes provoquant des échauffements localisés, des changements microstructuraux et de l'endommagement. La modélisation de tels phénomènes nécessite la connaissance précise des lois de comportement et d'endommagement sur une très large gamme de déformations, vitesses de déformation, et températures. Des essais de compression conduits en utilisant le simulateur thermomécanique GLEEBLE et le banc d'Hopkinson sur deux aciers, 42CrMo4 et 100Cr6, et un alliage d'aluminium 2017-T4, permettent de distinguer trois domaines de comportement : à froid, à mi-chaud et à chaud. Les phénomènes microstructuraux rencontrés dans chaque domaine sont analysés et reliés au comportement rhéologique observé. Après une comparaison de différents modèles de la littérature, un modèle empirique est proposé. Il rend efficacement compte de la compétition entre l'écrouissage et l'adoucissement et permet de refléter, via un couplage, les effets de température et de vitesse de déformation. Une caractérisation des mécanismes d'endommagement est conduite, à froid, à l'aide d'essais de traction in situ dans un Microscope Electronique à Balayage et de la méthode de corrélation d'images numériques. En température, des essais de traction sur éprouvettes axisymétriques entaillées sont utilisés. Deux modèles d'endommagement sont identifiés et analysés, le critère de rupture découplé de Johnson-Cook et le modèle micromécanique couplé de Gurson-Tvergaard-Needleman. Les différents modèles de comportement étudiés sont enfin utilisés et discutés dans le cadre de simulations Eléments Finis de la coupe orthogonale, du cisaillement sur éprouvette " chapeau " et du bipoinçonnement.

La spectrométrie gamma et neutronique forme un puissant outil de caractérisation géologique et géochimique des surfaces planétaires. Cette méthode permet d'aborder des questions capitales en planétologie: composition de la croûte, du manteau; glaces; volcanisme; processus d'altération... La majorité des neutrons et des photons gamma résultent de l'interaction du rayonnement cosmique avec la matière. Le premier chapitre présente la physique de ces interactions nucléaires dans les sols planétaires et dans les détecteurs. Nos travaux visent à optimiser les observations en précisant les performances instrumentales, et en isolant les relations entre la composition du sol et les flux de neutrons. Des simulations numériques utilisant le code GEANT du CERN supportent notre analyse. Le deuxième chapitre évalue les performances du spectromètre gamma à base de germanium pour MARS SURVEYOR 2001. Le résultat des simulations est confronté aux mesures de calibration; puis les performances sont calculées en configuration de vol. Le bruit de fond près de Mars est évalué à ~160 c/s. L'instrument offre une bonne sensibilité à: Fe, Mg, K, Si, Th, Cl et O. Il sera possible aussi de mesurer U, Ti, H, C, S, Ca et Al. Les lobes d'émission à la surface sont aussi calculés. Ces mesures permettront une meilleure compréhension de la surface martienne. Le dernier chapitre traite de l'émission des neutrons rapides lunaires [500 keV; 10 MeV]. La forte influence de l'oxygène est mise en évidence. Comme l'a observé LUNAR PROSPECTOR, le flux intégré montre une dépendance marquée avec la teneur du régolite en fer et en titane, permettant la cartographie. L'influence des autres éléments chimiques est quantifiée. Une formule mathématique simple est proposée pour estimer le flux intégré de neutrons suivant la composition du sol. Enfin, une étude des effets de l'hydrogène sur le flux de neutrons rapides est menée; nous examinons les possibilités de quantifier son abondance dans le sol par cette méthode.

La résolution numérique de l'équation de Schrödinger en domaine extérieur nécessite l'utilisation de conditions aux limites appropriées sur la frontière du domaine de calcul. Les conditions aux limites à utiliser sont directement reliées à la fonction de potentiel intervenant dans l'équation. Pour l'équation à potentiel nul, la condition aux limites exacte est connue, ainsi que des méthodes efficaces de discrétisation et d'implémentation numérique. L'objectif de cette thèse est d'étendre les méthodes mises en jeu à potentiel nul dans le cas d'un potentiel aussi général que possible, à l'image des situations physiques variées faisant intervenir un potentiel, linéaire ou non linéaire. Nous prenons le parti de renoncer à établir des conditions aux limites exactes, au profit d'une plus grande généralité de la méthode et d'une bonne adaptation à une implémentation numérique. En se basant sur le calcul pseudodifférentiel, on propose alors une recherche détaillée de méthodes permettant de prendre en compte le potentiel dans une condition aux limites artificielle (CLA). Cette thèse traite le cas de l'équation en dimension un ou deux avec potentiel linéaire ou non linéaire, ainsi que de l'équation stationnaire en dimension un. La construction de ces CLA repose sur l'analyse microlocale et le calcul symbolique associé aux opérateurs pseudodifférentiels fractionnaires. La discrétisation en temps est effectuée à l'aide de convolutions discrètes ou d'approximants de Padé, et la discrétisation en espace repose sur des éléments finis linéaires. On utilise la méthode de relaxation de Besse pour résoudre l'équation non linéaire. L'analyse mathématique des conditions construites dans cette thèse permet de démontrer dans certains cas des estimations a priori, sur le plan continu et sur le plan semi-discret. De nombreuses simulations numériques permettent de tester l'efficacité des conditions aux limites proposées et de les comparer entre elles.