The zero liquid discharge process is designed to eliminate the discharge of industrial wastewater into the environment, addressing significant environmental concerns that have emerged in recent decades. Falling film evaporator is one of the primary technologies employed in these systems. As observed, there is a noticeable lack of research numerically simulating these phenomena in a vertical pipe including falling film flow. In the present study, an effort has been made to address the existing research gap by conducting numerical simulation of particles deposition and fouling phenomena in a falling film flow in a three-dimensional vertical pipe using computational fluid dynamics. Continuous phases which contain two non-penetrating fluids, i.e., air and water, are simulated using volume of fluid multiphase model and the large eddy simulation turbulence model. Meanwhile, particles motion has been tracked utilizing discrete phase model. Water film enters the pipe with a thickness of 3.7 mm and an inlet velocity of 1 m/s, and particles with a density of 3110 kg/m3 are injected through the water inlet surface in various cases, each with different diameters. As observed, with an increase in particles diameter from 5 to 15 μm, the asymptotic total fouling mass has increased over 14 times, reaching from 0.016 to 0.18 kg/m2, and the total number of deposited particles has also increased. By increasing the fluids inlet velocity from 1 to 1.5 m/s, the total number of deposited particles has decreased throughout the pipe wall, and the total fouling mass has experienced a reduction of 55.8 %.