This thesis documents the hydrodynamics code Cholla and a numerical study of multiphase galactic winds. Cholla is a massively-parallel, GPU-based code designed for astrophysical simulations that is freely available to the astrophysics community. A static-mesh Eulerian code, Cholla is ideally suited to carrying out massive simulations (> 2048 ³ cells) that require very high resolution. The code incorporates state-of-the-art hydrodynamics algorithms including third-order spatial reconstruction, exact and linearized Riemann solvers, and unsplit integration algorithms that account for transverse fluxes on multidimensional grids. Operator-split radiative cooling and a dual-energy formalism for high mach number flows are also included. An extensive test suite demonstrates Cholla's superior ability to model shocks and discontinuities, while the GPU-native design makes the code extremely computationally efficient - speeds of 5-10 million cell updates per GPU-second are typical on current hardware for 3D simulations with all of the aforementioned physics. The latter half of this work comprises a comprehensive study of the mixing between a hot, supernova-driven wind and cooler clouds representative of those observed in multiphase galactic winds. Both adiabatic and radiatively-cooling clouds are investigated. The analytic theory of cloud-crushing is applied to the problem, and adiabatic turbulent clouds are found to be mixed with the hot wind on similar timescales as the classic spherical case (4-5 t_cc) with an appropriate rescaling of the cloud-crushing time. Radiatively cooling clouds survive considerably longer, and the differences in evolution between turbulent and spherical clouds cannot be reconciled with a simple rescaling. The rapid incorporation of low-density material into the hot wind implies efficient mass-loading of hot phases of galactic winds. At the same time, the extreme compression of high-density cloud material leads to long-lived but slow-moving clumps that are unlikely to escape the galaxy.

Large observational surveys have compiled substantial galaxy samples with an array of different properties across cosmic time. While we have a broad understanding of how galaxies grow and build their observable properties, the details of galaxy growth and evolution pose a fundamental challenge to galaxy evolution theories. Nonetheless, galaxy evolution is ultimately regulated by the properties of the gas reservoir. In this thesis I use numerical simulations to answer key questions related to the galactic gas reservoir, and galaxy mergers: a major transformational process. In Chapter 2 I present an analysis of 28 simulated L* galaxies to understand the physical processes that shape the massive gas reservoir surrounding galaxies (i.e. the circum-galactic medium; CGM). I show that (1) the gas and metal content of the CGM is driven by galaxy growth and the strength/presence of feedback processes, and (2) the ionisation and internal structures of the CGM are shaped by galactic outflows, and active galactic nucleus luminosity. Albeit dependent on internal galactic properties and the physical processes that shape them, the CGM remains greatly diverse, thus posing a challenge for observational surveys. As a follow-up to my study of normal L* galaxy gas halos, in Chapter 3 I present a theoretical study of the effect of galaxy mergers on the CGM. I demonstrate that galaxy mergers can leave a strong imprint on the CGM's gas and metal content, metallicity, and size. The merger can increase (1) the CGM's metallicity by 0.2-0.3 dex within 0.5 Gyr post-merge, and (2) the metal covering fractions by factors of 2-3. In spite of the increase in the CGM's metal content, the hard ionising field during the merger can drive a decline in the covering fractions of commonly observed ions. In Chapter 4 I shift focus to star formation, particularly the effects of galaxy mergers on star formation. While the effects of galaxy mergers have been proven observationally, theoretical predictions are limited to small binary merger suites and cosmological zoom-in studies. I present a statistical study of 27,691 post-merger galaxies from IllustrisTNG to quantify the effect of galaxy mergers on galactic star formation. I report a dependence in the merger-induced star formation rate (SFR) on mass ratio, stellar mass, gas fraction, and galaxy SFR. I also track the evolution of the effects of galaxy mergers demonstrating their decay over ~500 Myr. In Chapter 6, I leverage galactic scaling relations to extend my work on the effects of galaxy mergers to resolved scales. However, before using the simulated resolved scaling relations, I first examine their existence and robustness. In Chapter 5, I demonstrate the emergence of the kpc-scale star forming main sequence (rSFMS) in the FIRE-2 simulations. Nonetheless, the slope of the rSFMS is dependent on the (1) star formation tracer's timescale, and (2) observed resolution, which I propose is caused by the clumpiness of star formation. I develop a toy model that quantitatively captures the effects of clumpy star formation. I then illustrate how the model can be used to characterise the mass of star-forming clumps. Having demonstrated the existence and robustness of known scaling relations in numerical simulations, I explore the effects of galaxy mergers on resolved scales in Chapter 6. I generate synthetic observations for 1,927 post-mergers in IllustrisTNG and examine the radially-dependent merger-driven SFR enhancement, and metallicity suppression in post-mergers. Galaxy mergers preferentially boost star formation in the centres and suppress metallicities globally. The effects of the merger depends on galaxy properties such as stellar mass, SFR, mass ratio, and gas fraction. / Graduate