Most heat exchangers in wastewater treatment systems and the thermal desalination industry experience the deposition of unwanted materials on their surfaces. One of the primary causes is the supersaturation of salts with inverse solubility near the heat exchanger's heated surfaces. In this paper, a long-term simulation of the crystallization fouling of calcium carbonate by developing an efficient solver in the OpenFOAM software is proposed. The k-ω SST model was used to simulate turbulence, and the level-set model was adopted for high-precision and high-quality capturing of the fouling front inside the heat transfer tubes of a forced circulation crystallizer. The results demonstrated that the developed solver is a promising tool for accurately predicting moving boundary conjugate heat transfer. Upon validation of the results, the dynamic evolution of the fouling front and its sophisticated interactions with hydrodynamics were examined. The crystallization fouling was found to be under reaction control. Furthermore, the relatively slight supersaturation caused delayed front growth, which prolonged the fouling time to several months. Besides, findings demonstrated that the intensive competition between removal and deposition mechanisms at higher inlet velocities resulted in decreased asymptotic thermal resistance. Additionally, the velocity increase had a smaller impact on reducing the deposit's asymptotic resistance.