Cette thèse s'inscrit dans le cadre d'étude de la simulation de la propagation du son en écoulement. L'objectif de ces travaux est l'obtention de modèles approchés permettant une prise en compte aisée des zones de fortes variations de l'écoulement porteur (couche limite de paroi, couche de mélange...). Le modèle mathématique retenu pour l'étude est celui des équations de Galbrun. La première partie est consacrée à la propagation acoustique dans un tuyau mince bidimensionnel. Une analyse asymptotique qui s'apparente à une analyse basse fréquence est menée pour obtenir un problème approché original, faisant intervenir un terme intégral non local vis à vis de la coordonnée transverse. Du fait de son originalité, l'analyse de stabilité est complexe et nécessite une étude ad hoc. Cette approche nouvelle permet de retrouver des résultats sur la stabilité des écoulements incompressible, mais aussi d'en établir de nouveaux. Nous proposons ensuite une méthode de résolution numérique basée sur une expression quasi-explicite de la solution. La question de la prise en compte des couches limites de paroi fait l'objet de la deuxième partie. Nous considérons toujours un problème bidimensionnel à paroi plane. Les cas d'une paroi parfaitement rigide et d'une paroi sur laquelle on impose une condition d'impédance sont traités. Dans les deux cas nous remplaçons la couche limite par une condition aux limites approchée, au moyen d'une analyse asymptotique. Ces conditions font intervenir la résolution du problème limite du tube et l'analyse de stabilité repose sur les résultats de la première partie. Nous explorons ensuite les propriétés physiques et mathématiques de ces problèmes approchés.

Ce travail de thèse expérimental étudie les instabilités et la transition de la couche limite produite par un disque en rotation. Pour l'écoulement naturel (c.-à-d. sans forçage extérieur), les mesures des profils de vitesse moyenne, de spectres en fréquence et de moyennes de phase des séries temporelles de vitesse ont permis de distinguer différents régimes en fonction de la distance adimensionnelle R à l'axe du disque. Pour les faibles valeurs de R, les profils de vitesse moyenne suivent la solution de von Kármán. Pour des valeurs plus importantes de R, des écarts à cette solution analytique sont observés et augmentent avec R. Ces écarts sont dus à la croissance spatiale de modes instables de la couche limite (vortex .cross-flow.), et la mesure du taux de croissance spatiale de ces modes correspond bien aux prédictions théoriques de l'analyse de stabilité linéaire. Dans cet écoulement, la transition se produit vers R ≈ 530 et la turbulence pleinement développée s'installe vers R ≈ 600. Les profils dans la région pleinement turbulente suivent la loi logarithmique des couches limites turbulentes et les spectres de vitesse présentent une loi en puissance de type Kolmogorov. Pour étudier la réponse au forçage, un dispositif expérimental a été mis au point qui permet d'exciter des perturbations stationnaires (dans le référentiel du laboratoire) ou en rotation à une fréquence qui peut être réglée indépendamment de la fréquence de rotation du disque. La réponse de l'écoulement à ces deux types de forçage et avec deux formes différentes pour l'élément de forçage a été étudiée. Un forçage stationnaire produit un sillage qui décroît avec la distance à l'élément de forçage, en accord avec la théorie. Le forçage avec des éléments en rotation peut produire un paquet d'ondes amplifié qui, bien que non linéaire, suit des trajectoires proches de celles prédites par la théorie linéaire.

Dans cette thèse, nous étudions deux problèmes différentiels dépendant d'un paramètre ε et étu- dions l'asymptotique des solutions lorsque ce paramètre tend vers 0. Le premier problème est lié à l'équation de Helmholtz haute-fréquence. On construit un potentiel non captif ne satisfaisant pas l'hypothèse de refocalisation des rayons introduite par F. Castella. On montre que l'ensemble des tra jectoires hamiltoniennes (associées au potentiel construit) issues de l'origine et qui reviennent en 0 forme une sous-variété de dimension d − 1, où d est la dimension de l'espace. On montre alors que la solution de l'équation de Helmholtz converge vers une perturbation de la solution de Helmholtz avec condition de radiation à l'infini et coefficients figés en 0. Dans un second temps, nous étudions une équation de Navier-Stokes forcée par une source po- larisée fortement oscillante. On exhibe une famille de solutions exactes. On étudie alors la stabilité de cette famille lorsqu'on la perturbe à l'instant initial. On construit une solution approchée du pro- blème à l'aide d'une couche limite à l'instant initial (t=0). Ce développement montre en particulier que des interactions d'ondes, se propageant à des échelles différentes, peuvent se traduire au niveau macroscopique par une augmentation de la viscosité. Enfin, on justifie la convergence de la solution approchée vers la solution exacte à l'aide de méthodes d'énergie.

