O objetivo deste trabalho é simular numericamente o escoamento de um fluido viscoplástico através de um cilindro posicionado sobre duas placas planas paralelas, utilizando um algoritmo numérico baseado no método de elementos finitos. A modelagem mecânica do problema é baseada nas equações de conservação de massa e conservação da quantidade de movimento, acopladas ao modelo constitutivo recentemente proposto por Souza Mendes e Dutra [Souza Mendes e Dutra, 2004], por vezes mencionado apenas como fluido SMD ou simplesmente viscoplástico. A modelagem mecânica foi aproximada pelo método de Galerkin Mínimos-Quadrados (GLS), construído para superar a falha numérica presente no método de Galerkin clássico. De forma prévia avalia-se o escoamento de um fluido Newtoniano em uma cavidade forçada, sob condições de inércia desprezível (Re=1) e em condições inerciais onde o número de Reynolds é excursionado entre valores iguais a 100, 400 e 1000. Em seguida, avaliase sob esta mesma geometria, duas outras situações reológicas distintas, pseudoplasticidade e dilatância, nestes casos considerando-se o modelo constitutivo definido para um fluido Newtoniano generalizado (GNL), tendo também como principal objetivo à avaliação do código numérico utilizado no presente trabalho, frente os resultados obtidos na literatura. Por fim aborda-se a análise do escoamento do fluido viscoplástico sob a geometria proposta. Avalia-se nesta condição, variações dos valores do índice de potência n, o valor do número de salto J e principalmente as condições de inércia do escoamento, através do número de Reynolds característico para a função SMD. Em todas as situações analisadas, apresentam-se os resultados numéricos, através das isobáricas de pressão, isolinhas de velocidade e isoregiões de tensão (caracterizando regiões de maior e menor rigidez do escoamento), sempre com a principal finalidade de fornecer importantes informações sobre o fenômeno físico envolvido no escoamento do fluido viscoplástico proposto, além da caracterização da topologia presente em cada condição, podendo na maioria dos casos avaliados serem corroborados com a literatura disponível. / The objective of this text is to simulate numerically the draining of a viscoplastic fluid through a cylinder located between two plain parallel plates, using a numerical algorithm based in the method of finite elements. The modeling mechanics of the problem is based on the equations of mass conservation and movement conservation, connected to the constitutive model recently considered by Souza Mendes and Dutra [Souza Mendes and Dutra, 2004], for times mentioned only as SMD fluid or simply viscoplastic. The modeling mechanics was approached by the method of Galerkin Least-Square (GLS), constructed to surpass the present numerical imperfection of the classic Galerkin method. Previously the text perform an analysis of the Newtonian fluid draining through a lid-driven cavity, under conditions of worthless inertia (Re=1) and in inertial conditions where the Reynolds number is varied between values of 100, 400 and 1000. After that, over the same geometry, two other distinct rheology situations are evaluated, shear-thinning and shear-thickening, in these cases considering itself the defined constitutive generalized Newtonian fluid (GNL) model, also having as main objective to evaluate the numerical code used in this present work, comparing with the results gotten in literature. Finally it is approached analysis of the draining of the viscoplastic fluid under proposed mentioned geometry. It is evaluated in this condition, variations of the values of the power-law index n, value of the jump number J and the conditions of inertia of the draining, through the characteristic of the Reynolds number specified for the SMD function. In all the analyzed situations the numerical results are showing through the isobaric lines, velocity profiles and the stress regions (characterizing regions of greater and minor rigidity of the draining), always with the main purpose to supply important information on the involved physical phenomena in the draining of the considered viscoplastic fluid, beyond the characterization of the present topology in each condition, being able in the majority of the evaluated cases to be corroborated with available literature.

Apresenta-se, neste trabalho, um algoritmo para simular numericamente escoamentos compressíveis tridimensionais. As equações que governam o movimento do fluido são estabelecidas tanto para fluidos viscosos quanto para fluidos não viscosos. Para modelar numericamente o problema, emprega-se o esquema de Taylor-Galerkin de um passo, isto é, utiliza-se uma série de Taylor na discretização temporal do problema, enquanto que o fluido é discretizado no espaço aplicando o método de Galerkin, utilizando elementos finitos hexaédricos de oito nós. A integração das matrizes de elemento é efetuada analiticamente considerando o elemento não distorcido e usando um único ponto de integração para avaliar a matriz Jacobiana. Inclui-se um algoritmo para movimentação da malha, o qual é implementado no contexto da descrição Arbitrária Lagrangeana Euleriana (A.L.E.). O programa é codificado na linguagem FORTRAN e vetorizado para aproveitar as características dos processadores vetoriais existentes nos atuais supercomputadores. / An algorithm to simulate numerically 3-D compressible flows is proposed in this work. The governing equations are established both for viscous and non viscous fluids. The one step Taylor-Galerkin scheme is used in order to obtain a useful numerical model. Time integration is accomplished using Taylor series while spatial discretization is carried out applying the Galerkin method and employing the eight node hexahedron finite element. Element matrices integration is and analytically calculated considering undistorted elements and using one integration point to evaluate the Jacobian Matrix. A mesh movement algorithm is included in the context of the Arbitrary Lagrangian-Eulerian (A.L.E.) description. The code is written in FORTRAN language, and is vectorized in order to take advantage of fast vectorial processors existing in modern supercomputers.

