Le mémoire de maîtrise présente la dynamique de rotor d’une microturbine à gaz à configuration renversée permettant l’intégration de pales en céramique. Cette configuration permet d’utiliser des pales de turbine complètement en céramique car elles sont retenues en compression de l’extérieur par un anneau en composite qui est en rotation. Par contre, la configuration renversée nécessite une nouvelle conception de moyeu de turbine. En effet, l’anneau en composite est loin de l'axe de rotation. Par conséquent, il se dilate beaucoup dus aux efforts centrifuges ce qui entraîne un grand déplacement des pales. L’idée est donc de faire un moyeu flexible qui vient suivre le déplacement des pales pour appliquer une pression en dessous de ceux-ci, afin de les tenir en place par friction. La dynamique de rotor avec moyeu flexible est très peu documentée dans la littérature. Le projet de recherche consiste en la compréhension de la dynamique d’un moyeu flexible pour cibler les paramètres intervenant dans la vibration d’une turbine renversée. La recherche se fera via la conception d’un prototype permettant d’étudier l’effet de la flexibilité pour un concept donné. L’analyse de dynamique de rotor avec moyeu flexible est faite avec un modèle par éléments finis linéaire. Malgré le fait que des non-linéarités peuvent provenir des roulements à billes supportés par des films amortisseurs, un type de roulements couramment utilisé pour supporter des turbocompresseurs, il est démontré qu’un modèle linéarisé capte plutôt bien les tendances mesurées. Le but de ce projet est de déceler les paramètres problématiques d’un moyeu flexible et de découvrir les paramètres permettant une bonne dynamique.

Le développement du trafic aérien autour des agglomérations est limité notamment par les nuisances sonores qu'il engendre. Il est donc indispensable de parvenir à comprendre l'origine de ces phénomènes afin d'en réduire les effets. Les hélicoptères, qui sont amenés à survoler les agglomérations à basse altitude, font ainsi l'objet de nombreuses études. Le but principal de notre étude est de développer des méthodes numériques permettant d'évaluer le rôle des instabilités elliptiques de tourbillons dans la production du BWI. La première partie de cette thèse s'intéresse aux caractéristiques du bruit produit par l'interaction d'une pale avec un tourbillon instable. Pour ce faire, une simulation des instabilités elliptiques se développant dans une paire de tourbillons analytiques corotatifs parallèles rectilignes uniformes est d'abord réalisée. Ensuite, un outil numérique, basé sur la théorie d'Amiet, est développé pour obtenir les perturbations de pression sur pales issus d'une interaction entre une pale et le tourbillon instable. Un bon accord est trouvé entre ces résultats et les données expérimentales. Dans une deuxième partie nous nous intéressons à l'influence des caractéristiques géométriques du tourbillon d'extrémité de pale jeune isolé (coeur du tourbillon elliptique, trajectoire courbe) sur le développement des instabilités elliptiques. Pour cela le profil d'un tourbillon d'extrémité de pale a été calculé au moyen du code volumes finis. Une étude linéaire a ensuite été menée en considérant cette solution comme champ porteur, au moyen d'une méthode Euler linéarisée harmonique développée et validée au cours de cette thèse.

Avec les espaces de finitude, Ehrhard a exhibé une sémantique de la logique linéaire contenant une opération de différentiation. Dans ce cadre, l'interprétation des formules est décomposable en séries de Taylor. Cette étude a engendré des syntaxes différentielles. Cette thèse de sémantique dénotationnelle prolonge ce travail par une exploration de sémantiques vectorielles de la logique linéaire, et contribue à l'étude sémantique et syntaxique de la formule de Taylor. La première partie aborde la sémantique. Nous présentons l'interprétation des constructions de la logique linéaire dans les espaces vectoriels munis d'une topologie linéarisée, les espaces de Lefschetz. Nous définissons une notion intrinsèque d'espaces de finitude, les espaces de Lefschetz finitaires. Nous caractérisons les espaces de Lefschetz réflexifs complets à l'aide de bornologies linéaires. Enfin, nous montrons que la décomposition de Taylor reste valide dans ces espaces. La seconde partie porte sur les syntaxes différentielles. La formule de Taylor syntaxique traduit un terme en une superposition de termes différentiels qui sont autant de possibilités d'exécutions. Comme l'ont montré Ehrhard et Regnier, les termes issus de cette traduction vérifient une relation de cohérence. Nous introduisons une sémantique totale qui capture cette relation. Puis, nous construisons une extension vectorielle du lambda-calcul, le calcul barycentrique, interprété par cette sémantique totale. Enfin, dans le cadre des réseaux différentiels, nous présentons un algorithme non déterministe qui permet de décider si un ensemble fini de réseaux différentiels provient de la traduction d'un réseau de la logique linéaire par la formule de Taylor syntaxique.