L'objectif principal de cette thèse est de réaliser des études analytiques et numériques pour des problèmes de capillarité, de transfert de chaleur et de stabilité du procédé Bridgman démouillage. Pour le calcul de la forme du ménisque, sa surface sera donnée par l'équation de Young-Laplace décrivant l'équilibre sous la pression. Cette équation sera transformée en un système non linéaire d'équations différentielles. A partir d'études qualitatives et quantitatives de la solution, la dépendance de la forme du ménisque (convexe, concave, convexe-concave), de la différence de pression et d'autres paramètres du procédé, sera déterminée. Pour étudier la stabilité dynamique du système, l'épaisseur de l'espacement cristal-creuset et la position de l'interface liquide-solide sont des variables du problème et donc deux équations seront nécessaires, précisément, l'équation de Young-Laplace et le bilan thermique à l'interface liquide-solide. Par conséquence, ce travail est organisé comme suit: Des contributions récentes à la modélisation de certains problèmes de capillarité sont présentées dans le deuxième chapitre, à commencer par la formulation mathématique du problème capillaire régie par l'équation de Young-Laplace. Des études analytiques et numériques pour l'équation du ménisque sont élaborées pour le démouillage en microgravité et sur terre. Le troisième chapitre traite des contributions à la modélisation des problèmes de transfert de chaleur. Ainsi, les études analytiques et numériques pour l'équation non stationnaire de transfert de chaleur à une dimension sont effectuées afin de trouver des expressions analytiques de la distribution de la température et des gradients de température dans le liquide et dans le solide. L'équation de déplacement de l'interface liquide-solide est également obtenue du bilan énergétique à l'interface. Après quoi, l'effet de l'espacement cristal-creuset sur la courbure de l'interface liquide-solide est étudié pour un ensemble de paramètres représentatifs non-dimensionnels de la croissance de cristaux semi-conducteurs classiques. Une expression analytique pour la déflexion de l'interface, basée sur la théorie du flux de chaleur est rapportée. Afin de vérifier l'exactitude de la formule obtenue analytiquement et d'identifier ses limites de validité, l'équation de transfert de chaleur est résolue numériquement dans une symétrie axiale en 2D, pour un cas stationnaire et en utilisant le code d'éléments finis COMSOL Multiphysics 3.3. Le dernier chapitre est entièrement consacré à l'analyse de la stabilité. Tout d'abord, différents concepts de stabilité de Lyapunov qui peuvent survenir dans la croissance des cristaux: classique, uniforme, asymptotique et exponentielle d'un état d'équilibre; stabilité partielle de Lyapunov d'un état d'équilibre, et les mêmes types de stabilités de Lyapunov pour la solution temporelle sont présentés. Dans ce qui suit, après l'introduction de la notion de stabilité pratique sur une période de temps limitée, des études analytiques et numériques de la stabilité pratique sur une période de temps limitée du système non linéaire, des équations différentielles décrivant le déplacement d'interface liquide-solide et l'évolution de l'espacement cristal-creuset, pour des cristaux élaborés par le procédé Bridgman démouillage sous conditions terrestres sont développés. Enfin, les conclusions générales et perspectives de ce travail sont exposées.

L'objet de cette thèse est l'étude de l'espace de modules des G2-fibrés principaux sur une courbe complexe projective connexe lisse, où G2 désigne le groupe de Lie exceptionnel de plus petit rang. Le groupe G2 est caractérisé via trois approches différentes, la première étant celle où G2 est défini comme le groupe des automorphismes de l'algèbre complexe des octaves de Cayley. Les différentes réductions et extensions que peut admettre un G2-fibré principal sont étudiées ainsi que la relation entre la stabilité d'un G2-fibré principal et celle du fibré vectoriel qui lui est associé. L'espace de modules des G2-fibrés principaux semi-stables est analysé. Nous obtenons notamment une caractérisation de son lieu lisse, une décomposition explicite de son lieu singulier en trois composantes connexes et une analyse de l'espace de Verlinde de niveau 1 pour le groupe G2.

