Contraction of Surfactant-Laden Pores - Langmuir (ACS Publications)

Mar 14, 2018 - (25) Finally, the nonlinear system was solved simultaneously for velocity, pressure, surfactant concentration, and location of the free...
0 downloads 0 Views 577KB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

New Concepts at the Interface: Novel Viewpoints and Interpretations, Theory and Computations

Contraction of surfactant-laden pores Jiakai Lu, Diego Martín Martín Campana, and Carlos M Corvalan Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.8b00633 • Publication Date (Web): 14 Mar 2018 Downloaded from http://pubs.acs.org on March 22, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Contraction of surfactant-laden pores Jiakai Lu,† Diego M. Campana,‡ and Carlos M. Corvalan∗,† †Transport Phenomena Laboratory, Department of Food Science, Purdue University, West Lafayette, IN 47907, USA ‡Centro de Investigaciones y Transferencia de Entre R´ıos, CONICET-UNER. Facultad de Ingenier´ıa, Universidad Nacional de Entre R´ıos, Ruta 11 (Km. 10), 3100 (Oro Verde), Entre R´ıos, Rep´ ublica Argentina E-mail: [email protected]

Abstract The contraction of surfactant-laden pores at the microscale has implications for natural and technological processes ranging from the collapse of channels in lipid membranes to the stability of foams in the food processing industry. Despite its prevalence, our understanding of the mechanisms of pore contraction in the presence of surfactants remains unclear. These mechanisms have been challenging to study experimentally given the small length scale near the singularity and simulations capable of accurately characterizing the pore dynamics may help enhance our understanding of the process. Here, we use high-fidelity numerical simulations to gain insight into the fluid dynamics and interfacial phenomena underlying the contraction of viscous pores in the presence of an insoluble surfactant. Results show that surfactants accumulate on the advancing front of a collapsing pore due to the uneven deformation of the pore interface. Because of this accumulation, even a small amount of surfactant plays a major role in the way in which a collapsing pore approaches the singularity. keywords: pore, nanopore, surfactant, free-surface, DNA sequencing

1

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction The contraction of surfactant-laden pores at the micron and submicron scale is a critical phenomenon in nature, as in the collapse of pores in lipid membranes 1,2 , and central to technological systems ranging from the processing of foams in the food industry to the fabrication of pore-based biosensors for analysis of nucleic acids and proteins 3–8 . Currently, there is a good understanding of the fluid dynamics of both expanding 9–11 and contracting pores 12–14 for fluids with clean interface. But a recent theoretical and experimental study by Petit et al. 15 has highlighted the major role that the presence of surface-active species may play modulating the pore dynamics. Petit et al. 15 observed that opening pores nucleated in a soap film with strong surface elasticity expand at a velocity that is significantly smaller than the Taylor-Culick velocity expected from the balance between inertia and the equilibrium surface tension of the film 9,10 . This deviation from the theoretical Taylor-Culick velocity was attributed to local gradients of surface tension resulting from the transport of the surfaceactive species on the soap film. The experiments by Petit et al. 15 , as well as other previous works 16,17 , provide important insights into the influence of surfactants on the dynamics of opening pores. But the role of surfactants on the fluid dynamics of contracting pores is still largely unknown, in part due to the small length scale of the process near the singularity. In particular, how surface-active species disperse through the deforming pore interface during the contraction of small pores, and how the resulting changes in the local surface tension affect the velocity of contraction, remains an open question. Here, we address this question for the case of viscous pores covered with an small amount of an insoluble surfactant. We characterize the fluid dynamics and interfacial phenomena underlying the contraction of the pores using high-fidelity simulations that simultaneously solve the full Navier-Stokes system that governs the free-surface flow and the convection-diffusion equation that governs the surfactant transport. Results reveal that, contrary to viscous pores with clean interface, which contract at constant speed 14 , surfactant-laden pores contract at progressively decreasing speed. This new behavior, the 2

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

