Cette thèse est consacrée à l'étude théorique de quelques modèles d'évolution non-linéaires issus de la mécanique des fluides. Nous distinguons trois parties indépendantes. La première partie de la thèse traite essentiellement de l'existence des poches de tourbillon en rotation uniforme (appelées aussi V-states) pour un modèle quasi-géostrophique non visqueux. Notre étude est répartie sur deux chapitres où les poches présentent des structures topologiques différentes. Dans le premier chapitre nous étudions le cas simplement connexe et nous validons l'existence de ces structures dans un voisinage du tourbillon de Rankine en utilisant des techniques de bifurcation. Dans le deuxième chapitre nous abordons le cas doublement connexe où la poche admet un seul trou. Plus précisément, proche d'un anneau donné, nous décrivons cette famille par des branches dénombrables bifurquant de cet anneau à certaines valeurs explicites des vitesses angulaires liées aux fonctions de Bessel. Notre étude théorique a été complétée par des simulations numériques portant sur les V-states limites et un bon nombre de constatations ont été formulées ouvrant la porte à de nouvelles perspectives de recherche. La seconde partie concerne l'étude du problème de Cauchy pour le système de Boussinesq non visqueux 2D avec des données initiales de type Yudovich. Le problème est dans un certain sens critique à cause de quelques termes comportant la transformée de Riesz dans la formulation tourbillon-densité. Nous donnons une réponse positive pour une sous-classe comprenant les poches de tourbillon régulières et singulières. Dans la dernière partie nous analysons le problème de la limite incompressible pour les équations d'Euler isentropiques 2D associées à des données initiales très mal préparées et pour lesquelles les tourbillons ne sont pas forcément bornés mais appartiennent plutôt à des espaces de type ''BMO'' à poids. On utilise principalement deux ingrédients: d'un côté les estimations de Strichartz pour contrôler la partie acoustique. D'un autre côté, on se sert de la structure de transport compressible du tourbillon et on démontre une estimation de propagation linéaire dans l'esprit d'un travail récent de Bernicot et Keraani mené dans le cas incompressible. / In this dissertation, we are concerned with the study of some non-linear evolution models arising in fluid mechanics. We distinguish three independent parts. The first part of the thesis deals with the existence of the rotating vortex patches (called also V-states) for an inviscid quasi-geostrophic model. Our study is divided into two chapters dealing with different topological structures of the V-states. In the first chapter we study the simply connected case and we prove the existence of such structures in a neighborhood of the Rankine vortices by using the bifurcation theory. In the second chapter we discuss the doubly connected case where the patches admit only one hole. More precisely, close to a given annulus we describe this family by countable branches bifurcating from this annulus at some explicit angular velocities related to Bessel functions of the first kind. Our theoretical study was completed by numerical simulations on the limiting V-states and a number of interesting numerical observation were formulated opening new research perspectives. The second part of the thesis concerns the local well-posedness theory for the inviscid Boussinesq system with rough initial data. The problem is in some sense critical due to some terms involving Riesz transforms in the vorticity-density formulation. We give a positive answer for a special sub-class of Yudovich data including smooth and singular vortex patches. In the last part we address the problem of the incompressible limit for the 2D isentropic fluids associated to ill-prepared initial data and for which the vortices are not necessarily bounded and belong to some weighted BMO spaces. We mainly use two ingredients: On one hand, the Strichartz estimates to control the acoustic part and prove that it does not contribute for low Mach number. On the other hand, we use the transport compressible structure of the vorticity and we establish a linear propagation estimate in the spirit of a recent work of Bernicot and Keraani conducted in the incompressible case. The first part of the thesis deals with the existence of the rotating vortex patches (called also V-states) for an inviscid quasi-geostrophic model. Our study is divided into two chapters dealing with different topological structures of the V-states. In the first chapter we study the simply connected case and we prove the existence of such structures in a neighborhood of the Rankine vortices by using the bifurcation theory. In the second chapter we discuss the doubly connected case where the patches admit only one hole. More precisely, close to a given annulus we describe this family by countable branches bifurcating from this annulus at some explicit angular velocities related to Bessel functions of the first kind. Our theoretical study was completed by numerical simulations on the limiting V-states and a number of interesting numerical observation were formulated opening new research perspectives. The second part of the thesis concerns the local well-posedness theory for the inviscid Boussinesq system with rough initial data. The problem is in some sense critical due to some terms involving Riesz transforms in the vorticity-density formulation. We give a positive answer for a special sub-class of Yudovich data including smooth and singular vortex patches. In the last part we address the problem of the incompressible limit for the 2D isentropic fluids associated to ill-prepared initial data and for which the vortices are not necessarily bounded and belong to some weighted BMO spaces. We mainly use two ingredients: On one hand, the Strichartz estimates to control the acoustic part and prove that it does not contribute for low Mach number. On the other hand, we use the transport compressible structure of the vorticity and we establish a linear propagation estimate in the spirit of a recent work of Bernicot and Keraani conducted in the incompressible case.

