L'objet de cette thèse est d'étudier le comportement en temps long de solutions d'équations d'évolution non locales ainsi que la limite singulière d'équations et de systèmes d'équations aux dérivées partielles, où intervient un petit paramètre epsilon. Au Chapitre 1, nous considérons une équation de réaction-diffusion non locale avec conservation au cours du temps de l'intégrale en espace de la solution; cette équation a été initialement proposée par Rubinstein et Sternberg pour modéliser la séparation de phase dans un mélange binaire. Le problème de Neumann associé possède une fonctionnelle de Lyapunov, c'est-à-dire une fonctionnelle qui décroit selon les orbites. Après avoir prouvé que la solution est confinée dans une région invariante, nous étudions son comportement en temps long. Nous nous appuyons sur une inégalité de Lojasiewicz pour montrer qu'elle converge vers une solution stationnaire quand t tend vers l'infini. Nous évaluons également le taux de la convergence et calculons précisément la solution stationnaire limite en dimension un d'espace. Le Chapitre 2 est consacré à l'étude de l'équation différentielle non locale que l'on obtient en négligeant le terme de diffusion dans l'équation d'Allen-Cahn non locale étudiée au Chapitre 1. Sans le terme de diffusion, la solution ne peut pas être plus régulière que la fonction initiale. C'est la raison pour laquelle on ne peut pas appliquer la méthode du Chapitre 1 pour l'étude du comportement en temps long de la solution. Nous présentons une nouvelle méthode basée sur la théorie des réarrangements et sur l'étude du profil de la solution. Nous montrons que la solution est stable pour les temps grands et présentons une caractérisation détaillée de sa limite asymptotique quand t tend vers l'infini. Plus précisément, la fonction limite est une fonction en escalier, qui prend au plus deux valeurs, qui coïncident avec les points stables d'une équation différentielle associée. Nous montrons aussi par un contre-exemple non trivial que, quand une hypothèse sur la fonction initiale n'est pas satisfaite, la fonction limite peut prendre trois valeurs, qui correspondent aux points instable et stables de l'équation différentielle associée. Nous étudions au Chapitre 3 une équation différentielle ordinaire non locale qui a éte proposée par M. Nagayama. Une difficulté essentielle est que le dénominateur dans le terme de réaction non local peut s'annuler. Nous appliquons un théorème de point fixe lié a une application contractante pour démontrer que le problème à valeur initiale correspondant possède une solution unique qui reste connée dans un ensemble invariant. Ce problème possède une fonctionnelle de Lyapunov, qui est un ingrédient essentiel pour démontrer que la solution converge vers une solution stationnaire constante par morceaux quand t tend vers l'infini. Au Chapitre 4, nous considérons un modèle d'interface diffuse pour la croissance de tumeurs, où intervient une équation d'ordre quatre de type Cahn Hilliard. Après avoir introduit un modèle de champ de phase associé, on étudie formellement la limite singulière de la solution quand le coefficient du terme de réaction tend vers l'infini. Plus précisément, nous montrons que la solution converge vers la solution d'un problème à frontière libre. AMS subject classifications. 35K57, 35K50, 35K20, 35R35, 35R37, 35B40, 35B25.