simulations show, results from the continuous accumulation of the surface-active species on the advancing front of the pore. Because of this accumulation, even a tiny amount of the surfactant may play a major role in the way in which a collapsing pore approaches the singularity.

Free-surface Model and Methods A sufficiently small pore nucleated in a liquid sheet will contract, and eventually collapse, driven by capillary forces associated with its small radius of curvature 9 . Here, we consider the capillary contraction of a small circular pore nucleated in a viscous fluid sheet of thickness H. The fluid is a Newtonian liquid of constant density ρ and large viscosity µ, uniformly covered with an insoluble surfactant of concentration γˆ0 (Fig. 1). Recently, Savva and Bush 11 investigated the dynamics of opening pores using a onedimensional lubrication approximation that exploits the small thickness of the liquid sheet. This approximation has been shown to accurately predict pore expansion 13,14 but it cannot predict the evolution of closing pores because the thickness-to-radius ratio grows unbounded as the pore collapses. Therefore, to study the free-surface dynamics of collapsing pores, we have developed high-fidelity simulations that overcome the limitations of the lubrication model by solving the full Navier-Stokes governing equations. The model is described below in dimensionless form using H as lengthscale, γˆ0 as concentration scale, and V ≡ σ ˆ0 /µ, where σ ˆ0 is the equilibrium surface tension corresponding to the surfactant concentration γˆ0 , as velocity scale. The contraction of the pore is analyzed by solving the full axisymmetric Navier-Stokes equations Re(∂v/∂t + v · ∇v) = ∇ · T,

(1)

∇ · v = 0,

(2)

where Re ≡ ρV H/µ is the Reynolds number and where the influence of gravity is considered

3

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 24

negligible. The stress tensor T = −pI + [(∇v) + (∇v)T ], where p is pressure and v the fluid velocity. The fluids does not cross the free interface, which is ensured by imposing the kinematic boundary condition along the pore surface

n · (v − vs ) = 0,

(3)

where n is the unit normal vector to the free surface and vs is the velocity of the points on the interface. The driving surface tension σ is coupled to the fluid flow through the traction condition along the pore interface T · n = σ(2H)n + ∇s σ,

(4)

where H in the capillary stress term is the mean surface curvature, and ∇s ≡ (I − nn) · ∇ in the Marangoni stress term is the surface gradient operator 18 . The time evolution of the surfactant concentration γ along the pore interface is tracked by the transport equation

∂γ/∂t + γ(vs · n)(∇s · n) + ∇s · (γvs · t)t − P e−1 ∇2s γ = 0,

(5)

which accounts for surfactant transport due to surface dilatation (second term), convection (third term), and diffusion (fourth term) 19,20 . For realistic values of the surfactant diffusion coefficient D, the Peclet number P e ≡ HV /D  1. We use P e = 103 in our simulations unless stated; the results indicate that this value is sufficiently large that further increases do not significantly affect the pore dynamics for the ranges of parameters studied in this work. Following Hansen et al. 21 , Campana and Saita 22 and Dravid et al. 23 , we assume that for dilute surfactants the equation of state is linear

σ = 1 − β(γ − 1),

(6)

so that surface activity is represented by a single dimensionless parameter, the elasticity 4

ACS Paragon Plus Environment

Page 5 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

