O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Este trabalho tem o objetivo de realizar uma análise de sensibilidade dos parâmetros operacionais e de projeto de uma usina termelétrica. A análise de sensibilidade serve para indicar qual a importância dessas grandezas no resultado do sistema. A planta modelada é a da Usina Termelétrica Presidente Médici – Fase B, localizada em Candiota, RS, onde cada módulo é capaz de produzir até 160 MW de potência. Para esta planta é apresentado o equacionamento de cada componente do sistema. Destacam-se duas modelagens para a turbina a vapor, usadas para a simulação de seu comportamento em carga parcial, o da Elipse de Stodola e o de Schegliáiev. Ambos os modelos são comparados. Também são apresentados os modos de controle de carga por válvula de estrangulamento, por válvula de bocais e por pressão deslizante. A análise de sensibilidade é feita por meio de três métodos distintos, cujos resultados são comparados. O Método Diferencial utiliza as derivadas parciais para estimar a sensibilidade. O Método de Monte Carlo realiza uma série de avaliações do sistema, com os dados de entrada gerados randomicamente em cada avaliação e análise estatística dos resultados para avaliar a sensibilidade das respostas. É desenvolvido neste trabalho o Método da Transformada de Fourier, que gera os dados de forma senoidal e utiliza a transformada de Fourier para medir as sensibilidades. A medida da sensibilidade é feita através de dois índices. O índice de importância é a derivada parcial de um resultado em relação a um parâmetro, adimensionalizado pelas médias das variáveis, e o índice de sensibilidade é a composição fracional da variância de um resultado. Como conclusão, observa-se que os modelos para operação em carga parcial da turbina de Schegliáiev e de Stodola produzem praticamente os mesmos resultados. Os três métodos de análise de sensibilidade conduzem aos mesmos índices de sensibilidade no caso estudado. A ordenação da importância dos parâmetros estudados é idêntica quando se usa tanto o índice de importância quanto o de sensibilidade, mesmo seus valores sejam diferentes para cada parâmetros.

Problemas de transferência de calor e de massa em meios granulares são encontrados em inúmeras situações de interesse tecnológico. Particularmente, o processo de secagem de grãos envolve esses dois fenômenos. Nesta tese é desenvolvido um modelo matemático que descreve os balanços de energia e de massa para o ar e os grãos baseado nas equações de Navier-Stokes. O coeficiente convectivo de transferência de massa entre o grão e o ar é obtido a partir de dados experimentais para a secagem de grãos de soja. As equações governantes do modelo são resolvidas numericamente por um esquema em diferenças finitas. A análise de consistência, estabilidade e convergência também é realizada para o caso unidimensional. Uma solução analítica da equação diferencial parcial com a forma das equações de balanço da energia e da massa para o ar (caso unidimensional) é obtida, fazendo-se comparações entre dados numéricos e a solução analítica, para funções testes no termo fonte, com o objetivo de avaliar o esquema numérico utilizado. Resultados são apresentados fazendo comparações dos valores numéricos calculados através do modelo com dados experimentais encontrados na literatura. Avalia-se também a influência de parâmetros adimensionais envolvidos nos processos de transferência de calor e massa. Estes resultados contribuem para obter um melhor entendimento da transferência de calor e massa em meios granulares, cujas aplicações são encontradas em muitas situações de interesse prático. / Problems of heat and mass transfer in granular media are found in countless situations of technical interest. Particularly, the grain drying process involve those two phenomena. In the present work we develop a mathematical model that describes the energy and the mass balance for the air and the grain based on the Navier-Stokes equations. The convective mass transfer coefficient between the grain and the air is obtained from experimental data for the soy grain drying process. The governing equations of the model are approximated by means of a finite differences scheme. The analysis of consistency, stability and convergence is also made for the unidimensional case. An analytical solution for a partial differential equation with the form similar to the energy and mass balance equations for the air (in the unidimensional case) is obtained; comparisons are realized between the numerical data and the analytical solution for chosen test functions in the source term, with the objective of evaluating the numerical scheme used. Numerical results are presented and compared with the experimental data found in the literature. We evaluate the influence of dimensionless parameters involved in the heat and mass transfer processes as well. These results contribute to obtain a better understanding of the heat and mass transfer in a granular medium, whose applications are found in many situations of technical interest.

