Submitted by Caio Agra (caioagra@uol.com.br) on 2018-03-23T19:11:09Z No. of bitstreams: 1 Dissertacao Revisada_Versão Final_Caio Agra.pdf: 2374099 bytes, checksum: b0931ab1ef450034ad264d3d1b2db229 (MD5) / Approved for entry into archive by NUBIA OLIVEIRA (nubia.marilia@ufba.br) on 2018-03-26T14:01:46Z (GMT) No. of bitstreams: 1 Dissertacao Revisada_Versão Final_Caio Agra.pdf: 2374099 bytes, checksum: b0931ab1ef450034ad264d3d1b2db229 (MD5) / Made available in DSpace on 2018-03-26T14:01:46Z (GMT). No. of bitstreams: 1 Dissertacao Revisada_Versão Final_Caio Agra.pdf: 2374099 bytes, checksum: b0931ab1ef450034ad264d3d1b2db229 (MD5) / Dois tipos de conhecimentos com características distintas, porém complementares, emergem no momento de planejar e desenvolver cursos de formação profissional: os conhecimentos teóricos e os procedimentais necessitam ser considerados como objetos de ensino sistematizado, dimensões de conteúdos específicos pertinentes a cada área curricular em particular, sendo a mecânica dos fluidos uma área de interesse no presente estudo. Com o objetivo de explorar a evolução dos conhecimentos teóricos e procedimentais desses conteúdos, realizamos nossa pesquisa de natureza quali-quantitativa, sendo os dados coletados por meio de questionários com itens sobre ‘escoamento’, ‘vazão’, ‘pressão’, massa específica’, ‘medição de fluidos’, ‘viscosidade’, ‘temperatura’ e ‘energia’, e por entrevistas semi-estruturadas, ambos instrumentos aplicados à turma do Curso de Formação de Técnicos de Petróleo e Gás da Petrobras, na cidade de Salvador/BA (2014). O primeiro artigo dessa dissertação multipaper aborda a Teoria das Habilidades Dinâmicas, a tipologia dos conhecimentos teóricos e procedimentais e apresenta os resultados, análises e discussões a partir de testes estatísticos não-paramétricos. O segundo artigo desenvolve, valida e apresenta resultados a partir de um sistema categórico para análise do conteúdo das entrevistas feitas com alunos do Curso de Formação. As análises não revelaram uma prevalência genérica dos conhecimentos teóricos em relação aos procedimentais, mas sim do conteúdo específico ‘vazão’ sobre os demais aqui estudados. / Two types of knowledge with distinct but complementary characteristics emerge when planning and developing professional qualification training courses: theoretical and procedural knowledge need to be considered as objects of systematized teaching, to the specific dimensions of contents pertinent to each particular curricular area, being the fluid mechanics an area of interest in the present study. With the objective of exploring the evolution of the theoretical and procedural knowledge of these contents, we conducted our qualitative and quantitative research, and the data were collected through questionnaires of items about 'flow', 'flow rate', 'pressure', specific mass', 'fluid measurement', 'viscosity', 'temperature' and 'energy', and by semi-structured interviews, both instruments applied to the Petrobras Gas and Oil Technicians Training Course in the city of Salvador/BA (2014). The first article of this multipaper dissertation addresses the Dynamic Skills Theory, from the typology of theoretical and procedural knowledge, and presents results, analyzes and discussions based on non-parametric statistical tests. The second article develops, validates and presents results from a categorical system for analyzing the content of interviews with students of the professional qualification training courses. The analyzes did not reveal a generic prevalence of theoretical knowledge in relation to procedural knowledge, but of the specific content 'flow rate' over the others studied here.

