Here, we model the effect of non-uniform dynamical mass distributions and their associated gravitational fields on the stationary galactic superwind solution. We do this by considering an analogue injection of mass and energy from stellar winds and SNe. We give analytical formulae that establish when an outflow is possible and also characterize distinct flow regimes and enrichment scenarios. We also constraint the parameter space by giving approximate limits above which gravitation, self-gravitation and radiative cooling can inhibit the stationary flow. Later, we obtain analytical expressions for the free superwind hydrodynamical profiles. We find that the existence or inhibition of the superwind solution highly depends on the steepness and concentration of the dynamical mass and the mass and energy injection rates, not just in their effective total values. We compare our results with observational data and a recent numerical work that models artificial massive galaxies assuming a uniform distribution of their parameters. Finally, we put our results in the context of the mass-metallicity relationship to discuss observational evidence related to the selective loss of metals from the less massive galaxies and the recent discovery of high-metallicity dwarf spheroidal (dSph) and dwarf elliptical (dE) galaxies.

In this second part, we generalize the results of the previous paper. We present an analytic superwind solution considering extended gravitationally-interacting dark-matter and baryonic haloes. The incorporation of the latter is critical, since they can have a substantial effect on the hydrodynamics of superwinds generated by massive galaxies. Although the presence of extended and massive haloes does not change the limit for the closed-box enrichment of galaxies established in the first paper, they can trigger an earlier activation of the open-box enrichment scenario, since their gravitational potentials can contribute to the inhibition of the free superwind. Moreover, the incorporation of the extended haloes will also enhance the physical setting behind the superwind model, as we consider mass distributions with properties that emulate the results of recent simulations of $\Lambda$CDM haloes.