number β ≡ −(ˆ γ0 /ˆ σ0 )dˆ σ /dˆ γ = −dσ/dγ. The magnitude of the dimensionless interfacial tension for uniformly distributed surfactant at the initial concentration γˆ = γˆ0 is thus σ = 1 (or σ ˆ=σ ˆ0 in dimensional terms). A small elasticity number (β  1) can be understood as corresponding to a weak surfactant (dˆ σ /dˆ γ  1) or to a dilute surfactant (ˆ γ0  1). Here we assume dilute surfactant concentration. The pore is considered axisymmetric around the z-axis, initially at rest, and symmetry boundary conditions are applied on the plane of symmetry z = 0. The numerical scheme has been described previously in Lu et al. 14 , Lu and Corvalan 20 , Dravid et al. 23 and so only a brief summary is given. Following Lu and Corvalan 13 , the set of governing equations was discretized in space using the finite-element method, along with an arbitrary Lagrangian-Eulerian scheme in which the location of the free interface is traced using the method of spines introduced by Kistler and Scriven 24 . The time derivatives were discretized using the finite-difference method, with an Adam-Bashforth predictor and trapezoidal rule corrector in which the time steps were adaptively calculated using first-order continuation 25 . Finally, the non-linear system was solved simultaneously for velocity, pressure, surfactant concentration and location of the free interface using a full Newton’s method in which the entries of the Jacobian were calculated analytically as described in Kistler and Scriven 24 to enhance convergence. We have repeatedly verified this algorithm by contrasting its predictions against available analytical solutions and experimental data during the analysis of both viscous and inertial nanopores 13,14 , as well as surfactant-laden drops 20 , and fibers 26 .

Results and Discussion Below a threshold size, a pore punctuated or nucleated in a fuid sheet will contract driven by capillarity to decrease surface energy. The threshold size for capillary contraction is on the order the thickness of the sheet 13,27,28 As recently reviewed by Neel and Villermaux, 29

5

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

these small pores can be nucleated in liquid films using a variety of techniques including penetrating the liquid film with a sharp object, a projectile, 28 an air jet, 30 a spark, a laser beam 31 and an electron beam; 27 or by introducing bubbles, 7 oil droplets, 32 or hydrophobic particles 33 in the plane of the film. Studies of pore collapse have been driven by several areas relevant both in nature 1,2 and in industry. 3–8 In a pioneering work on the dynamics of small pores, Taylor and Michael 12 studied the collapse of pores nucleated in soap films, motivated by the stability of bubbles, foams and liquid sheets. And much of the recent impetus behind studies on the contraction of viscous pores – including the experimental studies by Storm et al., 27,34 Wu et al. 6 — have been motivated by the development of nanopore based sensors for sequencing necleic acids and analyzing epigenetic modifications. 3 These experiments have revealed that in the absence of surfactants viscous nanopores contract at constant speed when approaching collapse. As recently confirmed by computations in Lu et al. 14 , this constant terminal speed results from an eventual equilibrium between viscous and capillary forces, and thus corresponds to U =σ ˆ0 /(2µ) for pores with large viscosity µ and constant surface tension σ ˆ0 . Extending the work in Lu et al. 14 , here we present results from high-fidelity simulations of pore contraction that take into account the effect of surfactants. Particular attention is given to the influence of contaminants on the speed of contraction.

Influence of surfactants on the speed of contraction To gain preliminary insight into how the presence of contaminants affects the pore dynamics, we show in Figure 2 the computed pore velocity vm as a function of minimum pore radius rm for two viscous pores with elasticity numbers β = 5 × 10−4 (red line) and 1 × 10−3 (blue line). The pores are nucleated in a liquid sheet with small Reynolds number Re = 0.5. As a benchmark, and to compare our results to those expected for pores with constant surface tension, the figure also includes the computed pore speed for an otherwise identical pore with β = 0 (black line). The figure follows the contraction of the pores for about three

6

ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