Mes travaux de recherche sont liés à l'analyse asymptotique, l'approximation numérique et la théorie spectrale de problèmes elliptiques. J'allie les résultats théoriques et les simulations numériques pour préciser le comportement des solutions : la théorie permettant de proposer des méthodes numériques plus performantes et de prévoir certaines difficultés numériques, les simulations illustrant parfois des comportements plus fins que ceux démontrés jusque-là ou suggérant de nouvelles conjectures. Ce document se découpe en quatre chapitres, chacun correspondant à un thème de recherche. Le premier thème de recherche que je vais aborder concerne l'analyse mathématique de la supraconductivité qui était le sujet de ma thèse. Cette thématique a fait l'objet de collaborations avec F. Alouges, M. Dauge, S. Fournais, B. Helffer, D.~Martin, N. Popoff, N. Raymond et G. Vial. Notre objectif est de comprendre l'influence de la géométrie du matériau sur l'apparition de la supraconductivité. La première étape consiste à étudier le spectre de l'opérateur de Schrödinger avec champ magnétique et paramètre semi-classique dans les domaines à coins. Nous avons établi un développement asymptotique des modes propres et montré que les vecteurs propres avaient une structure double échelle, ce qui rend les simulations numériques très délicates. Nous avons proposé une approche basée sur la méthode d'éléments finis nodaux de haut degré et mis en évidence l'effet tunnel pour des domaines symétriques. Ces résultats ont ensuite permis d'étudier les minimiseurs de la fonctionnelle de Ginzburg-Landau et d'établir la localisation du paramètre d'ordre qui rend compte de la densité des électrons supraconducteurs, lorsqu'on abaisse progressivement le champ magnétique appliqué. Très peu d'études avaient été réalisées pour les domaines à coins. Nous en avons maintenant une compréhension assez précise en dimension 2. Récemment, nous avons commencé l'étude en dimension 3 dans le cadre de la thèse de N. Popoff, avec M. Dauge. La deuxième partie de ce document présente un modèle simplifié pour le transport quantique dans des diodes à effet tunnel résonant. Elle résulte de collaborations avec A. Faraj, F. Nier et Y. Patel. Ce dernier a réalisé une analyse asymptotique fine de systèmes de Schrödinger-Poisson stationnaires, non linéaires uni-dimensionnels dans un régime hors-équilibre. Nous avons proposé une adaptation numérique de cette analyse afin de déterminer rapidement des diagrammes courant-tension et de bifurcation et montré la pertinence de ce modèle réduit en le comparant à un modèle 1D de Schrödinger-Poisson avec traitement numérique complet des états résonnants. Dans le cadre du projet ANR jeunes chercheurs n° JCJC06-139561 Macadam, je me suis intéressée à l'analyse multi-échelle et numérique de problèmes elliptiques perturbés, en collaboration avec D. Brancherie, M. Dambrine, S. Tordeux, F. Hérau et G. Vial. Ce projet consiste à étudier l'influence de petites perturbations géométriques sur la solution de problèmes elliptiques. Les cas d'une inclusion isolée ou de plusieurs bien séparées ont été largement étudiés. Nous considérons plus précisément le cas où la distance entre deux inclusions tend vers zéro mais reste grande par rapport à leur taille caractéristique. Nous donnons un développement asymptotique multi-échelle complet de la solution de l'équation de Laplace dans la situation de deux inclusions. Nous présentons également quelques simulations numériques basées sur une méthode de superposition multi-échelle de la solution non perturbée et d'un profil (solution normalisée de l'équation de Laplace dans le domaine extérieur obtenu par blow-up de la perturbation). Nous étendons ces techniques aux équations de l'élasticité linéaire afin de prédire le comportement à rupture de certains matériaux présentant des micro-défauts. Nous avons également proposé des méthodes pour calculer effectivement les profils intervenant dans le développement asymptotique. Ceci a soulevé des questions mathématiques liées à la perte de coercivité provenant de conditions de Ventcel dégénérées. Le dernier chapitre propose quelques résultats sur les partitions minimales, en collaboration avec B. Helffer, T. Hoffmann-Ostenhof, C. Léna et G. Vial. Nous souhaitons comprendre le lien entre la $k$-partition minimale, pour laquelle la plus grande première valeur propre du Laplacien-Dirichlet sur les $k$ sous-domaines est minimale parmi les $k$-partitions, et les ensembles nodaux des vecteurs propres du Laplacien avec condition de Dirichlet. Nous nous sommes focalisés sur le cas $k=3$ pour lequel on ne connaît pas, en général, de partition optimale même pour des géométries très simples telles que le carré ou le disque. En se restreignant aux configurations symétriques, nous utilisons la méthode d'éléments finis pour exhiber des candidats aux 3-partitions minimales symétriques du disque, du carré ou d'autres géométries. Cette étude numérique nous a conduits à des problèmes d'isospectralité que nous avons résolus en utilisant le Hamiltonien de Aharonov-Bohm. L'introduction de cet opérateur pour résoudre une question théorique a ouvert une nouvelle piste numérique qui consiste à calculer les modes propres par une méthode d'éléments finis sur un revêtement à deux feuillets et d'étudier le comportement des lignes nodales en fonction du point singulier. Cela nous a permis de dégager de nouveaux candidats aux partitions minimales.

