____________________________________________________________________________________
Added Jan 10, 2019.
Solve for
Mathematica ✗
ClearAll[u, t, r, z, phi]; lap = Laplacian[u[r, phi, z, t], {r, phi, z}, "Cylindrical"]; pde = D[u[r, phi, z, t], {t, 2}] == c^2*lap; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, u[r, phi, z, t], {r, phi, z, t}], 60*10]];