Neste trabalho, foi utilizado a plataforma de simulação GEANT4 (Geometry and Tracking) para se avaliar a possibilidade de tratamento de esMago com a BNCT (BNCT em inglês, Boron Neutron Capture Therapy), sendo analisada a deposição de energia nas regiões que circundam o esMago. Como o trabalho apresentado é uma parte de um conjunto de estudos e pesquisas realizados sobre o assunto "Terapia de Captura de Nêutrons pelo Boro" (BNCT), direcionado ao câncer de esMago foi necessário partir de pontos básicos para o desenvolvimento do projeto. Para eonstrução dos volumes sensíveis foi levado em eonsicleração medidas aproximadas da anatomia da região do tronco, na altura da sa vértebra e também do esMago. Dentro do esMago, foi inserido um cilindro enriquecido com boro, com intuito de simular a regi ao a ser tratada. Após a construção, foi efetuada uma avaliação de qual biblioteea do GEANT4 relacionado a processos hadrônicos utilizar para obter resultados satisfatórios. Nas simulações realizadas, se utilizou um feixe de nêutrons monoenergético de 0,0253 eV. A fonte foi posicionada em quatro locais distintos, ou seja, três posições externas e uma interna ao tronco. Ao avaliar a deposição de energia desse feixe de nêutrons térmicos, percebe-se que o feixe, quando a fonte está posieionada no exterior do corpo, diverge praticamente para todas as regiões ocasionando assim uma deposição de energia em locais indesejados. No entanto, o mesmo efeito não ocorre quando a fonte encontra-se no interior do corpo, ou seja, a energia depositada, fica restrita a região do alvo. A simulação nos mostrou que o tratamento com fontes de nêutrons térmicos é realmente promissora, porém não existem fontes de nêutrons térmicos compactas para realizar o tratamento. No entanto, podemos pensar em conduzir nêutrons térmicos através de dutos com materiais reflexivos para nêutrons. / ln this work, we used the simulation platform Geant4 (Geometry and Tracking) to evaluate the possibility of cancer treatment of the esophagus with BNCT (Boron Neutron Capture Therapy). We analyzed the deposition of energy in the regions surrounding the esophagus. The present work constitutes a part of a series of studies and research on the subject "Neutron capture therapy with boron "(BNCT), directed to cancer of the esophagus. Therfore it was necessary to start from the basics for the development of the project. For the construction of sensitive volumes was taken into account approximate measures of the anatomy from the trunk, at the 8th. vertebra and the esophagus. Within the esophagus, a cylinder was inserted with enriched boron in order to simulate the region to be treated. After construction, a study was made in order to select the Geant4 library for the relvant hadronic processes and to get satisfactory results. ln the simulations performed, we used a beam of monoenergetic neutrons with 0.0253 eV. The source was placed at four separate locations, three positions outside and one inside of the trunk. To evaluate the deposition of energy of the thermal neutron beam, we find that the beam when the source is positioned outside the body, diverges in almost alI regions thereby causing a deposition of energy in unwanted tissues. However, the same effect does not occur when the source is inside the body, i.e. the deposited energy is restricted to the target region. The simulation has shown us that treatment with sources of thermal neutrons is realIy promising, but there are no sources of compact thermal neutrons to perform the treatment. However, one can think of conducting thermal neutrons through ducts with refiective materiaIs for neutrons.

