[Non fourni] / Acute Myeloid Leukemia (AML) is a cancer of white cells characterized by a quick proliferation of immature cells, that invade the circulating blood and become more present than mature blood cells. This thesis is devoted to the study of two mathematical models of AML. In the first model studied, the cell dynamics are represented by PDE’s for the phases G₀, G₁, S, G₂ and M. We also consider a new phase called Ğ₀, between the exit of the M phase and the beginning of the G₁ phase, which models the fast self-renewal effect of cancerous cells. Then, by analyzing the solutions of these PDE’s, the model has been transformed into a form of two coupled nonlinear systems involving distributed delays. An equilibrium analysis is done, the characteristic equation for the linearized system is obtained and a stability analysis is performed. The second model that we propose deals with a coupled model for healthy and cancerous cells dynamics in AML consisting of two stages of maturation for cancerous cells and three stages of maturation for healthy cells. The cell dynamics are modelled by nonlinear partial differential equations. Applying the method of characteristics enable us to reduce the PDE model to a nonlinear distributed delay system. For an equilibrium point of interest, necessary and sufficient conditions of local asymptotic stability are given. Finally, we derive stability conditions for both mathematical models by using a Lyapunov approach for the systems of PDEs that describe the cell dynamics.

Les systèmes décrits par les équations aux dérivées partielles, appartiennent à la classe des systèmes dynamiques impliquant des fonctions dépendantes de plusieurs variables, comme le temps et l'espace. Déjà fortement répandus pour la modélisation mathématique de phénomènes physiques et environnementaux, ces systèmes ont un rôle croissant dans les domaines de l'automatique. Cette expansion, provoquée par les avancées technologiques au niveau des capteurs facilitant l'acquisition de données et par les nouveaux enjeux environnementaux, incite au développement de nouvelles thématiques de recherche. L'une de ces thématiques, est l'étude des problèmes inverses et plus particulièrement l'identification paramétrique des équations aux dérivées partielles. Tout abord, une description détaillée des différentes classes de systèmes décrits par ces équations est présentée puis les problèmes d'identification qui leur sont associés sont soulevés. L'accent est mis sur l'estimation paramétrique des équations linéaires, homogènes ou non, et sur les équations linéaires à paramètres variant. Un point commun à ces problèmes d'identification réside dans le caractère bruité et échantillonné des mesures de la sortie. Pour ce faire, deux types d'outils principaux ont été élaborés. Certaines techniques de discrétisation spatio-temporelle ont été utilisées pour faire face au caractère échantillonné des données; les méthodes de variable instrumentale, pour traiter le problème lié à la présence de bruit de mesure. Les performances de ces méthodes ont été évaluées selon des protocoles de simulation numérique reproduisant des situations réalistes de phénomènes physique et environnementaux, comme la diffusion de polluant dans une rivière / A large variety of natural, industrial, and environmental systems involves phenomena that are continuous functions not only of time, but also of other independent variables, such as space coordinates. Typical examples are transportation phenomena of mass or energy, such as heat transmission and/or exchange, humidity diffusion or concentration distributions. These systems are intrinsically distributed parameter systems whose description usually requires the introduction of partial differential equations. There is a significant number of phenomena that can be simulated and explained by partial differential equations. Unfortunately all phenomena are not likely to be represented by a single equation. Also, it is necessary to model the largest possible number of behaviors to consider several classes of partial differential equations. The most common are linear equations, but the most representative are non-linear equations. The nonlinear equations can be formulated in many different ways, the interest in nonlinear equations with linear parameters varying is studied. The aim of the thesis is to develop new estimators to identify the systems described by these partial differential equations. These estimators must be adapted with the actual data obtained in experiments. It is therefore necessary to develop estimators that provide convergent estimates when one is in the presence of missing data and are robust to measurement noise. In this thesis, identification methods are proposed for partial differential equation parameter estimation. These methods involve the introduction of estimators based on the instrumental variable technique

