Dans la première partie, nous avons étudié l'effet du bruit sur la vitesse d'un front décrit par une équation de type Fisher-Kolmogorov. Ces équations interviennent souvent comme la limite d'un modèle aléatoire faisant intervenir N particules quand N devient grand. Elles ont beaucoup de solutions, mais c'est la vitesse marginalement stable v^* qui est sélectionnée pour une condition initiale localisée. Nous avons montré que si l'on prend en compte l'aspect discret du modèle microscopique en ajoutant un cut-off d'ordre 1/N dans la queue du front, alors, quelles que soient les conditions initiales, la vitesse de propagation v_N est proche de v^* et la différence v^* - v_N est d'ordre (log N)^(-2). Ces résultats peuvent s'appliquer au modèle aléatoire: grâce à des simulations faisant intervenir jusqu'à 10^14 particules, nous avons observé une correction de la vitesse compatible avec celle obtenue dans le modèle avec \textit{cut-off}. La méthode que nous avons employée permet également de retrouver les résultats de Bramson sur l'influence des conditions initiales sur la vitesse d'un front. La seconde partie est consacrée aux polymères dirigés dans un milieu aléatoire de largeur finie. La méthode des répliques permet de ramener le calcul des fluctuations de l'énergie libre d'un tel polymère à un problème de mécanique quantique avec n particules en interaction. Ce modèle peut être résolu grâce à l'Ansatz de Bethe, mais il faut extrapoler les solutions à des n non-entiers pour faire le lien avec l'énergie libre d'un polymère. Nous avons présenté une méthode qui nous a permis de calculer exactement les premiers cumulants de cette énergie libre. De plus, pour une dimension transversale périodique, on peut calculer tous ces cumulants dans la limite où la largeur du système devient grande et déterminer ainsi la distribution de l'énergie libre. Cette distribution est la même que celle obtenue dans le modèle ASEP et semble donc être une propriété universelle de l'équation KPZ.

Etude du signal de corrélations temporelles entre photons émis dans les bandes latérales du triplet de fluorescence. Ce signal constitue une analyse mixte de la lumière de fluorescence puisqu'il concerne a la fois ses caractéristiques spectrales et temporelles. Etude théorique de cet effet: le calcul perturbatif du champ diffuse par l'atome et la méthode utilisant les équations de bloch ont permis d'interpréter l'effet mis en évidence expérimentalement. Une autre approche théorique présentée est la méthode de l'atome habille.

Cette thèse est consacrée à l'étude mathématique de certains aspects de l'équation de Vlasov-Poisson, qui constitue un modèle cinétique classique en physique des plasmas. Dans un premier temps, nous nous intéressons à la justification rigoureuse d'approximations de l'équation de Vlasov-Poisson avec un champ magnétique extérieur intense, qui sont couramment utilisées, notamment lors des simulations numériques. Le but est de décrire certains régimes d'intérêt par des modèles asymptotiques, obtenus en faisant tendre un petit paramètre vers 0 (modélisant la physique du problème considéré) dans les équations originelles. Nous étudions pour commencer la limite quasineutre, c'est-à-dire la limite quand la longueur de Debye tend vers 0, pour l'équation de Vlasov-Poisson avec des électrons suivant une loi de Maxwell-Boltzmann. Dans la limite des plasmas froids, à l'aide de la méthode de l'entropie relative et de techniques de filtrage, nous montrons la convergence vers des équations hydrodynamiques compressibles telles que l'équation d'Euler isotherme. Nous nous intéressons ensuite à l'approximation "rayon de Larmor fini" en trois dimensions, qui permet de décrire le comportement turbulent d'un plasma soumis à un champ magnétique intense. Pour cette étude, qui peut en fait être interprétée comme une limite quasineutre anisotrope, nous montrons des résultats très différents selon la dynamique décrite. En effet, dans le cas de la dynamique avec des électrons sans masse, nous exhibons un effet stabilisant qui permet d'obtenir le même résultat que pour le système bidimensionnel, alors que pour la dynamique avec des ions lourds, nous mettons en évidence les conséquences d'instabilités de type multi-fluides. Dans un second temps, nous nous consacrons à l'étude mathématique du confinement d'un plasma de tokamak. Nous commençons par proposer un modèle hydrodynamique simplifié à deux températures et étudions la stabilité au sens de Lyapunov de deux états stationnaires permettant de modéliser l'équilibre du plasma. Nos résultats sont conformes à l'heuristique physique et mettent de surcroit en évidence qu'un fort gradient de température favorise la stabilité : cela pourrait fournir une explication aux modes de haut confinement (H-modes) dans les tokamaks. Pour finir, nous attaquons ce problème du point de vue de la théorie du contrôle et prouvons des résultats pour l'équation de Vlasov-Poisson en présence de champs extérieurs (typiquement un champ magnétique).

