This paper deals with isothermal Euler-Poisson system which is used to model collapse of self-gravitating Newtonian star. Density dependent viscosity term is added on the right-hand side of momentum equation and it has been proved that there exists stable shadow wave solution with unbounded density at the origin. This results is extended to the vanishing pressure case.
Nedeljkov et al. (Sun,) studied this question.