decades, from the initial dimensionless radius rm = 0.1 to a final radius rm ≈ 3 × 10−4 , or approximately 300 nm for a 100 µm initial pore radius. Note that length scales that are absent from our study would eventually become relevant as the pore radius rm → 0. For example, thermal fluctuations, which can provoke instabilities on the interface, become p relevant for length scales in the order of lT = kb T /σ, where kb T is the thermal energy for unit area. 35 The thermal length scale lT is typically less than 1 nm at room temperature. Despite the small Marangoni numbers, the results show a significant qualitative difference between the pores with β = 0 and those with β 6= 0. For the pore with β = 0, the interfacial tension remains constant on the surface of the pore until collapse (Eq. 6). Accordingly, the pore eventually attains a constant terminal speed, as observed in experiments for pores with constant surface tension during the late stage of contraction 6,27 . Furthermore, the terminal speed corresponds to the theoretical viscous-capillary velocity U = σ ˆ0 /(2µ), which translates to u = 1/2 in the dimensionless scale of Figure 2 (dashed line). For the pores with β 6= 0, the interfacial tension changes dynamically on the surface of the pores due to surfactant diffusion, surfactant convection or area dilatation (Eqs. 5 and 6). As a result, the pores do not attain a constant terminal speed; rather, they contract at continuously decreasing speed after the initial transient. By the time the pore with β = 1 × 10−3 reaches rm = 3 × 10−4 , the speed of contraction is about 20% slower than the viscous-capillary velocity u based on the equilibrium surface tension. The interfacial stress balance helps understand the qualitative difference between the two dynamics and leads, in the next section, to a simple scaling of the decreasing pore speed with pore size. Indeed, as discussed in Figure 2, viscous pores with clean interface contract with constant speed σ ˆ0 /(2µ) due to the eventual equilibrium between viscous and capillary stresses. But this equilibrium is affected by the presence of surfactants in two ways. First, surfactants modify the capillary stress by lowering the dynamic surface tension and, second, they introduce Marangoni stresses into the stress balance (Eq. 4). However, at low elasticity numbers, Marangoni stresses are not sufficient to displace the equilibrium (see below), and

7

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 24

thus the viscous-capillary equilibrium still applies provided that the constant surface tension σ ˆ0 is replaced by the dynamic surface tension. Therefore, for surfactants with low elasticity numbers the velocity of contraction can be estimated as σ ˆm (t)/(2µ), which in dimensionless form becomes vm (t) ≈ σm (t)/2,

(7)

where σm = σ ˆm /ˆ σ0 is the instantaneous surface tension at the tip of the pore. As shown in Figure 3, the computed pore velocities agree reasonably well with this estimation for pores with low elasticity numbers and a wide range of Peclet numbers 10 ≤ P e ≤ 103 .

Evolution of the dynamic surface tension Analyzing the mechanisms of surfactant transport, and the corresponding evolution of the surface tension σm (t) is important to understanding the new pore dynamics. To illustrate these mechanisms, Figure 4 shows the spatial and temporal evolution of the interfacial tension for the pore of Figure 2 with β = 1 × 10−3 . Results show a highly localized decrease of the surface tension within a narrow band on the advancing pore front. On this narrow band, the surface tension exhibits a minimum at the tip of the pore (z = 0), rises sharply as we move away from the tip, and finally levels out approaching the dimensionless initial value σ = 1 on the surface of the liquid sheet (Fig. 4a). Critically, results also describe how, as the pore contracts, the surface tension at the front of the pore decreases, reducing the capillary stress that drives the collapse, so that the pore velocity diminishes (Fig. 4b). As may be expected, the local decrease of the surface tension is caused by accumulation of surfactant around the tip of the pore. The appropriate scaling for the surfactant accumulation is readily estimated. As a pore collapses, the surface area of the pore front contracts, and this contraction is most severe at the tip where the local surface area decreases proportional to the pore circumference. As a result, the surfactant concentration γm −1 at the tip of the pore increases proportional to the inverse pore radius; that is, γm ∼ rm .

This simple scaling, along with the results in Figure 2, reveals a significant feature of the 8

ACS Paragon Plus Environment

Page 9 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

pore dynamics: even a modest amount of surfactant is eventually capable of unexpectedly large effects due to surfactant crowding. Results in Figure 5 show that as the pore of Figure −1 (dashed line) agrees well with the 2 with β = 1 × 10−3 contracts, the scaling γm ∼ rm