Orientador: José Luiz Gasche / Resumo: Originalmente, stents foram projetados para agir como barreiras mecânicas, impedindo a herniação de coils para a artéria-mãe. Recentemente, estudos mostraram que a atual geração de stents intracranianos auto-expansíveis altera a geometria local das artérias: um fenômeno com efeitos hemodinâmicos em parte incompreendidos. Nós realizamos simulações numéricas para avaliar a influência da remodelagem arterial induzida por stent sobre a hemodinâmica em aneurismas intracranianos. As simulações foram realizadas utilizando o software open-source OpenFOAM. O sangue foi modelado como fluido Newtoniano incompressível e as paredes arteriais foram consideradas rígidas. Para quantificar as alterações hemodinâmicas, avaliamos os parâmetros wall shear stress, WSS, e oscillatory shear index, OSI. Quatro geometrias reais de aneurismas intracranianos em bifurcações foram utilizadas. Em um aneurisma na bifurcação da artéria comunicante anterior (ACoA), um stent foi implantando, levando ao endireitamento das artérias que o receberam. Após o procedimento, os níveis de WSS e OSI aumentaram aproximadamente 60% e 25%, respectivamente. Em dois aneurismas em bifurcações da artéria cerebral média (MCA), dois stents foram implantados em uma configuração em “Y”, resultando em um endireitamento de ambas as artérias-filhas. O WSS máximo na superfície do aneurisma aumentou aproximadamente 5% em um dos casos e 22% no outro. Em outro aneurisma em uma bifurcação da MCA, um stent foi implantado, resultando no en... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Stents were first designed to act as mechanical barriers, preventing coil herniation into the parent artery. The current generation of self-expanding intracranial stents has recently been shown to change the local vascular geometry, a phenomenon with unclear hemodynamic effects. We carried out numerical simulations to assess the role of stent-induced vascular remodeling in modifying intraaneurysmal hemodynamics. Simulations were performed using the open-source software OpenFOAM. Blood was assumed to behave as an incompressible Newtonian fluid; vessel walls were assumed to be rigid. Wall shear stress, WSS, and oscillatory shear index, OSI, were evaluated to quantify the hemodynamic changes in the aneurysm sac. Four pre- and post-stent patient-specific geometries of intracranial bifurcation aneurysm were used. In one aneurysm at the anterior communicating artery (ACoA) bifurcation, a single stent was deployed, resulting in straightening of the host vessels. After stenting, WSS and OSI increased by approximately 60% and 25%, respectively. In two aneurysms at middle cerebral artery (MCA) bifurcations, two stents in a “Y” configuration were deployed, resulting in straightening of both daughter arteries. The maximum WSS on the aneurysm surface increased by approximately 5% in one case and 22% in the other. In another aneurysm at a bifurcation of the MCA, a single stent was deployed, resulting in straightening of the host vessels. After stenting, WSS and OSI reduced by approximately... (Complete abstract click electronic access below) / Mestre

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.

