La propagation d'une fissure chargée par un écoulement de fluide visqueux est un phénomène complexe où la compréhension des phénomènes mécaniques mis en jeu en pointe de fissures reste encore partielle. C'est le cas de la zone de décollement entre le solide et le fluide qui apparaît pour un certain choix de débit d'injection, de viscosité du fluide et de ténacité du matériau. Cette thèse propose une modélisation simplifiée de ce problème d'interaction fortement couplé. Dans un premier chapitre, on étudie un modèle simplifié unidimensionnel de film élastique collé sur un substrat rigide et on considère une injection de fluide visqueux entre le film et le substrat. On suppose que la propagation de la fissure est régie par la loi de Griffith. On néglige l'existence du retard possible entre le fluide et le solide et on choisit une loi de comportement non-linéaire pour le fluide visqueux. A partir d'une analyse asymptotique pour une faible viscosité, on établit une solution approchée du problème. On montre que le champ pression est singulier en pointe de fissure et on montre l'influence du débit d'injection sur la cinétique du trajet de fissuration. Dans le deuxième chapitre on propose de prendre en compte l'existence de la zone de décollement en modifiant la formulation du modèle et en le réécrivant sous la forme d'un problème d'optimisation en temps discret où les zones de décollement font partie des inconnues du problème. On valide la formulation proposée sur l'exemple analytique de l'écrasement d'une goutte par une barre rigide. On montre ensuite que cette formulation et l'algorithme lié à son implémentation sont capables de gérer l'évolution de l'écrasement de plusieurs gouttes de forme quelconque en capturant correctement les phase d'étalement des gouttes ainsi que de leur coalescence. On étend ensuite cette formulation au cas de l'écrasement d'une goutte par un film élastique. Dans le dernier chapitre, on examine la validité de l'hypothèse de lubrification utilisée en fracturation hydraulique. A l'aide de la méthode de développement asymptotique, on construit une équation de Reynolds régularisée avec des termes de gradient supérieur tenant compte de la variation spatiale de la hauteur des parois. On compare alors le comportement des champs de pression donnés par les équations de Reynolds classique et régularisée sur des exemples d'écoulement entre des conduits de formes multiples. / The crack evolution under a viscous fluid action is a complex phenomenon where the understanding of the mechanical phenomena near the crack tip is still largely limited. This is the case for the lag between the solid and the fluid front propagation which appears for some configurations of injection rate, fluid viscosity and material toughness. This thesis proposes a simplified model for this strongly coupled interaction problem.The first chapter studies a simplified one-dimensional model of a elastic film bonded to a rigid substrate. We consider a viscous fluid injection between the film and the substrate. The crack propagation is assumed to follow the Griffith's law. The existence of the lag is neglected and a non-linear behavior law is chosen for the viscous fluid. Using an asymptotic analysis, an approximate solution is established for the low viscosity case. It is shown that the pressure field diverges at the crack tip and that the kinetics of the crack is influenced by the injection rate. The second chapter proposes to take into account the existence of the lag by modifying the model formulation and rewriting it as a discrete time optimisation problem where the delamination zones are part of the unknowns of the problem. This formulation is validated for the analytical example of a drop crushed by a rigid bar. It is shown that this formulation and its implementation can manage the evolution of several drops of any shape and correctly captures the drops spreading and coalescence. This formulation is then extended to the case of a drop crushed by an elastic film. In the last chapter, the validity of the lubrication hypothesis is examinated. Using an asymptotic analysis, a regularized Reynolds equation is constructed with higher gradient terms taking into account the spatial variation of the walls height. A comparison between the pressure fields behaviour given by the classical and the regularized Reynolds equation is shown for different conducts.