Neste trabalho é apresentada uma solução analítica para a dispersão vertical turbulenta em uma camada limite estável, para uma fonte área instantânea. Esta solução é obtida pela aplicação da transformada de Laplace na variável temporal da equação de difusão e também foi usado um coeficiente de dispersão derivado das teorias de similaridade local e da difusão estatística. A transformada inversa é ,obtida através do esquema numérico de quadratura Gaussiana. São apresentadas várias simulações para diversas alturas da fonte área e os resultados obtidos apresentaram boa comparação com os resultados disponíveis na literatura. / In this work is presented an analytical solution for vertical poluent dispersion of a instantaneous area source in an atmospheric stable boundary layer. This solution is obtained by the aplication ofLaplace Transform on the time-dependent, one dimensional diffusion equation and using the dispersion coefficient derived from local similarity and statistical diffusion theories. The In verse o f the Transform is achieved through numerical f scheme of Gausian quadrature. Numerical simulations and comparisions are reported for several heights o f the source and the results attained a good concordance with other results.

O objetivo do presente trabalho é determinar numericamente o perfil de umedecimento em um meio poroso submetido a um processo de infiltração. A descrição do movimento transiente de água em solo não-expansivo é dada em termos da equação não-linear de Fokker-Planck ou da equação de Richards. As dificuldades na solução da equação são não somente devido aos termos não-lineares envolvidos mas também às condições de fronteira dependentes do tempo e ao ajuste dos dados experimentais para os parâmetros fisicos. Neste trabalho é considerado um esquema de diferença finita de três niveis para a simulação de infiltração em uma coluna homogênea, semi-infinita satisfazendo a condição de fluxo constante na superficie. O esquema proposto é implicito com erro de truncamento de segunda ordem no tempo e espaço. O termo convectivo é aproximado por quatro diferentes esquemas. A influência da discretização da derivada na fronteira também é analisada através de quatro esquemas. Os resultados obtidos na simulação numérica pelo esquema de três niveis com termo convectivo aproximado por diferença ascendente (upwind) e condição de fronteira discretizada com ponto ficticio apresentam melhor aproximação aos perfis de umedecimento experimentais que os obtidos por esquemas de dois niveis. / The purpose of the present work is to determine numerically the water content profile in a porous medium submitted to an infiltration process. The description of the transient water movement in nonswelling soil is given in terms of the non-linear Fokker-Planck equation or Richards equation. The difficulties in solving this equation are due not only the nonlinearities involved but also to time-dependent boundary conditions and to the fitting of the experimental data for the physical parameters. In this work, we consider a three-level finite difference scheme for the simulation of infiltration in an homogeneous semi-infinite column satisfying constant flux condition at the surface. This is an implicit scheme which is second-order accurate in time and space. The convective term is approximated by four different schemes. The influence of the discretization of the boundary derivative is analysed by four schemes. The results of the numerical simulation by the three-level scheme with convective term approximated by forward difference (upwind) discretizated with ficticious and point boundary exhibited condition a closer approximation to the experimental water content pro1ile than the ones obtained by two-level scheme.

El objetivo de este trabajo es presentar una formulación matemática para el estudio del movimiento en sistemas de aguas "poco profundas" y su implementación en el sistema HYDRO, lenguage orientado para Hidrodinimica computacional, el cual forma parte del Proyecto LORANE, desarrollado por el Curso de Pos-Graduación en Ingenier1a Civil de la U.F.R.G.S.- Porto Alegre (Brasil).

