Le passage de la méso-échelle, échelle des modèles de circulation générale GCM, de l'anglais "General Circulation Models", à la micro-échelle (échelle hydrologique), pour les précipitations, est un exercice assez complexe. Les champs de précipitations comme la plupart des champs géophysiques turbulents obéissent au concept d'invariance d'échelle, qui est une caractéristique principale des champs multifractals. Par ailleurs, il a été prouvé que le transfert d'énergie des grosses structures aux plus petites structures au sein d'un phénomène géophysique turbulent s'effectue de façon multiplicative (Kolmogorov, 1962; Mandelbrot, 1974 ...): un facteur aléatoire déterminant la fraction de flux transmis d'un gros tourbillon à un plus petit. Le travail que nous présentons ici s'inscrit dans le cadre du projet EDF "Prévisions saisonnières et Hydraulicité" dans la gestion de son parc hydroélectrique et a pour objectif la construction d'un modèle de désagrégation basé sur le principe d'invariance d'échelle des champs de précipitation, donc utilisant les propriétés des champs géophysiques mentionnées ci-dessus. Dans un premier temps, nous conduisons une analyse multifractale (Schertzer et Lovejoy , 1991) sur des séries pluviométriques de la France (243 séries pluviométriques au pas de temps de six minutes, constituées sur une dizaine d'années distribuées sur la France métropolitaine), ce qui nous permet de déduire les paramètres multifractals, dans le temps, dans l'espace ou dans le cas spatio-temporel. La seconde étape consiste à construire des cascades multifractales, à partir des valeurs saisonnières de pluies avec les paramètres déterminés dans la première étape. Le principe de cette deuxième partie consiste, à partir d'une prévision mensuelle sur des mailles de dimensions 243km×243km×32jours (correspondant à une anisotropie espace-temps de l'ordre de H=2/3 :x=y=t (3/2) ) voisines de celles des modèles de circulation générale (dimensions de l'ordre de 250km×250km×30jours) et à conduire la cascade multifractale, avec les paramètres multifractals préalablement déterminés, pour atteindre des valeurs de prévision sur des mailles de l'ordre de 1km×1km×1j. Les résultats obtenus devront faire l'objet d'un conditionnement orographique avant d'être comparés avec les valeurs réelles obtenues.

Dans cette étude nous considérons l'ébullition en paroi sous les angles de la modélisation et de la simulation numérique. Dans un premier temps nous proposons une revue bibliographique au sujet du régime d'ébullition nucléée à fort ux de chaleur pariétal et analysons plus particulièrement la compréhension du phénomène de crise d'ébullition de type caléfaction. Nous en déduisons une motivation pour l'étude de la dynamique de croissance de bulle au moyen de la simulation numérique. L'essentiel du travail concerne alors le développement d'un modèle de type champ de phase pour l'étude des écoulements liquide-vapeur avec changement de phase. Nous proposons une fermeture thermodynamique quasi-compressible dont les propriétés sont adaptées aux simulations envisagées. Le système d'équations du mouvement qui s'en déduit constitue une régularisation thermodynamiquement cohérente de la description discontinue du système diphasique, ce qui est l'avantage des modèles à interfaces diffuses. Nous démontrons que la formulation retenue permet de définir l'épaisseur de la zone régularisée indépendamment de la description thermodynamique des phases, ce qui est intéressant numériquement. Nous établissons la relation cinétique et analysons ainsi la modélisation champ de phase des mécanismes dissipatifs. Finalement nous étudions la résolution numérique du modèle à l'aide de simulations de transition de phase en configuration simple et de dynamique de bulle en isotherme.

