Next Article in Journal
A Direct-Forcing Immersed Boundary Method for Incompressible Flows Based on Physics-Informed Neural Network
Next Article in Special Issue
Rim Breakups of Impacting Drops on a Superhydrophobic Surface and a Superheated Surface
Previous Article in Journal
Numerical Investigation of Orifice Nearfield Flow Development in Oleo-Pneumatic Shock Absorbers
Previous Article in Special Issue
Experimental Correlation for Splashing Condition of Droplets on Solid Substrates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers

1
Graduate School of Engineering, Kobe University, 1-1 Rokkodai Nada, Kobe 657-8501, Japan
2
Department of Mechanical Engineering, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
3
Eindhoven Institute for Renewable Energy Systems (EIRES), Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
4
Graduate School of Societal Safety Sciences, Kansai University, 7-1 Hakubaicho, Takatsuki 569-1116, Japan
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(2), 55; https://doi.org/10.3390/fluids7020055
Submission received: 26 November 2021 / Revised: 18 January 2022 / Accepted: 18 January 2022 / Published: 25 January 2022
(This article belongs to the Special Issue Advances in Multiphase Flow Science and Technology)

Abstract

:
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, Γ , is almost zero for angles, θ , below a certain value (the stagnant-cap angle, θ c a p ), whereas it steeply increases and reaches a large value for θ > θ c a p (the stagnant-cap region). The increase in the initial surfactant concentration, C 0 , in the drop enhances the adsorption from the drop to the interface, which results in the increase in Γ and the decrease in θ c a p . Peaks appear in the predicted Marangoni stresses around θ c a p , 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 C in the vicinity of the interface for θ < θ c a p and the weak diffusion cannot compensate for the reduction in C by adsorption, which results in C at the interface smaller than C 0 . The pattern of the low C region is determined by the advection and does not smear out because of a small diffusive flux.

1. Introduction

Surface-active agents (surfactants) have various effects on the motion of a fluid particle, e.g., the drag coefficient of a spherical bubble under a fully contaminated condition corresponds to that of a solid particle [1,2], capillary waves formed on the surface of a bubble are dampened by the presence of surfactant [3,4], and so on. Boussinesq [5] assumed that the interface contaminated with surfactant behaves similar to a viscous membrane and modeled the increase in drag caused by the presence of surfactant as surface viscosity [1,6]. The drag model for clean spherical fluid particles of low Reynolds numbers, i.e., the Hadamard–Rybczynski model [7,8], can be easily extended to that for contaminated fluid particles by introducing the surface viscosity. On the other hand, Levich [9] succeeded in explaining the surfactant effect on the drag force by accounting for the surface tension gradient, i.e., the Marangoni stress, caused by non-uniform distribution of the surface tension due to the adsorption of surfactant to the interface. Numerical simulations of contaminated fluid particles have been carried out so far by taking into account the Marangoni stress [2,10,11,12,13,14,15,16,17,18,19,20], in contrast, measurements of the Marangoni stress and the interfacial surfactant concentration have rarely been carried out.
Reduction in the surface tension due to surfactant adsorption takes place for various surfactants such as Triton X-100, sodium dodecyl sulfate (SDS), alcohols such as 1-pentanol and so on. In the presence of interfacial flow, the distribution of surfactant at the interface is non-uniform due to advection. For example, the interfacial flow formed at the interface of a drop falling through a stagnant liquid due to gravity transfers surfactant molecules adsorbing to the interface from the nose side to the rear side of the drop, which results in a high surfactant concentration in the rear side. The non-uniform surfactant distribution also results in a non-uniform distribution of surface tension. The distribution of the surface tension can be calculated if the Marangoni stress is measured, and the interfacial surfactant concentration can also be computed from the surface tension by making use of the Langmuir–Szyszkowski isotherm [21]. The Marangoni stress can be evaluated if the velocity gradients at the interface and the viscosities of the two phases are known since the jump of the tangential viscous stress at the interface balances with the Marangoni stress. Hosokawa et al. [22] measured velocity fields about single contaminated drops of a glycerol–water solution falling through a silicone oil at low drop Reynolds numbers using SFV (spatiotemporal filter velocimetry) [23] to evaluate the distributions of the interfacial surfactant concentration (Figure 1). However, other field variables, e.g., the surfactant concentration in the drop and the pressure fields, required for full understanding of the phenomena, have not been obtained in the experiment yet.
In this study, the motions of contaminated drops of low Reynolds numbers are numerically simulated using the finite difference method under the same conditions as those in Hosokawa et al. [22] to discuss the distributions of the bulk and interfacial surfactant concentrations and the surfactant effects on the interface velocity and pressure in detail. In the following sections, the numerical method is first introduced. The validity of the numerical method is examined for clean drops, spherical particles and contaminated bubbles [10]. Then, the falling motions of contaminated drops are simulated, and the surfactant effects are discussed.

2. Numerical Method

2.1. Governing Equations

