O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual 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. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual 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. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual 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. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

O estudo de escoamentos de fluidos não-Newtonianos através de expansões desperta um grande interesse em pesquisadores nas diversas áreas da engenharia, devido a sua ampla aplicação em indústrias e no meio acadêmico. O objetivo principal desta Dissertação é simular problemas de escoamentos envolvendo fluidos viscoplásticos através de expansões axissimétricas abruptas. O modelo mecânico empregado é baseado nas equações de conservação de massa e de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), com a função de viscosidade de Herschel-Bulkley regularizada pela equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados, ou Galerkin Least-squares (GLS). Esse método (GLS) é usado a fim superar as dificuldades numéricas do modelo de Galerkin clássico: a condição de Babuška-Brezzi e a instabilidade inerente em regiões advectivas do escoamento. O método é construído adicionando termos de malha-dependentes a fim aumentar a estabilidade da formulação de Galerkin clássica sem danificar sua consistência. A formulação GLS é aplicada para estudar a influência do índice power-law, da tensão limite de cisalhamento e razão de aspecto na dinâmica do escoamento de fluidos de Herschel-Bulkley através de expansões axissimétricas abruptas de razão de aspecto 1:2 e 1:4. Os problemas que envolvem números de Reynolds desprezíveis, para uma escala do número de Herschel-Bulkley entre 0 e 100 e índice de comportamento entre 0,2 e 1,0 são apresentados. Os resultados são fisicamente detalhados e estão de acordo com a literatura. / The study of non-Newtonian fluid flows in expansions is of great interest for researchers in the several branches of engineering, due to their wide application both in industry and academy. The objective of this Dissertation is to simulate flow problems involving a viscoplastic fluid through an axisymmetric abrupt expansion. The mechanical model employed is based on the mass and momentum conservative equations for isochoric flows coupled with the Generalized Newtonian Liquid (GNL) constitutive equation, with the Papanastasiou-regularized Herschel-Bulkley viscosity function. The mechanical model is approximated by a stabilized finite element scheme, namely the Galerkin Least-squares method. This method (GLS) is used in order to overcome the numerical difficulties of the classical Galerkin method: the Babuška- Brezzi condition and the inherent instability in advective flow regions. The method is built adding mesh-dependent terms in order to increase the stability of the classical Galerkin formulation without damaging its consistency. The GLS formulation is applied to study the influence of power-law index, yield stress and aspect reason in the flow dynamics of Herschel- Bulkley fluids through an axisymmetric abrupt expansions of aspect reason 1:2 and 1:4. Problems involving negligible Reynolds numbers, for a Herschel-Bulkley number range between 0 and 100 and e power-law index range between 0.2 and 1.0 are presented. The results are physically comprehensive and are in accordance with the literature.

A grande maioria dos fluidos encontrados na natureza se comportam como fluidos não- Newtonianos o que torna o seu estudo muito importante para diversas áreas da engenharia. Este trabalho tem como objetivo simular o problema específico de escoamentos de fluidos viscoplásticos em expansões abruptas planares com razão de aspecto de 1:4. O problema em questão se mostra interessante pois em muitos sistemas industriais são apresentados a geometria proposta para estudo associada a fluidos viscoplásticos. Foi empregado o modelo mecânico multi-campos (mult-field) baseado nas equações de conservação de massa e balanço de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), associada à função de viscosidade de Herschel-Bulkley regularizada através da equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados (GLS). A fim de se pesquisar os fenômenos reológicos ali presentes é feito o estudo da influência do índice de power-law na topologia de um escoamento creeping flow (Re@0) para uma vasta faixa de números de Herschel-Bulkley variando entre 0.1 e 100. Os resultados mostraram-se satisfátorios, apresentando uma forte influencia do número de Herschel-Bulkley e do índice de power law na topologia e na dinâmica do escoamento. Obteve-se uma validação do trabalho comparando-se os resultados obtidos nesta dissertação com os obtidos em artigo de grande credibilidade e aceito por toda a academia científica. / The majority of fluids found in nature behave like non-Newtonian fluids what makes their study of great importance to various areas of engineering. The present work aims to simulate the specific problem of the viscoplastic fluids flows through in planar abrupt expansion with the aspect ratio of 1:4. The problem in question is interesting because in many industrial systems are the proposal to study geometry associated with viscoplastic fluid. It employed the multi-field mechanical model based on equations of conservation of mass and momentum balance for the isochoric flow coupled with the constitutive equation of the Generalized Newtonian Liquids (GNL), associated with the function of viscosity of Herschel-Bulkley regularized by the equation of Papanastasiou. The mechanical model is approximated by a stabilized finite element model, called Galerkin Least-Squares method (GLS). In order to research the rheological phenomena present is done studying the influence of power-law index in the topology of the creeping flow (Re@0) for a wide range of numbers of Herschel-Bulkley ranging between 0.1 and 100. The results have proved satisfactory, showing a strong influence on the number of Herschel-Bulkley and power law index of the topology and dynamics of the flow. Got to be a validation of the work comparing the results with those obtained in this thesis in an article of great credibility and accepted throughout the scientific academy.

