Main topic of this thesis are magnetic fluids. These are specific type of fluids which can simplistically be considered as Bingham fluids. Main issues regarding the magnetic fluids mentioned in this thesis are: rheological properties of the magnetic fluids, behaviour of the magnetic fluids and the use of the magnetic fluids in industrial applications. Main goal is to apply Binghams model on the the magnetic fluids assuming that this model can be applied with a good accuracy. Equations which describe behaviour of the Bingham fluids are constructed. Some assumptions which have to be respected to use this analytical equations for magnetic fluids are formulated. Flow of bingham fluid is analytically solved in some simplificated cases with consideration of laminar flow. Analytical results are confronted with numerical ones obtained from CFD software Fluent for the purpose of verification.

No description available.

<p dir="ltr">Due to advances in 3D integration and miniaturization of chips, the power density and number of hotspots within electronic packages have increased rapidly. A major bottleneck in the chip-to-coolant thermal resistance pathway is the interfacial resistance at solid-solid contacts and, therefore, thermal interface materials (TIMs) are employed to minimize interfacial thermal resistance. To improve heat dissipation, thermal grease (a type of TIM) is generally employed to reduce the overall thermal resistance from a heat-generating component to the heat sink. However, these materials degrade throughout their lifetime and the process is not well understood.</p><p dir="ltr">The first part of this dissertation focuses on investigating the degradation behavior of thermal greases using traditional and accelerated reliability techniques. The performance of a thermal grease often worsens with time due to the thermomechanical cycling driven by the coefficient of thermal expansion mismatch between the substrates via pumpout (material moves out of the interface) and dryout (phase separation of the composite material) phenomena. I isolate the effect of thermal cycling (from mechanical cycling) on the degradation of thermal greases by subjecting them to power cycling while holding the bond line thickness constant. In addition to thermocouples in the system, a high-resolution temperature mapping of the thermal grease is leveraged using an infrared microscope for evaluation of local degradation <i>in situ</i>. The results demonstrate a novel pathway for evaluating thermal grease performance by showcasing the importance of the viscosity-temperature hysteresis. However, traditional reliability testing methods such as thermal cycling have long testing periods, often of the order of days or months. Therefore, to accelerate the degradation analysis of thermal greases, I propose adding mechanical cycling while maintaining a constant heat flow rate. The reliability of thermal greases is investigated at different mechanical oscillation amplitudes and squeezing pressures using a novel custom-designed and machined experimental rig. The results uncover that the mechanical reliability of thermal greases depends on the ratio of elastic modulus to viscosity, with higher ratios being more desirable. Meanwhile, the thermal reliability depends upon the synergy of material properties with higher elastic modulus and higher thermal conductivity, resulting in a lesser increase in thermal resistance over the lifetime of thermal greases. </p><p dir="ltr">The second part of this dissertation focuses on the characterization of the rheology of the thermal greases and the associated uncertainty. Thermal greases have complex rheological properties that impact the performance over their lifetime. I perform rheological experiments on thermal greases and observe both stress relaxation and stress buildup regimes, which are not captured by steady shear-thinning models. Instead, a thixo-elasto-visco-plastic and a nonlinear-elasto-visco-plastic constitutive model characterizes each of the observed regimes. I use the models within a data-driven approach based on physics-informed neural networks (PINNs) to solve the inverse problem of determining the rheological model parameters from the dynamic response in experiments. Further, from a microscopic point of view, these rheological behaviors and associated uncertainties arise from the microstructure rearrangements due to particles' inhomogeneous mixing or separation/settling over time. However, this model calibration approach does not address parameter uncertainty arising due to epistemic (limited rheological data) and aleatoric (randomness of rheological experiments) sources. The last part of this dissertation addresses this limitation and quantifies uncertainties arising in the model calibration process. A hierarchical Bayesian inference methodology is used to obtain distributions of the rheological parameters. The uncertainty is further propagated to shear stress distributions and thermal resistances of thermal grease to demonstrate that the rheological models considered are suitable representations of the experimentally observed regimes. </p><p dir="ltr">Therefore, the current dissertation addresses the thermomechanical degradation behaviors and associated complex rheological characteristics of thermal greases. Understanding the degradation and rheology of thermal greases can help design better thermal greases which are degradation-resistant and hence can improve the reliability of electronic packages.</p>

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.