Neste trabalho apresentamos uma solução analítica para equações de transporte multigrupo de nêutrons unidimensional em reatores nucleares com fontes pulsadas, como por exemplo, reatores ADS (Accelerate Driven Systems). Este trabalho foi desenvolvido em geometria cartesiana para casos de espalhamento isotrópico e anisotrópico utilizando técnicas de transformações integrais e analisando o espaço transformado via aproximações por frequências Bn em uma placa homogênea. A solução final do problema consiste na decomposição do fluxo usando as soluções de casos monoenergéticos como base, formando assim uma solução hierárquica para casos multi-grupos através de soluções de ordem menores e de menos complexidade. A formulação do problema de anisotropia consiste em aproximar o núcleo do termo integral da equação do transporte via polinômios de Legendre. Foram obtidas soluções via inversão numérica da transformada de Laplace no tempo para o caso monoenergético e obtida uma solução analítica para o caso multigrupo dependente da solução monoenergetica. / In this work we present an analytical solution for transport equations in one-dimensional neutron in nuclear reactors using pulsed sources, for example, in ADS reactors (Accelerate or Driven Systems). This work was developed in cartesian geometry for cases of isotropic and anisotropic scattering, using integral transformation techniques and analyzing the transformed space via approximated frequencies by Bn in a homogeneous slab. The final solution of the proposed problem consists in a flux decomposition using the solutions of the one-group case as basis, forming this way a hierarchical solution for the multi-group cases through lower order solutions and with less complexity. The problem formulation for anisotropy is made by approximating the kernel of the integral term of the transport equation via Legendre polynomials. Solutions were obtained via numerical inversion of the Laplace transform in time in the one-group case and we obtained an analytical solution for the multi-group case that depends on the one-group solution.

Neste trabalho, o problema adjunto de transporte de partículas neutras unidimensional é pela primeira vez resolvido pelo Método de Ordenadas Discretas Analítico (ADO). A solução é derivada para problemas monoenergéticos em meios heterogêneos, considerando geometria plana infinita, com simetria azimutal, grau arbitrário de anisotropia e condições de contorno gerais. No caso especial de espalhamento isotrópico, é obtida pelo Método ADO uma solução em termos de autofunções explícitas, bem como um problema de autovalores na forma de perturbação de matriz diagonal. A solução ADO é verificada através de comparações com soluções presentes na literatura, mostrando exatidão superior. A solução é também testada através da abordagem do problema fonte-detector, a qual permite estabelecer uma análise comparativa dos resultados para a resposta de detectores posicionados no interior do domínio. Adicionalmente, a formulação é aplicada em um problema inverso de reconstrução de fontes isotrópicas de partículas, baseado na formulação fonte-detector, para problemas de transporte com geometria e materiais conhecidos. Finalmente, como o problema inverso é mal posto, é abordada a técnica de regularização de Tikhonov, uma modificação no processo de cálculo do problema inverso original, que torna o método de inversão menos sensível a ruídos nas medições. / In this study, the one-dimensional adjoint particle transport problem is solved for the rst time using the Analytical Discrete Ordinates Method (ADO). Monoenergetic problems in heterogeneous slabs with azimuthal symmetry, arbitrary degree of anisotropy and general boundary conditions are taken in consideration. In the particular case of isotropic scattering, the ADO provides a solution in terms of explicit eigenfunctions, as well as with an eigenvalue problem in the form of a rank-1 update of a diagonal matrix. The ADO solution for the adjoint equation is veri ed against solutions presented in the literature, showing to be more accurate. The solution is also tested through source-detector problem approach, which allows for a comparative analysis of the results to be established by computing the response of the detectors inside the domain. In addition, the ADO solution to the adjoint equation is applied in an inverse problem of isotropic source reconstruction, based on the source-detector problem, in domains with known geometry and materials. Finally, since the inverse problem is ill posed, the Tikhonov's regularization technique, a modi cation in the process of computing the inverse problem solution is applied in order to deal with noisy data.