Le but de cette thèse est de développer une méthode numérique simple, robuste, efficace et précise pour résoudre des problèmes d'ingénierie de grande taille à partir de la méthode Taylor Meshless (TMM) et fournir de nouvelles idées principales de TMM est d'utiliser comme fonctions de forme des polynômes d'ordre élevé qui sont des solutions approchées de l'EDP. Ainsi la discrétisation ne concerne que la frontière. Les coefficients de ces fonctions de forme sont obtenus en discrétisant les conditions aux limites par des procédures de collocation associées à la méthode des moindres carrés. TMM est alors une véritable méthode sans maillage sans processus d'intégration, les conditions aux limites étant obtenues par collocation. Les principales contributions de cette thèse sont les suivantes: 1) Basé sur TMM, un algorithme général et efficace a été développé pour résoudre des EDP elliptiques tridimensionnelles; 2) Trois techniques de couplage pour des résolutions par morceaux ont été discutées dans des cas de problèmes à grande échelle: la méthode de collocation par les moindres carrés et deux méthodes de couplage basées sur les multiplicateurs de Lagrange; 3) Une méthode numérique générale pour résoudre les EDP non-linéaires a été proposée en combinant la méthode de Newton, la TMM et la technique de différentiation automatique. 4) Pour résoudre des problèmes avec un bord non régulier, des solutions singulières satisfaisant l'équation de contrôle sont introduites comme des fonctions de forme complémentaires, ce qui fournit une base théorique pour la résolution de problèmes singuliers / Based on Taylor Meshless Method (TMM), the aim of this thesis is to develop a simple, robust, efficient and accurate numerical method which is capable of solving large scale engineering problems and to provide a new idea for the follow-up study on meshless methods. To this end, the influence of the key factors in TMM has been studied by solving three-dimensional and non-linear Partial Differential Equations (PDEs). The main idea of TMM is to use high order polynomials as shape functions which are approximated solutions of the PDE and the discretization concerns only the boundary. To solve the unknown coefficients, boundary conditions are accounted by collocation procedures associated with least-square method. TMM that needs only boundary collocation without integration process, is a true meshless method. The main contributions of this thesis are as following: 1) Based on TMM, a general and efficient algorithm has been developed for solving three-dimensional PDEs; 2) Three coupling techniques in piecewise resolutions have been discussed and tested in cases of large-scale problems, including least-square collocation method and two coupling methods based on Lagrange multipliers; 3) A general numerical method for solving non-linear PDEs has been proposed by combining Newton Method, TMM and Automatic Differentiation technique; 4) To apply TMM for solving problems with singularities, the singular solutions satisfying the control equation are introduced as complementary shape functions, which provides a theoretical basis for solving singular problems