The governing equations are written in the frame of reference moving with a single drop. The drop shape is assumed to be spherical, and the drop Reynolds number is given as a parameter from the experimental data [22]. The flow passing through the spherical drop at the prescribed Reynolds number is therefore obtained by solving the governing equations. The ( r , θ ) two-dimensional spherical coordinates are used since the flow has axial symmetry at low drop Reynolds numbers. The flow is assumed to be incompressible and isothermal. The Navier–Stokes equation for Newtonian fluids is therefore given by [24]
ρ k u r , k t + u k · u r , k u θ , k 2 r = p k r + μ k 2 u r , k 2 u r , k r 2 2 r 2 sin θ sin θ u θ , k θ
ρ k u θ , k t + u k · u θ , k + u r , k u θ , k r = 1 r p k θ + μ k 2 u θ , k + 2 r 2 u r , k θ u θ , k r 2 sin 2 θ
where u r and u θ are the velocity components in the r and θ directions, respectively, t is the time, p is the pressure, ρ is the density, μ is the viscosity, and k denotes either the dispersed phase ( k = D ) or the continuous phase ( k = C ). The differential operators, u · and 2 , are given in Appendix A. The following continuity equation holds in each phase:
1 r 2 r r 2 u r , k r + 1 r sin θ ( sin θ u θ , k ) θ = 0
The normal component of the velocity is continuous at the interface when there is no interfacial mass transfer. Hence the velocity boundary condition at the drop interface is given by
u r , D = u r , C = 0
and the following no-slip condition is used for the tangential velocity component:
u θ , D = u θ , C = u θ , i n t
where the subscript i n t denotes the interface. The tangential viscous stress, τ r θ , in each phase balances with the Marangoni stress at the interface:
τ r θ i n t = 1 a σ θ
where [ X ] i n t represents the jump of the quantity, X, at the interface, σ is the surface tension, and a is the drop radius.
In the present study, surfactant is assumed to be soluble in the dispersed phase, so that it is present inside the drop at the molar concentration, C. The transport equation of C is given by
C t + · u D C = · D D C
where D D is the diffusion coefficient. Adsorption and desorption of surfactant molecules take place at the interface and the magnitude of the adsorption–desorption flux depends on C at the interface, C S , and the molar concentration, Γ , of surfactant adsorbing to the interface. The transport equation of Γ is given by [2,10,25]
Γ t + S · u S Γ = S · D S S Γ + S ˙ Γ
where D S is the diffusion coefficient of surfactant at the interface, S is the surface gradient operator, and u S is the interface velocity. The S ˙ Γ is the net adsorption–desorption flux, which is evaluated by making use of the following Frumkin–Levich model [9,10,26]:
S ˙ Γ = k a C S ( Γ m a x Γ ) k d Γ
where k a is the adsorption rate constant, k d is the desorption rate constant, and Γ m a x is the saturation value of Γ . The diffusive flux of C at the interface balances with S ˙ Γ , i.e.,
D D C r r = a = S ˙ Γ
The surface tension, σ , of the contaminated interface depends on Γ as [21]
σ = σ 0 + R G T Γ m a x ln 1 Γ Γ m a x
where σ 0 is the surface tension of a clean interface, R G is the universal gas constant, and T is the temperature.
The following dimensionless forms of the governing equations are useful to make clear the relevant dimensionless groups:
u r , k t + u k · u r , k u θ , k 2 r = p k r + 1 R e k 2 u r , k 2 u r , k r 2 2 r 2 sin θ sin θ u θ , k θ
u θ , k t + u k · u θ , k + u r , k u θ , k r = 1 r p k θ + 1 R e k 2 u θ , k + 2 r 2 u r , k θ u θ , k r 2 sin 2 θ
C t + · u D C = 2 C P e D
Γ t + S · u S Γ = S 2 Γ P e S + S ˙ Γ
where u = u / U , p = p / ρ U 2 , C = C / C 0 , Γ = Γ / Γ m a x , r = r / d , t = U t / d , = d , S = d S , and C 0 is the initial concentration. The Reynolds numbers, R e k , the Peclet number, P e D , and the surface Peclet number, P e S , are defined by R e k = ρ k U d / μ k , P e D = U d / D D , and P e S = U d / D S , respectively, where d is the drop diameter, and U is the drop velocity. The momentum jump condition in the tangential direction is given by
τ r θ i n t = M a a σ θ
where the Marangoni number, M a , and the dimensionless surface tension, σ , are defined by M a = R G T Γ m a x / μ C U and σ = σ / R G T Γ m a x , respectively. The dimensionless radius is given by a = a / d = 1 / 2 . The interfacial surfactant flux in the dimensionless form is given by
S ˙ Γ = α C S ( 1 Γ ) Γ L a
where C S = C S / C 0 , L a is the Langmuir number defined by L a = k a C 0 / k d , and α is the ratio of the adsorption rate to the characteristic velocity, i.e., α = k a C 0 d / U . The boundary condition of C at the interface is
C r r = a = P e D K S ˙ Γ
where K is the dimensionless adsorption length defined by K = Γ m a x / C 0 d . The σ is given by
σ = σ 0 + ln ( 1 Γ )
where σ 0 = σ 0 / R G T Γ m a x .

2.2. Discretization Schemes

2.2.1. Velocity-Pressure Coupling

The present numerical method for obtaining the velocity and pressure fields is based on the fractional step method [27,28]. The Adams–Bashforth scheme is used for the advection terms of the Navier–Stokes equations, while the Crank–Nicolson scheme is used for the diffusion terms. Therefore,
2 u r , k m 2 u r , k m r 2 2 R e k Δ t u r , k m 2 r 2 sin θ sin θ u θ , k m 1 θ = S r , k
2 u θ , k m u θ , k m r 2 sin 2 θ 2 R e k Δ t u θ , k m + 2 r 2 u r , k m 1 θ = S θ , k
where Δ t is the time step size and
S r , k = D r , k n + 2 R e k u r , k n Δ t + 3 A r , k n A r , k n 1 2 + p k n r
D r , k = 2 u r , k 2 u r , k r 2 2 r 2 sin θ sin θ u θ , k θ
A r , k = u k · u r , k u θ , k 2 r
S θ , k = D θ , k n + 2 R e k u θ , k n Δ t + 3 A θ , k n A θ , k n 1 2 + 1 r p k n θ
D θ , k = 2 u θ , k + 2 r 2 u r , k θ u θ , k r 2 sin 2 θ
A θ , k = u k · u θ , k + u r , k u θ , k r
The superscript n denotes the time step, A and D are the advection and diffusion terms, respectively, and S is the source term in the Crank–Nicolson solver. Equations (20) and (21) are solved for u m , where the superscript m is the step number for the iteration of the Crank–Nicolson solver. Since the former and latter include the terms of u θ and u r , respectively, in the diffusion term, these terms are treated as source terms in each iteration. The converged solutions of the velocity components are denoted as u .
The pressure contribution is then subtracted from u to obtain the intermediate velocity u , i.e.,
u r , k = u r , k + Δ t p k n r
u θ , k = u θ , k + Δ t r p k n θ
The pressure at n + 1 is obtained by solving the following Poisson equation:
2 p k n + 1 = Θ k Δ t
where
Θ k = 1 r 2 r r 2 u r , k r + 1 r sin θ ( sin θ u θ , k ) θ
The velocity components are then updated as
u r , k n + 1 = u r , k Δ t p k n + 1 r
u θ , k n + 1 = u θ , k Δ t r p k n + 1 θ
The staggered variable arrangement is used, i.e., the velocity components are located at the faces of computational cells, whereas the scalar variables are located at the cell centers. The advection terms are discretized using the ENO (essentially non-oscillatory) scheme [29]. The second-order centered difference scheme is used for the diffusion terms.

2.2.2. Interface Velocity and Numerical Treatment at Pole