Die Dissertation beschäftigt sich mit den Auswirkungen des Klimawandels auf europäische Agrarmärkte im Jahre 2050, unter besonderer Berücksichtigung der Getreide- und Ölsaatenmärkte. Dazu werden die klimabedingten Änderungen der Pflanzenproduktivität des Vegetationsmodells LPJmL, welche auf fünf unterschiedlichen Klimamodellprojektionen basieren, in das Marktmodell ESIM implementiert. ESIM ist ein partielles Gleichgewichtsmodell, welches explizit Agrarmärkte der einzelnen EU-Mitgliedsstaaten simuliert. Zur Berücksichtigung der Unsicherheiten die der Klima-Einfluss-Modellierung zugrunde liegt, werden in dieser Arbeit zwei Ansätze berücksichtigt. Zunächst wird, mittels Gauss-Quadraturen, Stochastizitätin das Marktmodell implementiert, um die Unsicherheit bezüglich klimawandelbedingter steigender Ertragsvariabilität, zu berücksichtigen. Die zweite Methode verwendet die fünf individuellen Produktivitätsänderungen aus dem Vegetationsmodell, woraufhin eine Verteilung der Ergebnisse generiert wird. Darüber hinaus wird das Anpassungsverhalten der Landwirte in das Marktmodell integriert. Dies wird mittels der durch den Klimawandel veänderten Profitabilität der Ackerpflanzen berücksichtigt. Die Ergebnisse weisen darauf hin, dass die Pflanzenproduktivität innerhalb der EU, zumindest bis zum Jahre 2050, weitestgehend positiv vom Klimawandel beeinflusst wird. Die Stärke der Auswirkungen variiert jedoch stark zwischen den einzelnen Ackerpflanzen und Ländern, welche von den zugrundeliegenden Annahmen und Emissionszenarien abhängen. Diese Arbeit leistet einen Beitrag zur aktuellen Klimawandeldiskussion indem sie potentielle Schäden und Nutzen des Klimawandels auf den globalen und den europäischen Agrarsektor quantifizert. Darüber hinaus liefern die stochastische Simulation, sowie die multiplen Simualtionsläufe, ein realistisches Spektrum künftiger potentieller Auswirkungen des Klimawandels. / This study aims to assess potential economic effects of climate change on European agricultural markets at member state level by 2050, focusing on cereal and oilseed markets. The future scenarios include social as well as economic developments derived from two potential emission scenarios. In this modelling framework, crop simulation results of crop productivity changes from the dynamic vegetation model LPJmL, which are based on five individual climate projections, serve as inputs which are administered as a supply shock to the European Simulation Model (ESIM). ESIM is a partial equilibrium model depicting the agricultural sector of the EU in substantial detail. Changes in yields, production quantity and crop prices by the year 2050 are simulated. In order to account for the uncertainty inherent in climate impact assessments, two approaches are considered in this thesis. First, in order to account for climate change increased yield variability, stochasticity is implemented in ESIM, using the method of Gaussian Quadratures. The second method uses the five individual LPJmL outputs to generate a distribution of results. Further, a closely connected purpose of this study is to consider climate change induced adaptation of farmers to changes in the relative profitability of crops. Simulation results indicate, that agricultural productivity in most European countries is positively affected by climate change, at least until the year 2050. However, the degree of impacts vary among crop categories and countries and are also dependent on scenario assumptions. This thesis contributes to the current discussion about climate change impacts by quantifying the potential damages and benefits that may arise from climate change on EU member state level, as well as globally. Further, the stochastic and multiple simulation results based on different future climate and emission projections deliver a more realistic spectrum of potential impacts.