Nous étudions dans cette thèse le problème de gestion du risque physique en production électrique à l'horizon hebdomadaire. Dans un premier temps, nous nous intéressons à l'intégration de l'aléa d'apport hydraulique dans la gestion locale d'une vallée hydraulique. Cette approche est menée à l'aide de l'optimisation robuste et de règles de décision linéaires. Les résultats de multiples modes de simulation montrent que ces approches permettent de réduire notoirement les déversés en comparaison des modèles déterministes appliqués en exploitation, moyennant une faible augmentation du coût. La deuxième problématique traitée est la gestion active de la marge de production, définie comme l'écart entre l'offre totale et la demande totale, compte tenu des aléas affectant le système électrique. Il s'agit de déterminer quelles décisions optimales prendre, selon un certain critère économique, pour se couvrir contre un risque trop élevé de non-satisfaction de la demande dans au moins 99% des situations. Pour cela, une formulation inédite en boucle ouverte, basée sur le processus stochastique de marge de production et des contraintes en probabilité est proposée. Pour les besoins de cette formulation, nous générons des scénarios à l'aide de techniques plus réalistes qu'en exploitation. Enfin, une résolution moins anticipative est étudiée en appliquant l'heuristique «Programmation Stochastique avec Règles de Décision Constantes par Morceaux» introduite par Thénié et Vial. Les premiers résultats sont très encourageants en comparaison avec les modèles en boucle ouverte.

On a souvent recours à des assemblages métal-céramique pour des systèmes complexes couplant les propriétés, souvent opposées, des alliages métalliques et des céramiques. Les besoins industriels pour de hautes températures d'utilisation, en milieu corrosifs, restent toujours difficiles à satisfaire de façon fiable. Deux céramiques techniques industrielles ont été retenues pour cette étude : le carbure de silicium SiC et l'alumine Al2O3. Elles ont été assemblées à un superalliage de Nickel réfractaire (HAYNES™214®), selon le procédé de thermo-compression sous vide, par l'intermédiaire d'un joint métallique de faible épaisseur. La forte réactivité du SiC avec le nickel a conduit à rechercher une barrière de diffusion entre l'alliage et la céramique. Un joint d'argent ne formant aucun siliciure a permis d'éviter cette réactivité et a conduit ainsi à des tenues mécaniques en cisaillement significatives (45 MPa). Pour l'alumine (Al2O3), nous nous sommes orientés vers l'utilisation de joints de nickel, peu réactif avec cet oxyde. Cependant la différence de coefficient de dilatation thermique entre la céramique (Al2O3) et les métaux (Ni et HAYNES™214®) conduit à de fortes concentrations de contraintes résiduelles à l'interface. Une simulation par éléments finis, en utilisant le code ABAQUS, a permis de localiser le champ de contraintes résiduelles dans ce type d'assemblage. Une comparaison entre un calcul élasto-plastique et un calcul élasto-visco-plastique a été réalisée. La distribution des contraintes est sensible à l'épaisseur du joint de nickel ainsi qu'aux dimensions et formes respectives de l'alliage et de la céramique. Il a ainsi été montré que le calcul élasto-visco-plastique est indispensable pour caractériser la liaison Al2O3/Ni/HAYNES™214®, contrairement à des liaisons Al2O3/Ni/Al2O3 étudiées plus classiquement. De plus, les contraintes de tension dans la céramique sont beaucoup plus importantes pour l'assemblage alumine/alliage métallique que pour l'assemblage alumine/alumine. Ceci est retrouvé expérimentalement par les résultats de caractérisation mécanique (cisaillement, délamination) de ces deux types d'assemblage et par des mesures expérimentales de contraintes résiduelles par indentation et par diffraction des rayons X. La fragilisation de l'alumine près de l'interface est un phénomène connu dû essentiellement à la diffusion du métal de liaison et au frottement entre le joint métallique et l'alumine. Une corrélation entre la ténacité, l'intensité des contraintes résiduelles et le piégeage de charges électriques a été montrée. Enfin, l'étude de joints multicouches Cu-Ni-Cu a permis d'accommoder les contraintes et d'améliorer la résistance de l'assemblage Alumine/HAYNES™214®. De plus, la préoxydation des joints de cuivre côté alumine, reproduisant ainsi la méthode de l'eutectique, ouvre une voie quant à une utilisation des assemblages pour de hautes températures (> 1200 °C) et l'industrialisation du procédé pour de grandes séries.