On étudie dans cette thèse certains aspects de la régularité de l'ensemble singulier d'un minimiseur pour la fonctionnelle de Mumford-Shah. On se place principalement en dimension 3 même si certains résultats fonctionnent encore en dimension supérieure. Dans une première partie on étudie les minimiseurs globaux dans R^N et on montre que si (u;K) est un minimiseur global et que si K est un cône assez régulier, alors u (modulo les constantes) est une fonction homogène de degré 1/2 dans R^N\K. Ceci nous permet de lier l'existence d'un minimiseur global et le spectre du laplacien sphérique dans la sphère unité privée de K. Une conséquence est qu'un secteur angulaire stricte ne peut pas être l'ensemble singulier d'un minimiseur global de Mumford-Shah dans R^3. Dans la deuxième partie on montre un théorème de régularité au voisinage des cônes minimaux P, Y et T. On montre que si K est proche (en distance) d'un Y ou d'un T dans une certaine boule, alors K est l'image C^1,alpha d'un P, Y ou d'un T dans une boule légèrement plus petite, ce qui généralise un théorème de L. Ambrosio, N. Fusco et D. Pallara [AFP07]. Les techniques employées ne sont pas exclusives à la dimension 3 et devraient permettre de démontrer des résultats analogues en toute dimension pour un minimiseur de Mumford-Shah, dès lors qu'un résultat de régularité sur les ensembles presque minimaux existerait.

L'étude de l'évolution morphologique d'une surface de cuivre plongée dans un électrolyte liquide et soumise à une dissolution anodique a été abordée selon les trois angles suivants. La modélisation du problème et sa résolution analytique simplifiée ont permis d'introduire un outil de caractérisation et prédiction des instabilités morphologiques globales en fonction des constantes cinétiques associées a l'évolution du système. Un modèle numérique a permis de préciser les rôles respectifs de la réaction d'interface, de la diffusion en phase liquide et de certains paramètres comme le rayon de courbure. Enfin, l'évolution de la forme de l'interface cuivre-électrolyte a été caractérisée par microrugosimétrie à différentes échéances de temps pour des expériences de dissolution menées dans des conditions électrochimiques contrôlées. Si les résultats expérimentaux obtenus ont permis de valider globalement la démarche analytique, les résultats numériques ont montré qu'elle ne pouvait pas rendre compte de toute la complexité des comportements observes. Ainsi, il apparait que la simulation numérique peut constituer une aide précieuse pour la description des aspects morphologiques de l'évolution d'un système