Cette thèse traite de l'étude des équations de Navier-Stokes stationnaires incompressibles et, plus précisément, le comportement quand x→∞ de ses solutions. On étudie la situation dans différents types de domaines non bornés en supposant une condition de nullité à l'infini. On regarde d'abord la dimension 3, dans lequel on sait que si le terme de force décroît très vite à l'infini, le comportement asymptotique est donné par les solutions de Landau, qui sont homogènes de degré -1. On généralise donc ce résultat à des termes de force petits dont le comportement asymptotique est donné par un terme avec l'homogénéité correspondante, c'est-à-dire de degré -3. Pour cela, on trouve une condition nécessaire et suffisante qui est que la partie homogène du terme de force soit de moyenne nulle sur la sphère. Pour finir, on généralise ce résultat au cas d'un domaine extérieur. Dans le cas d'un demi-espace, on va plus loin en montrant que si le terme de force décroit assez à l'infini on obtient des solutions décroissant comme 1/|x|2 à l'infini et on trouve une expression explicite du terme dominant. On peut aussi montrer le même type de résultat que dans l'espace entier avec un terme de force en 1/|x|3 mais la condition de moyenne nulle sur la sphère disparaıt. Dans l'étude de la dimension 2 dans le plan tout entier, on se rend compte que les choses sont plus compliquées. D'abord, pour les solutions homogènes, on arrive à trouver les conditions pour que, si le terme de force est suffisamment petit, on obtienne l'existence de solution qui forment alors une famille à deux paramètres. Mais en leur imposant la restriction d'avoir un flux nul sur le cercle unité, on obtient une famille avec un paramètre seulement. Enfin on étudie les solutions non homogènes, mais pour cela on doit supposer certaines conditions de symétrie sur les données. On trouve alors, pour des termes de force décroissant très vite à l'infini, des solutions en 1/|x|3 et on obtient une formule explicite pour le terme principal de leur développement asymptotique. Ce résultat se généralise aussi au cas d'un domaine extérieur et pour finir, dans ce cadre symétrique, on trouve un résultat analogue au cas de la dimension 3 pour des termes de force qui décroissent en 1/|x|3 à l'infini / This thesis deals with the steady incompressible Navier-Stokes equations, more precisely with the asymptotic behavior of its solutions when |x| → ∞. We consider several types of unbounded domains and we assume that the velocity vanishes at infinity. We first look at the three dimensional case, for which we know that if the forcing term decays fast enough at infinity, the asymptotic behavior of the solutions is given by the Landau solutions that are homogeneous of degree -1. We generalize this result to small forcing terms whose asymptotic behavior at infinity is homogeneous of degree -3. To obtain solutions with an asymptotic behavior at infinity homogeneous of degree -1 we find a necessary and sufficient condition on the forcing : the homogeneous part of the forcing term must have zero mean over the unit sphere. Finally, we generalize this result to the case of an exterior domain. In the case of a half space, we prove that if the forcing term decays sufficiently fast at infinity, then we obtain solutions that decay as 1/|x|2 at infinity and we find an explicit formula for the dominant term in the expansion at infinity of the solution. We can also prove the same type of result as in the full space with forcing terms decaying like 1/|x|3 but the condition of zero mean over the sphere is not required any more. The case of the dimension two is much more difficult. We study first homogeneous solutions and find a family indexed on two real parameters. Imposing the restriction of having zero flux through the unit circle, we get a family of solutions with only one parameter. Finally we deal with non homogeneous solutions, but to do this we need to assume some symmetry conditions on the data. If the forcing term is small and decays sufficiently fast at infinity, we find solutions that decay like 1/|x|3 at infinity and we also obtain an explicit formula for the main term in their asymptotic expansion. We generalize this result to the case of an exterior domain and we also obtain, again under symmetry assumptions, an analogous result to the three dimensional case for forcing terms that decay like 1/|x|3 at infinity

