No description available.

The existence of hydrodynamically developed, laminar Viscoelastic Secondary Flows (VSFs) of non-Newtonian fluids in straight ducts of non-circular cross-section was proposed in the 1950's. VSFs have since been observed sporadically, and only once with a velocimetric technique. Using axial and transverse full flow-field velocity-position raster maps made with Rheological Nuclear Magnetic Resonance (Rheo-NMR), Newtonian and non-Newtonian fluid flows were quantified in Hagen-Poiseuille and Power Law contexts, over more than two orders of magnitude of flow rate, in ducts of circle, square, triangle, and pentagon cross-section. VSF was reliably and repeatedly observed to occur at between one part in 130 and one part in 600 of the primary axial flow velocity. Velocity measurements ranged from <10 µm/s to approximately 30 cm/s, suggesting a velocity dynamic range >3E4 without optimization. To obtain VSF flow direction information, a novel flow directional phantom was developed and characterized. Aqueous solutions of Polyethylene Oxide (PEO), Viscarin GP-109NF, Viscarin GP-209NF (V209), Hyaluronan (HA) in a Phosphate-Buffered Saline-like solvent, and an aqueous Polyethylene Glycol/PEO-based Boger fluid were investigated. Axial data was corroborated with related data gathered by an independent method. Basic simulations corroborated the VSF observations. Duct hydraulic diameters (>= 1.6 mm) approached the micro-channel regime. VSF detections in HA --- synovial fluid's principal component --- and V209 were novel, as were observations of some artifacts which were subsequently characterized and corrected. The detection of VSF in HA represents the first experimental evidence suggesting that its second normal stress (N_2) is comparable to that of better-characterized fluids. In the first application of a new VSF-based method, a particular Boger fluid's constant viscosity and, in the square duct, its lack of VSF were used with established criteria to suggest that the fluid's N_2 approached zero. The development of a rudimentary, but versatile and inexpensive home-built velocimetric spectrometer is detailed, as are several new components. An exhaustive VSF literature review is included. The remarkable transverse velocimetric ability of Rheo-NMR in both optically opaque and transparent system is highlighted, suggesting that perhaps the technique might represent, in both micro-channels and conventional ducts, the gold-standard in flow velocimetry.

Cette thèse est composée de deux parties indépendantes. La première est consacrée à l'étude de quelques problèmes elliptiques de type de Kirchhoff de la forme suivante : -M(ʃΩNul² dx) Δu = f(x, u) xЄΩ ; u(x) = o xЄƋΩ où Ω cRN, N ≥ 2, f une fonction de Carathéodory et M une fonction strictement positive et continue sur R+. Dans le cas où la fonction f est asymptotiquement linéaire à l’infini par rapport à l'inconnue u, on montre, en combinant une technique de troncature et la méthode variationnelle, que le problème admet au moins une solution positive quand la fonction M est non décroissante. Et si f(x, u) = |u|p-1 u + λg(x), où p >0, λ un paramètre réel et g une fonction de classe C1 et changeant de signe sur Ω, alors sous certaines hypothèses sur M, il existe deux réels positifs λ. et λ. tels que le problème admet des solutions positives si 0 < λ <λ. et n'admet pas de solutions positives si λ > λ.. Dans la deuxième partie, on étudie deux problèmes soulevés en dynamique des fluides. Le premier est une généralisation d'un modèle décrivant la propagation unidirectionnelle dispersive des ondes longues dans un milieu à deux fluides. En écrivant le problème sous la forme d'une équation de point fixe, on montre l'existence d'au moins une solution positive. On montre ensuite sa symétrie et son unicité. Le deuxième problème consiste à prouver l'existence de la vitesse, la pression et la température d'un fluide non newtonien, incompressible et non isotherme, occupant un domaine borné, en prenant en compte un terme de convection. L’originalité dans ce travail est que la viscosité du fluide ne dépend pas seulement de la vitesse mais aussi de la température et du module du tenseur des taux de déformations. En se basant sur la notion des opérateurs pseudo-monotones, le théorème de De Rham et celui de point fixe de Schauder, l'existence du triplet, (vitesse, pression, température) est démontré / This thesis consists of two independent parts. The first is devoted to the study of some elliptic problems of Kirchhoff-type in the following form : -M(ʃΩNul² dx) Δu = f(x, u) xЄΩ ; u(x) = o xЄƋΩ where Ω cRN, N ≥ 2, f is a Caratheodory function and M is a strictly positive and continuous function on R+. In the case where the function f is asymptotically linear at infinity with respect to the unknown u, we show, by combining a truncation technique and the variational method, that the problem admits a positive solution when the function M is nondecreasing. And if f(x, u) = |u|p-1 u + λg(x) where p> 0, λ a real parameter and g is a function of class C1 and changes the sign in Ω, then under some assumptions on M, there exist two positive real λ. and λ. such that the problem admits positive solutions if 0 < λ <λ., and no positive solutions if λ > λ.. In the second part, we study two problems arising in fluid dynamics. The first is a generalization of a model describing the unidirectional propagation of long waves in dispersive medium with two fluids. By writing the problem as a fixed point equation, we prove the existence of at least one positive solution. We then show its symmetry and uniqueness. The second problem is to prove the existence of the velocity, pressure and temperature of a non-Newtonian, incompressible and isothermal fluid, occupying a bounded domain, taking into account a convection term. The originality in this work is that the fluid viscosity depends not only on the velocity but also on the temperature and the modulus of deformation rate tensor. Based on the notion of pseudo-monotone operators, the De Rham theorem and the Schauder fixed point theorem, the existence of the triplet, (velocity, pressure, temperature) is shown

