De récents résultats sur l’étude des systèmes de Lur’e (commutés) à temps discret mettent en avant une fonction de Lyapunov de type Lur’e avancée, dont les lignes de niveau peuvent être non convexes et non connexes. Celles-ci soulèvent de larges questions pour les systèmes de Lur’e à temps discret obtenus par la discrétisation d’un système continu. Les contributions de cette thèse sont d’apporter des éléments de réponse à ces questions. Tout d’abord, le verrou des lignes de niveau non-connexes est levé en construisant à partir de celles-ci une suite décroissante d’ensembles connexes et bornés qui converge vers l’origine et qui contient le futur de la trajectoire à temps continu. Dans un second temps, le problème de la stabilisation conjointe d’un système de Lur’e à données échantillonnées avec un échantillonnage non-uniforme est traité. Quand la période d’échantillonnage est à choisir parmi un nombre fini de valeurs, il est montré que ce problème se traduit comme la stabilisation conjointe d’un système commuté de Lur’e avec des incertitudes bornées en norme. En associant de plus à chaque mode un critère quadratique, une stratégie de type min-switching permet de résoudre cette question à l’aide d’un problème d’optimisation sous contraintes LMI. Enfin, les propriétés de la stratégie de min-switching pour les systèmes de Lur’e commutés à temps discret sont étudiées. Une extension de la notion de consistance permet de prouver que cette stratégie est consistante vis-à-vis de majorants quadratiques modaux du critère de performance et ainsi de garantir l’intérêt de la stratégie d’échantillonnage non-uniforme développée / Recent studies dealing with discrete-time (switched) Lur’e systems involve an adapted Lur’e type function exhibiting possibly non-convex and disconnected level sets. These properties raise fundamental issues in the case of discrete-time Lur’e system obtained by the sampling of a continuous time one. This PhD thesis aims at answering these questions. The first contribution is to avoid the discrete-time disconnected level sets by a decreasing sequence of bounded and connected sets that converges to the origin and that contain the future of the continuous-time trajectory. The second contribution deals with the joint stabilization of a sampled-data Lur’e system with non-uniform sampling. When the sampling period belongs to a finite set of values, this problem is reformulated as the joint stabilization of a discrete-time Lur’e switched system with norm-bounded uncertain parameters. Futhermore, if a quadratic criterion is associated with each mode, a min-switching strategy combined with LMI constraints allow to provide a solution to this problem. Finally the property of consistency for discrete-time switched Lur’e systems is investigated. It is shown that the min-switching strategy is consistent with respect to quadratic upper bounds of the performances. This result is applied on the stabilization of Lur’e systems with non-uniform sampling.

Ce travail s'inscrit dans une thématique de recherche sur l'étude des opérateurs pseudo-différentiels sous représentations diffusives ; l'intégration fractionnaire est un exemple devenu classique d'opérateurs diffusifs.<br />Le première partie consiste en la mise en place des représentations diffusives à temps discret. Certains filtres non-relationnels, notamment les différences frationnaires, sont une agrégation continue de dynamiques purement amorties. Les représentations diffusives s'appliquent à toutes les discrétisations de l'intégration fractionnaire y compris celles pour lesquelles la fonction de transfert n'est pas connue analytiquement. Les filtres diffusifs peuvent être réalisés par un système de dimension infinie. Cette structure est un cadre adapté à l'approximation par un filtre relationnel, à l'analyse asymptotique aux temps longs et à l'élaboration d'un critère de dissipativité.<br />La deuxième partie consiste à appliquer ces outils pour l'étude des couplages formés de filtres diffusifs et de filtres rationnels positifs. L'application d'un critère de Nyquist prouve la stabilité énergétique. Ces couplages sont en fait la somme d'une partie entière et d'une partie diffusive, ce résultat de décomposition montre que certains couplages sont stables EBSB (entrée-bornée, sortie-bornée). La dissipativité de la réalisation diffusive ainsi que le lemme de Kalman-Yacubovich-Popov montrent notamment la stabilité interne de ces couplages ; une démonstration originale du caractère asymptotique de la stabilité interne est ainsi proposée. Les approches utilisées pour prouver ces stabiblités permettent une analyse asymptotique aux temps longs.