simulations (black symbols). However, because the surfactant concentration must remain finite as rm → 0, the scaling is necessarily an intermediate one and must yield to a different dynamics as the pore moves even closer to the singularity. Moreover, because diffusion mitigates the accumulation process, the surfactant concentration γm should also depend on the Peclet number. Indeed, comparison of the simulations for Peclet numbers P e = 10 (blue), 100 (red) and 1000 (black) reveals that there is a transition between the dilatation and the diffusion dominated regimes of surfactant transport. Results in Figure 5 show that as the Peclet number decreases, diffusion effects eventually becomes comparable to area-dilatation effects, and ultimately diffusion dominate the dynamics (P e ≈ 10). The foregoing results suggest that for realistic values of the Peclet number, for which surfactant transport is dominated by area dilatation, the normalized pore velocity s ≡ (u − vm )/u can be estimated in function of the elasticity number and instantaneous pore size as: −1 (s − β)/β ∼ rm .

(8)

This results from using the surface equation of state to substitute σm in Equation 7, and then introducing the scaling of γm with the inverse pore radius. Taking into account that this relationship applies to contaminants with low elasticity numbers, the scaling can be further simplified to −1 s/β ∼ rm .

(9)

This provides a simple scaling of the pore speed with pore size and surface activity. As illustrated in Figure 6, the scaling (dashed line) agrees well with results from the full numerical simulations (symbols) at the later times of contraction.

9

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Influence of interfacial stresses on the interfacial velocity field A close examination of the flow field developed during the contraction of the pores reveals that surfactants may also play an important role on the interfacial velocity profiles. In this subsection, we briefly discuss the separate influence of the capillary and Marangoni stresses on the interfacial velocity field near collapse. In the absence of surfactants the maximum interfacial velocity always occurs at the tip of the pore where the curvature is largest. The reason for this is that when the interfacial tension is constant the distribution of the capillary forces driving contraction is defined solely by the local curvature. In the presence of surfactants, however, the distribution of the capillary forces (n · T · n) = σ(∇ · n) at the interface is defined not only by the local curvature ∇ · n but also by the distribution of the local surface tension σ(x; t). To illustrate the effect of the evolving surface tension, we compare in Figure 7 the interfacial velocity profiles for the pores of Figure 6 with β = 1 × 10−4 , β = 5 × 10−4 , and β = 1 × 10−3 at an intermediate time and at the imminence of collapse (minimum pore radii rm = 0.095 and 3 × 10−4 respectively). Results show that because of the progressive decrease of the surface tension at the tip of the pore there is a qualitative change in the flow pattern when approaching collapse. The location of the maximum interfacial velocity switches from the tip of the pore at z = 0 (black lines) to two peaks about the midplane (blue lines). The rapid contraction of the surface area not only lowers the local surface tension around the tip of the pore but also creates steep surface tension gradients, particularly at high Peclet numbers (see Fig. 4). In turn, surface tension gradients induce Marangoni stresses, which act tangential to the interface (Eq. 4). These stresses cannot directly affect the interfacial shape but could modulate the velocity field near the interface through the tangential stress condition (n · T · t) = t · ∇s σ. Figure 8, illustrates the influence of the Marangoni stresses on the interfacial tangential velocity for the pores with the lowest (β = 1 × 10−4 ) and highest (β = 1 × 10−3 ) elasticity numbers in Figure 7. The Marangoni stresses are zero at the tip of the pore because of symmetry, then grow rapidly as we move away from the tip to reach 10

ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

a maximum value and finally decrease as we approach the surface of the liquid sheet where the surfactant concentrations remains essentially constant (Fig. 8a). Despite the different values of the Marangoni stresses, they have a negligible influence on the tangential interfacial velocity for the small elasticity numbers used in this work (Fig. 8b).

