L'objet de cette thèse est l'étude de quelques aspects de <br />la notion d'hyperbolicité, plus particulièrement de la <br />relation qui existe entre celle-ci et la nature bien posée <br />d'un problème de Cauchy obtenu à partir d'un système<br />d'équations aux dérivées partielles issu de la mécanique <br />des fluides ou la réalisation de la simulation numérique <br />d'un tel problème.<br /><br />Dans un premier temps, nous rappelons en quoi la notion de<br />linéarisation d'un système d'équations aux dérivées <br />partielles semble naturelle à l'étude de ce système et <br />comment, de l'étude de ces problèmes linéarisés, plus <br />précisément de leur nature bien posée c'est-à-dire de leur <br />stabilité, découle la notion d'hyperbolicité.<br /><br />Nous étudions ensuite le cas particulier d'un modèle à <br />quatre équations pour un écoulement bifluide comportant des <br />termes de diffusion pour les équations de quantité de <br />mouvement. Nous montrons alors que, bien que, pour ce <br />système, l'ajout des termes de diffusion n'entraîne pas <br />l'hyperbolicité du modèle obtenu, les problèmes de Cauchy <br />construits à partir de la linéarisation de ce système, <br />autour d'un état constant, sont désormais bien posés.<br /><br />Enfin, nous considérons le cas d'un modèle à cinq équations <br />pour un écoulement bifluide. Ce modèle ne nécessite pas de <br />loi de fermeture algébrique (équations d'état ou lois <br />tabulées) mais comporte une équation aux dérivées <br />partielles portant sur la pression. Le système ainsi <br />obtenu n'est pas hyperbolique mais les valeurs propres de <br />l'opérateur d'advection sont toutes réelles. La simulation<br />numérique d'un écoulement régi par ce modèle, pour le cas <br />test du robinet de Ransom, ne fait néanmoins pas apparaître <br />les instabilités numériques que la nature mal posée du<br />linéarisé nous faisait craindre et qui sont présentes dans <br />les simulations réalisées à partir du modèle isentropique <br />classique à quatre équations.

Nous utilisons trois discrétisations connues pour leur localisation fréquentielle et spatiale: les bases d'ondelettes, les paquets d'ondelettes et les bases de cosinus locaux. Nous avons construit et programmé deux algorithmes: --- pour l'équation parabolique non-linéaire $\Delta(u)+\1e^(c*u)=f$ avec $f$ présentant une singularité, notre algorithme calcule la compression optimale en dimension 1 et 2, avec résultats numériques pour la dimension 1. --- pour l'équation intégrale oscillante correspondant à la Combined Integral Field Equation qui est en rapport avec le problème de diffraction des ondes (Helmholtz) par un obstacle régulier 2D, lorsque la longueur d'onde diminue vers $0$. Les trois discrétisations ci-dessus sont testées, et nous étudions sa bonne compressibité dans une analyse précise des obstacles à la compression menée de manière asymptotique. Des résultats originaux, montrant que N degrés de liberté par longueur d'onde suffisent à hautes fréquences, ont été démontrés, et les matrices résultant de ce seuillage ont été étudiées, illustrations et preuves à l'appui.

Le calcul direct du rayonnement acoustique d'une couche de mélange est réalisé en résolvant les équations de Navier-Stokes dans un domaine de calcul qui permet l'accès au champ acoustique lointain de l'écoulement. L'outil que nous avons développé utilise la technique de la Simulation Numérique Directe alliée à des méthodes numériques précises pour assurer une bonne qualité du calcul acoustique direct. Des simulations 2D et 3D de l'évolution temporelle de couches de mélange isothermes et anisothermes sont effectuées en régime subsonique et notamment à faible nombre de Mach. L'approche temporelle adoptée permet d'identifier l'émission acoustique associée à chacune des phases de la transition d'une couche de mélange. Les résultats des calculs acoustiques directs révèlent d'une part que l'émission sonore prédominante provient dans tous les cas d'un processus d'appariement et que les effets 3D réduisent globalement l'émission acoustique de la couche de mélange. D'autre part, les principales conséquences des effets combinés de la température et de la compressibilité observées expérimentalement (augmentation/réduction de l'émission sonore à nombre de Mach faible/élevé) sont reproduites avec le modèle temporel. Les données des simulations sont également utilisées pour réaliser des prédictions de bruit basées sur l'analogie de Lighthill. La comparaison de ces estimations avec les solutions de référence obtenues par le calcul acoustique direct montre un excellent accord dans tous les cas. Dans les cas anisothermes, nous montrons la nécessité de tenir compte d'un terme source supplémentaire et l'analyse des contributions séparées de chaque terme source nous permet d'interpréter le mécanisme de renversement de tendance comme le résultat d'un changement de terme source prépondérant.