The author presents some generalized Burgers' equations for incompressible and isothermal flow of viscous non-Newtonian fluids based on the Cross model, the Carreau model, and the Power-Law model and some simple assumptions on the flows. The author numerically solves the traveling wave equations for the Cross model, the Carreau model, the Power-Law model by using industrial data. The author proves existence and uniqueness of solutions to the traveling wave equations of each of the three models. The author also provides numerical estimates of the shock thickness as well as maximum strain $\varepsilon_{11}$ for each of the fluids.

Die vorliegende Arbeit beschreibt eine experimentelle Untersuchung der Schichtbildung von nichtnewtonschen Flüssigkeiten im Tiefdruckverfahren auf nicht saugfähigen Substraten. Das fluiddynamisch bedingte „viscous fingering“ beim Farbspaltungsprozess soll mittels Hydrophobieren der Druckform gehemmt werden. Ziel ist es, möglichst homogene sowie wellenfreie Schichten zu erzeugen. Um ein direkt miteinander vergleichbares Druckergebnis zu erhalten, wird der Druckstoff parallel mit einer unbehandelten und hydrophobierten Form bedruckt. Als Druckstoff werden anstelle von Druckfarbe funktionale Materialien (vorzugsweise PEDOT:PSS) verwendet und variiert, wobei die elektrischen und geometrischen Schichteigenschaften, beispielsweise der elektrische Widerstand und die Rauheit, zur Ermittlung der gesetzten Ziele untersucht wurden. Hiermit und mittels Nutzung einer hydrophobierten Druckform kann eine deutliche Minderung der Wellenbildung (viscous fingering) bei vielen Druckstoffarten beobachtet werden. Die Minderung des viscous fingering im Farbspaltungsprozess und eine nahezu vollständige Leerung der hydrophobierten Tiefdruckform haben einen wesentlichen Nutzwert für den künftigen Einsatz nicht nur für die „gedruckte Elektronik“. / In this work is described experimental research about layer forming from non-Newtonian fluids in gravure printing on non-porous substrates. The viscous fingering, caused through fluid dynamics at splitting of printed material should be decreased by hydrophobic-surface modification of gravure printing form. The aim was to print wave-free homogenous layers. To achieve comparable results, modified and pure form were used simultaneously to print the same material. The printed material was mainly PEDOT:PSS and other, which is used in printed electronics. The properties (surface tension, viscosity) of printed materials were varied by additives. Printing conditions were varied too. The characteristic of printed layers were studied: resistivity, roughness, density, etc. The results shows decreasing of waviness, roughness and viscous fingering in final layer through use of hydrophobic gravure printing form, compared to print results with common printing form. This can be applied not only in the field of printed electronics.

