Mes travaux de recherche concernent les méthodes pour la conception optimale en mécanique des fluides numérique. Il s'agit d'un contexte difficile, dans la mesure où l'évaluation de la fonctionnelle coût pour l'optimisation s'appuie sur la résolution des équations aux dérivées partielles régissant l'écoulement, caractérisées par de fortes non- linéarités, anisotropies et éventuellement discontinuités. Plus spécifiquement, mes travaux ont porté sur le développement de méthodes de conception hiérarchiques, pour la modélisation et la paramétrisation, dans le but d'améliorer les propriétés de convergence; sur des méthodes de conception isogéométriques, permettant l'emploi d'une base unique d'ordre élevé pour la représentation géométrique et la simulation; sur des méthodes de conception robuste, pour prendre en compte les incertitudes et les erreurs durant la procédure de conception. Ces développements ont été appliqués à divers problèmes de conception optimale, principalement en aérodynamique compressible, mais aussi en hydrodynamique, pour des modélisations allant des équations d'Euler stationnaires aux équations de Navier-Stokes instationnaires avec fermeture turbulente. Les résultats obtenus ont notamment montré l'intérêt d'introduire un niveau de modélisation supplémentaire : l'utilisation de méta-modèles de type Processus Gaussiens, associés à une recherche des paramètres de conception maximisant l'amélioration probable, conduit à une optimisation globale pour un nombre de simulations significativement réduit (de l'ordre de quelques dizaines). De plus, cette approche peut permettre de prendre en compte les incertitudes et erreurs, en introduisant la notion d'observation bruitée. Le second résultat important concerne l'emploi des méthodes d'analyse isogéométrique, qui offrent un cadre conceptuel intégrant la géométrie et l'analyse. J'ai montré que cette approche permet le développement de schémas d'ordre élevé, de méthodes hiérarchiques pour la modélisation et l'optimisation, et facilite l'analyse de sensibilité et les couplages multidisciplinaires.

En électromagnétisme les structures périodiques suscitent un grand intérêt. Ces structures agissent ainsi comme des filtres fréquentiels et permettent la fabrication de méta-matériaux, composites et artificiels. Elles présentent des propriétés électromagnétiques inédites pour les matériaux naturels telles que des bandes interdites. On a ainsi pu fabriquer de nouveaux dispositifs permettant de guider, de focaliser ou de stopper la propagation. C'est par exemple utile pour éviter le couplage entre différents éléments rayonnants notamment via la caractérisation des ondes de surface qui se propagent à l'interface entre l'air et la structure périodique. Ce travail de thèse s'inscrit dans ce contexte et propose une description de la méthode des éléments finis dédiée à la caractérisation des structures périodiques. La modélisation numérique aboutit à des problèmes de valeurs propres de grandes tailles. Elle implique la résolution de systèmes linéaires composés de matrices creuses. Une méthode est abordée pour résoudre ce type de problème, en optimisant et combinant différents algorithmes. Avant d'aborder les différents aspects de la méthode développée, nous établissons une liste exhaustive de l'ensemble des méthodes qui existent en énonçant leurs avantages et leurs inconvénients. Nous constatons notamment que la méthode des éléments finis permet de traiter un large éventail de structures périodiques en trois dimensions sans limitation sur leur forme géométrique. Nous présentons alors les différentes formulations de cette méthode. Ensuite les aspects algorithmiques de la méthode sont détaillés. Nous montrons notamment qu'une analyse des paramètres de résolution permet de préciser les interprétations physiques des résultats obtenus. Finalement nous présentons les performances de notre outil sur des cas d'applications issus de la littérature et nous abordons la caractérisation des ondes de surface. Pour cela, l'étude d'un réseau d'antennes patchs insérées dans des cavités métalliques est conduite. Notons pour conclure que les études conduites au cours de cette thèse ont abouti à la production d'un code utilisable dans un environnement de calcul initialement présent à l'ONERA.

Cette thèse est divisée en deux parties, le fil directeur étant la contrôlabilité en temps optimal. Dans la première partie, après un rappel du principe du maximum de Pontryagin dans le cas des systèmes de dimension finie, nous mettrons en œuvre ce principe sur le cas d'un intégrateur non-holonome connu sous le nom de système de Brockett pour lequel nous imposons des contraintes sur l'état. La difficulté de cette étude provient du fait que l'on considère un problème de contrôle avec des contraintes sur l'état. Après cet exemple, nous nous intéressons à une extension du principe du maximum de Pontryagin au cas des systèmes de dimension infinie. Plus précisément, l'extension que nous considérons s'applique au cas de systèmes exactement contrôlables en tout temps. Typiquement, ce résultat s'applique à l'équation de Schrödinger avec contrôle interne. Pour de tels systèmes, sous une condition de contrôlabilité approchée, depuis un ensemble de temps non négligeable, nous montrons l'existence d'un contrôle bang-bang. Dans la seconde partie, nous étudions le problème de la nage à bas nombre de Reynolds. Une modélisation physique convenable nous permet de le formaliser comme un problème de contrôle. Nous obtenons alors un résultat de contrôlabilité sur ce problème. Plus précisément, nous montrons que quelque soit la forme du nageur, celui-ci peut se déformer légèrement pour suivre une trajectoire imposée. Nous étudions ensuite le cas d'un nageur à symétrie axiale. Les résultats de la première partie permettent alors la recherche d'un contrôle en temps optimal.

