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

Cette thèse s'intéresse à la conception de méthodes de discrétisation non-conforme pour un modèle de poromécanique. Le but de ce travail est de simplifier les couplages liant la géomécanique d'un milieu poreux à l'écoulement polyphasique compositionnel ayant cours en son sein tels qu'ils sont réalisés actuellement dans l'industrie pétrolière, en discrétisant sur un même maillage, typiquement non-conforme car à l'image de la lithologie, la mécanique et l'écoulement. La nouveauté consiste donc à traiter la mécanique par une méthode d'approximation non-conforme sur maillages généraux. Dans cette thèse, nous nous concentrons sur un modèle d'élasticité linéaire. Les difficultés inhérentes à son approximation non-conforme sont son manque de coercivité (se traduisant par la nécessité de satisfaire une inégalité de Korn sur un espace discret discontinu), ainsi que le phénomène de verrouillage numérique lorsque le matériau tend à devenir incompressible. Dans une première partie, nous construisons un espace d'approximation sur maillages généraux, s'apparentant à une extension de l'espace de Crouzeix-Raviart. Nous explicitons ses propriétés d'approximation et de conformité, et montrons que ce dernier est adapté à une discrétisation primale coercive et robuste au locking du modèle d'élasticité sur maillages généraux. La méthode proposée est moins coûteuse que son équivalent éléments finis (en termes de propriétés) P2. Nous nous intéressons dans une deuxième partie à l'approximation non-conforme d'un modèle couplé de poroélasticité. Nous étudions la convergence d'une famille de schémas numériques dont la discrétisation en espace utilise le formalisme des schémas Gradient, auquel appartient la méthode développée pour la mécanique. Nous prouvons la convergence de telles approximations vers la solution de régularité minimale du problème continu, indépendamment des paramètres physiques du système

Dans cette thèse, nous proposons un modèle pour le calcul des écoulements le long des failles. Ce modèle, baptisé modèle double interface permet de traiter deux difficultés majeures rencontrées lors de la modélisation des failles. Tout d'abord, l'utilisation d'un modèle interface, dans lequel les failles sont représentées par des éléments de dimension inférieure permet de s'affranchir du problème d'échelle spatiale. Ensuite, l'utilisation de deux interfaces pour représenter chaque faille permet de traiter naturellement les maillages non-conformes apparaissant dans ce type de problème. Les questions de failles non-planes et de réseaux de failles sont aussi abordées. Ce modèle est validé numériquement sur différents cas tests académiques et un cas synthétique inspiré du stockage du CO2 a aussi été réalisé. Finalement, une étude théorique a été menée afin de confirmer mathématiquement l'approche retenue. / In this thesis, we are interested in the modelisation of fluid flow along conductive faults. This model, so-called double interface model tackles two majors difficulties encountered when modelising faults. First of all, the use of an interface model, in which the faults are represented by lower dimension elements allows to treat the problem of space scale. Then, the use of two interfaces to modelise each fault allows to handle quite naturally the non-matching grid problem arising from this kind of problem. The question of non-planar fault and fault networks is also addressed. This model is then validated on several academic test cases and a synthetic case inspire by CO2 storage is also performed. Finally, a theoric study is also conducted in order to validate our approach.

Le problème de la mécanique des fluides, concernant les écoulements décollés derrière des obstacles immobiles ou en mouvement, est traité par l’étude de deux sujets: a) recherches sur l’existence de solutions stationnaires des equations d’Euler et Navier-Stokes pour grands nombres de Reynolds, au-delà des corps caractérisés par pointes ou singularités géométriques; b) analyse du sillage non stationnaire derrière une turbine à axe vertical (VAT). L’étude de ces deux différents régimes d’écoulements, concernant le phénomène du détachement derrière corps émoussés ou profils alaires à haut angle d’incidence, a permis la mise au point de plusieurs techniques analytiques et numériques basées sur le champ de vorticité. / The problem of the separated flows dynamics past obstacles at rest or moving bodies is addressed by means of the study of two topics a) investigation on the existence of some steady solutions of the Euler equations and of the Navier-Stokes equations at large Reynolds number, past bodies characterized by a cusp; b) analysis of the unsteady wake behind a Vertical Axis Turbine (VAT). The survey of such different flow regimes related to the separation phenomenon past bluff bodies or bodies at incidence allowed to devise several numerical and analytical techniques based on the evaluation of the vorticity field.