Ce travail de thèse a pour but l'optimisation de forme en forgeage 3D. Les problèmes à résoudre consistent à chercher la forme optimale du lopin initial ou des outils de préforme afin de minimiser une fonction coût F qui représente une mesure de non-qualité définie par les industriels. Ce sont souvent des problèmes multi optima, et le temps nécessaire pour une évaluation de la fonction coût est très élevé (de l'ordre de la journée). L'objectif de cette thèse est de construire un module d'optimisation automatique qui permet de localiser l'extremum global à un coût raisonnable (moins de 50 calculs de la fonction coût à chaque optimisation). La simulation du procédé est effectuée avec le logiciel éléments finis FORGE3®. Les formes axisymétrique des pièces initiales ou des outils de préforme (dans le cadre du forgeage multi-passes) sont paramétrées en utilisant des polygônes quadratiques ou des courbes Bsplines. Différentes fonctions coûts sont considérées, comme l'énergie totale de forgeage ou la mesure non-qualité de la surface (défaut de repli). Le gradient de ces fonctions coûts est obtenu par la méthode de l'Etat Adjoint combinée avec la méthode de différentiation semi-analytique. Dans ce travail, afin d'aborder une famille de procédés de forgeage plus vaste, ce calcul du gradient (initié dans la thèse de M. Laroussi) a été étendu aux paramètres de forme des outils de préformes dans le cadre du forgeage multi passes.Différents algorithmes d'optimisation ont été étudiés : un algorithme BFGS standard, un algorithme de type asymptotes mobiles, une stratégie d'évolution couplée avec une surface de réponse basée sur le Krigeage et deux nouveaux algorithmes hybrides proposés dans le cadre de ce travail. Cette approche hybride consiste à coupler un algorithme génétique avec une méthode de surface de réponse pour réduire le nombre d'évaluations de la fonction coût. Tous les algorithmes étudiés sont comparés sur deux problèmes caractéristiques de forgeage 3D, respectivement l'optimisation de la géométrie de la préforme et celle des outils de préforme. Les résultats obtenus montrent la faisabilité de l'optimisation de forme en forgeage 3D, c'est-à-dire l'obtention de résultats satisfaisant en moins de 50 simulations 3D et la robustesse des algorithmes à base de méta-modèle.

De nombreuses techniques de résolution d'équations aux dérivées partielles sans maillage ont été développées dans la dernière décennie, proposant une alternative attrayante lorsque les éléments finis atteignent leurs limites. Notre travail se concentre sur l'étude de l'approximation diffuse, de ses applications au lissage et a la résolution des équations différentielles : les éléments diffus. Cependant, les solutions proposées s'appliquent aussi à d'autres méthodes et de nombreux résultats numériques illustrent chaque développement théorique. Dans un premier temps, nous étudions les techniques d'approximation sans maillage et les comparons à l'aide d'une méthode hybride. L'approximation myope est construite en modifiant le critère de construction des splines et montre ainsi la différence de nature entre moindres carres glissants et interpolateurs radiaux (krigeage, splines). Nous développons ensuite l'approximation diffuse Hermite dont l'application au calcul de la courbure des surfaces triangulées aboutit à une technique de haute résolution. Un algorithme de reverse engineering de modèle CAO s'appuie sur cette nouvelle technique et met en évidence son potentiel. L'intégration numérique est un ingrédient essentiel pour une méthode sans maillage. L'analyse du patch test nous conduit a la définition d'une technique robuste et assurant de bons rangs de convergence. Nous l'appliquons à l'approximation Hermite pour construire un modèle de poutre 3d pour le forage pétrolier. Un algorithme de contact roche-structure original est propose pour ce problème.

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.

