Cette thèse est divisée en trois chapitres. Le premier explique comment utiliser la méthode «level-set» de manière rigoureuse pour faire la simulation de feux de forêt en utilisant comme modèle physique pour la propagation le modèle de l'ellipse de Richards. Le second présente un nouveau schéma semi-implicite avec une preuve de convergence pour la solution d'une équation de type Hamilton-Jacobi anisotrope. L'avantage principal de cette méthode est qu'elle permet de réutiliser des solutions à des problèmes «proches» pour accélérer le calcul. Une autre application de ce schéma est l'homogénéisation. Le troisième chapitre montre comment utiliser les méthodes numériques des deux premiers chapitres pour étudier l'influence de variations à petites échelles dans la vitesse du vent sur la propagation d'un feu de forêt à l'aide de la théorie de l'homogénéisation. / This thesis is divided in three chapters. The first explains how to use the level-set method in a rigorous way in the context of forest fire simulation when the physical propagation model for firespread is Richards' ellipse model. The second chapter presents a new semi-implicit scheme with a proof of convergence for the numerical solution of an anisotropic Hamilton-Jacobi partial differential equation. The advantage of this scheme is it allows the use of approximative solutions as initial conditions which reduces the computation time. The third chapter shows how to use the tools introduced in the first two chapters to study the influence of small-scale variations on the wind speed on firespread using the theory of homogenization.

Cette thèse a pour objet la simulation d'une zone de mélange turbulente issue de l'instabilité de Richtmyer-Meshkov à l'aide d'un modèle à fonction de densité de probabilité (PDF). Nous analysons plus particulièrement la prise en charge par le modèle PDF du transport de l'énergie cinétique turbulente dans la zone de mélange.Dans cette optique, nous commençons par mettre en avant le lien existant entre les statistiques en un point de l'écoulement et ses conditions initiales aux grandes échelles. Ce lien s'exprime à travers le principe de permanence des grandes échelles, et permet d'établir des prédictions pour certaines grandeurs de la zone de mélange, telles que son taux de croissance ou son anisotropie.Nous dérivons ensuite un modèle PDF de Langevin capable de restituer cette dépendance aux conditions initiales. Ce modèle est ensuite validé en le comparant à des résultats issus de simulations aux grandes échelles (LES).Enfin, une analyse asymptotique du modèle proposé permet d'éclairer notre compréhension du transport turbulent. Un régime de diffusion est mis en évidence, et l'expression du coefficient de diffusion associé à ce régime atteste l'influence de la permanence des grandes échelles sur le transport turbulent.Tout au long de cette thèse, nous nous sommes appuyés sur des résultats issus de simulations de Monte Carlo du modèle de Langevin. A cet effet, nous avons développé une méthode spécifique eulérienne et à l'avons comparé à des alternatives lagrangiennes. / The aim of the thesis is to simulate a turbulent mixing zone resulting from the Richtmyer-Meshkov instability using a probability density function (PDF) model. An emphasis is put on the analysis of the turbulent kinetic energy transport.To this end, we first highlight the link existing between the one-point statistics of the flow and its initial conditions at large scales. This link is expressed through the principle of permanence of large eddies, and allows to establish predictions for quantities of the mixing zone, such as its growth rate or its anisotropy.We then derive a Langevin PDF model which is able to reproduce this dependency of the statistics on the initial conditions. This model is then validated by comparing it against large eddy simulations (LES).Finally, an asymptotic analysis of the derived model helps to improve our understanding of the turbulent transport. A diffusion regime is identified, and the expression of the diffusion coefficient associated with this regime confirms the influence of the permanence of large eddies on the turbulent transport.Throughout this thesis, our numerical results were based on Monte Carlo simulations for the Langevin model. In this regard, we proceeded to the development of a specific Eulerian method and its comparison with Lagrangian counterparts.