Artérias contém curvas ou bifurcações, de tamanho pequeno a grande, que são comumente afetadas pela aterosclerose, doença que atinge milhares de pessoas. Pesquisas indicam que a predisposição ao aparecimento da lesão é conseqüência do comportamento do fluido, que apresenta características peculiares nestas regiões. No presente trabalho, realiza-se a simulação numérica do fluxo sangüíneo na bifurcação da artéria carótida. Para escoamentos completamente desenvolvidos, com altas taxas de cisalhamento, como ocorre na região da bifurcação da artéria carótida, o sangue assume comportamento de um fluido newtoniana incompressível com massa específica e viscosidade um pouco acima da água. Observa-se a ocorrência de uma interação fluido-parede na artéria; por ser o fluxo de natureza pulsátil, as paredes detentoras de propriedades elásticas respondem à este processo das fases de sístole e diástole do coração. O modelo utilizado considera o escoamento como bidimensional, newtoniana, incompressível, viscoso e com paredes rígidas e fornece informações válidas a respeito do escoamento em questão, detectando as zonas de recirculação e o aumento da tensão de cisalhamento na artéria carótida. O modelo matemático utilizado tem como base as equações de 1\avierStokes, em coordenadas generalizadas com condições de contorno específicas. Adotase o método de solução de diferenças finitas baseado no processo de integração temporal de Runge-Kutta de três estágios com aproximações espaciais e temporais de segunda ordem. Resultados numéricos são apresentados para o escoamento em três geometrias: duto simples, duto curvo e o interior da artéria carótida; eles comparam adequadamente com resultados numéricos e experimentais encontrados na literatura. / Arteries have curves or bifurcations, from small to big size in which can appear atherosclerosis, a disease common in many people. Researches indicate that the lesion surgement is a consequence of fluid behavior, which presents special characteristics in this region. The present work develops the numerical simulation of blood flow trough the carotid artery bifurcation. For fully developed flow, with high strain rates, common in the carotid artery bifurcation region, the blood behaves as an incompressible non-Newtonian fluid with density and viscosity little times greater than that of water. There appear ftuid and arteries wall interaction; as the flow is pulsatile, the walls answer the sistole/diastole of heart cycle. The model employed considers the flow as two-dimensional, Newtonian, incompressible, viscous and for rigid walls and give valuable information to the flow analysis, detecting the recirculation zones presence and shear strain increase in the carotid artery. The mathematical model used follows the Navier-Stokes equations, in generalized coordinates with special boundary conditions. The finite differences based on the second order Runge-Kutta time-stepping scheme of three-stages is adopted. Numerical results are presented for three geometries: simple duct, curved duct and the carotid artery; they compare very well with numerical/analytical results found in the literature.

O presente trabalho simula processos de conformação mecânica em domínios bidimensionais. Para tanto, considera-se o material como sendo um fluido não newtoniana e utiliza-se a formulação da mecânica dos fluidos com função de penalidade. O método empregado para realizar a discretização espacial é o método dos elementos finitos com integração reduzida do termo da função de penalidade. Para discretizar o tempo usa-se diferenças finitas com o esquema de Crank-Nicolson. E apresentado também o método Lagrangiano-Euleriano arbitrário para descrever o movimento da malha de elementos finitos. / Simulations o f metal forming processes in two dimensional domains are presented in this work. The material is considered as a non newtonian flu id and the flow formulation with the penalty function approach is used. Spatial discretization is accomplished using the finite element method with reduced integration for the pressure term. Finite differences is used for time integration with a Crank-Nicolson scheme. The arbitrary Lagrangian-Eulerian scheme is employed to describe the motion of the finite element mesh.