[pt] O presente trabalho apresenta um modelo simplificado baseado na teoria de lubrificação e conceito de viscosidade equivalente para estudar o escoamento dos fluidos de perfuração através de um espaço anular excêntrico com rotação da coluna de perfuração. Em razão do comportamento não Newtoniano destes fluidos, é significativa a influência da rotação da coluna na perda de carga, já que a viscosidade dos fluidos varia com a taxa de deformação imposta. Uma completa análise deste problema requer uma solução tridimensional da equação de conservação de quantidade de movimento de fluidos não Newtonianos, cuja solução é extremamente cara computacionalmente, tornando necessário o desenvolvimento de modelos simplificados que descrevam bem os principais fenômenos observados nestes escoamentos. No modelo proposto, as equações que governam o escoamento são simplificadas pela teoria da lubrificação em coordenadas cilíndricas. Desta forma, o escoamento tridimensional é descrito por uma única equação diferencial bidimensional para o campo de pressão. O comportamento mecânico do fluido é avaliado através do método da Viscosidade Newtoniana Equivalente. O modelo proposto foi validado a partir de soluções numéricas para o modelo completo, disponíveis na literatura, avaliando a precisão das simplificações adotadas. Os resultados mostram o intervalo no qual o modelo simplificado é preciso e apresentam o efeito das propriedades dos fluidos, da geometria do poço e da rotação da coluna na perda de carga do escoamento.O presente trabalho apresenta um modelo simplificado baseado na teoria de lubrificação e conceito de viscosidade equivalente para estudar o escoamento dos fluidos de perfuração através de um espaço anular excêntrico com rotação da coluna de perfuração. Em razão do comportamento não Newtoniano destes fluidos, é significativa a influência da rotação da coluna na perda de carga, já que a viscosidade dos fluidos varia com a taxa de deformação imposta. Uma completa análise deste problema requer uma solução tridimensional da equação de conservação de quantidade de movimento de fluidos não Newtonianos, cuja solução é extremamente cara computacionalmente, tornando necessário o desenvolvimento de modelos simplificados que descrevam bem os principais fenômenos observados nestes escoamentos. No modelo proposto, as equações que governam o escoamento são simplificadas pela teoria da lubrificação em coordenadas cilíndricas. Desta forma, o escoamento tridimensional é descrito por uma única equação diferencial bidimensional para o campo de pressão. O comportamento mecânico do fluido é avaliado através do método da Viscosidade Newtoniana Equivalente. O modelo proposto foi validado a partir de soluções numéricas para o modelo completo, disponíveis na literatura, avaliando a precisão das simplificações adotadas. Os resultados mostram o intervalo no qual o modelo simplificado é preciso e apresentam o efeito das propriedades dos fluidos, da geometria do poço e da rotação da coluna na perda de carga do escoamento.O presente trabalho apresenta um modelo simplificado baseado na teoria de lubrificação e conceito de viscosidade equivalente para estudar o escoamento dos fluidos de perfuração através de um espaço anular excêntrico com rotação da coluna de perfuração. Em razão do comportamento não Newtoniano destes fluidos, é significativa a influência da rotação da coluna na perda de carga, já que a viscosidade dos fluidos varia com a taxa de deformação imposta. Uma completa análise deste problema requer uma solução tridimensional da equação de conservação de quantidade de movimento de fluidos não Newtonianos, cuja solução é extremamente cara computacionalmente, tornando necessário o desenvolvimento de modelos simplificados que descrevam bem os principais fenômenos observados nestes escoamentos. No modelo proposto, as equações que governam o escoamento são simplificadas pela teoria da lubrificação em coordenadas cilíndricas. Desta forma, o escoamento tridimensional é descrito por uma única equação diferencial bidimensional para o campo de pressão. O comportamento mecânico do fluido é avaliado através do método da Viscosidade Newtoniana Equivalente. O modelo proposto foi validado a partir de soluções numéricas para o modelo completo, disponíveis na literatura, avaliando a precisão das simplificações adotadas. Os resultados mostram o intervalo no qual o modelo simplificado é preciso e apresentam o efeito das propriedades dos fluidos, da geometria do poço e da rotação da coluna na perda de carga do escoamento. / [en] In this work a simplified model based on the Lubrication Theory and on the equivalent viscosity is presented to study the drilling fluid flow dynamics through an eccentric annular space with rotation of the drilling column. As a result of the non-Newtonian behavior of the fluid, the rotation of the column has a significant impact on the pressure drop since the viscosity of the fluids is a function of the strain rate applied. A complete analysis of this problem requires a tridimensional solution of the equation of momentum conservation, which makes the solution computationally expensive. For such, simplified accurate models that describe the main phenomena observed have to be developed. On the proposed model, the flow dynamics equations are simplified by the Lubrication Theory in cylindrical coordinates. Therefore, the tridimensional flow is described by one single two dimensional equation for the pressure field. The mechanical behaviour of the fluid is evaluated through the Equivalent Newtonian Viscosity method. The proposed model, as well as the accuracy of the simplifications used, was validated by comparison with numerical solutions of the complete set of equations available on the literature. The results show the range in which the proposed simplified model is accurate and presents the effects of the fluids properties, well geometry and column rotation on the pressure loss.