Cette thèse s'intéresse à une étude variée des EDSRs. Une grande partie des résultats sont obtenus sous l'hypothèse d'une croissance de type quadratique du générateur en sa dernière variable. Un premier lien entre EDSRs quadratiques unidimensionnelles et théorie des jeux nous amène à développer des résultats avec générateurs convexes. La théorie du contrôle optimal nécessite quant à elle de traiter du cas multidimensionnel, dans lequel existence et unicité globales ne sont obtenues que pour des générateurs diagonalement quadratiques. Les résultats majeurs sur les EDSRs réfléchies (dont la solution est contrainte à rester dans un domaine) concernent des générateurs Lipschitziens. C'est dans ce cadre que nous développons un résultat de propagation du chaos, avec une contrainte portant sur la loi de la solution plutôt que sur sa trajectoire. Nous dressons enfin un pont entre EDSRs quadratiques et EDSRs réfléchies grâce aux EDSRs quadratiques de type champ moyen. Nous donnons plusieurs nouveaux résultats sur la possibilité de résoudre une équation quadratique dont le générateur dépend également de la moyenne des deux variables. / In this thesis, we are interested in studying variously Backward Stochastic Differential Equations. A large proportion of the results are obtained under the assumption that the driver is of quadratic growth in its last variable. A first link between one-dimensional quadratic BSDEs and game theory leads us to develop results with convex drivers. Optimal control theory requires as for it to deal with the multidimensional case, in which global existence and uniqueness are obtained only for diagonaly quadratic drivers. Major achievements in reflected BSDEs (whose solution is constrained to remain in a domain) are reached for Lipschitz drivers. We develop a result of chaos propagation in this setting, with a constraint on the law of the solution rather than on its path. We finaly build bridge between quadratic BSDEs and reflected BSDEs thanks to mean field quadratic BSDEs. We give several new results on solvability of a quadratic BSDE whose driver depends also on the mean of both variables.

Le but de cette thèse est d’étudier certaines équations aux dérivées partielles hyperboliques-dispersives. Une part importante est consacrée à l’analyse numérique et plus particulièrement à la convergence de schémas aux différences finies pour l’équation de Korteweg-de Vries et les systèmes abcd de Boussinesq. L’étude numérique suit les étapes classiques de consistance et de stabilité. Nous transposons au niveau discret la propriété de stabilité fort-faible des lois de conservations hyperboliques. Nous déterminons l’ordre de convergence des schémas et le quantifions en fonction de la régularité de Sobolev de la donnée initiale. Si nécessaire, nous régularisons la donnée initiale afin de toujours assurer les estimations de consistance. Une étape d’optimisation est alors nécessaire entre cette régularisation et l’ordre de convergence du schéma. Une seconde partie est consacrée à l’existence d’ondes progressives pour l’équation de Korteweg de Vries-Kuramoto-Sivashinsky. Par des méthodes classiques de systèmes dynamiques : système augmenté, fonction de Lyapunov, intégrale de Melnikov, par exemple, nous démontrons l’existence d’ondes oscillantes de petite amplitude. / The aim of this thesis is to study some hyperbolic-dispersive partial differential equations. A significant part is devoted to the numerical analysis and more precisely to the convergence of some finite difference schemes for the Korteweg-de Vries equation and abcd systems of Boussinesq. The numerical study follows the classical steps of consistency and stability. The main idea is to transpose at the discrete level the weak-strong stability property for hyperbolic conservation laws. We determine the convergence rate and we quantify it according to the Sobolev regularity of the initial datum. If necessary, we regularize the initial datum for the consistency estimates to be always valid. An optimization step is thus necessary between this regularization and the convergence rate of the scheme. A second part is devoted to the existence of traveling waves for the Korteweg-de Vries-Kuramoto-Sivashinsky equation. By classical methods of dynamical systems : extended systems, Lyapunov function, Melnikov integral, for instance, we prove the existence of oscillating small amplitude traveling waves.