Neste trabalho, discutem-se aspectos das transformações conformes na eletrodinâmica clássica com ênfase na invariância e nas leis de conservação. Inicialmente, abordaram-se aspectos gerais das transformações conformes e fez-se um resumo histórico da evolução dessas transformações. Procurou-se fazer uma apresentação didática, revisando-se a formulação Lagrangiana e o Teorema de Noether para campos aplicado à eletrodinâmica. Estudaram-se as transformações conformes no espaço plano, onde se mostrou que para dimensões maiores ou iguais a três o número de transformações é finito. A partir das equações de Maxwell em coordenadas curvilíneas, chegou-se à condição para que essas equações mantivessem sua forma cartesiana. Com essa condição, mostrou-se que a eletrodinâmica clássica é invariante para o grupo de transformações conformes. Foram discutidas as leis de conservação associadas à invariância conforme da eletrodinâmica clássica a partir do teorema de Noether. Das simetrias por translações no espaço-tempo, obtiveram-se as leis de conservação do momento linear e da energia. Das simetrias associadas às rotações, obtiveram-se seis quantidades conservadas: três delas ligadas à conservação do momento angular e, com relação às três restantes, observou-se, a partir de analogias com a mecânica, que estavam associadas ao movimento do centro de energia do campo. Para a interpretação da grandeza conservada por simetria de escala, verificou-se, também a partir de uma analogia mecânica, que essa simetria somente é verificada para partículas não massivas ou para partículas massivas a altas energias. Finalmente, para as transformações conformes especiais, verificou-se que as leis de conservação resultantes são consequências das leis anteriores de conservação para o campo eletromagnético, e neste caso, essa simetria também somente se manifesta para partículas de massa nula ou para altas energias. / In this work, aspects of conformal transformations in classical electrodynamics are discussed with emphasis on the invariance and conservation laws. Initially, a general view of conformal transformations was shown and a summary of the historical evolution of those transformations was presented. The work was approached didactically, and Noethers theorem based on the electrodynamics Lagrangian formulation was revised. The conformal transformations were studied in plane spaces and it was shown that, for dimensions greater than or equal to three, the number of transformations is finite. Starting from Maxwells equations in curvilinear coordinates, a condition for maintaining those equations in Cartesian form was established. With that condition, it was shown that the classical electrodynamics laws are invariant for the group of conformal transformations. The conservation laws associated with the conformal invariance of classical electrodynamics were discussed, based on Noethers theorem. From the space-time translation symmetry, the laws of conservation of linear momentum and of energy were obtained. From rotational symmetry, six conserved quantities were obtained: three of them associated with angular momentum and the remaining three, observed, starting from analogies with mechanics, were associated with the movement of the center of energy of the field. For the interpretation of the quantity conserved by scale symmetry, it was verified, also from a mechanical analogy, that that symmetry is only valid for null mass particles or for high energies. Finally, for the special conformal transformations, it was verified that the resultant laws of conservation are consequences of the previous laws, and in that case, symmetry is also valid only for particles of null mass or for high energies.

