Le résumé contient des caractères spéciaux

Il existe peu de codes commerciaux pour la simulation numérique des écoulements de liquides polymères. Les difficultés proviennent des propriétés intrinsèques des polymères, qui sont des fluides viscoélastiques non-newtoniens. Ceci implique un couplage entre la viscoélasticité du liquide et l'écoulement, couplage quantifié par le nombre de Weissenberg. D'un point de vue numérique, la source du problème est la perte de convergence des algorithmes lorsque ce nombre devient trop élevé. Cette thèse porte sur le développement de schémas numériques robustes pour la simulation de ces écoulements en considérant principalement le modèle de Giesekus. Nous nous sommes d'abord intéressés au problème de Stokes et nous avons fait l'étude d'une méthode de Galerkin discontinue moins coûteuse et plus robuste que la méthode "Interior Penalty" classique. Nous avons fait une analyse a priori et a posteriori et nous avons mis en évidence les relations entre cette méthode dG et les éléments finis non-conformes. Les résultats théoriques obtenus ont été validés numériquement. Par la suite, nous avons considéré le modèle à trois champs de Giesekus. La vitesse et la pression sont approchées par éléments finis non-conformes tandis que l'équation constitutive est traitée à l'aide d'éléments finis discontinus et d'un schéma décentré de type Lesaint-Raviart. L'analyse de ces schémas dans le cas quadrangulaire et triangulaire a été faite pour le problème de Stokes sous-jacent. Ces schémas ont ensuite été implémentés dans la librairie C++ Concha. Nous avons effectué des comparaisons avec des données expérimentales mettant en évidence le bon comportement du modèle de Giesekus mais aussi avec le code commercial Polyflow et une solution semi-analytique afin de valider nos schémas numériques. Nous avons obtenu des simulations réalistes pour des nombres de Weissenberg élevés sur des cas-tests populaires : écoulement autour d'un cylindre, contractions 4:1 et 4:1:4

Dans cette thèse, nous analysons une méthode de Galerkine discontinue (GD) et deux estimateurs d'erreur a posteriori pour l'équation d'advection-diffusion-réaction linéaire et stationnaire avec diffusion hétérogène. La méthode GD considérée, la méthode SWIP, utilise des moyennes pondérées dont les poids dépendent de la diffusion. L'analyse a priori montre que la convergence est optimale en le pas du maillage et robuste par rapport aux hétérogénéités de la diffusion, ce qui est confirmé par les tests numériques. Les deux estimateurs d'erreur a posteriori sont obtenus par une analyse par résidus et contrôlent la (semi-)norme d'énergie de l'erreur. L'analyse d'efficacité locale montre que presque tous les estimateurs sont indépendants des hétérogénéités. Le deuxième estimateur d'erreur est plus précis que le premier, mais son coût de calcul est légèrement plus élevé. Cet estimateur est basé sur la construction d'un flux H(div)-conforme dans l'espace de Raviart-Thomas-Nédéléc.

La première partie est consacrée à la résolution de systèmes linéaires. Le chapitre 1 expose des résultats théoriques et numériques sur les méthodes proposées par Altman et précise le lien avec les méthodes de Krylov. Le chapitre 2 utilise des techniques d'extrapolation introduites par Brezinski pour obtenir une estimation du vecteur erreur. Plusieurs méthodes de projection sont retrouvées et de nouvelles procédures d'accélération données. Dans la deuxième partie, une nouvelle stratégie inspirée de la méthode de Cauchy-Barzilai-Borwein permet de définir de nouveaux schémas résolvant des problèmes de point fixe. Des résultats numériques sur un problème de bifurcation et un théorème de convergence sont donnés. Les chapitres 4, 5 et 6 sont consacrés à l'accélération de l'algorithme EM utilisé pour calculer des estimateurs du maximum de vraisemblance. Une classe de schémas itératifs basés sur la stratégie précédente est présentée, un théorème de convergence et une application à un problème de tomographie sont donnés. La dernière partie, fruit d'un projet du cemracs 2003, traite d'un problème issu de la physique des plasmas : l'amélioration des Codes Particles in Cell à l'aide d'une reconstruction de la densité basée sur une méthode d'ondelettes et sa validation numérique.