COSTA, Marcelo Henrique de Araújo Santos. Fenômenos de Transporte em Meios Porosos e Interfaces Fractais. 2006. 103 f. Tese (Doutorado em Física) - Programa de Pós-Graduação em Física, Departamento de Física, Centro de Ciências, Universidade Federal do Ceará, Fortaleza, 2006. / Submitted by Edvander Pires (edvanderpires@gmail.com) on 2015-05-22T18:08:10Z No. of bitstreams: 1 2006_tese_mhascosta.pdf: 4361476 bytes, checksum: d1ad57b4bdfc1adb7ee1dc8cafc3333d (MD5) / Approved for entry into archive by Edvander Pires(edvanderpires@gmail.com) on 2015-05-22T20:04:44Z (GMT) No. of bitstreams: 1 2006_tese_mhascosta.pdf: 4361476 bytes, checksum: d1ad57b4bdfc1adb7ee1dc8cafc3333d (MD5) / Made available in DSpace on 2015-05-22T20:04:44Z (GMT). No. of bitstreams: 1 2006_tese_mhascosta.pdf: 4361476 bytes, checksum: d1ad57b4bdfc1adb7ee1dc8cafc3333d (MD5) Previous issue date: 2006 / In this work, we investigate different transport phenomena through irregular media by means of numerical simulations. Initially, we study the effect of the critical percolation disorder on pore networks under diffusion-reaction conditions. Our results indicate the existence of three distinct regimes of reactivity, determined by the dimensionless parameter E=D/(Kl^2), where D is the molecular diffusivity of the reagent, K is its chemical reaction coefficient, and l is the length scale of the pore. At low values of E, the flux of the reacting species penetrating the network follows the classical scaling behavior, namely F~LE^(1/2). At intermediate values of E, the influence of the fractal morphology of the percolating cluster results in an anomalous behavior, F~L^(A/2)E^B, with an exponent B=0.34. At high values of E, the flux of the reagent reaches a saturation limit, F_SAT, that scales with the system size as F_SAT=L^A, with an exponent A=1.89, corresponding to the fractal dimension of the sample-spanning cluster. In the second part of this work, we study how the irregularity of the geometry influences the sequential deactivation of an interface accessed by diffusion. By using the notion of active zone, we propose a conjecture which constitutes an extension of Makarov theorem. In the third part, we investigate the steady-state heat transport in a fluid flowing through a two-dimensional channel whose walls are irregular interfaces. Once more, we apply the notion of active zone to investigate the effect of the interface geometry on the heat exchange efficiency of the system for different conductive-convective conditions. Compared with the behavior of a channel with smooth interfaces and under conditions in which the mechanism of heat conduction dominates, the results indicate that the effect of roughness is almost negligible on the efficiency of the heat transport system. On the other hand, when the convection becomes dominant, the role of the interface roughness is to generally increase both the heat flux across the wall as well as the active length of heat exchange, when compared with the smooth channel. Finally, we show that this last behavior is closely related with the presence of recirculation zones in the reentrant regions of the fractal geometry. / Neste trabalho investigamos diversos fenômenos de transporte tendo lugar através de meios irregulares por meio de simulação computacional. Inicialmente, tratamos do efeito da desordem crítica em redes percolantes de poros sujeitas à difusão e reação química. Verificamos a existência de três regimes distintos, determinados pelo parâmetro adimensional E=D/(Kl^2), onde D é a difusão molecular, K o coeficiente de reação química e l um comprimento característico. Para valores baixos de E, o fluxo de reagente que penetra a rede obedece à relação de escala clássica, F~LE^(1/2). Para valores intermediários de E, a influência da morfologia fractal do agregado de percolação resulta em um regime anômalo, F~L^(A/2)E^B, com um expoente B=0.34. Para valores altos de E, o fluxo de reagente atinge um limite de saturação, F_SAT, e escala com o tamanho do sistema na forma F_SAT=L^A, onde A=1.89 corresponde à dimensão fractal do agregado incipiente de percolação. Em uma segunda etapa do trabalho, analisamos o efeito da geometria irregular na desativação seqüencial de uma interface acessada por difusão. Aplicando o conceito de zona ativa, propomos uma conjectura que se constitui numa extensão do teorema de Makarov. Na terceira parte deste trabalho, investigamos o transporte estacionário de calor no escoamento de um fluido através de um tubo bidimensional, cujas paredes são interfaces irregulares. Mais uma vez, utilizando o conceito de zona ativa, investigamos o efeito da geometria da interface na eficiência de troca térmica do sistema em diferentes condições difusivo-convectivas. Em condições nas quais o mecanismo de transporte dominante é a condução, a comparação entre os resultados dos tubos liso e rugosos indica que o efeito da rugosidade é quase desprezível sobre a eficiência de dispositivos de transporte de calor. Por outro lado, quando a convecção torna-se dominante, a rugosidade passa a ter um papel importante e, em geral, o fluxo de calor e o comprimento da zona ativa aumentam com a rugosidade da interface de troca. Finalmente, mostramos que esse último comportamento está relacionado com as zonas de recirculação, presentes nas reentrâncias da geometria fractal.

