A fast algorithm for particle simulations
Journal of Computational Physics
Numerical solution of partial differential equations
Numerical solution of partial differential equations
Influence of surfactant on rounded and pointed bubbles in two-dimensional stokes flow
SIAM Journal on Applied Mathematics
Journal of Computational Physics
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
Journal of Computational Physics
A level-set method for interfacial flows with surfactant
Journal of Computational Physics
Journal of Computational Physics
A front tracking algorithm for limited mass diffusion
Journal of Computational Physics
A front-tracking method for computation of interfacial flows with soluble surfactants
Journal of Computational Physics
An immersed boundary method for interfacial flows with insoluble surfactant
Journal of Computational Physics
A diffuse-interface method for two-phase flows with soluble surfactants
Journal of Computational Physics
Journal of Computational Physics
A continuum model of interfacial surfactant transport for particle methods
Journal of Computational Physics
An embedded boundary method for soluble surfactants with interface tracking for two-phase flows
Journal of Computational Physics
Hi-index | 31.47 |
We address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the physically representative limit of large bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a transition layer near the interface in which the surfactant concentration varies rapidly, and large gradients at the interface must be resolved accurately to evaluate the exchange of surfactant between the interface and bulk flow. We use the slenderness of the layer to develop a fast and accurate 'hybrid' numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the interfacial free boundary problem. The accuracy and efficiency of the method is assessed by comparison with a more 'traditional' numerical approach that uses finite differences on a curvilinear coordinate system exterior to the bubble, without the separate transition layer reduction. The traditional method implemented here features a novel fast calculation of fluid velocity off the interface.