Ce manuscrit de thèse porte sur l’analyse mathématique de modèles intégro-différentiels issus de la génétique des populations. Les deux modèles étudiés sont des équations de réaction-dispersion de type ∂tp(t,m) = UD[p](t,m) + f[p](t,m). Ils décrivent la dynamique de la distribution de la fitness (ou valeur sélective) dans une population asexuée sous l’effet des mutations et de la sélection représentées respectivement par les termes non locaux UD[p](t,m) et par f[p](t,m). La différence entre les deux modèles se situe au niveau du terme de mutation. En effet, dans le premier modèle, les effets des mutations sur la fitness ne dépendent pas de la fitness du parent, cela se traduit donc par un terme de convolution classique : D[p](t,m) =RR J(m−y)p(t,y)dy−p(t,m). Lorsqu’une mutation a lieu, la fonction J(m−y) représente la densité de probabilité pour un individu de fitness y d’avoir un descendant de fitness m. Le taux de mutation est donné par la constante U. Dans le second modèle, les effets des mutations sur la fitness dépendent aussi de la fitness du parent. Dans ce cas, un individu de fitness y a un descendant de fitness m avec la densité de probabilité Jy(m−y). Ce type de dépendance apparaît naturellement lorsque l’on suppose qu’il existe une fitness optimale (ou encore un optimum phénotypique). Pour chacun des deux modèles, nous établissons dans un premier temps des résultats d’existence et d’unicité ainsi que des propriétés de décroissance de la solution. Cette décroissance permet de définir la fonction génératrice des cumulants (CGF) associée à la distribution de fitness. La CGF est la solution d’une équation de transport non locale. Pour le premier modèle, l’étude de cette équation permet d’obtenir une solution analytique et donc d’obtenir une description complète de la distribution p(t,m) via ses moments. Nous étudions ensuite les états stationnaires pour chacun des deux modèles, et établissons des conditions suffisantes pour l’existence et la non-existence de phénomènes de concentration, correspondant à une accumulation d’individus de phénotypes optimaux. Nos résultats sont comparés à des sorties de modèles stochastiques individu-centrés représentant le même type de dynamiques évolutives. / This manuscript is devoted to the mathematical analysis of integro-differential models from population genetics. Both models are reaction-dispersion equations of the form ∂tp(t,m) = UD[p](t,m)+ f[p](t,m). They describe the dynamics of fitness distribution in an asexual population under the effect of mutation and selection. These two processes are represented by the nonlocal terms UD[p](t,m) and by f[p](t,m) respectively. The difference between the models rests on the mutation term. Indeed, in the first model, the mutation effects on fitness do not depend on the fitness of the parent. Thus, the mutation term is a standard convolution product: D[p](t,m) =RR J(m−y)p(t,y)dy −p(t,m). When a mutation occurs, the function J(m − y) represents the density of probability for an individual with fitness y to have an offspring with fitness m. The mutation rate is given by the constant U. In the second model, the mutation effects on fitness depend on the fitness of the parent. In this case, an individual with fitness y has an offspring with fitness m with a probability density Jy(m−y). This type of dependence naturally arises when the existence of an optimal fitness (or a phenotypic optimum) is assumed. For both models, we first establish existence and uniqueness results as well as decay properties of the solution. The decay property allows us to define the cumulant generating function (CGF). The CGF obeys a nonlocal transport equation. In the first model, we compute the analytical solution of this transport equation and thus, we obtain a complete description of the distribution p(t,m) through its moments. Then, we study the stationary states for both models, and establish sufficient conditions for the existence and non-existence of a concentration phenomenon corresponding to an accumulation of individuals with best possible phenotype. The results are compared to the results of stochastic individual based models which represent the same kind of evolutionary dynamics.