Cette thèse s'inscrit dans l'étude des modèles asymptotiques aux équations des ondes de surface dans le régime modulationnel. Le problème des ondes de surface consiste à décrire le mouvement - sous l'influence de la gravitation et éventuellement de tension de surface - d'un fluide dans un domaine délimité par la surface libre du fluide et par un fond fixe. Dans l'étude de ce problème, on s'intéresse en particulier aux ondes se propageant à la surface du fluide.Dans le régime modulationnel, on considère l'évolution des ondes de surface sous forme de paquets d'ondes de faible amplitude se propageant dans une direction. Il est bien connu que la motion de l'enveloppe du paquet d'onde sur une échelle de temps d'ordre t = O(1/ϵ²), où ϵ est un petit paramètre désignant l'amplitude, est décrite approximativement par des systèmes d'équations appelés systèmes de Benney-Roskes (BR) / Davey-Stewartson (DS). Ces systèmes sont donnés par une équation de type Schrödinger cubique couplée à une équation d'ondes. L'approximation classique de BR / DS est bien établie et a été largement étudiée au cours des dernières décennies. Récemment, David Lannes a introduit une version à "dispersion exacte" de ces systèmes. Contrairement aux équations de BR / DS standard, les systèmes à dispersion exacte préservent la relation de dispersion des équations des ondes de surface. On devrait obtenir ainsi une description plus riche du vrai comportement dynamique des ondes de surface que dans le cas de l'approximation classique.Le systèmes de BR / DS à dispersion exacte sont étudiés dans cette thèse. La première partie est consacrée à la déduction formelle des systèmes de BR / DS en tant que modèles asymptotiques aux équations des ondes de surface. Nous donnons en outre un résultat sur la consistance de cette approximation.Ensuite, nous étudions le problème de Cauchy pour le système de BR à dispersion exacte. En fait, afin de justifier la consistance de l'approximation de BR avec les équations exactes, on doit prouver que ce système est bien posé (en espace de Sobolev) sur une échelle de temps d'ordre O(1/ϵ). Ceci est un problème ouvert même dans le cas classique, du moins pour le système de dimension 1 + 2. De même, nous ne pouvons pas démontrer l'existence de solutions en temps long pour le système de BR à dispersion exacte, mais nous obtenons un théorème d'existence locale (t = O(1)) à condition que la tension de surface soit assez forte. Si nous nous restreignons au système de dimension 1+1, nous pouvons enlever la contrainte sur la tension de surface. L'idée de la preuve d'existence locale, qui est inspirée par un travail de Schochet-Weinstein, est d'écrire le système de BR comme un système symétrique hyperbolique quasi-linéaire perturbé par un terme dispersif ne contribuant pas à l'énergie du système. Ainsi, nous pouvons appliquer les méthodes standard de résolution des systèmes hyperboliques.En modifiant le terme non-linéaire du système de BR de dimension 1+1 sans changer l'ordre de consistance, nous obtenons un système qui est bien posé sur l'échelle de temps appropriée O(1/ϵ). Cependant, cette démarche ne peut pas être généralisée au cas de dimension 1+2.Dans le dernier chapitre de cette thèse, nous donnons quelques résultats sur les systèmes de Davey-Stewartson à dispersion exacte. Pour les systèmes de DS, il est suffisant de démontrer qu'ils sont bien posés localement afin de justifier leur consistance avec les équations des ondes de surface. La théorie d'existence de solutions est assez complète pour le système de DS classique. Dans le cas de dispersion exacte cependant, les équations paraissent mal posées généralement, si bien que l'existence locale ne peut être démontrée pour l'instant que pour quelques cas particuliers simples. / This thesis is concerned with asymptotic models to the water wave equations in the modulational regime. The water wave equations describe the motion - under the influence of gravity and possibly surface tension - of an inviscid fluid in a domain which is bounded by a fixed bottom from below and the free surface of the fluid from above. In the study of the water wave problem, one is in particular interested by waves propagating on the surface of the fluid.In the modulational regime, one considers the evolution of surface waves under the form of small amplitude wave packets traveling in one direction. It is well known that the evolution of the wave packet envelope on the long time scale t = O(1/ϵ²), where ϵ is a small parameter denoting the amplitude of the wave, is approximately governed by a set of equations known as the Benney-Roskes (BR) / Davey-Stewartson (DS) systems. These systems are essentially given by a cubic Schrödinger-type equation coupled to a wave equation. The classical BR / DS approximation is well established and has been largely studied in the past decades. Recently, David Lannes has introduced a "full dispersion" version of these systems. In contrast to the standard BR / DS equations, the full dispersion systems preserve the linear dispersion relation of the full water wave equations, and should therefore give a richer description of the original wave dynamics than the classical approximation.The full dispersion BR / DS systems are studied in this thesis. In the first part, we formally derive the full dispersion BR / DS approximation from the water wave equations both in the case of zero and positive surface tension. The formal derivation is completed by a consistency result.We then study well-posedness in Sobolev space of the full dispersion BR system. In order to justify consistency of the BR approximation with the full water wave equations, one needs to show that the BR system is well posed on a time scale of order O(1/ϵ). This is an open problem even in the classical case, at least for the 1 + 2 dimensional system. We also do not obtain well-posedness on the long time scale for the full dispersion BR system, but we can show that it is locally well-posed in the case of sufficiently strong surface tension, and additionally in the zero surface tension case if we restrict ourselves to the 1+1 dimensional system. The proof is inspired by a paper of Schochet-Weinstein, and is based on writing the full dispersion BR system as a quasilinear symmetric hyperbolic system with dispersive perturbation, where the dispersive terms do not contribute to the energy. We can therefore apply classical solution methods for hyperbolic systems.By modifying the nonlinear part of the 1+1 dimensional full dispersion BR system without changing consistency, we obtain a system that is well-posed on the appropriate O(1/ϵ) time scale. This approach however does not generalize to the 1+2 dimensional case.In the last chapter of the thesis, we give some results on the full dispersion DS systems, which are obtained as special limits of the full dispersion BR system. For these systems, it is sufficient to prove local well-posedness in order to show consistency with the water wave equations. For the standard DS systems, local well-posedness theory is quite complete. For the full dispersion systems, the analysis is complicated by some nonlocal operators and the equations seem to be generally ill-posed. There are however some simple cases where local well-posedness can be shown. We also discuss some modifications of the full dispersion DS system that might allow to solve it for a larger range of parameters.