Cette thèse présente une étude expérimentale du problème de mélange dans les fluides complexes, étude menée en partenariat avec l’entreprise Sanofi Pasteur. Le mélange est un acte des plus fréquents dans la vie courante et aussi dans l’activité industrielle. On trouve dans la littérature de nombreuses études s’intéressant aux cuves de mélange pour en améliorer les performances à partir d’observations faites à grande échelle. Par contre, à notre connaissance, il y a peu de recherche sur l’hydrodynamique du mélange dans les fluides complexes. Dans notre travail, on étudie des fluides non-Newtoniens formés de solutions diluées de polymères caractérisés par leurs propriétés rhéofluidifiante et viscoélastique. Il s’agit de solutions aqueuses de Polyacrylamide (PAA) ou de la gomme de Xanthan (XG). Afin d’identifier la différence de comportement avec les fluides Newtoniens, une étude expérimentale avec de l’eau est effectuée dans les mêmes conditions que celles pour les fluides non-Newtoniens. Cette étude a été menée, en premier, sur un modèle réduit d’une cuve de mélange de Sanofi Pasteur. Les résultats obtenus, non représentés dans ce mémoire de thèse, nous ont amenés à mettre en place une étude fondamentale de l’écoulement dans un mélangeur de géométrie plus simple. Il s’agit alors de pouvoir contrôler les conditions initiales et de s’affranchir des effets secondaires de l’agitation pour ne s’intéresser qu’au mélange. Pour cela, la géométrie retenue est celle d’un mélangeur en T avec deux entrées perpendiculaires. L’exploration en 2D des champs de vitesse et de concentration de scalaire dans cette jonction en T est assurée simultanément aux moyens des techniques optiques (PIV et PLIF). Les observations montrent un effet non négligeable sur l’hydrodynamique et le mélange lié à la présence de polymères dans l’écoulement. De plus, les résultats obtenus permettent de calculer la tension de Reynolds uv et les flux de masse vc et uc. Ils seront utilisés par la suite pour vérifier leur conformité avec le modèle k epsilon couramment utilisé dans l’industrie. / This thesis presents an experimental study of the mixing in complex fluids which is conducted in partnership with Sanofi Pasteur. The mixture is one of the most common act in everyday life and also in industrial activities. We found in the literature many studies focusing on the mixing tanks with objective to improve performance based on observation of large scale. By cons, in our knowledge, there is few or no research on the hydrodynamics of a mixture in complex fluides. In our work, we study non-Newtonian fluids formed of diluted solution of polymer which characterized by their viscoelastic and shear thinning properties. We used in this study aqueous solutions of polyacrylamide (PAA) or xanthan gum (XG). To identify the difference in behavior with Newtonian fluid, an experimental study with water is carried out under the same conditions as those non-Newtonian fluids. At first, this study was on a reduced mixing tank of Sanofi Pasteur. The results, which not shown in this thesis, led us to develop a fundamental study of flow in a mixer with a simple geometry. The objective is to be able to control the initial conditions and to avoid the side effects of agitation to focus on the mixture. For this, we chose a mixer in a T shape with two perpendicular inputs. Exploring 2D velocity and scalar concentration fields in this T-junction is provided simultaneously of optical techniques (PIV and PLIF). Observations show a significant effect on the hydodynamic and mixture related to the presence of polymers in the flow. In addition, results are used to calculate the Reynolds stress uv and the scalar flux vc and uc. They will be used to check their compliance with the k epsilon model that commonly used in industry.