Les environnements CAO fournissent des outils puissants pour la programmation graphique et la manipulation d’entités géométriques complexes. Dans cette thèse, nous proposons d’exploiter ce potentiel dans le domaine de la conception de robots parallèles. Ces robots sont considérés comme une alternative intéressante vis-à-vis de leurs homologues sériels dans différentes applications comme le « pick and place » et l’usinage. Cependant, leur utilisation industrielle est encore restreinte en raison d’un espace de travail limité, de modèles géométriques difficiles à résoudre et l’existence de configurations singulières délimitant leur domaine d’exploitation. L’analyse et la caractérisation de l’espace de travail jouent alors un rôle fondamental dans la phase de conception de robots manipulateurs parallèles. Dans ce travail de thèse, nous proposons des approches géométriques originales donnant lieu à un ensemble de méthodes et techniques basées CAO pour l’analyse et la caractérisation de l’espace de travail de robots parallèles plans et spatiaux. L’espace de travail est généré comme un solide dans l’environnement CAO à partir d’un paramétrage géométrique, d’esquisses et d’opérations élémentaires telles que le balayage hélicoïdal et l’intersection. Nous avons montré que ces méthodes constituent des outils pertinents et efficaces d’aide à la conception des mécanismes parallèles. Ils permettent également la résolution du problème géométrique direct et la génération de trajectoires libres de singularités. Plusieurs types de manipulateurs ont été considérés dans ce travail pour mettre en avant et illustrer les techniques CAO / Géométriques proposées : robots parallèles plans à 3 degrés de mobilité de type 3-RPR, 3-RRR, 3-PPR et 3-PRR, robots parallèles spatiaux à 6 degrés de mobilité de type ou 3-CRS ou 3-PRRS. / CAD environments provide very powerful tools for graphical programming and manipulation of complex geometric entities. In this thesis, we propose to exploit such potential in the design of parallel robots. These robots are considered an attractive and important alternative towards their serials counterparts in various applications, like “pick and place” and machining. However, their industrial applications are restricted due to limited workspace, complexity related to resolution of the direct geometric model, and in addition the existence of the singular configurations which bound their application field. The analysis and the characterization of the workspace therefore play an essential role in the design phase of parallel robotic manipulators. In this thesis, we suggest original geometric approaches giving rise to a set of methodologies and techniques based on the use of CAD in order to analyze and characterize the workspace of planar and spatial parallel robotic manipulators. Workspace is generated as a solid in CAD environment by using a parametric geometric model, sketches, and elementary operations such as helical scanning and performing then Boolean intersection operation. We have shown in this thesis, that the proposed methodologies represent relevant and efficient tools which assist designers of parallel mechanisms. Moreover, they allow us to solve the direct geometric problem and to plan singularity-free trajectories. Several types of robotic manipulators have been considered in this work to highlight and illustrate the proposed CAD / Geometric techniques : planar parallel manipulators having three degree of freedom such as 3-RPR, 3-RRR, 3-PPR, and 3-PRR, and spatial parallel robotic manipulators having six degree of freedom 3-CRS-type.

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille / This thesis aims to develop a method for numerical modelling of fluid flow through fractured porous media and for determination of their effective permeability by taking advantage of recent results based on formulation of the problem by Singular Integral Equations. In parallel, it was also an occasion to continue on the theoretical development and to obtain new results in this area. The governing equations for flow in such materials are reviewed first and mass conservation at the fracture intersections is expressed explicitly. Using the theory of potential, the general potential solutions are proposed in the form of a singular integral equation that describes the steady-state flow in and around several fractures embedded in an infinite porous matrix under a far-field pressure condition [136, 139]. These solutions represent the pressure field in the whole body as functions of the infiltration in the fractures, which fully take into account the fracture interaction and intersections. Closed-form solutions for the fundamental problem of fluid flow around a single fracture are derived, which are considered as the benchmark problems to validate the numerical solutions. In particular, the solution obtained for the case of an elliptical disc-shaped crack obeying to the Poiseuille's law has been compared to that obtained for ellipsoidal inclusions with Darcy's law [140].The numerical programs have been developed based on the singular integral equations method to resolve the general potential equations [132, 180]. These allow modeling the fluid flow through a porous medium containing a great number of fractures. Besides, this formulation of the problem also allows obtaining a semi-analytical infiltration solution over a single fracture depending on the matrice permeability, the fracture conductivity and the fracture geometry. This result is the important key to upscalling the effective permeability of a fractured porous medium by using different homogeneisation schemes. The results obtained by the self-consistent scheme have been in particular established. The multi-region approach can be used to extend the general potential solution written for the infinite domain to that for a finite domain [181]. A closed-form solution for flow in and around a single partially saturated fracture, surrounded by an infinite matrix subjected to a far-field condition, is also derived combining the solutions for a superconductive fracture and for an imprevious fracture. This solution is then employed to estimate the effective permeability of unsaturated fractured porous media [141].The effective permeability model is applied to study the hydromechanical behaviour of a fault zone constituted by a clay core surrounded by fractured zones in the context of CO2 geological storage. The pressure injection induces an overpressure in the reservoir that may affect the permeability of the fractured zones leading to complexe coupled hydromechanical phenomena. The simulation results allow evaluating the risk of leakage of the reservoir brine to higher aquifers as well as the risk of fault reactivation