O estudo de escoamentos de fluidos não-Newtonianos através de expansões desperta um grande interesse em pesquisadores nas diversas áreas da engenharia, devido a sua ampla aplicação em indústrias e no meio acadêmico. O objetivo principal desta Dissertação é simular problemas de escoamentos envolvendo fluidos viscoplásticos através de expansões axissimétricas abruptas. O modelo mecânico empregado é baseado nas equações de conservação de massa e de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), com a função de viscosidade de Herschel-Bulkley regularizada pela equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados, ou Galerkin Least-squares (GLS). Esse método (GLS) é usado a fim superar as dificuldades numéricas do modelo de Galerkin clássico: a condição de Babuška-Brezzi e a instabilidade inerente em regiões advectivas do escoamento. O método é construído adicionando termos de malha-dependentes a fim aumentar a estabilidade da formulação de Galerkin clássica sem danificar sua consistência. A formulação GLS é aplicada para estudar a influência do índice power-law, da tensão limite de cisalhamento e razão de aspecto na dinâmica do escoamento de fluidos de Herschel-Bulkley através de expansões axissimétricas abruptas de razão de aspecto 1:2 e 1:4. Os problemas que envolvem números de Reynolds desprezíveis, para uma escala do número de Herschel-Bulkley entre 0 e 100 e índice de comportamento entre 0,2 e 1,0 são apresentados. Os resultados são fisicamente detalhados e estão de acordo com a literatura. / The study of non-Newtonian fluid flows in expansions is of great interest for researchers in the several branches of engineering, due to their wide application both in industry and academy. The objective of this Dissertation is to simulate flow problems involving a viscoplastic fluid through an axisymmetric abrupt expansion. The mechanical model employed is based on the mass and momentum conservative equations for isochoric flows coupled with the Generalized Newtonian Liquid (GNL) constitutive equation, with the Papanastasiou-regularized Herschel-Bulkley viscosity function. The mechanical model is approximated by a stabilized finite element scheme, namely the Galerkin Least-squares method. This method (GLS) is used in order to overcome the numerical difficulties of the classical Galerkin method: the Babuška- Brezzi condition and the inherent instability in advective flow regions. The method is built adding mesh-dependent terms in order to increase the stability of the classical Galerkin formulation without damaging its consistency. The GLS formulation is applied to study the influence of power-law index, yield stress and aspect reason in the flow dynamics of Herschel- Bulkley fluids through an axisymmetric abrupt expansions of aspect reason 1:2 and 1:4. Problems involving negligible Reynolds numbers, for a Herschel-Bulkley number range between 0 and 100 and e power-law index range between 0.2 and 1.0 are presented. The results are physically comprehensive and are in accordance with the literature.

