A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

Les composites thermoplastiques hautes performances ont suscité l'intérêt de l'industrie aéronautique pour remplacer l'aluminium dans certaines pièces. Ils permettent la production de pièces moins lourdes et moins coûteuses. Parmi ces matériaux, le PEEK renforcé par des fibres de carbone courtes semble intéressant pour des applications nécessitant de hautes propriétés thermomécaniques. Ce matériau a déjà fait l'objet de travaux au Laboratoire ICA de l'IMT Mines Albi, dans le cadre de la thèse de Jérémy Crevel ainsi que les projets régionaux Midi-Pyrénées CINTHTE, CINTHTE2 et FUI INMAT2. En outre, des applications industrielles sont en cours de développement chez Liebherr Aerospace. Cependant, le comportement à long terme en température de ce matériau est mal connu. L'objectif de cette thèse est de déterminer la plage d'utilisation du matériau en fonction du temps, de la température et du chargement. Deux phénomènes intervenant à hautes températures et à long terme sont étudiés : le vieillissement et le fluage. L'étude du vieillissement sous air s'est concentrée sur l'évolution des propriétés mécaniques en traction après une longue exposition à haute température. L'influence du vieillissement sur les propriétés mécaniques usuelles (module, contrainte à rupture et allongement à rupture) ainsi que sur l'endommagement et la plasticité a été mesurée en dessous et au-dessus de la température de transition vitreuse et pour différentes orientations de fibres. L'évolution de la structure du matériau avec le vieillissement a été suivie pour faire le lien avec les propriétés mécaniques. L'étude du fluage s'est faite à des contraintes inférieures au seuil d’endommagement du matériau. Des essais à court terme et dynamiques fréquentiels ont été réalisés pour définir un modèle de comportement dépendant du temps et l’identifier sur le long terme. Ils ont en outre permis de prendre en compte la température, le vieillissement et l'anisotropie du matériau pour des durées longues. Ces deux études (vieillissement et fluage) ont été unifiées dans un critère de rupture pour prévoir le temps à rupture du matériau. Enfin, un modèle viscoélastique isotrope transverse et des indicateurs de rupture en fluage ont été implémentés dans Abaqus pour réaliser des calculs sur des pièces industrielles. / The high performance thermoplastic based composites sparked interest of the aeronautical industry to replace aluminum in some parts. These composites permit to produce lighter and less expensive parts. Among these materials, short carbon fiber reinforced PEEK is interesting for applications requiring high thermomechanical properties. This material has already been studied in some preliminary works at the ICA-Albi laboratory, especially in the PhD thesis of Jérémy Crevel and projects funded by the Midi-Pyrénées region CINTHTE, CINTHTE2 and FUI INMAT2. Moreover, some industrial parts are currently developed by Liebherr Aerospace. However, the long term and high temperature behavior is ill known. The aim of this PhD thesis is to determine the range of use for the material in function of time, temperature and loading. Two phenomena occurring at high temperature for long period of time are studied: ageing and creep. Ageing was studied under air atmosphere at high temperature. Mechanical properties evolution under tensile traction (elastic modulus, ultimate stress and stress at break) and damage and plasticity are measured below and above glass transition temperature and for different fiber orientations. Similarly, influence of ageing on the material structure is studied and linked to mechanical results. Creep study is done at stresses where material does undergo no damage. Short term and dynamic mechanical tests are carried out to define and identify a time dependent mechanical behaviour law over long period of time. Moreover, they permit to take into account effects of temperature, ageing and anisotropy. Both previous studies (ageing and creep) are unified in a failure criteria to forecast failure of the material. Finally, a transverse isotropic viscoelastic model and failure indices have been implemented into Abaqus to realize calculation of an industrial parts.

