Cette thèse porte sur la caractérisation numérique et la modélisation stochastique de l'accélération du fluide pour l'écoulement en canal à grand nombre de Reynolds. La motivation concerne l'observation et l'analyse des effets de l'intermittence liés aux interactions à longue portée à travers le canal. Dans la première partie, l'accélération est étudiée par simulation numérique directe pour trois différents nombres de Reynolds (180, 590 et 1000). La lognormalité de la norme de l'accélération est observée quelle que soit la distance à la paroi. Un profil universel de la norme de l'accélération est également recherché par analyse dimensionnelle. La seconde partie présente une modélisation stochastique de l'accélération basée sur la décomposition norme/orientation. Le modèle stochastique pour la norme s'appuie sur un processus de fragmentation afin de représenter les interactions à longue portée à travers le canal. Pour l'orientation, l'évolution vers l'isotropie lorsque la distance à la paroi augmente (observée par la DNS) est reproduite grâce à un modèle de marche aléatoire sur une sphère. Ces modèles ont été appliqués à l'approche LES-SSAM (Stochastic Subgrid Acceleration Model) introduite par Sabel'nikov, Chtab et Gorokhovski. Nos calculs montrent que les estimations de la vitesse moyenne, du spectre d'énergie, des contraintes de l'écoulement et de la non-gaussianité des statistiques de l'accélération peuvent être améliorées de façon significative par rapport à la LES classique. L'intérêt de l'approche LES-SSAM, donnant un accès vers la structure intermittente de sous-maille, est illustré dans la dernière partie, par l'étude du transport de particules inertielles ponctuelles par l'écoulement de canal. Cette étude commence par l'analyse par DNS de l'influence des structures de paroi sur la dynamique des particules

Dans une première partie, on s'intéresse à un système d'EDP stochastiques variant selon deux échelles de temps, et plus particulièrement à l'approximation de la composante lente à l'aide d'un schéma numérique efficace. On commence par montrer un principe de moyennisation, à savoir la convergence de la composante lente du système vers la solution d'une équation dite moyennée. Ensuite on prouve qu'un schéma numérique de type Euler fournit une bonne approximation d'un coefficient inconnu apparaissant dans cette équation moyennée. Finalement, on construit et on analyse un schéma de discrétisation du système à partir des résultats précédents, selon la méthodologie dite HMM (Heterogeneous Multiscale Method). On met en évidence l'ordre de convergence par rapport au paramètre d'échelle temporelle et aux différents paramètres du schéma numérique; on étudie les convergences au sens fort (approximation des trajectoires) et au sens faible (approximation des lois). Dans une seconde partie, on étudie une méthode d'approximation de solutions d'EDP paraboliques, en combinant une approche semi-lagrangienne et une discrétisation de type Monte-Carlo. On montre d'abord dans un cas simplifié que la variance dépend des pas de discrétisation; enfin on fournit des simulations numériques de solutions, afin de mettre en avant les applications possibles d'une telle méthode.

Les structures planaires de grandes tailles sont de plus en plus utilisées dans les applications des satellites et des radars. Deux grands types de ces structures à savoir les FSS et les Reflectarrays sont particulièrement les plus intéressants dans les domaines de la conception RF. Mais en raison de leur grande taille et de la complexité des cellules élémentaires, l'analyse complète de ces structures nécessite énormément de mémoire et des temps de calcul excessif. Par conséquent, les techniques classiques basées sur maillage linéaire soit ne parviennent pas à simuler de telles structures soit, exiger des ressources non disponibles à un concepteur d'antenne. Une technique appelée " technique par changement d'échelle " tente de résoudre ce problème par partitionnement de la géométrie du réseau par de nombreux domaines imbriqués définis à différents niveaux d'échelle du réseau. Le multi-pôle par changement d'échelle, appelé " Scale changing Network (SCN) ", modélise le couplage électromagnétique entre deux échelles successives, en résolvant une formulation intégral des équations de Maxwell par une technique basée sur la méthode des moments. La cascade de ces multi-pôles par changement d'échelle, permet le calcul de la matrice d'impédance de surface de la structure complète qui peut à son tour être utilisées pour calculer la diffraction en champ lointain. Comme le calcul des multi-pôles par changement d'échelle est mutuellement indépendant, les temps d'exécution peuvent être réduits de manière significative en parallélisant le calcul. Par ailleurs, la modification de la géométrie de la structure à une échelle donnée nécessite seulement le calcul de deux multi-pôles par changement d'échelle et ne requiert pas la simulation de toute la structure. Cette caractéristique fait de la SCT un outil de conception et d'optimisation très puissant. Des structures planaires uniformes et non uniformes excité par un cornet ont étés modélisés avec succès, avec des temps de calcul dé lais intéressants, employant les ressources normales de l'ordinateur.

