Using three-dimensional smoothed particle hydrodynamics simulations, we investigate the formation of multiple stellar populations (MSPs) in globular clusters (GCs) within the context of their parent galaxies. In our scenario, the second generation (2G) of stars originate from both asymptotic giant branch (AGB) polluters and pristine gas accreted from the host galaxy. Previous theoretical and numerical studies have demonstrated that this "AGB with dilution" model has the potential to alleviate several problems faced by the classical AGB scenario. However, the accretion of pristine gas onto the GC has yet to be investigated within the context of the parent galaxy. This paper presents the preliminary results from our original simulation code which models GC formation from giant molecular clouds in a host galaxy, and subsequent gas accretion onto the GC. By simulating the genesis of the 2G over a 370 Myr time frame, we demonstrate that the fraction of 2G stars is inextricably linked to the GC's environment. Our simulations rationalise the wide variety of abundance patterns, kinematics and 2G concentrations by altering the initial conditions of both the GC progenitor and the host galaxy itself. Most notably, we reproduce a positive correlation between the fraction of 2G stars and the initial mass of the cluster. We discuss the physical implications of our scenario and compare our simulations with observations of the Galactic GC 47 Tucanae (47 Tuc). Finally, we present scaling relations which encompass the wider GC population and serve as a reference for future observations.