“…We have introduced all the components of the mathematical model for the two-phase NavierStokes system involving surfactants and contact point dynamics, summarized as follows: a) Interface tracking using the domain-decomposition method in Section 2, involving the evolution equation (6) on each segment; b) Surfactant dynamics, which are governed by an advection-diffusion equation (10) in the local coordinates of each segment, as discussed in Section 3; c) Navier-Stokes equations (16, 17) using the pressure-correction method of Timmermans et al [63], given in Section 4; d) Contact-point boundary conditions (30) following Ren and E [49], to account for the contact point dynamics; e) Flow-interface coupling using an immersed boundary method [47], which amounts to computing (15), (22) and (23).…”