O principal objetivo desta tese é analisar a influência do parâmetro seção de choque macroscópica total na solução da aproximação SN da equação de transporte e os efeitos do parâmetro de condução-radiação no problema não-linear condutivo-radiativo acoplado. Inicialmente, é estudada a aproximação SN da equação de transporte estacionária, com um grupo de energia, numa placa plana homogênea, sem fonte externa. Esta aproximação corresponde a um sistema de N equações diferenciais ordinárias lineares de primeira ordem, cuja solução é a função exponencial da matriz do sistema, denominada matriz de transporte SN. Através da resolução de problemas e considerando a seção de choque macroscópica total como parâmetro de controle, é verificada numericamente a mudança de comportamento dos autovalores da matriz de transporte SN que determinam a passagem de soluções não-oscilatórias para soluções oscilatórias da aproximação SN· Estes pontos, referidos como pontos de bifurcação, são calculados para problemas com seção de choque diferencial de espalhamento com anisotropia de graus L = 1 e L = 2. Para pequenas ordens de quadratura N, são realizadas simulações para a análise computacional do polinômio característico da matriz de transporte SN e do seu número de condicionamento em termos da seção de choque macroscópica total. Este estudo é estendido para o problema com dois grupos de energia. Os resultados deste trabalho são comparados com resultados encontrados na literatura. Na parte final deste trabalho, é estudada a influência do parâmetro de condução-radiação na solução do problema acoplado de transferência radiativa e calor condutivo numa placa plana homogênea, com espalhamento anisotrópico, considerando condições de contorno de Dirichlet bem como reflexão especular e difusa. Simulações numéricas obtidas com o uso do método LTSN combinado com o método da decomposição são apresentadas. / The main objective of this thesis is to analyze the infiuence of the parameter total macroscopic cross section in the solution of the SN approximation to the transport equation and the effects of the conduction-radiation parameter in a coupled nonlinear conductive-radiative problem. lnitially, it is studied the SN approximation of the steadystate transport equation with one group of energy, in a homogeneous slab, without externai source. This approximation corresponds to a system of N linear ordinary differential equation of first order, whose solution is the exponential function of the matrix system, termed SN transport matrix. Through the resolution of problems and viewing the total macroscopic cross section as the control parameter, it is numerically verified the change of behavior of the eigenvalues of the SN transport matrix that determine the passage from non-oscillatory to oscillatory solutions of SN approximation. These points, refereed as bifurcation points, are calculated for problems with differential scattering cross section with anisotropy of degrees L = 1 and L = 2. For small orders o f quadrature N, simulations are accomplished for the computational analysis of the characteristic polynomial of SN transport matrix and of its conditioning number in terms of the total macroscopic cross section. This study is extended for the problem with two groups of energy. The results of this work are compareci with available results in the literature. In the final part of this work, it is studied the infiuence of the radiation-conduction parameter in the solution of coupled radiative transfer and conductive heat problem in a homogeneous slab, with anisotropic scattering, considering the Dirichlet's boundary conditions as well specular and diffuse refiection. Numerical simulations obtained with the use of the combined LTSN and the decomposition methods are reported.

No trabalho proposto são apresentados dois novos métodos para a obtenção de soluções de equações diferenciais parciais. O primeiro fornece soluções exatas para problemas difusivos transientes e o segundo mapeia as soluções obtidas em novas soluções para equações diferenciais parciais não-lineares. As soluções dos problemas difusivos são expressas como séries geométricas truncadas, enquanto os mapeamentos são obtidos através do emprego de transformações de Bäcklund. As principais características das formulações propostas são o caráter analítico das soluções obtidas e o baixo custo computacional requerido para efetuar as operações envolvidas. Simulações numéricas são apresentadas. / In this work two analytical methods for solving partial differential equations are proposed. The first method furnishes exact solutions for unsteady diffusion problems and the second one performs mappings which converts the solutions obtained into new exact solutions for nonlinear partial differential equations. The solutions for the diffusion problems are written as truncated geometric series and the mappings are obtained by means of Bäcklund transformations. The main features of the proposed formulations are the analytical character of the solutions obtained and the low computational cost demanded to carry out the calculations. Numerical results are reported.