Le travail essentiel de cette thèse réside dans l'élaboration d'un modèle mathématique permettant de calculer l'écoulement puis les pollutions thermiques ou chimiques en rivière. Les équations résolues sont les équations de Saint-Venant avec une hypothèse supplémentaire de "plafonnage" de la surface ("rivière à plafond"). La résolution utilise une méthode de différences finies avec un algorithme de pas fractionnaires. La discrétisation se fait sur un maillage curviligne à pas variable. Le modèle est conçu pour s'adapter à une géométrie quelconque sans îlot. Il a été testé de façon satisfaisante sur trois cas de figures présentant un degré de difficulté croissante.

AFIN D'OBTENIR LES PROPRIETES METALLURGIQUES OPTIMALES DU MATERIAU, LES CULASSES AUTOMOBILES EN ALLIAGE D'ALUMINIUM SUBISSENT UN TRAITEMENT THERMIQUE (TREMPE-REVENU) GENERANT D'IMPORTANTES CONTRAINTES RESIDUELLES. LEUR MODELISATION EST NECESSAIRE POUR POUVOIR PREDIRE LA TENUE EN FATIGUE POLYCYCLIQUE. LES CONTRAINTES RESIDUELLES ETANT PILOTEES PAR LA THERMIQUE DE TREMPE, IDENTIFIER L'HISTOIRE THERMIQUE DU REFROIDISSEMENT EST PRIMORDIAL. LES COEFFICIENTS D'ECHANGE H(T), MODELISANT LES TRANSFERTS THERMIQUES AVEC L'EXTERIEUR, DETERMINENT LE REFROIDISSEMENT. LEUR EVOLUTION EST DEFINIE PAR UNE FONCTION ANALYTIQUE SIMPLE DONT LES PARAMETRES SONT RECALES PAR OPTIMISATION A PARTIR DE DONNEES EXPERIMENTALES. UNE FOIS L'HISTOIRE THERMIQUE SIMULEE CORRECTEMENT, LES CONTRAINTES RESIDUELLES SONT OBTENUES SUITE A UN CALCUL MECANIQUE METTANT EN EVIDENCE, DANS LES REGIONS A RISQUE, DES ZONES DE TRACTION DEFAVORABLES A LA TENUE EN SERVICE. LA PRISE EN COMPTE DE CET ETAT INITIAL PERMET ALORS DE CARACTERISER LES ZONES CRITIQUES EXPERIMENTALEMENT OBSERVEES. PAR AILLEURS, AFIN DE RENDRE L'OUTIL OPERATIONNEL EN BE, IL EST NECESSAIRE DE PROPOSER UNE METHODE NUMERIQUE PERMETTANT DE DIMINUER LES DUREES DE SIMULATION POUR LA RESOLUTION D'UN PROBLEME D'EVOLUTION ET POUR LA RECHERCHE DIRECTE DE LA REPONSE STABILISEE D'UNE STRUCTURE SOUMISE A UN CHARGEMENT CYCLIQUE. L'APPROCHE DEVELOPPEE REPOSE SUR LA METHODE A GRAND INCREMENT DE TEMPS ET LA METHODE CYCLIQUE DIRECTE : L'EQUILIBRE GLOBAL, ECRIT SOUS SA FORME RESIDUELLE, EST RESOLU SUR UNE BASE D'ONDELETTES REDUITE. DES GAINS SIGNIFICATIFS EN TEMPS CPU SONT OBTENUS POUR LA RECHERCHE DU CYCLE STABILISE.