Les déplacements cellulaires, collectifs ou individuels, sont essentiels pour assurer des fonctions fondamentales de l'organisme (réponse immunitaire, morphogenèse), mais jouent également un rôle crucial dans le développement de certaines pathologies (invasion métastatique).Les processus cellulaires à l'origine du déplacement forment une activité complexe, auto-organisée et fortement multi-échelle en temps mais aussi en espace. Mettre en évidence des principes généraux de la migration est donc un enjeu majeur. Dans cette thèse, nous nous intéressons à la construction de modèles de migration individuelle qui prennent en compte ce caractère multi-échelle de manière minimale.Dans une première partie, nous nous intéressons à des modèles particulaires. Nous décrivons des processus intracellulaires clés de la migration de manière discrète au moyen de processus de population. Puis, par une renormalisation en grand nombre d'individus, taille infinitésimale et dynamique accélérée, nous obtenons des équations de dynamique continue et stochastique, permettant de faire le lien entre la dynamique intracellulaire et le déplacement macroscopique.Nous nous confrontons d'abord à la situation d'un leucocyte se déplaçant dans une artère, et développant des liaisons de différentes natures avec les molécules de la paroi, jusqu'à éventuellement s'arrêter. La dynamique de formation de liaisons est décrite par un processus stochastique de type Naissance et Mort avec Immigration. Ces liaisons correspondent à des forces de résistance au mouvement. Nous obtenons explicitement le temps d'arrêt moyen de la cellule.Puis, nous nous intéressons à la reptation cellulaire, qui se produit grâce à la formation d'excroissances au bord de la cellule, appelées protrusions, qui avancent sur le substrat et exercent des forces de traction. Nous modélisons cette dynamique au moyen d'un processus de population structurée par l'orientation de la protrusion. Le modèle continu limite obtenu peut être étudié pour la migration 1D, et donne lieu à une équation de Fokker-Planck sur la distribution de probabilité de la population de protrusion. L'étude d'une configuration stationnaire permet de mettre en avant une dichotomie entre un état non motile et un état de déplacement directionnel.Dans une seconde partie, nous construisons un modèle déterministe minimal de migration dans un domaine discoïdal non déformable. Nous nous basons sur l'idée selon laquelle les structures responsables de la migration renforcent la polarisation de la cellule, ce qui favorise en retour un déplacement directionnel. Cette boucle positive passe par le transport d'un marqueur moléculaire dont la répartition inhomogène caractérise un état polarisé.Le modèle comporte un problème de convection-diffusion sur la concentration en marqueur, où le champs d'advection correspond à la vitesse d'un fluide de Darcy modélisant le cytosquelette. Son caractère actif est porté par des termes de bord, ce qui fait l'originalité du modèle.Du point de vue analytique, le modèle 1D présente une dichotomie face à une masse critique. Dans les cas sous-critique et critique, il est possible de montrer l'existence globale de solutions faibles, ainsi que la convergence à taux explicite vers l'unique état stationnaire correspondant à un état non polarisé. Au delà de la masse critique et pour des masses intermédiaires, nous mettons en évidence deux états stationnaires supplémentaires correspondant à des profils polarisés. De plus, pour des conditions initiales assez asymétrique, nous démontrons l'apparition d'un blow-up en temps fini.Du point de vue numérique, des tests numériques en 2D sont effectués en volumes finis (Matlab) et éléments finis (FreeFem++). Ils permettent de mettre en évidence à nouveau des états motiles et non motiles. L'effet de perturbations stochastiques est étudié, permettant d'aborder des cas de réponse à des signaux extérieurs chimique (chimiotactisme) ou mécanique (obstacle). / Collective or individual cell displacements are essential in fundamental physiological processes (immune response, embryogenesis) as well as in pathological developments (tumor metastasis). The intracellular processes responsible for cell motion have a complex self-organized activity spanning different time and space scales. Highlighting general principles of migration is therefore a challenging task.In a first part, we build stochastic particular models of migration. To do so, we describe key intracellular processes as discrete in space by using stochastic population models. Then, by a renormalization in large population, infinitesimal size and accelerated dynamics, we obtain continuous stochastic equations for the dynamics of interest, allowing a relation between the intracellular dynamics and the macroscopic displacement.First, we study the case of a leukocyte carried by the blood flow and developing adhesive bonds with the artery wall, until an eventual stop. The binding dynamics is described by a stochastic Birth and Death with Immigration process. These bonds correspond to resistive forces to the motion. We obtain explicitly the mean stopping time of the cell.Then, we study the case of cell crawling, that happens by the formation of protrusions on the cell edge, that grow on the substrate and exert traction forces. We describe this dynamics by a structured population process, where the structure comes from the protrusions' orientations. The limiting continuous model can be analytically studied in the 1D migration case, and gives rise to a Fokker-Planck equation on the probability distribution for the protrusion density. For a stationary profile, we can show the existence of a dichotomy between a non motile state and a directional displacement state.In a second part, we build a deterministic minimal migration model in a discoïdal cell domain. We base our work on the idea such that the structures responsible for migration also reinforce cell polarisation, which favors in return a directional displacement. This positive feedback loop involves the convection of a molecular marker, whose inhomogeneous spatial repartition is characteristic of a polarised state.The model writes as a convection-diffusion problem for the marker's concentration, where the advection field is the velocity field of the Darcy fluid that describes the cytoskeleton. Its active character is carried by boundary terms, which makes the originality of the model.From the analytical point of vue, the 1D model shows a dichotomy depending on a critical mass for the marker. In the subcritical and critical cases, it is possible to show global existence of weak solutions, as well as a rate-explicit convergence of the solution towards the unique stationary profile, corresponding to a non-motile state. Above the critical mass, for intermediate values, we show the existence of two additional stationary solutions corresponding to polarised motile profiles. Moreover, for asymmetric enough initial profiles, we show the finite time apparition of a blowup.Studying a more complex model involving activation of the marker at the cell membrane permits to get rid of this singularity.From the numerical point of vue, numerical experiments are led in 2D either in finite volumes (Matlab) or finite elements (FreeFem++) discretizations. They allow to show both motile and non motile profiles. The effect of stochastic fluctuations in time and space are studied, leading to numerical simulations of cases of responses to an external signal, either chemical (chemotaxis) or mechanical (obstacles).

Cette thèse est divisée en deux parties. Dans la première partie, on considère le système de réaction-diffusion-advection (Pε), qui est un modèle d'haptotaxie, mécanisme lié à la dissémination de tumeurs cancéreuses. Le résultat principal concerne la convergence de la solution du systeme (Pε) vers la solution d'un problème à frontière libre (P0) qui est bien défini. Dans la seconde partie, on considère une classe générale d'équations elliptiques du type Hénon:−∆u = |x|^{α} f(u) dans Ω ⊂ R^N avec α > -2. On examine deux cas classiques : f(u) = e^u, |u|^{p−1} u et deux autres cas : f(u) = u^{p}_{+} puis f(u) nonlinéarité générale. En étudiant les solutions stables en dehors d'un ensemble compact (en particulier, solutions stables et solutions avec indice de Morse fini) avec différentes méthodes, on obtient des résultats de classification.

