In this paper, we construct a class of spherically symmetric and self-similar analytical solutions to the pressureless Navier–Stokes equations with density-dependent viscosity coefficients satisfying [Formula: see text], [Formula: see text] for all [Formula: see text]. Under the continuous density free boundary conditions imposed on the free surface, we investigate the large-time behavior of the solutions according to various [Formula: see text] and [Formula: see text]. When the time grows up, such solutions exhibit interesting information: Case (i) If the free surface initially moves inward, then the free surface infinitely approaches to the symmetry center and the fluid density blows up at the symmetry center, or the free surface tends to an equilibrium state; Case (ii) If the free surface initially moves outward, then the free surface infinitely expands outward and the fluid density decays and tends to zero almost everywhere away from the symmetry center, or the free surface tends to an equilibrium state. We also study the large-time behavior of the solutions for [Formula: see text] without any boundary conditions.