The u θ , i n t is updated at each time step after obtaining u k n + 1 . The momentum jump condition in the tangential direction is given by
r r u θ , C r = κ r r u θ , D r M a r σ θ
where κ is the viscosity ratio defined by κ = μ D / μ C . Applying the continuity of the velocity at the interface and the second-order one-sided finite difference to the velocity gradients for the configuration in Figure 2 yields
u θ , i n t = k = C , D c k a Δ r k n Δ r k w Δ r k w Δ r k n f k 1 Δ r k n Δ r k w f k 2 M a 1 a σ θ k = C , D c k Δ r k n Δ r k w Δ r k w Δ r k n Δ r k n Δ r k w
where f = u θ / r , c C = 1 , c D = κ , and the superscripts for f denote the cell indices. The surface tension gradient is evaluated using the second-order centered difference scheme.
Discretization of some derivatives in the governing equations needs a special treatment at the pole ( r = 0 ). In the evaluation of the divergence of u and the Laplacian of the velocity components and the pressure, no special treatments are required since the terms r ϕ and r ϕ / r vanish at the pole, where ϕ = u r , u θ or p. On the other hand, for instance, the finite difference approximation of the advection term u r u θ / r requires u θ at r = 0 . A numerical treatment similar to that in Mohseni and Colonius [30] and Sugioka and Komori [31] is used in this study to set u θ at r = 0 , i.e., u θ ( 0 , θ ) = ( u θ ( Δ r / 2 , θ ) + u θ ( Δ r / 2 , θ + π ) ) / 2 (Figure 3). Although 0 θ π , u θ ( Δ r / 2 , θ + π ) can be set to the value of u θ ( Δ r / 2 , θ π ) because of the axial symmetry.

2.2.3. Surfactant Transport Equations

The transport equation of Γ is solved using the second-order Runge–Kutta scheme:
Γ ( 1 ) = Γ n + Δ t 2 S · u S n + 1 Γ n + S 2 Γ n P e S + S ˙ Γ ( Γ n )
Γ n + 1 = Γ n + Δ t S · u S n + 1 Γ ( 1 ) + S 2 Γ ( 1 ) P e S + S ˙ Γ ( Γ ( 1 ) )
The Runge–Kutta scheme is also used for the transport equation of C as follows:
C ( 1 ) = C n + Δ t 2 · u D n + 1 C n + 2 C n P e D
C n + 1 = C n + Δ t · u D n + 1 C ( 1 ) + 2 C ( 1 ) P e D
The advection and diffusion terms are discretized with the second-order ENO and second-order centered difference schemes, respectively.
The bulk concentration, C S , at the interface is evaluated by applying the second-order one-sided finite difference to C / r ,
C S = P e D K S ˙ Γ + 1 Δ r D n Δ r D w Δ r D w Δ r D n C 1 Δ r D n Δ r D w C 2 1 Δ r D n Δ r D w Δ r D w Δ r D n Δ r D n Δ r D w
and is used to evaluate S ˙ Γ . The S ˙ Γ is then used to set the boundary condition, Equation (18). The u θ , i n t is updated for S σ ( Γ n + 1 ) .

2.2.4. Solution Procedure

The solution procedure of the numerical method is summarized in the following:
  • Set the dimensionless groups ( ρ D / ρ C , κ , R e C , L a , P e D , P e S , M a , α , K and σ 0 ); generate computational grids inside and outside the drop; and set the initial conditions for each field variable.
  • Evaluate the Marangoni stress, S σ ( Γ n ) , and set the interface velocity, u θ , i n t n , using the boundary condition (Equation (35)).
  • Solve Equations (20) and (21) to obtain u .
  • Compute the intermediate velocity, u using Equations (28) and (29).
  • Solve the Poisson equation, Equation (30), to obtain p n + 1 .
  • Update the velocity to u n + 1 using Equations (32) and (33).
  • Update the interface velocity, u θ , i n t n + 1 , for u n + 1 .
  • Compute the surfactant concentration, C S n , at the interface using Equation (40) with C n .
  • Solve Equations (36) and (37) to obtain Γ n + 1 , where the adsorption–desorption flux, S ˙ Γ , is evaluated using Equation (17).
  • Set the boundary condition of C at the interface, Equation (18), with S ˙ Γ .
  • Solve Equations (38) and (39) to obtain C n + 1 .
  • Return to Step 2.

3. Validation of Numerical Method

3.1. Clean Drop and Solid Sphere

The numerical method is validated for simple problems with well-accepted correlations, i.e., uniform flows about solid spheres and clean spherical drops. Figure 4 shows the computational domain and the grid system. The dimensionless particle diameter is unity, and the domain size is 40 in the radial direction. The outer boundary of the domain is split into either inflow or outflow. The zero-gradient condition is used for the velocity at the outflow boundary and the outlet pressure is set to zero. The numbers of the computational cells in the r and θ directions are 25 and 120, respectively. The cell size is uniform in the θ direction, whereas in the r direction the cell size decreases (increases) toward the drop interface (the outer boundary of the domain) as r increases. For clean drops, it is not necessary to use the very thin computational cells in the vicinity of the drop interface for capturing the velocity boundary layer at low and moderate Reynolds numbers. On the other hand, in the contaminated drop cases discussed in Section 4, the thin cells are required to capture concentration boundary layers of C forming at the interface. Cuenot et al. [10] estimated the dimensionless thickness of the concentration boundary layer as δ π / ( 4 P e C ) [9] and generated computational grids so as to have at least three computational cells within the boundary layer. The present grid is prepared to meet the same requirement. This grid will also be used in the contaminated drop simulations. The mesh dependence of the predictions was confirmed to be small enough at the present spatial resolution.
Figure 5 shows predicted u θ , i n t of clean drops at several drop Reynolds numbers, R e C , the range of which ( 0.1 R e C 20 ) covers R e C in the experiment of the contaminated drops [22], i.e., 0.5 R e C 0.73 . The viscosity ratio in the simulations is κ = 0.021 . The predicted u θ , i n t at R e C = 0.1 agrees well with the Hadamard–Rybczynski (HR) solution, u θ , i n t = sin θ / ( 2 ( 1 + κ ) ) . The slight difference between the prediction and the analytical solution for the Stokes flow is attributed to the presence of a weak inertial effect even at the small Reynolds number. The increase in R e C up to R e C = 1 makes the inertial effect obvious, i.e., u θ , i n t for θ < 3 π / 4 is larger than the HR solution. The deviation from the Stokes flow is increased by further increase in R e C .
Figure 6 shows predicted C D of clean spherical drops. The open circles are for κ = 0.021 and they are compared with the HR drag (the black dashed line)
C D = 8 R e C 2 + 3 κ 1 + κ
and the empirical correlation proposed by Myint et al. [32] (the black solid line)
C D = 8 R e C 2 + 3 κ 1 + κ 1 + 0.15 R e C 0.687
The C D at R e C = 0.1 agrees with both drag correlations. The predictions become larger than the HR drag as R e C increases due to the inertial effect. The trend of the drag curve of Equation (42) agrees with that of the data, the values of which are however smaller than those of the correlation since the correlation accounts for the increase in drag by shape deformation. The green and blue symbols in the figure are predictions for κ = 1 and 5, respectively. The simulations reproduce well the increase in C D by the factor of ( 2 + 3 κ ) / ( 1 + κ ) . Simulations of flows about spherical particles are also carried out as the limiting case of κ by directly applying the no-slip boundary condition to the particle surface. The predicted drag coefficients are shown in the figure by the red symbols and compared with the following empirical correlation proposed by Schiller and Naumann [33]:
C D = 24 R e C 1 + 0.15 R e C 0.687
The data agree well with the correlation. Thus, the present numerical method gives good predictions of the motion of drops for wide ranges of R e C and κ .

