In this paper, a quasi-unsteady mathematical model was presented to examine the crystallization fouling of the calcium sulfate aqueous solutions under upward subcooled flow boiling (SFB). An experimental set-up is also fabricated to collect experimental data for validating the model results. The developed model can axially and radially determine the temperature variation. An asymptotic predictive model was also used to calculate the heat transfer coefficients. Theoretical results of the bulk and surface temperatures at any desired location of the test section were coupled with the fouling formation equations to provide an integrated model for prediction of variable fouling resistance and deposit thickness along the heat transfer surface. Besides, variation of local heated surface and deposit surface temperatures as well as the pressure drop were theoretically evaluated over time while the deposit grows on the heated surface. The simulation results showed that reaction is the controlling mechanism when SFB is predominant.