The studies of disordered heterogeneous media and galaxy cosmology share a common goal: analyzing the disordered distribution of particles and/or building blocks at 'microscales' to predict physical properties of the medium at 'macroscales', whether it be a liquid, colloidal suspension, composite material, galaxy cluster, or entire Universe. The theory of disordered heterogeneous media provides an array of theoretical and computational techniques to characterize a wide class of complex material microstructures. In this work, we apply them to describe the disordered distributions of galaxies drawn from realistic simulations. We focus on the determination of lower-order correlation functions, 'void' and 'particle' nearest-neighbor functions, certain cluster statistics, pair-connectedness functions, percolation properties, and a scalar order metric to quantify the degree of order. Compared to analogous homogeneous Poisson and typical disordered systems, the cosmological simulations exhibit enhanced large-scale clustering and longer tails in the void and particle nearest-neighbor functions, due to the presence of quasi-long-range correlations imprinted by early Universe physics, with a minimum particle separation far below the mean nearest-neighbor distance. On large scales, the system appears 'hyperuniform', as a result of primordial density fluctuations, whilst on the smallest scales, the system becomes almost 'antihyperuniform', as evidenced by its number variance. Additionally, via a finite scaling analysis, we compute the percolation threshold of the galaxy catalogs, finding this to be significantly lower than for Poisson realizations (at reduced density ηc = 0.25 in our fiducial analysis compared to ηc = 0.34), with strong dependence on the mean density; this is consistent with the observation that the galaxy distribution contains voids of up to 50% larger radius. However, the two sets of simulations appear to share the same fractal dimension on scales much larger than the average inter-galaxy separation, implying that they lie in the same universality class. We also show that the distribution of galaxies are a highly correlated disordered system (relative to the uncorrelated Poisson distribution), as measured by the τ order metric. Finally, we consider the ability of large-scale clustering statistics to constrain cosmological parameters, such as the Universe's expansion rate, using simulation-based inference. Both the nearest-neighbor distribution and pair-connectedness function (which includes contributions from correlation functions of all order) are found to considerably tighten bounds on the amplitude of quantum-mechanical fluctuations from inflation at a level equivalent to observing twenty-five times more galaxies. This provides a useful alternative to the standard three-particle correlation, and can be computed in a substantially reduced time, though likely requires simulation-based modeling. This work provides the first application of such techniques to cosmology, providing both a novel system to test ...