3.2. Contaminated Bubbles

Cuenot et al. [10] carried out numerical simulations of contaminated spherical bubbles in stagnant liquid at R e C = 100 , P e C = 1 × 10 5 and K = 1 . They changed M a , α and L a to deal with different situations, i.e., the stagnant-cap regime and the fully covered regime. Case 1 ( M a = 61 , α = 0.001 , L a = 0.112 ), Case 3 ( M a = 61 , α = 0.001 , L a = 0.0112 ) and Case 4 ( M a = 610 , α = 0.001 , L a = 0.112 ) in their study are simulated with the present numerical method. The flows inside bubbles are not solved; instead, the tangential viscous stress in the gas phase is set to zero in Equation (6). Figure 7 shows predicted u θ , i n t and Γ , which are compared with the predictions by Cuenot et al. [10]. The u θ , i n t increases with increasing θ , whereas it steeply decreases to zero at θ π / 4 and 5 π / 8 in Cases 1 and 3, respectively. On the other hand, the gradients of Γ are large at these θ . The steep decrease in u θ , i n t is due to the Marangoni stress caused by the gradient of Γ . In comparison with Case 3, Γ in the contaminated region is larger in Case 1 and the attenuation of u θ , i n t is more significant because of the larger L a . In Case 4, u θ , i n t is almost zero, in other words the interface is immobile, due to the large M a , and Γ 0.09 even at the bubble nose. Thus, being similar to the predictions of Cuenot et al. [10], the effects of L a and M a on u θ , i n t and Γ are clearly observed.

4. Simulations of Contaminated Drops

The main difference in the simulations of contaminated drops carried out in this section is the value of L a , which ranges from 0 to 10 2 , leading to varying surface concentrations. The fluid properties of a drop and a bulk liquid are the same as those in the experiment [22], i.e., those for a glycerol–water solution of 54 wt.% and a silicone oil. They are summarized in Table 1. The drop diameter is 8.3 mm. Under this condition, R e C is less than unity. Triton X-100 is used for surfactant. The adsorption–desorption properties and the diffusion coefficients of Triton X-100 are shown in Table 2 (see Appendix B for evaluation methods). The C 0 in a drop is initially set either to 0 (clean liquid), 0.0020 , 0.0050 , 0.010 or 0.10 mol/m 3 . These conditions are referred to as Cases 1, 2, 3, 4 and 5, respectively. The concentrations are less than the critical micellar concentration (0.28 mol/m 3 ). The measurements were carried out at 250 mm below the nozzle tip. The dimensionless time when drops reached the measurement region was 30. The relevant dimensionless groups in each case are summarized in Table 3. These values of the dimensionless groups are used as the input parameters of the numerical simulations. See Hosokawa et al. [22] for further details of the experimental conditions.
Figure 8a shows flow fields in Case 1 ( C 0 = 0 mol/m 3 ). The left half is the flow field visualized by integrating the grayscale values in high-speed video images while tracking the drop; the flow structure is therefore observed from the frame of reference moving with the drop. The liquid in the outside of the drop flows from the nose ( θ = 0 ) side to the rear ( θ = π ) side and no separation of the flow takes place. Inside the drop, an internal circulation is formed. Hence the interface is mobile. The right half of the figure shows the predicted velocity vectors inside and outside the drop. The number of the velocity vectors are reduced for visualization purposes. The flow characteristics inside and outside the drop agree with the measurement.
Figure 9 shows Γ at t = 10 , 20 and 30. In Case 2 ( C 0 = 0.0020 mol/m 3 ), Γ increases in the rear region of the drop as t increases. The contaminated area also increases, but is not so large even at t = 30 . In Case 3 ( C 0 = 0.0050 mol/m 3 ) and 4 ( C 0 = 0.010 mol/m 3 ), the trend of the temporal change in Γ is similar to that in Case 1, whereas the contaminated area becomes much wider as C 0 increases. The Γ in Case 5 ( C 0 = 0.10 mol/m 3 ) is much larger than the other cases, and the whole interface is contaminated. The presence of surfactant reduces the region of the internal circulation both in the experiment and simulation of Case 3 as shown in Figure 8b (the flow field of Case 2, whose flow characteristics are in-between those of Cases 1 and 3, was omitted for saving the space). The magnitude of the drop velocity in the vicinity of the drop rear is very small. The increase in C 0 from 0.0050 mol/m 3 to 0.010 mol/m 3 (Case 4) makes the circulation region smaller (Figure 8c). The interface velocity in the rear half, θ > π / 2 , is close to zero, so that the formation of velocity boundary layer along the interface can be clearly observed in the continuous phase. In Case 5 (Figure 8d), no internal circulation is formed inside the drop. In the continuous phase, the grayscale values near the interface in the measurement are very low (almost black), which implies that the velocities there are very small. The prediction agrees with the measurement; the interface velocity is almost zero at the whole interface; in other words, the interface is fully contaminated and immobile, so that the velocity field outside the drop seems to be similar to that about a solid sphere.
Figure 10a shows the distribution of C in Case 2. The t increases from the left to the right. The rightmost figures are magnified views of the regions labeled by (i) and (ii). The C near the interface in the rear region decreases at t = 5 . The low C region develops toward the centerline of the drop as t increases due to the advection. However C near the rear does not decrease so much. Because of the low diffusion coefficient, the pattern of the low C region is determined by the advection and does not smear out. Figure 11 shows the time evolution of C S in Case 2. The surfactant flux from the bulk to the interface decreases C in the vicinity of the interface and the weak diffusion cannot compensate for the reduction in C by adsorption. The C S at t = 5 is therefore smaller than the initial value, C S = 1 . The θ for the minimum C S corresponds to that for steep velocity decrease; in other words, the stagnant-cap angle, θ c a p . For θ < θ c a p , Γ is very small, so that the adsorption is dominant in the surfactant flux. The increase in C S with increasing θ for θ > θ c a p implies that the desorption increases C S in the stagnant-cap region. In the stagnant-cap region, Γ approaches a value for the local adsorption–desorption equilibrium, Γ 1 / ( 1 + ( C S ) 1 ) . The u θ , i n t is very small in that region, but is still non-zero. The weak interface compression effect ( Γ S · u S ) appears and slightly increases Γ . However, some excess amount of surfactant is desorbed from the interface to the bulk. Thus, C S has an increasing trend in the high θ region. At larger t , C S decreases moderately for θ < θ c a p , whereas for θ > θ c a p C S increases and becomes larger than unity.
Figure 10b shows the time evolution of C in Case 4. As shown in the magnified view of C near the interface (iii), a low C region is formed around the drop equator. The interface area is more than half covered with surfactant at this moment (Figure 9c). The formation of the low C region is therefore the result of surfactant adsorption from the bulk to the interface. The thickness of the concentration boundary layer is, as expected, very thin because of the high P e D . In Case 5, the whole interface is covered by surfactant as shown in Figure 9d, which results in the formation of the thin concentration boundary layer of low C in the nose region (Figure 10c(v)). As it is similar to Case 2, the interface compression effect increases Γ in the rear half of the drop and the desorption becomes dominant near the rear, resulting in the formation of boundary layer of C larger than unity in the rear region (the magnified view (iv)).
Figure 12 shows the pressure distribution in Case 3, in which p D and p C are subtracted by the pressure, p 0 , at the drop center and the far field pressure, p = 0 , respectively. The θ c a p represents the cap angle evaluated from the velocity distribution and θ c a p 0.6 π rad. In the pressure distribution, high-pressure spots are formed in the vicinity of the interface at θ slightly smaller than θ c a p . Comparing p with the velocity field clearly shows that these high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region in the rear half of the drop and the attenuation of the tangential velocity in the continuous phase.
Figure 13a shows the interface velocity, u θ , i n t , at t = 30 . The lines and symbols are the predictions and the measured data [22], respectively. The predicted velocity of the clean drop (Case 1) agrees fairly with the data, except for the rear region. The lower values at large θ in the data are due to slight contamination in the rear region, the cause of which would be contamination brought into the drop by seeding particles. The retardation effect of surfactant in u θ , i n t is clearly seen in both data and predictions of Cases 2–5. In Case 2, the predicted velocity profile agrees with the data. The u θ , i n t decreases from θ π / 2 as θ increases and becomes very small for θ > 3 π / 4 , i.e., θ c a p 3 π / 4 . The increase in C 0 decreases θ c a p , i.e., θ c a p are 0.6 π and 0.45 π rad in Case 3 and 4, respectively. Although the predicted velocities are somewhat larger than the data, the dependence of u θ , i n t on C 0 is reproduced. In Case 5, u θ , i n t 0 at any θ as already seen in Figure 8.
The distributions of Γ are shown in Figure 13b. The predicted Γ of Case 2 is almost zero for θ < θ c a p , whereas it steeply increases and reaches a large value for θ > θ c a p , which is larger than L a / ( 1 + L a ) = 0.67 for the absorption-desorption equilibrium for C 0 since C S in this region is larger than unity (Figure 11). The measured Γ shows a similar profile, while the following differences are present: the measured Γ is smoother around θ = θ c a p than predicted and approaches 0.67 . The increase in C 0 (Cases 3 and 4) increases Γ , whereas the drop nose is still almost clean. On the other hand, Γ takes a large value even at the nose in Case 5 since the adsorption is much more dominant than the desorption as represented by the large L a .
As shown in Figure 13c, peaks appear in the predicted Marangoni stresses in Cases 2, 3 and 4. However, in Case 5, there are no peaks. This is due to the smaller Γ gradient. Despite a smaller Γ gradient, the Marangoni stress immobilizes the interface almost completely since the Marangoni stress acts on the whole interface. In this situation, the Marangoni stress agrees very well with the viscous stress acting on the solid sphere given by the Stokes solution (the black dotted line in the figure), i.e., u θ / r r = a = 3 sin θ . The smaller peaks in the data are attributed to a smoother distribution of Γ .
The peaked distribution of the Marangoni stress also causes similar peaks in the pressure distribution as shown in Figure 13d. The locations of the strong Marangoni stress correspond to those of the steep velocity reduction (the cap angles). The peaks of the pressure however locate at slightly smaller θ . This trend was already observed in Figure 12. The pressure distribution in Case 5 is smooth as in the clean drop case and agrees well with the Stokes solution given by p C = 3 cos θ / ( 2 R e C ) .
In the distributions of C S shown in Figure 13d, θ for the minimum C S in Cases 2–4 roughly correspond to θ c a p . The C S in these cases exhibit the increasing trend in the stagnant-cap region. On the other hand, on the completely immobilized interface in Case 5, the profile of C S monotonically increases with increasing θ .
The measured Γ is much smoother than predicted, which suggests that in reality P e S is much lower than that used in the simulation or that there is some other transport process much more effective for interface transport. A brief discussion on this point is given in Appendix C.

