Описан подход к численному моделированию двухфазных флюидопотоков, основанный на методе фазового поля, в котором фазы задаются функцией концентрации. Эта функция является гладкой, так что межфазная граница заменяется достаточно тонким слоем, в котором значения функции концентрации непрерывно изменяются. Такое представление позволяет устойчиво вычислять силы поверхностного натяжения и проводить учет угла смачивания при использовании метода конечных разностей для численного моделирования флюидопотоков и переноса фазы. Работоспособность подхода иллюстрируется на ряде тестовых примеров.
