Quasi-1D unsteady bubbly cavitating nozzle flows are considered by employing a homogeneous bubbly liquid flow model, where the non-linear dynamics of cavitating bubbles is described by a modified Rayleigh-Plesset equation. The various damping mechanisms are considered by a single damping coefficient lumping them together in the form of viscous dissipation and by assuming a polytropic law for the expansion and compression of the gas. The complete system of equations, by appropriate uncoupling, are then reduced to two evolution equations, one for the flow speed and the other for the bubble radius when all damping mechanisms are considered by a single damping coefficient. The evolution equations for the bubble radius and flow speed are then perturbed with respect to flow unsteadiness resulting in a coupled system of linear partial differential equations (PDEs) for the radius and flow speed perturbations. This system of coupled linear PDEs is then cast into an eigenvalue problem and the exact solution of the eigenvalue problem is found by normal mode analysis in the inlet region of the nozzle. Results show that the steady-state cavitating nozzle flow solutions are stable only for perturbations with very small wave numbers. The stable regions of the stability diagram for the inlet region of the nozzle are seen to be broadened by the effect of turbulent wall shear stress.