Les contributions que nous avons choisi de présenter dans ce mémoire ont trait en général à la commande, à la stabilité et à l'observation des systèmes dynamiques hybrides. Ce mémoire retrace le cheminement de mes recherches sur les dix dernières années au CRAN au sein du projet systèmes dynamiques hybrides et complexes mais également aux travers de projets et collaborations locales (CPER : projet SC2), nationales (GdR SDH, AS 192, ANR ArHyCo) et européennes (Hycon). L'extension du principe du minimum et l'établissement de conditions nécessaires pour la résolution d'un problème de commande optimale pour les SDH ont permis de soulever et de mettre en évidence des difficultés liées à la résolution de ce type de problème. Deux freins principaux à la mise en oeuvre d'algorithmes généraux efficaces ont été identifiés et proviennent d'une part, de l'explosion combinatoire engendrée par la dynamique discrète couplée aux dynamiques continues et d'autre part, de l'existence de trajectoires singulières conduisant nécessairement à des solutions sous optimales. Une discussion sur les méthodes de résolution directes et indirectes montre que même dans des situations simplifiées, l'existence de trajectoires singulières pose un réel problème et les algorithmes utilisant le principe du minimum sont inadaptés sans traitement particulier pour gérer ces situations. Nous montrons sur un exemple qu'il peut être préférable d'utiliser des méthodes directes sous une formulation de tirs multiples. Pour la classe des systèmes affines commutés qui représente une classe de systèmes technologiques très répandus, les points de fonctionnement rendent particulièrement évident le rôle joué par les trajectoires singulières. On constate qu'elles correspondent aux solutions optimales au sens de Fillipov du système commuté. Nous avons proposé une méthode permettant d'effectuer la synthèse de ces trajectoires optimales pour des critères de type temps optimal ou quadratique en temps infini applicable sur des systèmes de dimension faible et en général inférieure ou égale à 3. Une seconde méthode de synthèse à base de commande prédictive utilise les degrés de liberté sur la commande autour du point de fonctionnement pour poursuivre un cycle particulier défini par un critère. L' algorithme rapide d'optimisation des solutions tient compte des discontinuités issues des instants de commutations et utilise les fonctions de sensibilité vis à vis de ces instants. La généralisation aux SDH du calcul des fonctions de sensibilité lorsque le système hybride est le siège de discontinuités de champs et de saut sur l'état, est un résultat utile qui permet d'utiliser des outils nécessitant une différentiation. Le calcul rapide des cycles limites et l'étude de leur stabilité est une première application possible, l'optimisation est très clairement la cible principale. L'étude et l'analyse du comportement asymptotique des systèmes linéaires commutés via la formulation d'un principe d'invariance de Lasalle nous amène à un constat similaire à celui fait sur la commande optimale : Le comportement asymptotique est identique à celui que l'on obtient en considérant les solutions de Fillipov du système commuté. Nous avons montré en considérant plusieurs hypothèses sur les lois de commutation que la séparation est opérée si on restreint la loi de commutation en considérant par exemple des hypothèses fortes de temps d'activation minimum sur les modes. Ces résultats d'analyse sont primordiaux dans la compréhension de la complexité de la dynamique issue de la commutation. La synthèse d'observateurs commutés pour la classe des systèmes linéaires commutés que nous avons proposée, et l'identification algébrique des lois ne permettant pas d'observer le système ont été déduites du principe d'invariance mentionné à l'item précédent. La mise en évidence de la décroissance de la vitesse de convergence de l'erreur d'observation en fonction de la vitesse de commutation est une conséquence de cette caractérisation algébrique. Les travaux réalisés ont permis de soulever de très nombreux problèmes spécifiques à la classe des systèmes commutés (Robustesse vis à vis d'incertitudes paramétriques, extension au cas discret et échantillonné, synthèse de lois préservant l'observabilité, vitesse de convergence de l'observateur, ...) et de dégager des pistes de recherche concernant leur résolution.

La compréhension des mécanismes de transition La compréhension des mécanismes de transition d'écoulements décollés laminaires trouve un intérêt particulier de par sa présence au sein de nombreuses applications industrielles. Notamment, la prédiction précise de la première bifurcation associée à un tel écoulement et la construction de modèles réduits, pourrait autoriser des développements d'outils de contrôle efficaces, permettant une augmentation des performances aérodynamiques des profils considérés. Dans ce cadre, nous proposons à travers ce mémoire une analyse de stabilité originale, prenant en compte le caractère fortement non parallèle de l'écoulement. En particulier, une étude de stabilité linéaire "globale", sur une configuration académique d'écoulement décollé de plaque plane, révèle l'influence des modes globaux bidimensionnels et tridimensionnels, dans la dynamique en espace et en temps de la perturbation, en fonction du type de forçage appliqué. En outre, une telle étude identifie des mécanismes non détectés par les approches classiques, comme une très large amplification transitoire 2D de la perturbation le long de la couche de mélange, et un mécanisme centrifuge pouvant conduire à une tridimensionnalisation du décollement. Enfin, la décomposition de la perturbation dans la base de modes globaux s'avère une réduction de modèle pertinente, dans l'objectif d'un contrôle efficace.

La plupart des écoulements en dynamique des fluides génèrent ou rencontrent des phénomènes instationnaires. Dans le domaine de l'aérospatial, au cours de la mise au point de certains moteurs de fusées d'altitude, des instationnarités basses fréquences ont été observées. Ces phénomènes, liés à l'interaction entre l'onde de choc et la couche limite turbulente qui se développe sur les parois, peuvent engendrer des déformations de la tuyère ou des mouvements inopportuns par rapport à son système d'attache. De nombreux travaux ont ´et´e entrepris afin d'en déterminer l'origine. L'´etude de la stabilité entre dans ce cadre. Dans ce contexte, le propos de cette thèse est de développer un code de stabilité globale adapté aux écoulements en régime compressible et en géométrie curviligne afin d'étudier la stabilité de l'interaction au sein de la tuyère. La résolution du problème est assuré par un schéma aux différences finies de type DRP. Le code est validé à travers plusieurs cas-tests incompressibles, compressibles et en maillage curviligne. Ces configurations permettent également de souligner les avantages du schéma de discrétisation retenu par rapport `a d'autres solutions plus classiques utilisées généralement dans ce type de problèmes.

Face aux problèmes vibratoires en usinage, deux approches ont été développées. L'une consiste à prédire les vibrations grâce à des modélisations de plus en plus élaborées des opérations d'usinage. La deuxième est une approche plus expérimentale, d'instrumentation des machines dans le but de surveiller et de comprendre le comportement de l'outil pendant l'usinage. Dans le cadre de cette deuxième approche, de nombreuses technologies de capteurs sont utilisées. Le but de ces travaux est d'étudier le comportement de l'outil grâce à des capteurs de déplacement sans contact qui mesurent directement les vibrations de l'outil pendant l'opération d'usinage. Ces capteurs d'abords utilisés en fraisage puis en tournage permettent l'étude des différents aspects de l'interaction entre l'outil et la pièce. Le premier point étudié est l'investigation du phénomène de la coupe, avec la possibilité de tracer en trois dimensions ou dans le plan la trajectoire de l'outil. Le deuxième point est l'étude de la stabilité ou de l'instabilité de l'usinage. Par la suite, les mesures de déplacement sont utilisées pour déterminer les efforts de coupe, que ce soit en fraisage ou en tournage. Les efforts appliqués sur l'outil vont provoquer sa déformation qui sera mesurée par les capteurs sans contact. Dans le cas du fraisage, une méthode de calcul spécifique à partir de la Fonction de Réponse en Fréquence (FRF) de l'outil a été mise en œuvre pour calculer ces efforts de coupe pour des vitesses de rotation importantes. Le dernier aspect traité est la reconstruction de l'état de surface de la pièce usinée. Les déplacements mesurés de l'outil sont superposés à la trajectoire théorique de l'outil pour déterminer l'aspect et le profil de la surface réellement usinée lors de l'opération.

Ce travail s'intéresse à la résolution du système de Maxwell dans le domaine temporel sur des maillages héxaédriques orthogonaux localement raffinés de manière conforme ou non-conforme. Une méthode de Galerkin discontinu, reposant sur une approximation centrée pour le calcul des intégrales de surface et un schéma saute-mouton d'ordre 2 pour l'intégration temporelle est présentée. On définit ainsi une classe de schémas non-diffusifs : un équivalent discret de l'énergie électromagnétique est conservée.