Ce travail présente une technique de simulation de feux de forêt qui utilise la méthode Level-Set. On utilise une équation aux dérivées partielles pour déformer une surface sur laquelle est imbriqué notre front de flamme. Les bases mathématiques de la méthode Level-set sont présentées. On explique ensuite une méthode de réinitialisation permettant de traiter de manière robuste des données réelles et de diminuer le temps de calcul. On étudie ensuite l’effet de la présence d’obstacles dans le domaine de propagation du feu. Finalement, la question de la recherche du point d’ignition d’un incendie est abordée. / This work presents a forest fire simulation model which uses the Level-Set method. We use a partial differential equation to deform a surface on which our flame front is inscribed. The mathematical foundations of the Level-set method are presented. We then explain a reinitialization method that allows us to treat in a robust way real data and to reduce the calculation time. The effect of the presence of barriers in the fire propagation domain is also studied. Finally, we make an attempt to find the ignition point of a forest fire.

[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

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.

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.

Cette thèse est consacrée au développement d'outils de traitement pour l'Imagerie par Résonance Magnétique du Tenseur de Diffusion (IRM-TD). Cette technique d'IRM récente est d'une grande importance pour comprendre le fonctionnement du cerveau ou pour améliorer le diagnostic de pathologies neurologiques. Nous proposons des méthodes de traitement basées sur la géométrie Riemannienne, les équations aux dérivées partielles et les techniques de propagation de front. La première partie de ce travail est théorique. Après des rappels sur le système nerveux humain, l'IRM et la géométrie différentielle, nous étudions l'espace des lois normales multivariées. L'introduction d'une structure Riemannienne sur cet espace nous permet de définir des statistiques et des schémas numériques intrinsèques qui sont à la base des algorithmes proposés dans la seconde partie. Les propriétés de cet espace sont importantes pour l'IRM-TD car les tenseurs de diffusion sont les matrices de covariance de lois normales modélisant la diffusion des molécules d'eau en chaque voxel du milieu imagé. La seconde partie est méthodologique. Nous y introduisons des approches originales pour l'estimation et la régularisation d'IRM-TD. Puis nous montrons comment évaluer le degré de connectivité entre aires corticales et introduisons un modèle statistique d'évolution de surface permettant de segmenter ces images. Finalement, nous proposons une méthode de recalage non-rigide. La dernière partie de cette thèse est consacrée à l'analyse des connexions entre le cortex cérébral et les noyaux gris centraux, impliquées dans des tâches motrices, et à l'étude du réseau anatomo-fonctionnel du cortex visuel humain..

La tomographie optique diffuse de fluorescence permet de développer des systèmes d'imagerie fonctionnelle. Elle permet le suivi in vivo de la biodistribution de marqueurs spécifiques fluorescents, dans le but de développer de nouvelles thérapies, en particulier dans le domaine de l'étude de cancers. Le travail de cette thèse est appliqué au problème mathématique de la tomographie, qui peut être formulé en termes d'évaluation de paramètres physiques dans un ensemble d'équations aux dérivées partielles (EDPs).<br />La Méthode des Eléments Finis (MEF) a été choisie ici pour résoudre ces EDPs, car elle permet de ne pas avoir besoin d'hypothèses sur la géométrie et les inhomogénéités du système. Mais elle consomme beaucoup de temps et de mémoire de calcul, principalement en raison des grandes dimensions des matrices impliquées.<br />Notre objectif est d'accélérer la résolution du système.<br />Pour cela, dans un premier temps, une analyse multirésolution a été choisie : les matrices apparaissant dans la version discrétisée des EDPs sont projetées sur une base orthonormale d'ondelettes de Haar ; on obtient alors des approximations des matrices, de taille réduite, contenant l'essentiel de l'information. En utilisant ces approximations à la place des matrices originales, le temps de calcul pour la résolution des EDPs est alors sensiblement réduit.<br />Dans un deuxième temps, la méthode de Galerkin-ondelettes est choisie à la place de la MEF standard : les polynômes d'interpolation classiquement utilisés par la MEF sont directement remplacés par des fonctions d'ondelettes de Daubechies. On s'affranchit alors de l'étape du maillage, coûteuse en temps de calcul, ce qui permet de diminuer également du temps de calcul tout en gardant une qualité comparable de la solution.