Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Nous présentons une nouvelle méthode de type frontière immergée (immersed boundary method, ou méthode IB) pour le calcul d'écoulements visqueux incompressibles en géométries irrégulières. Dans les méthodes IB , la frontière irrégulière de la géométrie n'est pas alignée avec la grille de calcul, et le point crucial de leur développement demeure le traitement numérique des cellules fluides qui sont coupées par la frontière irrégulière, appelées cut-cells. La partie dédiée à la résolution des équations de Navier-Stokes de notre méthode IB, appelée méthode LS-STAG , repose sur la méthode MAC pour grilles cartésiennes décalées, et sur l'utilisation d'une fonction de distance signée (la fonction level-set ) pour représenter précisément les frontières irrégulières du domaine. L'examen discret des lois globales de conservation de l'écoulement (masse, quantité de mouvement et énergie cinétique) a permis de bâtir une discrétisation unifiée des équations de Navier-Stokes dans les cellules cartésiennes et les cut-cells . Cette discrétisation a notamment la propriété de préserver la structure à 5 points du stencil original et conduit à une méthode extrêmement efficace sur le plan du temps de calcul en comparaison à un solveur non-structuré. La précision de la méthode est évaluée pour l'écoulement de Taylor-Couette et sa robustesse éprouvée par l'étude de divers écoulements instationnaires, notamment autour d'objets profilés. Le champ d'application de notre solveur Newtonien s'étend au cas d'écoulements en présence de géométries mobiles, et la méthode LS-STAG s'avère être un outil prometteur puisqu'affranchie des étapes systématiques (et coûteuses) de remaillage du domaine. Finalement, la première application d'une méthode IB au calcul d'écoulements de fluides viscoélastiques est présentée. La discrétisation de la loi constitutive est basée sur la méthode LS-STAG et sur l'utilisation d'un arrangement totalement décalé des variables dans tout le domaine assurant le couplage fort requis entre les variables hydrodynamiques et les composantes du tenseur des contraintes élastiques. La méthode est appliquée au fluide d'Oldroyd-B en écoulement dans une contraction plane 4:1 à coins arrondis. / This thesis concerns the development of a new Cartesian grid / immersed boundary (IB) method for the computation of incompressible viscous flows in two-dimensional irregular geometries. In IB methods, the computational grid is not aligned with the irregular boundary, and of upmost importance for accuracy and stability is the discretization in cells which are cut by the boundary, the so-called ``cut-cells''. In this thesis, we present a new IB method, called the LS-STAG method, which is based on the MAC method for staggered Cartesian grids and where the irregular boundary is sharply represented by its level-set function. This implicit representation of the immersed boundary enables us to calculate efficiently the geometry parameters of the cut-cells. We have achieved a novel discretization of the fluxes in the cut-cells by enforcing the strict conservation of total mass, momentum and kinetic energy at the discrete level. Our discretization in the cut-cells is consistent with the MAC discretization used in Cartesian fluid cells, and has the ability to preserve the 5-point Cartesian structure of the stencil, resulting in a highly computationally efficient method. The accuracy and robustness of our method is assessed on canonical flows at low to moderate Reynolds number~: Taylor Couette flow, flows past a circular cylinder, including the case where the cylinder has forced oscillatory rotations. We extend the \em LS-STAG \em method to the handling of moving immersed boundaries and present some results for the transversely oscillating cylinder flow in a free-stream. Finally, we present the first IB method that handles flows of viscoelastic fluids. The discretization of the constitutive law equation is based on the \em LS-STAG \em method and on the use of a fully staggered arrangement of unknowns, which ensures a strong coupling between all flow variables in the whole domain. The resulting method is applied to the flow of an Oldroyd-B fluid in a 4:1 planar contraction with rounded corner.