L’objectif de cette étude est de simuler l’érosion d’un sol cohésif sous l’effet d’un écoulement incompressible. Le modèle élaboré décrit une vitesse d’érosion interfaciale qui dépend de la contrainte de cisaillement de l’écoulement. La modélisation numérique proposée est une approche eulérienne, où une méthode de pénalisation de domaines est utilisée pour résoudre les équations de Navier-Stokes autour d’un obstacle. L’interface eau/sol est décrite par une fonction Level Set couplée à une loi d’érosion à seuil.L’approximation numérique est basée sur un schéma DDFV (Discrete Duality Finite Volume) autorisant des raffinements locaux sur maillages non-conformes et non-structurés. L’approche par pénalisation a mis en évidence une couche limite d'inconsistance à l'interface fluide/solide lors du calcul de la contrainte de cisaillement. Deux approches sont proposées pour estimer précisément la contrainte de ce problème à frontière libre. La pertinence du modèle à prédire l’érosion interfaciale du sol est confirmée par la présentation de plusieurs résultats de simulation, qui offrent une meilleure évaluation et compréhension des phénomènes d'érosion / This study focuses on the numerical modelling of the interfacial erosion occurring at a cohesive soil undergoing an incompressible flow process. The model assumes that the erosion velocity is driven by a fluid shear stress at the water/soil interface. The numerical modelling is based on the eulerian approach: a penalization procedure is used to compute Navier-Stokes equations around soil obstacle, with a fictitious domain method, in order to avoid body- fitted unstructured meshes. The water/soil interface’s evolution is described by a Level Set function coupled to a threshold erosion law.Because we use adaptive mesh refinement, we develop a Discrete Duality Finite Volume scheme (DDFV), which allows non-conforming and non-structured meshes. The penalization method, used to take into account a free velocity in the soil with non-body-fitted mesh, introduces an inaccurate shear stress at the interface. We propose two approaches to compute accurately the erosion velocity of this free boundary problem. The ability of the model to predict the interfacial erosion of soils is confirmed by presenting several simulations that provide better evaluation and comprehension of erosion phenomena.

Nous nous intéressons à l'étude d'un problème d'analyse des transferts de chaleur qui modélise une opération de soudage. L'approche que nous considérons ne s'occupe que de la partie solide de la plaque. Elle consiste à résoudre un problème à frontière libre. Pour cela, nous proposons une formulation en optimisation de forme. Le problème d'état est gouverné par un opérateur qui, pour certaines données, n'est pas coercif. Cela complique l'étude de la continuité du problème d'état. Nous surmontons cette difficulté en utilisant le degré topologique de Leray-Shauder, ainsi nous montrons l'existence d'un domaine optimal. Ensuite, nous considérons une discrétisation de ce problème basée sur les éléments finis linéaires. Nous prouvons alors que le problème discret admet une solution et nous montrons qu'une sous-suite des solutions de ce problème convergence vers la solution du problème continu. Enfin, nous présentons des résultats numériques réalisés par deux méthodes : la méthode déterministe basée sur le calcul du gradient de forme, et les algorithmes génétiques combinés avec la logique floue et le calcul parallèle. Ainsi une étude comparative de ces deux méthodes aux niveaux qualitatif et quantitatif a été présentée.

