We consider a set of equations modelling contact line dynamics. The model consists of the Navier-Stokes free surface flow with local slip type boundary conditions and gradient surface tension coupled with a mesoscopic local surface model (nonlinear degenerated equations) describes the surface tension variations. The dynamical contact angle and the local surface tension variations are unknowns of the model. We present some mathematical and numerical results.