Dans cette thèse, on s'intéresse à la contrôlabilité de différentes équations aux dérivées partielles. La première partie est consacrée à la méthode de Lebeau-Robbiano pour le contrôle des équations paraboliques linéaires. On étend tout d'abord cette méthode à des opérateurs elliptiques non-autoadjoints, montrant une inégalité spectrale ainsi que la contrôlabilité de l'équation parabolique associée. On prouve ensuite ces deux propriétés pour un modèle de transmission à travers une interface, pour lequel la condition de transmission implique une diffusion tangentielle. La preuve repose sur une inégalité de Carleman, uniforme par rapport au petit paramètre représentant l'épaisseur de l'interface. Dans la deuxième partie, on analyse les propriétés de certains systèmes d'équations aux dérivées partielles linéaires couplées par des termes d'ordre zéro. Après avoir étudié la stabilisation de deux équations d'ondes, dont une seulement est amortie, on montre la contrôlabilité en temps grand d'un système similaire au moyen d'un seul contrôle, sous des conditions géométriques optimales sur les zones de contrôle et de couplage. Par des méthodes d'analyse microlocale, on obtient de plus la contrôlabilité de systèmes d'ondes en cascade, ainsi que l'expression exacte du temps minimal de contrôle. On déduit de ces résultats la contrôlabilité des systèmes paraboliques associés, dans des situations où les zones de contrôle et de couplage sont disjointes. Enfin, dans la troisième partie, on étudie la contrôlabilité uniforme de perturbations visqueuses de lois de conservation scalaires, dans la limite de viscosité évanescente. On montre la contrôlabilité exacte globale aux états constants au moyen de contrôles uniformément bornés lorsque la viscosité tend vers zéro.

La tomographie des temps de première arrivée cherche à estimer un modèle de vitesse de propagation des ondes sismiques à partir des temps de première arrivée pointés sur les sismogrammes. Le modèle de vitesse obtenu peut alors permettre une interprétation structurale du milieu ou bien servir de modèle initial pour d'autres traitements de l'imagerie sismique. Les domaines d'application de cette méthode s'étendent, à des échelles différentes, de la géotechnique à la sismologie en passant par la géophysique pétrolière. Le savoir-faire du géophysicien joue un rôle important dans la difficile résolution du problème tomographique non-linéaire et mal posé. De nombreuses recherches ont entrepris de faciliter et d'améliorer cette résolution par des approches mathématique ou physique. Dans le cadre de ce travail, nous souhaitons développer une approche pragmatique, c'est-à-dire que nous considérons que le problème tomographique doit être résolu par un algorithme interactif dont les paramètres de réglage sont clairement définis. L'aspect interactif de l'algorithme facilite l'acquisition du savoir-faire tomographique car il permet de réaliser, dans un temps raisonnable, de nombreuses simulations pour des paramétrisations différentes. Le but poursuivi dans cette thèse est de définir, pour le cas spécifique de la tomographie des temps de première arrivée, un algorithme qui réponde au mieux à ces critères. Les algorithmes de tomographie des temps de première arrivée classiquement mis en oeuvre aujourd'hui ne répondent pas à nos critères d'une approche pragmatique. En effet, leur implémentation ne permet pas d'exploiter l'architecture parallèle des supercalculateurs actuels pour réduire les temps de calcul. De plus, leur mise en oeuvre nécessite une paramétrisation rendue complexe du fait de la résolution du système linéaire tomographique. Toutes ces limitations pratiques sont liées à la formulation même de l'algorithme à partir de la méthode de Gauss- Newton. Cette thèse repose sur l'idée de formuler la résolution du problème tomographique à partir de la méthode de plus grande descente pour s'affranchir de ces limitations. L'étape clé de cette formulation réside dans le calcul du gradient de la fonction coût par rapport aux paramètres du modèle. Nous utilisons la méthode de l'état adjoint et une méthode définie à partir d'un tracé de rais a posteriori pour calculer ce gradient. Ces deux méthodes se distinguent par leur formulation, respectivement non-linéaire et linéarisée, et par leur mise en oeuvre pratique. Nous définissons ensuite clairement la paramétrisation du nouvel algorithme de tomographie et validons sur un supercalculateur ses propriétés pratiques : une parallélisation directe et efficace, une occupation mémoire indépendante du nombre de données observées et une mise en œuvre simple. Finalement, nous présentons des résultats de tomographie pour des acquisitions de type sismique réfraction, 2-D et 3-D, synthétiques et réelles, marines et terrestres, qui valident le bon comportement de l'algorithme, en termes de résultats obtenus et de stabilité. La réalisation d'un grand nombre de simulations a été rendue possible par la rapidité d'exécution de l'algorithme, de l'ordre de quelques minutes en 2-D.