Fluidos Não-newtonianos estão presentes em grande quantidade na natureza e em muitas aplicações industriais. Tais fluidos apresentam uma variedade de comportamentos específicos. Dentre eles podemos citar: elasticidade, plasticidade, viscosidade variando com a taxa de deformação, tixotropia. Esta tese de doutorado apresenta um recente modelo para Fluidos Estruturados Tixotrópicos (de Souza Mendes, 2009), no qual baseasse nas tradicionais equações de balanço de momentum e massa, em uma equação constitutiva de Maxwell-b modificada e uma equação de evolução para o nível de estrutura do material. Na equação constitutiva de Maxwell-b modificada, tempo de relaxação e viscosidade estrutural são funções da estrutura do material. A equação evolutiva para o parâmetro de estruturação prevê efeitos de tixotropia, dependentes do tempo de equilíbrio do material. As equações que constituem o modelo mecânico para Fluidos Estruturados Tixotrópicos são aproximadas numericamente através do Método de Elementos Finitos de Galerkin Mínimos-quadrados (GLS). Tal método adicionou termos dependentes de malha ao método clássico de Galerkin, de forma a aumentar a estabilidade da formulação original, proporcionando compatibilizar os sub-espaços de velocidade e pressão. São investigados escoamentos creeping flow de Fluidos Estruturados Tixotrópicos em dois problemas clássicos de engenharia: em torno de um cilindro confinado entre duas placas e em uma expansão planar abrupta. Parâmetros relevantes ao modelo são variados em faixas pertinentes ao problema físico simulado e os resultados demonstraram uma boa capacidade de previsão. / Non-Newtonian fluids are present in large quantities in nature and in many industrial applications. Such fluids exhibit a variety of specific behaviors. Among them we can mention: elasticity, plasticity, viscosity varying with the strain rate, thixotropy, etc. This thesis presents a new model for thixotropic Structured Fluids (de Souza Mendes, 2009), in which were a basis in traditional balance equations of mass and momentum, the constitutive equation of Maxwell-type and an evolution equation for the level structure of the material. In the modified constitutive equation of Maxwell-b, relaxation time and structural viscosity are functions of material structure. The structure parameter evolutive equation provides effects of thixotropy, dependent on the equilibration time of the material. The equations which constitute the mechanical model for Structured Thixotropic fluids are approximated numerically using the Galerkin Least-squares Finite Element Method (GLS). Such a method added mesh dependent terms to the classical method of Galerkin in order to increase the stability of the original formulation, providing compatibility between sub-spaces of velocity and pressure. Are investigated Structured Thixotropic Fluids creeping flow in two classical problems in engineering: around a cylinder confined between two plates and in a planar abrupt expansion. Relevant parameters to the model are varied in pertinent ranges to physical problem simulated and the results demonstrated a good predictive capability.

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnologico / Made available in DSpace on 2016-01-08T21:07:33Z (GMT). No. of bitstreams: 1 103916.pdf: 3688429 bytes, checksum: 7eb414d2824d20a2c8544d4868002676 (MD5) Previous issue date: 1996 / Estudo dos fundamentos dos esquemas de alta resolução TVD e NVD. Revisa os principais aspectos matemáticos de leis de conservação hiperbólicas cujas soluções podem envolver descontinuidades. As características dissipativas e oscilatórias dos principais esquemas numéricos, usados no desenvolvimento dos esquemas de alta resolução, são analisadas e testadas na equação de Burgers invícida. Os esquemas TVD são testados para as Eq. de Euler nos problemas de tubo de choque e bocal divergente. Os esquemas NVD são testados para as equações de Navier-Stokes no problema de convecção forçada em uma cavidade quadrada. É utilizado o método de volumes finitos e o sistema de equações é resolvido segregadamente ou via características locais (simultâneamente). Os esquemas NVD foram utilizados para gerar esquemas de interpolação de alta resolução, que foram implementados num procedimento geral com o uso de limitadores. Os resultados obtidos são comparados com soluções analíticas, benchmark e entre si.

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.

Fluidos Não-newtonianos estão presentes em grande quantidade na natureza e em muitas aplicações industriais. Tais fluidos apresentam uma variedade de comportamentos específicos. Dentre eles podemos citar: elasticidade, plasticidade, viscosidade variando com a taxa de deformação, tixotropia. Esta tese de doutorado apresenta um recente modelo para Fluidos Estruturados Tixotrópicos (de Souza Mendes, 2009), no qual baseasse nas tradicionais equações de balanço de momentum e massa, em uma equação constitutiva de Maxwell-b modificada e uma equação de evolução para o nível de estrutura do material. Na equação constitutiva de Maxwell-b modificada, tempo de relaxação e viscosidade estrutural são funções da estrutura do material. A equação evolutiva para o parâmetro de estruturação prevê efeitos de tixotropia, dependentes do tempo de equilíbrio do material. As equações que constituem o modelo mecânico para Fluidos Estruturados Tixotrópicos são aproximadas numericamente através do Método de Elementos Finitos de Galerkin Mínimos-quadrados (GLS). Tal método adicionou termos dependentes de malha ao método clássico de Galerkin, de forma a aumentar a estabilidade da formulação original, proporcionando compatibilizar os sub-espaços de velocidade e pressão. São investigados escoamentos creeping flow de Fluidos Estruturados Tixotrópicos em dois problemas clássicos de engenharia: em torno de um cilindro confinado entre duas placas e em uma expansão planar abrupta. Parâmetros relevantes ao modelo são variados em faixas pertinentes ao problema físico simulado e os resultados demonstraram uma boa capacidade de previsão. / Non-Newtonian fluids are present in large quantities in nature and in many industrial applications. Such fluids exhibit a variety of specific behaviors. Among them we can mention: elasticity, plasticity, viscosity varying with the strain rate, thixotropy, etc. This thesis presents a new model for thixotropic Structured Fluids (de Souza Mendes, 2009), in which were a basis in traditional balance equations of mass and momentum, the constitutive equation of Maxwell-type and an evolution equation for the level structure of the material. In the modified constitutive equation of Maxwell-b, relaxation time and structural viscosity are functions of material structure. The structure parameter evolutive equation provides effects of thixotropy, dependent on the equilibration time of the material. The equations which constitute the mechanical model for Structured Thixotropic fluids are approximated numerically using the Galerkin Least-squares Finite Element Method (GLS). Such a method added mesh dependent terms to the classical method of Galerkin in order to increase the stability of the original formulation, providing compatibility between sub-spaces of velocity and pressure. Are investigated Structured Thixotropic Fluids creeping flow in two classical problems in engineering: around a cylinder confined between two plates and in a planar abrupt expansion. Relevant parameters to the model are varied in pertinent ranges to physical problem simulated and the results demonstrated a good predictive capability.