Ce travail de thèse est consacré à la modélisation et à la simulation numérique des transitions de phase liquide-vapeur. L'étude effectuée se découpe en deux randes parties: une première où on étudie les phénomènes de transition de phase avec une loi d'état de type Van Der Waals (perte de monotonie de la loi d'état) et une deuxième partie où on choisit une approche alternative avec deux loi d'états. La première partie consiste à étudier les critères visqueux classiques de sélection des solutions du système d'équations utilisé lorsque la loi d'état n'est pas monotone. Les critères classiques ne sélectionnant pas des solutions a priori physiques, un critère plus récent est introduit: le critère visco-capillaire. L'utilisation de ce critère avec un solveur de Riemann exact (sous la contrainte de trouver le zéro d'une fonction non linéaire) permet d'obtenir des résultats mais avec un coût de calcul trop élevé. Une approche alternative est alors envisagée avec deux lois d'états (une pour chaque phase). A l'aide d'un procédé de minimisation de l'action hamiltonienne, un modèle bifluide de changement de phase est proposé. Celui-ci respecte alors le second principe de la thermodynamique. Deux sous-systèmes en sont déduits à l'aide d'un procédé de retour à l'équilibre: mécanique dans un premier temp puis mécanique et thermodynamique dans un deuxième temps. Malgré la faible hyperbolicité du dernier sous-système obtenu, des schémas numériques stables basés sur une méthode de splitting sont proposés. On montre alors que le système ainsi obtenu est naturellement capable de nucléer des bulles de vapeur dans du liquide.

La circulation de l'information devient une préoccupation prioritaire des sociétés modernes. C'est pourquoi, la production audiovisuelle destinée aux télévisions et aux entreprises doit considérablement se développer dans les prochaines années. Une concurrence sévère nait dès maintenant de cette situation. Elle oppose les entreprises artisanales, qui jusqu'à ce jour occupaient le marché, et les sociétés industrielles qui se créent pour répondre à la demande prévisible. La vie et le savoir-faire d'un secteur d'activité sont menacé. Les petites entreprises, qui constituent le potentiel de production de dessin animé de notre pays, doivent s'adapter à une nouvelle situation internationale. En effet, la Grande-Bretagne, le Japon et les USA ont pris une place importante dans le domaine de la série télévisée de dessins animés. Ils ont conquis le marché en proposant leurs séries à la fois plus vite et moins cher. Contrairement à une idée répandue, ils n *utilisent pas encore systématiquement les techniques informatiques. Ils ont porté 1 *économie sur les frais de personnels et la simplicité de réalisation en mettant en oeuvre des mét~odes classiques d'industrialisation d'une fabrication: parcellisation du travail, normes de qualité souvent respectées à leur niveau minimum ... Mais la force commerciale et l'expérience acquises leur permettent désormais de proposer des productions de qualité croissante et d'intégrer progressivement des techniques modernes. dessin animé. Cette base de travail a été établie en collaboration étroite avec des professionnels du domaine. Nous l'avons complétée par une esquisse des apports actuels ou potentiels de l'informatique. Nous terminons en abordant les problèmes techniques que cé travail contribue à résoudre. Les principales références qui ont servi a établir ce chapitre sont: [BRI84], [MAR70], [BCT84]. Cette description ne concerne que le dessin animé classique le cartoon. Nous n'avons pas explicité ici les mille et une façons de faire du cinéma d'animation. Citons pour mémoire l'animation de papier découpé, de marionnettes, la pixilation. Face à la demande fortement croissante de dessins animés, destinés en particulier aux programmes de télévision, il devient nécessaire d'améliorer la productivité dans ce domaine. Outre des méthodes industrielles classiques, comme la parcellisation des tâches, des apports de nature différente peuvent être espérés de l'informatique: gestion de l'information (flux, mises à jours, cohérences entre documents), automatisation de tâches répétitives, outils interactifs de définition numérique d'images et de mouvements. Des recherches tendent à faire de l'ordinateur un outil nouveau de créa ti on d *images. animées. Not a mme nt, des études s*ont faites pour définir des objets avec des caractéristiques: élasticité, poids, ... et les manipuler en temps réel sur un écran d'ordinateur (formes du langage Logo, [LUC84]). D'autres études ont plutôt porté sur l'imitation du processus traditionnel de production de dessin animé. Des travaux ont été entrepris dans ce sens depuis plus de dix ans. Les plus grands studios: Hanna-Barbera, Walt Disn~y ont investi dans ces recherches. Aujourd'hui les systèmes informatiques d'aide à la production de dessins animés commencent seulement à sortir des laboratoires pour être utilisés en production [WAL81] [DIS84]. La fabrication des dessins animés est une activité créatrice nécessitant un grand travail manuel. Les trois points clés qui déterminent la concurrence entre les studios sont le coût, les délais et les volumes de production. Il sera intéressant de gagner sur chacun de ces points.