The water treatment process with ozone is influenced by various operating parameters and environmental factors that can impact its efficiency. In this study, experiments were conducted using a Venturi pumping frame to investigate the effects of three controllable variables: oxygen flow height, applied voltage level, and water flow rate. The tests aimed to develop a mathematical model that accurately represents the relationship between these input variables and the resulting ozone concentration in the treated water. The experimental data was analyzed using the MODDE 5.0 software, a specialized application for statistical modeling and design of experiments. By fitting the data to appropriate model equations, a mathematical model was obtained that quantifies the influence of each variable and their interactions on the ozone concentration response. To optimize the process performance, a particle swarm optimization (PSO) algorithm was employed to extract the best-fit parameters for the mathematical model. PSO is a computational technique inspired by the social behavior of bird flocks or fish schools, utilizing a population of candidate solutions that evolve iteratively to converge on the global optimum solution. In this case, PSO searched for the model parameter values that minimized the error between predicted and experimentally measured ozone concentrations, rapidly converging to an accurate solution. The optimized mathematical model enables predicting the ozone concentration under any combination of oxygen flow height, voltage, and water flow rate within the experimental range. This predictive capability facilitates identifying the optimum operating conditions to maximize ozone concentration, thereby enhancing the efficiency of the water treatment process. The model serves as a valuable tool for process control, monitoring, and optimization, ensuring consistent treatment quality while minimizing resource consumption and operational costs.