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 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.

Este trabalho aborda o estudo de modelos constitutivos a fim de analisar o comportamento mecânico de um biopolímero derivado de óleo de mamona (Ricinus communis). A importância deste trabalho se dá ao fato de: servir como roteiro de ensaio de materiais poliméricos; completar caracterização do material estudado além do regime viscoelástico, bem como aplicação da análise do comportamento do material além do regime viscoelástico linear. A aplicação dos modelos avaliados via Método dos Elementos Finitos (MEF) exigiu o levantamento das propriedades e do comportamento mecânico do biopolímero, bem como a preparação dos ensaios, com as devidas especificações sugeridas pela ASTM (American Society for Testing and Materials) para cada tipo de ensaio. Assim, foram obtidos os corposde- prova (CDP) para a realização dos ensaios mecânicos quase-estáticos de tração, compressão e flexão (monotônicos e cíclicos), bem como, ensaio para análise dinâmicomecânica (Dynamic Mechanical Analysis - DMA) para a obtenção de propriedades associada às parcelas viscosas do material. Com estes foram obtidos parâmetros inerentes a cada tipo de ensaio, permeados com uma cuidadosa revisão bibliográfica sobre as implicações para a realização de cada ensaio. Uma outra etapa foi a verificação de modelos de plastificação baseados na teoria de von Mises, Drucker-Prager e de viscoelasticidade implementados no programa comercial ABAQUS®. Sendo que o objetivo maior foi a comparação dos resultados experimentais obtidos com os modelos em MEF. Os resultados para tração monotônica obtiveram erro de 0,61%, os resultados de compressão monotônica apresentaram erro de 1,5% e as simulações para os carregamentos cíclicos de tração e compressão não foram tão bons quanto os anteriores. Por isso se justifica a criação do modelo fenomenológico utilizando os parâmetros do material que foram identificados. / This work develops the study involving constitutive models with the objective to analyze the mechanical behavior of a biopolymer obtained from the Castor Oil Polyurethane (Ricinus communis). The importance of this work is due to the fact that: serves as a test script for the use of polymer materials; complete the characterization of the material studied beyond the viscoelastic realm, as well as the application of the analysis of the behavior of the material beyond the linear viscoelastic domain. The application of the evaluated models by way of the Finite Element Method (FEM), required the gathering of the properties and mechanical behavior of the biopolymer, as well as the preparation of the tests, with the required specifications suggested by the ASTM (American Society for Testing and Materials) for each type of test. Therefore, there were obtained the body-tests in order to run the mechanical quasi-static traction tests, tensile, compression and flexion (monotonic and cyclical) dynamicmechanical analysis (DMA) in order to obtain the properties associated with the viscous parts of the material. With these the inherent parameters were obtained for each specimen, permeated with a careful biographical revision of the implications for the experimenting of each specimen. Another phase was the verification of laminated models based on the theory by Von Mises, Drucker-Prager and of the viscoelasticity implemented in the commercial program ABAQUS®. The main objective being the comparison of the experimental results obtained with the models in FEM. The results for the monotonic traction obtained a error of 0,61%, the results of monotonic comparison presented an error of 1,5% and the simulations for the cyclical charges of traction and compression were not as good as the former. Therefore it justifies the creation of the phenomenological model utilizing the parameters of the materials that were identified.