L'objet de ce travail a principalement consisté en l'étude et la mise en oeuvre d'une méthode de calcul des gradients des fonctions aérodynamiques par rapport à des paramètres géométriques pour un système aéroélastique soumis à un écoulement lointain stationnaire. Dans un premier temps, une méthodologie de calcul de l'équilibre aéroélastique statique a tout d'abord été développée. Dans ce cadre, le comportement du fluide peut être modélisé par les équations d'Euler ou par les équations de Navier-Stokes moyennées (RANS). Celles-ci sont numériquement résolues par elsA - code de simulation numérique pour la mécanique des fluides développé à l'ONERA. Le comportement de la structure est, quant à lui, prédit par la théorie des poutres et les équations d'Euler-Bernoulli. Le chargement aérodynamique est transmis à la structure par l'intermédiaire de la matrice des coefficients d'influence également appelée matrice de flexibilité. Seuls les efforts de torsion et de flexion sont transmis de manière consistante à la structure, dont seuls les mouvements induits de torsion et de flexion sont calculés sous l'hypothèse des petits déplacements. La déformation résultante sur le maillage du domaine fluide est prédite analytiquement par analogie avec la mécanique du solide. Enfin, le système aéroélastique couplé est résolu selon un processus itératif inspiré de la méthode du point fixe. Dans un deuxième temps, un cadre de calcul, pour le système aéroélastique décrit précédemment, des gradients des fonctions d'intérêt (objectif et contraintes) par rapport à un vecteur de paramètres géométriques de la forme solide a été mis en oeuvre. Les gradients peuvent être calculés par la méthode de l'équation linéarisée discrète ou par la méthode du vecteur adjoint discret. Ces méthodes reposent sur la résolution de systèmes linéaires couplés, effectuée, dans le cadre de cette étude, par un processus itératif doublement retardé. Pour finir, ces développements ont été appliqués au calcul des gradients des coefficients aérodynamiques de traînée et de portance par rapport à un ensemble de paramètres de forme pour trois configurations aérodynamiques de complexité croissante: équations d'Euler résolues sur un maillage multibloc coïncident, équations RANS résolues sur un maillage monobloc, et, finalement, équations RANS résolues sur un maillage multibloc non-coïncident. La validité des résultats a été établie par comparaison aux gradients calculés par différences finies. Une dernière partie du travail a été consacrée à l'évaluation des performances de quatre modèles réduits non physiques dans le cadre d'un processus d'optimisation de forme d'une configuration bidimensionnelle de turbomachine.