Sur une variété kählérienne compacte connexe de dimension $2m$, $\omega $ étant la forme de Kähler, $\Omega $ une forme volume donnée dans $[\omega ]^m$ et $k$ un entier $1

Cette thèse est scindée en deux parties. Dans un premier temps, nous présentons deux schémas volumes finis pour la discrétisation des problèmes de diffusion-convection-réaction sur les surfaces mouvantes. Le premier schéma présente une extension du schéma volumes finis avec flux à deux points sur les surfaces mouvantes. Le deuxième développe une méthode de type O-méthode. Cette dernière consiste à construire à partir des inconnus au centre des mailles, des solutions linéaires autour des nœuds de maillage qui intègrent la continuité des flux aux interfaces de mailles. La méthode permet aussi la construction des décentrages amont d'ordre 2 et ainsi, offre au procédé l'ordre 2 de convergence sur tout maillage non dégénéré. Ensuite, nous modélisons l'écoulement du couplage filme mince-surfactant (surface active agent) sur les surfaces mouvantes et simulons à l'aide des schémas volumes finis précédemment définis. Ici, l'utilisation du calcul tensoriel et de la méthode d'approximation par la lubrification permettent de réduire les équations de Navier-Stokes caractérisant le mouvement du filme mince en dimension 3 en un système d'équations définies sur la surface courbe mouvante dont l'inconnu est la hauteur du fluide. Le surfactant supposé insoluble est modélisé par une équation de diffusion convection à la l'interface fluid-air. Nous simulons l'ensemble en utilisant une méthode dite de capture d'interface (Interface tracking method) dérivée des volumes finis définis plus haut. Plusieurs exemples illustrent à suffisance l'efficacité et la précision des différentes méthodes.

Cette thèse porte sur les équations d'évolution et s'articule autour de trois parties. Dans la première partie, on se propose de se concentrer sur le critère oscillatoire de certaines équations différentielles. Des résultats classiques sur les fonctions presque-périodiques sont rassemblés dans le premier chapitre. Le deuxième chapitre de cette thèse a pour objectif de prouver l'existence d'une solution presque-périodique de Besicovitch d'une équation différentielle de second ordre sur un espace de Hilbert. L'approche utilisée se base sur un formalisme variationnel. La deuxième partie de cette thèse traite le comportement asymptotique des problèmes de Cauchy dans le cas non autonome. Les semi-groupes et les familles d'évolution étant les outils principaux utilisés dans cette partie, le troisième chapitre introduit des résultats importants de cette théorie, notamment ceux permettant de caractériser la stabilité des semi-groupes et des familles d'évolution périodiques. Dans le quatrième chapitre de cette contribution, on prouve, en utilisant une approche basée sur les semi-groupes, un résultat liant la bornitude de solutions de problèmes de Cauchy périodiques et la stabilité exponentielle uniforme des familles d'évolution issues de ces problèmes. Dans une troisième partie, on focalise l'attention sur quelques résultats sur la dichotomie exponentielle comme une propriété liée au comportement asymptotique des systèmes différentiels. Quelques résultats connus sont, par suite, réunis au cinquième chapitre qui introduit brièvement la notion de dichotomie exponentielle. Dans un dernier chapitre, une caractérisation de la dichotomie exponentielle d'une famille d'évolution en termes de bornitude des solutions de problèmes de Cauchy opératoriels correspondants sera démontrée.

Dans cette thèse on s'intéresse à deux problèmes faisant intervenir des formes normales de champs de vecteurs et des phénomènes exponentiellement petits. Dans le premier chapitre on démontre tout d'abord deux théorèmes de normalisation avec restes exponentiellement petits pour des champs de vecteurs analytiques au voisinage d'un point d'équilibre, dans le cas non autonome périodique. Le premier théorème de normalisation permet de construire une quasi-variété invariante à un exponentiellement petit près, tandis que le deuxième met le champ de vecteur sous la forme normale de Elphick-Tirapegui-Brachet-Coullet-Iooss à un exponentiellement petit près. Dans le deuxième chapitre on travaille près d'un point d'équilibre d'une famille de systèmes hamiltoniens au voisinage d'une résonance 0²iw. On démontre l'existence d'une famille d'orbites périodiques entourant l'équilibre puis l'existence d'orbites homoclines à plusieurs boucles à chacune de ces orbites périodiques, aussi proche de cet équilibre que l'on veut à l'exception de l'équilibre lui-même. La démonstration est basée sur la preuve d'un théorème de forme normale hamiltonien inspiré des formes normales de Elphick-Tirapegui-Brachet-Coullet-Iooss ainsi que sur une normalisation locale hamiltonienne s'appuyant sur un résultat de Moser. On obtient ensuite le résultat grâce à des arguments géométriques liés à la petite dimension et à un théorème KAM qui permet de confiner les boucles. Pour le même problème dans le cadre d'un champ de vecteurs réversible non hamiltonien, l'apparition d'exponentiellement petits lors de la perturbation de l'orbite homocline de la forme normale empêche la démonstration de l'existence d'orbites homoclines à des orbites périodiques de taille exponentiellement petite. Le même phénomène apparait ici mais l'obstacle est contourné grâce à des arguments géométriques spécifiques aux système Hamiltoniens.

