Cette thèse est consacrée à l'étude des équations et systèmes de réaction-diffusion dans des milieux hétérogènes. Elle est divisée en deux parties. La première est dédiée à l'étude des équations de réaction-diffusion dans des milieux périodiques. Nous nous intéressons en particulier aux équations posées dans des domaines qui ne sont pas l'espace entier $\mathbb{R}^{N}$, mais des domaines périodiques, avec des "obstacles". Dans un premier chapitre, nous étudions l'effet de la géométrie du domaine sur la vitesse d'invasion des solutions. Après avoir dérivé une formule de type Freidlin-Gartner, nous construisons des domaines où la vitesse d'invasion est strictement inférieure à la vitesse critique des fronts. Nous donnons également des critères géométriques qui garantissent l'existence de directions où l'invasion se produit à la vitesse critique. Dans le chapitre suivant, nous donnons des conditions nécessaires et suffisantes pour garantir que l'invasion ait lieu, après quoi nous construisons des domaines où des phénomènes intermédiaires (blocage, invasion orientée) se produisent. La deuxième partie de cette thèse est consacrée à l'étude de modèles décrivant l'influence de lignes à diffusion rapide (une route, par exemple) sur la propagation d'espèces invasives. Il a en effet été observé que certaines espèces, dont le moustique-tigre, envahissent plus rapidement que prévu certaines zones proches du réseau routier. Nous étudions deux modèles : le premier décrit l'influence d'une route courbe sur la propagation. Nous nous intéressons en particulier au cas de deux routes non-parallèles. Le second modèle décrit l'influence d'une route sur une niche écologique, en présence d'un changement climatique. Le résultat principal est que l'effet de la route est ambivalent : si la niche est stationnaire, alors l'effet de la route est délétère. Cependant, si la niche se déplace, suite à un changement climatique, nous montrons que la route peut permettre à une population de survivre. Pour étudier ce second modèle, nous développons une notion de valeur propre principale généralisée pour des systèmes de type KPP, et nous dérivons une inégalité de Harnack, qui est nouvelle pour ce type de systèmes. / This thesis is dedicated to the study of reaction-diffusion equations and systems in heterogeneous media. It is divided into two parts. The first one is devoted to the study of reaction-diffusion equations in periodic media. We pay a particular attention to equations set on domains that are not the whole space $\mathbb{R}^{N}$, but periodic domains, with "obstacles". In a first chapter, we study how the geometry of the domain can influence the speed of invasion of solutions. After establishing a Freidlin-Gartner type formula, we construct domains where the speed of invasion is strictly less than the critical speed of fronts. We also give geometric criteria to ensure the existence of directions where the invasion occurs with the critical speed. In the second chapter, we give necessary and sufficient conditions to ensure that invasion occurs, and we construct domains where intermediate phenomena (blocking, oriented invasion) occur. The second part of this thesis is dedicated to the study of models describing the influence of lines with fast diffusion (a road, for instance) on the propagation of invasive species. Indeed, it was observed that some species, such as the tiger mosquito, invade faster than expected some areas along the road-network. We study two models : the first one describes the influence of a curved road on the propagation. We study in particular the case of two non-parallel roads. The second model describes the influence of a road on an ecological niche, in presence of climate change. The main result is that the effect of the road is ambivalent: if the niche is stationary, then effect of the road is deleterious. However, if the niche moves, because of a shifting climate, the road can actually help the population to persist. To study this model, we introduce a notion of generalized principal eigenvalue for KPP-type systems, and we derive a Harnack inequality, that is new for this type of systems.

