Dans ce travail, nous nous intéressons à des problèmes issus de la Mécanique des Fluides et plus particulièrement au cas des aérosols (ou sprays, c'est-à-dire un ensemble de particules en suspension dans un fluide environnant). Les phénomènes physiques mis en jeu sont modélisés par des équations aux dérivées partielles (EDP). La phase continue (fluide environnant) est décrite par des équations issues de la mécanique des milieux continus de type Navier-Stokes ou Euler. La phase dispersée est décrite par une équation cinétique de type Boltzmann.<br /><br />Le premier résultat que nous présentons est consacré à l'étude mathématique d'un couplage entre une équation cinétique de type Vlasov et les équations d'Euler isentropiques. Ces équations modélisent un spray fin. Nous démontrons l'existence en temps petit d'une solution régulière pour le couplage Vlasov-Euler isentropique.<br /><br />Ensuite, nous présentons les équations précises relatives à la modélisation des collisions, coalescences et fragmentations dans un spray.<br /><br />Nous décrivons par la suite la simulation numérique du couplage fluide-cinétique dans un code industriel (Commissariat à l'Énergie Atomique), en particulier l'ajout des phénomènes de collisions.<br /><br />Un deuxième modèle de fragmentation est également présenté. Ce modèle est plus pertinent dans les cas où les particules de la phase dispersée ont un grand nombre de Weber.<br /><br />Enfin, nous présentons un résultat concernant une estimation explicite de trou spectral pour l'opérateur de Boltzmann avec potentiels durs linéarisé, et pour l'opérateur de Landau avec potentiels durs linéarisé.

Un grand nombre de systèmes physiques sont le siège d'interactions à longue portée : systèmes auto gravitants, plasmas, interactions entre vortex... et partagent de ce fait certaines propriétés. Dans cette thèse, un autre type de système expérimental est envisagé : des atomes froids ; dans ce cas, ce sont les échanges de photons qui peuvent induire des interactions effectives à longue portée. La dynamique de ces systèmes à longue portée est décrite sur une certaine échelle de temps par une équation de type Vlasov, ou Vlasov-Fokker-Planck. Le but de cette thèse est d'étudier le comportement hors équilibre de plusieurs systèmes de particules comportant en général des interactions à longue portée, d'un point de vue théorique, numérique et expérimental. Dans une première partie, nous étudions dans le cadre de l'équation de Vlasov la dynamique d'un système de particules au voisinage d'un état stationnaire inhomogène. Nous montrons que si un amortissement de type Landau apparaît aux temps courts, une relaxation vers un état stationnaire en loi de puissance domine toujours aux temps longs. Nous testons et validons ensuite nos prédictions par des simulations numériques du modèle HMF (archétype des systèmes à longue portée). Nous nous intéressons ensuite aux oscillations de respiration et du centre de masse d'un système de particules en interaction. En supposant une invariance de la forme de la distribution des particules, nous obtenons deux équations qui décrivent approximativement l'évolution de ces modes pour une grande gamme de systèmes (longue/courte portée, avec/sans thermostat, ...). Pour finir, nous présentons l'utilisation des résultats précédemment obtenus pour explorer un régime d'instabilité dans un piège magnéto-optique, et la possible existence de l'analogue d'un régime auto-gravitant 1d.

Une des difficultés fondamentales de la simulation numérique de plasmas magnétisés est l'existence d'échelles de temps et d'espace mettant en jeu plusieurs ordres de grandeurs très différents. Afin de réaliser des simulations numériques efficaces de ces phénomènes physiques, il est essentiel de développer des modèles et des méthodes numériques adaptés à ces problèmes. <br />A ce jour, la notion de convergence 2-échelles introduite par G. Allaire et G. Nguetseng est un des outils permettant de dériver rigoureusement des limites multi-échelles, ce qui nous permet d'obtenir des modèles limites qu'il est possible de discrétiser avec une méthode numérique usuelle : nous parlons alors d'une méthode numérique 2-échelles. <br />L'objectif de cette thèse est de développer une méthode semi-lagrangienne 2 échelles sur un modèle de type Vlasov gyrocinétique afin de simuler un plasma soumis à un champ magnétique fort du même type que ceux utilisés pour le projet ITER. Cependant, comme les phénomènes physiques à simuler sont assez complexes et comme nous ne savons que peu de choses sur le comportement d'une méthode numérique 2-échelles sur un modèle non-linéaire, il convient de procéder par étapes avant de développer une telle méthode sur un modèle gyrocinétique. <br />Dans une première partie, nous construisons une méthode de volumes finis 2-échelles sur les équations d'Euler 1D isentropiques faiblement compressibles. Bien que ce modèle soit assez différent d'un modèle de type Vlasov, il n'en est pas moins un cadre de travail relativement simple pour étudier le comportement d'une méthode numérique 2-échelles face à un modèle non-linéaire. <br />Dans une seconde partie, nous nous basons sur le modèle limite développé par E. Frénod, F. Salvarani et E. Sonnendrücker afin de construire une méthode semi-lagrangienne 2-échelles pour simuler des faisceaux de particules en géométrie axisymétrique. Même si le modèle de Vlasov axisymétrique utilisé est différent d'un modèle gyrocinétique, il constitue un contexte idéal pour établir les bases d'une méthode semi-lagrangienne 2 échelles.<br />Enfin, dans une troisième partie, nous utilisons la convergence 2-échelles afin d'améliorer les résultats de convergence faible-* établis par M. Bostan en 2007, et nous proposons une méthode semi-lagrangienne en avant permettant de valider numériquement ces résultats.

Nuclear spinodal instabilities are investigated in non-relativistic and relativistic stochastic mean-field approaches for charge asymmetric and charge symmetric nuclear matter. Quantum statistical effect on the growth of instabilities are calculated in non-relativistic approach. Due to quantal effects, in both symmetric and asymmetric matter, dominant unstable modes shift towards longer wavelengths and modes with wave numbers larger than the Fermi momentum are strongly suppressed. As a result of quantum statistical effects, in particular at lower temperatures, amplitude of density fluctuations grows larger than those calculated in semi-classical approximation. Relativistic calculations in the semi-classical limit are compared with the results of non-relativistic calculations based on Skyrme-type effective interactions under similar conditions. A qualitative difference appears in the unstable response of the system: the system exhibits most unstable behavior at higher baryon densities around $rho_{B}=0.4 rho_{0}$ in the relativistic approach while most unstable behavior occurs at lower baryon densities around $rho_{B}=0.2 rho_{0}$ in the non-relativistic calculations.

The nuclear matter liquid-gas phase transition is expected to be a signal of nuclear spinodal instabilities as a result of density fluctuations. Nuclear spinodal instabilities in symmetric nuclear matter are studied within a stochastic relativistic density-dependent model in semi-classical approximation. We use two parameterization for the Lagrange density, DDME1 and TW sets. The early growth of density fluctuations is investigated by employing relativistic Vlasov equation based on QHD and discussed the cluster size of the condensations from the early growth of density correlation functions. Expectations are that hot nuclear matter behaves unstable around &rho / b &asymp / &rho / 0/4 (below the saturation density) and at low temperatures. We therefore present our results at low temperature T=1 MeV and at higher temperature T=5 MeV, and also at a lower initial baryon density &rho / b = 0.2 &rho / 0 and a higher value &rho / b = 0.4 &rho / 0 where unstable behavior is within them. Calculations in density-dependent model are compared with the other calculations obtained in a relativistic non-linear model and in a Skyrme type nonivrelativistic model. Our results are consistent with them. Qualitatively similar results show that the physics of the quantities are model-independent. The size of clusterization is estimated in two ways, by using half-wavelength of the most unstable mode and from the width of correlation function at half maximum. Furthermore, the average speed of condensing fragments during the initial phase of spinodal decomposition are determined by using the current density correlation functions.

Les résultats présentés dans cette thèse portent sur l'approximation adaptative de deux problèmes de transport non-linéaire : le système de Vlasov-Poisson et les lois de conservation scalaires. Pour le premier, et dans une approche semi-lagrangienne, on a proposé un schéma adaptatif original à base d'éléments finis hiérarchiques où l'évolution des maillages est réalisée par une étape de prédiction très simple suivie d'une étape de correction plus classique. En introduisant la notion de courbure totale pour étendre la semi-norme W2,1(R2) aux fonctions affines par morceaux, on a alors établi une estimation d'erreur a priori prouvant la convergence de ce schéma en distance L∞, et donné des éléments de preuve concernant sa complexité optimale. Les lois de conservations scalaire ne pouvant être approchées en distance L∞, on a considéré leur analyse en distance uniforme de Hausdorff, moins répandue bien que plus géométrique. Après avoir montré que les solutions de ces équations étaient stables vis-à-vis de cette distance, on a établi un résultat d'approximation adaptative d'ordre élevé.

Cette thèse traite de deux sujets: la résolubilité forte d'équations différentielles stochastiques à dérive hölderienne et bruit hypoelliptique et la simulation de processus progressifs-rétrogrades découplés de McKean-Vlasov. Dans le premier cas, on montre qu'un système hypoelliptique, composé d'une composante diffusive et d'une composante totalement dégénérée, est fortement résoluble lorsque l'exposant de la régularité Hölder de la dérive par rapport à la composante dégénérée est strictement supérieur à 2/3. Ce travail étend au cadre dégénéré les travaux antérieurs de Zvonkin (1974), Veretennikov (1980) et Krylov et Röckner (2005). L'apparition d'un seuil critique pour l'exposant peut-être vue comme le prix à payer pour la dégénérescence. La preuve repose sur des résultats de régularité de la solution de l'EDP associée, qui est dégénérée, et est basée sur une méthode parametrix. Dans le second cas, on propose un algorithme basé sur les méthodes de cubature pour la simulation de processus progessifs-rétrogrades découplés de McKean-Vlasov apparaissant dans des problèmes de contrôle dans un environnement de type champ moyen. Cet algorithme se divise en deux parties. Une première étape de construction d'un arbre de particules, à dynamique déterministe, approchant la loi de la composante progressive. Cet arbre peut être paramétré de manière à obtenir n'importe quel ordre d'approximation (en terme de pas de discrétisation de l'intervalle). Une seconde étape, conditionnelle à l'arbre, permettant l'approximation de la composante rétrograde. Deux schémas explicites sont proposés permettant un ordre d'approximation de 1 et 2.

Cette thèse traite de deux sujets: la résolubilité forte d'équations différentielles stochastiques à dérive hölderienne et bruit hypoelliptique et la simulation de processus progressifs-rétrogrades découplés de McKean-Vlasov. Dans le premier cas, on montre qu'un système hypoelliptique, composé d'une composante diffusive et d'une composante totalement dégénérée, est fortement résoluble lorsque l'exposant de la régularité Hölder de la dérive par rapport à la composante dégénérée est strictement supérieur à 2/3. Ce travail étend au cadre dégénéré les travaux antérieurs de Zvonkin (1974), Veretennikov (1980) et Krylov et Röckner (2005). L'apparition d'un seuil critique pour l'exposant peut-être vue comme le prix à payer pour la dégénérescence. La preuve repose sur des résultats de régularité de la solution de l'EDP associée, qui est dégénérée, et est basée sur une méthode parametrix. Dans le second cas, on propose un algorithme basé sur les méthodes de cubature pour la simulation de processus progessifs-rétrogrades découplés de McKean-Vlasov apparaissant dans des problèmes de contrôle dans un environnement de type champ moyen. Cet algorithme se divise en deux parties. Une première étape de construction d'un arbre de particules, à dynamique déterministe, approchant la loi de la composante progressive. Cet arbre peut être paramétré de manière à obtenir n'importe quel ordre d'approximation (en terme de pas de discrétisation de l'intervalle). Une seconde étape, conditionnelle à l'arbre, permettant l'approximation de la composante rétrograde. Deux schémas explicites sont proposés permettant un ordre d'approximation de 1 et 2.

Cette thèse traite de la solution numérique de deux types de problèmes stochastiques. Premièrement, nous nous intéressons aux EDS fortement oscillantes, c'est-à-dire, les systèmes composés de variables ergodiques évoluant rapidement par rapport aux autres. Nous proposons un algorithme basé sur des résultats d'homogénéisation. Il est défini par un schéma d'Euler appliqué aux variables lentes couplé avec un estimateur à pas décroissant pour approcher la limite ergodique des variables rapides. Nous prouvons la convergence forte de l'algorithme et montrons que son erreur normalisée satisfait un résultat du type théorème limite centrale généralisé. Nous proposons également une version extrapolée de l'algorithme ayant une meilleure complexité asymptotique en satisfaisant les mêmes propriétés que la version originale. Ensuite, nous étudions la solution des EDS de type McKean-Vlasov (EDSPR-MKV) associées à la solution de certains problèmes de contrôle sous un environnement formé d'un grand nombre de particules ayant des interactions du type champ-moyen. D'abord, nous présentons un nouvel algorithme, basé sur la méthode de cubature sur l'espace de Wiener, pour approcher faiblement la solution d'une EDS du type McKean-Vlasov. Il est déterministe et peut être paramétré pour atteindre tout ordre de convergence souhaité. Puis, en utilisant ce nouvel algorithme, nous construisons deux schémas pour résoudre les EDSPR-MKV découplées et nous montrons que ces schémas ont des convergences d'ordres un et deux. Enfin, nous considérons le problème de réduction de la complexité de la méthode présentée tout en respectant la vitesse de convergence énoncée.

The Vlasov-Poisson system is most commonly used to model the movement of charged particles in a plasma or of stars in a galaxy. It consists of a kinetic equation known as the Vlasov equation coupled with a force determined by the Poisson equation. The system in Euclidean space is well-known and has been extensively studied under various assumptions. In this paper, we derive the Vlasov-Poisson equations assuming the particles exist only on the 2-sphere, then take an in-depth look at particles which initially lie along a great circle of the sphere. We show that any great circle is an invariant set of the equations of motion and prove that the total energy, number of particles, and entropy of the system are conserved for circular initial distributions. / Graduate