In questo lavoro vengono presentati dei modelli di dispersione analitici, per la qualità dell'aria, basati sulle equazioni di diffusione ed avvezione multidimensionali, in geometria Cartesiana. Questi modelli sono validi per un contaminante passivo e per condizioni di turbolenza omogenea e uniforme e velocità dei vento uniforme, che equivale ad emissioni a quote elevate che si verificano in condizioni di stabilità intermedia e in assenza di di una forte forza di "buoyancy". Le soluzioni sono ottenute con l'uso dei metodo della Trasformata Integralle Generalizzata e sono analitiche nel senso che nessuna approssimazione e fatta nella derivazione. L'idea di base dei metodo si sviluppa nelle fasi seguenti: costruzione di un paio di trasformate attraverso la soluzione di un problema di Sturm-Liouville; applica.zione della trasformata e della soluzione dei sistema di equazioni differenziali ordinarie risultanti; applicazione della trasformata inversa. Sono presentati le simulazioni numeriche e i confronti con risultati disponibili in letteratura. In queste simulazioni i coefficienti di dispersione sono stati considerati dipendenti sia dali a distanza dali a sorgente (caso bidimensionale e tridimensionale) sia dai loro valore asintotico (caso tridimensionale). / Neste trabalho são apresentados modelos de dispersão analíticos de qualidade do ar, baseados nas equações de difusão e advecção estacionárias bidimensional e tridimensional, na geometria cartesiana. Estes modelos são válidos para um contaminante passivo, sujeito a situações de turbulência homogênea e com velocidades do vento médio uniformes, ou seja, para abandonos elevados que ocorrem em condições de estabilidade intermediária e na ausência de fortes empuxos. As soluções são alcançadas com o uso do método da Transformada Integral Generalizada e são analíticas no sentido de que nenhuma aproximação é feita ao longo de sua derivação. A idéia básica do método consiste: elaboração de um par de transformadas pela solução do problema de Sturm-Liouville; aplicação da transformada integral e solução do sistema de equações diferenciais ordinárias resultantes e aplicação da transformada inversa. Simulações numéricas e comparações com resultados disponíveis na literatura são apresentadas. Nestas simulações foram considerados coeficientes de dispersão dependentes da distância da fonte (caso bidimensional e tridimensional) e o valor assintótico do coeficiente de dispersão (caso bidimensional). / In this work, analytical dispersion models for air quality are reportedo The models are based upon the multidimensional steady-state diffusion-advection equation, in the cartesian geometryo They are valid for a passive contaminant and for physical conditions of homogeneous turbulence and uniform mean wind speed, that means, for elevated contaminant releases that occur in intermediate stability conditions without strong buoyancyo The solutions are obtained using the Generalized Integral Transform Technique; they are analytical in the sense that no approximation is made along their derivationo The main idea consists of the steps: construction of a pair of transformation through the solution of a Sturm-Liouville problem, application of the integral transform and solution of result system of ordinary differential equations followed by the application of the inverse transformo Numerical simulations and comparisons with available results in the literature are presentedo In these simulations the dispersion coefficients were assumed to be dependent on the distance o f the source (two and three dimensional cases) as well its asymptotic bebaviour (two dimensional case).

Uma grande parte dos medidores de vazão instalados em todo mundo são placas de orifício. Existem motivos para esta grande aceitação: o seu baixo custo e a boa precisão destes medidores, isto quando se trata de escoamentos monofásicos. Para estes escoamentos uma grande quantidade de dados técnicos estão disponíveis na literatura. Entretanto quando se trata de escoamentos multifásicos a quantidade de informações técnicas disponíveis é muita pequena. Além disto é bem conhecido que as informações obtidas para os escoamentos monofásicos não pedem ser transferidas aos multifásicos sem comprometer a exatidão das medidas. O objetivo deste trabalho é obter informações a respeito do comportamento de escoamentos bifásicos água-ar, ao passarem por uma placa de orifício de bordos chanfrados. Este estudo está baseado na resposta dinâmica das flutuações no escoamento (flutuação do diferencial de pressão) e a posição da veia contraída na placa de orifício. Estas informações, permitirão no futuro estabelecer uma correlação entre estes escoamentos e fornecer subsídios para futuros trabalhos que possibilitem a redução nos erros cometidos ao se usar placas de orifício na medida de escoamentos bifásicos. / Orifice plate is one of the most used flow meter, all over the world. There are good reasons for the great usage: low costs and good uncertainty, which is true for single-phase flow situations, only. Great amount of data are also available in the technical literature. For multi-phase flow situations, thougt, the amount of data is not so. It is well known, from the experience, that technical informatios coming from a singlephase flow situation can not be used straighforward in multi-phase flow, without compromising the uncertainty of the measurement. The purpose of this work is provide more information about the orifice plate in two-phase flow mixture, air-water, to better understand the flow behaviour. The study is basead on the dynamic response of the flow fluctuations (differential pressure fluctuations) and on the position of the vena contracts in the sharp orifice plate. These informations might give good correlation between both types of flow and, in the future, reduce the uncertainty of the measurements in two-phase flow situations.

Este trabalho teve como motivação buscar as causas e apresentar alternativas de soluções que busquem sanar ou reduzir o problema de desgaste localizado, produzido pela erosão das cinzas, nos tubos do economizador das caldeiras de uma usina termoelétrica a carvão. Para este fim foi construído um modelo em escala de parte da caldeira, reproduzindo o duto curvo de retorno dos gases e, na saída do mesmo, um banco de barras de seção circular, simulando o banco de tubos do economizador. Medições de velocidade, utilizando uma sonda tipo tubo de Pitot e anemômetro de fio quente, medições do campo de pressão nas paredes e visualização do escoamento utilizando fios de lã foram aplicados para a análise do escoamento nas condições existentes na usina e nas configurações modificadas, propostas segundo critérios de controle da camada limite. A análise dos dados, comparada com medições de desgaste localizado na caldeira, levaram à identificação do processo hidrodinâmico existente e os efeitos resultantes nas modificações propostas.