5. Conclusions

The motion of contaminated spherical drops falling through a stagnant liquid were numerically simulated using the finite difference method to investigate the effects of surfactant on the flow inside and outside the drops. The numerical method was first validated through simulations of spherical particles, clean drops and contaminated bubbles. The predicted drag coefficients of the spherical particles and clean drops agreed well with available drag correlations. The interface velocity and the interfacial surfactant concentration of contaminated bubbles also agreed with those predicted in the literature [10]. The numerical simulations of contaminated drops were then carried out. The numerical conditions were the same as those in the experiment carried out in Hosokawa et al. [22], i.e., the fluid properties for the drop and the bulk liquid were for a glycerol–water solution and a silicone oil, respectively, and the surfactant was Triton X-100. The drop Reynolds numbers were less than unity and the drop shape kept spherical. The surfactant concentration, C 0 , was varied below the critical micellar concentration; C 0 = 0.0020 , 0.0050 , 0.010 and 0.10 mol/m 3 in Cases 1, 2, 3, 4 and 5, respectively. The characteristics of the field variables can be summarized as follows:
  • The predicted interfacial surfactant concentration, Γ , is almost zero for the angle, θ , smaller than a certain value ( θ c a p ), whereas it steeply increases and reaches a large value for θ > θ c a p . The interface velocity, u θ , i n t , decreases as θ increases and becomes very small for θ > θ c a p . The increase in C 0 increases Γ and decreases θ c a p . The presence of surfactant attenuates the internal circulation and the increase in C 0 makes the circulation region smaller.
  • The surfactant flux from the bulk to the interface decreases C in the vicinity of the interface and the weak diffusion cannot compensate for the reduction in C by adsorption. The bulk concentration, C S , at the interface therefore tends to be smaller than C 0 for θ < θ c a p . The pattern of the low C region is determined by the advection and does not smear out because of a small diffusive flux.
  • Peaks appear in the predicted Marangoni stresses in Cases 2–4, while in Case 5 no peaks develop due to a smaller Γ gradient. The peaked distribution of the Marangoni stress also causes similar peaks in the pressure distribution. The locations of the strong Marangoni stress correspond to θ c a p , whereas the peaks of the pressure appear at slightly smaller θ . The high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region in the rear half of the drop and the attenuation of the tangential velocity in the continuous phase.

Author Contributions