La théorie des jeux à champ moyen fut introduite en 2006 par Jean-Michel Lasry et Pierre-Louis Lions. Elle permet l'étude de la théorie des jeux dans certaines configurations où le nombre de joueurs est trop grand pour espérer une résolution pratique. Nous étudions la théorie des jeux à champ moyen sur les graphes en nous appuyant sur les travaux d'Olivier Guéant que nous étendrons à des formes plus générales d'Hilbertien. Nous étudierons aussi les liens qui existent entres les K-moyennes et les jeux à champ moyen ce qui permettra en principe de proposer de nouveaux algorithmes pour les K-moyennes grâce aux techniques de résolution numérique propres aux jeux à champ moyen. Enfin nous étudierons un jeu à champ moyen à savoir le problème "d'heure de début d'une réunion" en l'étendant à des situations où les agents peuvent choisir entre deux réunions. Nous étudierons de manière analytique et numérique l'existence et la multiplicité des solutions de ce problème. / The mean field game theory was introduced in 2006 by Jean-Michel Lasry and Pierre-Louis Lions. It allows us to study the game theory in some situations where the number of players is too high to be able to be solved in practice. We will study the mean field game theory on graphs by learning from the studies of Oliver Guéant which we will extend to more generalized forms of Hilbertian. We will also study the links between the K-means and the mean field game theory. In principle, this will offer us new algorithms for solving the K-means thanks to the techniques of numerical resolutions of the mean field games. Findly, we will study a mean field game called the "starting time of a meeting". We will extend it to situations where the players can choose between two meetings. We will study analytically and numerically the existence and multiplicity of the solutions to this problem.

Les modèles structure-fonction de la croissance des plantes (FSPMs) combinent la description du fonctionnement biophysique et du développement architectural des plantes. On peut distinguer deux grandes familles de FSPM : d'une part les modèles décrivant finement la structure de la plante au niveau de l'organe et d'autre part les modèles à plus grande échelle qui s'intéressent directement à la forme du houppier. La paramétrisation du premier type de modèle est souvent difficile car elle nécessite des données expérimentales très riches. A l'inverse, les modèles à plus grande échelle mettent généralement en œuvre des lois empiriques qui ne permettent pas de décrire la plasticité de la croissance, et l'adaptation de la plante à des conditions environnementales différentes.Pour répondre à ces problématiques, nous nous tournons vers un nouveau paradigme : Motivé par le succès du concept de la densité spatiale dans les modèles en écologie des populations, cette thèse caractérise la distribution spatiale de feuillage dans les plantes par la densité de surface foliaire , ce qui permet une description locale ouvrant la voie à une prise en compte de la plasticité des plantes, tout en ne décrivant pas chaque feuille individuellement, ce qui permet de modéliser des vieux et grands arbres, dont le nombre de feuilles est sinon trop lourd à gérer du point de vue des calculs. Cette thèse présente des modèles dynamiques de croissance développés spécifiquement pour les plantes agricoles et les arbres. Nous explorons des approches mathématiques différentes en temps discrète et continue, tout en examinant d'un œil critique leurs aptitudes conceptuelles ainsi que des possibilités de simplifications et de solutions analytiques dans l'optique de l'accélération des simulations.La densité foliaire permet le calcul de l'interception de lumière par la loi de Beer-Lambert et la production de biomasse grâce au concept d'efficience d'utilisation de la lumière. Le mécanisme central qui est considéré pour les différentes approches développées dans cette thèse est celui de l'expansion locale de la surface foliaire dans la direction du gradient de lumière. Par ce concept téléonomique, nous faisons l'hypothèse que la plante cherche par son développement à optimiser la productivité de la surface foliaire pour la production de biomasse. Ce principe induit ainsi un développement horizontal et vertical du feuillage vers l'extérieur du houppier. Le développement horizontal cesse quand on s'approche trop de plantes voisines, leur ombrage diminuant le gradient de lumière et donc l'expansion de densité de surface foliaire dans ces directions. Le modèle de production de biomasse est également généralisé pour une prise en compte explicite de la teneur en eau du sol en introduisant une composante hydraulique permettant de décrire l'équilibre mécaniste entre le potentiel hydrique dans les feuilles et la transpiration par la régulation stomatale. Finalement, nous prenons en compte l'allocation de biomasse produite à d'autres compartiments de la plante tels que les racines et le bois selon la théorie du « pipe model ».Les résultats des modèles sont comparés à un large jeu de données expérimentales sur des plantations à différentes densités et conditions environnementales. Celui-ci montre de remarquables capacités d'une part à prévoir les variables biométriques importantes (hauteur, diamètre du tronc) ainsi que certaines relations d'allométrie, et d'autre part à générer des formes de houppier en accord avec les formes observées, ceci pour les différents scénarios de compétition et comme propriété émergente du modèle. Ainsi, cette thèse démontre le potentiel du concept de densité de surface foliaire en modélisation de la croissance des plantes, par sa capacité à reproduire les comportements locaux et l'adaptation à des conditions environnementales variées sans compromettre l'efficacité et la robustesse. / Functional-structural plant growth models (FSPMs) have emerged as the synthesis of mechanistic process-based models, and geometry-focussed architectural models. In terms of spatial scale, these models can essentially be divided into small-scale models featuring a topologic­al architecture – often facing data-demanding parametrisations, parameter sensitivity, as well as computational heaviness, which imposes problematic limits to the age and size of in­dividuals than can be simulated – and large-scale models based on a description of crown shape in terms of rigid structures such as empirical crown envelopes – commonly strug­gling to allow for spatial variability and plasticity in crown structure and shape in response to local biotic or abiotic growth conditions.In response to these limitations, and motivated not least by the success-story of spatial density approaches in theoretical populations ecology, the spatial distribution of foliage in plants in this thesis is characterised in terms of spatial leaf density, which allows for a com­pletely local description that is a priori unrestricted in terms of plasticity, while being robust and computationally efficient. The thesis presents dynamic growth models specific­ally developed for crops and trees, exploring different mathematical frameworks in continu­ous and discrete time, while critically discussing their conceptual suitability and exploring analyt­ical simplifications and solutions to accelerate simulations.The law of Beer-Lambert on the passing of light though an absorbing medium allows to infer the local light conditions based on which local biomass production can be computed via a radiation use efficiency. A key unifying mechanism of the different models is the local ex­pansion of leaf density in the direction of the light gradient, which coincides with the direc­tion most promising with regard to future biomass productivity. This aspect falls into the line of teleonomic and optimization-oriented plant growth models, and allows to set aside the otherwise complex modelling of branching processes. The principle induces an expans­ive horizontal and upward-directed motion of foliage. Moreover, it mechanistically accounts for a slow-down of the horizontal expansion as soon as a neighbouring competitor's crown is reached, since the appropriate region is already shaded, implying a corresponding adapta­tion of the light gradient. This automatically results in narrower crowns in scenarios of in­creased competition, ultimately decreasing biomass production and future growth due to lesser amount of intercepted light. In an extension, the impact of water availability is incor­porated into the previously light-only dependency of biomass production by means of a novel hydraulic model describing the mechanistic balancing of leaf water potential and tran­spiration in the context of stomatal control. The allocation of produced biomass to other plant compartments such as roots and above-ground wood, e.g. by means of the pipe model theory, is readily coupled to leaf density dynamics.Simulation results are compared against a variety of empirical observations, ranging from long-term forest inventory data to laser-recorded spatial data, covering multiple abi­otic environmental conditions and growth resources as well as stand densities and thus de­grees of competition. The models generate a series of complex emergent properties includ­ing the realistic prediction of biometric growth parameters, the spontaneous adaptability and plasticity of crown morphologies in different competitive scenarios, the empirically documented insensitivity of height to stand density, the accurate deceleration of height growth, as well as popular allometric relationships – altogether demonstrating the potential of leaf density based approaches for efficient and robust plant growth modelling.