Este trabalho define um ambiente que auxilia o usuário a determinar as raízes complexas de um polinômio, dados apenas seus coeficientes complexos. Para tanto é definida uma interface que habilita a comunicação e a apresentação de resultados de forma mais expressiva do que as habitualmente usadas, permitindo que o usuário perceba mais rapidamente as informações que a ferramenta lhe apresenta. Para capacitar o usuário a determinar as soluções o sistema faz uso de uma abordagem baseada em imagens, obtidas a partir da interpretação do polinômio sendo estudado como se fosse uma superfície originada por uma função de dual variáveis. Através do uso de imagens o usuário pode orientar o sistema para determinar as soluções de seu interesse particular, e pode adicionalmente perceber de forma muito clara detalhes do polinômio que não são apresentados por outras ferramentas. São definidas operações que facilitam a manipulação das imagens obtidas, podendo-se efetuar uma serie de computações que eventualmente dispensam os métodos de determinação numérica de raízes, obtendo-se de uma forma alternativa possíveis aproximações para as soluções. Esta abordagem apoiada em imagens é muito útil na determinação de raízes em regiões onde a avaliação do polinômio se mostra instável ou sujeita a erros grosseiros, como por exemplo raízes múltiplas. Quando são geradas imagens pode ser percebida claramente a presença de erros de avaliação, permitindo que o usuário tome a atitude que julgar adequada. Outros sistemas não permitem que tais instabilidades sejam notadas, fornecendo respostas que podem lido ter qualquer nexo com a realidade. Para uso na ferramenta que implementa a solução proposta, existe uma serie de comandos definidos que capacitam o usuário a fazer quaisquer operações de seu interesse. Estes comandos foram já testados no período da implementação do sistema e permitem que sejam feitas operações sobre o polinômio, as imagens geradas e as raízes determinadas. Em todo o projeto do sistema existe a preocupação de tornar claro o que está acontecendo a cada momento com o polinômio ou com os dados dele obtidos, para que o usuário não se canse em demasia, seja interpretando dados do sistema ou procurando decidir o que fazer em seguida. Assim, o use da ferramenta d bastante natural, e tenta atingir a intuição daqueles que a irão usar. Exemplos são apresentados, bem como descrições dos comandos utilizados na interface. De especial interesse são os algoritmos usados na construção das imagens, pois estes permitem uma grande economia nas avaliações do polinômio quando da geração das cenas, reduzindo sobremaneira o tempo de espera do usuário. Quanto as imagens, são fornecidos três tipos de representações: imagens coloridas em duas dimensões, curvas de nível e curvas de sinal. O uso de três dimensões nas imagens foi descartado, e as razões para tanto estão descritas no texto. Todos os tipos de imagens permitem que se determine e identifique com facilidade as soluções desejadas, e são construídas a partir de uma transformação sobre o polinômio original. Esta transformação é explicada em detalhe no transcorrer do texto. / This work describes an environment that helps the user to find the complex solutions of a polynomial, given its coefficients. To do so, a quite expressive user interface is defined providing easy communication between man and machine and also allowing users to perceive much more rapidly the informations given by the tool. The system uses an image-based interface, obtained from a two-variable surface defined by the polynomial. By using these pictures, it is possible to direct the tool, finding any solutions of special meaning to the user. In addition, some details and features of the polynomial can be very easily seen. To manipulate images a set of operations was defined, and this set has so many features that the use of some operations can even avoid the effort of numerically computing the zeros of the polynomial, also reducing the numerical error embedded in these computations. This image-based approach is very useful when finding domains where the evaluation of the polynomial is unstable or prone to errors, e. g., near multiple roots. In pictures, any evaluation problems are easily detected, and the user can decide how to take the answers given by the system. Other tools usually do not provide any means to detect such behavioral oscillations, and can give nonsense answers. There is a set of commands embedded on the tool that implements the proposed solution, and they allow the user to perform any operations of his (her) interest. Such commands were tested at implementation time, and all were useful at some moment. They perform operations on the polynomial, on the pictures and on the detected zeros. All over the design phase there was the care of making clear what is happening to the polynomial or any data, so that the user can always work without being stressed by trying to find out what is happening. So, we try to build a natural approach to this interfacing problem, trying to reach user's intuition. Examples and descriptions of the commands are given, and also there is a description of a quite important algorithm, the one that builds the pictures themselves. This algorithm needs quite few evaluations of the function to generate pictures, thus reducing the waiting times. For the images, there are three presentation types: color images in two dimensions, height curves or signal curves. Three-dimensional pictures were discarded, and the reasons to such decision are described on the text. All kinds of images provide easy and simple identification of solutions. The images are produced from a simple transformation of the starting polynomial. This transformation is also explained in detail.