A Left Ventricular Assist Device (LVAD) is a mechanical pump that helps patients with heart failure conditions. This rotary pump works in parallel to the ailing heart and provides an alternative path for blood flow from the weak left ventricle to the aorta. The LVAD is controlled by the power supplied to the pump motor. An increase in the pump motor power increases the pump speed and the pump flow. The LVAD is typically controlled at a fixed setting of pump power. This basically means that the controller does not react to any change in the activity level of the patient. An important engineering challenge is to develop an LVAD feedback controller that can automatically adjusts its pump motor power so that the resulting pump flow matches the physiological demand of the patient. To this end, the development of a mathematical model that can be used to accurately simulate the interaction between the cardiovascular system of the patient and the LVAD is essential for the controller design. The use of such a dynamic model helps engineers and physicians in testing their theories, assessing the effectiveness of prescribed treatments, and understanding in depth the characteristics of this coupled bio-mechanical system. The first contribution of this dissertation is the development of a pump power-based model for the cardiovascular-LVAD system. Previously, the mathematical models in the literature assume availability of the pump speed as an independent control variable. In reality, however, the device is controlled by pump motor power which, in turn, produces the rotational pump speed. The nonlinear relationship between the supplied power and the speed is derived, and interesting observations about the pump speed signal are documented. The second contribution is the development of a feedback controller for patients using an LVAD as either a destination therapy or a bridge to transplant device. The main objective of designing this controller is to provide a physiological demand of the patient equivalent of that of a healthy individual. Since the device is implanted for a long period of time, this objective is chosen to allow the patient to live a life as close to normal as possible. The third contribution is an analysis of the aortic valve dynamics under the support of an LVAD. The aortic valve may experiences a permanent closure when the LVAD pump power is increased too much. The permanent closure of the aortic valve can be very harmful to the patients using the device as a bridge to recovery treatments. The analysis illustrates the various changes in the hemodynamic variables of the patient as a result of aortic valve closing. The results establish the relationship between the activity level and the heart failure severity with respect to the duration of the aortic valve opening.

No description available.

Cette thèse porte sur les relations entre la limite supérieure de la zone sismogène et la construction et la structuration du prisme de Nankai, au Sud-Est du Japon. La détermination des contraintes dans le prisme à partir des données d’ovalisation de puits montrent une permutation des contraintes au dessus d’un hors séquence, appelé splay fault, situé à la limite de la zone sismogène. Le régime de contrainte est extensif dans cette zone et décrochant dans le reste du prisme. L’étude de la compaction des sédiments à partir des données de porosité corrigée de l’effet des argiles permet de quantifier l’érosion récente des sédiments de pentes, et montre l’existence d’une érosion plus ancienne au dessus des chevauchements majeurs. L’étude des styles de déformations des sédiments montre que l’érosion peut rendre compte de la répartition des structures de déformation dans le prisme. Enfin, la modélisation mécanique de la marge montre que l’activité de la splay fault, très faible depuis le dernier million d’années, et la présence d’une zone d’extension en arrière de cette faille s’expliquent par l’arrêt de l’accrétion et un décollement faible dans le prisme frontal. / This study is focused on the upper limit of the seismogenic zone in the Kumano transect of the Nankai margin in relation with the construction of the wedge. The first step was to determine the stress orientations and amplitude within the wedge from borehole breakouts analysis. We show the striking presence of a trench normal extension zone behind the splay fault in a wedge dominated by a strike-slip stress regime. Locally an extensional stress regime is also observed in the slope sediment of the outer part of the wedge. In a second part of this work, we used porosity corrected from water bound to clay minerals to study the compaction state of the sediments and to quantify the amount of erosion in the slope sediments. Results also highlight the existence of past erosion related to the activity of the splay fault. In a third part, the observation of deformation style within sediments showed that erosion could explain the distribution of dilatant and compactive structures within the wedge. Finally, we show, based on mechanical modeling, that the splay fault slow slip rate over the last million year, and the presence of an extension zone landward of the splay fault, could both be explained by a pause in the accretion and a weak décollement beneath the outer wedge.