Nos travaux concernent l'homogénéisation du transfert de chaleur dans un milieu poreux périodique qui modélise la géométrie d'un cœur de réacteur nucléaire à caloporteur gaz. Cette géométrie est constituée d'un milieu solide traversé par plusieurs longs et minces cylindres parallèles dont le diamètre est du même ordre que la période. La chaleur est transportée par conduction dans la partie solide du domaine et par conduction, convection et rayonnement dans la partie fluide (les cylindres). Le rayonnement est modélisé par une condition non-locale sur les parois des cylindres. C'est une analyse stationnaire qui correspond à un fonctionnement nominal du cœur, et aussi non-stationnaire qui correspond à un arrêt nor- mal du cœur. Pour obtenir le problème homogénéisé nous utilisons d'abord une analyse formelle par développement asymptotique à deux échelles. La justification mathématique de nos résultats est basée sur la méthode de convergence à deux échelles. Une caractéristique de ce travail en dimension 3 est qu'il combine l'analyse asymptotique par homogénéisation avec une analyse asymptotique par réduction de la dimension de l'espace 3D en 2D pour remédier à la non-périodicité de la condition de rayonnement suivant la direction axiale des cylindres. Une deuxième caractéristique de ce travail est l'étude de ce transfert de chaleur lorsqu'il contient une source thermique oscillante au niveau microscopique et un échange thermique entre les parties fluide et solide du cœur, dans un tel contexte, notre analyse numérique montre une contribution non-négligeable du correcteur dit d'ordre 2 qui nous aide à reproduire les gradients qui apparaissent entre la zone de la source thermique et la partie fluide (les cylindres).

Dans ce travail, nous développons un cadre de calcul dédié à la simulation d'écoulements à plusieurs fluides. Nous présentons des validations et vérifications de ces méthodes sur des problèmes de capture d'interfaces et de simulations de bulles visqueuses. Nous montrons ensuite que ce cadre de calcul est adapté à la simulation d'objet rigides en écoulement. Puis, nous étendons ces méthodes à la simulation d'objets déformables simulant le comportement des globules rouges : les vésicules. Nous validons aussi ces simulations. Enfin nous appliquons les précédents modèles à des problèmes ouverts de microfluidique tels que la séparation d'une suspension dans une bifurcation microfluidique et la rhéologie en milieu confiné.

Cette thèse est consacrée à l'étude de la dynamique des gaz à faible nombre de Mach. Le modèle étudié provient des équations de Navier-Stokes complètes lorsque le nombre de Mach tend vers zéro. On cherche à montrer que le problème de Cauchy correspondant est bien posé. Les cas visqueux et non visqueux sont tous deux considérés. Les coefficients physiques peuvent dépendre de la densité (ou de la température) inconnue. En articulier, nous prenons en compte les effets de onductivité thermique et on autorise de grandes variations d'entropie. Rappelons qu'en absence de diffusion thermique, la limite à faible nombre de Mach implique la condition d'incompressibilité. Dans le cadre étudié ici, en introduisant un nouveau champ de vitesses à divergence nulle, le système devient un couplage non linéaire entre une équation quasi-parabolique pour la densité et un système de type Navier-Stokes (ou Euler) pour la vitesse et la pression. \\\\ Pour le cas avec viscosité, on établit le résultat classique, à savoir qu'il existe une solution forte existant localement (resp. globalement) en temps pour des données initiales grandes (resp. petites). On considère ici le problème de Cauchy avec données initiales dans des espaces de Besov critiques. Lorsque les coefficients physiques du système vérifient une relation spéciale, le système se simplifie considérablement, et on peut alors établir qu'il existe des solutions faibles globales en temps à énergie finie. Par un argument d'unicité fort-faible, on en déduit que les solutions fortes à énergie finie existent pour tous les temps positifs en dimension deux. \\\\ Pour le cas sans viscosité, on montre d'abord le caractère bien posé dans des espaces de Besov limites, qui s'injectent dans l'espace des fonctions lipschitziennes. Des critères de prolongement et des estimations du temps de vie sont établis. Si l'on suppose la donnée initiale à énergie finie dans l'espace de Besov limite à exposant de Lebesgue infini, on a également un résultat d'existence locale. En dimension deux, le temps de vie tend vers l'infini quand la densité tend vers une constante positive. \\\\ Des estimations de produits et de commutateurs, ainsi que des estimations a priori pour les équations paraboliques et pour le système de Stokes (ou d'Euler) à coefficients variables, se trouvent dans l'annexe.Ces estimations reposent sur la théorie de Littlewood-Paley et le calcul paradifférentiel.