Ce travail, qui a été réalisé dans le cadre d'un projet européen, avait pour objectif de faire la simulation numérique de l'injection d'un alliage d'aluminium (A356) à l'état semi-solide. Pour atteindre cet objectif, nous avons travaillé sur deux principaux axes qui sont l'expérience et la simulation numérique. Sur le plan expérimental, l'objectif était de trouver une loi de comportement et ses paramètres pour caractériser le comportement de notre alliage à l'état semi-solide. Cette étude nous a conduit d'une part à passer en revue les différentes classes de loi disponibles dans la bibliographie et à utiliser une loi de comportement viscoplastique, et d'autre part à réaliser un nouveau test rhéologique, le test de l'écoulement de Stephan. Les résultats expérimentaux ont été satisfaisants et l'identification des paramètres a pu être faite, en utilisant une stratégie de type essai-erreur. Ceci nous a permis de trouver pour notre cas des paramètres en concordance avec ceux trouvés dans la bibliographie et d'obtenir un accord raisonnable entre simulation numérique et mesures sur l'ensemble des essais effectués. Sur le plan purement numérique, nous avons dans un premier temps adapté le code de calcul R3 à la mise en forme à l'état semi-solide. Ce travail nous a conduit à y introduire un module de contact matière/matière et à développer le module de contact matière/outil existant. Nous avons alors été confrontés au problème de préconditionnement de la matrice hessienne, dont le conditionnement est fortement dégradé par l'ajout de la condition de non-interpénétration du contact matière/matière. Ce problème a été résolu grâce à l'utilisation du préconditionneur par factorisation incomplète de Crout. Nous avons ensuite travaillé sur la formulation eulérienne lagrangienne arbitraire et sur la procédure de remaillage automatique. Malheureusement, ces travaux ne nous ont pas permis de pouvoir simuler jusqu'au bout le remplissage de pièces industrielles. Nous sommes alors passés dans un deuxième temps à une formulation eulérienne en utilisant le code de calcul Rem3D, initialement développé pour l'injection de thermoplastiques. Nous y avons introduit les termes d'inertie, ainsi que notre loi de comportement. Ce dernier code de calcul nous a alors permis de simuler l'injection de notre alliage d'aluminium à l'état semi-solide. Les résultats obtenus montrent une bonne concordance avec le procédé industriel et le logiciel peut constituer un outil précieux d'aide à la conception.

Les explosions sphériques entraînent des perturbations importantes de l'interface entre les produits de détonation et l'air. Ces instabilités jouent un rôle dominant dans la détermination du volume de la "boule de feu". Un calcul sphérique unidimensionnel classique conduit un volume de sphère très inférieur à celui mesuré expérimentalement. De plus, des réactions de post-combustion peuvent avoir lieu dans la zone de mélange, libérant une énergie deux fois supérieur à celle de la détonation, déjà considérable. À une échelle suffisamment petite, on distingue les longueurs d'onde des instabilités et les tailles de jets, mais à une échelle plus globale, on observe une couche de mélange où la forme précise de l'interface n'est plus visible. Les deux phases (produits détonation et de l'air) s'interpénètrent, et par conséquent, l'interface devient une zone de mélange. Pour calculer correctement chacune des instabilités, une approche multidimensionnelle semble s'imposer. Cependant, un grand nombre de cellules est nécessaire pour calculer une structure unique de la zone de mélange. En outre, pour une instabilité isolé, le maillage entraînent des instabilités parasites qui dépendent fortement de la viscosité numérique du schéma utilisé. L'approche multidimensionnelle, basée sur la simulation numérique directe, présente donc des difficultés. En réalité, nous ne voulons pas calculer la forme exacte des instabilités de l'interface, mais seulement l'épaisseur de la couche de mélange et les champs de concentrations des phases dans celle-ci. Ainsi, une approche unidimensionnelle peut être suffisante. L'objectif est d'écrire un modèle unidimensionnel décrivant le phénomène d'interpénétration. Trois modèles ont alors été construits à partir du modèle diphasique de l'Baer et Nunziato (1986). Nous obtenons des résultats intéressants avec les deux premiers sur des problématiques d'épaississement d'interface, mais ils sont insuffisants. Le dernier modèle, qui dérive des deux premiers, a été validé sur des tests d'explosions sphériques.

L'objectif est d'améliorer la stabilité et la précision de la discrétisation spatiale de type volumes finis sur des maillages non structurés. La thèse fournit une analyse générale de la reconstruction des polynômes de degré k en maillage non structuré et présente plusieurs algorithmes permettant de reconstruire des polynômes sur de petits voisinages compacts. Une étude théorique de la stabilité établit des principes pour concevoir des méthodes de reconstruction stables. Une étude théorique de la précision caractérise les erreurs induites par le maillage non structuré à l'aide de l'approche de l'équation modifiée. L'étude formule également des algorithmes de limitation en maillage non structuré basés sur une approche géométrique. Toutes les études théoriques sont complétées par des expériences numériques. Les calculs LES d'un écoulement subsonique au-dessus d'une cavité et d'un jet supersonique permettent de valider et comparer plusieurs options de discrétisation spatiale.

