This will be easier in the next release, but for now something like the following may work:
NMinimize[{Norm[{u0 - uif[x, y], v0 - vif[x, y]}], SignedRegionDistance[MeshRegion[mesh]][{x, y}] <= 0}, {x, y}]
Thank you for the quick response! Just tried your solution and it seems to work fine.
