A flat-flux expansion on axially symmetric zones has been used, together with a uniform double PI expansion for the boundary angular fluxes, to construct an interface-current method for the calculation of two-dimensional hexagonal cells. Collision and escape probabilities are computed in a volume-preserving, cylindrical cell model, while one of three separate PSS models is used for the calculation of the transmission probabilities. When only one zone has been cylinderized, then use of the heterogeneous PSS model is equivalent to an exact calculation in the actual two-dimensional cell geometry (without cylinderization). Comparison between the different approximations and a Monte Carlo calculation are presented for a typical undermoderated assembly.