I will present a new numerical scheme to model surface tension for an interface represented by a level-set function. In contrast with previous schemes, the method conserves fluid momentum and recovers Laplace’s equilibrium exactly. It is formally consistent and does not require the introduction of an arbitrary interface thickness, as is classically done when approximating surface-to-volume operators using Dirac functions. Variable surface tension is naturally taken into account by the scheme and accurate solutions are obtained for thermocapillary Marangoni flows.