A simulação computacional de escoamentos incompressíveis multifásicos tem avançado continuamente e é uma área extremamente importante em Dinâmica de Fluidos Computacional (DFC) por suas várias aplicações na indústria, em medicina e em biologia, apenas para citar alguns exemplos. Apresentamos modelos matemáticos e métodos numéricos tendo em vista simulações computacionais de fluidos bifásicos newtonianos e viscoelásticos (não newtonianos), em seus regimes transiente e estacionário de escoamento. Os ingredientes principais requeridos são o Modelo de Um Fluido e o Método da Fronteira Imersa em malhas adaptativas, usados em conjunto com os métodos da Projeção de Chorin-Temam e de Uzawa. Tais metodologias são obtidas a partir de equações a derivadas parciais simples as quais, naturalmente, são resolvidas em malhas adaptativas empregando métodos multinível-multigrid. Em certas ocasiões, entretanto, para escoamentos modelados pelas equações de Navier-Stokes (e.g. em problemas onde temos altos saltos de massa específica), tem-se problemas de convergência no escopo destes métodos. Além disso, no caso de escoamentos estacionários, resolver as equações de Stokes em sua forma discreta por tais métodos não é uma tarefa fácil. Verificamos que zeros na diagonal do sistema linear resultante impedem que métodos de relaxação usuais sejam empregados. As dificuldades mencionadas acima motivaram-nos a pesquisar por, a propor e a desenvolver alternativas à metodologia multinível-multigrid. No presente trabalho, propomos métodos para obter explicitamente as matrizes que representam os sistemas lineares oriundos da discretização daquelas equações a derivadas parciais simples que são a base dos métodos de Projeção e de Uzawa. Ter em mãos estas representações matriciais é vantajoso pois com elas podemos caracterizar tais sistemas lineares em termos das propriedades de seus raios espectrais, suas definições e simetria. Muito pouco (ou nada) se sabe efetivamente sobre estes sistemas lineares associados a discretizações em malhas compostas bloco-estruturadas. É importante salientarmos que, além disso, ganhamos acesso ao uso de bibliotecas numéricas externas, como o PETSc, com seus pré-condicionadores e métodos numéricos, seriais e paralelos, para resolver sistemas lineares. Infraestrutura para nossos desenvolvimentos foi propiciada pelo código denominado ``AMR2D\'\', um código doméstico para problemas em DFC que vem sendo cuidado ao longo dos anos pelos grupos de pesquisa em DFC do IME-USP e da FEMEC-UFU. Estendemos este código, adicionando módulos para escoamentos viscoelásticos e para escoamentos estacionários modelados pelas equações de Stokes. Além disso, melhoramos de maneira notável as rotinas de cálculo de valores fantasmas. Tais melhorias permitiram a implementação do Método dos Gradientes Bi-Conjugados, baseada em visitas retalho-a-retalho e varreduras da estrutura hierárquica nível-a-nível, essencial à implementação do Método de Uzawa. / Numerical simulation of incompressible multiphase flows has continuously of advanced and is an extremely important area in Computational Fluid Dynamics (CFD) because its several applications in industry, in medicine, and in biology, just to mention a few of them. We present mathematical models and numerical methods having in sight the computational simulation of two-phase Newtonian and viscoelastic fluids (non-Newtonian fluids), in the transient and stationary flow regimes. The main ingredients required are the One-fluid Model and the Immersed Boundary Method on dynamic, adaptive meshes, in concert with Chorin-Temam Projection and the Uzawa methods. These methodologies are built from simple linear partial differential equations which, most naturally, are solved on adaptive grids employing mutilevel-multigrid methods. On certain occasions, however, for transient flows modeled by the Navier-Stokes equations (e.g. in problems where we have high density jumps), one has convergence problems within the scope of these methods. Also, in the case of stationary flows, solving the discrete Stokes equations by those methods represents no straight forward task. It turns out that zeros in the diagonal of the resulting linear systems coming from the discrete equations prevent the usual relaxation methods from being used. Those difficulties, mentioned above, motivated us to search for, to propose, and to develop alternatives to the multilevel-multigrid methodology. In the present work, we propose methods to explicitly obtain the matrices that represent the linear systems arising from the discretization of those simple linear partial differential equations which form the basis of the Projection and Uzawa methods. Possessing these matrix representations is on our advantage to perform a characterization of those linear systems in terms of their spectral, definition, and symmetry properties. Very little is known about those for adaptive mesh discretizations. We highlight also that we gain access to the use of external numerical libraries, such as PETSc, with their preconditioners and numerical methods, both in serial and parallel versions, to solve linear systems. Infrastructure for our developments was offered by the code named ``AMR2D\'\' - an in-house CFD code, nurtured through the years by IME-USP and FEMEC-UFU CFD research groups. We were able to extend that code by adding a viscoelastic and a stationary Stokes solver modules, and improving remarkably the patchwise-based algorithm for computing ghost values. Those improvements proved to be essential to allow for the implementation of a patchwise Bi-Conjugate Gradient Method which ``powers\'\' Uzawa Method.

