Added functions for div(phi*vi*v) and div(phi*tau)