On étudie la pulvérisation numérique de jets et de nappes liquides. Dans la première partie une revue bibliographique du problème présente la théorie des instabilités de Kelvin-Helmholtz et différents résultats expérimentaux antérieurs. Dans le seconde partie on présente différentes méthodes numériques : technique de suivi d'interface (Volume of Fluid), de calculs des efforts de tension de surface (CSF) et de résolution des équations de Navier-Stokes (Mac Cormack). La dernière partie décrit une méthode particulaire (SPH) permettant de traiter des interactions d'écoulements liquides et gazeux. Les principaux phénomènes interfaciaux, tels que la tension de surface et la vaporisation sont pris en compte. Des calculs 2D et 3D de déformation de jets liquides dans les écoulements gazeux supersoniques sont analyses et comparés avec des résultats expérimentaux.

La génération de fréquences discrètes dans l'écoulement autour d'un profil d'aile est étudiée dans cette thèse au moyen de simulations numériques non-linéaires et d'études de stabilité globale. À cette fin, un code numérique mettant en œuvre une nouvelle technique pour accéder à la dynamique linéaire, directe et adjointe, a d'abord été développé et ensuite appliqué à l'écoulement autour d'un profil. Les simulations non-linéaires confirment l'apparition de fréquences discrètes dans le spectre sonore, et pour le cas considéré, l'analyse de stabilité globale de la dynamique linéaire de l'écoulement moyen montre que celle-ci est stable. Cependant, la réponse de l'écoulement à des perturbations incidentes révèle de fortes croissances transitoires amenant à l'établissement de cycles de rétroaction aéroacoustique. Ces cycles comprennent la croissance des instabilités hydrodynamiques dans les couches limites sur l'intrados et l'extrados ainsi que leurs interactions avec les radiations acoustiques du bord de fuite. Les processus associés à l'apparition des fréquences discrètes dans le spectre sonore ainsi que les cycles de rétroaction sont ensuite mis en relation avec les modes globaux les moins stables: d'un coté, la structure spatiale des modes directs montre la croissance des instabilités hydrodynamiques sur l'extrados et la zone du sillage proche du bord de fuite; d'un autre coté, les modes adjoints associés présentent l'intrados comme la zone la plus réceptive à des perturbations externes. Finalement, l'analyse de la région dite du wavemaker indique, en accord avec les expériences, le rôle fondamental de la couche limite sur l'intrados.

L'objet de cette thèse a été de développer un code de calcul pour la simulation d'écoulements diphasiques de fluides non miscibles, incompressibles et isothermes afin de l'appliquer au phénomène de fragmentation d'un jet liquide, et plus particulièrement à la désintégration assistée d'une nappe de liquide par deux écoulements d'air portés à haute vitesse.<br /><br />Les choix des hypothèses physiques et des schémas numériques associés ont été fait afin de respecter au mieux la physique complexe de brisure d'un jet liquide. La résolution des équations de Navier-Stokes incompressibles est faite de façon directe par le biais d'une méthode de projection. La méthode naissante et prometteuse Level-Set assure quant à elle le suivi numérique de la surface de séparation entre deux fluides non miscibles. Enfin, la méthode Ghost Fluid permet un traitement correct des discontinuités à la traversée de l'interface en préservant au niveau discret les conditions de saut entre équilibre des forces de capillarité, de pression et de viscosité. Les bonnes aptitudes de tels schémas numériques sont montrées à travers une batterie de cas tests académiques.<br /><br />Les mécanismes physiques mis en jeu lors de l'atomisation primaire d'un jet liquide ont fait l'objet de nombreuses études théoriques et expérimentales. Les instabilités se développant à la surface du liquide sont multiples et clairement tridimensionnelles. Actuellement et malgré beaucoup d'efforts de recherche, aucune théorie ou modèle ne rend compte rigoureusement de ce phénomène. Or, il est montré que l'outil numérique peut apporter une lumière nouvelle. Par exemple, l'impact de la couche limite gazeuse sur la fréquence d'oscillation est étudié. Nous avons également pu retrouver par la simulation la dynamique ligamentaire intrinsèque à la désintégration assistée d'un jet liquide.