L'analyse statique a pour but de vérifier qu'un programme a le comportement souhaité c.à.d. satisfait des propriétés de sûreté. Toutefois, inférer les propriétés vérifiées par un programme est un problème difficile : le théorème de Rice énonce que toute propriété non triviale d'un langage de programmation Turing-complet est indécidable. Afin de contourner cette difficulté, les analyses statiques effectuent des approximations des comportements possibles du programme. La théorie de l'interprétation abstraite permet de donner un cadre formel à ces approximations. Cette théorie, introduite par Cousot & Cousot propose un cadre d'approximation basé sur la notion de treillis, de connexion de Galois et de calculs de points fixes par itération. Ce cadre permet de définir la qualité des approximations effectuées et notamment la notion de meilleure approximation. À l'opposé, les notions quantitatives n'apparaissent pas naturellement dans ce cadre. Nous nous sommes donc posés la question de l'inférence, par analyse statique, de propriétés s'exprimant de manière quantitative (telles que l'utilisation de la mémoire ou le temps d'exécution). / Static analysis aims to verify that programs behave correctly i.e. satisfy safety properties. However, generating properties verified by a program is a difficult problem : Rice’s theorem states that any non-trivial property about the language recognized by a Turing machine is undecidable. In order to avoid this difficulty, static analyses approximate the possible behaviours of the program. Abtract interpretation theory defines a formal framework for approximating programs. This theory, introduced by Cousot & Cousot is based on the mathematical structure of lattices, Galois connections and iterative fixpoints calculus. This framework defines the notion of correct approximation and allows for qualitatively compare approximations. On the contrary, it is not suitable for handling quantitative properties (such as memory usage and execution time).

Le travail présenté dans cette thèse vise à contribuer à l'analyse des phénomènes magnéto-mécanique et magnéto-thermique à l'aide de la méthode des éléments finis. Généralement, les problèmes magnéto-mécanique et magnéto-thermique sont constitués de plusieurs sous-problèmes dont la nature physique est différente. En utilisant la stratégie du couplage faible , ces sous-problèmes peut être calculés séparément en utilisant la méthode des éléments finis et donc sur leurs propres maillages afin d'assurer la précision. Pour obtenir une solution précise pour l'ensemble du problème multi-physique, il est très important d'assurer la transmission d'informations entre les sous-problèmes. Dans ce travail , nous étudions les méthodes de la projection. Les formulations de la projection sont données pour l'espaces L2 (scalaire et vectoriel), mais aussi pour H(grad), H(rot) et H(div ) afin d'améliorer la précision sur les dérivés. Une méthode de Petrov-Galerkin est présenté pour remplir l'espace test avec une base bi-orthogonale, qui permet de réduire le coût de calcul des projections L2 Ritz-Galerkin. Les techniques d'implémentation sont également exprimées en détails au niveau de la précision, la rapidité et la simplicité de réalisation. Ensuite avec les formulations énergétiques des champs électromagnétiques, des applications de la projection sont démontrées. La précision et l'efficacité de la projection sont données dans les problèmes multi -physiques. La thèse se termine avec les conclusions et des possibilités pour l'avenir. / The work presented in this thesis aims to contribute to finite element analysis of magneto-mechanical and magneto-thermal phenomena. In general, magneto-mechanical and magneto-thermal problems are made up of subproblems of which the physical nature differs. Using weak-coupling strategies, these subproblems can be calculated separately using finite element methods and thus on their own meshes in order to ensure precision. To obtain a precise solution for the entire problem, it is crucial to ensure the transmission of information between the subproblems. In this work, we study field projection methods on overlapping domains. Field projection formulations are given for classical L2 space, as well as for H(grad), H(curl) and H(div) in order to obtain increased projection accuracy for the distributional derivatives. A Petrov-Galerkin method is presented to fill the test space using a bi-orthogonal basis, in order to reduce the computation cost of L2 or L2 Ritz-Galerkin projections. Practical implementation techniques are also discussed in details for the consideration of accuracy, speed and simplicity of realization. Subsequently, with energy-conserving formulations for electromagnetic fields, applications of mesh-to-mesh projections are demonstrated. The accuracy and efficiency of the presented projection methods are given through multi-physics problems. The thesis closes with some conclusions and possibilities for future work.