Cette thèse s'inscrit dans le cadre général de l'étude de la dynamique de populations. Elle porte sur la modélisation et l'analyse mathématique de l'hématopoïèse, le processus de production et de régulation des cellules sanguines. La population de cellules est perçue comme un milieu continu avec une structuration en âge et en espace. Nous avons commencé par analyser des modèles d'équations différentielles et aux différences à retard discret et distribué. Ces modèles à retard permettent de mettre en évidence des comportements particuliers tels que l'existence de solutions périodiques. Ensuite, nous avons pris en compte l'aspect spatial et la diffusion des cellules dans ces modèles, tout en sachant que la structuration en espace, dans le cas de l'hématopoïèse, a été très peu abordée par le passé. Un nouveau modèle a été obtenu du point de vue mathématique. Une étude d'existence d'ondes progressives est effectuée lorsque le domaine est non borné et lorsque le domaine est borné une étude de stabilité des états stationnaires ainsi que de l'existence d'une bifurcation de Hopf est réalisée / This thesis focuses on the study of population dynamics. It is devoted to the mathematical analysis and modeling of hematopoiesis, which is the process leading to the production and regulation of blood cells. The cell's population is seen as a continuous medium structured in age and space. We analyzed models of differential-difference system with discrete- and distributed -delay. These models can exhibit specific behaviors such as the existence of periodic solutions. Then we consider a space structuration and the diffusion of cells in such models, knowing that the space structure has not been widely studied in the case of hematopoiesis. A new model is obtained from the mathematical point of view. We studied the existence of traveling waves when the domain is unbounded. When the domain is bounded, the stability of stationary solutions and the existence of a Hopf bifurcation are obtained