The objective of this work is the study of phase transitions on the context of Groupoids and their C*-Algebras. The main result of this dissertation is due to Klaus Thomsen in [Tho17], which investigates the connection between conformal measures in the classical formalism and KMS-states in the quantum formalism. The phase transition in the quantum setting is a consequence of this connection between both formalisms and the fact that on the classical setting it was known examples of continuous potentials that show the phenomena of phase transition. The potential used was introduced by Hofbauer [Hof77], an example that shows, dierently from potential of summable variations, potentials only continuous can exhibit phase transition. / O objetivo deste trabalho é o estudo do fenômeno de transição de fase no contexto de Grupóides e suas C*-álgebras. O resultado principal é devido a Klaus Thomsen em [Tho17], que explora a conexão entre medidas conformes no formalismo clássico e estados KMS do contexto quântico. A transição de fase no caso quântico é consequência desta ligação entre os dois formalismos e do fato de que no setting clássico eram conhecidos exemplos de potenciais contínuos que apresentam o fenômeno de transição de fase. O potencial utilizado é aquele introduzido por Hofbauer [Hof77], um exemplo que mostra que, diferentemente de potenciais de variação somável, potenciais apenas contínuos podem apresentar transição de fase.

Cette thèse propose l’étude d’une méthode numérique permettant de simuler un plasma. On considère un ensemble de particules, dont le mouvement est régi par l’équation de Vlasov, et qui est sensible aux forces électromagnétiques, qui proviennent des équations de Maxwell. La résolution numérique des équations de Vlasov-Maxwell est réalisée par une méthode Particle In Cell (PIC). La résolution des équations de Maxwell nécessite un maillage suffisamment fin afin de modéliser correctement les problémes multi-échelles que nous souhaitons traiter. Cependant, mailler finement tout le domaine de calcul a un coût. La nouveauté de cette thèse est de proposer un solveur PIC sur des maillages cartésiens localement raffinés, des maillages non conformes, afin de garantir la bonne modélisation du phénomène physique et d’éviter une trop forte pénalisation des temps de calcul.Nous utilisons une méthode Galerkin Discontinue en domaine temporelle (GDDT), qui offre l’avantage d’être d'une grande flexibilité dans le choix du maillage et qui est une méthode d’ordre élevé. Un point fondamental dans l’étude des solveurs PIC concerne le respect de la conservation de la charge. Nous proposons deux approches afin de traiter cet aspect. La première concerne les méthodes utilisant un système de Maxwell augmenté, dont la nouveauté a été de les étendre aux maillages non conformes. La seconde approche repose sur une méthode originale de pré-traitement du calcul du terme source de courant. / This thesis deals with the study of a numerical method to simulate a plasma. We consider a set of particles whose displacement is governed by the Vlasov equation and which creates an electromagnetic field thanks to Maxwell equations. The numerical resolution of the Vlasov-Maxwell system is performed by a Particle In Cell (PIC) method. The resolution of Maxwell equations needs a sufficiently fine mesh to correctly simulate the multi scaled problems that we have to face. Yet, a uniform fine mesh of the whole domain has a prohibitive cost. The novelty of this thesis is a PIC solver on locally refined Cartesian meshes : non conforming meshes, to guarantee the good modeling of the physical phenomena and to avoid too large CPU time. We use the Discontinuous Galerkin in Time Domain (DGTD) method which has the advantage of a great flexibility in the choice of the mesh and which is a high order method. A fundamental point in the study of PIC solvers is the respect of the charge conserving law. We propose two approaches to tackle this point. The first one deals with augmented Maxwell systems, that we have adapted to non conforming meshes. The second one deals with an original method of preprocessing of the calculation of the current source term.

Submitted by Cássia Santos (cassia.bcufg@gmail.com) on 2016-03-09T12:52:41Z No. of bitstreams: 2 Dissertação - Ricardo Ferreira da Cunha - 2014.pdf: 3971057 bytes, checksum: b27175e5abd2d9d25ecf057c50549a54 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2016-03-10T11:54:28Z (GMT) No. of bitstreams: 2 Dissertação - Ricardo Ferreira da Cunha - 2014.pdf: 3971057 bytes, checksum: b27175e5abd2d9d25ecf057c50549a54 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2016-03-10T11:54:28Z (GMT). No. of bitstreams: 2 Dissertação - Ricardo Ferreira da Cunha - 2014.pdf: 3971057 bytes, checksum: b27175e5abd2d9d25ecf057c50549a54 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2014-07-11 / In this work we describe soem properties of the complex numbers and analytical functions. As an application, we use the conformal transformations to caculate analytical solutions for electromagnetic problems defined in nonconventional domain. The conformal transformation are of excepcional importance in solving boundary value problems in electromagnetic theory. The conformal transformation have the property of modifying the geometry, but preserving the pysical quantities. / Neste trabalho descrevemos algumas propriedades dos números complexos e das funções analíticas complexas. Como aplicação, usamos as transformações conformes para calcular de forma analítica a solução de problemas eletromagnéticos definidos em domínios não convencionais. As transformações conformes são de excepcional importância na resolução de problemas de valores de contorno na teoria eletromagnética. As transformações conformes têm a propriedade de modificar a geometria, preservando as grandezas físicas.