La propension des systèmes informatiques à subir des défaillances matérielles est à l'origine d'une recherche abondante afin de concevoir des systèmes dits tolérants aux pannes. Le procédé couramment retenu consiste à procéder à des réplications, donnant alors naissance à ce que l'on nomme un système distribué. La question se pose alors de savoir si l'on peut garantir que les multiples copies sont cohérentes entre elles. Ainsi, la recherche d'un accord devient-elle un problème à résoudre, à portée paradigmatique : le Consensus. Or, la complexité des algorithmes de Consensus rend la tache ardue : il n'est donc pas rare que l'on commette des erreurs lors de leur conception. De là découle l'idée, développée depuis plus de trente ans, de recourir à des procédés de vérification mécanique des algorithmes et de leurs preuves de correction. Ces procédés prennent place parmi ce que l'on désigne usuellement comme étant des méthodes formelles. C'est à la croisée des recherches en algorithmique distribuée et en méthodes formelles que se situent nos travaux. Plus spécifiquement, il s'agit de faire usage d'un logiciel de certification formelle, Isabelle/HOL, afin de garantir l'exactitude des preuves de correction d'algorithmes de Consensus exprimés dans un cadre formel uniforme du nom de Heard-Of, proposé en 2009 par Charron-Bost et Schiper. Nous montrons que, du fait de leur expression dans un même cadre formel, et du fait de leur proximité, suivant les cas, soit de conception (nombre de rondes, recours à des mécanismes de vote, ...) soit de forme syntaxique, soit d'hypothèses de fonctionnement (synchronisme partiel, ...), ces algorithmes présentent des preuves dont une part conséquente d?arguments sont communs. Cela permet de copier certains d'entre eux d'une preuve à l'autre, afin de réduire l'effort de certification : ces arguments peuvent alors être automatiquement évalués par la machine pour chacun d'entre eux, l'utilisateur n'ayant à intervenir que là où celle-ci est en peine, c'est-à-dire lorsque les différences algorithmiques induisent qu'il faille réviser les détails de l'argumentation. L'exposé que nous faisons de la certification que nous avons effectuée pour six algorithmes distribués dédiés à la résolution du problème du Consensus illustre cette démarche. Par conséquent, nous présentons d'abord les portions communes des démonstrations, puis détaillons ce qui est propre à chacune, l'objectif n'étant pas de permettre une lecture linéaire de chaque démonstration mais de mettre en évidence notre proposition / Computer systems fail. Whatever the reason of these failures, it has been a widespread approach to try and increase the faults-tolerance of a computer system through its replication. The resulting system is said to be a distributed one, in which replicas have to be kept consistent with each others. Hence, reaching agreement, and Consensus in particular, becomes the problem to solve - indeed, the paradigm. Solving Consensus (under various assumptions) is a hard task : algorithms designed on this purpose are subtle and proving their being correct is error-prone - whenever they are, which occasionally appears not to be the case. For more that thirty years, researchers interested in what is called Formal Methods have been working on mechanizing the verification process, in order to increase confidence in the correctness of (distributed) algorithms. The work we present here is at the intersection of distributed algorithms and formal methods. We use the Isabelle/HOL software to certify the correctness proof of various Consensus algorithms expressed in a uniform framework based on the Heard-Of Model, introduced by Charron-Bost and Schiper in 2009. Expressed in a common model, these algorithms, which, depending on the case, share some common mecanisms (number of steps, intermediate votes, ...), some elements of syntax, or types of assumptions (partial synchronism...), can be proved using some common arguments. As a consequence, the certification effort can be reduced by copying some intermediate lemmas from one proof to another and let the computer automatically parse them until some manual adaptation is required. This lead to the idea of certifying the correctness of multiple algorithms all together, instead of proving them one after the other, as one would do on paper in a traditional way. The effort of translation in the formal language of the proof assistant is then possibly reduced. Of course, each proof will also contain specific arguments, which will have to be isolated and translated into the software. Here, we illustrate this proposition through the presentation of formal certificates of correctness for six Consensus algorithms. As a consequence, on should not expect to find here a comprehensive linear presentation of each proof : we first show the arguments shared by multiple proofs, followed by those which are specific to each o them