Le calcul de gradients de fonctions aérodynamiques par rapport à des paramètres géométriques de la forme solide est une sous-discipline de la simulation numérique en mécanique des fluides. Cette dérivation par rapport à deux variables dépendantes, le champ aérodynamique et le maillage, liés par les équations discrètes de la mécanique des fluides, nécessite l'inversion d'un système linéaire dont la matrice est la matrice jacobienne des équations de la mécanique des fluides discrétisées par rapport au champ aérodynamique (méthode de l'équation linéarisée) ou la transposée de cette matrice jacobienne (méthode adjointe). La précision du calcul de cette matrice jacobienne fait débat lorsque les équations de la mécanique des fluides sont modélisées avec l'approche RANS. L'objectif de cette thèse est de déterminer le degré d'exactitude de la linéarisation d'un flux visqueux discret et des équations discrètes de certains modèles de turbulence, nécessaire à l'obtention de gradients précis de fonctions métiers des concepteurs de turbomachines par rapport à des paramètres géométriques d'une aube. L'écriture de linéarisations approchées du flux visqueux (avec ou sans approche dite de "couche mince") et de deux modèles de turbulence (modèle algébrique de Michel et al. et modèle à deux équations de transport k-e de Launder-Sharma) est détaillée. Pour le modèle de Michel et al., plusieurs approximations pour la linéarisation des équations du modèle ont été testées et comparées. Des résultats (valeurs de gradients de fonctions aérodynamiques, sensibilités d'écoulement pour la méthode linéarisée) sont présentés pour la tuyère de Déléry dite C, l'aile ONERA M6 et deux configurations d'aube isolée de turbine. Des recommandations sont formulées pour le calcul de gradients pour des configurations de machine tournante avec modélisation RANS.

L'imagerie des paramètres physiques du sous-sol à partir d'enregistrements sismiques de surface constitue un problème inverse non linéaire. L'inversion des formes d'onde élastique est une méthode d'imagerie quantitative multiparamètres de diffractants, nécessitant au préalable la connaissance précise d'un macromodèle de vitesse. <br /><br />Le problème direct associé, la propagation des ondes élastiques, est résolu dans le domaine fréquentiel, permettant la prise en compte efficace d'acquisitions multisources et multirécepteurs, par une méthode numérique de différences finies modélisant la propagation de tous les types d'onde (ondes de volume, de surface, diffractées, réfractées ...). Le stencil de différences finies choisi permet simule précisément la surface libre et la propagation des ondes de surface. <br /><br />Cette inversion linéarisée s'appuie sur une méthode de gradient, qui minimise une fonction coût contenant les différences entre données observées et calculées. Deux paramètres sismiques sont reconstruits à partir de sismogrammes verticaux et horizontaux. L'inversion est effectuée des basses fréquences vers les hautes fréquences, introduisant des longueurs d'onde de plus en plus courtes dans les images des paramètres. Ces dernières sont sensibles au choix des approximations physiques effectuées pour calculer le gradient de la fonction coût (approximation de Born ou de Rytov), au dispositif d'acquisition, au préconditionnement des données et au choix des paramètres inversés. <br /><br />Cette méthode est ensuite appliquée à des milieux contenant une surface libre. La surface libre est une interface très contrastée qui donne lieu à des ondes de surface très énergétiques dans les sismogrammes. Les milieux de subsurface sont complexes, les ondes qui s'y propagent subissent des réflexions ou diffractions multiples. Lorsque le milieu contient deux anomalies dont le contraste en vitesse vaut $20 \%$, l'inversion les localise correctement dans un macromodèle connu. Les images sont améliorées et l'amplitude des anomalies est très bien reconstruite lorsque le nombre de fréquences inversées augmente et les données sont sélectionnées des faibles déports vers les grands déports. Avec une anomalie fortement contrastée comme une cavité vide, l'inversion retrouve correctement la position, la forme et la taille de l'objet mais son amplitude est sous-estimée.<br /><br />L'application à des données réelles verticales de subsurface acquises dans un milieu complexe contenant une cavité maçonnée montre que le milieu hétérogène issu de l'inversion ajuste mieux les données et permet de bien reproduire les ondes inversées. Néanmoins, la cavité n'est pas imagée.

Le calcul de quantités dépendants du temps pour des systèmes quantiques est limité le scaling exponentiel des méthodes exactes. Néanmoins, ces quantités présentes un intérêt scientifique important. Un compromis, entre précision et coût, est trouvé par les méthodes quasi-classiques. Dans ces méthodes, la densité thermique exacte est combinée à des trajectoires approximant la dynamique quantique. Durant ma thèse, j’ai développé et appliqué une méthode quasi-classique : PIM (Phase Integration Methode) qui combine des algorithmes MC et MD pour calculer des fonction de corrélation. Le Chapitre 2 décrit les méthodes quasi-classiques ainsi que les approximations qui permettent d’en tirer les fonctions de corrélations dépendants du temps.Le Chapitre 3 illustre comment PIM est adapté au calcul de la densité de Wigner qui est une quantité clé pour les méthodes quasi-classique. À travers le calcul de cette quantité, PIM est capable de capturer des corrélations entre différents degrés de liberté. Dans le Chapitre 4, on montre comment PIM est adapté au calcul de spectres infrarouge. La comparaison des résultats avec d’autres méthodes montre que PIM est une méthode précise pour les systèmes à basse dimensionnalité. Les spectres de OH et CH4 confirment que PIM ne souffrent pas de problèmes intrinsèques comme CMD ou RPMD et peut être appliqué à des systèmes à plus haute dimensionnalité. Le Chapitre 5 présente la méthodologie pour calculer des constantes de vitesse à l’aide de PIM. Les résultats sont bons jusqu’à 300 K mais pas en dessous. Le travail futur se concentrera sur le calcul de la fonction de corrélation de Kubo flux-side pour remédier à ce problème. / Simulation of time-dependent quantities for quantum systems is limited by the exponential scaling of exact methods. However, the calculation of these quantities is key in many problems. A reasonable compromise among accuracy and cost is done by the quasi-classical methods for computing time correlation functions. In these methods, the thermal density is combined with trajectories that approximate quantum dynamics. In my thesis, I develop and apply quasi-classical methods for vibrational spectroscopy. The focus is on the Phase Integration Method. PIM is based on combining MD and MC algorithms to compute appropriate quantities. Chapter 2 is devoted to a general description of the quasi-classical methods. We introduce the different approximations used to compute quantum time correlation functions. Chapter 3 illustrates how PIM is adapted to the calculation of the Wigner density, which is a key quantity in quasi-classical methods. Via this quantity, we show that PIM is able to capture quantum correlation effects among different degrees of freedom. Chapter 4 focuses on the adaptation of PIM for the infrared spectroscopy. Comparison of our results, show that PIM is accurate for low dimensional models. OH and CH4 spectrum confirms that our approach does not suffer from the pathologies such as CMD and RPMD but also that it can treat systems with a larger number of degrees of freedom reliably. Chapter 5 presents the methodology used to calculate rate constants with PIM. The results are in good agreement with the exact reference until 300 K. Future work will focus on using the Kubo flux side correlation function.

Un programme d'inversion des formes d'ondes sismiques a été développé pour estimer les propriétés d'un milieu stratifié élastique dans le domaine temps-distance. Ce programme est capable de prendre en compte des données haute~ fréquences (0-50 Hz) et des milieux finement stratifiés, ce qui en fait une méthode particulièrement adaptée à l'inversion de profils de sismique marine, y compris ceux enregistrés par des sismographes fond de mer. Le problème inverse est résolu par une tchnique de moindres carrés généralisée. Cette technique nécessite de pouvoir résoudre rapidement le problème direct et calculer efficacement la matrice des dérivées de Fréchet qui représente la sensibilité au premier ordre des champs d'ondes par rapport à de faibles modifications des paramètres élastiques du modèle. Le problème direct est résolu en utilisant une technique classique de matrices de réflexion-transmission généralisée combinée avec une intégration en nombres d'ondes discrets. Les dérivées de Fréchet sont obtenues sous forme analytique à l'aide des fonctions de Green du milieu non perturbé. Nous proposons deux approches pour calculer ces dérivées. La première consiste à les calculer directement dans le domaine temps-distance en représentant les interfaces par des arrangements réguliers de points diffractants. La deuxième approche est une nouvelle méthode de calcul basée sur une formulation analytique qui permet d'évaluer avec rapidité et précision les dérivées de Fréchet pour une perturbation des paramètres élastiques d'une couche plane d'extension infinie. La méthode d'inversion proprement dite utilise une technique de gradients conjugués dans le domaine temps-distance, en intégrant toutes les techniques développées précédemment. La méthode est illustrée par des exemples d'inversion de données synthétiques réalistes, mais aussi par des enregistrements réels de type OBS de la campagne Malis en Mer Ligure. La robustesse de la méthode a été étudiée en imposant aux données ou au milieu de propagation certaines erreurs violant les hypothèses de base de l'algorithme. La méthode permet par ailleurs de considérer un grand nombre de géométries source-récepteur et peut s'appliquer en particulier à des enregistrements de profils sismiques verticaux, ainsi qu'à des modèles présentant une forte variabilité des propriétés élastiques suivant la profondeur.

L'objectif de cette thèse est de contribuer à l'optimisation de problèmes dynamiques en présence de retard. Le point de vue qui nous intéressera est celui de Pontryagin qui dans son ouvrage publié en 1962 a donné les conditions nécessaires d'existence de solutions pour ce type de problème. Warga dans son ouvrage publié en 1972 a fait un catalogue des solutions possible, Li et al. ont étudié le cas de contrôle périodique. Notre méthode de démonstration est directement inspirée de la démonstration de P. Michel du cas des systèmes gouvernés par des équations différentielles ordinaires. La principale difficulté pour cette approche est l'utilisation de la résolvante de l'équation différentielle fonctionnelle linéarisée de l'équation différentielle fonctionnelle d'évolution qui gouverne le système. Nous traitons aussi de condition d'Euler-Lagrange dans le cadre d'un problème de calcul variationnel avec retard. / In this thesis, we have attempted to contribute to the optimization of dynamical problems with delay in state space. We are specifically interested in the viewpoint of Pontryagin who outlined in his book published in 1962 the necessary conditions required for solving such problems. In his work published in 1972, Warga catalogued the possible solutions. Li and al. analyzed the case of periodic control. We will treat an optimal control problem governed by a Delay Functional Differential Equation. Our method is close to the one of P. Michel on dynamical system governed by Ordinary Differential Equations. The main problem ariving out in this approach is the use of the resolvent of the Delay Functional Differential Equation. We also consider with Euler-Lagrange condition in the framework of variational problems with delay.