Computational fluid dynamics (CFD) and population balance models (PBM) were coupled together for the first time to simulate the synthesis of zirconia nanoparticles in a continuous hydrothermal flow synthesis (CHFS) system with a self-designed confined impinging jet mixing (CJM) reactor. The hydrodynamic and thermodynamic behaviors within the CJM reactor strongly influenced the formation of the ZrO2nanoparticles. Crucial parameters, such as velocities, temperatures, mixing conditions, and reaction rates, were analyzed under various supercritical conditions. Temperature and velocity measurements as functions of distance were also investigated. Normal particle size distribution (PSD) patterns were observed in all cases. The mean particle sizes in this study were calculated and compared using PBM aggregation analysis.