The time discretization was performed via the explicit Euler time scheme. A variable time-step method was chosen to ensure numerical stability, and controlled by the time step restrictions based on a Courant-Freidich-Lewis(CFL) number derived from the momentum equation and the surface tension[1]}. We enforced CFL numbers below 0.2 for both momentum and interface driven time-step restrictions.

Дискретизация по времени была выполнена с помощью явной схемы Эйлера. Был выбран метод переменного шага, чтобы обеспечить численную устойчивость, контролируемую ограничениями шага времени на основе числа Куранта-Фрайдиха-Льюиса (CFL), полученного из уравнения движения и поверхностного натяжения[1]}. Мы поддерживали значения CFL ниже 0.2 для ограничений шага времени, связанных с импульсом и интерфейсом.