Les piles à combustible à membrane échangeuse de protons (PEMFC) sont sensibles aux polluantsde l'hydrogène et de l'air. Cette étude s'est focalisée sur l'impact du monoxyde de carbone (CO) et dusulfure d'hydrogène (H2S), deux polluants majeurs dans l'hydrogène (H2), suivant une approchecombinant expériences et modélisation.Le volet expérimental a consisté à étudier l'effet de la concentration des polluants individuels et enmélange et des chargements en catalyseurs, pour différents modes de fonctionnement. Cette étude amis en lumière un impact sur les deux électrodes (anode et cathode) dû à la distribution hétérogènedes polluants à la surface de l'anode et à la désactivation de la partie de la cathode en regard. Deplus, dans le cas d'un empoisonnement par H2S, cette étude a montré que la tension de cellule atteintun état quasi-stationnaire, en mode galvanostatique, ce qui n'avait jamais été mis en évidence dans lalittérature.Dans l'approche de modélisation multi-échelles, le couplage de l'électrochimie et de la fluidique ainsique le développement de différentes " briques " du modèle ont permis de perfectionner la descriptiondes phénomènes physico-chimiques. Le modèle permet maintenant de simuler le fonctionnementd'une cellule de pile à combustible dans les conditions opératoires réelles, en intégrant les cinétiquesd'empoisonnement du platine par CO et H2S.Enfin, la comparaison des données expérimentales et des simulations a montré des résultats trèssatisfaisants appuyant certains arguments pour l'interprétation de l'impact des impuretés de H2.

L'objet de cette thèse est l'étude mathématique de modèles probabilistes pour la génération et la propagation d'un potentiel d'action dans les neurones et plus généralement de modèles aléatoires pour les systèmes excitables. En effet, nous souhaitons étudier l'influence du bruit sur certains systèmes excitables multi-échelles possédant une composante spatiale, que ce soit le bruit contenu intrinsèquement dans le système ou le bruit provenant du milieu. Ci-dessous, nous décrivons d'abord le contenu mathématique de la thèse. Nous abordons ensuite la situation physiologique décrite par les modèles que nous considérons. Pour étudier le bruit intrinsèque, nous considérons des processus de Markov déterministes par morceaux à valeurs dans des espaces de Hilbert ("Hilbert-valued PDMP"). Nous nous sommes intéressés à l'aspect multi-échelles de ces processus et à leur comportement en temps long. Dans un premier temps, nous étudions le cas où la composante rapide est une composante discrète du PDMP. Nous démontrons un théorème limite lorsque la composante rapide est infiniment accélérée. Ainsi, nous obtenons la convergence d'une classe de "Hilbert-valued PDMP" contenant plusieurs échelles de temps vers des modèles dits moyennés qui sont, dans certains cas, aussi des PDMP. Nous étudions ensuite les fluctuations du modèle multi-échelles autour du modèle moyenné en montrant que celles-ci sont gaussiennes à travers la preuve d'un théorème de type "central limit". Dans un deuxième temps, nous abordons le cas où la composante rapide est elle-même un PDMP. Cela requiert de connaître la mesure invariante d'un PDMP à valeurs dans un espace de Hilbert. Nous montrons, sous certaines conditions, qu'il existe une unique mesure invariante et la convergence exponentielle du processus vers cette mesure. Pour des PDMP dits diagonaux, la mesure invariante est explicitée. Ces résultats nous permettent d'obtenir un théorème de moyennisation pour des PDMP "rapides" couplés à des chaînes de Markov à temps continu "lentes". Pour étudier le bruit externe, nous considérons des systèmes d'équations aux dérivées partielles stochastiques (EDPS) conduites par des bruits colorés. Sur des domaines bornés de $\mathbb{R}^2$ ou $\mathbb{R}^3$, nous menons l'analyse numérique d'un schéma de type différences finies en temps et éléments finis en espace. Pour une classe d'EDPS linéaires, nous étudions l'erreur de convergence forte de notre schéma. Nous prouvons que l'ordre de convergence forte est deux fois moindre que l'ordre de convergence faible. Par des simulations, nous montrons l'émergence de phénomènes d'ondes ré-entrantes dues à la présence du bruit dans des domaines de dimension deux pour les modèles de Barkley et Mitchell-Schaeffer.