Les matériaux composites occupent une place de plus en plus importante dans l'industrie, en particulier aéronautique. Les composites à matrice thermoplastique suscitent aujourd'hui dans ce secteur un très fort engouement par rapport à leurs concurrents thermodurcissables. Le présent travail de thèse se concentre sur le soudage de matériaux composites à matrice thermoplastique. Le procédé de soudage ultrasons consiste à dissiper un travail mécanique au niveau de l'interface. La température augmente alors localement et entraîne un écoulement dans la zone soudée. Le procédé dit "continu", proposé par EADS IW, permet aujourd'hui d'envisager l'assemblage de structures aéronautiques à l'échelle industrielle par la technique ultrason. Bien évidemment, une bonne compréhension du procédé est un préalable indispensable à une telle application aéronautique. Afin d'étudier la qualité de l'adhésion, la présente étude se situe à l'échelle dite mésoscopique, de quelques millimètres. La problématique consiste à modéliser puis simuler numériquement la thermo-mécanique de l'écoulement à l'interface lors du soudage. L'objectif est de mieux comprendre les liens entre les différents paramètres du procédé, les phénomènes physiques dans la zone soudée et la qualité de la soudure. Après identification des phénomènes physiques prépondérants, la mise en équations de cet écoulement révèle la co-existence de phénomènes rapides (vibration) et lents (écoulement). Une méthode d'homogénéisation temporelle par développements asymptotiques permet alors de modéliser le procédé à l'aide de trois problèmes aux limites couplés. La résolution numérique de ces problèmes nécessite l'utilisation de méthodes particulières permettant de gérer les couplages multiphysiques et les grandes déformations de la géométrie. Finalement, l'outil numérique développé permet une analyse de l'influence des paramètres procédés sur la qualité de l'écoulement et du soudage.

De très nombreux problèmes rencontrés dans l'industrie sont liés aux évolutions que subissent les interfaces. Le contact, le frottement et l'adhésion y jouent un rôle fondamental. Ce travail est une contribution au développement d'un modèle qui couple adhésion et frottement et qui autorise la cicatrisation des liens adhésifs lorsque les corps sont remis en contact après avoir été séparés. Ce modèle est appelé modèle d'adhésion cicatrisante. Dans le cadre de ce mémoire, les applications concernent plus particulièrement l'étude du contact verre/élastomère. La formation du col d'adhésion, les phénomènes « jump-on » et « jump-off » et la formation et la propagation des ondes de Schallamach sont abordés. Le travail propose une modélisation à l'échelle des surfaces, issue de considérations thermodynamiques et qui simule les effets des interactions de type van der Waals. Outre la possible cicatrisation des liens lors de l'approche des corps, l'originalité du modèle est que les processus d'endommagement et de cicatrisation des liens adhésifs sont tous deux liés à la vitesse de sollicitation au niveau de l'interface. Selon les applications envisagées, l'approche permet ainsi de considérer deux sources possibles de dissipations. L'une est surfacique et liée au comportement de l'interface. La seconde est volumique et liée au comportement des corps. La formulation dynamique est donnée dans le cadre de la mécanique non régulière. Les méthodes numériques correspondantes sont mises en œuvre dans la plateforme LMGC90 (Montpellier - Marseille). Le modèle est testé sur des Benchmarks et validé sur la simulation d'une expérience de tribologie de contact verre/élastomère.