Conclusion In conclusion, results from our numerical study reveal that surfactant-laden pores collapse at progressively decreasing speed due to surfactant accumulation on the advancing front of the pore. As a consequence of this accumulation even a small amount of contaminant can drastically reduce the speed of contraction. These findings provide new mechanistic insights into the flow physics of surfactant-laden pores at low elasticity numbers, and suggest interesting directions for future studies. For example, considering the strong influence that modest concentrations of surfactants have on the speed of contraction, it would be of interest to investigate the use of small pores as analytical tools for interfacial rheological characterization of dilute surfactants. In addition, although our direct numerical simulations solve the full Navier-Stokes system that governs the free-surface flow along with the full convection-diffusion equation that governs the interfacial surfactant transport, the results are still limited by simplifying constitutive assumptions, including negligible surfactant solubility and linear surface equation of state. We expect that our results would motivate the further development of simulations that incorporate soluble surfactants. Due to the surfactant accumulation at the pore front, it is conceivable that solubility may play an important role in the pore dynamics. Insoluble surfactants cannot easily leave the concentrated front region leading to a large drops of the interfacial tension; but surfactant solubility might mitigate this effect by enabling transfer of surfactant from the concentrated region to the bulk liquid 22 . Similarly, although contamination of the pore interface by a surfactant monolayer is a common situation, a more

11

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

sophisticated surface equation of state would enable the study of high surfactant concentrations with micelles formation. In these conditions, the surface tension in the concetrated region may eventually take a constant equilibrium value as the pore contracts, potentially stabilizing the pore velocity.

Acknowledgement This work was partially supported by the USDA National Institute of Food and Agriculture, Hatch projects 1008409 and 199888. D. C. would like to thank the following institutions for financial support: Consejo Nacional de Investigaciones Cient´ıficas y T´ecnicas of Argentina (CONICET) and Universidad Nacional de Entre R´ıos (UNER).

References (1) Chanturiya, A.; Chernomordik, L. V.; Zimmerberg, J. Flickering fusion pores comparable with initial exocytotic pores occur in protein-free phospholipid bilayers. Proceedings of the National Academy of Sciences 1997, 94, 14423–14428. (2) Photos, P. J.; Bermudez, H.; Aranda-Espinoza, H.; Shillcock, J.; Discher, D. E. Nuclear pores and membrane holes: generic models for confined chains and entropic barriers in pore stabilization. Soft Matter 2007, 3, 364–371. (3) Dekker, C. Solid-state nanopores. Nature Nanotechnology 2007, 2, 209–215. (4) Iqbal, S. M.; Akin, D.; Bashir, R. Solid-state nanopore channels with DNA selectivity. Nature nanotechnology 2007, 2, 243–248. (5) Iqbal, S. M.; Bashir, R. Nanopores: sensing and fundamental biological interactions; Springer Science & Business Media, 2011.

12

ACS Paragon Plus Environment

Page 12 of 24

Page 13 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(6) Wu, S.; Park, S. R.; Ling, X. S. Lithography-free formation of nanopores in plastic membranes using laser heating. Nano letters 2006, 6, 2571–2576. (7) Lhuissier, H.; Villermaux, E. Effervescentatomization in two dimensions. Journal of Fluid Mechanics 2013, 714, 361–392. (8) Timounay, Y.; Lorenceau, E.; Rouyer, F. Opening and retraction of particulate soap films. EPL (Europhysics Letters) 2015, 111, 26001. (9) Taylor, G. I. The Dynamics of Thin Sheets of Fluid. III. Disintegration of Fluid Sheets. Proceedings of the Royal Society of London A 1959, 253, 313–321. (10) Culick, F. E. C. Comments on ruptured soap film. Journal of Applied Physics 1960, 31, 1128–1129. (11) Savva, N.; Bush, J. W. M. Viscous sheet retraction. Journal of Fluid Mechanics 2009, 626, 211 – 240. (12) Taylor, G. I.; Michael, D. H. On making holes in a sheet of fluid. Journal of Fluid Mechanics 1973, 4, 625. (13) Lu, J.; Corvalan, C. M. Free-surface dynamics of small pores. Chemical Engineering Science 2015, 132, 93–98. (14) Lu, J.; Yu, J.; Corvalan, C. M. Universal scaling law for the collapse of viscous nanopores. Langmuir 2015, 31, 8618–8622. (15) Petit, P.; Le Merrer, M.; Biance, A.-L. Holes and cracks in rigid foam films. Journal of Fluid Mechanics 2015, 774, R3. (16) McEntee, W. R.; Mysels, K. J. Bursting of soap films. I. An experimental study. The Journal of Physical Chemistry 1969, 73, 3018–3028.

