Find in Library

Numerical simulations of contaminated spherical drops falling through a stagnant liquid at low Reynolds numbers are carried out using the finite difference method. The numerical results are used to describe the behavior of the surfactant concentrations and to understand the surfactant effects on the fluid motions in detail. The predicted interfacial surfactant concentration, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi mathvariant="sans-serif">Γ</mi></semantics></math></inline-formula>, is almost zero for angles, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi>θ</mi></semantics></math></inline-formula>, below a certain value (the stagnant-cap angle, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>θ</mi><mrow><mi>c</mi><mi>a</mi><mi>p</mi></mrow></msub></semantics></math></inline-formula>), whereas it steeply increases and reaches a large value for <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>θ</mi><mo>></mo><msub><mi>θ</mi><mrow><mi>c</mi><mi>a</mi><mi>p</mi></mrow></msub></mrow></semantics></math></inline-formula> (the stagnant-cap region). The increase in the initial surfactant concentration, <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>C</mi><mn>0</mn></msub></semantics></math></inline-formula>, in the drop enhances the adsorption from the drop to the interface, which results in the increase in <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi mathvariant="sans-serif">Γ</mi></semantics></math></inline-formula> and the decrease in <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>θ</mi><mrow><mi>c</mi><mi>a</mi><mi>p</mi></mrow></msub></semantics></math></inline-formula>. Peaks appear in the predicted Marangoni stresses around <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>θ</mi><mrow><mi>c</mi><mi>a</mi><mi>p</mi></mrow></msub></semantics></math></inline-formula>, which causes similar peaks in the pressure distribution. The high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region and the attenuation of the tangential velocity in the continuous phase. The surfactant flux from the bulk to the interface decreases <i>C</i> in the vicinity of the interface for <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>θ</mi><mo><</mo><msub><mi>θ</mi><mrow><mi>c</mi><mi>a</mi><mi>p</mi></mrow></msub></mrow></semantics></math></inline-formula> and the weak diffusion cannot compensate for the reduction in <i>C</i> by adsorption, which results in <i>C</i> at the interface smaller than <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>C</mi><mn>0</mn></msub></semantics></math></inline-formula>. The pattern of the low <i>C</i> region is determined by the advection and does not smear out because of a small diffusive flux.