is that the contact angle no longer is a constraint imposed onto the system, but that the contact angle is a self-selecting variable. Using this approach has a couple of advantages over existing models. Because there is no enforcement of a model that relates the contact angle to contact line speed, this model can reproduce static contact angle hysteresis, without artificially fixing the position of the contact line [@fang2008; @dupont2010]. Additionally, a variable contact angle model can describe flow over a chemically patterned surface [@dupuis2004; @wang2008]. While the Generalized Navier Boundary Condition has been implemented for continuum simulations using a diffuse interface Cahn Hillard method [@qian2006; @ren2007; @wang2008; @cai2015], using a Volume Of Fluid approach has the advantage that many less grid points are needed to resolve the interface [@ding2007]. In contrast to the body force term derived by @mahady2015, the model presented in this work describes wetting of a dry surface. Opposed to @deganello2011, the model presented here is a analytical derivation of the line tension force.

In our work, we uncover the following findings: when looking at the static behavior of droplets, it is found there is a good match between the droplet shape resolved in the simulations 