The polar interface optical (IO) and surface optical (SO) phonon modes and the corresponding Froehlich electron phonon-interaction Hamiltonian in a freestanding multi-layer wurtzite cylindrical quantum wire (QWR) are derived and studied by employing the transfer matrix method in the dielectric continuum approximation and Loudon's uniaxial crystal model. A numerical calculation of a freestanding wurtzite GaN/AlN QWR is performed. The results reveal that for a relatively large azimuthal quantum number m or wave-number kz in the free z-direction, there exist two branches of IO phonon modes localized at the interface, and only one branch of SO mode localized at the surface in the system. The degenerating behaviours of the IO and SO phonon modes in the wurtzite QWR have also been clearly observed for a small kz or m. The limiting frequency properties of the IO and SO modes for large kz and m have been explained reasonably from the mathematical and physical viewpoints. The calculations of electron-phonon coupling functions show that the high-frequency IO phonon branch and SO mode play a more important role in the electron phonon interaction.