Lors de l’infection par un pathogène intracellulaire, l’organisme déclenche une réponse immunitaire spécifique dont les acteurs principaux sont les lymphocytes T-CD8. Ces cellules sont responsables de l’éradication de ce type d’infections et de la constitution du répertoire immunitaire de l’individu. Les processus qui composent la réponse immunitaire se répartissent sur plusieurs échelles physiques inter-connectées (échelle intracellulaire, échelle d’une cellule, échelle de la population de cellules). La réponse immunitaire est donc un processus complexe, pour lequel il est difficile d’observer ou de mesurer les liens entre les différents phénomènes mis en jeu. Nous proposons trois modèles mathématiques multi-échelles de la réponse immunitaire, construits avec des formalismes différents mais liés par une même idée : faire dépendre le comportement des cellules TCD8 de leur contenu intracellulaire. Pour chaque modèle, nous présentons, si possible, sa construction à partir des hypothèses biologiques sélectionnées, son étude mathématique et la capacité du modèle à reproduire la réponse immunitaire au travers de simulations numériques. Les modèles que nous proposons reproduisent qualitativement et quantitativement la réponse immunitaire T-CD8 et constituent ainsi de bons outils préliminaires pour la compréhension de ce phénomène biologique. / Upon infection by an intracellular pathogen, the organism triggers a specific immune response,mainly driven by the CD8 T cells. These cells are responsible for the eradication of this type of infections and the constitution of the immune repertoire of the individual. The immune response is constituted by many processes which act over several interconnected physical scales (intracellular scale, single cell scale, cell population scale). This biological phenomenon is therefore a complex process, for which it is difficult to observe or measure the links between the different processes involved. We propose three multiscale mathematical models of the CD8 immune response, built with different formalisms but related by the same idea : to make the behavior of the CD8 T cells depend on their intracellular content. For each model, we present, if possible, its construction process based on selected biological hypothesis, its mathematical study and its ability to reproduce the immune response using numerical simulations. The models we propose succesfully reproduce qualitatively and quantitatively the CD8 immune response and thus constitute useful tools to further investigate this biological phenomenon.

Cette thèse est consacrée à l’étude des équations aux dérivées partielles stochastiques de type parabolique. Dans la première partie nous démontrons de nouveaux résultats concernant l’existence et l’unicité de solutions variationnelles globales et locales à des problèmes avec des conditions aux bords de type Neumann pour une classe d’équations aux dérivées partielles stochastiques non-autonomes. Les équations que nous considérons sont définies sur des domaines non bornés de l’espace euclidien qui satisfont à certaines conditions géométriques, et sont dirigées par un bruit multiplicatif dérivé d’un processus de Wiener fractionnaire infini-dimensionnel caractérisé par une suite de paramètres de Hurst H = (Hi) i ∈ N+ ⊂ (1/2,1). Ces paramètres sont en fait soumis à d’autres contraintes intimement liées à la nature de la non-linéarité dans le terme stochastique des équations, et au choix des espaces fonctionnels dans lesquels le problème à résoudre est bien posé. Notre méthode de preuve repose essentiellement sur des arguments d’injections compactes. Dans la seconde partie, nous étudions la possibilité de l’explosion de solutions d’une classe d’équations aux dérivées partielles stochastiques semi-linéaire avec des conditions aux bords de type Dirichlet, perturbées par un mélange d’un mouvement brownien et d’un mouvement brownien fractionnaire et dirigées par une classe d’opérateurs différentiels non autonomes contenant des processus de diffusions et des processus de Lévy. Notre but est de comprendre l’influence de la partie stochastique et de l’opérateur différentiel sur le comportement d’explosion des solutions. En particulier, nous donnons des expressions explicites pour des bornes inférieures et supérieures du temps de l’explosion de la solution, et des conditions suffisantes pour l’existence d’une solution globale positive. Nous estimons également la probabilité d’une explosion en temps fini et la loi d’une borne supérieur du temps d’explosion de la solution / This thesis is concerned with stochastic partial differential equations of parabolic type. In the first part we prove new results regarding the existence and the uniqueness of global and local variational solutions to a Neumann initial-boundary value problem for a class of non-autonomous stochastic parabolic partial differential equations. The equations we consider are defined on unbounded open domains in Euclidean space satisfying certain geometric conditions, and are driven by a multiplicative noise derived from an infinite-dimensional fractional Wiener process characterized by a sequence of Hurst parameters H = (Hi) i ∈ N+ ⊂ (1/2,1). These parameters are in fact subject to further constraints that are intimately tied up with the nature of the nonlinearity in the stochastic term of the equations, and with the choice of the functional spaces in which the problem at hand is well-posed. Our method of proof rests on compactness arguments in an essential way. The second part is devoted to the study of the blowup behavior of solutions to semilinear stochastic partial differential equations with Dirichlet boundary conditions driven by a class of differential operators including (not necessarily symmetric) Lévy processes and diffusion processes, and perturbed by a mixture of Brownian and fractional Brownian motions. Our aim is to understand the influence of the stochastic part and that of the differential operator on the blowup behavior of the solutions. In particular we derive explicit expressions for an upper and a lower bound of the blowup time of the solution and provide a sufficient condition for the existence of global positive solutions. Furthermore, we give estimates of the probability of finite time blowup and for the tail probabilities of an upper bound for the blowup time of the solutions