Cette thèse porte sur l'étude théorique et numérique des méthodes rapides pour résoudre les équations de l'élastodynamique 3D en domaine fréquentiel, et se place dans le cadre d'une collaboration avec la société Shell en vue d'optimiser la convergence des problèmes d'inversion sismique. La méthode repose sur l'utilisation des éléments finis de frontière (BEM) pour la discrétisation et sur les techniques de matrices hiérarchiques (H-matrices) pour l'accélération de la résolution du système linéaire. Dans le cadre de cette thèse on a développé un solveur direct pour les BEMs en utilisant une factorisation LU et un stockage hiérarchique. Si le concept des H-matrices est simple à comprendre, sa mise en oeuvre requiert des développements algorithmiques importants tels que la gestion de la multiplication de matrices représentées par des structures différentes (compressées ou non) qui ne comprend pas mois de 27 sous-cas. Un autre point délicat est l'utilisation des méthodes d'approximations par matrices compressées (de rang faible) dans le cadre des problèmes vectoriels. Une étude algorithmique a donc été faite pour mettre en oeuvre la méthode des H-matrices. Nous avons par ailleurs estimé théoriquement le rang faible attendu pour les noyaux oscillants, ce qui constitue une nouveauté, et montré que la méthode est utilisable en élastodynamique. En outre on a étudié l'influence des divers paramètres de la méthode en acoustique et en élastodynamique 3D, à fin de calibrer leur valeurs numériques optimales. Dans le cadre de la collaboration avec Shell, un cas test spécifique a été étudié. Il s'agit d'un problème de propagation d'une onde sismique dans un demi-espace élastique soumis à une force ponctuelle en surface. Enfin le solveur direct développé a été intégré au code COFFEE développé a POEMS (environ 25000 lignes en Fortran 90) / This thesis focuses on the theoretical and numerical study of fast methods to solve the equations of 3D elastodynamics in frequency-domain. We use the Boundary Element Method (BEM) as discretization technique, in association with the hierarchical matrices (H-matrices) technique for the fast solution of the resulting linear system. The BEM is based on a boundary integral formulation which requires the discretization of the only domain boundaries. Thus, this method is well suited to treat seismic wave propagation problems. A major drawback of classical BEM is that it results in dense matrices, which leads to high memory requirement (O (N 2 ), if N is the number of degrees of freedom) and computational costs.Therefore, the simulation of realistic problems is limited by the number of degrees of freedom. Several fast BEMs have been developed to improve the computational efficiency. We propose a fast H-matrix based direct BEM solver.

