M Motaghian1, R Shirsavar2, M Erfanifam2, M Sabouhi2, E van der Linden1, H A Stone3, D Bonn4, Mehdi Habibi1. 1. Physics and Physical Chemistry of Foods , Wageningen University , Wageningen 6708 PB , Gelderland , The Netherlands. 2. Department of Physics, Faculty of Science , University of Zanjan , Zanjan 45371-38791 , Zanjan , Iran. 3. Department of Mechanical and Aerospace Engineering , Princeton University , Princeton , New Jersey 08544 , United States. 4. Institute of Physics, van der Waals-Zeeman Institute , University of Amsterdam , Science Park 904 , Amsterdam 1098 XH , North Holland , The Netherlands.
Abstract
We study the spreading of a droplet of surfactant solution on a thin suspended soap film as a function of dynamic surface tension and volume of the droplet. Radial growth of the leading edge (R) shows power-law dependence on time with exponents ranging roughly from 0.1 to 1 for different surface tension differences (Δσ) between the film and the droplet. When the surface tension of the droplet is lower than the surface tension of the film (Δσ > 0), we observe rapid spreading of the droplet with R ≈ tα, where α (0.4 < α < 1) is highly dependent on Δσ. Balance arguments assuming the spreading process is driven by Marangoni stresses versus inertial stresses yield α = 2/3. When the surface tension difference does not favor spreading (Δσ < 0), spreading still occurs but is slow with 0.1 < α < 0.2. This phenomenon could be used for stretching droplets in 2D and modifying thin suspended films.
We study the spreading of a droplet of surfactant solution on a thin suspended soap film as a function of dynamic surface tension and volume of the droplet. Radial growth of the leading edge (R) shows power-law dependence on time with exponents ranging roughly from 0.1 to 1 for different surface tension differences (Δσ) between the film and the droplet. When the surface tension of the droplet is lower than the surface tension of the film (Δσ > 0), we observe rapid spreading of the droplet with R ≈ tα, where α (0.4 < α < 1) is highly dependent on Δσ. Balance arguments assuming the spreading process is driven by Marangoni stresses versus inertial stresses yield α = 2/3. When the surface tension difference does not favor spreading (Δσ < 0), spreading still occurs but is slow with 0.1 < α < 0.2. This phenomenon could be used for stretching droplets in 2D and modifying thin suspended films.
The
dynamics of the spreading of a droplet over a solid or liquid
surface is a phenomenon that has attracted much attention in the past
decades.[1,2] Beside its relevance to many fields of technology,[3−5] studying the dynamics of spreading can provide understanding of
the various forces acting on an interface.It has been shown
that depositing a droplet on a solid surface
will initially result in the formation of a precursor film[6,7] with a molecular thickness,[8] which decreases
the spreading rate of the main droplet. However, Fraaije and Cazabat[9] showed that spreading a droplet on a liquid substrate
exhibits faster dynamics, and an even faster spreading has been observed
in a drop-on-drop geometry.[10] In these
previous studies, the radius of the spreading front increases as a
function of time according to power law. The exponent of this power-law
growth depends on physicochemical properties of both the spreading
liquid and the substrate, for example, surface energy, roughness,
viscosity,[11] miscibility,[12] immiscibility,[13] and the depth
of the substrate for liquid substrates.[9,14,15] When capillary forces are driving the spreading of
a drop on a solid substrate and inertial effects are negligible, a
balance between viscous forces and capillary forces predicts a slow
time evolution of the spreading front, R ≈ t1/10, known as Tanner’s law.[16−18] When two different liquids are brought into contact with each other,
it is often the difference in the surface tensions (Marangoni stress)
that drives the spreading.[19] For spreading
on a thin liquid substrate where the thickness of the film is much
smaller than the radius of the spreading front, the lubrication approximation
applies. Jensen and Grotberg[20] have shown
that the when viscous effects dominate inertial effects, the width
of a planar monolayer strip containing a fixed mass of the surfactant
spreads in time as t1/3, and the radius
of an axisymmetric monolayer droplet evolves as t1/4. Numerical studies[21] also
confirmed the exponent of 1/4 for spreading of a droplet on a rigid
substrate.Although spreading has been studied under various
conditions, spreading
a droplet on a soap film of a few micrometers thickness, suspended
in air, has not been studied yet to the best of our knowledge. Here,
we study a specific case of the droplet spreading on a liquid substrate:
we deposit droplets of surfactant solution on a suspended thin soap
film and document the spreading dynamics as a function of surface
tension and droplet size using high-speed imaging. In contrast to
our experimental system, thin liquid substrates studied in previous
works[15] were at least a few millimeters
thick and in contact with a solid boundary, which causes significant
viscous dissipation because of the no-slip boundary condition.
Experimental Section
To make a horizontal soap film, we used a concentric cylinder made
of stainless steel with inner and outer diameters of 39.5 and 40.0
mm, respectively. The concentric cylinder was dipped into sodium dodecyl
sulfate (SDS) (Merck) solution to produce a horizontal soap film.
Then, a droplet of ammonium lauryl sulfate (ALS) (Fluka) solution
or SDS solution was deposited on the soap film using a capillary tube.
The average volume of deposited droplets was 6 ± 0.5 μL,
although in one series of experiments, we varied the deposited volume
to see how spreading depends on the droplet size. To vary the droplet
size, capillary tubes of different diameters were used; for very small
droplet sizes, the liquid was sprayed very close to the soap film.
SDS solutions with molar concentrations of 0.005 M were used to make
the initial soap films, and SDS solutions with concentrations of 0.005,
0.008, 0.01, and 0.02 M and ALS solutions with concentrations of 0.001,
0.002, 0.003, 0.005, 0.01, 0.04, 0.1, and 0.3 M were used for the
droplets deposited on the soap film. The critical micelle concentrations
(cmc’s) of ALS and SDS solutions at 25 °C are 0.0065[22] and 0.0082 M,[23] respectively.
For the deposited droplets, the concentrations of solutions corresponded
to a range from 0.15 to 46 and 0.6 to 2.4 cmc for ALS and SDS, respectively.
Dynamic surface tensions of solutions were measured via the bubble
pressure method using a Krüss BP50 bubble surface tensiometer.
The thickness of the initial soap film was estimated indirectly from
the speed of a traveling mechanical wave on the soap film using the
equation for wave speed on an elastic sheet under tension (see Supporting Information). After depositing an
ALS droplet, the evolution of the deposited drop on the soap film
was recorded using a high-speed camera (Phantom) at a rate of 4000
frames per second. To illuminate the surface of the soap film, we
used a white light source with the incident and reflected angles of
∼45°. The radius of the leading edge was measured by pixel
counting on the frames.
Results and Discussion
We focus on the spreading of a droplet on a suspended liquid film
of a few micron thickness. This is different from the other work on
Marangoni spreading, where the minimum thickness of liquid substrates[15] was typically about a few millimeters and where
the substrates are in contact with a solid boundary. Therefore, we
expect the viscous dissipation in our system to be much smaller, which
can lead to faster spreading. Furthermore, having a suspended soap
film allows us to confine the liquid–liquid interface solely
to a circular rim.To investigate the position of the droplet
and the soap film during
the course of spreading, we first made qualitative observations. Figure shows an example
of spreading a droplet of 0.3 M ALS on a 0.005 M SDS soap film. In Figure , spreading of a
droplet of ALS 0.4 M on an SDS 0.005 M soap film is shown, in which
a fluorescent dye was added to the SDS solution. The sequence of images
reveals that, after deposition of the droplet, a nonfluorescent circular
film is formed in the center, gradually growing and pushing the fluorescent
SDS soap film toward the edge of the container. The absence of fluorescence
in the spreading film indicates that the ALS droplet is hardly contaminated
by the initial soap film.
Figure 1
Spreading a droplet of 0.3 M ALS on a soap film
of 0.005 M SDS.
Images from left to right show the instant of deposition and subsequent
stretching of the droplet. From left to right, the images are taken
at (a) 5, (b) 12, and (c) 24 ms after deposition.
Figure 2
Spreading
of a droplet containing ALS 0.4 M on a soap film of SDS
0.005 M. Fluorescent dye (blue) is added to the initial SDS soap film.
From left to right: the droplet (black circle) is stretched toward
the edges of the container. Two liquids do not mix during spreading.
From left to right, images are taken at (a) 5, (b) 15, and (c) 30
ms after deposition of the droplet.
Spreading a droplet of 0.3 M ALS on a soap film
of 0.005 M SDS.
Images from left to right show the instant of deposition and subsequent
stretching of the droplet. From left to right, the images are taken
at (a) 5, (b) 12, and (c) 24 ms after deposition.Spreading
of a droplet containing ALS 0.4 M on a soap film of SDS
0.005 M. Fluorescent dye (blue) is added to the initial SDS soap film.
From left to right: the droplet (black circle) is stretched toward
the edges of the container. Two liquids do not mix during spreading.
From left to right, images are taken at (a) 5, (b) 15, and (c) 30
ms after deposition of the droplet.This qualitative observation is confirmed by calculating the ratio
of advection rate to diffusion rate of our system. This ratio, the
Peclet number, is given by Pe = lu/D, where l is the characteristic
length of the system, u is the average flow velocity,
and D is the diffusion coefficient. By considering
the size of the droplet (∼1 mm) as the characteristic length,
the lowest spreading speed (∼0.1 m/s), and diffusion coefficient
of D ≈ 10–9 m2 s–1 (for a typical surfactant in water[24]), we arrive at a Peclet number significantly
greater than 1, confirming that diffusion at the boundaries is negligible
during the course of fast spreading.
Effect
of Surface Tension
Figure shows the dynamic
surface tension for different concentrations of ALS and SDS as a function
of the surface age measured using a bubble pressure surface tensiometer.
To determine the effective surface tension difference (Δσ)
between the droplets and the soap film in our experiments, we use
the time-dependent data of Figure . For the ALS droplet, we consider a characteristic
time scale given by τ = L/u, where u is the average speed of the spreading
film and L is the maximum radius of the spreading
film (∼1 cm). For the SDS substrate, in each experiment there
were at least a few seconds between making the soap film and depositing
the ALS droplet; therefore, the surface tension of the SDS at 5 s
(Figure ) is used
to calculate Δσ. Taking this into account, the surface
tension difference between droplets (ALS) and the soap film (SDS)
covers a wide range, from −19.4 mN/m for ALS 0.001 M to 6.3
mN/m for ALS 0.3 M.
Figure 3
Dynamic surface tension as a function of surface age for
(a) different
concentrations of ALS solutions and 0.005 M SDS (open symbols) and
(b) different concentrations of SDS solutions. Details of the dynamic
surface tension measurements are explained in Supporting Information.
Dynamic surface tension as a function of surface age for
(a) different
concentrations of ALS solutions and 0.005 M SDS (open symbols) and
(b) different concentrations of SDS solutions. Details of the dynamic
surface tension measurements are explained in Supporting Information.In order to identify how surface tension regulates spreading, droplets
with different concentrations of ALS and SDS (different surface tensions)
were deposited on the soap film containing SDS 0.005 M. After the
moment of deposition, the radius (R) of the leading
front was measured as a function of time using high-speed imaging.
The results are shown in Figure a and b for ALS and SDS droplets spreading on an SDS
film, respectively. For the experiments where the radius of the leading
front went beyond 1 cm, a cross-over regime is observed (open symbols
in Figure ), probably
due to boundary effects (the radius of the cylinder is 2 cm). After
discarding the data points at and beyond the cross-over region, a
clear power-law behavior for R versus time is observed
(R ∝ tα).
By fitting power-law functions to the closed-symbol data sets, we
find exponents α ranging from 0.11 to 0.9 for different surface
tension differences (Figure a). Power-law growth of spreading fronts has been observed
in many other spreading phenomena.[9,12,15,25]
Figure 4
Spreading radius versus
time for droplets of ALS (a) and SDS (b)
with different concentrations on a soap film of 0.005 M SDS. Open
symbols show the cross-over regime because of the boundary effects.
Error bars are smaller than the size of the symbols.
Figure 5
(a) Exponent α as a function of surface tension difference
for the deposited droplets with different concentrations of ALS and
SDS on the 0.005 M SDS soap film. α is calculated by a power-law
fit to data in Figure . (b) Average spreading speed calculated by linear fits to the solid
data points in Figure . For the sake of clarity, only several typical error bars for the
velocity and Δσ are shown.
Spreading radius versus
time for droplets of ALS (a) and SDS (b)
with different concentrations on a soap film of 0.005 M SDS. Open
symbols show the cross-over regime because of the boundary effects.
Error bars are smaller than the size of the symbols.(a) Exponent α as a function of surface tension difference
for the deposited droplets with different concentrations of ALS and
SDS on the 0.005 M SDS soap film. α is calculated by a power-law
fit to data in Figure . (b) Average spreading speed calculated by linear fits to the solid
data points in Figure . For the sake of clarity, only several typical error bars for the
velocity and Δσ are shown.For droplets with low concentrations of ALS, where the surface
tension of the droplets is larger than or equal to the surface tension
of the soap film (Δσ ≤ 0), a slow spreading regime
was observed. Here, the surface tension difference is not favorable
for the spreading as the new stretched liquid film has a higher surface
tension than the original film. For deposited droplets of ALS 0.001,
0.002, and 0.003 M, for which Δσ = −19.4, −10.7,
and −3 mN/m, respectively, we obtain power-law exponents of
0.12 ± 0.04, 0.14 ± 0.04, and 0.19 ± 0.04, respectively
(Figure a).Figure shows the
situation when the droplet and soap film are of the same composition
and concentration: SDS 0.005 M. After deposition of the droplet, it
forms a lens and spreads slowly (black squares in Figure b) with a power-law time dependence
and an exponent of 0.15 ± 0.04. This is in contrast with an earlier
experimental work by Aarts et al.[26] for
a three-dimensional configuration where coalescence of two droplets
of the same liquid was studied and a linear dependence of radius on
time (for early times) was measured, although the current work considers
a system with surfactant molecules at the surface and a different
surface geometry.
Figure 6
Droplet of 0.005 M SDS is deposited on the same soap film.
After
deposition, a liquid lens is formed that spreads slowly in time. R as a function of time for this experiment is shown by
black squares in Figure b.
Droplet of 0.005 M SDS is deposited on the same soap film.
After
deposition, a liquid lens is formed that spreads slowly in time. R as a function of time for this experiment is shown by
black squares in Figure b.When a drop of SDS spreads slowly
on a film of the same liquid
(Δσ = 0), there is a clear energetic reason because of
lowering the total surface area. This slow spreading with small exponents
may be explained by the regime of Tanner’s law. In Tanner’s
regime, the driving force of spreading is because of capillary pressure
induced by the curvature at the rim. We note that Tanner’s
law applies for spreading on a rigid boundary. If the surface exhibits
a slip, then, one might expect a slightly different result.[27] We have an almost similar situation for spreading
a droplet on the same liquid film, where large curvatures are observed
at the rim of the liquid lens. The exponent predicted by Tanner’s
law is shown by a horizontal dashed line on Figure a; it is in agreement with the spreading
exponents observed for (Δσ = 0). For droplets with larger
surface tension with respect to the film (Δσ < 0),
we still see slow spreading dynamics with similar exponents (α
< 0.2). At first glance, it seems that the spreading should not
occur because of higher surface tension of the final film. However,
if we consider the very slow dynamics, we can expect that the surface
tension difference between the droplet and the film rapidly decreases
and reaches a uniform surface tension, which is lower than the original
surface tension of the film surface, shortly after the coalescence.
In particular, in the slow spreading regime, the spreading speed is
about 1 mm/s, so that it takes several seconds to spread about 1 cm.
The data in Figure shows that the surface tension decreases considerably over this
time scale, so we anticipate that the surface tension everywhere becomes
uniform in the slow spreading regime. Therefore, even in case of Δσ
< 0, the slow spreading continues with the Tanner-like dynamics,
as discussed above.For the droplets with lower surface tension
than the soap film
(Δσ > 0), the exponent α strongly increases with
increasing Δσ (Figure a). Here, the positive surface tension difference favors
droplet spreading. Because of very fast stretching of the droplet,
we expect a gradient of surface tension from the outside film toward
the center of the spreading droplet.We estimate the average
spreading speed by fitting a linear function
to the data in Figure a and b. The results are shown as a function of Δσ in Figure b, ranging from 0.003
m/s for negative Δσ to about 0.7 m/s for Δσ
= 6.3 mN/m. From the average spreading speed, we determine an experimental
time scale for the spreading. Using this time scale, we obtain the
dynamical surface tensions in Figure .In our system, we observe a range of exponents
from α ≈
0.1 to linear growth of the front radius (α ≈ 1). A similar
range of spreading exponents was reported previously by Rafaï
et al.[28] in a completely different system
consisting of a surfactant solution (Trisiloxan) spread on a solid
substrate, where by increasing the concentration of surfactant, α
increased from 0.2 to 1. This wide range of exponents was associated
with Marangoni effects.
Effect of the Droplet Volume
In order
to study the effect of the volume of the droplets on the spreading
process, droplets with different diameters (D) were
deposited on an SDS 0.005M film. We used ALS 0.04 M droplets for fast
spreading and SDS 0.005 M droplets for slow spreading. Figure shows the radius of spreading
as a function of time for droplets with different volumes deposited
on the soap film. We fit the time-dependencies with a power law and
find that the power-law exponent is roughly the same for all droplet
volumes in each regime (α ≈ 0.9 for fast spreading and
α ≈ 0.11 for the slow regime), as shown in Figure a. Only for the smallest droplet
of ALS, created by spraying, an exponent α = 0.8 ± 0.05
was obtained, which is a little smaller than that of the other experiments
in the fast regime. For the fast spreading experiments, because the
exponent α is close to unity, we can estimate the spreading
speed from the slope of linear fits to the data points. In Figure b, the spreading
speed is reported as a function of the droplet diameter. By increasing
the size of the droplets, the average speed of spreading increases
slowly before reaching a plateau. The capillary length of the liquid
(, where σ = 41.2
mN/m is the dynamic
surface tension of the ALS 0.04 M) determines where the plateau starts.
For the droplets with diameters smaller than the capillary length
of the system, the average velocity grows linearly with the increasing
radius of the droplet (slope 1 in the log–log plot of Figure b). After that, the
speed remains constant at 0.5 ± 0.15 m/s. Lc is indicated in Figure b as a vertical dashed line. According to our results
for the fast spreading, the volume of the spreading droplet does not
affect the spreading exponent but the spreading speed can be influenced
by the droplet volume when the diameter of the droplet is smaller
than the capillary length of the liquid.
Figure 7
Time evolution of spreading
droplets with different volumes on
SDS soap films. (a) Spreading of 0.04 M ALS droplets with different
volumes (0.02, 1.7, 2.2, 3.5, 4.3, 5, and 11.3 μL) on a 0.005
M SDS soap film (fast spreading). (b) Spreading of 0.005 M SDS droplets
with different volumes (0.04,0.23, 2.42, 5.41, and 9.49 μL)
on a 0.005 M SDS soap film (slow spreading). Error bars are smaller
than the size of the symbols.
Figure 8
(a) Exponent
α calculated by a power-law fit to the data
in Figure as a function
of the diameter of the droplets (D). Upper and lower
horizontal dashed lines represent the average value of α for
rapid and slow spreading, respectively. (b) Average spreading speed
as a function of the droplet size with logarithmic scales for fast
(0.04 M ALS droplet on a soap film of 0.005 M SDS) and slow spreading
(0.005 M SDS droplet on 0.005 M SDS soap film). Vertical dashed line
illustrates the capillary length of the system. Blue line with slope
one is shown as a guide to the eye.
Time evolution of spreading
droplets with different volumes on
SDS soap films. (a) Spreading of 0.04 M ALS droplets with different
volumes (0.02, 1.7, 2.2, 3.5, 4.3, 5, and 11.3 μL) on a 0.005
M SDS soap film (fast spreading). (b) Spreading of 0.005 M SDS droplets
with different volumes (0.04,0.23, 2.42, 5.41, and 9.49 μL)
on a 0.005 M SDS soap film (slow spreading). Error bars are smaller
than the size of the symbols.(a) Exponent
α calculated by a power-law fit to the data
in Figure as a function
of the diameter of the droplets (D). Upper and lower
horizontal dashed lines represent the average value of α for
rapid and slow spreading, respectively. (b) Average spreading speed
as a function of the droplet size with logarithmic scales for fast
(0.04 M ALS droplet on a soap film of 0.005 M SDS) and slow spreading
(0.005 M SDS droplet on 0.005 M SDS soap film). Vertical dashed line
illustrates the capillary length of the system. Blue line with slope
one is shown as a guide to the eye.For the slow spreading, although the power law exponents are significantly
smaller than unity (α ≈ 0.11), we still are able to linearly
fit the time evolution data of Figure and estimate a spreading speed. The results are shown
by square symbols in Figure b and indicate that the average spreading speed remains approximately
constant with an increase in the volume of the droplet.
Scaling Arguments
Bringing two miscible
liquids into contact with each other will generate a gradient of surface
tension and induces a Marangoni-driven flow. The driving force in
such flows is the tangential stress associated with gradients of surface
tension at the interface which should be balanced by viscous stresses.
Based on scaling the Marangoni stresses for viscosity-limited spreading
on a deep bath (distance as a function of time), a spreading exponent
of 3/4 can be derived.[29−31]The Reynolds number (Re =
ρuh/η, with h as the
film thickness) of the system in the rapid spreading regime (Δσ
> 0) is always larger than one, indicating that the inertial terms
should be considered in the balance of stresses. In the first approximation,
we ignore the viscous dissipation and balance the surface tension
gradient by the inertial term: |∇σ| ≈ ρu2, where u ≈ Ṙ with a dot denoting the time derivative. The surface
tension gradient term scales with the dynamic surface tension difference
divided by the radius of the spreading front (|∇σ| ≈
Δσ/R), so we arrive at Δσ/R ≈ ρṘ2 and
hence, R ≈ (Δσ/ρ)1/3t2/3. The exponent 2/3 is within the
range of exponents we observed for the fast spreading experiments,
as shown in Figure a by a dashed line marked by “inertial”.
Conclusions
Droplets will spread when deposited on
a soap film. Spreading evolves
as a function of time in a power-law manner with exponents ranging
from about 0.1 to about 1 depending on the dynamic surface tension
differences between the film and the droplet. Small spreading exponents
with slow spreading dynamics occur when the surface tension of the
droplet is greater than or equal to the surface tension of the film
so that the tension difference does not favor spreading. This slow
spreading regime resembles Tanner’s spreading on a solid substrate
with exponent 0.1, as driven by capillary pressure at the rim. By
increasing the surface tension difference, the spreading exponent
increases toward 1. In this regime, the spreading speed increases
significantly with surface tension difference. Here, the driving mechanism
is the Marangoni stress. By scaling the Marangoni stress with inertial
terms, we can derive spreading exponents of α = 2/3. Although
this exponent is in the range of experimentally observed exponents
for the fast spreading regime, however it does not explain the very
fast spreading dynamics with exponents about one. This might be because
of the fact that in our scaling argument, we have simplified all the
kinetic effects and surfactant absorption dynamics in a surface tension
gradient. We also find that the spreading speed depends on the size
of the droplet in the fast spreading regime. The spreading speed depends
linearly on the droplet size for drops smaller than the capillary
length of the liquid and remains constant for larger droplet sizes;
however, the spreading exponent is not influenced by the droplet size.
Our experimental results may help in understanding of the various
forces acting on a droplet deposited on a liquid interface. Using
a more detailed theoretical model in which the surfactant absorption
kinetics is considered may elucidate all our experimental observations.
In addition, our methodology can be used for stretching complex liquid
droplets in 2D and therefore may have applications in modifying thin
suspended films. For example, using this approach, the effect of elongational
viscosity on stretching a complex liquid droplet in 2D can be studied.