.

Un robot manipulateur parallèle est composé de deux solides, une plate-forme et une base, reliés par des chaînes cinématiques articulées et motorisées permettant leur mouvement relatif. Le problème du modèle géométrique direct des robots parallèles est essentiel. Il consiste a déterminer la position et l'orientation de la plate-forme par rapport a un repère de base lorsque les variables articulaires sont fixées. Ce problème revient a considérer le placement d'un solide quelconque sur six sphères de rayons et de centres fixés. Il présente deux facettes, une géométrique et l'autre algébrique, menant a diverses méthodes à mettre en oeuvre. Le modèle géométrique direct est connu pour admettre au plus quarante solutions complexes. Nous présentons une approche consistant a utiliser des capteurs additionnels pour simplifier et résoudre ce problème. Nous utilisons ou développons toute une gamme de méthodes symboliques et parfois numériques permettant la résolution dans ce contexte précis. Les résultats obtenus permettent de définir des architectures de robot ayant un nombre minimum de solutions pour la résolution du modèle géométrique direct, Voire même des solutions explicites ou uniques.

De nombreux systèmes embarqués implémentent des applications de traitement du signal, notamment lors de communications. Certains de ces traitements sont effectués par des filtres linéaires, qu'il est donc nécessaire de mettre en oeuvre numériquement sur ces cibles. Les systèmes embarqués sont sujets à diverses contraintes qu'il faut optimiser tout en conservant des systèmes fiables en terme de performance et de précision. L'arithmétique virgule fixe est généralement préférée à l'arithmétique flottante pour des systèmes embarqués de traitement du signal, entre autres car elle est moins coûteuse, disponibles dans tous les systèmes, permet d'utiliser des largeurs arbitraires sur des cibles matérielles et est généralement suffisante en terme de précision pour les applications de traitement du signal. Le calcul en virgule fixe nécessite d'aligner les positions des virgules pour ainsi rendre cohérent des calculs à base de nombres entiers. Cela implique des quantifications et l'enjeu est donc de minimiser la répercussion de ces arrondis sur le résultat final, en proposant une garantie sur l'erreur sur la sortie. Une méthodologie a été proposée durant cette thèse qui, à partir d'un algorithme de filtre linéaire, utilise une méthode analytique pour implémenter cet algorithme en virgule fixe et générer du code. Cette méthodologie considère à la fois les implémentations logicielles, et les implémentations matérielles qui impliquent la résolution d'un problème d'optimisation. Un outil, nommé FiPoGen, a été développé pour mettre en oeuvre les méthodes proposées et fournir automatiquement un code virgule fixe implémentant un filtre donné avec garantie sur l'erreur sur la sortie. / Embedded systems implement signal processing systems, such as linear filters, for example for communication through networks. These devices are subject to various constraints, such as power consumption, time-To-Market, area consumption, and so on, that is necessary to optimize while guaranteeing reliability and accuracy of the implemented systems. Fixed-Point arithmetic is generally used instead of floating-Point arithmetic for signal processing embedded systems because it is less expensive, all devices support fixed-Point numbers (as they are implemented only using integers), allows arbitrary word-Lengths in hardware implementation and is enough accurate for signal processing programs. Fixed-Point computations need the operands of an operation to be aligned together with the same position of the binary point. This leads to quantification errors and the goal is to minimize these round-Off effects onto the final result, by proposing a guarantee on the output error. During this thesis, a methodology has been proposed which implements a given algorithm in fixed-Point using analytical approach, and generates some codes. This methodology consider both software and hardware targets. The hardware approach is realized solving a word-Length optimization problem. A tool, named FiPoGen, has been developed that realizes this methodology and automatically yields fixed-Point code corresponding to a given filter algorithm with a guarantee on the output error.

