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.