A l'intérieur d'un solide, l'endommagement est provoqué par des mouvements microscopiques. Nous décidons de tenir compte de ces mouvements dans le principe des puissances virtuelles. Les vitesses microscopiques étant reliées à la vitesse d'endommagement, la puissance des efforts intérieurs que nous retenons dépend de cette vitesse et aussi de la vitesse du gradient de l'endommagement pour décrire les interactions. L'utilisation de cette approche permet de représenter le comportement non linéaire du béton, qui est en grande partie dû à l'endommagement, par des modèles adoucissants. Les applications aux calculs de structures donnent des résultats objectifs surmontant la difficulté liée à la dépendance du maillage lors de la discrétisation par éléments finis. La prise en compte du gradient de l'endommagement dans la formulation permet également de rendre compte de l'effet d'échelle de structures. Le phénomène unilatéral lié à la refermeture des microfissures, qui conduit à la restauration de la raideur lors du passage d'un chargement endommageable en traction à un chargement en compression, est pris en compte par l'utilisation de deux variables scalaires d'endommagement et par une partition du tenseur des déformations dans l'expression de l'énergie libre. Les modèles que nous présentons sont implantés dans un code de calcul par éléments finis et les résultats obtenus permettent une analyse prévisionnelle des structures.

Le calcul par éléments finis des structures en élastoplasticité classique conduit à des dépenses informatiques très élevées, sans qu'on puisse estimer de façon correcte les erreurs commises. Dans notre travail, nous proposons une nouvelle méthode d'estimation a priori des erreurs. Nous construisons d'abord les solutions exactes de problèmes mécaniques par une méthode inverse. Ensuite, nous utilisons ces solutions exactes pour évaluer l'erreur commise lors d'une analyse numérique, avec des discrétisations spatiale et temporelle données. Enfin, nous utilisons cette base d'exemples de deux manières: 1) la recherche de règles simples de remaillage quasi optimal; 2) la construction, par des techniques d'apprentissage automatique, de règles générales donnant a priori la discrétisation nécessaire pour respecter une erreur fixée.

La simulation numérique est aujourd'hui un outils majeur dans la conception des objets aérodynamiques, que ce soit dans l'aéronautique, l'automobile, l'industrie navale, etc... Un des défis majeurs pour repousser les limites des codes de simulation est d'améliorer leur précision, tout en utilisant une quantité fixe de ressources (puissance et/ou temps de calcul). Cet objectif peut être atteint par deux approches différentes, soit en construisant une discrétisation fournissant sur un maillage donné une solution d'ordre très élevé, soit en construisant un schéma compact et massivement parallélisable, de manière à minimiser le temps de calcul en distribuant le problème sur un grand nombre de processeurs. Dans cette thèse, nous tentons de rassembler ces deux approches par le développement et l'implémentation de Schéma Distribuant le Résidu (RDS) d'ordre très élevé et de compacité maximale. Ce manuscrit commence par un rappel des principaux résultats mathématiques concernant les Lois de Conservation hyperboliques (CLs). Le but de cette première partie est de mettre en évidence les propriétés des solutions analytiques que nous cherchons à approcher, de manière à injecter ces propriétés dans celles de la solution discrète recherchée. Nous décrivons ensuite les trois étapes principales de la construction d'un schéma RD d'ordre très élevé : \begin{itemize} \item la représentation polynomiale d'ordre très élevé de la solution sur des polygones et des polyèdres; \item la description de méthodes distribuant le résidu de faible ordre, compactes et conservatives, consistantes avec une représentation polynomiale des données de très haut degré. Parmi elles, une attention particulière est donnée à la plus simple, issue d'une généralisation du schéma de Lax-Friedrichs (LxF); \item la mise en place d'une procédure préservant la positivité qui transforme tout schéma stable et linéaire, en un schéma non linéaire d'ordre très élevé, capturant les chocs de manière non oscillante. \end{itemize} Dans le manuscrit, nous montrons que les schémas obtenus par cette procédure sont consistants avec la CL considérée, qu'ils sont stables en norme $\L^{\infty}$ et qu'ils ont la bonne erreur de troncature. Même si tous ces développements théoriques ne sont démontrés que dans le cas de CL scalaires, des remarques au sujet des problèmes vectoriels sont faites dès que cela est possible. Malheureusement, lorsqu'on considère le schéma LxF, le problème algébrique non linéaire associé à la recherche de la solution stationnaire est en général mal posé. En particulier, on observe l'apparition de modes parasites de haute fréquence dans les régions de faible gradient. Ceux-ci sont éliminés grâce à un terme supplémentaire de stabilisation dont les effets et l'évaluation numérique sont précisément détaillés. Enfin, nous nous intéressons à une discrétisation correcte des conditions limites pour le schéma d'ordre élevé proposé. Cette théorie est ensuite illustrée sur des cas test scalaires bidimensionnels simples. Afin de montrer la généralité de notre approche, des maillages composés uniquement de triangles et des maillages hybrides, composés de triangles et de quadrangles, sont utilisés. Les résultats obtenus par ces tests confirment ce qui est attendu par la théorie et mettent en avant certains avantages des maillages hybrides. Nous considérons ensuite des solutions bidimensionnelles des équations d'Euler de la dynamique des gaz. Les résultats sont assez bons, mais on perd les pentes de convergence attendues dès que des conditions limite de paroi sont utilisées. Ce problème nécessite encore d'être étudié. Nous présentons alors l'implémentation parallèle du schéma. Celle-ci est analysée et illustrée à travers des cas test tridimensionnel de grande taille. Du fait de la relative nouveauté et de la complexité des problèmes tridimensionels, seuls des remarques qualitatives sont faites pour ces cas test : le comportement global semble être bon, mais plus de travail est encore nécessaire pour définir les propriétés du schémas en trois dimensions. Enfin, nous présentons une extension possible du schéma aux équations de Navier-Stokes dans laquelle les termes visqueux sont traités par une formulation de type Galerkin. La consistance de cette formulation avec les équations de Navier-Stokes est démontrée et quelques remarques au sujet de la précision du schéma sont soulevées. La méthode est validé sur une couche limite de Blasius pour laquelle nous obtenons des résultats satisfaisants. Ce travail offre une meilleure compréhension des propriétés générales des schémas RD d'ordre très élevé et soulève de nouvelles questions pour des améliorations futures. Ces améliorations devrait faire des schémas RD une alternative attractive aux discrétisations classiques FV ou ENO/WENO, aussi bien qu'aux schémas Galerkin Discontinu d'ordre très élevé, de plus en plus populaires.

