ABSTRACTA numerical model, based on the open source code OpenFOAM, was developed to predict jet regimes and total dissolved gas downstream of spillways. The model utilizes the volume of fluid method to track the interface between air and water. A detached eddy simulation model is used for turbulence closure. Transport and dissolution of bubbles are predicted using an Eulerian approach. A bubble number density equation was implemented to predict bubble size changes caused by dissolution and compression. Total dissolved gas was computed using a transport equation that includes the mass transfer between bubbles and water. The model simultaneously captured spillway jet regimes and distribution of total dissolved gas in a spillway sectional model of McNary Dam. Model parameters, gas volume fraction and bubble size at the entrainment region, were calibrated to match total dissolved gas measured in the field under different dam operations.