L'incidence rasante est un problème spécifique, qui apparaît notamment lorsqu'une antenne est placée sur un mât (télécommunications, défense...) ou sur la côte (surveillance environnementale ou militaire de l'espace maritime). Elle rend la modélisation du problème de diffraction difficile, de par le faible niveau de rétrodiffusion et l'importance de phénomènes complexes comme la diffusion multiple. La question reste importante même si l'écho est très faible, puisqu'il est potentiellement suffisant pour perturber le bon fonctionnement de systèmes antennaires microondes sur un navire. Il porte de plus des informations intéressantes sur l'état de la mer, comme cela a été démontré aux bandes HF et VHF. Un modèle rigoureux de diffraction tridimensionnelle précédemment développé est étendu au calcul des quatre polarisations fondamentales (polarisations Horizontale et Verticale des ondes incidente et réfléchie). Il permet désormais de prendre en compte la conductivité finie de la surface, point crucial dans le cas d'une polarisation incidente Verticale. L'opérateur hyper-singulier impliqué dans l'équation intégrale discrétisée par la méthode des moments est étudié pour évaluer la précision des calculs numériques.Les méthodes approchées de diffraction permettent des calculs numériques beaucoup plus rapides, et sont donc en pratique incontournables. Le modèle rigoureux est donc utilisé, en conjonction avec des données expérimentales, pour servir de référence permettant d'étudier la précision, en incidence rasante et dans le cas de la surface de la mer, de ces méthodes approchées. Nous étudions en particulier la méthode à deux échelles GOSSA, et proposons une correction à son comportement aux angles rasants.Le mouvement de la surface de la mer crée un décalage de fréquence radar dans l'onde rétrodiffusée (effet Doppler), décalage mesuré expérimentalement et que l'algorithme de méthode des moments permet de simuler. Nous étudions par des simulations bidimensionnelles l'évolution du décalage Doppler micro onde avec l'incidence, et l'influence des nonlinéarités de la surface de la mer. Le comportement limite en incidence rasante est précisé, et les contributions respectives des phénomènes électromagnétiques et hydrodynamiques discutées.

Une méthode de frontières immergées est développée pour la simulation d'écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d'ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d'une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l'adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d'un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L'étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l'écoulement principal et présence ou non d'un swirl). Les spécificités de mise en œuvre de la simulation d'un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d'instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L'existence d'un cycle limite acoustique est souligné et une formule permettant d'anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d'assurer une reproduction plus précise de la dynamique de flamme.