Les polluants et les matières qui transitent dans les cours d'eau peuvent constituer des stocks dans la colonne sédimentaire. Ces dépôts, soumis aux processus d'autoconsolidation, sont sujets au vieillissement et à l'enfouissement, qui augmentent leur temps de résidence dans la colonne sédimentaire par augmentation de leur résistance à l'érosion. Les flux de polluants peuvent être ramenés aux flux de matières sédimentaires, dont l'érosion est décrite par des lois à seuil. Nous proposons une modélisation du seuil d'érosion d'un sédiment cohésif. Ce modèle considère les propri étés cohésives des matières qui composent les couches sédimentaires, en s'appuyant sur la force attractive électrodynamique de contact de Van der Waals. Cette force, que nous avons paramétrée par la granulométrie et la porosité, vient s'ajouter au bilan des forces généralement considéré pour modéliser le seuil d'érosion. Des expériences d'érosion artificielle ont été menées dans un canal à érosion, avec des sédiments naturels prélevés dans la partie en aval du Rhône, pour en évaluer le seuil d'érosion, la granulométrie et la porosité. Les résultats des expériences alimentent la littérature en données expérimentales sur l'érosion de sédiments cohésifs consolidés. Le modèle d'érosion ainsi développé est appliqué dans les conditions expérimentales. Il donne un bon accord avec les mesures. Nous illustrons sa prise en compte dans un modèle de consolidation, permettant de construire des diagrammes du type Mohr pour l'érosion d'un sédiment cohésif faiblement consolidé.

Pour étudier l'influence de la gravité sur la croissance par T.H.M. du semi-conducteur CdTe dopé indium (programme EURECA), nous avons effectué une analyse comparative d'expérience spatiales et terrestres, complétée par une simulation numérique de l'hydrodynamique de la phase fluide utilisant un code 3D. Notre deuxième objectif a été de prépérer des monocristaux de CdTe dopé vanadium par la technique Bridgman présentant de bonnes propriérés photoréfractives. Pour une meilleure compréhension du mécanisme de l'effet photoréfractif, les matériaux ont été caractérisés par R.P.E., D.C.M. et P.I.C.T.S. afin de déterminer le degré de valence de V et la signature des niveaux vanadium dans la bande interdite. Des mesures optiques (absorption, photoconductivité...) et d'effet photoréfractif ont été effectuées. Enfin, pour évaluer la pertinence spatiale de cette croissance, nous avons étudié par simulation numérique l'influence de la gravité sur la configuration hydrodynamique de la phase fluide.

L'étude physique des enveloppes externes de la Terre est fortement influencée par leur stratification. Cette dernière permet la propagation d'ondes anisotropes appelées ondes de gravité internes. L'interaction de ces ondes avec la topographie gouverne la génération et la dissipation de ces ondes, et joue un rôle non négligeable dans le bilan énergétique global de la planète.<br />Nous avons mis au point une technique de visualisation des ondes internes par strioscopie synthétique, et mené parallèlement des expériences à grande échelle sur leur réflexion critique. Nous avons également élaboré un nouveau type d'excitateur d'ondes internes planes.<br />Ces expériences nous ont permis de valider une théorie analytique d'émission par des corps oscillants et de l'appliquer à la génération d'ondes par un talus continental. Nous avons également confirmé la présence d'effets nonlinéaires lors de la réflexion critique. Enfin, nous avons étudié l'influence de la force de Coriolis sur ces phénomènes.

De nombreuses hypothèses sont classiquement utilisées pour décrire le comportement du fluide dans un contact lubrifié : film continu, viscosité constante dans l'épaisseur, film mince, fluide newtonien... Or, certaines s'avèrent erronées dès lors que l'on s'intéresse aux contacts Elasto- HydroDynamiques fortement glissants ou à la répartition du lubrifiant en sortie de contact. Une approche numérique originale, basée sur un retour aux équations de la mécanique des fluides générale et prenant en compte le couplage fluide/solide et les effets thermiques sont proposés ici dans le but d'apporter des éléments physiques supplémentaires aux modélisations usuelles. Dans un premier temps, l'influence des effets thermiques sur l'évolution du frottement dans les contacts Thermo-EHD est mise en évidence. La présence d'un minimum de frottement pour le cas du glissement pur est expliquée par l'analyse des transferts thermiques entre le lubrifiant et les solides. L'origine des modifications locales d'épaisseur de film observées et l'existence même d'une épaisseur de film lubrifiant pour les cas de vitesse d'entraînement nulle sont alors reliées à la présence d'un fort gradient de viscosité dans l'épaisseur de film. Une comparaison qualitative avec des éléments expérimentaux de la littérature est réalisée, validant les tendances obtenues. Dans un second temps, l'écoulement à surface libre du lubrifiant en périphérie du contact est étudié expérimentalement puis numériquement par une méthode à interface diffuse. Le rôle des effets capillaires est analysé et les résultats numériques confrontés à des résultats issus de la littérature. Un bon accord est obtenu tant qualitativement que quantitativement. Validé par l'étude numérique diphasique (air/lubrifiant) réalisée, un modèle analytique simplifié est alors développé, prédisant une loi de répartition du lubrifiant en sortie de contact. La zone de sortie des contacts EHD est ensuite traitée par un modèle de cavitation vaporeuse et la prise en compte nécessaire de l'air environnant est discutée. Enfin, une première modélisation tridimensionnelle de l'écoulement à surface libre du lubrifiant autour d'un contact ponctuel est réalisée mettant en avant l'influence des effets capillaires et la faisabilité d'une telle approche.