Dispositivos de reação miniaturizados tendem a ganhar espaço na indústria de processos químicos por elevarem o transporte de massa e de calor e a segurança dos processos. Para compreender o papel de cada elemento constituinte de um microrreator sobre seu campo de velocidades e fenômenos de mistura, foram simulados um conjunto de dois canais, uma junção em T, 30 canais em formato de serpentina e um microchip completo através da fluidodinâmica computacional. A seção transversal destes microdispositivos têm dimensões entre 100 e 300 µm, enquanto o comprimento dos canais varia de 3000 a 25190 µm. Os modelos computacionais foram discretizados por malhas hexaédricas e os campos de velocidade em estado estacionário foram calculados para vazões de alimentação entre 12,5 e 2000 µL min-1, considerando regime laminar. A mistura foi avaliada pela injeção de traçadores não-reativos e distribuição das respectivas frações mássicas. As simulações foram validadas usando microvelocimetria por imagens de partículas. Os campos de velocidade possuem magnitudes significativas apesar das dimensões reduzidas e baixas vazões de operação dos sistemas. As imagens experimentais do escoamento evidenciaram o formato parabólico do campo de velocidades e o deslocamento de seu ponto máximo nas regiões curvas causado pela força centrífuga, como estimado pelo modelo computacional. Tal força, em conjunto com as forças viscosas na parede, gera fluxos secundários no escoamento. A distribuição de traçadores não-reativos evidenciou a importância dos fluxos secundários para promover mistura na direção ortogonal ao escoamento principal, ocorrendo sob o regime estacionário nas vazões analisadas. O estudo aqui realizado evidencia o emprego da fluidodinâmica computacional como ferramenta para melhor compreensão da fluidodinâmica e como apoio ao projeto de microdispositivos. / Miniaturized reaction vessels are drawing attention of chemical industries because they promote better mass and heat transfer and also enhance process safety. To understand the relevance of each element of a microreactor on the velocity field of the equipment and the corresponding mixing processes, several microdevices were simulated using computational fluid dynamics: an assembly of two channels, a T-junction, 30 channels in a serpentine assembly and a full microreactor. The cross section of the devices is 100 - 300 µm wide and the length of the channels varies between 3000 and 25190 µm. Computational domains were discretized using hexahedral meshes and steady-state velocity fields were computed considering laminar flow for flow rates between 12,5 and 2000 µL min-1. Mixing was evaluated by injecting inert tracers and monitoring their distribution. Simulations were validated against experimental micro particle image velocimetry data. Velocities throughout the devices are relatively high despite the small dimensions of the cross sections and small flow rates. Experimental images of the flow elucidated the parabolic shape of the velocity profile and its distortion on curved segments caused by centrifugal forces, matching predictions of the computational model. Tracer maps indicated secondary flows play an important role in mixing stream perpendicular to the main flow direction. This study emphasizes the use of computational fluid dynamics as a tool for understating flow throughout microdevices and supporting their design.

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.