The evolution of droplets during the printing process is modeled using the volume of fluid (VOF) method, which involves solving the Navier–Stokes and continuity equations for incompressible flow with multiple immiscible phases on a finite volume grid. An indicator function tracks the interfaces and calculates surface tension forces. A grid independence study confirmed the convergence and efficacy of the solutions. The computational model agreed well with experimental data, accurately capturing the impact, spreading, and recoiling of droplets on a solid surface. Additionally, the model validated the interaction of droplets with hydrophilic and hydrophobic surfaces for both constant and dynamic contact angles. Key non-dimensional numbers (Re, We, Oh) were considered to study the interplay of forces during droplet impact on a solid surface. The final print quality is influenced by droplet dynamics, governed by body forces (surface tension, gravity), contact angle, dissipative forces due to motion, and material properties. Computational studies provide insights into the overall process performance and final print quality under various process conditions and material properties.