Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper presents a mathematical model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. This mathematical model was incorporated in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. The numerical model developed with GFSSP has been validated against two sets of experimental data and has been shown to predict both propellant subcooling and helium consumption to within 30% in most cases. The test data used for the model validation were taken in 1-g, but the mass diffusion model was developed to be applied in both 1-g and micro-g environments.