K.H.: conceptualization, investigation, methodology, validation, visualization, writing—original draft preparation, writing—review and editing; Y.M.: validation, visualization, writing—original draft preparation; M.J.A.v.d.L.: investigation, validation, visualization, writing—review and editing; N.G.D.: investigation, validation, writing—review and editing; S.H.: investigation, visualization, writing—review and editing; A.T.: conceptualization, investigation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by JSPS KAKENHI, Grant Numbers JP 20K04267 and JP 19H02067.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Differential Operators

The differential operators used in the governing equations are defined by
u · ( ) = u r ( ) r + u θ r ( ) θ
2 ( ) = 1 r 2 r r 2 ( ) r + 1 r sin θ θ 1 r ( ) θ sin θ

Appendix B. Surfactant Property

The isotherm is given in terms of C 0 :
σ = σ 0 R G T Γ m a x ln 1 + C 0 χ
where χ = k d / k a At large C 0 , this equation can be approximated as
σ σ 0 R G T Γ m a x ln C 0 χ
The Γ m a x is therefore calculated by fitting the expression
Γ m a x = 1 R G T σ ln C 0 T
to data of σ ( C 0 ) , which can be measured using a sessile drop method (Figure A1).
Figure A1. Surface tension [22]. Solid line: Equation (A3) and σ ( C M C ) .
Figure A1. Surface tension [22]. Solid line: Equation (A3) and σ ( C M C ) .
Fluids 07 00055 g0a1
Substituting Γ m a x and the interfacial surfactant concentration
Γ e q = C 0 Γ m a x χ + C 0
in the adsorption–desorption equilibrium state into Equation (11) yields
σ ( C 0 ) = σ 0 + R G T Γ m a x ln χ χ + C 0
Fitting this functional form to the σ data determines χ .
By eliminating the advection and diffusion terms in Equation (8) for a static interface with a uniform Γ distribution, Equation (8) can be immediately integrated, i.e.
Γ ( t ) = C 0 Γ m a x χ + C 0 1 e ( χ + C 0 ) k a t
where t is the time measured from the instant at which a drop intrudes into the bulk liquid. Substituting thus obtained Γ ( t ) into Equation (A3) yields
σ ( t ) = σ 0 + R G T Γ m a x ln χ + C 0 e ( χ + C 0 ) k a t χ + C 0
where only k a is the unknown constant. The least-square fitting to the data of σ ( t ) gives k a .
The change of σ in time at a stationary interface is given by [34]
σ ( t ) = σ ( Γ e q ) + R G T Γ e q 2 2 C 0 π D D t
The D D can therefore be calculated by applying the sessile drop method. Since there is no knowledge of D S of Triton X-100 at the interface between the glycerol–water solution and silicone oil, it is assumed that D S = D D [2,10].

Appendix C. Contaminated Drops of Lower Peclet Numbers

The measured Γ is much smoother than predicted as shown in Figure 13, which leads to a speculation that in reality the diffusion of surfactant is much larger; in other words, the Peclet number is much lower than that in the simulation. As described in Appendix B, D D was evaluated from the data of σ obtained by the sessile drop method (Equation (A10)). By considering the measurement uncertainties in D D we increased the value of D D in the simulation by three times. The P e D for the increased D D is given in Table A1. The simulations of the contaminated drops were carried out with these P e D and the assumption of P e S = P e D . The predicted Γ profiles were still much sharper than those in the measurement although the reduction in the Peclet numbers slightly affected the results. By further decreasing P e S down to the order of O ( 10 ) (Table A1), we obtained better agreements between the predictions and the data as shown in Figure A2. The peaks in the pressure and the Marangoni stress observed in Figure 13 are mitigated.
The simulation thus implies the possibility of lower P e in reality. However, according to values of D D and D S measured for different contaminated systems [35,36], the assumption, D D = D S , is considered to be acceptable, and P e S of O ( 10 ) seems too small. Further studies are required to make clearer the cause of the deviation from the measurement.
Table A1. Reduced Peclet numbers.
Table A1. Reduced Peclet numbers.
Case2345
C 0  [mol/m 3 ]0.00200.00500.0100.10
P e D × 10 5 3.53.02.72.6
P e S × 10 1 7.06.15.45.2
Figure A2. (a) u θ , i n t , (b) Γ , (c) M a σ / ( a θ ) , (d) p C and (e) C S predicted at lower Peclet numbers (see Table A1). The legends are the same as those in Figure 13.
Figure A2. (a) u θ , i n t , (b) Γ , (c) M a σ / ( a θ ) , (d) p C and (e) C S predicted at lower Peclet numbers (see Table A1). The legends are the same as those in Figure 13.
Fluids 07 00055 g0a2