La thèse est consacrée au développement de nouvelles méthodes de modélisations mathématiques en biologie et en médecine, du type “off-lattice" modèles hybrides discret-continus, et de leurs applications à l’hématopoïèse et aux maladies sanguines telles la leucémie et l’anémie. Dans cette approche, les cellules biologiques sont considérées comme des objets discrets alors que les réseaux intracellulaire et extracellulaire sont décrits avec des modèles continus régis par des équations aux dérivées partielles et des équations différentielles ordinaires. Les cellules interagissent mécaniquement et biochimiquement entre elles et avec le milieu environnant. Elles peuvent se diviser, mourir par apoptose ou se différencier. Le comportement des cellules est déterminé par le réseau de régulation intracellulaire et influencé par le contrôle local des cellules voisines ou par la régulation globale d’autres organes. Dans la première partie de la thèse, les modèles hybrides du type “off-lattice" dynamiques sont introduits. Des exemples de modèles, spécifiques aux processus biologiques, qui décrivent au sein de chaque cellule la concurrence entre la prolifération et l’apoptose, la prolifération et la différenciation et entre le cycle cellulaire et de l’état de repos sont étudiés. L’émergence des structures biologiques est étudiée avec les modèles hybrides. L’application à la modélisation des filamente de bactéries est illustrée. Dans le chapitre suivant, les modèle hybrides sont appliqués afin de modéliser l’érythropoïèse ou production de globules rouges dans la moelle osseuse. Le modèle inclut des cellules sanguines immatures appelées progéniteurs érythroïdes, qui peuvent s’auto-renouveler, se différencier ou mourir par apoptose, des cellules plus matures appelées les réticulocytes, qui influent les progéniteurs érythroïdes par le facteur de croissance Fas-ligand, et des macrophages, qui sont présents dans les îlots érythroblastiques in vivo. Les régulations intracellulaire et extracellulaire par les protéines et les facteurs de croissance sont précisées et les rétrocontrôles par les hormones érythropoïétine et glucocorticoïdes sont pris en compte. Le rôle des macrophages pour stabiliser les îlots érythroblastiques est montré. La comparaison des résultats de modélisation avec les expériences sur l’anémie chez les souris est effectuée. Le quatrième chapitre est consacré à la modélisation et au traitement de la leucémie. L’érythroleucémie, un sous-type de leucémie myéloblastique aigüe (LAM), se développe à cause de la différenciation insuffisante des progéniteurs érythroïdes et de leur auto-renouvellement excessif. Un modèle de type “Physiologically Based Pharmacokinetics-Pharmacodynamic” du traitement de la leucémie par AraC et un modèle de traitement chronothérapeutique de la leucémie sont examinés. La comparaison avec les données cliniques sur le nombre de blast dans le sang est effectuée. Le dernier chapitre traite du passage d’un modèle hybride à un modèle continu dans le cas 1D. Un théorème de convergence est prouvé. Les simulations numériques confirment un bon accord entre ces deux approches. / This dissertation is devoted to the development of new methods of mathematical modeling in biology and medicine, off-lattice discrete-continuous hybrid models, and their applications to modelling of hematopoiesis and blood disorders, such as leukemia and anemia. In this approach, biological cells are considered as discrete objects while intracellular and extracellular networks are described with continuous models, ordinary or partial differential equations. Cells interact mechanically and biochemically between each other and with the surrounding medium. They can divide, die by apoptosis or differentiate. Their fate is determined by intracellular regulation and influenced by local control from the surrounding cells or by global regulation from other organs. In the first part of the thesis, hybrid models with off-lattice cell dynamics are introduced. Model examples specific for biological processes and describing competition between cell proliferation and apoptosis, proliferation and differentiation and between cell cycling and quiescent state are investigated. Biological pattern formation with hybrid models is discussed. Application to bacteria filament is illustrated. In the next chapter, hybrid model are applied in order to model erythropoiesis, red blood cell production in the bone marrow. The model includes immature blood cells, erythroid progenitors, which can self-renew, differentiate or die by apoptosis, more mature cells, reticulocytes, which influence erythroid progenitors by means of growth factor Fas-ligand, and macrophages, which are present in erythroblastic islands in vivo. Intracellular and extracellular regulation by proteins and growth factors are specified and the feedback by the hormones erythropoietin and glucocorticoids is taken into account. The role of macrophages to stabilize erythroblastic islands is shown. Comparison of modelling with experiments on anemia in mice is carried out. The following chapter is devoted to leukemia modelling and treatment. Erythroleukemia, a subtype of Acute Myeloblastic Leukemia (AML), develops due to insufficient differentiation of erythroid progenitors and their excessive slef-renewal. A Physiologically Based Pharmacokinetics-Pharmacodynamics (PBPKPD) model of leukemia treatment with AraC drug and chronotherapeutic treatments of leukemia are examined. Comparison with clinical data on blast count in blood is carried out. The last chapter deals with the passage from a hybrid model to a continuous model in the 1D case. A convergence theorem is proved. Numerical simulations confirm a good agreement between these approaches.