La thermographie infrarouge est une méthode largement employée pour la caractérisation des propriétés thermophysiques des matériaux. L’avènement des diodes laser pratiques, peu onéreuses et aux multiples caractéristiques, étendent les possibilités métrologiques des caméras infrarouges et mettent à disposition un ensemble de nouveaux outils puissants pour la caractérisation thermique et le contrôle non desturctif. Cependant, un lot de nouvelles difficultés doit être surmonté, comme le traitement d’une grande quantité de données bruitées et la faible sensibilité de ces données aux paramètres recherchés. Cela oblige de revisiter les méthodes de traitement du signal existantes, d’adopter de nouveaux outils mathématiques sophistiqués pour la compression de données et le traitement d’informations pertinentes. Les nouvelles stratégies consistent à utiliser des transformations orthogonales du signal comme outils de compression préalable de données, de réduction et maîtrise du bruit de mesure. L’analyse de sensibilité, basée sur l’étude locale des corrélations entre les dérivées partielles du signal expérimental, complète ces nouvelles approches. L'analogie avec la théorie dans l'espace de Fourier a permis d'apporter de nouveaux éléments de réponse pour mieux cerner la «physique» des approches modales.La réponse au point source impulsionnel a été revisitée de manière numérique et expérimentale. En utilisant la séparabilité des champs de température nous avons proposé une nouvelle méthode d'inversion basée sur une double décomposition en valeurs singulières du signal expérimental. Cette méthode par rapport aux précédentes, permet de tenir compte de la diffusion bi ou tridimensionnelle et offre ainsi une meilleure exploitation du contenu spatial des images infrarouges. Des exemples numériques et expérimentaux nous ont permis de valider dans une première approche cette nouvelle méthode d'estimation pour la caractérisation de diffusivités thermiques longitudinales. Des applications dans le domaine du contrôle non destructif des matériaux sont également proposées. Une ancienne problématique qui consiste à retrouver les champs de température initiaux à partir de données bruitées a été abordée sous un nouveau jour. La nécessité de connaitre les diffusivités thermiques du matériau orthotrope et la prise en compte des transferts souvent tridimensionnels sont complexes à gérer. L'application de la double décomposition en valeurs singulières a permis d'obtenir des résultats intéressants compte tenu de la simplicité de la méthode. En effet, les méthodes modales sont basées sur des approches statistiques de traitement d'une grande quantité de données, censément plus robustes quant au bruit de mesure, comme cela a pu être observé. / Infrared thermography is a widely used method for characterization of thermophysical properties of materials. The advent of the laser diodes, which are handy, inexpensive, with a broad spectrum of characteristics, extend metrological possibilities of infrared cameras and provide a combination of new powerful tools for thermal characterization and non destructive evaluation. However, this new dynamic has also brought numerous difficulties that must be overcome, such as high volume noisy data processing and low sensitivity to estimated parameters of such data. This requires revisiting the existing methods of signal processing, adopting new sophisticated mathematical tools for data compression and processing of relevant information.New strategies consist in using orthogonal transforms of the signal as a prior data compression tools, which allow noise reduction and control over it. Correlation analysis, based on the local cerrelation study between partial derivatives of the experimental signal, completes these new strategies. A theoretical analogy in Fourier space has been performed in order to better understand the «physical» meaning of modal approaches.The response to the instantaneous point source of heat, has been revisited both numerically and experimentally. By using separable temperature fields, a new inversion technique based on a double singular value decomposition of experimental signal has been introduced. In comparison with previous methods, it takes into account two or three-dimensional heat diffusion and therefore offers a better exploitation of the spatial content of infrared images. Numerical and experimental examples have allowed us to validate in the first approach our new estimation method of longitudinal thermal diffusivities. Non destructive testing applications based on the new technique have also been introduced.An old issue, which consists in determining the initial temperature field from noisy data, has been approached in a new light. The necessity to know the thermal diffusivities of an orthotropic medium and the need to take into account often three-dimensional heat transfer, are complicated issues. The implementation of the double singular value decomposition allowed us to achieve interesting results according to its ease of use. Indeed, modal approaches are statistical methods based on high volume data processing, supposedly robust as to the measurement noise.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