References

  1. Clift, R.; Grace, J.R.; Weber, M.E. Bubbles, Drops, and Particles; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  2. Takagi, S.; Uda, T.; Watanabe, Y.; Matsumoto, Y. Behavior of a rising bubble in water with surfactant dissolution (1st report, steady behavior). Trans. JSME 2003, 69, 2192–2199. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  3. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Butterworth-Heinemann: Oxford, UK, 1987. [Google Scholar]
  4. Aoki, J.; Hayashi, K.; Hosokawa, S.; Tomiyama, A. Effects of surfactant on mass transfer from single carbon dioxide bubbles in vertical pipes. Chem. Eng. Technol. 2015, 38, 1955–1964. [Google Scholar] [CrossRef]
  5. Boussinesq, J. Vitesse de la chute lente, devenue uniforme, d’une goutte liquide spherique. dans un fluid visqueux de poids specifique moindre. C. R. Acad. Sci. 1913, 156, 1124–1130. [Google Scholar]
  6. Scriven, L.E. Dynamics of a fluid interface. Chem. Eng. Sci. 1960, 12, 98–108. [Google Scholar] [CrossRef]
  7. Hadamard, J.S. Mouvement permanent tent d’une sphere liquide et visquese dans un liquide visqueux. Comptes Rendus I’Acad. Des Sci. 1911, 152, 1735–1738. [Google Scholar]
  8. Rybczynski, W. On the translatory of a fluid sphere in a viscous medium. Bull. Acad. Sci. Kracow Ser. 1911, 3, 40–46. [Google Scholar]
  9. Levich, V.G. Physicochemical Hydrodynamics; Prentice Hall: Englewood Cliffs, NJ, USA, 1962. [Google Scholar]
  10. Cuenot, B.; Magnaudet, J.; Spennato, B. The effects of slightly soluble surfactants on the flow around a spherical bubble. J. Fluid Mech. 1997, 339, 25–53. [Google Scholar] [CrossRef]
  11. Takagi, S.; Ogasawara, T.; Fukuta, M.; Matsumoto, Y. Surfactant effect on the bubble motions and bubbly flow structures in a vertical channel. Fluid Dyn. Res. 2009, 41, 065003. [Google Scholar] [CrossRef]
  12. Muradoglu, M.; Tryggvason, G. A front-tracking method for computation of interfacial flows with soluble surfactants. J. Comput. Phys. 2008, 227, 2238–2262. [Google Scholar] [CrossRef]
  13. Hayashi, K.; Tomiyama, A. Effects of surfactant on terminal velocity of a Taylor bubble in a vertical pipe. Int. J. Multiph. Flow 2012, 39, 78–87. [Google Scholar] [CrossRef]
  14. Hayashi, K.; Tomiyama, A. Effects of numerical treatment of viscous and surface tension forces. J. Comput. Multiph. Flows 2014, 6, 111–126. [Google Scholar] [CrossRef] [Green Version]
  15. Engberg, R.F.; Wegener, M.; Kenig, E.Y. The influence of Marangoni convection on fluid dynamics of oscillating single rising droplets. Chem. Eng. Sci. 2014, 117, 114–124. [Google Scholar] [CrossRef]
  16. Atasi, O.; Haut, B.; Pedrono, A.; Scheid, B.; Legendre, D. Influence of soluble surfactants and deformation on the dynamics of centered bubbles in cylindrical microchannels. Langmuir 2018, 34, 10048–10062. [Google Scholar] [CrossRef] [PubMed]
  17. Shin, S.; Chergui, J.; Juric, D.; Kahouadji, L.; Matar, O.; Craster, R. A hybrid interface tracking-level set technique for multiphase flow with soluble surfactant. J. Comput. Phys. 2018, 359, 409–435. [Google Scholar] [CrossRef] [Green Version]
  18. Piedfert, A.; Lalanne, B.; Masbernat, O.; Risso, F. Numerical simulations of a rising drop with shape oscillations in the presence of surfactants. Phys. Rev. Fluids 2018, 3, 103605. [Google Scholar] [CrossRef] [Green Version]
  19. Batchvarov, A.; Kahouadji, L.; Magnini, M.; Constante-Amores, C.; Shin, S.; Chergui, J.; Juric, D.; Craster, R.; Matar, O. Effect of surfactant on elongated bubbles in capillary tubes at high Reynolds number. Phys. Rev. Fluids 2020, 5, 093605. [Google Scholar] [CrossRef]
  20. Ubal, S.; Brown, N.; Lu, J.; Corvalan, C. Active motion of contaminated microbubbles. Chem. Eng. Sci. 2021, 238, 116574. [Google Scholar] [CrossRef]
  21. Dukhin, S.S.; Kretzschmar, G.; Miller, R. Dynamics of Adsorption at Liquid Interfaces; Elsevier Science: Amsterdam, The Netherlands, 1995. [Google Scholar]
  22. Hosokawa, S.; Masukura, Y.; Hayashi, K.; Tomiyama, A. Experimental evaluation of Marangoni stress and surfactant concentration at interface of contaminated single spherical drop using spatiotemporal filter velocimetry. Int. J. Multiph. Flow 2017, 97, 157–167. [Google Scholar] [CrossRef] [Green Version]
  23. Hosokawa, S.; Tomiyama, A. Spatial Filter Velocimetry based on Time-Series Particle Images. Exp. Fluids 2012, 52, 1361–1372. [Google Scholar] [CrossRef]
  24. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge University Press: Cambridge, UK, 1967. [Google Scholar]
  25. Stone, H.A. A simple derivation of the time-dependent convective-diffusion equation for surfactant transport along a deforming interface. Phys. Fluids 1990, 2, 111–112. [Google Scholar] [CrossRef]
  26. Frumkin, A.; Levich, V.G. On surfactants and interfacial motion. Zhurnal Fiz. Khimii 1947, 21, 1183–1204. (In Russian) [Google Scholar]
  27. Kim, J.; Moin, P. Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations. J. Comput. Phys. 1985, 323, 308–323. [Google Scholar] [CrossRef]
  28. Choi, H.; Moin, P. Effects of the computational time step on numerical solutions of turbulent flow. J. Comput. Phys. 1994, 113, 1–4. [Google Scholar] [CrossRef]
  29. Shu, C.; Osher, S. Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 1988, 77, 439–471. [Google Scholar] [CrossRef] [Green Version]
  30. Mohseni, K.; Colonius, T. Numerical treatment of polar coordinate singularities. J. Comput. Phys. 2000, 157, 787–795. [Google Scholar] [CrossRef] [Green Version]
  31. Sugioka, K.; Komori, S. Drag and lift forces acting on a spherical water droplet in homogeneous linear shear air flow. J. Fluid Mech. 2007, 570, 155–175. [Google Scholar] [CrossRef]
  32. Myint, W.; Hosokawa, S.; Tomiyama, A. Terminal velocity of single drops in stagnant liquids. J. Fluid Sci. Technol. 2006, 1, 72–81. [Google Scholar] [CrossRef] [Green Version]
  33. Schiller, L.; Naumann, A. Über die grundlegende Berechnungen bei der Schwerkraft-aufbereitung. Z. Des Vereines Dtsch. Ingenieure 1933, 77, 318–320. [Google Scholar]
  34. Ward, A.F.H.; Tordai, L. Time-dependence of Boundary Tensions of Solutions: I. the role of diffusion in time-effects. J. Chem. Phys. 1946, 14, 453–461. [Google Scholar] [CrossRef]
  35. Valkovska, D.S.; Danov, K.D. Determination of Bulk and Surface Diffusion Coefficients from Experimental Data for Thin Liquid Film Drainage. J. Colloid Interface Sci. 2000, 223, 314–316. [Google Scholar] [CrossRef] [Green Version]
  36. Morgan, C.E.; Breward, C.J.W.; Griffiths, I.M.; Howell, P.D.; Penfold, J.; Thomas, R.K.; Tucker, I.; Petkov, J.T.; Webster, J.R.P. Kinetics of Surfactant Desorption at an Air-Solution Interface. J. Colloid Interface Sci. 2012, 28, 17339–17348. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Single drop falling through stagnant liquid [22], where d is the drop diameter, a is the drop radius, r and θ are the radial and azimuthal coordinates, U is the drop velocity, and g is the magnitude of the gravitational acceleration.