Dans cette thèse, on s'est intéressé à la dynamique des poches de tourbillon pour des équations issues de la mécanique des fluides posées dans le plan. La thèse est composée de trois partie indépendantes. Un des objectifs est d'établir l'existence des tourbillons uniformément concentrés et rigides, c’est-à-dire, qui ne se déforment pas lors de l'évolution. Nous analysons deux configurations liées à la nature topologique du support: poches simplement et doublement connexes. Nos solutions sont obtenues via des techniques de bifurcations et d'analyse complexe. Le deuxième objectif est d'obtenir des précisions sur la structure globale du diagramme de bifurcation et sa réponse vis-à-vis des petites perturbations dans le modèle. Plus précisément, dans le deuxième chapitre on prouve l'existence de V-states doublement connexes dans un voisinage de l'anneau pour le modèle des surfaces quasi-géostrophique. On montre que l'on peut construire des branches de solutions qui sont des anneaux perturbés pour certaines valeurs explicites de vitesses angulaires qui sont liées aux fonctions hypergéométriques de Gauss et aux fonctions de Bessel. Le troisième chapitre porte sur l'étude de la structure du diagramme de bifurcation dans le cas doublement connexes pour l'équation d'Euler. Numériquement, près d'un cas dégénéré, les deux branches issues des deux vitesses angulaires possibles semblaient se rejoindre pour former un lacet. Nous avons prouvé analytiquement ce résultat. Le quatrième chapitre porte sur le modèle shallow water quasi-géostrophique. Dans une première partie, on prouve l'existence de V-states simplement connexes dans un voisinage du tourbillon de Rankine pour un nombre dénombrable de vitesses angulaires liées aux fonctions de Bessel modifiées. La deuxième partie porte sur la réponse du diagramme de bifurcation lorsque l'on fait varier un paramètre du modèle. On montre en particulier qu'une singularité présente lors d'un cas limite est éclatée. Notre étude analytique a été complétée par des simulations numériques portant sur les V-states limites pour les symétries deux et trois. / In this dissertation, we are concerned with the vortex dynamics for some equations arising in fluid mechanics. We distinguish three independent parts. One of the objectives is to prove the existence of uniformly concentrated rigid vortices, they do not change their shapes during the motion. We examine two configurations related to the topological nature of the support: simply and doubly connected vortex patches. Our solutions are obtained using bifurcation arguments and complex analysis tools. The second objective is to obtain some precisions on the global structure of the bifurcation diagram and its response to small perturbations. More precisely, in the second chapter we prove the existence of doubly connected V-states in a neighborhood of the annulus for the surface quasi-geostrophic model. We check that we can construct some branches of solutions which are perturbated annulus at some angular velocities related to hypergeometric Gauss functions and Bessel functions. The goal of the third chapter is to study the structure of the bifurcation diagram in the doubly connected case for Euler equations. Numerically, close to a degenerate case, the two branches of solutions come from the two angular velocities seems to merge to form a loop. We prove analytically this result. In the last chapter, we focus on the shallow quasi-geostrophic model. In the first part, we prove the existence of the simply V-states in a neighborhood of the Rankine Vortices for a countable number of angular velocities related to modified Bessel functions. In the second part, we study the reaction of the diagram bifurcation for small perturbations of the parameter. In particular, we prove that some singularities are broken due to a resonance phenomenon. Our analytical study is completed by numerical simulations on the limiting V-states for the two and three fold symetries.