Les problèmes de transmission avec couche mince sont délicats à approcher numériquement, en raison de la nécessité de construire des maillages à l’échelle de la couche mince. Il est courant d’éviter ces difficultés en usant de problèmes avec conditions aux limites approchées — dites d’impédance. Si l’approximation des problèmes de transmission par des problèmes d’impédance s’avère performante dans le cas de domaines réguliers, elle l’est beaucoup moins lorsque ceux-ci comportent des coins ou arêtes. L’objet de cette thèse est de proposer de nouvelles conditions d’impédance, plus performantes, afin de corriger cette perte de performance. Pour cela, les développements asymptotiques des différents problèmes-modèles sont construits et étudiés afin de localiser avec précision l’origine de la perte, en lien avec les profils singuliers associés aux coins et arêtes. De nouvelles conditions d’impédance sont construites, de type Robin multi-échelle ou Venctel. D’abord étudiées en dimension 2, elles sont ensuite généralisées à certaines situations en dimension 3. Des simulations viennent confirmer l’efficience des méthodes théoriques. / Transmission problems with thin layer are delicate to approximate numerically, because of the necessity to build meshes on the scale of the thin layer. It is common to avoid these difficulties by using problems with approximate boundary conditions — also called impedance conditions. Whereas the approximation of transmission problems by impedance problems turns out to be successful in the case of smooth domains, the situation is less satisfactory in the presence of corners and edges. The goal of this thesis is to propose new impedance conditions, more efficient, to correct this lack of performance. For that purpose, the asymptotic expansions of the various models -problems are built and studied to locate exactly the origin of the loss, in connection with the singular profiles associated to corners and edges. New impedance conditions are built, of multi-scale Robin or Venctel types. At first studied in dimension 2, they are then generalized in certain situations in dimension 3. Simulations have been carried out to confirm the efficiency of the theoretical methods to some.