Figure 1. Single drop falling through stagnant liquid [22], where d is the drop diameter, a is the drop radius, r and θ are the radial and azimuthal coordinates, U is the drop velocity, and g is the magnitude of the gravitational acceleration.
Fluids 07 00055 g001
Figure 2. Interface velocity.
Figure 2. Interface velocity.
Fluids 07 00055 g002
Figure 3. Numerical treatment at pole.
Figure 3. Numerical treatment at pole.
Fluids 07 00055 g003
Figure 4. Computational grid. (a) whole domain, (b) drop, (c) continuous phase around drop.
Figure 4. Computational grid. (a) whole domain, (b) drop, (c) continuous phase around drop.
Fluids 07 00055 g004
Figure 5. Interface velocities of clean drops at several Reynolds numbers. HR: Hadamard-Rybczynski solution.
Figure 5. Interface velocities of clean drops at several Reynolds numbers. HR: Hadamard-Rybczynski solution.
Fluids 07 00055 g005
Figure 6. Drag coefficients, C D , of clean drops and solid spheres. HR: Hadamard–Rybczynski [7,8] (Equation (41)), WM: Win Myint [32] (Equation (42)).
Figure 6. Drag coefficients, C D , of clean drops and solid spheres. HR: Hadamard–Rybczynski [7,8] (Equation (41)), WM: Win Myint [32] (Equation (42)).
Fluids 07 00055 g006
Figure 7. Interface velocities, u θ , i n t , and interfacial surfactant concentrations, Γ of contaminated bubbles. Steady-state solutions for Cases 1, 3 and 4 in Cuenot et al. [10].
Figure 7. Interface velocities, u θ , i n t , and interfacial surfactant concentrations, Γ of contaminated bubbles. Steady-state solutions for Cases 1, 3 and 4 in Cuenot et al. [10].
Fluids 07 00055 g007
Figure 8. Comparisons of flow fields between predictions and experimental data [22]. (a) Case 1 (clean), (b) Case 3, (c) Case 4, (d) Case 5. Right: predicted; Left: measured (Reprinted from International Journal of Multiphase Flow 97, S. Hosokawa, Y. Masukura, K. Hayashi, A. Tomiyama, Experimental evaluation of Marangoni stress and surfactant concentration at interface of contaminated single spherical drop using spatiotemporal filter velocimetry, 157–167, Copyright (2017), with permission from Elsevier).
Figure 8. Comparisons of flow fields between predictions and experimental data [22]. (a) Case 1 (clean), (b) Case 3, (c) Case 4, (d) Case 5. Right: predicted; Left: measured (Reprinted from International Journal of Multiphase Flow 97, S. Hosokawa, Y. Masukura, K. Hayashi, A. Tomiyama, Experimental evaluation of Marangoni stress and surfactant concentration at interface of contaminated single spherical drop using spatiotemporal filter velocimetry, 157–167, Copyright (2017), with permission from Elsevier).
Fluids 07 00055 g008
Figure 9. Time evolutions of Γ . (a) Case 2; (b) Case 3; (c) Case 4; (d) Case 5.
Figure 9. Time evolutions of Γ . (a) Case 2; (b) Case 3; (c) Case 4; (d) Case 5.
Fluids 07 00055 g009
Figure 10. Time evolutions of C . (a) Case 2; (b) Case 4; (c) Case 5. The time t increases from left to right, whereas the rightmost figures are magnified views of the distribution at t = 30 .
Figure 10. Time evolutions of C . (a) Case 2; (b) Case 4; (c) Case 5. The time t increases from left to right, whereas the rightmost figures are magnified views of the distribution at t = 30 .
Fluids 07 00055 g010
Figure 11. Surfactant concentration, C S , at interface in Case 2.
Figure 11. Surfactant concentration, C S , at interface in Case 2.
Fluids 07 00055 g011
Figure 12. Pressure distributions inside and outside drop in Case 3 (right). The velocity field is also shown for comparing the high-pressure location with the cap angle, θ c a p 0.6 π rad (left).
Figure 12. Pressure distributions inside and outside drop in Case 3 (right). The velocity field is also shown for comparing the high-pressure location with the cap angle, θ c a p 0.6 π rad (left).
Fluids 07 00055 g012
Figure 13. Predictions of (a) u θ , i n t , (b) Γ , (c) M a σ / ( a θ ) , (d) p C and (e) C S at interface. The symbols and lines are measured [22] and predicted values, respectively. The black dotted line in σ / ( a θ ) represents the dimensionless velocity gradient, u θ / r r = a = 3 sin θ , given by the Stokes solution for solid sphere. The black dotted line in p C is also the Stokes solution given by p C = 3 cos θ / ( 2 R e C ) .
Figure 13. Predictions of (a) u θ , i n t , (b) Γ , (c) M a σ / ( a θ ) , (d) p C and (e) C S at interface. The symbols and lines are measured [22] and predicted values, respectively. The black dotted line in σ / ( a θ ) represents the dimensionless velocity gradient, u θ / r r = a = 3 sin θ , given by the Stokes solution for solid sphere. The black dotted line in p C is also the Stokes solution given by p C = 3 cos θ / ( 2 R e C ) .
Fluids 07 00055 g013
Table 1. Fluid properties of glycerol–water solution of 54 wt.% and silicone oil at 298 K [22].
Table 1. Fluid properties of glycerol–water solution of 54 wt.% and silicone oil at 298 K [22].
PhaseDispersedContinuous
ρ [kg/m 3 ]1132967
μ [Pa·s]0.006130.299
σ 0 [N/m]0.0343
Table 2. Properties of Triton X-100 [22].
Table 2. Properties of Triton X-100 [22].
k a [m 3 /mol·s] k d [1/s] Γ max [mol/m 2 ] D D [m 2 /s] D S [m 2 /s]
9.40.00961.69 × 10 6 2 × 10 10 2 × 10 10
Table 3. Dimensionless groups.
Table 3. Dimensionless groups.
Case12345
C 0  [mol/m 3 ]00.00200.00500.0100.10
ρ D / ρ C 1.171.171.171.171.17
κ 0.0210.0210.0210.0210.021
R e D 4239343029
R e C 0.730.680.590.520.50
L a 02.04.99.898
P e D × 10 6 -1.10.910.800.77
P e S × 10 6 -1.10.910.800.77
M a -0.550.640.720.75
α -0.00620.0180.0400.42
K-0.100.0410.0200.0020
σ 0 8.28.28.28.28.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hayashi, K.; Motoki, Y.; van der Linden, M.J.A.; Deen, N.G.; Hosokawa, S.; Tomiyama, A. Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers. Fluids 2022, 7, 55. https://doi.org/10.3390/fluids7020055

AMA Style

Hayashi K, Motoki Y, van der Linden MJA, Deen NG, Hosokawa S, Tomiyama A. Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers. Fluids. 2022; 7(2):55. https://doi.org/10.3390/fluids7020055

Chicago/Turabian Style

Hayashi, Kosuke, Yuya Motoki, Matheus J. A. van der Linden, Niels G. Deen, Shigeo Hosokawa, and Akio Tomiyama. 2022. "Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers" Fluids 7, no. 2: 55. https://doi.org/10.3390/fluids7020055

APA Style

Hayashi, K., Motoki, Y., van der Linden, M. J. A., Deen, N. G., Hosokawa, S., & Tomiyama, A. (2022). Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers. Fluids, 7(2), 55. https://doi.org/10.3390/fluids7020055

Article Metrics

Back to TopTop