A simulação computacional de escoamentos incompressíveis multifásicos tem avançado continuamente e é uma área extremamente importante em Dinâmica de Fluidos Computacional (DFC) por suas várias aplicações na indústria, em medicina e em biologia, apenas para citar alguns exemplos. Apresentamos modelos matemáticos e métodos numéricos tendo em vista simulações computacionais de fluidos bifásicos newtonianos e viscoelásticos (não newtonianos), em seus regimes transiente e estacionário de escoamento. Os ingredientes principais requeridos são o Modelo de Um Fluido e o Método da Fronteira Imersa em malhas adaptativas, usados em conjunto com os métodos da Projeção de Chorin-Temam e de Uzawa. Tais metodologias são obtidas a partir de equações a derivadas parciais simples as quais, naturalmente, são resolvidas em malhas adaptativas empregando métodos multinível-multigrid. Em certas ocasiões, entretanto, para escoamentos modelados pelas equações de Navier-Stokes (e.g. em problemas onde temos altos saltos de massa específica), tem-se problemas de convergência no escopo destes métodos. Além disso, no caso de escoamentos estacionários, resolver as equações de Stokes em sua forma discreta por tais métodos não é uma tarefa fácil. Verificamos que zeros na diagonal do sistema linear resultante impedem que métodos de relaxação usuais sejam empregados. As dificuldades mencionadas acima motivaram-nos a pesquisar por, a propor e a desenvolver alternativas à metodologia multinível-multigrid. No presente trabalho, propomos métodos para obter explicitamente as matrizes que representam os sistemas lineares oriundos da discretização daquelas equações a derivadas parciais simples que são a base dos métodos de Projeção e de Uzawa. Ter em mãos estas representações matriciais é vantajoso pois com elas podemos caracterizar tais sistemas lineares em termos das propriedades de seus raios espectrais, suas definições e simetria. Muito pouco (ou nada) se sabe efetivamente sobre estes sistemas lineares associados a discretizações em malhas compostas bloco-estruturadas. É importante salientarmos que, além disso, ganhamos acesso ao uso de bibliotecas numéricas externas, como o PETSc, com seus pré-condicionadores e métodos numéricos, seriais e paralelos, para resolver sistemas lineares. Infraestrutura para nossos desenvolvimentos foi propiciada pelo código denominado ``AMR2D\'\', um código doméstico para problemas em DFC que vem sendo cuidado ao longo dos anos pelos grupos de pesquisa em DFC do IME-USP e da FEMEC-UFU. Estendemos este código, adicionando módulos para escoamentos viscoelásticos e para escoamentos estacionários modelados pelas equações de Stokes. Além disso, melhoramos de maneira notável as rotinas de cálculo de valores fantasmas. Tais melhorias permitiram a implementação do Método dos Gradientes Bi-Conjugados, baseada em visitas retalho-a-retalho e varreduras da estrutura hierárquica nível-a-nível, essencial à implementação do Método de Uzawa. / Numerical simulation of incompressible multiphase flows has continuously of advanced and is an extremely important area in Computational Fluid Dynamics (CFD) because its several applications in industry, in medicine, and in biology, just to mention a few of them. We present mathematical models and numerical methods having in sight the computational simulation of two-phase Newtonian and viscoelastic fluids (non-Newtonian fluids), in the transient and stationary flow regimes. The main ingredients required are the One-fluid Model and the Immersed Boundary Method on dynamic, adaptive meshes, in concert with Chorin-Temam Projection and the Uzawa methods. These methodologies are built from simple linear partial differential equations which, most naturally, are solved on adaptive grids employing mutilevel-multigrid methods. On certain occasions, however, for transient flows modeled by the Navier-Stokes equations (e.g. in problems where we have high density jumps), one has convergence problems within the scope of these methods. Also, in the case of stationary flows, solving the discrete Stokes equations by those methods represents no straight forward task. It turns out that zeros in the diagonal of the resulting linear systems coming from the discrete equations prevent the usual relaxation methods from being used. Those difficulties, mentioned above, motivated us to search for, to propose, and to develop alternatives to the multilevel-multigrid methodology. In the present work, we propose methods to explicitly obtain the matrices that represent the linear systems arising from the discretization of those simple linear partial differential equations which form the basis of the Projection and Uzawa methods. Possessing these matrix representations is on our advantage to perform a characterization of those linear systems in terms of their spectral, definition, and symmetry properties. Very little is known about those for adaptive mesh discretizations. We highlight also that we gain access to the use of external numerical libraries, such as PETSc, with their preconditioners and numerical methods, both in serial and parallel versions, to solve linear systems. Infrastructure for our developments was offered by the code named ``AMR2D\'\' - an in-house CFD code, nurtured through the years by IME-USP and FEMEC-UFU CFD research groups. We were able to extend that code by adding a viscoelastic and a stationary Stokes solver modules, and improving remarkably the patchwise-based algorithm for computing ghost values. Those improvements proved to be essential to allow for the implementation of a patchwise Bi-Conjugate Gradient Method which ``powers\'\' Uzawa Method.