Neste trabalho, usamos a transformada de Laplace para desenvolver expressões para as condições de contorno tipo albedo para uma e duas regiões refletoras. Nós apresentamos a aplicação da condição de contorno tipo albedo de maneira não convencional. Na prática, os meios multiplicativos dos reatores nucleares são normalmente circundados por materiais refletores, usados para reduzir a fuga de nêutrons. No intuito de retirar a região refletora dos cálculos, introduzimos um coeficiente de reflexão ou parâmetro albedo. Usamos este parâmetro para resolver numericamente a equação da difusão monoenergética e as equações da difusão multigrupo com dois grupos de energia, tanto com fonte fixa quanto com fonte de fissão, pelo método de diferenças finitas. Para os casos de fonte fixa, encontramos os fluxos de nêutrons sem albedo e comparamos com albedo para uma região e duas regiões e verificamos a precisão e a redução no tempo computacional. Já para os casos de fonte de fissão (problema de autovalor), encontramos, sem o parâmetro albedo, os fluxos de nêutrons, os fatores de multiplicação efetivos (Keff), e a potência gerada por região. Comparamos os resultados com albedo para uma região e duas regiões e verificamos a precisão e redução no tempo computacional. A extensão para mais regiões torna-se possível seguindo os passos aqui utilizados, ainda que haja, em contrapartida, um esforço algébrico crescente com o aumento de regiões. / In this dissertation we use the Laplace transform to derive expressions for nonstandard albedo boundary conditions for one and two non-multiplying regions at the ends of onedimensional domains. In practice, the fuel regions of reactor cores are surrounded by refletor regions that reduce neutron leakage. In order to exclude the refletor regions from the calculations, we introduce a reflection coefficient or albedo. We use the present albedo boundary conditions to solve numerically slab-geometry monoenergetic and multigroup diffusion equations using the conventional finite difference method. Numerical results are generated for fixed source and eigenvalue diffusion problems in slab geometry.

Neste trabalho, foi utilizado a plataforma de simulação GEANT4 (Geometry and Tracking) para se avaliar a possibilidade de tratamento de esMago com a BNCT (BNCT em inglês, Boron Neutron Capture Therapy), sendo analisada a deposição de energia nas regiões que circundam o esMago. Como o trabalho apresentado é uma parte de um conjunto de estudos e pesquisas realizados sobre o assunto "Terapia de Captura de Nêutrons pelo Boro" (BNCT), direcionado ao câncer de esMago foi necessário partir de pontos básicos para o desenvolvimento do projeto. Para eonstrução dos volumes sensíveis foi levado em eonsicleração medidas aproximadas da anatomia da região do tronco, na altura da sa vértebra e também do esMago. Dentro do esMago, foi inserido um cilindro enriquecido com boro, com intuito de simular a regi ao a ser tratada. Após a construção, foi efetuada uma avaliação de qual biblioteea do GEANT4 relacionado a processos hadrônicos utilizar para obter resultados satisfatórios. Nas simulações realizadas, se utilizou um feixe de nêutrons monoenergético de 0,0253 eV. A fonte foi posicionada em quatro locais distintos, ou seja, três posições externas e uma interna ao tronco. Ao avaliar a deposição de energia desse feixe de nêutrons térmicos, percebe-se que o feixe, quando a fonte está posieionada no exterior do corpo, diverge praticamente para todas as regiões ocasionando assim uma deposição de energia em locais indesejados. No entanto, o mesmo efeito não ocorre quando a fonte encontra-se no interior do corpo, ou seja, a energia depositada, fica restrita a região do alvo. A simulação nos mostrou que o tratamento com fontes de nêutrons térmicos é realmente promissora, porém não existem fontes de nêutrons térmicos compactas para realizar o tratamento. No entanto, podemos pensar em conduzir nêutrons térmicos através de dutos com materiais reflexivos para nêutrons. / ln this work, we used the simulation platform Geant4 (Geometry and Tracking) to evaluate the possibility of cancer treatment of the esophagus with BNCT (Boron Neutron Capture Therapy). We analyzed the deposition of energy in the regions surrounding the esophagus. The present work constitutes a part of a series of studies and research on the subject "Neutron capture therapy with boron "(BNCT), directed to cancer of the esophagus. Therfore it was necessary to start from the basics for the development of the project. For the construction of sensitive volumes was taken into account approximate measures of the anatomy from the trunk, at the 8th. vertebra and the esophagus. Within the esophagus, a cylinder was inserted with enriched boron in order to simulate the region to be treated. After construction, a study was made in order to select the Geant4 library for the relvant hadronic processes and to get satisfactory results. ln the simulations performed, we used a beam of monoenergetic neutrons with 0.0253 eV. The source was placed at four separate locations, three positions outside and one inside of the trunk. To evaluate the deposition of energy of the thermal neutron beam, we find that the beam when the source is positioned outside the body, diverges in almost alI regions thereby causing a deposition of energy in unwanted tissues. However, the same effect does not occur when the source is inside the body, i.e. the deposited energy is restricted to the target region. The simulation has shown us that treatment with sources of thermal neutrons is realIy promising, but there are no sources of compact thermal neutrons to perform the treatment. However, one can think of conducting thermal neutrons through ducts with refiective materiaIs for neutrons.