Cette thèse porte sur l'étude de méthodes numériques pour l'analyse de deux modèles stochastiques apparaissant dans le contexte de la condensation de Bose-Einstein. Ceux-ci constituent deux généralisations de l'équation de Gross-Pitaevskii. Cette équation aux dérivées partielles déterministe modélise la dynamique de la fonction d'onde d'un condensat de Bose-Einstein piégé par un potentiel extérieur confinant.Le premier modèle étudié permet de modéliser les fluctuations de l'intensité du potentiel confinant et prend la forme d'une équation aux dérivées partielles stochastiques. Celles-ci conduisent en pratique à un échauffement du condensat, et parfois mêmeà son effondrement. Nous proposons dans un premier chapitre la construction d'un schéma numérique pour la résolution de ce modèle. Il est fondé sur une discrétisation spectrale en espace, et une discrétisation temporelle de type Crank-Nicolson. Nous démontrons que le schéma proposé converge fortement en probabilité à l'ordre au moins 1 en temps, et nous présentons des simulations numériques illustrant ce résultat. Le deuxième chapitre est consacré à l'étude théorique et numérique de la dynamique d'une solution stationnaire (pour l'équation déterministe) de type vortex. Nous étudions l'influence des perturbations aléatoires du potentiel sur la solution, et montrons que la solution perturbée garde les symétries de la solution stationnaire pour des temps au moins de l'ordre du carré de l'inverse de l'intensité des fluctuations. Ces résultats sont illustrés par des simulations numériques exploitant une méthode de Monte-Carlo adaptée à la simulation d'événements rares.Le deuxième modèle permet de modéliser les effets de la température sur la dynamique d'un condensat. Lorsque celle-ci n'est pas nulle, la condensation n'est pas complète et le condensat interagit avec les particules non condensées. Ces interactions sont d'un grand intérêt pour comprendre la dynamique de transition de phase et analyser les phénomènes de brisure de symétrie associés, comme la formation spontanée de vortex. Nous nous sommes intéressés dans les chapitres 3 et 4 à des questions relatives à la simulation de la distribution des solutions de cette équation en temps long. Le troisième chapitre est consacré à la construction d'une méthode d’échantillonnage sans biais pour des mesures connues à une constante multiplicative près. C'est une méthode de Monte Carlo par chaînes de Markov qui a la particularité de permettre un échantillonnage non-réversible basé sur une équation de type Langevin sur-amortie. Elle constitue une extension de Metropolis-Adjusted Langevin Algorithm (MALA). Le quatrième chapitre est quant à lui consacré à l'étude numérique de dynamiques métastables liées à la nucléation de vortex dans des condensats en rotation. Un intégrateur numérique pour la dynamique étudiée est proposé, ainsi qu'une méthode de Monte-Carlo adaptée à la simulation d'événements rares correspondant aux changements de configurations métastables. Cette dernière est basée sur l'algorithme Adaptive Multilevel Splitting (AMS). / This thesis is devoted to the numerical study of two stochastic models arising in Bose-Einstein condensation physics. They constitute two generalisations of the Gross-Pitaevskii Equation. This deterministic partial differential equation model the wave function dynamics of a Bose-Einstein condensate trapped in an external confining potential. The first chapter contains a simple presentation of the Bose-Einstein condensation phenomenon and of the experimental methods used to construct such systems.The first model considered enables to model the fluctuations of the confining potential intensity, and takes the form of a stochastic partial differential equation. In practice, these fluctuations lead to heating of the condensate and possibly to its collapse. In the second chapter we propose to build a numerical scheme to solve this model. It is based on a spectral space discretisation and a Crank-Nicolson discretisation in space. We show that the proposed scheme converges strongly at order at least one in probability. We also present numerical simulations to illustrate this result. The third chapter is devoted to the numerical and theoretical study of the dynamics of a stationary solution (for the deterministic equation) of vortex type. We study the influence of random disturbances of the confining potential on the solution. We show that the disturbed solution conserves the symmetries of the stationary solution for times up to at least the square of the inverse of the fluctuations intensity. These results are illustrated with numerical simulations based on a Monte-Carlo method suited to rare events estimation.The second model can be used to model the effects of the temperature on the dynamics of a Bose-Einstein condensate. In the case of finite temperature, the Bose-Einstein condensation is not complete and the condensate interacts with the non-condensed particles. These interactions are interesting to understand the dynamics of the phase transition and analyse the phenomena of symmetry breaking associated, like the spontaneous nucleation of vortices We have studied in the fourth and the fifth chapters some questions linked to the long time simulation of this model solutions. The fourth chapter is devoted to the construction of an unbiased sampling method of measures known up to a multiplicative constant. The distinctive feature of this Markov-Chain Monte-Carlo algorithm is that it enables to perform an unbiased non-reversible sampling based on an overdamped Langevin equation. It constitutes a generalization of the Metropolis-Adjusted Langevin Algorithm (MALA). The fifth chapter is devoted to the numerical study of metastable dynamics linked to the nucleation of vortices in rotating Bose-Einstein condensates. A numerical integrator and a suited Monte-Carlo methods for the simulation of metastable dynamics are proposed. This Monte-Carlo method is based on the Adaptive Multilevel Splitting (AMS) algorithm.