Cette thèse est divisée en deux parties reliées par un même fil conducteur. Elle porte sur l'étude théorique et l'application de modèles mathématiques décrivant des dynamiques de population où les individus se reproduisent et meurent à des taux dépendant de leur âge et d'un trait phénotypique. Le trait est fixé durant la vie de l'individu. Il est modifié au fil des générations par des mutations apparaissant lors de la reproduction. On modélise la sélection naturelle en introduisant un taux de mortalité densité-dépendant décrivant la compétition pour les ressources.Dans une première partie, nous nous intéressons au comportement en temps long d'une équation aux dérivées partielles de sélection-mutation structurée en âge décrivant une grande population d'individus. En étudiant les propriétés spectrales d'une famille d'opérateurs positifs sur un espace de mesures, nous montrons l'existence de mesures stationnaires pouvant admettre des masses de Dirac en les traits maximisant la fitness. Lorsque ces mesures admettent une densité continue, nous montrons la convergence des solutions vers cet (unique) état stationnaire.La seconde partie de cette thèse est motivée par un problème issu de la biologie du vieillissement. Nous voulons comprendre l'apparition et le maintien au cours de l'évolution d'un marqueur de sénescence observé chez l'espèce Drosophila melanogaster. Pour cela, nous introduisons un modèle individu-centré décrivant la dynamique d'une population structurée par l'âge et par le trait d'histoire de vie suivant : l'âge de fin de reproduction et celui où la mortalité devient non-nulle. Nous modélisons également l'effet Lansing, qui est l’effet suivant lequel « la descendance de parent jeune vit plus longtemps que celle de parents vieux » . Nous montrons, sous des hypothèses de grande population et de mutations rares, que l'évolution amène ces deux traits à coïncider. Pour cela, nous sommes amenés à étendre l'équation canonique de la dynamique adaptative à une situation où le gradient de fitness n'admet pas des propriétés de régularité suffisantes. L'évolution du trait n'est plus décrite par la trajectoire (unique) d'une équation différentielle ordinaire mais par un ensemble de trajectoires solutions d'une inclusion différentielle. / This thesis is divided into two parts connected by the same thread. It concerns the theoretical study and the application of mathematical models describing population dynamics. The individuals reproduce and die at rates which depend on age a and phenotypic trait. The trait is fixed duringthe life of the individual. It is modified over generations by mutations appearing during reproduction. Natural selection is modeled by introducing a density-dependent mortality rate describing competition for resources.In the first part, we study the long-term behavior of a selection-mutation partial differential equation with age structure describing such a large population. By studying the spectral properties of a family of positive operators on a measures space, we show the existence of stationary measures that can admit Dirac masses in traits maximizing fitness. When these measures admit a continuous density, we show the convergence of the solutions towards this (unique) stationary state.The second part of this thesis is motivated by a problem from the biology of aging. We want to understand the appearance and maintenance during evolution of a senescence marker observed in the species Drosophila melanogaster. For this, we introduce an individual-based model describing the dynamics of a population structured by age and by the following life history trait: the age of reproduction ending and the one where the mortality becomes non-zero. We also model the Lansing effect, which is the effect through which the “progeny of old parents do not live as long as those of young parents”. We show, under large population and rare mutation assumptions, that the evolution brings these two traits to coincide. For this, we are led to extend the canonical equation of adaptive dynamics to a situation where the fitness gradient does not admit sufficient regularity properties. The evolution of the trait is no longer described by the (unique) trajectory of an ordinary differential equation but by a set of trajectories solutions of a differential inclusion.

Ce travail présente quelques résultats sur les systèmes de particules en interaction pour l'interprétation probabiliste des équations aux dérivées partielles, avec des applications à des questions de dynamique moléculaire et de chimie quantique. On présente notamment une méthode particulaire permettant d'analyser le processus de la force biaisante adaptative, utilisé en dynamique moléculaire pour le calcul de différences d'énergies libres. On étudie également la sensibilité de dynamiques stochastiques par rapport à un paramètre, en vue du calcul des forces dans l'approximation de Born-Oppenheimer pour rechercher l'état quantique fondamental de molécules. Enfin, on présente un schéma numérique basé sur un système de particules pour résoudre des lois de conservation scalaires, avec un terme de diffusion anormale se traduisant par une dynamique de sauts sur les particules / This work presents some results on stochastically interacting particle systems and probabilistic interpretations of partial differential equations with applications to molecular dynamics and quantum chemistry. We present a particle method allowing to analyze the adaptive biasing force process, used in molecular dynamics for the computation of free energy differences. We also study the sensitivity of stochastic dynamics with respect to some parameter, aiming at the computation of forces in the Born-Oppenheimer approximation for determining the fundamental quantum state of molecules. Finally, we present a numerical scheme based on a particle system for the resolution of scalar conservation laws with an anomalous diffusion term, corresponding to a jump dynamics on the particles

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 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

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.