Dans une démarche de modélisation hydrogéologique, la méthode des éléments analytiques présente une alternative intéressante par rapport aux méthodes numériques usuellement employées pour l'hydrodynamique souterraine. C'est une méthode aux frontières qui a été développée spécifiquement pour la modélisation des écoulements souterrains. Après avoir identifié sa place parmi les autres méthodes existantes et discuté de ses avantages et inconvénients, il s'agit d'approfondir cette comparaison sur un cas d'étude réel. C'est pourquoi, cette méthode est comparée à la méthode des différences finies à travers une méthodologie envisageant toutes les étapes de la modélisation ; de la création du modèle conceptuel jusqu'aux résultats de simulation. Cette comparaison fait apparaître l'importance de l'étape de la discrétisation et l'importance cruciale qu'elle représente pour l'utilisateur sur la qualité générale de la modélisation. Après avoir identifié les approximations inhérentes à la méthode des éléments analytiques, une analyse mathématique de problèmes simplifiés et localisés sur quelques éléments analytiques permet de quantifier ces erreurs et d'y extraire une analyse de leur comportement numérique. Cela permet d'introduire par la même occasion des éléments de théorie de la méthode des éléments analytiques. La représentation objet et vectorielle des éléments analytiques permet de bien valoriser les interactions possibles avec les SIG. En particulier, il s'agit d'élaborer des algorithmes de discrétisation adaptés aux spécificités et aux besoins de la méthode, mais aussi de mettre en valeur des critères d'évaluation de la qualité de cette discrétisation. La synthèse de ces outils et méthodes permet à l'utilisateur de disposer d'un outil d'aide à la discrétisation pour réduire au mieux ses effets sur la qualité de modélisation. Il s'agit ensuite d'étudier le comportement numérique des éléments analytiques adaptés à la simulation en régime transitoire. Pour cela, des éléments analytiques 1D sont d'abord développés pour faciliter cette étude. Leur comportement numérique reste similaire tout en réduisant les erreurs numériques qui sont induites par le calcul des éléments analytiques euxmême. Cette analyse permet de mettre en évidence une méthode de résolution adaptée au régime transitoire qui est ensuite transposée au problème 2D. Elle y expose en particulier une méthode pour estimer le champ d'influence de l'élément transitoire, utilisant un critère de discrétisation spatial qui doit être corrélé avec le pas de discrétisation temporel.