A grande maioria dos fluidos encontrados na natureza se comportam como fluidos não- Newtonianos o que torna o seu estudo muito importante para diversas áreas da engenharia. Este trabalho tem como objetivo simular o problema específico de escoamentos de fluidos viscoplásticos em expansões abruptas planares com razão de aspecto de 1:4. O problema em questão se mostra interessante pois em muitos sistemas industriais são apresentados a geometria proposta para estudo associada a fluidos viscoplásticos. Foi empregado o modelo mecânico multi-campos (mult-field) baseado nas equações de conservação de massa e balanço de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), associada à função de viscosidade de Herschel-Bulkley regularizada através da equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados (GLS). A fim de se pesquisar os fenômenos reológicos ali presentes é feito o estudo da influência do índice de power-law na topologia de um escoamento creeping flow (Re@0) para uma vasta faixa de números de Herschel-Bulkley variando entre 0.1 e 100. Os resultados mostraram-se satisfátorios, apresentando uma forte influencia do número de Herschel-Bulkley e do índice de power law na topologia e na dinâmica do escoamento. Obteve-se uma validação do trabalho comparando-se os resultados obtidos nesta dissertação com os obtidos em artigo de grande credibilidade e aceito por toda a academia científica. / The majority of fluids found in nature behave like non-Newtonian fluids what makes their study of great importance to various areas of engineering. The present work aims to simulate the specific problem of the viscoplastic fluids flows through in planar abrupt expansion with the aspect ratio of 1:4. The problem in question is interesting because in many industrial systems are the proposal to study geometry associated with viscoplastic fluid. It employed the multi-field mechanical model based on equations of conservation of mass and momentum balance for the isochoric flow coupled with the constitutive equation of the Generalized Newtonian Liquids (GNL), associated with the function of viscosity of Herschel-Bulkley regularized by the equation of Papanastasiou. The mechanical model is approximated by a stabilized finite element model, called Galerkin Least-Squares method (GLS). In order to research the rheological phenomena present is done studying the influence of power-law index in the topology of the creeping flow (Re@0) for a wide range of numbers of Herschel-Bulkley ranging between 0.1 and 100. The results have proved satisfactory, showing a strong influence on the number of Herschel-Bulkley and power law index of the topology and dynamics of the flow. Got to be a validation of the work comparing the results with those obtained in this thesis in an article of great credibility and accepted throughout the scientific academy.

O estudo de escoamentos de fluidos não-Newtonianos através de expansões desperta um grande interesse em pesquisadores nas diversas áreas da engenharia, devido a sua ampla aplicação em indústrias e no meio acadêmico. O objetivo principal desta Dissertação é simular problemas de escoamentos envolvendo fluidos viscoplásticos através de expansões axissimétricas abruptas. O modelo mecânico empregado é baseado nas equações de conservação de massa e de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), com a função de viscosidade de Herschel-Bulkley regularizada pela equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados, ou Galerkin Least-squares (GLS). Esse método (GLS) é usado a fim superar as dificuldades numéricas do modelo de Galerkin clássico: a condição de Babuška-Brezzi e a instabilidade inerente em regiões advectivas do escoamento. O método é construído adicionando termos de malha-dependentes a fim aumentar a estabilidade da formulação de Galerkin clássica sem danificar sua consistência. A formulação GLS é aplicada para estudar a influência do índice power-law, da tensão limite de cisalhamento e razão de aspecto na dinâmica do escoamento de fluidos de Herschel-Bulkley através de expansões axissimétricas abruptas de razão de aspecto 1:2 e 1:4. Os problemas que envolvem números de Reynolds desprezíveis, para uma escala do número de Herschel-Bulkley entre 0 e 100 e índice de comportamento entre 0,2 e 1,0 são apresentados. Os resultados são fisicamente detalhados e estão de acordo com a literatura. / The study of non-Newtonian fluid flows in expansions is of great interest for researchers in the several branches of engineering, due to their wide application both in industry and academy. The objective of this Dissertation is to simulate flow problems involving a viscoplastic fluid through an axisymmetric abrupt expansion. The mechanical model employed is based on the mass and momentum conservative equations for isochoric flows coupled with the Generalized Newtonian Liquid (GNL) constitutive equation, with the Papanastasiou-regularized Herschel-Bulkley viscosity function. The mechanical model is approximated by a stabilized finite element scheme, namely the Galerkin Least-squares method. This method (GLS) is used in order to overcome the numerical difficulties of the classical Galerkin method: the Babuška- Brezzi condition and the inherent instability in advective flow regions. The method is built adding mesh-dependent terms in order to increase the stability of the classical Galerkin formulation without damaging its consistency. The GLS formulation is applied to study the influence of power-law index, yield stress and aspect reason in the flow dynamics of Herschel- Bulkley fluids through an axisymmetric abrupt expansions of aspect reason 1:2 and 1:4. Problems involving negligible Reynolds numbers, for a Herschel-Bulkley number range between 0 and 100 and e power-law index range between 0.2 and 1.0 are presented. The results are physically comprehensive and are in accordance with the literature.

