A multiphase Volume-of-fluid (VOF) model was developed to gain further insights into the reactive flow parameters and electrical capacitance tomography (ECT) measurements on the remediation of hazardous organic pollutants. Low ozone bubble frequencies were obtained for high surface tension fluids, and the liquid viscosity affected the ozone bubbling frequency. The VOF model indicated that the increase of inlet gas velocity enriched the ozone bubble detachment and concomitantly generated larger ozone bubbles, decreasing the detoxification rates. VOF mappings and ECT visualizations of gas-liquid unveiled preferential routes and highlighted the attenuation of the axisymmetric behavior of the ozonation bubble column under high-interaction regimes.