Ce travail de thèse aborde deux sujets : (i) L'utilisation d'une nouvelle méthode numérique pour l'évaluation des options sur un panier d'actifs, (ii) Le risque de liquidité, la modélisation du carnet d'ordres et la microstructure de marché. Premier thème : Un algorithme glouton et ses applications pour résoudre des équations aux dérivées partielles. L'exemple typique en finance est l'évaluation d'une option sur un panier d'actifs, laquelle peut être obtenue en résolvant l'EDP de Black-Scholes ayant comme dimension le nombre d'actifs considérés. Nous proposons d'étudier un algorithme qui a été proposé et étudié récemment dans [ACKM06, BLM09] pour résoudre des problèmes en grande dimension et essayer de contourner la malédiction de la dimension. L'idée est de représenter la solution comme une somme de produits tensoriels et de calculer itérativement les termes de cette somme en utilisant un algorithme glouton. La résolution des EDP en grande dimension est fortement liée à la représentation des fonctions en grande dimension. Dans le Chapitre 1, nous décrivons différentes approches pour représenter des fonctions en grande dimension et nous introduisons les problèmes en grande dimension en finance qui sont traités dans ce travail de thèse. La méthode sélectionnée dans ce manuscrit est une méthode d'approximation non-linéaire appelée Proper Generalized Decomposition (PGD). Le Chapitre 2 montre l'application de cette méthode pour l'approximation de la solution d'une EDP linéaire (le problème de Poisson) et pour l'approximation d'une fonction de carré intégrable par une somme des produits tensoriels. Un étude numérique de ce dernier problème est présenté dans le Chapitre 3. Le problème de Poisson et celui de l'approximation d'une fonction de carré intégrable serviront de base dans le Chapitre 4 pour résoudre l'équation de Black-Scholes en utilisant l'approche PGD. Dans des exemples numériques, nous avons obtenu des résultats jusqu'en dimension 10. Outre l'approximation de la solution de l'équation de Black-Scholes, nous proposons une méthode de réduction de variance des méthodes Monte Carlo classiques pour évaluer des options financières. Second thème : Risque de liquidité, modélisation du carnet d'ordres, microstructure de marché. Le risque de liquidité et la microstructure de marché sont devenus des sujets très importants dans les mathématiques financières. La dérégulation des marchés financiers et la compétition entre eux pour attirer plus d'investisseurs constituent une des raisons possibles. Dans ce travail, nous étudions comment utiliser cette information pour exécuter de façon optimale la vente ou l'achat des ordres. Les ordres peuvent seulement être placés dans une grille des prix. A chaque instant, le nombre d'ordres en attente d'achat (ou vente) pour chaque prix est enregistré. Dans [AFS10], Alfonsi, Fruth et Schied ont proposé un modèle simple du carnet d'ordres. Dans ce modèle, il est possible de trouver explicitement la stratégie optimale pour acheter (ou vendre) une quantité donnée d'actions avant une maturité. L'idée est de diviser l'ordre d'achat (ou de vente) dans d'autres ordres plus petits afin de trouver l'équilibre entre l'acquisition des nouveaux ordres et leur prix. Ce travail de thèse se concentre sur une extension du modèle du carnet d'ordres introduit par Alfonsi, Fruth et Schied. Ici, l'originalité est de permettre à la profondeur du carnet d'ordres de dépendre du temps, ce qui représente une nouvelle caractéristique du carnet d'ordres qui a été illustré par [JJ88, GM92, HH95, KW96]. Dans ce cadre, nous résolvons le problème de l'exécution optimale pour des stratégies discrètes et continues. Ceci nous donne, en particulier, des conditions suffisantes pour exclure les manipulations des prix au sens de Huberman et Stanzl [HS04] ou de Transaction-Triggered Price Manipulation (voir Alfonsi, Schied et Slynko) / This work is organized in two themes : (i) A novel numerical method to price options on manyassets, (ii) The liquidity risk, the limit order book modeling and the market microstructure.First theme : Greedy algorithms and applications for solving partial differential equations in high dimension Many problems of interest for various applications (material sciences, finance, etc) involve high-dimensional partial differential equations (PDEs). The typical example in finance is the pricing of a basket option, which can be obtained by solving the Black-Scholes PDE with dimension the number of underlying assets. We propose to investigate an algorithm which has been recently proposed and analyzed in [ACKM06, BLM09] to solve such problems and try to circumvent the curse of dimensionality. The idea is to represent the solution as a sum of tensor products and to compute iteratively the terms of this sum using a greedy algorithm. The resolution of high dimensional partial differential equations is highly related to the representation of high dimensional functions. In Chapter 1, we describe various linear approaches existing in literature to represent high dimensional functions and we introduce the high dimensional problems in finance that we will address in this work. The method studied in this manuscript is a non-linear approximation method called the Proper Generalized Decomposition. Chapter 2 shows the application of this method to approximate the so-lution of a linear PDE (the Poisson problem) and also to approximate a square integrable function by a sum of tensor products. A numerical study of this last problem is presented in Chapter 3. The Poisson problem and the approximation of a square integrable function will serve as basis in Chapter 4for solving the Black-Scholes equation using the PGD approach. In numerical experiments, we obtain results for up to 10 underlyings. Second theme : Liquidity risk, limit order book modeling and market microstructure. Liquidity risk and market microstructure have become in the past years an important topic in mathematical finance. One possible reason is the deregulation of markets and the competition between them to try to attract as many investors as possible. Thus, quotation rules are changing and, in general, more information is available. In particular, it is possible to know at each time the awaiting orders on some stocks and to have a record of all the past transactions. In this work we study how to use this information to optimally execute buy or sell orders, which is linked to the traders' behaviour that want to minimize their trading cost. In [AFS10], Alfonsi, Fruth and Schied have proposed a simple LOB model. In this model, it is possible to explicitly derive the optimal strategy for buying (or selling) a given amount of shares before a given deadline. Basically, one has to split the large buy (or sell) order into smaller ones in order to find the best trade-off between attracting new orders and the price of the orders. Here, we focus on an extension of the Limit Order Book (LOB) model with general shape introduced by Alfonsi, Fruth and Schied. The additional feature is a time-varying LOB depth that represents a new feature of the LOB highlighted in [JJ88, GM92, HH95, KW96]. We solve the optimal execution problem in this framework for both discrete and continuous time strategies. This gives in particular sufficient conditions to exclude Price Manipulations in the sense of Huberman and Stanzl [HS04] or Transaction-Triggered Price Manipulations (see Alfonsi, Schied and Slynko). The seconditions give interesting qualitative insights on how market makers may create price manipulations

