Este trabalho apresenta resultados numéricos para o escoamento incompressível axial laminar em feixes de barras. Para discretização das equações de transporte é empregado o método de volume-de-controle aplicado a um sistema de coordenadas generalizadas. Os resultados reproduzem primeiramente o escoamento em dutos de várias formas através do mapeamento do campo físico (X,Y) no domínio computacional (n). Soluções numéricas e analíticas são comparadas com a finalidade de se qualificar o programa computacional desenvolvido. Em seguida, são mostradas soluções para o escoamento em feixes de barras dispostos em arranjos triangular e quadrado. Também nestes casos os resultados apresentam concordância satisfatória com soluções analíticas existentes.