Ce travail s’inscrit dans un contexte de contrôle de la pollution des ressources en eau. On s’intéresse plus particulièrement à l’impact des engrais d’origine agricole sur la qualité de l’eau, en alliant modélisation économique et hydrogéologique. Pour cela, nous définissons d’une part un objectif économique spatio-temporel prenant en compte le compromis entre l’utilisation d’engrais et les coûts de dépollution. D’autre part, nous décrivons le transport du polluant dans le sous-sol (3D en espace) par un système non linéaire d’équations aux dérivées partielles couplées de type parabolique (réaction-convection-dispersion) et elliptique dans un domaine borné. Nous prouvons l’existence globale d’une solution au problème de contrôle optimal. L’unicité est quant à elle démontrée par analyse asymptotique pour le problème effectif tenant compte de la faible concentration d’engrais en sous-sol. Nous établissons les conditions nécessaires d’optimalité et le problème adjoint associé à notre modèle. Quelques exemples analytiques sont donnés et illustrés. Nous élargissons ces résultats au cadre de la théorie des jeux, où plusieurs joueurs interviennent, et prouvons notamment l’existence d’un équilibre de Nash. Enfin, ce travail est illustré par des résultats numériques (2D en espace), obtenus en couplant un schéma de type Éléments Finis Mixtes avec un algorithme de gradient conjugué non linéaire. / This work is devoted to water ressources pollution control. We especially focus on the impact of agricultural fertilizer on water quality, by combining economical and hydrogeological modeling. We define, on one hand, the spatio-temporal objective, taking into account the trade off between fertilizer use and the cleaning costs. On an other hand, we describe the pollutant transport in the underground (3D in space) by a nonlinear system coupling a parabolic partial differential equation (reaction-advection-dispersion) with an elliptic one in a bounded domain. We prove the global existence of the solution of the optimal control problem. The uniqueness is proved by asymptotic analysis for the effective problem taking into account the low concentration fertilizer. We define the optimal necessary conditions and the adjoint problem associated to the model. Some analytical results are provided and illustrated. We extend these results within the framework of game theory, where several players are involved, and we prove the existence of a Nash equilibrium. Finally, this work is illustrated by numerical results (2D in space), produced by coupling a Mixed Finite Element scheme with a nonlinear conjugate gradient algorithm.

L'objet de cette thèse est l'étude, sur les variétés riemanniennes compactes $(V_n,g)$ de dimension $n>4$, de l'équation aux dérivées partielles elliptique de quatrième ordre $$(E)\; \Delta^2u+\nabla [a(x)\nabla u] +h(x)u= f(x)|u|^(N-2)u$$ où $a$, $h$, $f$ sont fonction $C^\infty $, avec $f(x)$ fonction constante ou partout positive et $N=(2n\over((n-4)))$ est l'exposant critique. En utilisant la méthode variationnelle on prouve dans le théorème principal que l'équation $(E)$ admet une solution $C^((5,\alpha))(V)$ $0<\alpha<1$ non nulle si une certaine condition qui dépend de la meilleure constante dans les inclusion de Sobolev ($H_2\subset L_(2n\over(n-4))$) est satisfaite. De plus on montre que si $a$ et $h$ sont des fonctions constantes bien précisées la solution de l'équation est positive et $C^\infty(V)$. Lorsque $n\geq 6$, on donne aussi des applications du théorème principal. Dans la dernière partie de cette thèse sur une variété riemannienne compacte à bord de dimension $n$, $(\overline(W)_n,g )$ nous nous intéressons au problème : $$ (P_N) \; \left\lbrace \begin(array)(c) \Delta^2 v+\nabla [a(x)\nabla u] +h(x) v= f(x)|v |^(N-2)v \; \hbox(sur)\; W \\ \Delta v =\delta \, , \, v = \eta \;\hbox(sur) \;\partial W \end(array)\right.$$ avec $\delta$,$\eta$,$f$ fonctions $C^\infty (\overline (W))$ avec $f(x)$ fonction partout positive et on démontre l'existence d'une solution non triviale pour le problème $(P_N)$.