Les performances d’un système de reconnaissance de formes dépendent en bonne partie de la qualité de l’image segmentée. Malgré les progrès effectués, une segmentation complète (c’est-à-dire avec des contours entiers) ne peut pas être toujours atteinte. Dans un premier temps, nous nous plaçons dans le cas où seulement certaines parties de la forme entière sont disponibles. D’abord, afin d’assurer l’invariance des parties de formes aux transformations géométriques, nous développons une méthode d’optimisation de l’analyse procustéenne qui consiste à retrouver les points extrémités optimaux qui minimisent la distance de Procutse. Ensuite, nous proposons une approche de reconnaissance de parties de formes et une approche de reconnaissance de formes partielles. Ces deux approches, basées-contour, sont fondées sur un recalage robuste entre les parties de formes. La méthode de recalage proposée consiste à optimiser une mesure de similarité basée sur les géodésiques dans l’espace de formes. Ainsi, nous exploitons le résidu du recalage pour définir une nouvelle métrique pour la reconnaissance de parties de formes. Puis, nous décrivons une stratégie de combinaison avec cette même métrique pour la reconnaissance de formes partielles. Par la suite, nous proposons d’utiliser la distance géodésique proposée pour la reconnaissance des parties de formes dans la définition d’une métrique globale pour la reconnaissance de formes entières. Les tests de reconnaissance (classification et recherche) sont effectués sur des parties requêtes et des formes entières de la base d’images MPEG-7, puis sur des images réelles segmentées. Les résultats expérimentaux montrent la supériorité de nos méthodes par rapport aux autres méthodes de l’état de l’art. / The quality of the segmentation process directly affects the performance of the shape recognition. Despite the progress that has been made, it is often unreachable to segment the entire object (i.e. closed contour). In fact, only some parts/fragments of objects can be detected. We first develop a new alignment method based on Procrustes analysis in order to ensure invariance of shape parts to geometric transformations (translation, rotation and scale factor). The proposed method consists in finding optimal extremities which minimize the Procrustes distance. Then, we propose a shape part recognition approach and a partial shape recognition approach. These two contour-based approaches are based on matching between shape parts to compare. This matching process consists in establishing a robust registration between shape parts based on geodesics in the shape space. Thus, we exploit the registration residual to define a novel distance for shape part recognition. Later, for partial shape recognition, we describe a geodesics-based combining strategy with the same distance. As well, we propose to use the geodesics distance proposed for shape part recognition to define a global distance for entire shape recognition. Experiments are carried out on parts of shapes and entire shapes of theMPEG-7 database, then on parts issued from segmented real images. The obtained results demonstrate the effectiveness of our proposed recognition schemes. The proposed approaches are shown to significantly outperform previous works for classification and retrieval applications.

On s'interesse aux problèmes elliptiques avec coéficients variables à travers des interfaces intérieures. La solution et ses dérivées normales peuvent subir des variations significatives à travers les frontières intérieures. On présente une méthode compacte aux différences finies sur des maillages adaptés de type octree conçues pour une résolution en parallèle. L'idée principale est de minimiser l'erreur de troncature sur la discretisation locale, en fonction de la configuration du maillage, en rapprochant une convergence à l'ordre deux. On montrera des cas 2D et 3D des résultat liés à des applications concrètes. / We consider problems governed by a linear elliptic equation with varying coéficients across internal interfaces. The solution and its normal derivative can undergo significant variations through these internal boundaries. We present a compact finite-difference scheme on a tree-based adaptive grid that can be efficiently solved using a natively parallel data structure. The main idea is to optimize the truncation error of the discretization scheme as a function of the local grid configuration to achieve second order accuracy. Numerical illustrations relevant for actual applications are presented in two and three-dimensional configurations.