Ce travail se situe dans le cadre de la modélisation numérique par éléments finis du procédé de forgeage à chaud des métaux. Il a pour but l'optimisation du logiciel de simulation numérique tridimensionnelle Forge 3. Il s'oriente autour de deux points bloquants pour une utilisation de ce logiciel: des problèmes de précision sur la surface de contact entre la pièce et l'outil et un temps de calcul élevé. Dans la première partie, on étudie le problème du contact unilatéral. En partant de la méthode de pénalisation existant, on améliore l'algorithme de contact par la mise en place d'un algorithme implicite. On valide ces développements sur plusieurs cas de forgeage où l'on montre l'apport du nouvel algorithme en terme de précision de gestion de la surface entre la pièce et l'outil. Dans la seconde partie, on traite du problème de temps de calcul. Malgré la mise en place d'une méthode de calcul parallèle les délais de calcul pour une simulation tridimensionnelle restent très élevés. Nos travaux ont donc consisté à remplacer le solveur de type résidu minimal implémenté dans Forge 3 par une méthode asymptotiquement optimale basée sur l'utilisation de multigrilles. Nous avons développé une méthode originale utilisant la méthode de résidu minimal comme lisseur et une méthode de raffinement de maillage pour la construction des différents niveaux. Nous avons optimisé cet algorithme en introduisant en particulier un préconditionneur utilisant une factorisation de Crout incomplète pour les résolutions grille grossière. Les résultats obtenus en bigrille sur différents cas de forgeage sont excellents et les quelques résultats sur un calcul avec trois niveaux de maillage sont très encourageants.

Cette thèse présente des solutions analytiques en conjugaison de phase acoustique unidimensionnelle dans un milieu actif. Le conjugateur dans ces études est en contact direct avec des milieux passifs et la zone active peut occuper la totalité ou une partie du conjugateur. <br /><br />La base de la théorie développée dans cette thèse est une analyse directe par une méthode multi-échelle. Les résultats sont obtenues sous une formes simple à utiliser : le seuil du mode supercritique, les taux de croissance et sa relation avec les répartitions spatiales de contrainte dans le conjugateur sont donnés sous forme analytique. <br /><br />Une approche du problème de focalisation est présentée par transposition des résultats 1D à la symétrie sphérique. Toutes les solutions analytiques sont comparées avec des résultats numériques. Enfin le problème de la conjugaison de phase en 2D est initié numériquement.

Ce travail porte sur la production de son par interaction entre un écoulement et un résonateur : le couplage est affuté de manière importante par le détail géométrique de la lumière comme en témoigne la facture de la flûte à bec, dans laquelle des chanfreins sont taillés à la sortie du canal. Ce point de facture est considéré comme étant crucial, tant pour le timbre et la stabilité de l'instrument que pour le confort de jeu de l'instrumentiste. Des visualisations de naissances de jets à la caméra rapide complétées par des simulations numériques révèlent que la géométrie de la sortie du canal a une influence sur le comportement du jet faible par rapport au nombre de Reynolds de l'écoulement, et indiquent que le rôle joué par les chanfreins se situe dans l'interaction entre le jet et le champ acoustique issu du résonateur de l'instrument. Un modèle analytique simple de réceptivité d'un jet soumis à une perturbation transverse est développé, tenant compte de la géométrie de la lumière. Ce modèle est étudié dans le cas d'un oscillateur simple. Le comportement d'un instrument dont la géométrie du système excitateur peut être changée est étudié à la bouche artificielle, et l'adaptation du contrôle exercé par un instrumentiste est considérée. Enfin, l'étude d'une famille de flûtes révèle que le facteur adapte la géométrie de l'excitateur dans chaque instrument dans deux buts : d'une part le spectre des notes des flûtes les plus graves est enrichi. D'autre part, les pressions de jeu sont du même ordre de grandeur pour tous les instruments.

Cette HDR, présente les travaux de simulation et de modélisation de la combustion turbulente des écoulements diphasiques. la simulation numérique directe est la base de ces travaux. Un code de DNS a été couplé avec un solveur Lagrangien pour décrire une phase liquide dispersé s'évaporant dans le domaine. Tous les aspects présents dans les chambres de combustion sont abordés : dispersion, évaporation, mélange et combustion.