“…For square ducts, the secondary motion in fully developed laminar flow of viscoelastic fluids was first investigated numerically by Wheeler and Whissler [32] and Townsend et al [33], both numerically and experimentally. More recently, Gervang and Larsen [34], Debbaut et al [35] and Siline and Leonov [36] showed that the magnitude of the secondary-flow velocities are typically two orders of magnitude lower than that of the streamwise velocity (the secondary motion actually varies with the fourth power of the axial velocity at low shear rates and linearly with the axial velocity at high shear rates according to Refs. [35,36]).…”