Dans une démarche de modélisation hydrogéologique, la méthode des éléments analytiques présente une alternative intéressante par rapport aux méthodes numériques usuellement employées pour l'hydrodynamique souterraine. C'est une méthode aux frontières qui a été développée spécifiquement pour la modélisation des écoulements souterrains. Après avoir identifié sa place parmi les autres méthodes existantes et discuté de ses avantages et inconvénients, il s'agit d'approfondir cette comparaison sur un cas d'étude réel. C'est pourquoi, cette méthode est comparée à la méthode des différences finies à travers une méthodologie envisageant toutes les étapes de la modélisation ; de la création du modèle conceptuel jusqu'aux résultats de simulation. Cette comparaison fait apparaître l'importance de l'étape de la discrétisation et l'importance cruciale qu'elle représente pour l'utilisateur sur la qualité générale de la modélisation. Après avoir identifié les approximations inhérentes à la méthode des éléments analytiques, une analyse mathématique de problèmes simplifiés et localisés sur quelques éléments analytiques permet de quantifier ces erreurs et d'y extraire une analyse de leur comportement numérique. Cela permet d'introduire par la même occasion des éléments de théorie de la méthode des éléments analytiques. La représentation objet et vectorielle des éléments analytiques permet de bien valoriser les interactions possibles avec les SIG. En particulier, il s'agit d'élaborer des algorithmes de discrétisation adaptés aux spécificités et aux besoins de la méthode, mais aussi de mettre en valeur des critères d'évaluation de la qualité de cette discrétisation. La synthèse de ces outils et méthodes permet à l'utilisateur de disposer d'un outil d'aide à la discrétisation pour réduire au mieux ses effets sur la qualité de modélisation. Il s'agit ensuite d'étudier le comportement numérique des éléments analytiques adaptés à la simulation en régime transitoire. Pour cela, des éléments analytiques 1D sont d'abord développés pour faciliter cette étude. Leur comportement numérique reste similaire tout en réduisant les erreurs numériques qui sont induites par le calcul des éléments analytiques euxmême. Cette analyse permet de mettre en évidence une méthode de résolution adaptée au régime transitoire qui est ensuite transposée au problème 2D. Elle y expose en particulier une méthode pour estimer le champ d'influence de l'élément transitoire, utilisant un critère de discrétisation spatial qui doit être corrélé avec le pas de discrétisation temporel.

Cette étude porte sur les apports de l'automatique dans les systèmes intégrés sur Silicium pour la synthèse de radio-fréquences de l'ordre du Giga-Hertz avec une pureté spectrale optimale. La base des architectures qui réalisent cette fonction repose sur celle des boucles à verrouillage de phase (PLL). Ce travail se positionne par rapport aux di érents modèles existants et apporte des améliorations pour l'étude du bruit de phase et de la robustesse de ces systèmes non linéaires avec entre autres un outil pour l'analyse de la stabilité et de la robustesse semi globale. Ces travaux ont permis l'optimisation d'une nouvelle architecture fractionnaire réalisée et testée. L'approche systémique propre à l'automatique a aussi permis d'améliorer les performances des architectures numériques grâce à un observateur et un contrôleur robuste de type Hinfini. Leur implémentation sur Silicium est rendue possible grâce à des algorithmes qui visent à réduire la surface du circuit.

Cette thèse se divise en deux parties indépendantes. Dans la première, nous considérons un modèle microscopique individu-centré pour décrire une population structurée par traits et âges. Nous étudions l'écologie de ce système (problèmes de dynamique de populations) dans une asymptotique de grandes populations. Sous certaines renormalisations, le processus microscopique converge par la solution à valeurs mesures d'une équation d'évolution déterministe. Un théorème central limite et les déviations exponentielles associées à cette convergence sont étudiés. Nous appliquons ensuite ces résultats pour établir des généralisations aux populations structurées par âge de modèles d'évolution tirés de la récente théorie des dynamiques adaptatives. Ces derniers modélisent l'évolution de la structure en traits sur des grandes échelles de temps et sous les hypothèses de mutations rares (éventuellement petites) et de grandes populations. Dans la seconde partie de la thèse, nous considérons des équations aux dérivées partielles de McKean-Vlasov et de Navier-Stokes 2D avec conditions initiales aléatoires. La loi des solutions, qui sont alors des variables aléatoires, est appelée solution statistique. En nous basant sur une approche probabiliste de ces équations aux dérivées partielles, nous proposons de nouvelles approximations particulaires stochastiques avec ondelettes pour les moments d'ordre 1 des solutions statistiques, et nous étudions leurs vitesses de convergence.

Dans cette thèse, nous étudions l'approximation et la stabilisation de certaines équations d'évolution, en utilisant la théorie des semi-groups et l'analyse spectrale. Cette thèse est divisée en deux parties principales. Dans la première partie, comme dans [3, 4], nous considérons l'approximation des équations d'évolution du deuxième ordre modélisant les vibrations de structures élastiques. Il est bien connu que le système approché par éléments finis ou différences finies n'est pas uniformément exponentiellement ou polynomialement stable par rapport au paramètre de discrétisation, même si le système continu a cette propriété. Dans la première partie, notre objectif est d'amortir les modes parasites à haute fréquence en introduisant des termes de viscosité numérique dans le schéma d'approximation. Avec ces termes de viscosité, nous montrons la décroissance exponentielle ou polynomiale du schéma discret lorsque le problème continu a une telle décroissance et quand le spectre de l'opérateur spatial associé au problème conservatif satisfait la condition du gap généralisée. En utilisant le Théorème de Trotter-Kato, nous montrons la convergence de la solution discrète vers la solution continue. Quelques exemples sont également présentés.