No presente trabalho são apresentadas simulações bi e tridimensionais de escoamentos em canais com expansão brusca, usando dinâmica de fluidos computacional (CFD). O objetivo é analisar as zonas de recirculação formadas à jusante da expansão para escoamentos laminares. As configurações geométricas são baseadas em trabalhos experimentais com razão de expansão (ER) igual a 1,9423 e razão de aspecto (AR) igual a 36 e com ER=2 e AR=8. Resultados de escoamentos bidimensionais (assumindo hipótese de bidimensionalidade no plano central do canal) e tridimensionais com AR=36 são apresentados para valores de Reynolds abaixo de 2500. Também são apresentados resultados tridimensionais com AR=8 para números de Reynolds abaixo de 400. O conjunto de equações diferenciais parciais (continuidade e quantidade de movimento) foi resolvido usando os métodos numéricos de elementos (Petrov-Galerkin compensado) e volumes finitos (upwind de 1 e 2 ordem) com um algoritmo não permanente, obtendo soluções estáveis permanentes com t?8. Os estudos de malha de volumes finitos mostraram que o elemento quadrilateral uniforme apresenta resultados mais coerentes com os experimentais. Comparando os resultados das simulações bidimensionais e tridimensionais (AR=36) com os experimentais da literatura, os resultados bidimensionais para x1/s apresentam melhor concordância com os experimentais. No entanto, os resultados tridimensionais apresentaram melhor concordância com os resultados experimentais para as outras zonas de recirculação. Portanto, as diferenças encontradas entre os resultados experimentais e numéricos não são devidas aos efeitos tridimensionais não capturados pelo modelo numérico bidimensional. As discordâncias podem estar sendo causadas por fenômenos não relatados nos trabalhos experimentais.