“…This spectral-type approach is based on eigenfunction expansions yielding to solutions where the most features are: (i) an automatic and straightforward global error control and, (ii) an only mild cost increase in overall computational effort for multidimensional situations. Due to its hybrid nature, this scheme has been well indicated for benchmarking purposes and for the validation of different numerical methods in many classes of problems such as, non-linear heat and¯uid¯ow problems, including the Navier±Stokes equations (Pe Ârez Guerrero and Lima et al, 1997;Quaresma and Cotta, 1997), the laminar and turbulent boundary layer equations in duct¯ows (Cotta and Carvalho, 1991;Carvalho et al, 1993;Machado and Cotta, 1995;Figueira da Silva and Cotta, 1996;Cotta and Pimentel, 1998) and, convection-diffusion and eigenvalue problems (Campos Silva et al, 1992;Mikhailov and Cotta, 1994), as well as, laminar boundary layer equations in non-Newtonian¯ows in channels (Magno et al, 1999).…”