Ce mémoire porte sur l'étude d'un système de gestion d'énergie électrique dans un système multi-sources soumis à des perturbations exogènes. L'application visée est l'alimentation d'une propulsion hybride diesel/électrique équipée d'un système d'absorption des pulsations de couple. Les perturbations exogènes considérées peuvent être transitoires ou persistantes. Une perturbation transitoire correspond à une variation rapide du couple de charge, due par exemple à une accélération ou une décélération du véhicule. Une perturbation persistante provient du système de compensation des pulsations de couple générées par le moteur thermique. Le premier objectif du contrôle est de maintenir constante la tension du bus continu. Le deuxième objectif est d'absorber dans un système de stockage rapide constitué de super condensateur ces perturbations qui peuvent à terme provoquer une usure prématurée de la batterie. Le troisième objectif est de compenser l'auto-décharge dans le super condensateur en maintenant constante sa tension nominale. Les deux sources (batterie et super condensateur) sont reliées au bus continu par l'intermédiaire de deux convertisseurs boost DC/DC. La commande consiste à piloter les rapports cycliques de chaque convertisseur. C'est un système non linéaire où la commande est multiplicative de l'état. L'approche classique consistant à résoudre les équations Francis-Byrnes-Isidori ne s'applique pas directement dans ce cas où la sortie et la matrice d'interconnection dépendent de la commande. De plus, si cette approche est bien adaptée au rejet de perturbations persistantes, elle montre ces limites pour le rejet de perturbations non persistantes combiné à des objectifs de régulation. Notre approche a consisté à écrire le système sous un formalisme Port-Controlled Hamiltonian et à s'affranchir de la contrainte de la dépendance de la matrice d'interconnection avec la commande en utilisant la théorie des perturbations singulières. La commande du système dégénéré peut ensuite être calculée par une approche passive. Les performances de cette commande ont été testées en simulation et à l'aide d'un banc d'essai expérimental. Les résultats montrent l'efficacité du système d'absorption des différents types de perturbation tout en respectant les deux objectifs de régulation. / This thesis presents the research of energy management in a battery/ultracapacitor hybrid energy storage system with exogenous disturbance in hybrid electric vehicular application. Transient and harmonic persistent disturbances are the two kinds of disturbances considered in this thesis. The former is due to the transient load power demand during acceleration and deceleration, and the latter is introduced from the process of the internal combustion engine torque ripples compensation. Our control objective is to absorb the disturbances causing battery wear via the ultracapacitor, and meanwhile, to maintain a constant DC voltage and to compensate the self-discharge in the ultracapacitor to maintain it operating at the nominal state of charge. The object system is nonlinear due to the multiplicative relation between the input and the state. The traditional approach to solve Francis-Byrnes-Isidori equations cannot be directly applied in this case since the interconnect matrix depends on the control input. Besides, even if this approach is well suited to the rejection of persistent disturbances, it shows the limits for the case of non-persistent disturbances which is also our object. Our contributed control method is realized through a cascade control structure based on the singular perturbation theory. The ultracapacitor current with the fastest motion rate is controlled in the inner fast loop through which we impose the desired dynamic to the system. The reduced system controlled in the outer slow loop is a Hamiltonian system and the controller is designed via interconnection and damping assignment. Simulations and experiments have been carried out to evaluate the control performance. A contrast of the system responses with and without the control algorithm shows that, with the control algorithm, the ultracapacitor effectively absorbs the disturbances; and verifies the effectiveness of the control algorithm.

