Numerical simulations have unexplored potential in the study of droplet impact on non-uniform wettability surfaces. In this work, we compare numerical and experimental results to investigate the application potential of a Volume of Fluid method utilized in OpenFOAM®. The approach implements the Kistler model for the dynamic contact angle of impacting droplets. We begin with an investigation on the influence of the most important solver parameters in order to optimize the computational setup and reach the best compromise between computational cost and solution errors, as assessed in comparison to experimental results. Next, we verify the accuracy of the predictions for droplet impact on uniformly hydrophilic or superhydrophobic surfaces. Benchmarking the maximal spreading factor, contact and spreading times, as well as contact-line behavior, we show strong agreement between the present numerical results and the models of Pasandideh-Fard et al. (1996) and Clanét et al. (2004). Lastly, we demonstrate the capability of the model to accurately predict outcome behaviors of droplets striking distributed-wettability surfaces, which introduce 3-D outcome characteristics, even in orthogonal impact. The model successfully predicts droplet splitting and vectoring, as reported in the experiments of Schutzius et al. (2014). Finally, we demonstrate a configuration wherein a spreading droplet becomes arrested within a disc of higher wettability than its surrounding domain. The main contribution of the present work is a numerical model capable of accurately simulating droplet impact on spatially non-uniform wettability patterns of any foreseeable design.