Laminar, fully developed flow through single- and double-trapezoidal (or hexagonal) ducts is modeled using a finite-difference method. A coordinate transformation is employed to map the irregular flow cross-section onto a rectangular computational domain. Both H1 and T thermal boundary conditions are considered as they represent the fundamental limiting conditions in most practical applications. Solutions for velocity and temperature variations are obtained for a wide range of duct aspect ratios and with four different trapezoidal angles. The friction factor and Nusselt number results show a strong dependence on duct geometry (aspect ratio gamma and trapezoidal angle theta). The variations off Re, Nu(H1), and Nu(T) With duct aspect ratio for each theta-valued duct are presented in the form of polynomials in gamma. These equations describe the computed numerical values within +/-2% for single-trapezoidal and within +/-1.5% for hexagonal ducts and are of much importance to the design of compact heat exchangers.