Cette étude est consacrée à l'optimisation géométrique du contrôle actif dans les gaines de ventilation. La première partie de ce travail concerne l'étude du contrôle hybride (actif+passif) dans les gaines de ventilation et de ses avantages suivant la géométrie de la gaine et la présence ou non de revêtements absorbants. Les effets du contrôle passif du bruit (dissipatif, réflexif et diffusif) dans une gaine de ventilation améliorent l'efficacité du contrôle actif. Un contrôle hybride (actif + passif) présente donc des avantages, même aux basses fréquences. La seconde partie de ce mémoire s'intéresse au placement des microphones d'erreur et des sources secondaires. Pour une excitation harmonique, une méthode de calcul par programmation linéaire et entière détermine un nombre suffisant de microphones d'erreur ainsi qu'une position optimale des microphones d'erreur et des sources secondaires. Cette méthode de placement est appliquée pour la réduction du bruit dans les gaines de ventilation.

L'étude consiste à décrire les phénomènes thermiques et aérauliques dans le bâtiment à l'aide d'un outil simplifié : la méthode zonale. Il s'agit d'une méthode tridimensionnelle basée sur le partitionnement en un petit nombre de sous-volumes, intermédiaire entre les modèles à un noeud et les maillages fins. On écrit des bilans de masse et d'énergie dans chaque sous-volume tandis que les échanges dans les interfaces sont déterminés par des lois reliant les débits aux différences de pression. L'aspect modulaire de la méthode facilite son implémentation dans un environnement orienté objet et le logiciel SPARK, adapté à la résolution de gros systèmes d'équations non-linéaires est utilisé à cet effet. Les résultats sont validés par rapport à différentes références expérimentales et numériques. Une étude paramétrique détermine les coefficients empiriques judicieux ainsi que les caractéristiques d'un maillage optimal. Un autre atout d'un environnement objet réside dans les possibilités de couplage. On traitera successivement l'exemple d'un modèle de description du confort, celui des transferts par conduction en tenant compte des effets tridimensionnels ainsi qu'un modèle de transferts de masse. La simulation des effets d'une source de chaleur donne des résultats conformes aux constatations expérimentales dans l'ensemble du volume. Enfin, le cas de la convection mixte est traité en prenant en compte la conservation de l'énergie cinétique dans l'écoulement ; les résultats correspondent à ceux obtenus avec un modèle de champ.

Le travail présenté dans cette thèse porte sur l'étude numérique d'écoulements réactifs en régime transsonique. Le modèle mathématique pour résoudre ces écoulements est constitué des équations d'Euler multi-espèces écrites pour un mélange de gaz parfaits auxquelles nous ajoutons les termes modélisant les effets diffusifs et les effets dus à la combustion. Pour résoudre numériquement le système d'équations ainsi défini, nous utilisons une formulation mixte éléments finis - volumes finis basée sur la méthodologie MUSCL et utilisant des fonctions de flux numériques décentrés. Pour le traitement des conditions aux limites, on utilise des flux numériques adaptés aux écoulements multi-espèces qui traitent les frontières où l'écoulement est subsonique ou supersonique. A l'aide de ce schéma on étudie différents problèmes monodimensionnels de détonation stable et instable ainsi que des problèmes bidimensionnels d'interaction réactive de jets et de flammes de diffusion. La difficulté des cas tests étudiés avec des rapports de pression parfois supérieurs à 30 prouve la robustesse de la méthode. De plus l'utilisation de la méthode MUSCL associée à des maillages fins obtenus par raffinement statique ou dynamique conduit à des solutions numériques précises et sans oscillations. Pour s'affranchir des difficultés liées à la disparité des échelles temporelles qui apparaissent dans ces écoulements, on propose différents schémas explicites et implicites par décomposition des opérateurs qui permettent (tout en conservant une solution instationnaire de bonne qualité) des gains supérieurs à 5 par rapport à l'explicite.

