Spelling suggestions: "subject:"equations aux dérivés partielle"" "subject:"equations aux dérivés partiellement""
61 |
Modélisation et simulation de l'effet Leidenfrost dans les micro-gouttesDenis, Roland 26 November 2012 (has links) (PDF)
L'effet Leidenfrost répresente un cas particulier de caléfaction : lorsqu'une goutte de liquide est déposée sur une surface dont la température est très supérieure à la température d'ébullition du liquide, ce dernier s'évapore avant de toucher la surface et la vapeur ainsi créée forme un coussin sous la goutte qui la maintient en sustentation et l'isole de la plaque chauffante. Ce travail de thèse concerne la modélisation et la simulation de ce phénomène complexe. Dans une première partie, nous étudions un modèle avec interface raide basée sur les équations de Navier-Stokes enrichies avec des termes interfaciaux prenant en compte le changement de phase et la tension de surface. La simulation d'une couche uniforme de liquide sur un film de vapeur nous ramène à un cas unidimensionnel pour lequel on utilise la méthode ALE (Arbitrary Lagrangian Eulerian) afin de gérer la hauteur variable de chaque phase. La discrétisation du modèle est validée sur un cas test. Dans une seconde partie, on utilise la méthode de capture d'interface Level-Set dans laquelle la frontière liquide/gaz est réprésentée par la ligne de niveau zéro d'une fonction. Cette interface est artificiellement épaissie et les quantités thermodynamiques y sont régularisées. La tension de surface et le changement de phase sont alors introduits sous forme de termes volumiques dans nos équations. L'hypothèse d'incompressibilité de chaque phase pure nous amène alors à un fluide généralisé dont la compressibilité se manifeste uniquement dans la zone interfaciale, là où se produit le changement de phase. La troisième partie est consacrée à la discrétisation de ce modèle pour l'étude tridimensionnelle d'une goutte d'eau, immobile et symétrique par rotation, se ramenant ainsi à un problème bi-dimensionnel axisymétrique. La méthode Level-Set nécessite des choix numériques particuliers qui sont alors explicités : schéma d'advection peu diffusif, redistanciation par résolution d'une équation de Hamilton-Jacobi et correction globale du volume de la goutte, prenant en compte le changement de phase. Un algorithme de projection de type Chorin est également utilisé afin de satisfaire la contrainte sur la compressibilité de notre fluide. On présentera également un nouveau schéma aux différences finies à stencil compact pour l'approximation du gradient. La dernière partie présente et compare nos résultats numériques avec plusieurs courbes théoriques, décrivant chacune l'évolution de certains paramètres de la goutte de liquide : son volume, son rayon et la hauteur de la couche de vapeur.
|
62 |
Contrôle non destructif via des sondes courants de Foucault : nouvelles approchesJiang, Zixian 23 January 2014 (has links) (PDF)
L'objectif principal de cette thèse est de proposer et de tester quelques méthodes de l'optimisation de forme afin d'identifier et de reconstruire des dépôts qui couvrent la paroi extérieure d'un tube conducteur dans un générateur de vapeurs en utilisant des signaux courant de Foucault. Ce problème est motivé par des applications industrielles en contrôle non-destructive dans le secteur de l'énergie nucléaire. En fait, des dépôts peuvent obstruer le passage de circuit de refroidissement entre les tubes et les plaques entretoises qui les soutiennent, ce qui entraînerait une baisse de productivité et mettrait la structure en danger. On considère dans un premier temps un cas axisymétrique dans le cadre duquel on construit un modèle 2-D par des équations aux dérivées partielles pour le courant de Foucault, ce qui nous permet ensuite de reproduire des mesures d'impédances qui correspondent en réalité les signaux courants de Foucault. Pour réduire le coût de calculs de la simulation numérique, on tronque le domaine du problème en posant des conditions aux bords artificielles basées sur des études sur le comportement de la solution, notamment sur un calcul semi-analytique de l'opérateur Dirichlet-to-Neumann dans la direction axiale. Pour la partie identification et reconstruction, on distingue deux sortes de dépôts et établit pour chacun une méthode d'inversion spécifique. Le premier cas concernent des dépôts dont la conductivité est relativement faible (d'environs 1.e4 S/m). On utilise la méthode d'optimisation de forme qui consiste à exprimer explicitement la dérivée des mesures d'impédance par rapport au domaine du dépôt en fonction d'une déformation et à représenter le gradient d'un fonctionnel de coût à minimiser par l'intermédiaire d'un état adjoint proprement défini. Motivée par la présence des dépôts et des plaques de maintient non-axisymétriques, on fait aussi une extension en 3-D des méthodes précédentes. Pour le deuxième cas, des dépôts sont fortement conducteurs et sous forme de couche mince d'une épaisseur de l'ordre de micron. La méthode adaptée à la première sorte de dépôts devient ici trop coûteuse à cause du degré de liberté très élevé des éléments finis sur un maillage extrêmement raffiné à l'échelle de la couche mince. Pour relever cette difficulté, les études sont portées sur plusieurs modèles asymptotiques qui remplace la couche mince par des conditions d'interface sur la surface du tube portante du dépôt. Le nom de modèle asymptotique vient du fait que les conditions d'interface effectives sont obtenues par le développement asymptotique de la solution en fonction d'un paramètre caractérisant la conductivité et l'épaisseur de la couche mince. La validation numérique a permis de retenir un modèle asymptotique qui est à la fois suffisamment précis et simple à inverser. L'inversion (recherche de l'épaisseur du dépôt) consiste alors à rechercher des paramètres dans les conditions d'interface (non standard). Pour les deux cas, la validation et des exemples numériques sont proposés pour le modèle direct et l'inversion.
|
63 |
Contrôle d'équations de Schrödinger et d'équations paraboliques dégénérées singulièresMorancey, Morgan 27 November 2013 (has links) (PDF)
Ce mémoire présente les travaux réalisés au cours de ma thèse sur le contrôle d'équations aux dérivées partielles. La première partie est consacrée à l'étude d'équations de Schrödinger bilinéaires unidimensionnelles autour de deux axes : la non contrôlabilité en temps petit avec des contrôles petits et la contrôlabilité simultanée. On établit un cadre pour lequel, bien que la vitesse de propagation du système soit infinie, la contrôlabilité exacte locale avec des contrôles petits est vérifiée si et seulement si le temps est suffisamment grand. Ces résultats, basés sur la coercivité d'une forme quadratique associée au développement de l'état à l'ordre deux, sont étendus dans le contexte de la contrôlabilité simultanée. On montre alors, en utilisant la méthode du retour de J.-M.~Coron, des résultats de contrôle exact local simultané pour deux ou trois équations, à phase globale et/ou à retard global près. La trajectoire de référence utilisée est construite via des résultats de contrôle partiel. En utilisant un argument de perturbation, on étend cette idée pour montrer la contrôlabilité exacte globale d'un nombre quelconque d'équations sans hypothèse sur le potentiel. Dans la deuxième partie, on prend en compte dans le modèle un terme supplémentaire quadratique en le contrôle. Ce terme, dit de polarisabilité, généralement négligé, présente un intérêt physique dans la modélisation, mais aussi mathématique dans le cas où le terme bilinéaire est insuffisant pour conclure à la contrôlabilité. En dimension quelconque, on construit des contrôles explicites réalisant la contrôlabilité approchée de l'état fondamental. En adaptant conjointement l'argument de perturbation précédent et certains résultats du contrôle bilinéaire, on prouve la contrôlabilité globale exacte de l'équation de Schrödinger avec polarisabilité 1D. La dernière partie de ce mémoire est consacrée à l'étude de la continuation unique pour un opérateur de type Grushin sur un rectangle 2D. Cet opérateur présente une singularité et une dégénérescence sur un segment séparant le domaine en deux composantes. On donne une condition nécessaire et suffisante sur le coefficient du potentiel singulier pour obtenir la continuation unique.
|
64 |
Modélisation, analyse mathématique et simulations numériques de quelques problèmes aux dérivées partielles multi-échellesRambaud, Amélie 05 December 2011 (has links) (PDF)
Nous étudions plusieurs aspects d'équations aux dérivées partielles multi-échelles. Pour trois exemples, la présence de multiples échelles, spatiales ou temporelles, motive un travail de modélisation mathématique ou constitue un enjeu de discrétisation. La première partie est consacrée à la construction et l'étude d'un système multicouche de type Saint-Venant pour décrire un fluide à surface libre (océan). Son obtention s'appuie sur l'analyse des échelles spatiales, précisément l'hypothèse " eau peu profonde ". Nous justifions nos équations à partir du modèle primitif et montrons un résultat d'existence locale de solution. Puis nous proposons un schéma volumes finis et des simulations numériques. Nous étudions ensuite un problème hyperbolique de relaxation, inspiré de la théorie cinétique des gaz. Nous construisons un schéma numérique via une stratégie préservant l'asymptotique : nous montrons sa convergence pour toute valeur du paramètre de relaxation, ainsi que sa consistance avec le problème à l'équilibre local. Des estimations d'erreurs sont établies et des simulations numériques sont présentées. Enfin, nous étudions un problème d'écoulement sanguin dans une artère avec stent, modélisé par un système de Stokes dans un domaine contenant une petite rugosité périodique (géométrie double échelle). Pour éviter une discrétisation coûteuse du domaine rugueux (l'artère stentée), nous formulons un ansatz de développement de la solution type Chapman-Enskog, et obtenons une loi de paroi implicite sur le bord du domaine lisse (artère seule). Nous montrons des estimations d'erreurs et des simulations numériques
|
65 |
Contributions à l'étude de l'instant de défaut d'un processus de Lévy en observation complète et incomplète / Contributions to the study of default time of a Lévy process in complete observation and in incomplete ObservationNgom, Waly 06 July 2016 (has links)
Dans nos travaux, nous avons considéré un processus de Lévy X avec une composante brownienne non nulle et dont la partie à sauts est un processus de Poisson composé. Nous avons supposé que la valeur d'une entreprise est modélisée par un processus stochastique de la forme V = Vo exp X et que cette entreprise est mise à défaut dès lors que sa valeur passe sous un certain seuil b déterminé de façon exogène et qui donc, est une donnée du problème. L'instant de défaut T est alors de la forme Tx pour x= ln(Vo) ln((b) où x> 0, Tx = inf{t 2:0: X, 2:x}. Dans un premier temps, nous supposons que des agents observant la valeur V des actifs de la firme souhaitent connaître le comportement de l'instant de défaut. Dans ce modèle, au chapitre 2, nous avons étudié d'une part la régularité de la densité de la loi de l'instant de défaut. D'autre part, nous avons étudié la loi conjointe de l'instant de défaut, de l'overshoot et de l'undershoot. Au chapitre 3, nous avons obtenu une équation à valeurs mesures dont le quadriplet formé par la variable aléatoire X,, le su premum du processus X à l'instant t, le supremum du processus X au dernier instant de saut avant l'instant t et le dernier instant de saut à l'instant t est solution au seris faible, puis une équation dont ce quadriplet est une solution forte. Dans un second temps, au chapitre 4, nous avons supposé que des investisseurs souhaitant détenir une part de cette entreprise ne disposent pas de l'information complète. Ils n'observent pas la valeur des actifs de la firme V, mais sa valeur bruitée. Leur information est modélisée par la filtration Ç = (Ç,, t 2: 0) engendrée par cette observation. Dans ce modèle, nous avons montré que la loi conditionnelle de l'instant de défaut sachant la tribu Ç, admet une densité par rapport à la mesure de Lebesgue et obtenu une équation de Volttera dont cette densité est solution. Cette connaissance permet aux investisseurs de prévoir au vu de leur information, quand est-ce que l'instant de défaut va intervenir après l'instant t. Nous avons complété ce travail par des simulations numériques. / In this Ph.D thesis, we consider a jump-diffusion process which the diffusion part is a drifted Brownian motion and the jump part is a compound Poisson process. We assume that a firm value is modelling by a stochastic process V = V0 exp-X. This firm goes to default whenever its value is below a specified tlrreshold b which is exo genously determined. For x = ln(Vo) - ln(b) > 0, the default time is of the form Tx = inf{t 2:0: X, 2: x}. First, we suppose that agents observe perfectly the firm value. In this mode, we sho wed in chapter 2 that the density of the default time is continuons, then study the joint law of the default time, overshoot an undershoot. We obtained in chapter 3 a valued measure differentia equation which the solution is the quadruplet formed by the random variable X,, the running supremum x; of X at time t, the supremum of X at the last jump time before t and the last jump time before t. Secondly, we assume that investors wishing detain a part of the firm can not observe the firm value. They observe a noisy value of the firm and their information is madel ling by the filtration g = (9,,t 2: 0) generated by their observation. In this mode, we have shown that the conditional density of Tx with respect to Ç has a density which is solution of one stochastic integral-differentia equation The knowledge of this density allows investors to predict the default time after time t. This second part is the chapter 4.
|
66 |
Flot de Ricci sans borne supérieure sur la courbure et géométrie de certains espaces métriques / Ricci flow without upper bounds on the curvature and the geometry of some metric spaces.Richard, Thomas 21 September 2012 (has links)
Le flot de Ricci, introduit par Hamilton au début des années 80, a montré sa valeur pour étudier la topologie et la géométrie des variétés riemanniennes lisses. Il a ainsi permis de démontrer la conjecture de Poincaré (Perelman, 2003) et le théorème de la sphère différentiable (Brendle et Schoen, 2008). Cette thèse s'intéresse aux applications du flot de Ricci à des espaces métriques à courbure minorée peu lisses. On définit en particulier ce que signifie pour un flot de Ricci d'avoir pour condition initiale un espace métrique. Dans le Chapitre 2, on présente certains travaux de Simon permettant de construire un flot de Ricci pour certains espaces métriques de dimension 3. On démontre aussi deux applications de cette construction : un théorème de finitude en dimension 3 et une preuve alternative d'un théorème de Cheeger et Colding en dimension 3. Dans le Chapitre 3, on s'intéresse à la dimension 2. On montre que pour les surfaces singulières à courbure minorée (au sens d'Alexandrov), on peut définir un flot de Ricci et que celui-ci est unique. Ceci permet de montrer que l'application qui à une surface associe son flot de Ricci est continue par rapport aux perturbations Gromov-Hausdorff de la condition initiale. Le Chapitre 4 généralise une partie de ces méthodes en dimension quelconque. On doit y considérer des conditions de courbure autres que les usuelles minorations de la courbure de Ricci ou de la courbure sectionnelle. Les méthodes mises en place permettent de construire un flot de Ricci pour certains espaces métriques non effondrés limites de variétés dont l'opérateur de courbure est minoré. On montre aussi que sous certaines hypothèses de non-effondrement, les variétés à opérateur de courbure presque positif portent une métrique à opérateur de courbure positif ou nul. / The Ricci flow was introduced by Hamilton in the beginning of the 90's. It has been a valuable tool to study the topology and the geometry of smooth Riemannian manifolds. For example, it was essential in the of the Poincaré conjecture (Perelman, 2003) and of the differentiable sphere theorem (Brendle and Schoen, 2008). In this thesis, we are interested in the applications of Ricci flow to metric spaces with curvature bounded from below which are not smooth. We define what it means for a Ricci flow to admit a metric space as initial condition. In Chapter 2, we present some works of Simon which allow to build a Ricci flow for some metric spaces of dimension 3. We also give two applications of this result : a finiteness theorem in dimension 3 and an alternative of a theorem of Cheeger and Colding in dimension 3. In Chapter 3, we treat the special case of dimension 2. We show that for singular surfaces whose curvature is boded from below (in the sense of Alexandrov), we can define a Ricci and it is unique. This allow to show that for surfaces with curvature bounded from below, the application which maps a surface to its Ricci flow is continuous with respect to Gromov-Hausdorff perturbations of the initial condition. Chapter 4 generalizes some of these methods in higher dimension. Here one needs to consider other conditions on the curvature than the usual "Ricci curvature bounded from below" and "sectional curvature bounded from below". The methods used there allow us to build a Ricci flow for some non-collapsed metric spaces which are limits of manifolds whose curvature operator is bounded from below. We also show that under some non-collapsing assumptions manifolds with almost non-negative curvature operator admit metrics with non-negative curvature operator.
|
67 |
Modélisation, analyse et simulation numérique de solides combinant plasticité, rupture et dissipation visqueuse / Modeling, analysis and numerical simulation of solids combining plasticity, fracture and viscous dissipationJakabčin, Lukáš 22 September 2014 (has links)
Dans cette thèse nous nous intéressons à la modélisation, analyse mathématique et simulation numérique d'une classe de modèles combinant différents phénomènes dissipatifs liés à la plasticité, rupture et dissipation visqueuse.Tout d'abord, nous construisons des modèles d'évolution contenant plasticité, viscoplasticité, écrouissage cinématique linéaire et rupture. En particulier, nous montrons une inégalité thermodynamique de type Clausius-Duhem pour nos modèles. Ensuite, nous montrons l'existence d'évolutions pour deux modèles: celui d'élasto-visco-plasticité avec la rupture approchée via la fonctionnelle Ambrosio-Tortorelli et celui d'élasto-viscoplasticité avec écrouissage cinématique linéaire et rupture approchée basée sur l'utilisation de la fonctionnelle d'Ambrosio-Tortorelli avec un r-Laplacien. Enfin, nous étudions numériquement nos modèles en fonction de différents paramètres mécaniques. Nous proposons aussi une extension de la méthode numérique de backtracking aux matériaux à mémoire. Au final, nous effectuons des comparaisons numériques entre un de nos modèles et l'expérience géophysique de plasticine de Peltzer et Tapponnier qui modélise la propagation des failles dans la crôute terrestre. / In this work, we are interested in modeling, mathematical analysis and numerical simulation of a class of models that combine several mecanisms of dissipation: plasticity, fracture and viscous dissipation. Firslty, we construct evolution models containing plasticity, viscoplasticity, linear kinematic hardening and fracture. In particular, we show for our models a Clausius-Duhem like thermodynamical inequality. Then, we prove an existence result for evolutions for an elasto-visco-plastic model with regularized fracture using the Ambrosio-Tortorelli functional and for an elasto-viscoplastic model with kinematic hardening and fractures regularized with the modified r-Laplacian Ambrosio-Tortorelli functional. Finally, we study from a numerical point of view our models in function of various mecanical parameters. We also propose an extension of the backtracking algorithm for materials with memory. In the end, we test numerically one of our models on a geophysical Peltzer and Tapponnier's experiment of plasticine that models failure propagation in the Earth crust.
|
68 |
Analytical properties of viscosity solutions for integro-differential equations : image visualization and restoration by curvature motions / Propriétés analytiques des solutions de viscosité des équations integro-différentielles : visualisation et restauration d'images par mouvements de courbureCiomaga, Adina 29 April 2011 (has links)
Le manuscrit est constitué de deux parties indépendantes.Propriétés des Solutions de Viscosité des Equations Integro-Différentielles.Nous considérons des équations intégro-différentielles elliptiques et paraboliques non-linéaires (EID), où les termes non-locaux sont associés à des processus de Lévy. Ce travail est motivé par l'étude du Comportement en temps long des solutions de viscosité des EID, dans le cas périodique. Le résultat classique nous dit que la solution u(¢, t ) du problème de Dirichlet pour EID se comporte comme ?t Åv(x)Åo(1) quand t !1, où v est la solution du problème ergodique stationaire qui correspond à une unique constante ergodique ?.En général, l'étude du comportement asymptotique est basé sur deux arguments: la régularité de solutions et le principe de maximumfort.Dans un premier temps, nous étudions le Principe de Maximum Fort pour les solutions de viscosité semicontinues des équations intégro-différentielles non-linéaires. Nous l'utilisons ensuite pour déduire un résultat de comparaison fort entre sous et sur-solutions des équations intégro-différentielles, qui va assurer l'unicité des solutions du problème ergodique à une constante additive près. De plus, pour des équationssuper-quadratiques le principe de maximum fort et en conséquence le comportement en temps grand exige la régularité Lipschitzienne.Dans une deuxième partie, nous établissons de nouvelles estimations Hölderiennes et Lipschitziennes pour les solutions de viscosité d'une large classe d'équations intégro-différentielles non-linéaires, par la méthode classique de Ishii-Lions. Les résultats de régularité aident de plus à la résolution du problème ergodique et sont utilisés pour fournir existence des solutions périodiques des EID.Nos résultats s'appliquent à une nouvelle classe d'équations non-locales que nous appelons équations intégro-différentielles mixtes. Ces équations sont particulièrement intéressantes, car elles sont dégénérées à la fois dans le terme local et non-local, mais leur comportement global est conduit par l'interaction locale - non-locale, par exemple la diffusion fractionnaire peut donner l'ellipticité dans une direction et la diffusion classique dans la direction orthogonale.Visualisation et Restauration d'Images par Mouvements de CourbureLe rôle de la courbure dans la perception visuelle remonte à 1954, et on le doit à Attneave. Des arguments neurologiques expliquent que le cerveau humain ne pourrait pas possiblement utiliser toutes les informations fournies par des états de simulation. Mais en réalité on enregistre des régions où la couleur change brusquement (des contours) et en outre les angles et les extremas de courbure. Pourtant, un calcul direct de courbures sur une image est impossible. Nous montrons comment les courbures peuvent être précisément évaluées, à résolution sous-pixelique par un calcul sur les lignes de niveau après leur lissage indépendant.Pour cela, nous construisons un algorithme que nous appelons Level Lines (Affine) Shortening, simulant une évolution sous-pixelique d'une image par mouvement de courbure moyenne ou affine. Aussi bien dans le cadre analytique que numérique, LLS (respectivement LLAS) extrait toutes les lignes de niveau d'une image, lisse indépendamment et simultanément toutes ces lignes de niveau par Curve Shortening(CS) (respectivement Affine Shortening (AS)) et reconstruit une nouvelle image. Nousmontrons que LL(A)S calcule explicitement une solution de viscosité pour le le Mouvement de Courbure Moyenne (respectivement Mouvement par Courbure Affine), ce qui donne une équivalence avec le mouvement géométrique.Basé sur le raccourcissement de lignes de niveau simultané, nous fournissons un outil de visualisation précis des courbures d'une image, que nous appelons un Microscope de Courbure d'Image. En tant que application, nous donnons quelques exemples explicatifs de visualisation et restauration d'image : du bruit, des artefacts JPEG, de l'aliasing seront atténués par un mouvement de courbure sous-pixelique / The present dissertation has two independent parts.Viscosity solutions theory for nonlinear Integro-Differential EquationsWe consider nonlinear elliptic and parabolic Partial Integro-Differential Equations (PIDES), where the nonlocal terms are associated to jump Lévy processes. The present work is motivated by the study of the Long Time Behavior of Viscosity Solutions for Nonlocal PDEs, in the periodic setting. The typical result states that the solution u(¢, t ) of the initial value problem for parabolic PIDEs behaves like ?t Å v(x) Å o(1) as t ! 1, where v is a solution of the stationary ergodic problem corresponding to the unique ergodic constant ?. In general, the study of the asymptotic behavior relies on two main ingredients: regularity of solutions and the strong maximum principle.We first establish Strong Maximum Principle results for semi-continuous viscosity solutions of fully nonlinear PIDEs. This will be used to derive Strong Comparison results of viscosity sub and super-solutions, which ensure the up to constants uniqueness of solutions of the ergodic problem, and subsequently, the convergence result. Moreover, for super-quadratic equations the strong maximum principle and accordingly the large time behavior require Lipschitz regularity.We then give Lipschitz estimates of viscosity solutions for a large class of nonlocal equations, by the classical Ishii-Lions's method. Regularity results help in addition solving the ergodic problem and are used to provide existence of periodic solutions of PIDEs. In both cases, we deal with a new class of nonlocal equations that we term mixed integrodifferential equations. These equations are particularly interesting, as they are degenerate both in the local and nonlocal term, but their overall behavior is driven by the local-nonlocal interaction, e.g. the fractional diffusion may give the ellipticity in one direction and the classical diffusion in the complementary one.Image Visualization and Restoration by CurvatureMotionsThe role of curvatures in visual perception goes back to 1954 and is due to Attneave. It can be argued on neurological grounds that the human brain could not possible use all the information provided by states of simulation. But actually human brain registers regions where color changes abruptly (contours), and furthermore angles and peaks of curvature. Yet, a direct computation of curvatures on a raw image is impossible. We show how curvatures can be accurately estimated, at subpixel resolution, by a direct computation on level lines after their independent smoothing.To performthis programme, we build an image processing algorithm, termed Level Lines (Affine) Shortening, simulating a sub-pixel evolution of an image by mean curvature motion or by affine curvature motion. Both in the analytical and numerical framework, LL(A)S first extracts all the level lines of an image, then independently and simultaneously smooths all of its level lines by curve shortening (CS) (respectively affine shortening (AS)) and eventually reconstructs, at each time, a new image from the evolved level lines.We justify that the Level Lines Shortening computes explicitly a viscosity solution for the Mean CurvatureMotion and hence is equivalent with the clasical, geometric Curve Shortening.Based on simultaneous level lines shortening, we provide an accurate visualization tool of image curvatures, that we call an Image CurvatureMicroscope. As an application we give some illustrative examples of image visualization and restoration: noise, JPEG artifacts, and aliasing will be shown to be nicely smoothed out by the subpixel curvature motion.
|
69 |
Quelques problèmes liés à l'erreur statistique en homogénéisation stochastique / Some problems related to statistical error in stochastic homogenizationMinvielle, William 25 September 2015 (has links)
Le travail de cette thèse a porté sur le développement de techniques numériques pour l'homogénéisation d'équations dont les coefficients présentent des hétérogénéités aléatoires à petite échelle. Les difficultés liées à la résolution de telles équations aux dérivées partielles peuvent être résolues grâce à la théorie de l'homogénéisation stochastique. On substitue alors la résolution d'une équation dont les coefficients sont aléatoires et oscillants à l'échelle la plus fine du problème par la résolution d'une équation à coefficients constants. Cependant, une difficulté subsiste : le calcul de ces coefficients dits homogénéisés sont définis par une moyenne ergodique, que l'on ne peut atteindre en pratique. Seuls des approximations aléatoires de ces quantités déterministes sont calculables, et l'erreur commise lors de l'approximation est importante. Ces questions sont développées en détail dans le Chapitre 1 qui tient lieu d'introduction. L'objet du Chapitre 2 de cette thèse est de réduire l'erreur de cette approximation dans un cas nonlinéaire, en réduisant la variance de l'estimateur par la méthode des variables antithétiques. Dans le Chapitre 3, on montre comment obtenir une meilleure réduction de variance par la méthode des vari- ables de contrôle. Cette approche repose sur un modèle approché, disponible dans le cas étudié. Elle est plus invasive et moins générique, on l'étudie dans un cas linéaire. Dans le Chapitre 4, à nouveau dans un cas linéaire, on introduit une méthode de sélection pour réduire l'erreur commise. Enfin, le Chapitre 5 porte sur l'analyse d'un problème in- verse, où l'on recherche des paramètres à l'échelle la plus fine, ne connaissant que quelques quantités macroscopiques, par exemple les coefficients homogénéisés du modèle / In this thesis, we design numerical techniques to address the homogenization of equations the coefficients of which exhibit small scale random heterogeneities. Solving such elliptic partial differential equations is prohibitively expensive. One may use stochastic homogenization theory to reduce the complexity of this task. We then substitute the random, fine scale oscillating coefficients of the equation with constant homogenized coefficients. These coefficients are defined through an ergodic average inaccessible to practical computation. Only random approximations thereof are available. The error committed in this approximation is significant. These issues are detailed in the introductory Chapter 1. In Chapter 2, we show how to reduce the error in this approximation, in a nonlinear case, by using an antithetic variable estimator that has a smaller variance than the standard Monte Carlo estimator. In Chapter 3, in a linear case, we show how to obtain an even better variance reduction with the control variate method. Such a method is based on a surrogate model. In Chapter 4, we use a selection method to reduce the global error. Chapter 5 is devoted to the analysis of an inverse problem, wherein we seek parameters at the fine scale whilst only being provided with a handful of macroscopic quantities, among which the homogenized coefficients
|
70 |
Analyse asymptotique d'équations intégro-différentielles : modèles d'évolution et de dynamique des populations / Asymptotic Analysis of Integro-differential Equations : populations dynamics and evolutionary modelsPatout, Florian 27 September 2019 (has links)
Cette thèse est consacrée à l’étude de phénomènes de propagation et de concentration dans des modèles d’équations intégro-différentielles venant de la écologie. On étudie certaines équations de réaction-diffusion non locales apparaissant en dynamique de populations, ainsi que des modèles représentant l’évolution Darwinienne avec un mode de reproduction sexué.Dans une première partie, nous étudions la propagation spatiale pour une équation de réaction-diffusion ou la dispersion opère via un noyau de convolution à queue lourde. Nous mesurons de manière précise l’accélération du front de propagation de la solution. Nous proposons également une échelle adaptée pour mesurer les «petites» mutations. Dans les deux cas nous utilisons le formalisme des équations de Hamilton-Jacobi.Dans un second temps nous étudions un modèle de génétique quantitative, avec un mode de reproduction sexuée. Un petit paramètre mesure la déviation entre le trait des descendants est la moyenne des traits des parents. Dans le régime où ce paramètre est petit nous étudions l’existence de solutions stationnaires, puis le problème de Cauchy lié à ce modèle. Les solutions se concentrent autour des optima de sélection, sous la forme de perturbations de distributions Gaussiennes avec petite variance fixée par le paramètre. Notre analyse généralise le cas linéaire de la reproduction asexuée en utilisant des outils d’analyse perturbative. Enfin dans une dernière partie nous fournissons des simulations numériques et des méthodes mathématiques pour étudier la dynamique interne des équilibres dans le régime de petite variance, pour les deux modes de reproduction : asexué et sexué. / This manuscript tackles propagation and concentration phenomena in different integro-differential equations with a background in ecology. We study non local reaction-diffusion equations from population dynamics, and models for Darwinian evolution with a sexual or asexual mode of reproduction, with a preference for the former.In a first part, we study spatial propagation for a reaction diffusion equation where dispersion acts through a fat tailed kernel. We measure accurately the acceleration of the propagation front of the population. We propose as well a scaling well adapted to “small mutations” when we consider the model in the context of adaptative dynamics. This scaling is very natural following the previous spatial investigation. In both cases we look at the long time behavior and we use the Hamilton-Jacobi framework. Then we turn our attention towards a quantitative genetics model, with a sexual mode of reproduction, imposed by the “infinitesimal operator”. In this non-linear setting, a small parameter tunes the deviation between the phenotypic trait of the offspring and the mean of the traits of the parents. In the regime where this parameter is small, we prove existence of stationary solutions, and their local uniqueness. We also provide an example of non-uniqueness in the case where the selection function admits several extrema. We prove that the solution concentrates around the points of minimum of the selection function. The analysis is carried by the small perturbations of special profiles : Gaussian distributions with small variance fixed by the parameter.We then study the stability of the Cauchy problem associated to the previous model. This time we prove that at all times, for a well prepared initial data, the solutions is arbitrary close to a Gaussian distribution with small variance. The proof follows the framework of the previous : we use perturbative analysis tools, but this time an even more precise description of the correctors is needed and we linearize the equation to obtain it. In a final part we show numerical simulations and different mathematical approaches to study inside dynamics of phenotypic lineages in the regime of small variance, with a moving environement.
|
Page generated in 0.1467 seconds