L'objet de cette thèse est d'étudier le comportement en temps long de solutions d'équations d'évolution non locales ainsi que la limite singulière d'équations et de systèmes d'équations aux dérivées partielles, où intervient un petit paramètre epsilon. Au Chapitre 1, nous considérons une équation de réaction-diffusion non locale avec conservation au cours du temps de l'intégrale en espace de la solution; cette équation a été initialement proposée par Rubinstein et Sternberg pour modéliser la séparation de phase dans un mélange binaire. Le problème de Neumann associé possède une fonctionnelle de Lyapunov, c'est-à-dire une fonctionnelle qui décroit selon les orbites. Après avoir prouvé que la solution est confinée dans une région invariante, nous étudions son comportement en temps long. Nous nous appuyons sur une inégalité de Lojasiewicz pour montrer qu'elle converge vers une solution stationnaire quand t tend vers l'infini. Nous évaluons également le taux de la convergence et calculons précisément la solution stationnaire limite en dimension un d'espace. Le Chapitre 2 est consacré à l'étude de l'équation différentielle non locale que l'on obtient en négligeant le terme de diffusion dans l'équation d'Allen-Cahn non locale étudiée au Chapitre 1. Sans le terme de diffusion, la solution ne peut pas être plus régulière que la fonction initiale. C'est la raison pour laquelle on ne peut pas appliquer la méthode du Chapitre 1 pour l'étude du comportement en temps long de la solution. Nous présentons une nouvelle méthode basée sur la théorie des réarrangements et sur l'étude du profil de la solution. Nous montrons que la solution est stable pour les temps grands et présentons une caractérisation détaillée de sa limite asymptotique quand t tend vers l'infini. Plus précisément, la fonction limite est une fonction en escalier, qui prend au plus deux valeurs, qui coïncident avec les points stables d'une équation différentielle associée. Nous montrons aussi par un contre-exemple non trivial que, quand une hypothèse sur la fonction initiale n'est pas satisfaite, la fonction limite peut prendre trois valeurs, qui correspondent aux points instable et stables de l'équation différentielle associée. Nous étudions au Chapitre 3 une équation différentielle ordinaire non locale qui a éte proposée par M. Nagayama. Une difficulté essentielle est que le dénominateur dans le terme de réaction non local peut s'annuler. Nous appliquons un théorème de point fixe lié a une application contractante pour démontrer que le problème à valeur initiale correspondant possède une solution unique qui reste connée dans un ensemble invariant. Ce problème possède une fonctionnelle de Lyapunov, qui est un ingrédient essentiel pour démontrer que la solution converge vers une solution stationnaire constante par morceaux quand t tend vers l'infini. Au Chapitre 4, nous considérons un modèle d'interface diffuse pour la croissance de tumeurs, où intervient une équation d'ordre quatre de type Cahn Hilliard. Après avoir introduit un modèle de champ de phase associé, on étudie formellement la limite singulière de la solution quand le coefficient du terme de réaction tend vers l'infini. Plus précisément, nous montrons que la solution converge vers la solution d'un problème à frontière libre. AMS subject classifications. 35K57, 35K50, 35K20, 35R35, 35R37, 35B40, 35B25. / The aim of this thesis is to study the large time behavior of solutions of nonlocal evolution equations and to also study the singular limit of equations and systems of parabolic partial differential equations involving a small parameter epsilon. In Chapter 1, we consider a nonlocal reaction-diffusion equation with mass conservation, which was originally proposed by Rubinstein and Sternberg as a model for phase separation in a binary mixture. The corresponding Neumann problem possesses a Lyapunov functional, namely a functional which decreases in time along solution orbits. After having proved that the solution is conned in an invariant region, we study its large time behavior and apply a Lojasiewicz inequality to show that it converges to a stationary solution as t tends to infinity. We also evaluate the rate of convergence and precisely compute the limiting stationary solution in one space dimension. Chapter 2 is devoted to the study of a nonlocal evolution equation which one obtains by neglecting the diffusion term in the nonlocal Allen-Cahn equation studied in Chapter 1. Without the diffusion term, the solution can not be expected to be more regular than the initial function. Moreover, because of the absence of the diusion term, the method of Chapter 1 can not be applied to study the large time behavior of the solution. We present a new method based up on rearrangement theory and the study of the solution profile. We show that the solution stabilizes for large times and give a detailed characterization of its asymptotic limit as t tends to infinity. More precisely, it turns out that the limiting function is a step function, which takes at most two values, which are stable points of a corresponding ordinary dierential equation. We also show by means of a nontrivial counterexample that, when a certain hypothesis on the initial function does not hold, the limiting function may take three values. One of them is the unstable point and the two others are the stable points of the ordinary dierential equation. We study in Chapter 3 a nonlocal ordinary dierential equation which has been proposed by M. Nagayama. The nonlocal term involves a denominator which may vanish. We apply a contraction fixed point theorem to prove the existence of a unique solution which stays confined in an invariant region. We also show that the corresponding initial value problem possesses a Lyapunov functional and prove that the solution stabilizes for large times to a step function, which takes at most two values. In Chapter 4, we consider a diffuse-interface tumor-growth model which involves a fourth order Cahn-Hilliard type equation. Introducing a related phase-field model, we formally study the singular limit of the solution as the reaction coecient tends to infinity. More precisely, we show that the solution converges to the solution of a moving boundary problem. AMS subject classifications. 35K57, 35K50, 35K20, 35R35, 35R37, 35B40, 35B25.