Dans cette thèse nous essayons d'améliorer l'évaluation de filtres numériques en nous concentrant sur la précision de calcul nécessaire.Ce travail est réalisé dans le contexte d'un générateur de code matériel/logiciel fiable pour des filtres numériques linéaires, en particulier filtres à Réponse Impulsionnelle Infinie (IIR). Avec ce travail, nous mettons en avant les problèmes liés à l'implémentation de filtres linéaires en arithmétique Virgule Fixe tout en prenant en compte la précision finie des calculs nécessaires à la transformation de filtres vers code. Ce point est important dans le cadre de filtres utilisés dans des systèmes embarqués critique comme les véhicules autonomes. Nous fournissons une nouvelle méthodologie pour l'analyse d'erreur lors de l'étude d'algorithmes de filtres linéaires du point de vue de l'arithmétique des ordinateurs. Au cœur de cette méthodologie se trouve le calcul fiable de la mesure Worst Case Peak Gain d'un filtre qui est la norme l1 de sa réponse impulsionnelle. L'analyse d'erreur proposée est basée sur la combinaison de techniques telles que l'analyse d'erreur en Virgule Flottante, l'arithmétique d'intervalles et les implémentations multi-précisions. Cette thèse expose également la problématique de compromis entre les coûts matériel (e.g. la surface) et la précision de calcul lors de l'implémentation de filtres numériques sur FPGA. Nous fournissons des briques de bases algorithmiques pour une solution automatique de ce problème. Finalement, nous intégrons nos approches dans un générateur de code pour les filtres au code open-source afin de permettre l'implémentation automatique et fiable de tout algorithme de filtre linéaire numérique. / In this thesis we develop approaches for improvement of the numerical behavior of digital filters with focus on the impact of accuracy of the computations. This work is done in the context of a reliable hardware/software code generator for Linear Time-Invariant (LTI) digital filters, in particular with Infinite Impulse Response (IIR). With this work we consider problems related to the implementation of LTI filters in Fixed-Point arithmetic while taking into account finite precision of the computations necessary for the transformation from filter to code. This point is important in the context of filters used in embedded critical systems such as autonomous vehicles. We provide a new methodology for the error analysis when linear filter algorithms are investigated from a computer arithmetic aspect. In the heart of this methodology lies the reliable evaluation of the Worst-Case Peak Gain measure of a filter, which is the l1 norm of its impulse response. The proposed error analysis is based on a combination of techniques such as rigorous Floating-Point error analysis, interval arithmetic and multiple precision implementations. This thesis also investigates the problematic of compromise between hardware cost (e.g. area) and the precision of computations during the implementation on FPGA. We provide basic brick algorithms for an automatic solution of this problem. Finally, we integrate our approaches into an open-source unifying framework to enable automatic and reliable implementation of any LTI digital filter algorithm.

Le sujet principal de cette thèse porte sur l'approximation polynômiale des fonctions aléatoires au moyen de la projection L2 aléatoire discrète, et son application aux problèmes inverses pour les équations aux dérivées partielles avec des données aléatoires. Les motivations proviennent de l'approximation paramétrique de la solution de modèles aux dérivées partielles. La thèse se compose de deux parties, avec un chapitre d'introduction qui résume les techniques modernes de l'approximation polynômiale des fonctions de variables aléatoires. La première partie, du chapitre 1 au chapitre 4, contient l'analyse théorique de la projection L2 aléatoire discrète pour résoudre le problème direct, par exemple, pour rapprocher les moments d'une fonction aléatoire à partir de ses observations, ou pour calculer la solution à un modèle numérique avec des coefficients stochastiques. La stabilité et l'optimalité de l'erreur d'approximation évaluée dans la norme L2 pondérée sont traités. Dans la dernière partie de la thèse, composé des chapitres 5 et 6, la méthodologie développée précédemment pour le problème direct est appliqué aux problèmes inverses pour les équations aux dérivées partielles à coefficients stochastiques. La méthode de factorisation est appliquée dans le cadre de la tomographie par impédance électrique, d'abord dans le cas de coefficient inhomogène, puis dans le cas de coefficient constante par morceaux, à valeurs dans chaque région affectée par l'incertitude. Enfin, dans le chapitre 6 les variantes de la méthode de factorisation proposées dans le chapitre précédent sont accélérés en utilisant les techniques qui ont été présentés dans la première partie de la thèse.