Dans le présent travail de thèse, nous souhaitons approfondir l’étude des systèmes à commutation pour des problèmes aux dérivées partielles en explorant de nouvelles pistes d’investigation, incluant notamment la question de la synthèse de contrôle garanti par décomposition de l’espace des états, la synthèse de contrôle nécessitant la réduction de modèle, le contrôle des différentes sources d’erreur sur des quantités d’intérêt, et la mesure des incertitudes sur les états et les paramètres du modèle. Nous envisageons l’utilisation de méthodes de calcul ensemblistes associées à des méthodes de réduction de modèle, ainsi que l’utilisation d’observateurs d’état pour l’estimation en ligne du système. / In this thesis, we focus on switched control systems described by partial differential equations, and investigate the issues of guaranteed control of such systems using state-space decomposition methods. The use of state-space decomposition methods requires model order reduction, control of the different sources of error for quantities of interest, and measure of uncertainties on the states and parameters of the system. We are considering using set-based computation methods, in association with model order reduction techniques, along with the use of state-observers for on-line estimation of the system.

L’objet de cette thèse est l’étude des propriétés de stabilité et contrôle pour des systèmes linéaires écrits à l’aide d’équations aux dérivées partielles (EDP) ou d’équations à retard. Nous souhaitons exploiter dans cette thèse les liens qui existent entre ces deux classes de systèmes de dimension infinie afin de développer une nouvelle approche permettant leur analyse. En effet dans plusieurs applications, il est possible de choisir l’un ou l’autre de ces deux types de systèmes pour modéliser la dynamique considérée. Par exemple, les phénomènes de congestion dans un réseau routier peuvent être modélisés à l’aide d’EDP de type transport [JKC], mais aussi par un modèle à retard distribué [MMN] ou encore à retard discret [SN]. On peut également renvoyer aux travaux de Krstic [K] sur la formulation d’un système à retard comme un système EDP. Ces deux classes de systèmes sont des cas particuliers de systèmes de dimension infinie, et contrairement aux cas de systèmes de dimension finie, on parle de fonctions d’état plutôt que vecteur d’état. Cela implique que l’analyse associée est plus délicate et fait appel à des outils dédiés. Dans le cadre de la thèse, l’étudiant se focalisera sur les approches basées sur une extension du théorème de Lyapunov pour les systèmes de dimension infinie utilisant des fonctionnelles spécifiques. Comme pour la modélisation, l’analyse de stabilité des systèmes à retard ou de type EDP peut être menée à l’aide de fonctionnelles de Lyapunov très similaires. Nous souhaitons que cette thèse tire parti des travaux existants dans les deux communautés sur les systèmes à retards et de type EDP pour développer une approche novatrice et unifiée pour l’analyse et le contrôle de systèmes de dimension infinie. Pour cela, le candidat s’appuiera sur ses acquis en automatique et en mathématiques ainsi que sur l’expertise des deux encadrants. Plusieurs contributions sont attendues durant la thèse. Dans un premier temps, il sera question d’étendre des résultats récents [SG1,2] développés pour l’analyse de stabilité des systèmes à retards au cas de systèmes régis par des EDP. Ces premiers résultats auront vocation à servir de base pour l’étude de la synthèse de commandes robustes dans le cadre d’applications telles que le contrôle de trafic routier [MMN], le contrôle de vibration [RBPA], etc… Cette thèse en automatique requiert plusieurs compétences parmi lesquelles des connaissances sur la théorie de Lyapunov pour les systèmes avec ou sans retard, sur les inégalités matricielles linéaires tout en s’appuyant sur les outils de mathématiques appliquées pour l’étude des équations aux dérivées partielles (algèbre linéaire, analyse fonctionnelle, espaces de Hilbert, de Sobolev). / The purpose of this thesis is the study of stability and control properties for linear systems described by partial differential equations (PDE) or delay differential equations. We wish to use in this thesis the relationship between these two classes of infinite-dimensional systems in view of developing a new paradigm for their analysis. Indeed, in many applications, it is possible to choose one or the other of these two classes of systems to model the dynamics of the system under consideration. For example, traffic flow can be modeled using PDE type of transportation [JKC], but also by a distributed delay model [SMP] or discrete delay [SN]. We may also refer to the work of Krstic [K] on the formulation of a delay system as an PDE system. These two classes of systems are special cases of infinite dimensional systems, unlike the case of finite-dimensional systems, we better called state functions rather than the state vector. This implies that the analysis is more delicate and refers to the use of dedicated tools. As part of the thesis, the student will focus on approaches based on an extension of Lyapunov theorem for infinite dimensional systems using specific functional. As for the modeling process, the stability analysis of delayed or PDE type systems can be conducted using very similar Lyapunov functionals. We hope that this thesis builds on existing work in the two communities on delay systems and PDE to develop an innovative and unified approach to the analysis and control of infinite dimensional systems. To do so, the candidate will build on its skills in automatic and mathematics as well as the on from expertise of both supervisors. Several contributions are expected during the thesis . Initially, we aim at extedning recent results [SG13,14] developed in the context of the stability analysis of delay systems to the case of systems governed by PDE. These first results will provide the basis for the design of robust control laws for various applications including traffic control, vibration control, etc ... Cette thèse portera sur l’étude des propriétés de stabilité et de contrôle des systèmes linéaires de dimension infinie, plus particulièrement écrits à l’aide d’EDP ou d’équations à retard. L’intérêt naturel pour l’étude de cette classe de systèmes à la frontière entre mathématiques appliquées et automatique connaît un succès grandissant de part la large gamme d’applications en contrôle pouvant être décrites par ces modèles : en ingénierie, biologie, informatique… L’émulation scientifique entre systèmes à retard et systèmes de type EDP permettra en outre à cette thèse de tirer parti des méthodes et outils propres à chacun des ces domaines. This PhD proposal in automatic control requires several skills including knowledge on Lyapunov theory for systems with or without delay , on linear matrix inequalities while relying on mathematical tools applied in the study of partial differential equations ( linear algebra functional analysis , Hilbert spaces , Sobolev) .

