Numerical methods for multi-scale interfacial flows