A two-dimensional theoretical model of a flat-plate dialyzer with ultrafiltration operation was developed in this study. The two-dimensional velocity profile of the fluid and the concentration distribution in the membrane of the dialysis system with ultrafiltration operation were also derived. The Crank-Nicolson method was used to numerically solve the two-dimensional theoretical model. The influences of the retentate phase flow rate, dialysate phase flow rate, ultrafiltration flow rate and channel thickness ratio on the concentration distribution and the mass transfer rate were illustrated. A considerable mass-transfer efficiency improvement was obtained by employing the flat-plate dialyzer with ultrafiltration operation in comparison to the pure dialysis system without ultrafiltration operation. An experiment involving urea separation by the flat-plate dialyzer with ultrafiltration operation was also performed to confirm the accuracy of the proposed mathematical model.