Cette thèse est divisée en deux parties. Dans la première partie, on considère le système de réaction-diffusion-advection (Pε), qui est un modèle d'haptotaxie, mécanisme lié à la dissémination de tumeurs cancéreuses. Le résultat principal concerne la convergence de la solution du systeme (Pε) vers la solution d'un problème à frontière libre (P0) qui est bien défini. Dans la seconde partie, on considère une classe générale d'équations elliptiques du type Hénon:−∆u = |x|^{α} f(u) dans Ω ⊂ R^N avec α > -2. On examine deux cas classiques : f(u) = e^u, |u|^{p−1} u et deux autres cas : f(u) = u^{p}_{+} puis f(u) nonlinéarité générale. En étudiant les solutions stables en dehors d'un ensemble compact (en particulier, solutions stables et solutions avec indice de Morse fini) avec différentes méthodes, on obtient des résultats de classification. / This thesis is divided into two main parts. In the first part, we consider an example of reaction-diffusion-taxis system (Pε), which is a haptotaxis model - a mechanism about the spread of cancer cells. The main result concerns the convergence of the solution of System (Pε) to the solution of a free boundary problem (P0), where system (P0) is well-posed. In the second part, we consider a general class of Hénon type elliptic equations : −∆u = |x|^{α} f(u) in Ω ⊂ R^Nwith α > −2. We investigate two classical cases f(u) = e^u, |u|^{p−1} u and two others cases f(u) = u^{p}_{+} , f(u) is a general function. By studying the solutions which are stable outside a compact set (in particular, stable solutions and finite Morse index solutions) with different methods, we establish some classification results.

Dans le cadre de la thèse des modèles physico-mathématiques pour des microsystèmes électromécaniques avec une permittivité générale sont développés et analysés par des méthodes mathématiques modernes du domaine des équations aux dérivées partielles. En particulier ces systèmes sont à frontière libre et pour conséquence difficiles à traiter. Des méthodes numériques ont été développées pour valider les résultats analytiques obtenus. / In the framework of this thesis physical/mathematical models for microelectromechanical systems with general permittivity have been developed and analysed with modern mathematical methods from the domain of partial differential equations. In particular these systems are moving boundary problems and thus difficult to handle. Numerical methods have been developed in order to validate the obtained analytical results.

Cette thèse est centrée autour de l'étude théorique et de l'analyse numérique des équations paraboliques non linéaires avec divers conditions aux limites. La première partie est consacrée aux équations paraboliques dégénérées mêlant des phénomènes non-linéaires de diffusion et de transport. Nous définissons des notions de solutions entropiques adaptées pour chacune des conditions aux limites (flux nul, Robin, Dirichlet). La difficulté principale dans l'étude de ces problèmes est due au manque de régularité du flux pariétal pour traiter les termes de bords. Ceci pose un problème pour la preuve d'unicité. Pour y remédier, nous tirons profit du fait que ces résultats de régularités sur le bord sont plus faciles à obtenir pour le problème stationnaire et particulièrement en dimension un d'espace. Ainsi par la méthode de comparaison "fort-faible" nous arrivons à déduire l'unicité avec le choix d'une fonction test non symétrique et en utilisant la théorie des semi-groupes non linéaires. L'existence de solution se démontre en deux étapes, combinant la méthode de régularisation parabolique et les approximations de Galerkin. Nous développons ensuite une approche directe en construisant des solutions approchées par un schéma de volumes finis implicite en temps. Dans les deux cas, on combine les estimations dans les espaces fonctionnels bien choisis avec des arguments de compacité faible ou forte et diverses astuces permettant de passer à la limite dans des termes non linéaires. Notamment, nous introduisons une nouvelle notion de solution appelée solution processus intégrale dont l'objectif, dans le cadre de notre étude, est de pallier à la difficulté de prouver la convergence vers une solution entropique d'un schéma volumes finis pour le problème de flux nul au bord. La deuxième partie de cette thèse traite d'un problème à frontière libre décrivant la propagation d'un front de combustion et l'évolution de la température dans un milieu hétérogène. Il s'agit d'un système d'équations couplées constitué de l'équation de la chaleur bidimensionnelle et d'une équation de type Hamilton-Jacobi. L'objectif de cette partie est de construire un schéma numérique pour ce problème en combinant des discrétisations du type éléments finis avec les différences finies. Ceci nous permet notamment de vérifier la convergence de la solution numérique vers une solution onde pour un temps long. Dans un premier temps, nous nous intéressons à l'étude d'un problème unidimensionnel. Très vite, nous nous heurtons à un problème de stabilité du schéma. Cela est dû au problème de prise en compte de la condition de Neumann au bord. Par une technique de changement d'inconnue et d'approximation nous remédions à ce problème. Ensuite, nous adaptons cette technique pour la résolution du problème bidimensionnel. A l'aide d'un changement de variables, nous obtenons un domaine fixe facile pour la discrétisation. La monotonie du schéma obtenu est prouvée sous une hypothèse supplémentaire de propagation monotone qui exige que la frontière libre se déplace dans les directions d'un cône prescrit à l'avance.