Les utilisateurs s'attendent à ce que les logiciels soient à la fois continuellement disponibles et mis-à-jour à la volée. L'introduction de nouvelles fonctionnalités et l'intégration de correctifs ne doit pas perturber la bonne marche d'une application. Les développeurs bénéficient aussi de la modification dynamique des logiciels en adaptant le code à de nouvelle situations ou bien an exécutant des analyses dynamiques à des fins de monitorage, de déboguage et d’optimisation. Les langages de programmation doivent donc fournir aux développeurs des mécanismes qui améliorent l'adaptabilité logicielle. Dans le même temps, l'augmentation de la taille et de la complexité des logiciels requièrent des applications faites de nombreux composants développés par différents groupes. Dans ce contexte chaque composant doit être convenablement protégé des autres. Les langages de programmation doivent donc aussi fournir des mécanismes qui améliorent l'encapsulation. Le besoin de supporter des adaptations dynamiques d'un côté, et le besoin d'avoir des composants convenablement confinés d'un autre, poussent la conception des langages de programmation dans deux directions opposées: soit plus de flexibilité, soit plus de contrôle. Cette dissertation étudie cette tension dans le contexte des langages orientée-objet dynamiquement typés via deux mécanismes qui promeuvent l'adaptabilité logicielle: la réflexion et les méthodes d’extension. Pour chacun de ces mécanismes, nous démontrons leur utilité, leurs effets néfastes sur l'encapsulation et proposons des solutions qui conservent leurs forces tout en maintenant l'encapsulation. / Users expect modern software to be both continually available and updated on the fly. Introduction of new features and integration of bug fixes should not trouble the smooth running of an application. Developers also benefit from dynamic modifications of software by adapting code to new unanticipated situations or by running dynamic analyses to get feedback about the behavior of running applications for monitoring, debugging and optimization purposes. Programming languages should thus provides developers with mechanisms that improve software adaptability. At the same time, the increasing size and complexity of software call for applications made up of many interacting components developed by different parties. In this context, each component needs to be properly protected from other ones. Programming languages should thus provide developers with mechanisms that improve encapsulation. The need for dynamic adaptations on the one hand and the need for properly confined components on the other hand pulls programming language design in two opposite directions: either more flexi-bility or more control.This dissertation studies this tension in the context of dynamically-typed object-oriented languages with two language mechanisms that promote software adaptability: reflection and extension meth-ods. For both mechanisms, we demonstrate their usefulness, their harmful effect on encapsulation and propose solutions to retain their power while maintaining encapsulation.

Safety-critical systems depend on the fact that their software components provide services that behave correctly (i.e. satisfy their requirements). Additionally, in many cases, these systems have to be adapted or reconfigured in case of failures or when changes in requirements or in quality of service occur. When these changes appear at the software level, they can be handled by the notion of substitution. Indeed, the software component of the source system can be substituted by another software component to build a new target system. In the case of safety-critical systems, it is mandatory that this operation enforces that the new target system behaves correctly by preserving the safety properties of the source system during and after the substitution operation. In this thesis, the studied systems are modeled as state-transition systems. In order to model system substitution, the Event-B method has been selected as it is well suited to model such state-transition systems and it provides the benefits of refinement, proof and the availability of a strong tooling with the Rodin Platform. This thesis provides a generic model for system substitution that entails different situations like cold start and warm start as well as the possibility of system degradation, upgrade or equivalence substitutions. This proposal is first used to formalize substitution in the case of discrete systems applied to web services compensation and allowed modeling correct compensation. Then, it is also used for systems characterized by continuous behaviors like hybrid systems. To model continuous behaviors with Event-B, the Theory plug-in for Rodin is investigated and proved successful for modeling hybrid systems. Afterwards, a correct substitution mechanism for systems with continuous behaviors is proposed. A safety envelope for the output of the system is taken as the safety requirement. Finally, the proposed approach is generalized, enabling the derivation of the previously defined models for web services compensation through refinement, and the reuse of proofs across system models.

L'objet de cette thèse est le développement de méthodes variationnelles pour l'étude des solutions presque-périodiques au sens de H. Bohr et au sens de Besicovitch de quelques classes d'équations différentielles ordinaires du second ordre à retard. Pour cela on utilise le Calcul des Variations en Moyenne Temporelle. Dans un premier temps on étudie une classe d'équations différentielles du type neutre, puis une classe d'équations différentielles à retard fini, enfin on s'intéresse à une classe d'équations différentielles à retard infini.

Cette HDR est dédiée à l'étude du management de l'innovation. Elle est orientée autour de trois axes : la complémentarité entre innovations technologique et organisationnelle, l'organisation de la propriété industrielle et une réflexion sur les méthodologies qualitatives en sciences de gestion.