O objetivo deste trabalho é estudar os efeitos não lineares que resultam na dinâmica de um dímero rígido e amortecido, provido de rotação deslizando sobre um substrato periódico unidimensional. Simulações numéricas são realizadas para dinâmicas com amortecimento nos graus de liberdade translacional e rotacional, para lançamentos do dímero com velocidades translacionais finitas (estado transiente); com o dímero sujeito a uma força externa aplicada na coordenada do centro de massa e com temperaturas finitas (estado estacionário). As equações de movimento descritas em coordenadas de centro de massa e rotacional são caracterizadas por um acoplamento roto-translacional, este que é ativado em regimes distintos durante o deslizamento do dímero, resultando em transferências de energia entre suas coordenadas. A motivação deste estudo é tentar compreender os efeitos rotacionais emergentes na dinâmica do menor objeto possível que seja capaz de rotar, de forma a possibilitar simples contribuições para o entendimento da origem do atrito em escalas nanométricas. No estado de transiente foi possível obter soluções analíticas para os diferentes regimes da velocidade do centro de massa do dímero, bem como uma compreensão extensiva de sua dinâmica, realizada através de um estudo sobre o comportamento geral de sua inclinação angular e com comparações instrutivas com trabalhos anteriores realizados com um dímero vibrante. No estado estacionário foi possível determinar o comportamento no limiar de deslizamento (atrito estático) através do seu estudo como função da inclinação angular inicial. Três diferentes dinâmicas surgem neste estado, uma relacionada com um adátomo, outra com caráter pendular, e a última em que o dímero é capaz de deslizar realizando rotações completas. Nas três dinâmicas surgem efeitos não lineares, sendo que nas duas últimas acontece devido ao acoplamento rototranslacional, em que a natureza das dinâmicas são qualitativamente compreendidas através de uma análise do potencial efetivo como função da inclinação angular. Argumentos e considerações concisas são desenvolvidas com o âmbito de se formular simples expressões aproximadas que predizem com eficácia alguns dos limites que foram previstos. Também foi estudado o efeito da inclusão de temperaturas finitas no sistema, resultando numa atenuação no efeito não linear. / The aim of this work is to study the effects of rotation of the non linear friction of a damped dimer sliding on a 1D periodic substrate. Numerical simulations are performed with: a damping in the translational and rotational coordinate, throwing the dimer with a finite initial translational velocity (transient state); with the dimer subjected to an external force applied in the center of mass coordinate and with finite temperatures (steady state). The equations of motion in terms of center of mass and rotational coordinate show a rototranslational coupling, whose is activated for distinct regimes while the dimer is sliding, resulting in an energy transfer between the coordinates. The motivation of this work is to understand the rotational effects that emerges from the dynamics of the smallest object that can rotate, in order to achieve simple contributions to the understanding of the friction origin in nanometric scale. The first system was the transient state, and we were able to obtain analytical solutions of the center of mass velocity while the dimer is sliding in the three distinct regimes that arise. As well as an extensive comprehension of this dynamics, performed by studying the whole behavior of the angular inclination and instructive comparisons with the vibrating dimer. In the steady state system, an analysis of the behavior of the dimer at the threshold sliding (static friction) as function of the angular inclination allowed a complete understanding of the threshold effects. Three distinct dynamics arise from this state, one related to the adatom, another with a pendulous nature, and the last with the dimer performing complete rotations while is sliding. Non linear effects arise from the thee dynamics, for the two former cases the effects arise from roto-translational coupling, which are qualitatively understood by analysing the effective potential as function of the angular inclinations. Reasonable and proper arguments are developed aiming to construct simple and approximated expressions that predict with good efficiency some of the limit parameters values, that describes the transitions between different dynamics. Also, numerical simulations performed with finite temperature shows an attenuation in the non linear effects.

Neste trabalho estudamos um sistema de equações diferenciais parabólicas que modelam um processo de difusão-reação em duas dimensões da mistura molecular e reação química irreverssível de um só passo entre duas espécies químicas A e B para formar um produto P. Apresentamos resultados analíticos e computacionais relacionados à existência e unicidade da solução, assim como estimativas do erro local e global utilizando elementos finitos. Para os resultados analíticos usamos a teoria de semigrupos e o principio do m´aximo, e a simulação numérica é feita usando diferenças finitas centrais e o esquema simplificado de Ruge-Kutta. As estimativas do erro local para o problema semi-discretizado são estabelecidas usando normas de Sobolev, e para estimar o erro global usamos shadowing finito a posteriori. Os resultados computacionais obtidos mostram que o comportamento da solução está dentro do esperado e concorda com resultados da referências. Assim mesmo as estimativas do erro local e global são obtidas para pequenos intervalos de tempo e assumindo suficiente regularidade sobre a velocidade do fluído no qual realiza-se o processo. Destacamos que a estimativa do erro global usando shadowing finito é obtida sob hipóteses a posteriori sobre o operador do problema e o forte controle da velocidade numa vizinhança suficientemente pequena.