Cette thèse est consacrée à l'amélioration des capacités de prédiction des émissions polluantes (CO, NOx . . . ) des foyers de combustion de turboréacteurs. L'étude, exclusivement numérique, repose sur des simulations aux grandes échelles (ou LES pour Large-Eddy Simulation) basées sur des méthodes de tabulation de la chimie détaillée. L'approche PCM-FPI (pour Presumed Conditional Moments - Flame Prolongation of ILDM) a été étendue à la chimie des oxydes d'azote via la modification de la variable d'avancement. Différentes validations sur des configurations laminaires simples puis des flammes de laboratoire (Cabra, Sandia) sont présentées. Les résultats en terme de structure de flamme et de champs d'espèces chimiques sont confrontés aux mesures expérimentales. Le rôle du formaldéhyde comme marqueur de la zone réactionnelle est illustré à l'aide de calculs de flammes laminaires puis confirmé par un calcul 3D LES. Une analyse des spécificités de l'implantation de ce type de modèle sur des machines à architecture massivement parallèle est ensuite menée. Diverses modifications de la structure de la table et des méthodes d'interpolation sont réalisées, servant de base à une étude de sensibilité de maillage appliquée à la flamme Sandia D. Les difficultés relatives à la prédiction du NO dans les flammes turbulentes sont exposées : divers modèles de sous-maille sont alors employés et comparés.

Cette thèse porte sur la simulation des jets de gouttes générés par des pulvérisateurs essence haute pression, pulvérisateurs qui sont un point clef des systèmes de combustion automobile de la présente et future génération devant diminuer les émissions de CO2 et de polluants. Dans un premier temps les jets de gouttes (" sprays ") sont simulés par simulation moyennée. Les résultats de simulation d'un jet donnant des résultats en moyenne satisfaisant, l'interaction de jets en injecteurs multi-trous est alors simulée. Les résultats sont cohérents par rapport aux mesures d'entraînement d'air. La simulation permettant d'avoir accès au champ complet 3D, le mécanisme d'interaction jet à jet et de développement instationnaire du spray est décrit en détail. La formation d'un mouvement descendant au centre du spray et celle d'un point d'arrêt central sont trouvés. Finalement, Ces résultats sont étendus au cas surchauffé, cas où la pression dans la chambre est inférieure à la pression de vapeur saturante. Un modèle simple semi-empirique est proposé pour tenir compte de la modification des conditions proches de la buse d'injection. Le modèle prédit correctement les tendances des variations de paramètres et capture la forme générale du spray qui se referme sur lui-même. La seconde grande partie est consacrée au développement d'un modèle de spray par l'approche des grandes échelles (SGE), limité ici aux cas non évaporant. Il comprend la modélisation de sous-maille de la dispersion turbulente, des collisions-coalescence et des termes d'échange de quantité de mouvement de sous-maille. L'effet du choix du modèle de sous-maille pour la viscosité turbulente de sous-maille est montré, le choix retenu étant le modèle de Smagorinski dynamique. Afin d'améliorer la représentativité cruciale des conditions d'injections, un couplage faible est réalisé à partir de résultats de simulations existantes de l'écoulement interne aux buses. Les fonctions densité de probabilité simple et jointes extraits des résultats de simulations sont validés par rapport aux mesures PDA en situation pseudo-stationnaire et la pénétration liquide et la forme du spray est comparée aux visualisations par ombroscopie. Enfin, différentes zones caractéristiques sont identifiées et des longueurs sont notées pour les cas d'injection à 100 et 200bar.

La combustion stratifiée, qui consiste à brûler un mélange carburant/oxydant inhomogène, et la combustion diluée, consistant à ajouter une quantité limitée de gaz brûlés, sont deux technologies utilisées dans les moteurs à piston pour réduire leur consommation. Cette thèse est dédiée à l'étude de l'allumage dans ces deux types de milieux en régimes laminaire et turbulent. Un nouveau schéma cinétique pour la combustion propane/air a été dérivé et combiné à deux approches de modélisation différentes : la chimie complexe et une approche de chimie tabulée de type FPI. Dans le cas laminaire, les deux approches de modélisation donnent des résultats similaires et un modèle simple a mis en évidence l'importance de la dynamique des gaz frais et des gaz brûlés sur le développement du noyau. Dans le cas turbulent, plusieurs techniques d'analyse ont montré la dépendance de la vitesse absolue de la flamme au champ de vitesse moyen et la décorrélation des fluctuations locales de richesse.