Les guides d'ondes optiques sont utilisés dans plusieurs dispositifs de télécommunications pour le transport d'informations numériques. Nous les retrouvons aussi dans les équipements médicaux ou d'ingénierie, comme par exemple les coupleurs, les interféromètres optiques, les convertisseurs de polarisation, etc. Le transport des informations dans ces structures demeure une question capitale et donc les méthodes utilisées doivent être capables de véhiculer ces informations avec un minimum de perte. Dans cette thèse, nous présentons une méthode numérique d'approximation qui permet la simulation du transport des ondes dans une direction choisie avec le maximum de précision et à un coût abordable. Cette méthode est basée sur une approximation au sens de Padé d'opérateurs pseudo-différentiels associés à la fonction algébrique y/1 — x . Essentiellement, nous approchons une équation différentielle ordinaire (EDO) d'ordre 2 par une EDO d'ordre 1 dont les solutions ont la même phase que celle des solutions de l'EDO d'ordre 2 qui se propagent dans la direction choisie. Nous avons évalué les capacités de cette méthode sur de nombreux problèmes bidimensionnels et tridimensionnels. Elle s'est révélée à la fois robuste et précise

Les travaux de cette thèse s'inscrivent dans le domaine des équations aux dérivées partielles et sont liés à la problématique des changements d'échelle dans le contexte de la cinétique des gaz. En effet, sachant qu'il existe plusieurs niveaux de description pour un gaz, on cherche à relier les différentes échelles associées dans un cadre où une part d'aléa intervient. Dans une première partie, on établit la dérivation rigoureuse de l'équation de Boltzmann linéaire sans cut-off en partant d'un système de particules interagissant via un potentiel à portée infinie en partant d'un équilibre perturbé.La deuxième partie traite du passage d'un modèle BGK stochastique avec champ fort à une loi de conservation scalaire avec forçage stochastique. D'abord, on établit l'existence d'une solution au modèle BGK considéré. Sous une hypothèse additionnelle, on prouve alors la convergence vers une formulation cinétique associée à la loi de conservation avec forçage stochastique.Au cours de la troisième partie, on quantifie dans le cas à vitesses discrètes le défaut de régularité dans les lemmes de moyenne et on établit un lemme de moyenne stochastique dans ce même cas. On applique alors le résultat au cadre de l'approximation de Rosseland pour établir la limite diffusive associée à ce modèle.Enfin, on s'intéresse à l'étude numérique du modèle de Uchiyama de particules carrées à quatre vitesses en dimension deux. Après avoir adapté les méthodes de simulation développées dans le cas des sphères dures, on effectue une étude statistique des limites à différentes échelles de ce modèle. On rejette alors l'hypothèse d'un mouvement Brownien fractionnaire comme limite diffusive / The work of this thesis belongs to the field of partial differential equations and is linked to the problematic of scale changes in the context of kinetic of gas. Indeed, knowing that there exists different scales of description for a gas, we want to link these different associated scales in a context where some randomness acts, in initial data and/or distributed on all the time interval. In a first part, we establish the rigorous derivation of the linear Boltzmann equation without cut-off starting from a particle system interacting via a potential of infinite range starting from a perturbed equilibrium. The second part deals with the passage from a stochastic BGK model with high-field scaling to a scalar conservation law with stochastic forcing. First, we establish the existence of a solution to the considered BGK model. Under an additional assumption, we prove then the convergence to a kinetic formulation associated to the conservation law with stochastic forcing. In the third part, first we quantify in the case of discrete velocities the defect of regularity in the averaging lemmas. Then, we establish a stochastic averaging lemma in that same case. We apply then the result to the context of Rosseland approximation to establish the diffusive limit associated to this model.Finally, we are interested into the numerical study of Uchiyama's model of square particles with four velocities in dimension two. After adapting the methods of simulation which were developed in the case of hard spheres, we carry out a statistical study of the limits at different scales of this model. We reject the hypothesis of a fractional Brownian motion as diffusive limit