O presente trabalho é dedicado ao estudo de métodos de simulação para ciclos de Rankine. O trabalho é iniciado com a modelagem de um ciclo de Rankine simples e segue evoluindo para configurações mais complexas tal como o ciclo de Rankine com reaquecimento e regeneração. São adotadas as considerações mais convencionais da prática de projeto de centrais termelétricas cujos sistema térmicos baseiam-se no ciclo de Rankine, incluindo-se queda de pressão em tubulações do circuito além de outras perdas. Em seguida, são estabelecidas as expressões matemáticas que possibilitam a determinação das propriedades termodinâmicas da água em seus mais diversos estados ao longo do ciclo. Por último, são desenvolvidos métodos de simulação, chamados neste trabalho de Substituição Sucessiva e Bloco Único, que caracterizam-se pela resolução simultânea do conjunto de equações algébricas dos ciclos elaborados. As simulações são efetuadas através de programas escritos na linguagem Fortran. Os métodos de simulação são aplicados para a obtenção dos resultados considerados mais importantes na análise de sistemas térmicos de potência, tais como rendimento térmico do ciclo, título na saída da turbina, vazões mássicas pelo sistema, potência nas bombas e calor trocado no gerador de vapor e no condensador Na maioria das simulações, estes resultados apresentam-se como funções da: (1) potência elétrica requerida, eficiência isentrópica e pressões na turbina; (2) eficiência térmica, pressão e temperatura no gerador de vapor; (3) pressão e grau de subresfriamento do líquido saturado no condensador e (4) eficiência isentrópica das bombas. São obtidos os mesmos resultados para os métodos de simulação utilizados. O método da Substituição Sucessiva apresentou menor tempo computacional, principalmente para configurações de ciclo mais complexas. Uma aplicação alternativa do método de Bloco Único demonstrou ser inconveniente para ciclos de configurações mais complexas devido ao elevado tempo computacional, quando todas as equações de cálculo das propriedades termodinâmicas são incluídas no sistema de equações a ser resolvido. Melhores rendimentos térmicos e título na saída da turbina foram obtidos para configurações de ciclo de Rankine com reaquecimento e regeneração.

A paralelização de aplicaçõpes envolvendo a solução de problemas definidos sob o escopo da Dinâmica dos Fluidos Computacional normalmente é obtida via paralelismo de dados, onde o domínio da aplicação é dividido entre os diversos processadores, bem como a manutenção do balancecamento durante a execução é um problema complexo e diversas heurísticas têm sido desenvolvidas. Aplicações onde a simulação é dividida em diversas fases sobre partes diferentes do domínio acrescentam uma dificuldade maior ao particionamento, ao se buscar a distirbuição equlibrada das cargas em todas as fases. este trabalho descreve a implementação de mecanismos de particionamento e balanceamento de carga em problemas multi-fase sobre clusters de PCs. Inicialmente é apresentada a aplicação desenvolvida, um modelo de circulação e transporte de susbtâncias sobre corpos hídricos 2D e 3 D, que pode ser utilizado para modelar qualquer corpo hídrico a partir da descrição de sua geometria, batimetria e condições de contorno. Todo o desenvolvimento e testes do modelo foi feito utilizando como caso de estudo o domínio do Lago Guaíba, em Porto Alegre. Após, são descritas as principais heurísticas de particionamento de domínio de aplicações multi-fase em clusters, bem como mecanismos para balanceamento de carga para este tipo de aplicação. Ao final, é apresentada a solução proposta e desenvolvida, bem como os resultados obtidos com a mesma.