La première partie de cette thèse est consacrée à l'étude d'interfaces critiques bidimensionnelles par des méthodes d'évolutions de Schramm-Loewner (SLE) et de théories conformes. Nous étudions en particulier le cas de SLE(2) qui est la limite d'échelle des marches à boucles effacées. La solution explicite du problème d'enroulement sur des domaines doublement connexes est discutée. Nous établissons une généralisation de la formule de Schramm pour SLE(2) dans la géométrie doublement connexe et étendons la solution au cas de conditions mixtes Dirichlet-Neumann. L'analyse par la théorie conforme permet l'identification de l'opérateur de changement des conditions aux bords. De plus, à partir de l'étude des lignes de discontinuité du champ gaussien libre sur des domaines doublement connexes nous mettons en évidence une relation entre SLE(4) et les ponts browniens.<br /><br />Le sujet de la seconde partie est l'étude de deux exemples de systèmes désordonnés de basse dimension. D'un coté nous établissons les propriétés de localisation et spectrales d'un hamiltonien aléatoire unidimensionnel qui interpole entre les cas du modèle de Halperin et le modèle supersymétrique désordonné. Un lien avec la diffusion unidimensionnelle dans un potentiel aléatoire permet d'étudier la modification de la dynamique ultra-lente de Sinai en présence d'absorbeurs. De l'autre côté nous analysons la transition vitreuse d'ARN pour des séquences aléatoires à l'aide de la théorie des champs de Lässig-Wiese-David. L'application au cas d'ARN soumis à une force extérieure conduit à la prédiction de la caractéristique force-extension pour des séquences hétérogènes. L'étude de la phase vitreuse nous amène à considérer un modèle hiérarchique combinatoire dont nous déterminons les exposants et lois d'échelle exactes ainsi que les corrections de taille finie.

Cette thèse étudie la classification des théories conformes à 2d à l'aide de symétries quantiques de diagrammes. Les fonctions de partition d'un système conforme - l'invariante modulaire ou celles provenant de l'introduction de lignes de défauts - s'expriment en fonction d'un ensemble de coefficients qui forment des nimreps de certaines algèbres. Ces coefficients définissent les diverses structures d'une classe d'algèbres de Hopf, dites faibles, et peuvent être codés par un ensemble de graphes. Le chapitre 1 présente les connaissances actuelles sur ce sujet. Dans le chapitre 2 sont introduites l'algèbre de Hopf faible et ses structures, notamment l'algèbre des symétries quantiques d'Ocneanu, qui joue un rôle important dans l'étude des systèmes conformes à 2d. Nous analysons en détails ces structures pour le diagramme A3 du modèle affin su(2). Le chapitre 3 est dédié à la présentation d'une réalisation de l'algèbre des symétries quantiques d'Ocneanu, construite comme un quotient du carré tensoriel de l'algèbre d'un graphe G (de type ADE pour le modèle affin su(2)). Cette réalisation permet d'obtenir un algorithme simple permettant le calcul des fonctions de partition du modèle conforme associé. Notre construction se prête naturellement à une généralisation aux cas affins su(n), pour n > 2, pour lesquels peu de résultats étaient connus. Dans le chapitre 4, nous traitons explicitement tous les cas du type su(2) ainsi que trois exemples choisis du type su(3).