In 2D, the default interpolation of the velocity is quadratic, and this works well for all free surface problems. In the presence of surface tension, a quadratic interpolation of the coordinates (which is linear by default) can help.

In 3D, the default interpolation of the velocity is linear. Using the mini-element interpolation can help get better results on coarse meshes.