A grande maioria dos fluidos encontrados na natureza se comportam como fluidos não- Newtonianos o que torna o seu estudo muito importante para diversas áreas da engenharia. Este trabalho tem como objetivo simular o problema específico de escoamentos de fluidos viscoplásticos em expansões abruptas planares com razão de aspecto de 1:4. O problema em questão se mostra interessante pois em muitos sistemas industriais são apresentados a geometria proposta para estudo associada a fluidos viscoplásticos. Foi empregado o modelo mecânico multi-campos (mult-field) baseado nas equações de conservação de massa e balanço de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), associada à função de viscosidade de Herschel-Bulkley regularizada através da equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados (GLS). A fim de se pesquisar os fenômenos reológicos ali presentes é feito o estudo da influência do índice de power-law na topologia de um escoamento creeping flow (Re@0) para uma vasta faixa de números de Herschel-Bulkley variando entre 0.1 e 100. Os resultados mostraram-se satisfátorios, apresentando uma forte influencia do número de Herschel-Bulkley e do índice de power law na topologia e na dinâmica do escoamento. Obteve-se uma validação do trabalho comparando-se os resultados obtidos nesta dissertação com os obtidos em artigo de grande credibilidade e aceito por toda a academia científica. / The majority of fluids found in nature behave like non-Newtonian fluids what makes their study of great importance to various areas of engineering. The present work aims to simulate the specific problem of the viscoplastic fluids flows through in planar abrupt expansion with the aspect ratio of 1:4. The problem in question is interesting because in many industrial systems are the proposal to study geometry associated with viscoplastic fluid. It employed the multi-field mechanical model based on equations of conservation of mass and momentum balance for the isochoric flow coupled with the constitutive equation of the Generalized Newtonian Liquids (GNL), associated with the function of viscosity of Herschel-Bulkley regularized by the equation of Papanastasiou. The mechanical model is approximated by a stabilized finite element model, called Galerkin Least-Squares method (GLS). In order to research the rheological phenomena present is done studying the influence of power-law index in the topology of the creeping flow (Re@0) for a wide range of numbers of Herschel-Bulkley ranging between 0.1 and 100. The results have proved satisfactory, showing a strong influence on the number of Herschel-Bulkley and power law index of the topology and dynamics of the flow. Got to be a validation of the work comparing the results with those obtained in this thesis in an article of great credibility and accepted throughout the scientific academy.

Esta dissertação investiga numericamente a influência dos efeitos viscosos e elásticos em escoamentos de materiais viscoplásticos no interior de uma cavidade dirigida. O modelo mecânico empregado é constituído pelas equações de conservação de massa e pelo princípio da quantidade de movimento linear, para fluidos incompressíveis, acoplado à equação constitutiva. Esta equação modifica o modelo viscoelástico de Oldroyd-B de modo a acomodar que os tempos de relaxação e retardo do material, bem como sua viscosidade viscoplástica, dependam das mudanças de sua microestrutura. A aproximação numérica do modelo emprega o método multi-campos de Galerkin mínimos-quadrados em termos do tensor de tensão extra, do vetor velocidade e do campo de pressão. Os resultados objetivam a determinção do tamanho e localização das regiões aparentemente não-escoadas do material, bem como sua deformação elástica, intensidade de tensão, e a sua vorticidade no interior da cavidade. Os resultados claramente indicam que o padrão do escoamento é fortemente influenciado pela variação dos efeitos elásticos (variação do tempo de relaxação adimensional, θ0 * ), viscosos (variação do índice de power-law, n) e cinemáticos (variação da velocidade adimensional, U* , do escoamento) no interior da cavidade. / This dissertation investigated numerically the influence of viscous and elastic effects on flows of viscoplastic materials within a lid-driven cavity. The mechanical model used is made up of mass and momentum balance equations, coupled with the constitutive equation. This equation modifies the viscoelastic Oldroyd-B model to accommodate both relaxation and retardation times, and viscosity function, dependent on the microstructure changes. Numerical approximations of the model make use a three-field Galerkin least squares method in terms of the extra stress tensor, velocity vector and pressure field. Computations focus on the determination of the size and position of apparently unyielded regions as well as the elastic deformation, stress intensity, and the vorticity within of the cavity. Results clearly indicate that the flow pattern is strongly influenced by the elastic (variation of the dimensionless relaxation time, θ0 * ), viscous (variation of the power-law index, n) and kinematic (variation of the dimensionless flow velocity, U* ) effects within the cavity.