Nesta tese desenvolveu-se um modelo estocástico para simular o transporte de nêutrons em um meio heterogêneo, considerando espectros de nêutrons contínuos e as propriedades nucleares com a sua dependência contínua de energia. Este modelo foi implementado utilizando o método Monte Carlo para a propagação dos nêutrons nos diferentes meios. Devido `a limitação com respeito ao numero de nêutrons que pode ser simulado em tempo de processamento computacional aceitável introduziu-se o volume de controle variável junto `as condições de contornos (pseudo-)periódicas para contornar este problema. A escolha pelo Monte Carlo físico clássico deve-se ao fato de poder decompor em constituintes mais simples o problema de resolver uma equação de transporte. Os constituintes podem ser tratados separadamente, estes são a propagação e a interação, respeitando as leis de conservação de energia e momento, e as relações de probabilidade que determinam a respectiva interação. Está-se consciente do fato que o problema abordado nesta tese ´e longe de ser comparável com a construção de um reator nuclear, porém nesta discussão o alvo principal era desenvolver o modelo Monte Carlo, implementar o código computacional numa linguagem que permite extensões de forma modular. Este estudo permitiu uma análise detalhada da influência da energia sobre a população de nêutrons e seu impacto sobre o ciclo de vida de nêutrons. Dos resultados obtidos, mesmo para um arranjo geométrico simples, pode-se concluir a necessidade de considerar a dependência de energia, ou seja, um fator de multiplicação efetivo espectral deve ser introduzido para cada grupo de energia separadamente. / This thesis has developed a stochastic model to simulate the neutrons transport in a heterogeneous environment, considering continuous neutron spectra and the nuclear properties with its continuous dependence on energy. This model was implemented using Monte Carlo method for the propagation of neutrons in different environment. Due to restrictions with respect to the number of neutrons that can be simulated in reasonable computational processing time introduced the variable control volume along the (pseudo-) periodic boundary conditions in order to overcome this problem. The choice of class physical Monte Carlo is due to the fact that it can decompose into simpler constituents the problem of solves a transport equation. The components may be treated separately, these are the propagation and interaction while respecting the laws of energy conservation and momentum, and the relationships that determine the probability of their interaction. We are aware of the fact that the problem approached in this thesis is far from being comparable to building a nuclear reactor, but this discussion the main target was to develop the Monte Carlo model, implement the code in a computer language that allows extensions of modular way. This study allowed a detailed analysis of the influence of energy on the neutron population and its impact on the life cycle of neutrons. From the results, even for a simple geometrical arrangement, we can conclude the need to consider the energy dependence, is a spectral effective multiplication factor should be introduced each energy group separately.

