В рамках метода функционала плотности в гидродинамике многофазных многокомпонентных смесей показана возможность моделирования поверхностно-активных веществ на границе фаз типа "жидкость-жидкость". Суть метода основана на непрерывном описании многофазных смесей за счет введения в энтропию или свободную энергию слагаемых, квадратичных по градиентам плотностей компонентов. Приведен вывод определяющих уравнений. Численно решены модельные гидродинамические задачи, демонстрирующие ряд характерных явлений: аккумуляция поверхностно-активных веществ на межфазной границе, соответствующее уменьшение поверхностного натяжения, эффекты Гиббса и Марангони.