Le sujet principal de cette thèse porte sur l'approximation polynômiale des fonctions aléatoires au moyen de la projection L2 aléatoire discrète, et son application aux problèmes inverses pour les équations aux dérivées partielles avec des données aléatoires. Les motivations proviennent de l'approximation paramétrique de la solution de modèles aux dérivées partielles. La thèse se compose de deux parties, avec un chapitre d'introduction qui résume les techniques modernes de l'approximation polynômiale des fonctions de variables aléatoires. La première partie, du chapitre 1 au chapitre 4, contient l'analyse théorique de la projection L2 aléatoire discrète pour résoudre le problème direct, par exemple, pour rapprocher les moments d'une fonction aléatoire à partir de ses observations, ou pour calculer la solution à un modèle numérique avec des coefficients stochastiques. La stabilité et l'optimalité de l'erreur d'approximation évaluée dans la norme L2 pondérée sont traités. Dans la dernière partie de la thèse, composé des chapitres 5 et 6, la méthodologie développée précédemment pour le problème direct est appliqué aux problèmes inverses pour les équations aux dérivées partielles à coefficients stochastiques. La méthode de factorisation est appliquée dans le cadre de la tomographie par impédance électrique, d'abord dans le cas de coefficient inhomogène, puis dans le cas de coefficient constante par morceaux, à valeurs dans chaque région affectée par l'incertitude. Enfin, dans le chapitre 6 les variantes de la méthode de factorisation proposées dans le chapitre précédent sont accélérés en utilisant les techniques qui ont été présentés dans la première partie de la thèse.