La perméabilisation des cellules à l’aide d’impulsions électriques intenses, appelée électroperméabilisation, est un phénomène biologique impliqué dans des thérapies anticancéreuses récentes. Elle permet, par exemple, d’améliorer l’efficacité d’une chimiothérapie en diminuant les effets secondaires, d’effectuer des transferts de gènes, ou encore de procéder à l’ablation de tumeurs. Les mécanismes de l’électroperméabilisation restent cependant encore méconnus, et l’hypothèse majoritairement admise par la communauté de formation de pores à la surface des membranes cellulaires est en contradiction avec certains résultats expérimentaux.Le travail de modélisation proposé dans cette thèse est basé sur une approche différente des modèles d’électroporation existants. Au lieu de proposer des lois sur les propriétés des membranes à partir d’hypothèses à l’échelle moléculaire, nous établissons des lois ad hoc pour les décrire, en se basant uniquement sur les informations expérimentales disponibles. Aussi, afin de rester au plus prèsde ces dernières et faciliter la phase de calibration à venir, nous avons ajouté un modèle de transport et de diffusion de molécules dans la cellule. Une autre spécificité de notre modèle est que nous faisons la distinction entre l’état conducteur et l’état perméable des membranes.Des méthodes numériques spécifiques ainsi qu’un code en 3D et parallèle en C++ ont été écrits et validés pour résoudre les équations aux dérivées partielles de ces différents modèles. Nous validons le travail de modélisation en montrant que les simulations reproduisent qualitativement les comportements observés in vitro. / Cell permeabilization by intense electric pulses, called electropermeabilization, is a biological phenomenon involved in recent anticancer therapies. It allows, for example, to increase the efficacy of chemotherapies still reducing their side effects, to improve gene transfer, or to proceed tumor ablation. However, mechanisms of electropermeabilization are not clearly explained yet, and the mostly adopted hypothesis of the formation of pores at the membrane surface is in contradiction with several experimental results.This thesis modeling work is based on a different approach than existing electroporation models. Instead of deriving equations on membranes properties from hypothesis at the molecular scale, we prefer to write ad hoc laws to describe them, based on available experimental data only. Moreover, to be as close as possible to these data, and to ease the forthcoming work of parameter calibration, we added to our model equations of transport and diffusion of molecules in the cell. Another important feature of our model is that we differentiate the conductive state of membranes from their permeable state.Numerical methods, as well as a 3D parallel C++ code were written and validated in order to solve the partial differential equations of our models. The modeling work was validated by showing qualitative match between our simulations and the behaviours that are observed in vitro

Nous analysons trois modèles non-locaux décrivant la dynamique évolutive d’un trait phénotypique continu soumis à l’action conjointe des mutations et de la sélection. Nous établissons l’existence et l’unicité des solutions du problème de Cauchy, et donnons la description du comportement en temps long de la solution. Dans le premier travail nous étudions l’équation du réplicateur-mutateur en domaine non borné et généralisons aux cas des valeurs sélectives confinantes les résultats connus dans le cas harmonique. À savoir, l’existence d’une unique solution globale, régulière, convergeant en temps long vers un profil universel ; pour cela, nous employons des techniques de décomposition spectrale d’opérateurs de Schrödinger. Le deuxième travail traite d’un modèle dont la valeur sélective est densité-dépendante. Afin de montrer le caractère bien posé de l’équation, nous combinons deux approches. La première est basée sur l’étude de la fonction génératrice des cumulants, satisfaisant une équation de transport non locale et permettant d’obtenir implicitement le trait moyen. La deuxième exploite un changement de variable (formule d’Avron-Herbst), permettant d’écrire la solution en termes du trait moyen et de la solution de l’équation de la chaleur avec même donnée initiale. Finalement, nous étudions un modèle dont le taux de mutation est proportionnel à la valeur moyenne du trait. Nous établissons un lien bijectif entre ce dernier modèle et le deuxième, permettant ainsi de décrire finement la dynamique de la solution. Nous montrons en particulier la croissance exponentielle du trait moyen. / We analyze three non-local models describing the evolutionary dynamics of a continuous phenotypic trait undergoing the joint action of mutations and selection. We establish the existence and uniqueness of the solutions to the Cauchy problem, and give a description of the long-time behaviour of the solution. In the first work we study the replicator-mutator equation in the unbounded domain and generalize to cases of selective values confining the known results in the harmonic case. Namely, the existence of a unique global regular solution, converging towards a universal profile; for this, we use spectral decomposition techniques of Schrödinger operators. In the second work, we discuss a model whose fitness value is density-dependent. In order to show the well-posedness of the equation, we combine two approaches. The first is based on the study of the cumulant generating functions, satisfying a non-local transport equation and making it possible to implicitly obtain the average trait. The second uses a change of variable (Avron-Herbst formula), allowing the solution to be written in terms of the average trait and the solution of the heat equation with the same initial data. Finally, we study a model whose mutation rate is proportional to the average value of the trait. We establish a bijective link between this last model and the second, thus making it possible to describe the dynamics of the solution in detail. In particular, we show the exponential growth of the average trait.