13

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(17) Frankel, S.; Mysels, K. J. Bursting of soap films. II. Theoretical considerations. The Journal of Physical Chemistry 1969, 73, 3028–3038. (18) Slattery, J. C. Interfacial Transport Phenomena; Springer, New York, 1990. (19) Wong, H.; Rumschitzki, D.; Maldarelli, C. On the surfactant mass balance at a deforming fluid interface. Physics of Fluids (1994-present) 1996, 8, 3203–3204. (20) Lu, J.; Corvalan, C. M. Coalescence of viscous drops with surfactants. Chemical Engineering Science 2012, 78, 9–13. (21) Hansen, S.; Peters, G.; Meijer, H. The effect of surfactant on the stability of a fluid filament embedded in a viscous fluid. Journal of Fluid Mechanics 1999, 382, 331–349. (22) Campana, D.; Saita, F. Numerical analysis of the Rayleigh instability in capillary tubes: The influence of surfactant solubility. Physics of Fluids 2006, 18, 022104. (23) Dravid, V.; Songsermpong, S.; Xue, Z.; Corvalan, C. M.; Sojka, P. E. Two-dimensional modeling of the effects of insoluble surfactant on the breakup of a liquid filament. Chemical engineering science 2006, 61, 3577–3585. (24) Kistler, S. F.; Scriven, L. E. Coating Flow: Computational Analysis of Polymer Processing; Applied Science Publishers, New York, 1983. (25) Corvalan, C. M.; Saita, F. A. Automatic stepsize control in continuation procedures. Computers & Chemical Engineering 1991, 15, 729–739. (26) Campana, D. M.; Ubal, S.; Giavedoni, M. D.; Saita, F. A. Dip Coating of Fibers in the Visco-inertial Regime: Numerical Analysis. Industrial & Engineering Chemistry Research 2013, 52, 12646–12653. (27) Storm, A. J.; Chen, J. H.; Ling, X. S.; Zandbergen, H. W.; Dekker, C. Electron-beaminduced deformations of SiO2 nanostructures. Journal of Applied Physics 2005, 98 . 14

ACS Paragon Plus Environment

Page 14 of 24

Page 15 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(28) Courbin, L.; Stone, H. A. Impact, puncturing, and the self-healing of soap films. Physics of Fluids 2006, 18, 091105. (29) N´eel, B.; Villermaux, E. The spontaneous puncture of thick liquid films. Journal of Fluid Mechanics 2018, 838, 192–221. (30) Lhuissier, H.; Brunet, P.; Dorbolo, S. Blowing a liquid curtain. Journal of Fluid Mechanics 2016, 795, 784–807. (31) Wedershoven, H.; Berendsen, C.; Zeegers, J.; Darhuber, A. Infrared-laser-induced thermocapillary deformation and destabilization of thin liquid films on moving substrates. Physical Review Applied 2015, 3, 024005. (32) Dombrowski, N.; Fraser, R. A photographic investigation into the disintegration of liquid sheets. Phil. Trans. R. Soc. Lond. A 1954, 247, 101–130. (33) Garrett, P. Defoaming: theory and industrial applications; CRC Press, 1992; Vol. 45. (34) Storm, A. J.; Chen, J. H.; Ling, X. S.; Zandbergen, H. W.; Dekker, C. Fabrication of solid-state nanopores with single-nanometre precision. Nature Materials 2003, 2, 537–540. (35) Eggers, J. Dynamics of liquid nanojets. Physical Review Letters 2002, 89, 084502.

