Neste trabalho são enfocadas aproximações de elementos finitos para o escoamento de fluidos com tensão de escoamento através de uma expansão planar abrupta. O modelo mecânico combina as equações de balanço de massa e momentum com o modelo viscoplástico de Bingham, regularizado pela equação proposta por Papanastasiou, que descreve através de uma única equação as regiões materiais escoadas e não escoadas. Uma metodologia de Galerkin mínimos-quadrados (GLS) multicampos, em termos de tensão-extra, velocidade e pressão é empregada na aproximação dos escoamentos. Este método é construído para contornar a necessidade de satisfazer as condições de compatibilidade entre os subespaços de elementos finitos tanto de pressão e velocidade quanto de tensão e velocidade. Além disso, devido à escolha apropriada de seus parâmetros de estabilização – resíduos dependentes da malha, resultantes da minimização das equações de Euler-Lagrange do problema, que são adicionados à formulação clássica de Galerkin, aumentando sua estabilidade sem prejudicar sua consistência – a metodologia é capaz de manter-se estável e precisa em escoamentos sujeitos a altos números de Bingham e de Reynolds. As simulações numéricas realizadas consideram escoamentos do fluido de Bingham através de uma expansão planar abrupta 1:4. Para escoamentos lentos, investigou-se o efeito da tensão de escoamento na dinâmica de materiais viscoplásticos, variando-se o número de Bingham de 0,2 até 100. Na sequência, foram considerados os efeitos de inércia através da variação do número de Reynolds entre 0 e 50. Os resultados numéricos mostraram-se capazes de caracterizar com precisão a morfologia das superfícies de escoamento para altos valores do número de Bingham em escoamentos com inércia
