Abstract: We consider the free boundary value problem to 3D spherically symmetric compressible isentropic Navier–Stokes–Poisson equations for self-gravitating gaseous stars with γ-law pressure density function for 6/5 < γ ≤ 4/3 .For stress-free boundary condition and zero flow density continuously across the free boundary, the global existence of spherically symmetric weak solutions is shown, and the regularity and long time behavior of global solution are investigated for spherically symmetric initial data with the total mass smaller than a critical mass.