Atualmente diversos produtos são fabricados por meio de injeção de polímeros, num processo denominado moldagem por injeção: material fundido é injetado em um molde no qual resfria e endurece. Contudo, ao contrário de outros processos de produção, a qualidade da peça criada por meio de moldagem por injeção não depende apenas do material e da sua forma geométrica, mas também da maneira na qual o material é processado durante a moldagem. Por esse motivo, o uso de modelagem matemática e simulações numéricas tem aumentado consideravelmente como maneira de auxiliar o processo de produção e tem-se tornado uma ferramenta indispensável. Desta forma, este projeto tem o propósito de simular o escoamento de fluidos durante a fase de preenchimento do processo de moldagem por injeção, utilizando o modelo 21/2-dimensional, composto por uma equação bidimensional para a pressão, conhecida como equação de Hele-Shaw, e uma equação tridimensional para a temperatura do fluido. Um modelo bidimensional para a temperatura é também desenvolvido e apresentado. Este projeto de doutorado propõe duas estratégias numéricas para a solução da equação de Hele-Shaw. A primeira delas é baseada em uma formulação euleriana do método Smoothed Particle Hydrodynamics, onde os pontos utilizados na discretização não se movem, e não há utilização de malhas. A segunda estratégia é baseada na criação de malhas dinamicamente construídas na região do molde que já encontra-se parcialmente cheio de fluido e subseqüente aplicação do método Control Volume Finite Element Method. Uma estratégia dinâmica do método semi lagrangeano é apresentada e aplicada à solução da equação bidimensional da temperatura. O projeto também pretende investigar três novas abordagens para o tratamento da superfície livre. Duas delas são baseadas na técnica Volume of Fluid e uma delas é uma adaptação meshless do método Front-Tracking. O comportamento não newtoniano do fluido é caracterizado por uma família de modelos de viscosidade. Testes numéricos indicando a confiabilidade das metodologias propostas são conduzidos / Currently, several plastic products are manufactured by polymer injection, in a process named injection molding: molten material is injected into a thin mold where it cools and solidifies. However, unlike other manufacturing processes, the quality of injection-molded parts depends not only on the material and shape of the part, but also on how the material is processed throughout the molding. For this reason, the use of mathematical modelling and numerical simulations has been increasing in order to assist in the manufacturing process, and it has become an essential tool. Therefore, this Sc.D. project has the purpose of simulating the fluid flow during the filling stage of the injection molding process, using the 21/2-dimensional model, compounded by a two-dimensional equation for the pressure field (also known as Hele-Shaw equation) and a three-dimensional equation for the temperature of the fluid. A simpler two-dimensional model for the temperature field is also derived and presented. This project proposes two novel numerical strategies for the solution of Hele-Shaw equation. The first one is based on an Eulerian formulation of the Smoothed Particle Hydrodynamics method, where the particles used in the discretization do not move along as the simulation evolves, thereby avoing the use of meshes. In the second strategy, local active dual patches are constructed on-the-fly for each active point to form a dynamic virtual mesh of active elements that evolves with the moving interface, then the Control Volume Finite Element Method is applied for the pressure field approximation. A dynamic approach of the semi-Lagrangian scheme is applied to the solution of the two-dimensional temperature equation. The project also assesses three new approaches for the treatment of the free surface of the fluid flow. Two of them are based on the Volume of Fluid technique and one of them is a meshless adaptation of the Front-Tracking method. The non-Newtonian behavior is characterized by a family of generalized viscosity models. Supporting numerical tests and performance studies, which assess the accuracy and the reliability of the proposed methodologies, are conducted