CETTE THESE RASSEMBLE DIVERSES CONTRIBUTIONS MATHEMATIQUES ET NUMERIQUES A LA CHIMIE QUANTIQUE. LE CHAPITRE 1 EST CONSACRE A UNE PRESENTATION DE L'ESPRIT ET DES MODELES DE LA CHIMIE QUANTIQUE. LE CHAPITRE 2 TRAITE DE LA CONVERGENCE D'ALGORITHMES POUR LA RESOLUTION DES EQUATIONS DE HARTREE-FOCK. LES CHAPITRES SUIVANTS PORTENT SUR DES PROBLEMES SPECIFIQUES AUX SYSTEMES MOLECULAIRES IN SITU, C'EST-A-DIRE EN INTERACTION AVEC UN ENVIRONNEMENT EXTERIEUR. UNE PREMIERE APPROCHE POUR SIMULER LES EFFETS D'ENVIRONNEMENT CONSISTE A TRAITER L'INTERACTION ENTRE LE SYSTEME MOLECULAIRE ET LE MILIEU EXTERIEUR COMME UNE PERTURBATION. AU CHAPITRE 3, ON ETEND LA THEORIE DES PERTURBATIONS DES OPERATEURS LINEAIRES AU CADRE NON LINEAIRE DU MODELE DE HARTREE-FOCK. L'INTERACTION D'UN SYSTEME MOLECULAIRE AVEC UN ENVIRONNEMENT EST SOUVENT UN PROCESSUS DYNAMIQUE. C'EST LE CAS BIEN EVIDEMMENT DES QU'ON ETUDIE UNE REACTION CHIMIQUE. LE CHAPITRE 4 CONSISTE EN L'ANALYSE MATHEMATIQUE D'UNE DES APPROXIMATIONS DE L'EQUATION DE SCHODINGER DEPENDANT DU TEMPS QUI DECRIT LA DYNAMIQUE DU SYSTEME : LE MODELE DE HARTREE-FOCK NON ADIABATIQUE. LA QUASI-TOTALITE DES REACTIONS CHIMIQUES INTERESSANT L'INDUSTRIE OU LES SCIENCES DE LA VIE SE DEROULENT EN PHASE LIQUIDE, OU LES EFFETS DE SOLVANTS JOUENT UN ROLE DETERMINANT. LES CHAPITRES 5, 6 ET 7 CONCERNENT LA RESOLUTION NUMERIQUE DES MODELES DE CONTINUUM QUI SONT LES MODELES DE SOLVATATION OFFRANT A L'HEURE ACTUELLE LE MEILLEUR COMPROMIS ENTRE QUALITE DES RESULTATS ET TEMPS DE CALCUL.

La première partie du mémoire est consacrée à une étude bibliographique. Après un rapide aperçu historique des travaux portant sur les milieux à microstructure, une étude de la théorie des matériaux simples de W. Noll est présentée, ainsi qu'un exposé des travaux de G. Capriz sur les matériaux à microstructure. On introduit ensuite les notions d'univers d'observation "profond", de microplacement et de macroplacement dans cet univers d'un ensemble de particules. Les propriétés de ces placements conditionnent les définitions des milieux monoparticulaires, multiparticulaires, multiconstituants ou multiphasiques que l'on donne. On décrit les propriétés cinématiques et dynamiques de ces milieux. On propose l'introduction d'un champ de fonctions de phase pour traiter les problèmes de changements de phases. Deux approches complémentaires pour obtenir les équations d'équilibre d'un milieu multiparticulaire à microstructure sont ensuite mises en oeuvre, en s'intéressant plus particulièrement aux milieux de Cosserat multiparticulaires (distribution continue d'orientation de trièdres au même point de l'espace physique). Le chapitre IV consiste en un exposé complet d'une méthode de changement d'échelle apparentée aux outils classiques de l'homogénéisation que nous nommons globalisation. Cette méthode permet, par défocalisation d'un milieu à microstructure "classique", de déduire les équations d'équilibre macro des équations d'équilibre du milieu initial. L'intérêt réside dans l'identification physique des grandeurs constitutives du milieu macro à partir des grandeurs constitutives du milieu micro. La démarche est menée à terme dans le cadre général des microstructures sur lesquelles l'action du groupe des rotations de IR3 est transitive. La fin du mémoire apporte quelques indications sur l'ébauche d'une méthode de détermination du comportement d'un milieu multiparticulaire à partir de la globalisation des principales équations de bilan du milieu micro.