La modélisation du trafic routier dans un réseau dense et de grande étendue nécessite un grand nombre de données, si bien qu'une modélisation par arcs est impossible en pratique. Pour simplifier le problème, une idée est d'agréger les tronçons du réseau en un continuum sur lequel le trafic routier s'écoule comme un fluide surfacique. Cette modélisation est qualifiée de bidimensionnelle. Même si la structure géométrique détaillée du réseau est perdue, une telle modélisation évite la description très fine du trafic sur un réseau dans lequel les points de mesure ne sont pas en nombre suffisant pour permettre une évaluation exhaustive de l'état du trafic. Une série de travaux commencée dans les années 1980 a permis de dégager quelques concepts importants. Cependant, ces travaux n'ont pas résolu les problèmes de modélisation fondamentaux : comment déduire et modéliser des comportements globaux à partir de comportements locaux (flux sur un axe, interactions aux intersections) ? Deux modèles bidimensionnels de trafic sont développés. Le premier modèle est statique. Le trafic s'écoule dans des directions de propagation privilégiées (orthogonales). Le modèle prend en compte l'équilibre entre l'offre de déplacement du réseau et une demande élastique de déplacement des usagers. Les principales sorties sont constituées en chaque point et pour chaque destination par les débits directionnels et les coûts de déplacement. Le deuxième modèle est dynamique. L'écoulement du trafic est décrit au niveau de cellules élémentaires du réseau dans lesquelles on définit les notions d'offre et de demande. À partir d'une loi comportementale obtenue sur un réseau discret, on écrit l'équation dynamique de conservation du trafic routier en tout point d'un réseau anisotrope