15

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 24

Figures

z ρ, μ r

H

σ0 Figure 1: Pore sketch. Small circular pore nucleated in a Newtonian liquid sheet with density ρ, viscosity µ, and thickness H. The initial surface tension σ0 is uniform, and corresponds to the initial surfactant concentration γ0 .

16

ACS Paragon Plus Environment

Page 17 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

0.7 0.5

0.3

-3

-2

10

10

-1

10

Figure 2: Effect of surfactant on the contraction of a small viscous pore. Evolution of the velocity of contraction vm with the minimum pore radius rm for pores with β = 0 (black line), β = 5 × 10−4 (red line), and β = 1 × 10−3 (blue line). The pore with β = 0 eventually approaches the theoretical viscous-capillary velocity U = σ0 /(2µ) (dashed line).

17

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 24

0.5

0.4

0.3 0.3

0.4

0.5

Figure 3: Effect of the local surface tension on the velocity of contraction. Pore velocity varies approximately as vm = σm /2 (solid line). The symbols indicate results from simulations with β = 5 × 10−4 (green), 1 × 10−3 (blue), 3 × 10−3 (red) and 5 × 10−3 (black), and Peclet numbers Pe = 10 (left triangle), 1 × 102 (right triangle) and 1 × 103 (up triangle).

18

ACS Paragon Plus Environment

Page 19 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 4: Evolution of the dynamic surface tension. Distribution of the instantaneous surface tension σ on the surface of the pore of Figure 2 with β = 1×10−3 . (a) Cross sectional fluid velocity and radial distribution of the surface tension for pore radius rm = 3.4 × 10−4 . (b) Axial distribution of the surface tension for pore radii rm = 2.2 × 10−3 , 6 × 10−4 , and 3.4 × 10−4 (from top to bottom).

19

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

2

Page 20 of 24

-1

10

101

0

10

-3

-2

10

10

-1

10

Figure 5: Scaling of surfactant concentration with pore radius. Surfactant concen−1 tration at the pore tip eventually follows the scaling γm ∼ rm for moderate and large Peclet −3 numbers. Here, β = 1 × 10 and P e = 10 (blue), 100 (red), 1 × 103 (black), and 1 × 104 (black). The results for P e = 1 × 103 and 1 × 104 overlap.

20

ACS Paragon Plus Environment

Page 21 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

1

100

10

1000

2000 3000

Figure 6: Scaling of pore velocity with pore radius. Scaling of normalized velocity s with inverse pore radius for pores with β = 5 × 10−4 (green), 1 × 10−3 (red), 3 × 10−3 (blue), and 5 × 10−3 (black).

21

ACS Paragon Plus Environment

ACS Paragon Plus Environment

22

0 -0.5-0.25 0 0.25 0.5

0.3

0.6

Figure 7: Influence of elasticity number on interfacial velocity profiles. Interfacial radial velocity for contracting pores with β = 1 × 10−4 (left), 5 × 10−4 (middle), and 1 × 10−3 (right). Results correspond to pore radius rm = 9.5 × 10−2 (black) and 3 × 10−4 (blue).

0 -0.5-0.25 0 0.25 0.5

0.3

0.3

0 -0.5-0.25 0 0.25 0.5

0.6

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

0.6

Langmuir Page 22 of 24

Page 23 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

10 5 0 0.4 0.2 00

0.05

0.1

Figure 8: Marangoni stresses and tangential interfacial flow. (a) Axial variation of tangential Marangoni stress on the surface of the pores of Figure 7 with β = 1 × 10−4 (blue) and 1 × 10−3 (black) at pore radius rm = 3 × 10−4 . (b) The corresponding tangential interfacial velocities vt overlap. 23 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For Table of Contents Only Contraction of surfactant-laden pores. Jiakai Lu, Diego M. Campana, and Carlos M. Corvalan

24

ACS Paragon Plus Environment

Page 24 of 24