La modélisation et la simulation numérique des interfaces sont au coeur de nombreuses applications en mécanique des fluides et des solides, telles que la biologie cellulaire (déformation des globules rouges dans le sang), l'ingénierie pétrolière et la sismique (modélisation de réservoirs, présence de failles, propagation des ondes), l'aérospatiale (problème de rupture, de chocs) ou encore le génie civil. Cette thèse porte sur l'approximation des problèmes d'interface et d'interphase en mécanique des fluides par la méthode NXFEM, qui permet de prendre en compte de façon précise une discontinuité non alignée avec le maillage. Nous nous sommes d'abord intéressés au développement de la méthode NXFEM pour des éléments finis non-conformes pour prendre en compte une interface séparant deux milieux. Nous avons proposé deux approches pour les équations de Darcy et de Stokes. La première consiste à modifier les fonctions de base de Crouzeix-Raviart sur les cellules coupées et la deuxième consiste à rajouter des termes de stabilisation sur les arêtes coupées. Les résultats théoriques obtenus ont été ensuite validés numériquement. Par la suite, nous avons étudié la modélisation asymptotique et l'approximation numérique des problèmes d'interphase, faisant apparaître une couche mince. Nous avons considéré d'abord les équations de Darcy en présence d'une faille et, en passant à la limite dans la formulation faible, nous avons obtenu un modèle asymptotique où la faille est décrite par une interface, avec des conditions de transmission adéquates. Pour ce problème limite, nous avons développé une méthode numérique basée sur NXFEM avec éléments finis conformes, consistante et stable. Des tests numériques, incluant une comparaison avec la littérature, ont été réalisés. La modélisation asymptotique a été étendue aux équations de Stokes, pour lesquelles nous avons justifié le modèle limite obtenu. Enfin, nous nous sommes intéressés à la modélisation de la membrane d'un globule rouge par un fluide non-newtonien viscoélastique de Giesekus, afin d'appréhender la rhéologie du sang. Pour un problème d'interphase composé de deux fluides newtoniens (l'extérieur et l'intérieur du globule) et d'un liquide de Giesekus (la membrane du globule), nous avons dérivé formellement le problème limite, dans lequel les équations dans la membrane sont remplacées par des conditions de transmission sur une interface. / Numerical modelling and simulation of interfaces in fluid and solid mechanics are at the heart of many applications, such as cell biology (deformation of red blood cells), petroleum engineering and seismic (reservoir modelling, presence of faults, wave propagation), aerospace and civil engineering etc. This thesis focuses on the approximation of interface and interphase problems in fluid mechanics by means of the NXFEM method, which takes into account discontinuities on non-aligned meshes.We have first focused on the development of NXFEM for nonconforming finite elements in order to take into account the interface between two media. Two approaches have been proposed, for Darcy and Stokes equations. The first approach consists in modifying the basis functions of Crouzeix-Raviart on the cut cells and the second approach consists in adding some stabilization terms on each part of a cut edge. We have studied them from a theoretical and a numerical point of view. Then we have studied the asymptotic modelling and numerical approximation of interphase problems, involving a thin layer between two media. We have first considered the Darcy equations in the presence of a highly permeable fracture. By passing to the limit in the weak formulation, we have obtained an asymptotic model where the 2D fracture is described by an interface with adequate transmission conditions. A numerical method based on NXFEM with conforming finite elements has been developed for this limit problem, and its consistency and uniform stability have been proved. Numerical tests including a comparison with the literature have been presented. The asymptotic modelling has been finally extended to Stokes equations, for which we have justified the limit problem. Finally, we have considered the mechanical behaviour of red blood cells in order to better understand blood rheology. The last part of the thesis is devoted to the modelling of the membrane of a red blood cell by a non-Newtonian viscoelastic liquid, described by the Giesekus model. For an interphase problem composed of two Newtonian fluids (the exterior and the interior of the red blood cell) and a Giesekus liquid (the membrane), we formally derived the limit problem where the equations in the membrane are replaced by transmission conditions on an interface.