Os sistemas de aquecimento através da energia solar são uma fonte interessante de energia acessível e distribuída. No Brasil, a utilização de coletores solares planos já é bem conhecida e difundida. Os coletores solares de tubo evacuado até pouco tempo tinham um custo considerado elevado para o nosso mercado, mas atualmente esta tecnologia está cada vez mais difundida no mercado nacional com custos menores. Para se entender se o custo adicional deste coletor quando comparado aos tradicionais coletores planos é vantajoso, são necessários estudos comparativos para determinar as reais vantagens e quando as mesmas podem proporcionar uma maior economia financeira à instalação. Esta dissertação compara as tecnologias atuais de coletores planos com os de tubo evacuado disponíveis no mercado nacional, mostrando comparações em diferentes regiões do país e analisando se as vantagens encontradas são interessantes do ponto de vista do custo da instalação. Para realizar estas analises foram desenvolvidos testes para caracterização dos tubos destes coletores, de forma a determinar os parâmetros de desempenho dos mesmos. Também foi preparada uma bancada de testes para os coletores solares e foram realizados ensaios de comportamento térmico. Conhecendo os resultados foi possível adequar e utilizar um software de simulação desenvolvido no Laboratório de Energia Solar da Universidade Federal do Rio Grande do Sul. Este software visa a simulação de sistemas de aquecimento de água com diferentes configurações em qualquer região do país. Com a implementação dos coletores de tubo evacuado, e com simulações do seu desempenho em diferentes regiões, é possível comparar o seu comportamento em relação ao coletor plano. Nesta dissertação as curvas de desempenho foram comparadas em diferentes situações, onde foi possível verificar em que momentos o coletor de tubos pode ser melhor que o plano. Também foram realizadas comparações entre a área ocupada e o custo destes coletores para instalações em 3 diferentes regiões do país. Com estas comparações ficou evidente que para climas mais frios o coletor de tubos evacuados pode apresentar melhor desempenho que os coletores planos, mas o seu custo ainda é muito maior que os coletores planos nacionais. / The solar water heating systems are a interesting source of accessible and distributed energy. In Brazil, the flat plate solar collector use is well known. The vaccum tube solar collectors until little time ago had an elevated cost for our market, but, today, the technology is available in our national market with lower costs. And to understand if this costs, that are higher than the flat plane collectors are worth it, some studies are necessary to understand the real benefits of this technology. This work aims to compare both technologies available in the national market, showing comparisons in different regions of the country and studying if the advantages are interesting in the economical point of view. To achieve that some tests were developed in order to obtain the performance parameters of this collectors. Also an external bench was developed to test this collectors according to the Brazilian standard. With the results and knowing the behavior of this collectors it was possible to program the simulation software that was developed in the Laboratório de Energia Solar of the Universidade Federal do Rio Grande do Sul. This software aims at the simulation of water heating systems in different regions of the country. With the implementation of the vacuum collectors in the software, it was possible to simulate different climate conditions over the country and compare the behavior with the flat plate collectors. In this work the performance curves were compared in different climate situations, so that was possible to understand when the vaccum collector can be better. Also some comparisons where made between the occupied area and the costs in three different regions of the country. With the results of this comparisons it is evident that for cold climates the vaccum collector can obtain a better performance than the flat plate collectors, but its cost is still higher.

O principal objetivo dessa dissertação consiste nas simulações numéricas de um escoamento de um fluido viscoplástico usando o método de elementos finitos. A modelagem mecânica empregada está baseada nas equações conservativas de massa e momentum acopladas ao modelo constitutivo recentemente proposto por Souza Mendes e Dutra – aqui chamado de SMD. A modelagem mecânica foi aproximada via método de Galerkin Mínimos-Quadrados (GLS), que é construído para superar a falha numérica do método clássico de Galerkin para formulações mistas sem danificar a sua consistência. Primeiro, a fim de checar o código computacional, um escoamento em uma cavidade forçada. Em seguida, dois escoamentos de fluidos viscoplásticos foram considerados: um fluido SMD escoando através de uma expansão seguida de uma contração abrupta e outro em torno de um cilindro circular. Os resultados numéricos gerados forneceram informações físicas dos fenômenos de viscoplasticidade estudados e se mostraram de acordo com a literatura. / The main objective of this paper consists of numerical simulations of viscoplastic fluid flows using the finite element method. The mechanical modeling employed is based on the mass and momentum conservative equations coupled with the constitutive model recently proposed by Souza Mendes and Dutra - herein only named by SMD. The mechanical modeling was approximated via the Galerkin Least-Squares method (GLS) which is built to overcome the numerical shortcoming of the classical Galerkin method for mixed formulations without damaging its consistency. First, in order to check the computational code, a Newtonian leaky cavity flow. Next, two viscoplastic fluid flows has been considered: a SMD fluid flow through a sudden expansion followed by a contraction geometry and around a circular cylinder. The generated numerical results gave some physical information on the viscoplastic phenomena studied and are in accordance with the literature.