A maioria dos líquidos encontrados na natureza são não-Newtonianos e o estudo do seu comportamento tem uma importância significante em diferentes áreas da engenharia. Entre eles, uma larga classe de materiais que exibem pequena ou nenhuma deformação quando sujeitos a um nível de tensões inferiores a uma tensão limite de escoamento – chamado de comportamento viscoplástico. A presente Dissertação tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos viscoplásticos não-lineares em uma cavidade forçada. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo viscoplástico recentemente introduzido por Souza Mendes e Dutra – SMD – e é aproximado por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão-extra, velocidade e pressão. As condições de compatibilidade entre os subespaços de elementos finitos para tensão-extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não-Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da UFRGS. Em diversos trabalhos encontrados na literatura, a superfície de escoamento do material é definida como a região onde o módulo da tensão-extra é igual à tensão limite de escoamento. É mostrado nesta Dissertação que esta metodologia pode conduzir à alguns erros, dado ao grande aumento experimentado pela taxa de cisalhamento em uma pequena faixa de tensões próximas à tensão limite de escoamento. Assim, foi adotada outra metodologia, definindo a superfície de escoamento como a linha onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade Newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o número de salto, J, o coeficiente de power-law, n, e a vazão adimensional, U*, são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos viscoplásticos. Os resultados obtidos estão de acordo com a literatura e atestam a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found on the nature and the study of their behavior has a significant importance on different areas of engineering. Among them, there is a wide class of materials that exhibits little or no deformation when subjected to a stress level behind an apparent yield stress – called the viscoplastic behavior. The present thesis aimed to a numerical study of two dimensional steady state laminar flows of non-linear viscoplastic fluids in a lid-driven cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the recently introduced Souza Mendes and Dutra – SMD – viscoplastic model and has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the extra-stress, velocity and pressure fields. In this way, the compatibility conditions between the extra-stressvelocity and pressure-velocity (Babuška-Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of UFRGS. In several works found on the literature, the yield surface of the material is defined as the region where the stress modulus is equal to the yield stress. Is shown in this work that this methodology can lead to some errors, due to the large strain rate increasing in a small range of values of stress on the vicinity of the yield stress. Therefore, it was adopted another approach, defining the yield surface as the line where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the jump number, J, the the power-law coefficient, n,and the non-dimensional flow rate, U*, are ranged in order to evaluate how they the influence on the viscoplastic fluid dynamics have been investigated. All results found were in accordance with the affine literature and attests the good stability features of the formulation.

A maioria dos líquidos encontrados na natureza são não-Newtonianos e o estudo do seu comportamento tem uma importância significante em diferentes áreas da engenharia. Entre eles, uma larga classe de materiais que exibem pequena ou nenhuma deformação quando sujeitos a um nível de tensões inferiores a uma tensão limite de escoamento – chamado de comportamento viscoplástico. A presente Dissertação tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos viscoplásticos não-lineares em uma cavidade forçada. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo viscoplástico recentemente introduzido por Souza Mendes e Dutra – SMD – e é aproximado por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão-extra, velocidade e pressão. As condições de compatibilidade entre os subespaços de elementos finitos para tensão-extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não-Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da UFRGS. Em diversos trabalhos encontrados na literatura, a superfície de escoamento do material é definida como a região onde o módulo da tensão-extra é igual à tensão limite de escoamento. É mostrado nesta Dissertação que esta metodologia pode conduzir à alguns erros, dado ao grande aumento experimentado pela taxa de cisalhamento em uma pequena faixa de tensões próximas à tensão limite de escoamento. Assim, foi adotada outra metodologia, definindo a superfície de escoamento como a linha onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade Newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o número de salto, J, o coeficiente de power-law, n, e a vazão adimensional, U*, são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos viscoplásticos. Os resultados obtidos estão de acordo com a literatura e atestam a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found on the nature and the study of their behavior has a significant importance on different areas of engineering. Among them, there is a wide class of materials that exhibits little or no deformation when subjected to a stress level behind an apparent yield stress – called the viscoplastic behavior. The present thesis aimed to a numerical study of two dimensional steady state laminar flows of non-linear viscoplastic fluids in a lid-driven cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the recently introduced Souza Mendes and Dutra – SMD – viscoplastic model and has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the extra-stress, velocity and pressure fields. In this way, the compatibility conditions between the extra-stressvelocity and pressure-velocity (Babuška-Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of UFRGS. In several works found on the literature, the yield surface of the material is defined as the region where the stress modulus is equal to the yield stress. Is shown in this work that this methodology can lead to some errors, due to the large strain rate increasing in a small range of values of stress on the vicinity of the yield stress. Therefore, it was adopted another approach, defining the yield surface as the line where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the jump number, J, the the power-law coefficient, n,and the non-dimensional flow rate, U*, are ranged in order to evaluate how they the influence on the viscoplastic fluid dynamics have been investigated. All results found were in accordance with the affine literature and attests the good stability features of the formulation.

Prolonger la durée de vie d’un parc nucléaire impose de garantir la pérennité des enceintes de confinement, troisième et dernière barrière de protection contre les potentielles émissions radioactives. Ces structures sont attentivement auscultées depuis la construction grâce à un dispositif de mesures spécifiques incluant des extensomètres noyés dans le béton. De conception intrinsèquement fiable et bien que largement éprouvés, ces capteurs de déformation commencent à présenter des taux de défaillance significatifs, nécessitant, de fait, le recours à des chaînes de mesure de substitution. En ne pouvant être mises en œuvre qu’au voisinage de la surface du béton, ces solutions ne permettent qu’un accès indirect à la déformation initialement observée. Des doutes peuvent dès lors être émis sur la représentativité des mesures et par conséquent sur la qualité de la surveillance.Pour analyser l'influence de la localisation de la mesure, le développement et l'exploitation de deux modèles numériques éléments finis complémentaires sont proposés. Il s’agira principalement, pour l’un, de quantifier les possibles écarts de comportement sur le long terme et en épreuve liés aux gradients de température et d’humidité ; et pour l’autre, la variabilité liée au caractère hétérogène et particulier de la géométrie des enceintes. Suite à la calibration des modèles de comportement thermo-hydro-mécanique, s’appuyant sur des résultats d'essais matériaux en laboratoire, une prévision à l'échelle de la structure, utilisant des chargements sur ouvrage reconstruits à partir de données historiques ou de synthèses fondées sur le retour d’expérience du parc nucléaire français exploité par EDF, est menée. L'historique de déformation est ainsi reconstruit en tout point de l'enceinte pour en déduire un modèle analytique des variations spatiales dans l'enceinte des déformations afin de faciliter son exploitation.En mettant en évidence de faibles écarts sur le long terme entre position à cœur et en surface, il est d'abord démontré que des chaînes de mesure de substitution sont utilisables avec un bon niveau de confiance en dehors des zones singulières. Par ailleurs, l’usage de fonctions de correction basées sur les transferts physiques régnant dans la structure permet de minimiser l’influence des changements de position de mesure dans l’épaisseur. / Extending the life of nuclear plant needs to ensure the sustainability of containment, third and last barrier against potential radioactive emissions. These structures are closely monitored from construction through a specific measuring system including strain gauges embedded in the concrete. Inherently reliable, though widely proven technology, these strain sensors start presenting significant failure rates, requiring the use of substitution measure chains. These solutions, which can only be placed near surface of concrete, allow only an indirect access to the initially observed strain. Doubts may therefore be issued on the representativeness of measurements and therefore on the quality of monitoring.To analyze the influence of the location of the measurement, the development and exploitation of two finite elements numerical models are proposed. This will be mainly for one to quantify the possible misbehavior in the long term and during tests related to temperature and moisture gradients; and for the other, the variability related to heterogeneous and particular nature of the geometry of nuclear power plant. As a result of the calibration of the thermo-hydro-mechanical behavior models, based on laboratory materials test, a forecast throughout the structure is conducted, using reconstructed loads from historical data or syntheses based on feedback from the french nuclear power plant operated by EDF. The history of deformation is thus reconstructed in any enclosure point to settle an analytical model of the spatial variations of the enclosure deformation to facilitate its operations.Highlighting small long-term differences between embedded and surface position, it is first shown that alternative measuring systems can be used with a high confident level outside singular position. Furthermore, the use of correction functions based on physical transfers prevailing in the structure makes it possible to minimize the influence of position measurement of changes in thickness.

A maioria dos líquidos encontrados na natureza são não-Newtonianos e o estudo do seu comportamento tem uma importância significante em diferentes áreas da engenharia. Entre eles, uma larga classe de materiais que exibem pequena ou nenhuma deformação quando sujeitos a um nível de tensões inferiores a uma tensão limite de escoamento – chamado de comportamento viscoplástico. A presente Dissertação tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos viscoplásticos não-lineares em uma cavidade forçada. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo viscoplástico recentemente introduzido por Souza Mendes e Dutra – SMD – e é aproximado por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão-extra, velocidade e pressão. As condições de compatibilidade entre os subespaços de elementos finitos para tensão-extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não-Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da UFRGS. Em diversos trabalhos encontrados na literatura, a superfície de escoamento do material é definida como a região onde o módulo da tensão-extra é igual à tensão limite de escoamento. É mostrado nesta Dissertação que esta metodologia pode conduzir à alguns erros, dado ao grande aumento experimentado pela taxa de cisalhamento em uma pequena faixa de tensões próximas à tensão limite de escoamento. Assim, foi adotada outra metodologia, definindo a superfície de escoamento como a linha onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade Newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o número de salto, J, o coeficiente de power-law, n, e a vazão adimensional, U*, são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos viscoplásticos. Os resultados obtidos estão de acordo com a literatura e atestam a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found on the nature and the study of their behavior has a significant importance on different areas of engineering. Among them, there is a wide class of materials that exhibits little or no deformation when subjected to a stress level behind an apparent yield stress – called the viscoplastic behavior. The present thesis aimed to a numerical study of two dimensional steady state laminar flows of non-linear viscoplastic fluids in a lid-driven cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the recently introduced Souza Mendes and Dutra – SMD – viscoplastic model and has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the extra-stress, velocity and pressure fields. In this way, the compatibility conditions between the extra-stressvelocity and pressure-velocity (Babuška-Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of UFRGS. In several works found on the literature, the yield surface of the material is defined as the region where the stress modulus is equal to the yield stress. Is shown in this work that this methodology can lead to some errors, due to the large strain rate increasing in a small range of values of stress on the vicinity of the yield stress. Therefore, it was adopted another approach, defining the yield surface as the line where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the jump number, J, the the power-law coefficient, n,and the non-dimensional flow rate, U*, are ranged in order to evaluate how they the influence on the viscoplastic fluid dynamics have been investigated. All results found were in accordance with the affine literature and attests the good stability features of the formulation.

La localisation de la déformation est un facteur essentiel à l’existence de la tectonique des plaques. Cependant, elle est le résultat d’interactions non-linéaires qui relient entre elles différentes échelles spatiales et temporelles que nous ne savons pas encore capturer dans un modèle rhéologique unique. Ce travail de thèse tente de combler en partie la fracture entre l’échelle pétrologique et celle de la mécanique à l’échelle crustale. Il cherche ainsi à mesurer l’impact des circulations de fluides sur la dynamique et la cinématique des dômes métamorphiques continentaux par le biais de modèles numériques d’échelle crus- tale. Pour répondre à cette question, nous avons cherché à établir un modèle empirique des écoulements de fluide et des interactions fluide-roche valide à l’échelle de la croûte. A cette échelle, il est important de capturer d’une part l’effet des circulations de fluides sur la densité qui détermine les forces de volume disponibles comme moteur de l’exhumation, mais aussi leur capacité à limiter ou accélérer la chute de résistance des bandes de cisaillement qui déterminent les forces de surface ralentissant l’exhumation. Ce modèle empirique est basé sur une observation au premier ordre du comportement métamorphique des roches le long d’un trajet P − T rétrograde : les roches métamorphiques de haut-grade (faciès amphibolite/granulite) sont exhumées vers la surface sans se rétromorphoser, excepté dans les zones où la déformation est localisée. Ce travail consiste donc à paramétrer cette observation de manière à pouvoir l’implémenter dans un code explicite thermomécanique 2D par le biais d’un couplage hydro-thermo-mécanique prenant en compte les échanges d’eau libre/eau liée avec la roche. Il introduit aussi la prise en compte simplifiée de l’hydrothermalisme lors de l’exhumation. Les résultats montrent que cette paramétrisation permet de mieux rendre compte des observations de terrain dans une gamme de paramètres validée par les mesures de laboratoire. Les résultats obtenus en limitant la rétromorphose des roches par la disponibilité de l’eau sont significativement différents de ceux obtenus dans les études publiées qui ne prenaient pas en compte ce paramètre. La thèse démontre notamment qu’il n’y a pas besoin d’hétérogénéités préexistantes pour obtenir des dômes asymétriques. Elle montre aussi que la chimie initiale de la roche peut exacerber ou diminuer la localisation de la déformation liée à la rétromorphose limitée par l’eau et contrôler ainsi la forme des dômes métamorphiques. / Strain localization is essential to the existence of plate tectonic on Earth. Yet, it is the result of non-linear interaction across several temporal and spatial scales and to date no unique rheological model exists that captures the effect of all these interactions at crustal scale. The work presented here is an attempt to fill part of the gap between out- crop scale and crustal scale models, between petrological scale and crustal scale. The central question is how to measure the effects of fluids and fluid-rock interactions on the dynamics and the kinematics of continental metamorphic core complexes by the mean of crustal scale numerical models. To answer this question, we have tried to esta- blish an empirical model of fluid flow and fluid-rock interaction, which could be valid at the scale of the crust. At that specific scale, it is important to capture the effects of free fluids both on the density of rocks, because it determines the body forces that are driving the exhumation of the crust, and on the strength of rocks and particularly on shear zones which constitute the ’hand-break’, the surface tractions, which resist this exhumation. This empirical model is based on a first order observation related to the behaviour of metamorphic rocks along the retrograde P − T path of their exhumation. High-grade metamorphic rocks (amphibolite/granulite) are indeed exhumed to the sur- face with very little retromorphose except those located within the shear bands where deformation localises. This works first consisted in parameterising this observation in a manner that could be implemented into a 2D thermo-mechanical code as an hydro- thermo-mechanical coupling which could account for free-water/rock-water exchange. In second part, a simplified scheme is introduced to account for hydrothermal cooling of the crust in the late stage of exhumation. The results show that this parameterization al- lows to better account for field observation in a range of parameters that are compatible with laboratory experiments. Limiting the retromorphose of rocks by the availability of free water is shown to produce significantly different metamorphic domes kinematics as compared to prior studies. We indeed demonstrate that it is not necessary to introduce heterogeneity in crustal composition to form asymmetric structures but that the effect of water is modulated by the chemistry of the protolith rocks which influences strain localization.