Made available in DSpace on 2015-03-04T18:50:19Z (GMT). No. of bitstreams: 1 tese cristiane.pdf: 2779864 bytes, checksum: 78fe9d6461cd9dd3cd23fe85acdd548c (MD5) Previous issue date: 2010-08-30 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Nesta tese apresenta-se uma nova formulação mista estabilizada regularizada de elementos finitos nas variáveis primais com interpolação contínua para a velocidade e descontínua para a pressão, para problemas de fluidos viscoplásticos.Esta formulação tem base nos métodos de lagrangeano aumentado-regularização e de estabilização por mínimos quadrados.Discutem-se as dificuldades de resolução dos problemas viscoplásticos especialmente aqueles modelados pela relação de Bingham que prevê uma descontinuidade na forma de restrição de desigualdade. São apresentadas criticamente a teoria de lagrangeano aumentado, os modelos regularizados e o método estabilizado via mínimos quadrados. São feitas análises matemáticas para a nova formulação em termos de estabilidade, existência e unicidade de solução. Ordens de convergência ótimas são obtidas matematicamente, superando os métodos clássicos que também apresentam limitação para o valor de tensão limite. Resultados numéricos são apresentados confirmando a teoria aqui desenvolvida e mostrando a robustez do novo método para resolver problemas onde a tensão limite é muito elevada.

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.

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.

Made available in DSpace on 2016-12-12T20:25:11Z (GMT). No. of bitstreams: 1 Thiago Andre Carniel.pdf: 5594070 bytes, checksum: 5e24f59f581ec79c2edb07d9afdd502e (MD5) Previous issue date: 2013-07-12 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This work presents the development of a numerical formulation for the transient analysis of spatial truss structures considering geometrical and material nonlinearities. The equilibrium problem is formulated using the Principle of the Virtual Power in a current Lagrangian description. The displacement and velocity fields are discretized employing the Finite Element Method and the nonlinear equilibrium equations are solved by the full Newton- Raphson procedure. In order to obtain quadratic convergence, the algorithmic tangent matrix expression is developed and presented. Next, the transient problem is formulated and solved using the Newmark integration algorithm. The material model adopted is viscoelastic-viscoplastic coupled to Lemaitre s damage. The viscoelastic behavior is represented by the generalized Kelvin-Voigt rheological model and the viscoplastic contribution is introduced via Perzyna s equation. The evolution laws of the material model are integrated numerically using the implicit Euler scheme, resulting in the corresponding return mapping. The algorithmic tangent modulus for this material model is derived and validated by means of numerical differentiation using the complex variable approach. The constitutive parameters for the model are obtained by minimization of the difference between experimental and numerical curves, aiming at polymeric materials. Since non-uniqueness is expected for the minimization problem, the gradient free Particle Swarm Optimization method (PSO) is adopted for its solution. Finally, the results obtained and some particularities of the proposed model are discussed. / Este trabalho tem por objetivo desenvolver uma formulação numérica para a análise transiente de estruturas treliçadas tridimensionais considerando não linearidade geométrica e material. O problema de equilíbrio é apresentado através do Princípio das Potências Virtuais segundo uma formulação Lagrangeana Corrente. O campo de deslocamentos e velocidades é discretizado através do Método dos Elementos Finitos. As equações não lineares de equilíbrio são resolvidas pelo método de Newton-Raphson, sendo apresentado o desenvolvimento da matriz de rigidez tangente consistente. Posteriormente, o problema transiente é formulado, sendo as equações do movimento solucionadas pelo método de integração implícito de Newmark. O modelo constitutivo proposto é viscoelástico-viscoplástico acoplado com o modelo de dano de Lemaitre. A formulação viscoelástica é representada pelo modelo generalizado de Kelvin-Voigt e a contribuição viscoplástica é dada pela equação de Perzyna. As leis de evolução do modelo material são integradas numericamente através do algoritmo implícito de Euler, resultando no algoritmo de retorno. Para obter a taxa de convergência quadrática do método iterativo de Newton-Raphson durante o equilíbrio estrutural, é apresentado o cálculo analítico do módulo tangente algoritmo consistente, o qual é validado através do método da derivada complexa. Os parâmetros materiais são encontrados através do ajuste entre as curvas experimental e teórica, visando a aplicação para materiais poliméricos. Para o ajuste das curvas utilizou-se o método de Otimização por Nuvem de Partículas Particle Swarm Optimization (PSO). Ao final é realizada uma discussão dos resultados encontrados e comentado as particularidades do modelo proposto.

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.

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.

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.