Este trabalho objetiva o desenvolvimento e a implementação computacional de aproximações Galerkin mínimos-quadrados (GLS) para escoamentos não Newtonianos inelásticos utilizando os modelos de Líquido Newtoniano Generalizado (GNL) e Fluido Quasi- Newtoniano sensível ao tipo de escoamento. Apesar de escoamentos de fluidos não Newtonianos serem cruciais em problemas de engenharia, eles representam, ainda hoje, um assunto em aberto, devido à dificuldade de se criar um modelo matemático para o comportamento não-linear das funções materiais medidas em laboratório e também à complexidade do problema gerado quando da aproximação dos modelos mecânicos disponíveis por métodos numéricos usuais. As equações clássicas de GNL possuem funções de viscosidade capazes de prever comportamentos pseudoplástico, viscoplástico ou dilatante, representado ajustes de curvas de tensão de cisalhamento versus taxa de cisalhamento em escoamentos viscométricos. Os modelos sensíveis ao tipo de escoamento são capazes de representar o comportamento diferenciado em regiões de cisalhamento e extensão, utilizando um parâmetro cinemático de classificação de escoamentos e curvas de escoamento em regimes viscométrico e extensional. As formulações multi-campos em tensão extra, pressão e velocidade (τ-p-u), ou ainda com a taxa de deformação como variável primal, são alternativas para a aproximação numérica de escoamentos não Newtonianos que possuem diversas vantagens. Elas tornam desnecessário o pós-processamento do campo de velocidade para a obtenção dos campos de tensão, representando uma aproximação mais precisa especialmente para elementos de baixa ordem, desejáveis do ponto de vista computacional. Além disso, permitem a intuitiva extensão para a implementação de modelos viscoelásticos diferenciais. No entanto, a compatibilização dos sub-espaços funcionais das variáveis é uma dificuldade na implementação numérica, além das instabilidades numéricas intrínsecas aos operadores advectivos dos modelos. Como forma de superar essas dificuldades, uma estratégia de estabilização do tipo GLS é utilizada, a fim de permitir a implementação numérica das formulações em estudo, gerando a estabilidade necessária para várias combinações de subespaços funcionais. É apresentado um estudo dos princípios da Mecânica do Contínuo, de teoria constitutiva, e do método de elementos finitos, com ênfase nos métodos estabilizados do tipo GLS para formulações u-p, τ-p-u e D-p-u. São apresentados resultados numéricos validando o código computacional desenvolvido, e também investigando escoamentos de fluidos Newtonianos, viscoplásticos, pseudoplásticos e quasi-Newtonianos em geometrias como cavidades e contrações. O método se mostra estável e os resultados fisicamente realistas. No entanto, o modelo sensível ao tipo de escoamento demonstra sérios problemas de convergência nos casos nos quais a não-linearidade material é acentuada. / This work aims the development and computational implementation of Galerkin leastsquares (GLS) approximations for non Newtonian purely viscous flows employing Generalized Newtonian Liquid (GNL) and Quasi-Newtonian flow type sensitive models. Despite non Newtonian flows are crucial in many problems in engineering, they represent an open subject, due to the difficulty of building a mathematical model for the non-linear behavior of material functions obtained experimentally, and also due to the complexity of the problem generated by the approximation of the available mechanical models via usual numerical methods. In the classical GNL equations there are viscosity functions that are able to predict pseudoplastic, viscoplastic and dilatant behavior, representing the curve fitting of shear stress versus shear rate in viscometric flows. The flow type sensitive models are able to predict the different behavior in regions of shearing and extension, employing a kinematic parameter to classify the flows and the flow curves in viscometric and extensional flows. Multi-field formulations in extra-stress, pressure and velocity (τ-p-u), or even with the strain rate as a primal variable, are alternatives to the numerical approximation of non Newtonian flows that possess many advantages. They render unnecessary the post-processing of the velocity field to obtain the stress fields, representing a more accurate for low-order elements, eligible from the numerical standpoint. Besides, they allow the intuitive extension for the implementation of viscoelastic differential models. However, the compatibilization of functional sub-spaces of the variables is a difficulty in the numerical implementation, besides the numerical instability intrinsic to the advective operators in the models. As a path to come over these difficulties, a stabilization strategy of GLS type is employed, so as to allow the numerical implementation of the formulations, achieving the necessary stability for various combinations of functional sub-spaces. A study of the principles of Continuum Mechanics, constitutive theory and the finite element method is presented, with emphasis in the GLS methods for u-p, τ-p-u and D-p-u formulations. Numerical results are presented, validating the computational code developed, and investigating flows of Newtonian, viscoplastic, pseudoplastic and quasi-Newtonian fluids different geometry, as cavities and contractions. The methods are noticeable stable and physically comprehensive. Although, the flow type sensitive model shows serious convergence problems in the cases when the nonlinearity is too pronounced.