volume 13 number 35 21 september 2017 pages 5811–6022 …
TRANSCRIPT
ISSN 1744-6848
PAPER John E. Dolbow, Eliot Fried et al. Influence of surface tension in the surfactant-driven fracture of closely-packed particulate monolayers
Volume 13 Number 35 21 September 2017 Pages 5811–6022
5832 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
Cite this: SoftMatter, 2017,
Influence of surface tension in the surfactant-driven fracture of closely-packed particulate monolayers
Christian Peco, a Wei Chen,b Yingjie Liu, a M. M. Bandi, c
John E. Dolbow *a and Eliot Fried *b
A phase-field model is used to capture the surfactant-driven formation of fracture patterns in particulate
monolayers. The model is intended for the regime of closely-packed systems in which the mechanical
response of the monolayer can be approximated as that of a linearly elastic solid. The model
approximates the loss in tensile strength of the monolayer with increasing surfactant concentration
through the evolution of a damage field. Initial-boundary value problems are constructed and spatially
discretized with finite element approximations to the displacement and surfactant damage fields.
A comparison between model-based simulations and existing experimental observations indicates a
qualitative match in both the fracture patterns and temporal scaling of the fracture process. The
importance of surface tension differences is quantified by means of a dimensionless parameter,
revealing thresholds that separate different regimes of fracture. These findings are supported by newly
performed experiments that validate the model and demonstrate the strong sensitivity of the fracture
pattern to differences in surface tension.
1 Introduction
When a densely packed monolayer of hydrophobic particles is placed on the surface of a liquid, the particles interact through capillary bridges,1 leading to the formation of particle rafts. The macroscopic properties of these rafts reflect an interplay between fluid and solid mechanics,2–4 giving rise to novel physics. This interplay is relevant to a wide range of applica- tions, from the synthesis of ‘‘liquid marbles’’5 to the design of drug delivery systems6 to the stabilization of drops.7
The interest in particle rafts has driven researchers to investigate their mechanical properties.4,8 It is now known that densely packed monolayers exhibit a two-dimensional linearly elastic solid response. The interaction of the particles through the capillary bridges allows the monolayer to support both tension and compression. Moreover, the mechanical properties of such monolayers depend proportionally on the surface tension of the liquid layer. The introduction of a controlled amount of surfactant generates a surface tension gradient,
producing Marangoni forces9 and causing the surfactant to spread, fracturing the monolayer. Studies of surfactant-driven fracture have examined the role of viscosity and the initial packing fraction on the evolution of cracks in closely and loosely-packed systems, respectively.10,11 Surprisingly, the potential importance of differences in the surface tension of the surfactant and the underlying liquid has not been explored. The magni- tude of this difference is interesting because it determines the Marangoni force exerted on the particulate monolayer and it is the main force driving the fracture process. Modulating the surface tension difference also provides a way to probe mechanical properties that are difficult to measure directly, such as the critical failure stress. Finally, the surface tension difference can easily be controlled in the laboratory by modifying the composition of the surfactant or the underlying liquid.
In this article, we focus on closely-packed systems. We develop a phase-field model that takes into account a two-way coupling between the flow of surfactant and the motion of the monolayer. Through model-based simulations and accompanying experiments, we demonstrate that surface tension differences play a vital role in the overall fracture response of particle rafts.
The general setup for a surfactant injection experiment, illustrated in Fig. 1, consists of a circular Petri dish containing a liquid layer onto which hydrophobic particles are deposited. The surfactant is introduced with a needle near the center of the dish. The surface tension of the liquid–vapor interface
a Department of Civil and Environmental Engineering, Duke University, Durham,
NC 27708, USA. E-mail: [email protected] b Mathematics, Mechanics, and Materials Unit, OIST Graduate University,
Onna-son, Okinawa, 904-0495, Japan. E-mail: [email protected] c Collective Interactions Unit, OIST Graduate University, Onna-son, Okinawa,
904-0495, Japan
DOI: 10.1039/c7sm01245d
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5833
decreases where the surfactant is present. Marangoni forces then cause the surfactant to spread over the surface of the liquid and through the monolayer. The subsequent response is sensitive to the ratio of the fraction of the area of the liquid– vapor interface that is occupied by particles, which we refer to as the packing fraction and denote by f. As f is increased, the properties of the liquid–vapor interface change from liquid to solid.11
Consistent with out interest in closely-packed systems, we consider situations in which the packing fraction is high (0.7 r f r 0.9). We model the particle laden liquid–vapor interface as a continuous two-dimensional linearly elastic solid, capable of supporting both tension and compression.4,8 The mechanical properties of such a particulate monolayer12 are basic to understanding its response to surfactant-driven stresses. An estimate of those properties is provided by Vella et al.,8 based on experimental measurements and geometrical arguments. Surface wave experiments have been used to characterize the stretching and bending stiffnesses of particulate monolayers,13
which are found to present a soft granular character with nonclassical response under fluid-driven deformation.14
Experiments show that when such a system is stimulated by the localized introduction of a surfactant, an advancing front forms and fractures the monolayer.10 This type of surfactant- induced effect has also been observed in other systems, such as agar gels.15 The resulting fracture patterns, illustrated in Fig. 1, are reminiscent of those observed in classical brittle materials but with significant differences. Bandi et al.11 suggested that the fracture patterns can be very sensitive to variations in the initial distribution of particles. Additionally, in conventional elastic solids, crack branching is mostly associated with inertial effects and dynamic instabilities (e.g., bifurcations which occur as the crack tip velocity approaches 60% of the Rayleigh wave speed16). In contrast, crack branching has been observed in particulate mono- layers at crack speeds as low as 0.2% of the shear wave speed.10
Regarding the time scales of the fracture process, the most obvious distinction is that crack tip velocities do not appear to be influenced by the elastic properties of the monolayer. Rather, experimental observations suggest that the fracture time scale is mostly governed by variations in the viscosity of the underlying liquid10 and the packing fraction.11 These observations raise a number of questions regarding the fluid-driven fracture of elastic media in general and in monolayers in particular.
The mechanisms underlying the surfactant-driven fracture of particulate monolayers are challenging to study experimentally.
There are practical limits to the range of particle sizes that can be used. Sufficiently small particles fall below the current resolution limit of imaging equipment when a full field of view is needed. Particles that are too large have too much inertia to move significantly in response to surfactant flow. Modeling and simulation efforts can address these concerns and provide detailed insight concerning the basic mechanisms and sensi- tivities. However, computational studies of these systems are in the early stage of development. Previous attempts to model particulate monolayers have focused on loosely-packed systems, which admit several simplifying assumptions. For example, Bandi et al.,11 developed a discrete-element method to examine the influence of the initial packing fraction on the number of fractures in loosely-packed systems. In that work, a one-way coupling in which surfactant flow influences the motion of the particles, but not vice versa, as assumed. While simulations based on this approach accurately reproduce the limiting number of cracks that develop, the underlying assumptions limit its applicability to situations where the packing fraction is low.
The numerical simulation of closely-packed systems is challenging due to the large number of particles and the complexity of the physics involved in fracture. In this work, we propose a model based on a phase-field method which makes it possible to smoothly represent transitions between the damaged and undamaged zones of the monolayer as the surfactant advances. Phase-field models are suitable for this kind of problem since they avoid the need to track the propa- gating front, while allowing for a simple and powerful way to introduce the essential physics. We take as a starting point the work of Miehe et al.17 and Borden et al.18 for phase-field regularizations of the Griffith model for fracture. Based on these approaches, we present a new model that includes several features that are needed to characterize the fracture of particulate monolayers. In particular, our model incorporates the distribution of particles, the force on the monolayer due to the presence of surfactant, and the viscosity of the underlying liquid. Our objective is to develop a model capable of capturing the salient features of this system, namely the sensitivity of fracture patterns to differences in surface tension and the temporal scaling of the fracture process.
This paper is organized as follows. In Section 2, we present a phase-field approach for modeling the fracture mechanics of particulate monolayers, placing emphasis on the fundamental nature of the terms used to describe each contribution to the physics, after which we propose a governing free energy per unit area from which the model is derived. In Section 3, we describe the materials and experimental methodology used to explore the fidelity of the computational model. In Section 4, we present numerical results which demonstrate the capability of the model to reproduce different cases of surfactant-driven fracture. We then present a phase diagram which delineates different fracture regimes as a function of the surface tension difference and fracture resistance of the monolayer. Finally, in Section 5, we summarize our main results and propose direc- tions for further research.
Fig. 1 Left, center: States before and after the introduction of surfactant. The surfactant is injected at the center, decreasing the surface tension and generating an advancing front at which fracture initiates. Right: Micrograph showing the microstructure of the closely-packed particulate monolayer.
Paper Soft Matter
2 Phase-field model
We model the surfactant-driven fracture of particulate mono- layers by adapting recently developed phase-field approaches to fracture.17–19 The model is designed to capture how a drop of surfactant introduced at the center of the domain can spread to form cracks or fissures in the particulate monolayer (Fig. 1). Accordingly, our model incorporates a number of postulates that are based on experimental observations, as detailed below. Central among these is the assumption that the fractured zones are completely filled with surfactant, without appreciable pene- tration of the surfactant beyond the contours of the cracks. On this basis, we use a single ‘‘surfactant damage’’ field as an indicator function for the surfactant concentration and for the damage to the monolayer.
2.1 General considerations
Consider a monolayer of particles floating on the surface of a liquid layer. Let the two-dimensional region occupied by that monolayer be denoted by O. In Fig. 2(a), the dark portion of O represents the intact portion of the monolayer and the light portion of O represents the portion of the monolayer damaged by the surfactant. The state of the monolayer is described by two independent variables, its vector displacement field u and a scalar surfactant damage field d. As indicated in Fig. 2(a), d takes values in [0,1], with d = 1 representing completely damaged material. For closely-packed monolayers, experimental observations suggest that fracture occurs even for small values of the (infinitesimal) strain tensor
e ¼ 1
2 ruþ ðruÞT
: (1)
The microstructure of the monolayer is described by a scalar packing fraction f that takes values in [0,1] and depends on position x in O, as illustrated in Fig. 2(b). This field characterizes the ratio of the subset of the surface of the liquid layer that is
covered by particles to the total area of that surface and is thought to influence the overall fracture pattern, crack kinking, and branching.11 We assume that the particles are rigid, so that any local contraction or expansion of the monolayer is accom- modated solely by variations of the surface area not covered by particles. Thus, given an initial packing fraction f0, the packing fraction f depends on the trace, tre, of the strain e through
fðeÞ ¼ f0
1þ tre : (2)
We further assume that surfactant damage acts only to degrade the tensile resilience of the monolayer and that crack propagation is prohibited under compression. This is achieved by employing a spectral decomposition17,20 of e into positive and negative components e+ and e. It is then possible to define tensile and compressive strain-energy densities W+ and W through
WðeÞ ¼ E
2ð1þ nÞ ej j 2 þ nE ð1 nÞ2 treð Þ2 (3)
where E 4 0 and 0 o n o 1 are the Young’s modulus and Poisson’s ratio of the undamaged monolayer and the positive and negative trace operations tr+ and tr are defined in accord with
trþe ¼ tre; if tre 0;
0 otherwise;
0 otherwise:
( (5)
We assume that the free-energy density of the monolayer is a function c of e, d, and rd with the form
cðe;d;rdÞ¼ ð1dÞ2WþðeÞ|fflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflffl} tensile
þxðfðeÞÞWðeÞ|fflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflffl} compressive
þQðd;rdÞ|fflfflfflfflfflffl{zfflfflfflfflfflffl} fracture
:
As in conventional phase-field models for brittle fracture,17,18
the tensile contribution to c decays quadratically with the surfactant damage d. The remaining contributions of the energy are, however, nonstandard and therefore require further discussion.
The compressive contribution to c accounts for increases in compressive energy that accompany increases in the packing fraction through the jamming factor x. This factor penalizes compression beyond a jamming threshold fj and prevents packing fractions from exceeding a maximum value fm. The variation of x with f—and, thus, with reference to (2), the initial packing fraction f0 and strain e—is illustrated in Fig. 3. Its value starts at x = 1 for 0 r f r fj and increases monotonically for fj r fo fm, exhibiting a vertical asymptote as f - fm. The thresholds fj = 0.84 and fm = 0.9 correspond to a random close-packing in two space dimensions21 and the maximum packing density for two-dimensional discs, respectively. Numerically, it is more convenient to use expressions that monotonically increase the compressive contribution without becoming unbounded. In this work, we use an expression with the
Fig. 2 (a) Schematic of a fractured particulate monolayer. The fractured zone is represented as a smooth transition zone across which the damage field d takes values between zero and unity. (b) The microstructure of the elastic domain O is captured by an additional packing fraction field f representing the solid to fluid ratio of the medium.
Soft Matter Paper
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5835
form 1.0 + f (1.0 + tanh((f 0.5(fm fj))/l)), with f being the factor of amplification and l being the length scale controlling the width of the regularization.
For the contribution to c associated with fracture, we choose a modified version
Qðd;rdÞ ¼ GðdÞ 2
(7)
of the expression proposed by Miehe et al.,17 with G defined by
GðdÞ ¼ G0 gs gf þ 1 gs
gf
gf
d
(8)
where gf and gs are the surface tension of the liquid layer and the surfactant, respectively, G0 4 0 is a constant that represents the fracture toughness of the undamaged monolayer, and l4 0 is a constant proportional to the characteristic thickness of a layer between damaged and undamaged material. This modifica- tion accounts for a reduction in the fracture toughness with increasing surfactant concentration. Accordingly, it is convenient to introduce
Gr ¼ Gð1Þ ¼ G0 gs gf
(9)
as a measure of the reduced fracture toughness of the monolayer. Finally, the surfactant contribution to c accounts for the
interplay between the monolayer and the surfactant (Fig. 4) and is assumed to be of the form
Fðe; dÞ ¼ F0 gf gsð ÞðtreÞd2
2 ; (10)
2.2 Governing equations
Following Silva et al.19 (but neglecting inertia), the governing equations of the phase-field model consist of the macroforce balance
div @cðe; d;rdÞ
8< : (12)
where b Z 0 represents the kinetic modulus that controls the rate at which cracks can propagate through the monolayer. The alternative on the right-hand side of (12) embodies the requirement that, consistent with experimental observations, cracks that form in the monolayer never heal. This requirement takes the form of the constraint
: d Z 0, and a reactive micro-
force pr is needed to ensure satisfaction of that constraint. In particular, pr vanishes for
: d 4 0 and is determined by the
left-hand side of (12) for : d = 0.
The kinetic modulus b incorporates two effects. First, it captures the capacity of the surfactant to penetrate the parti- culate monolayer and thereby generate damage, which is influenced by f. As f increase, the higher density of particles impede surfactant spreading. Second, it captures the resistance of the underlying liquid layer to rearrangements of the parti- cles, a resistance which is directly related to the viscosity mf
of the liquid comprising that layer. For the foregoing reasons, we assume that b increases with f and mf in accord with the relation
b ¼ b0 1 fm
1 f ; b0 ¼
; (13)
where C 4 0 is a dimensionless constant and b0 represents the kinetic modulus for maximal packing f = fm.
Since the particular choice of the kinetic modulus b defined in (13) is positive, damage only increases when the left-hand side of (12) is positive. Accordingly, the evolution eqn (12) for d becomes
1 fðeÞ b0 1 fmð Þ div
@cðe; d;rdÞ @rd
hhi ¼ 0; h 0;
h; h4 0;
( (15)
denotes its Macaulay bracket. With reference to the right-hand side of the definition (6) of
the free-energy density c and using (7)–(10) and (13) with reference to the parameters listed in Table 1, we find the
Fig. 3 Variation in the jamming factor x with the packing fraction f. The jamming factor tends towards N as f surpasses the jamming threshold fj = 0.84 and approaches the maximum packing fraction fm = 0.9.
Fig. 4 Representation of the surfactant force driving the system. A collection of particles, initially in equilibrium at a constant surface tension gf, are reached by the surfactant. The gradient between its lower surface tension gs and the surrounding media generates a Marangoni flow which pushes the particulate monolayer outward.
Paper Soft Matter
governing eqn (11) and (14) can be written as
div ð1 dÞ2@WþðeÞ @e
2 I
¼ 0
1 fðeÞ b0 1 fmð Þ G0l 1 1 gs
gf
d
G0d
respectively, and where I denotes the (two-dimensional) identity tensor. We consider (16) and (17) subject to the vanishing displacement condition
u = 0 (18)
on the boundary qO of O and, letting n denote a unit normal on qO, the natural boundary condition
rdn = 0. (19)
f(x,0) = f0(x), d(x,0) = d0(x), (20)
with f0 representing the initial distribution of particles and d0 the initial damage.
2.3 Model characterization: dimensionless number v
Introducing characteristic measures L and T of length and time, we define dimensionless counterparts x* and t* of x and t by
x ¼ x
T : (21)
rd ¼ Lrd; Dd ¼ L2Dd; and @d
@t ¼ T _d (22)
of the gradient, Laplacian, and partial time-derivative, dimen- sionless parameters,
G ¼ Gr
E ; L ¼ L
; g ¼ gs gf ; and l ¼ l
L ;
and dimensionless strain-energy densities
W ðeÞ ¼ 1
2ð1þ nÞ ej j 2 þ n ð1 nÞ2 treð Þ2; (24)
we arrive at dimensionless counterparts
div ð1 dÞ2@Wþ ðeÞ
@e þ F
G
g l 1 1 gð Þdð ÞDd G
g l
þ 2ð1 dÞLWþðeÞ LFðtreÞd ¼ @d
@t
(26)
of the governing eqn (16) and (17). While we have not undertaken an exhaustive suite of simula-
tions spanning the complete range of parameter space associated with the six dimensionless variables in (23), the extensive testing that we have conducted indicates that the fracture pattern is largely governed by the dimensionless driving force F* and the dimensionless fracture toughness G* of the monolayer. In particular, the number of fractures in the final configuration appears to be dictated by the ratio
w ¼ Fð Þ2
G : (27)
In Section 4.2, we discuss the threshold levels of the dimen- sionless parameter w that delineate different regimes, ranging from no fractures at all to a configuration with multiple branches.
The expression (27) for w can be further simplified by invoking the approximate scaling of Young’s modulus with the surface tension as derived by Vella et al.8 (i.e., E p gf/%rp). The expression for w can therefore be rewritten to isolate the influence of the surface tensions, yielding
w / F0rp 2 gf gsð Þ2
G0gs : (28)
Daniels et al.22 proposed a similar relationship between the elastic and Marangoni energies to explain fracture patterns in agar gels. In that work, the number of crack branches was found to scale linearly with the difference in surface tension. As described in Section 4, our experimental observations and
Table 1 List of parameters in the model
Parameter Description
fj Jamming packing fraction of monolayer fm Maximum packing fraction of monolayer %rp Average particle radius gf Surface tension of liquid layer gs Surface tension of surfactant G0 Fracture toughness of undamaged monolayer F0 Surfactant force constant E Young’s modulus of monolayer n Poisson’s ratio of monolayer mf Dynamic viscosity of monolayer b0 Kinetic parameter
Soft Matter Paper
model-based simulations suggest the system considered here exhibits a much stronger sensitivity to surface tension differences.
3 Experiments 3.1 Materials
All videos from experimental work were acquired using a Phantom V641 high-speed camera equipped with an AF Nikkor 50 mm f/1.8 D lens. Videos were saved using PCC software provided by Phantom. A 5 W white LED board used to illuminate the mixture was supplied by VANCO (series #33342). Dispersion of silica microballoons (SIG MFG) was performed using the Active Motif Q120AM (120 W, 20 kHz) ultrasonicator equipped with a CL18 3.2 mm probe. Oleic acid (surfactant) was supplied by WAKO Chemicals and the acetone used as a cosolvent was supplied by Nacalai Tesque. The surfactant was dropped using Terumo 2.5 mL syringes (SS-02SZ) equipped with a 25 G (0.50 16 mm) needle. Surface tension measurements of mixtures were obtained using the KSV NIMA LB Small trough (33473).
3.2 Methodology
Fig. 5 shows a schematic of the experimental setup. A clean low- form cylindrical glass vessel 17 cm in diameter and 9 cm in height filled to a height of 2.5 cm with Milli-Q water (or mixed with a cosolvent) was placed atop an LED light panel. Micro- balloons were carefully weighed to either 0.150 g or 0.250 g, introduced at the air–water interface, and dispersed by means of ultrasonication at 25% power to form a particulate mono- layer. Using a 2.5 mL syringe fitted with a 25 G metal tip, a single drop of surfactant (oleic acid with or without solvent) was introduced to the approximate center of the vessel at time t = 0 s. The water immiscible surfactant was observed to spread and push the microballoon particles radially outwards resulting in compaction of the particulate monolayers. Observations were captured using a Phantom high-speed camera at different frame rates. Data analysis was performed by first converting
raw data files (.cine) to multipage image files (.tif) and these files were subsequently analyzed using ImageJ. In cases where the surfactant or the bulk water phase was mixed with a solvent to reduce surface tension, acetone was added carefully and aliquots of the liquids were quickly transferred to an LB trough for surface tension measurements.
4 Results and discussion
We discretize the model described in Section 2 by recasting the evolution equations in a variational form and approximating the displacement and surfactant damage fields with finite element basis functions. Temporal integration is performed with a second-order accurate, unconditionally stable backward Euler algorithm (BDF2) with adaptive time stepping. The coupled system of nonlinear equations are solved with the Newton method in a staggered manner (i.e., alternating the updates for the displacement and damage fields).
4.1 Pattern comparison
Initial-boundary value problems are constructed over circular domains representative of the petri dishes used in the experi- ments. The experiments have a stochastic aspect corresponding to the initial particle locations. This effect is modeled at the continuum scale by using initial packing fields with spatial variability. Specifically, we numerically construct initial packing fraction fields over the circular domain that have a mean value of f = 0.795 with a 6% variation, corresponding to uniform random distributions with minimum and maximum packing fractions of 0.75 and 0.84, respectively. The initial drop of surfactant is approximated by a sufficiently small (on the order of l), empty circular region (f0 = 0) concentric with the center of the dish. Table 2 provides the reference values of the parameters used for the simulations in this section. All simulations were per- formed using MOOSE (Multiphysics Object-Oriented Simulation Environment), a finite-element framework developed primarily by Idaho National Laboratory.23
A first qualitative assessment of the model is shown in Fig. 6. In this example, the parameter w is adjusted to reproduce a three-branch pattern, which is a commonly observed configu- ration in the experiments for our given set of parameters. The result indicates that the model captures the general fracture process and reproduces the final three-branch configuration. Additionally, the cracks are sensitive to the packing fraction structure, which gives rise to more subtle features. These features include the particular angle in which the three- branch pattern is generated, the kinking and bending of cracks when gradients in the particle density are encountered, and the bifurcation of a crack when it encounters a high density region directly in line with its propagation. In this work, stochastic features of the experiments are approximated through random variations in the initial packing fraction field. The specific details of the numerically generated fracture patterns thus stem from the initial conditions and are sensitive to the particular random seed used in the simulations. As a result, our
Fig. 5 Setup: a clean low-form cylindrical vessel containing distilled water was placed on an LED tablet. A camera suspended vertically above records particulate monolayer compaction dynamics at the air–water interface when surfactant was introduced using a 25 G steel needle.
Paper Soft Matter
5838 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
simulation are designed to capture only the generic features of the experiments. Nevertheless, the model reproduces the number of cracks observed in the experiments and, at a qualitative level, the features of crack kinking and the for- mation of secondary branches.
We now provide a detailed description of a representative simulation in which a three-branch fracture pattern arises. Fig. 7 shows the results from a simulation of this common pattern, obtained with the parameters given in Table 2 and w = 0.0125 (see Table 3 for the correspondent dimensionless quantities). The columns in Fig. 7 show the evolution of three fields, namely the surfactant damage, the tensile energy, and the packing fraction, along with a comparison with snapshots of a matching experiment at comparable moments in time. The distribution of d closely matches the fractured portion of the particle raft in the experiment. The tensile strain-energy density field W+, which is representative of the processes that underpin crack propagation, transitions from an initially circular tensile distribution due to the surfactant force to one with broken symmetry depending on the available energy in the system. The tensile strain-energy density drives the three- branches until they eventually arrest due to an increase in the compressive energy. The field f depicts the evolution of the
initial packing fraction. As the fractured region expands, the monolayer contracts (to approximately f = 0.87) and the mobility of the surfactant in the fractured region (f = 0.3) increases.
4.2 Phase diagram for fracture pattern characterization
We now study whether the dimensionless parameter w can be used to construct a phase diagram for the fracture patterns that are produced. In so doing, we eliminate the random aspect of the previous simulations and consider initial packing fractions that are spatially uniform with f0 = 0.80. As illustrated in Fig. 8, the results allow us to identify threshold levels of w that delineate regions of a phase diagram for the fracture patterns obtained through simulations. The threshold w = 0.005 defines a limit below which fracture is not observed. As w is increased,
Table 2 Dimensional variables8,10
Property Description Value Unit
rd Radius of Petri dish 8.5 102 m f0 Range of initial packing fraction (7.5–8.4) 101 None fm Packing fraction of jammed monolayer 9.0 101 None G0 Fracture toughness of undamaged monolayer 5.7 (104–102) N m1
F0 Surfactant force constant 1.0 (101–102) m1
gf Surface tension of liquid layer (water) 7.2 102 N m1
gs Surface tension of surfactant (oleic acid) 4.0 102 N m1
mf Range of dynamic viscosity of liquid layer 1.0 (103–101) Pa s %rp Mean particle radius 1.0 104 m E Young’s modulus of monolayer 3.178 103 Pa n Poisson’s ratio of monolayer 1
ffiffiffi 3 p
None
b0 Kinetic parameter 1.0 Pa s l Characteristic length scale 1.0 102 m
Fig. 6 Comparison between the result of a simulation (left) and an experimental image (right) for a three-branch fracture pattern. The model captures the main traits of the final configuration, as represented by the number of fractures, and more detailed features corresponding to varia- tions in the distribution of particles, such as secondary fracturing, crack bending, and branching. The particular orientation of the fracture pattern shown in the simulation result is arbitrary, as it is sensitive to the random seed used to construct the initial packing field.
Fig. 7 Numerical fields obtained for a simulation of a three-branch fracture pattern. From left to right, each column shows the initial (top) and final (bottom) states of surfactant damage d, tensile strain-energy density W+ and f fields. The rightmost column shows the reference experiment, obtained by adding oleic acid to a 0.25 g mass monolayer of particles laying on a 18% salt saturation water solution.
Table 3 Dimensionless parameters
L* 8.5 102
b*0 2.7 (108–106) g* 5.6 101
l* 1.5 101
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5839
the fracture pattern starts to develop as a single crack, and then into three-branch shapes with further increases. We note that the three-branch configuration is quite ubiquitous in nature, as the presence of acute angles and the compression generated by the three cracks stabilizes the system. For larger values of w, a fourth and a fifth branch are gradually realized. Larger values of w are required to produce larger number of crack branches, resulting in unstable fracture patterns that cannot be clearly delineated in well defined regions of our phase diagram.
To validate the zones and thresholds predicted by the simulations, we first calibrated our model to experiments using oleic acid spreading over a pure water layer. The resulting fracture pattern is shown in the bottom of the center column of Fig. 9, which presents images of the final configurations for three different experiments, along with simulation results for comparable values of w. Numerical simulations indicate that such a pattern can be obtained using w = 0.018, which we label as wref. As indicated by expression (28), changing the surface tension difference corresponds to modulating the value of w in the model. Accordingly, we conducted two additional experi- ments. In the first of these experiments, the surface tension of the liquid layer was reduced from 72 to 56 mN m1 (colsolvent). In the second, the surface tension of the surfactant was reduced from 40 to 28 mN m1 (acetone). These changes in surface tension correspond approximately to modifying factors of 1/3 and 2, respectively, on the reference wref = 0018. We denote the corresponding levels of w as wdown and wup.
As illustrated in Fig. 9, simulations based on the values of wup and wdown yield fracture patterns that are remarkably consistent with the corresponding experiments. Simulations based on doubling w from wref to wup result in an increase in the number of crack branches, from three to five. Conversely, simulations based on a reduction of w from wref to wdown result in a decrease in the number of crack branches, from three to two. In each case, a qualitative match with the corresponding experiment is clear. As w is increased, the model presents a gradual change between patterns, from the formation of a
single crack to a configuration with multiple branches. The model nicely reproduces the central damage zone and multiple crack branches observed in experiments for the range (104–101) of w in this study.
This result is particularly important because the conven- tional view has been that the packing fraction distribution is the primary variable governing the morphological features of the crack pattern. Our results suggest that, although the randomness of the initial packing fraction can underpin features such as crack deflections and arrest, the branching and general patterns observed in the experiments correspond to a much more fundamental principle related to the interplay between the surfactant force and the fracture toughness. In other words, the difference in surface tension between the surfactant and underlying liquid are central to determining the crack patterns that emerge.
4.3 Temporal scaling
We now investigate the extent to which the model captures the speed of crack propagation. We examine the overall time required to fracture the particulate monolayer as well as the temporal scaling of the process. The snapshots in Fig. 10 show two fracture patterns that evolve at very different rates, but for which the final configurations are visually indistinguishable. For comparison purposes, we show our results using a dimensional version of the model (i.e., viscosity in Pa s). Previous experimental observa- tions have indicated that, for the same initial packing fractions, the geometry of the fracture pattern is insensitive to the viscosity of the underlying liquid layer.10 In the experiments, the liquid viscosity can be adjusted, for example, by increasing the ratio of glycerol in the glycerol–water mix based liquid layer, raising the viscosity from 103 Pa s to 101 Pa s. This is consistent with work showing that the rate at which the process occurs decreases by the same order of magnitude.10 We can capture this effect by adjusting the mobility implicitly defined in (13), which is propor- tional to the viscosity mf of the liquid layer.
Fig. 8 Phase diagram for fracture pattern characterization. Threshold levels of the dimensionless number w = (F*)2/G* delineate various patterns. In the region below w = 0.005 (zone i), the initially fractured region does not grow with time. As w increases, the fracture patterns gradually change from the single branch configuration (zone ii) to the three (zone iii), four (zone iv) and five (zone v) branch configurations.
Fig. 9 Comparison between fracture patterns from experiments (top) and simulations (bottom). In all cases, the final configuration is shown. Simula- tion results were obtained for values of the dimensionless parameter w (wdown = 0.006, wref = 0018, wup = 0.035), achieved by modulating the surface tension of the liquid layer and the surfactant used in each of the experiments, consistent with (28).
Paper Soft Matter
5840 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
Additionally, experimental measurements indicate that the fracture proceeds in two stages. Crack growth rates that scale with the 3/4 and higher powers of t are observed in the early stages of the fracture process, followed by a gradual flattening as the crack advances and eventually arrests.10 Arrest may occur in response to the increase in packing around the perimeter of the domain, which gradually makes it more difficult to displace the particles of the monolayer. As shown in Fig. 11, our model qualitatively reproduces this behavior with time, in terms of the normalized fractured area A* (namely the ratio of the area of the surface of the liquid layer that is exposed to the area of O). As argued by Bandi et al.,11 if the radially receding front of R(t) of a particle raft scales according to R(t) B t3/4, then the particle free area proceeds as R(t)2 B (t3/4)2 B t3/2. The growth rates in the early stages are slightly below t3/2, which suggests the presence of the monolayer slightly retards the expansion of the surfactant from what would be expected over a pure liquid surface.24 Finally, we recognize a shift between the results from experiment and simulation, with the initial expansion process
being noticeably faster in the experiment. This discrepancy at early times could be related to the rapid expansion in the drop of surfactant, a process that is not captured by the model. The subtleties of the physics at early times are likely sensitive to the volume and height of the drop of surfactant at the time of injection. Questions related to these subtleties are a subject for future work.
Following the rapid initial expansion, the dynamics slow down considerably once the particulate rafts jam together. As shown in the inset to Fig. 11, the late stage dynamics exhibit approximately logarithmically slow relaxation over a decade in time. This kind of dynamics, also known as ‘‘creep’’ or ‘‘aging’’, is observed in a variety of phenomena including granular compaction,25 electron glasses,26 polymer relaxation,27 and super- conductors,28 and is now considered a hallmark of non- exponential, slow relaxation in amorphous media.29 Although the observation times in our model and experiments do not extend over several decades, they are consistent with the expected behavior of slow relaxation in frustrated, amorphous media.
5 Conclusions
We have presented a phase-field model that captures the main features of the surfactant-driven fracture of closely-packed particulate monolayers. The numerical simulations and experi- mental results indicate that there is a competition between the spreading force of the surfactant and the fracture toughness of the monolayer that determines the number of fractures. The model gives rise to a dimensionless parameter that can be written in terms of the surface tensions, and which allows for a straightforward comparison with experimental conditions. Experiments were conducted to validate the model and delineate the regimes of fracture pattern as a function of that parameter.
Our model rests on a number of simplifying assumptions. In particular, we used the damage to the particulate monolayer as an indicator function for the surfactant concentration. As such, our model precludes investigations of the extent to which the surfactant can penetrate into the particulate monolayer network, ahead of (or behind) the crack front. We have also assumed that the strains exhibited by the monolayer are sufficiently small to justify modeling it as a linearly elastic solid. This assumption ceases to be reasonable if the surface tension difference is increased above some critical threshold. Moreover, its applic- ability is limited to systems in which the initial packing is low enough to allow the particulate layer to exhibit both loosely- packed and jammed behaviors as the surfactant spreads.11 In either case, transitions between linear and nonlinear response, or between fluid and solid behavior, are worthy of further study. Finally, the interplay between liquid and solid constituents in the monolayer may endow it with a viscoelastic behavior that could also be explored. As the rate of surfactant transport is a function of the Marangoni stress, modulating the surface tension difference might also be used to explore any rate depen- dence in the monolayer. Future work will focus on enhancing the
Fig. 10 Identical fracture patterns are produced at different speeds due to changes in the viscosity mf of the underlying liquid layer.
Fig. 11 Comparison between simulation and experiment for the temporal scaling for the normalized fractured area A* (the ratio of the portion of the area of the liquid layer that is exposed to the area of O), indicating two regimes of crack evolution. For reference, a t3/2 curve is provided (red dashed line). The inset presents a zoom of the latter stages using semi-log axes, revealing a slow logarithmic relaxation of the fracture pattern.
Soft Matter Paper
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5841
model to incorporate these effects and to allow for detailed exploration of their consequences and significance.
Acknowledgements
Christian Peco, Yingjie Liu, and John E. Dolbow are grateful for the support of NSF grant CMMI-1537306, to Duke University. Wei Chen, M. M. Bandi, and Eliot Fried gratefully acknowledge support from the Okinawa Institute of Science and Technology Graduate University with subsidy funding from the Cabinet Office, Government of Japan. John E. Dolbow would also like to acknowledge support from the Okinawa Institute of Science and Technology during a sabbatical stay in late 2014, during which time this work began.
References
1 A. Varshney, A. Sane, S. Ghosh and S. Bhattacharya, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2012, 86, 31402.
2 S. U. Pickering, J. Chem. Soc., Trans., 1907, 91, 2001–2021. 3 B. P. Binks and P. D. I. Fletcher, Langmuir, 2001, 17,
4708–4710. 4 P. Cicuta, E. J. Stancik and G. G. Fuller, Phys. Rev. Lett., 2003,
90, 236101. 5 P. Aussillous and D. Quere, Nature, 2001, 411, 924–927. 6 N. Tsapis, D. Bennett, B. Jackson, D. A. Weitz and D. A. Edwards,
Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 12001–12005. 7 A. B. Subramaniam, M. Abkarian and H. A. Stone, Nat. Mater.,
2005, 4, 553–556. 8 D. Vella, P. Aussillous and L. Mahadevan, Europhys. Lett.,
2004, 68, 212–218. 9 L. Scriven and C. Sterling, Nature, 1960, 187, 186–188.
10 D. Vella, H. Y. Kim, P. Aussillous and L. Mahadevan, Phys. Rev. Lett., 2006, 96, 1–4.
11 M. M. Bandi, T. Tallinen and L. Mahadevan, Europhys. Lett., 2011, 96, 36008.
12 P. Cicuta and D. Vella, Phys. Rev. Lett., 2009, 102, 138302. 13 C. Planchette, E. Lorenceau and A.-L. Biance, Soft Matter,
2012, 8, 2444–2451. 14 C. W. MacMinn, E. R. Dufresne and J. S. Wettlaufer,
Phys. Rev. X, 2015, 5, 11020. 15 C. Spandagos, T. B. Goudoulas, P. F. Luckham and O. K.
Matar, Langmuir, 2012, 28, 7197–7211. 16 E. H. Yoffe, Philos. Mag., 1951, 42, 739–750. 17 C. Miehe, M. Hofacker and F. Welschinger, Comput. Methods
Appl. Mech. Eng., 2010, 199, 2765–2778. 18 M. J. Borden, C. V. Verhoosel, M. A. Scott, T. J. R. Hughes
and C. M. Landis, Comput. Methods Appl. Mech. Eng., 2012, 217-220, 77–95.
19 M. N. da Silva Jr., F. P. Duda and E. Fried, J. Mech. Phys. Solids, 2013, 61, 2178–2195.
20 C. Miehe, Comput. Struct., 1998, 66, 37–43. 21 C. S. O’Hern, L. E. Silbert, A. J. Liu and S. R. Nagel, Phys. Rev.
E: Stat., Nonlinear, Soft Matter Phys., 2003, 68, 11306. 22 K. E. Daniels, S. Mukhopadhyay, P. J. Houseworth and
R. P. Behringer, Phys. Rev. Lett., 2007, 99, 124501. 23 D. Gaston, C. Newman, G. Hansen and D. Lebrun-Grandie,
Nucl. Eng. Des., 2009, 239, 1768–1778. 24 O. E. Jensen, J. Fluid Mech., 1995, 293, 349–378. 25 E. Ben-Naim, J. B. Knight, E. R. Nowak, H. M. Jaeger and
S. R. Nagel, Phys. D, 1998, 123, 380–385. 26 A. Amir, Y. Oreg and Y. Imry, Phys. Rev. Lett., 2009,
103, 126403. 27 J. D. Ferry, Viscoelastic Properties of Polymers, John Wiley and
Sons, 1980. 28 A. Gurevich and H. Kupfer, Phys. Rev. B: Condens. Matter
Mater. Phys., 1993, 48, 6477–6487. 29 A. Amir, Y. Oreg and Y. Imry, Proc. Natl. Acad. Sci. U. S. A.,
2012, 109, 1850–1855.
PAPER John E. Dolbow, Eliot Fried et al. Influence of surface tension in the surfactant-driven fracture of closely-packed particulate monolayers
Volume 13 Number 35 21 September 2017 Pages 5811–6022
5832 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
Cite this: SoftMatter, 2017,
Influence of surface tension in the surfactant-driven fracture of closely-packed particulate monolayers
Christian Peco, a Wei Chen,b Yingjie Liu, a M. M. Bandi, c
John E. Dolbow *a and Eliot Fried *b
A phase-field model is used to capture the surfactant-driven formation of fracture patterns in particulate
monolayers. The model is intended for the regime of closely-packed systems in which the mechanical
response of the monolayer can be approximated as that of a linearly elastic solid. The model
approximates the loss in tensile strength of the monolayer with increasing surfactant concentration
through the evolution of a damage field. Initial-boundary value problems are constructed and spatially
discretized with finite element approximations to the displacement and surfactant damage fields.
A comparison between model-based simulations and existing experimental observations indicates a
qualitative match in both the fracture patterns and temporal scaling of the fracture process. The
importance of surface tension differences is quantified by means of a dimensionless parameter,
revealing thresholds that separate different regimes of fracture. These findings are supported by newly
performed experiments that validate the model and demonstrate the strong sensitivity of the fracture
pattern to differences in surface tension.
1 Introduction
When a densely packed monolayer of hydrophobic particles is placed on the surface of a liquid, the particles interact through capillary bridges,1 leading to the formation of particle rafts. The macroscopic properties of these rafts reflect an interplay between fluid and solid mechanics,2–4 giving rise to novel physics. This interplay is relevant to a wide range of applica- tions, from the synthesis of ‘‘liquid marbles’’5 to the design of drug delivery systems6 to the stabilization of drops.7
The interest in particle rafts has driven researchers to investigate their mechanical properties.4,8 It is now known that densely packed monolayers exhibit a two-dimensional linearly elastic solid response. The interaction of the particles through the capillary bridges allows the monolayer to support both tension and compression. Moreover, the mechanical properties of such monolayers depend proportionally on the surface tension of the liquid layer. The introduction of a controlled amount of surfactant generates a surface tension gradient,
producing Marangoni forces9 and causing the surfactant to spread, fracturing the monolayer. Studies of surfactant-driven fracture have examined the role of viscosity and the initial packing fraction on the evolution of cracks in closely and loosely-packed systems, respectively.10,11 Surprisingly, the potential importance of differences in the surface tension of the surfactant and the underlying liquid has not been explored. The magni- tude of this difference is interesting because it determines the Marangoni force exerted on the particulate monolayer and it is the main force driving the fracture process. Modulating the surface tension difference also provides a way to probe mechanical properties that are difficult to measure directly, such as the critical failure stress. Finally, the surface tension difference can easily be controlled in the laboratory by modifying the composition of the surfactant or the underlying liquid.
In this article, we focus on closely-packed systems. We develop a phase-field model that takes into account a two-way coupling between the flow of surfactant and the motion of the monolayer. Through model-based simulations and accompanying experiments, we demonstrate that surface tension differences play a vital role in the overall fracture response of particle rafts.
The general setup for a surfactant injection experiment, illustrated in Fig. 1, consists of a circular Petri dish containing a liquid layer onto which hydrophobic particles are deposited. The surfactant is introduced with a needle near the center of the dish. The surface tension of the liquid–vapor interface
a Department of Civil and Environmental Engineering, Duke University, Durham,
NC 27708, USA. E-mail: [email protected] b Mathematics, Mechanics, and Materials Unit, OIST Graduate University,
Onna-son, Okinawa, 904-0495, Japan. E-mail: [email protected] c Collective Interactions Unit, OIST Graduate University, Onna-son, Okinawa,
904-0495, Japan
DOI: 10.1039/c7sm01245d
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5833
decreases where the surfactant is present. Marangoni forces then cause the surfactant to spread over the surface of the liquid and through the monolayer. The subsequent response is sensitive to the ratio of the fraction of the area of the liquid– vapor interface that is occupied by particles, which we refer to as the packing fraction and denote by f. As f is increased, the properties of the liquid–vapor interface change from liquid to solid.11
Consistent with out interest in closely-packed systems, we consider situations in which the packing fraction is high (0.7 r f r 0.9). We model the particle laden liquid–vapor interface as a continuous two-dimensional linearly elastic solid, capable of supporting both tension and compression.4,8 The mechanical properties of such a particulate monolayer12 are basic to understanding its response to surfactant-driven stresses. An estimate of those properties is provided by Vella et al.,8 based on experimental measurements and geometrical arguments. Surface wave experiments have been used to characterize the stretching and bending stiffnesses of particulate monolayers,13
which are found to present a soft granular character with nonclassical response under fluid-driven deformation.14
Experiments show that when such a system is stimulated by the localized introduction of a surfactant, an advancing front forms and fractures the monolayer.10 This type of surfactant- induced effect has also been observed in other systems, such as agar gels.15 The resulting fracture patterns, illustrated in Fig. 1, are reminiscent of those observed in classical brittle materials but with significant differences. Bandi et al.11 suggested that the fracture patterns can be very sensitive to variations in the initial distribution of particles. Additionally, in conventional elastic solids, crack branching is mostly associated with inertial effects and dynamic instabilities (e.g., bifurcations which occur as the crack tip velocity approaches 60% of the Rayleigh wave speed16). In contrast, crack branching has been observed in particulate mono- layers at crack speeds as low as 0.2% of the shear wave speed.10
Regarding the time scales of the fracture process, the most obvious distinction is that crack tip velocities do not appear to be influenced by the elastic properties of the monolayer. Rather, experimental observations suggest that the fracture time scale is mostly governed by variations in the viscosity of the underlying liquid10 and the packing fraction.11 These observations raise a number of questions regarding the fluid-driven fracture of elastic media in general and in monolayers in particular.
The mechanisms underlying the surfactant-driven fracture of particulate monolayers are challenging to study experimentally.
There are practical limits to the range of particle sizes that can be used. Sufficiently small particles fall below the current resolution limit of imaging equipment when a full field of view is needed. Particles that are too large have too much inertia to move significantly in response to surfactant flow. Modeling and simulation efforts can address these concerns and provide detailed insight concerning the basic mechanisms and sensi- tivities. However, computational studies of these systems are in the early stage of development. Previous attempts to model particulate monolayers have focused on loosely-packed systems, which admit several simplifying assumptions. For example, Bandi et al.,11 developed a discrete-element method to examine the influence of the initial packing fraction on the number of fractures in loosely-packed systems. In that work, a one-way coupling in which surfactant flow influences the motion of the particles, but not vice versa, as assumed. While simulations based on this approach accurately reproduce the limiting number of cracks that develop, the underlying assumptions limit its applicability to situations where the packing fraction is low.
The numerical simulation of closely-packed systems is challenging due to the large number of particles and the complexity of the physics involved in fracture. In this work, we propose a model based on a phase-field method which makes it possible to smoothly represent transitions between the damaged and undamaged zones of the monolayer as the surfactant advances. Phase-field models are suitable for this kind of problem since they avoid the need to track the propa- gating front, while allowing for a simple and powerful way to introduce the essential physics. We take as a starting point the work of Miehe et al.17 and Borden et al.18 for phase-field regularizations of the Griffith model for fracture. Based on these approaches, we present a new model that includes several features that are needed to characterize the fracture of particulate monolayers. In particular, our model incorporates the distribution of particles, the force on the monolayer due to the presence of surfactant, and the viscosity of the underlying liquid. Our objective is to develop a model capable of capturing the salient features of this system, namely the sensitivity of fracture patterns to differences in surface tension and the temporal scaling of the fracture process.
This paper is organized as follows. In Section 2, we present a phase-field approach for modeling the fracture mechanics of particulate monolayers, placing emphasis on the fundamental nature of the terms used to describe each contribution to the physics, after which we propose a governing free energy per unit area from which the model is derived. In Section 3, we describe the materials and experimental methodology used to explore the fidelity of the computational model. In Section 4, we present numerical results which demonstrate the capability of the model to reproduce different cases of surfactant-driven fracture. We then present a phase diagram which delineates different fracture regimes as a function of the surface tension difference and fracture resistance of the monolayer. Finally, in Section 5, we summarize our main results and propose direc- tions for further research.
Fig. 1 Left, center: States before and after the introduction of surfactant. The surfactant is injected at the center, decreasing the surface tension and generating an advancing front at which fracture initiates. Right: Micrograph showing the microstructure of the closely-packed particulate monolayer.
Paper Soft Matter
2 Phase-field model
We model the surfactant-driven fracture of particulate mono- layers by adapting recently developed phase-field approaches to fracture.17–19 The model is designed to capture how a drop of surfactant introduced at the center of the domain can spread to form cracks or fissures in the particulate monolayer (Fig. 1). Accordingly, our model incorporates a number of postulates that are based on experimental observations, as detailed below. Central among these is the assumption that the fractured zones are completely filled with surfactant, without appreciable pene- tration of the surfactant beyond the contours of the cracks. On this basis, we use a single ‘‘surfactant damage’’ field as an indicator function for the surfactant concentration and for the damage to the monolayer.
2.1 General considerations
Consider a monolayer of particles floating on the surface of a liquid layer. Let the two-dimensional region occupied by that monolayer be denoted by O. In Fig. 2(a), the dark portion of O represents the intact portion of the monolayer and the light portion of O represents the portion of the monolayer damaged by the surfactant. The state of the monolayer is described by two independent variables, its vector displacement field u and a scalar surfactant damage field d. As indicated in Fig. 2(a), d takes values in [0,1], with d = 1 representing completely damaged material. For closely-packed monolayers, experimental observations suggest that fracture occurs even for small values of the (infinitesimal) strain tensor
e ¼ 1
2 ruþ ðruÞT
: (1)
The microstructure of the monolayer is described by a scalar packing fraction f that takes values in [0,1] and depends on position x in O, as illustrated in Fig. 2(b). This field characterizes the ratio of the subset of the surface of the liquid layer that is
covered by particles to the total area of that surface and is thought to influence the overall fracture pattern, crack kinking, and branching.11 We assume that the particles are rigid, so that any local contraction or expansion of the monolayer is accom- modated solely by variations of the surface area not covered by particles. Thus, given an initial packing fraction f0, the packing fraction f depends on the trace, tre, of the strain e through
fðeÞ ¼ f0
1þ tre : (2)
We further assume that surfactant damage acts only to degrade the tensile resilience of the monolayer and that crack propagation is prohibited under compression. This is achieved by employing a spectral decomposition17,20 of e into positive and negative components e+ and e. It is then possible to define tensile and compressive strain-energy densities W+ and W through
WðeÞ ¼ E
2ð1þ nÞ ej j 2 þ nE ð1 nÞ2 treð Þ2 (3)
where E 4 0 and 0 o n o 1 are the Young’s modulus and Poisson’s ratio of the undamaged monolayer and the positive and negative trace operations tr+ and tr are defined in accord with
trþe ¼ tre; if tre 0;
0 otherwise;
0 otherwise:
( (5)
We assume that the free-energy density of the monolayer is a function c of e, d, and rd with the form
cðe;d;rdÞ¼ ð1dÞ2WþðeÞ|fflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflffl} tensile
þxðfðeÞÞWðeÞ|fflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflffl} compressive
þQðd;rdÞ|fflfflfflfflfflffl{zfflfflfflfflfflffl} fracture
:
As in conventional phase-field models for brittle fracture,17,18
the tensile contribution to c decays quadratically with the surfactant damage d. The remaining contributions of the energy are, however, nonstandard and therefore require further discussion.
The compressive contribution to c accounts for increases in compressive energy that accompany increases in the packing fraction through the jamming factor x. This factor penalizes compression beyond a jamming threshold fj and prevents packing fractions from exceeding a maximum value fm. The variation of x with f—and, thus, with reference to (2), the initial packing fraction f0 and strain e—is illustrated in Fig. 3. Its value starts at x = 1 for 0 r f r fj and increases monotonically for fj r fo fm, exhibiting a vertical asymptote as f - fm. The thresholds fj = 0.84 and fm = 0.9 correspond to a random close-packing in two space dimensions21 and the maximum packing density for two-dimensional discs, respectively. Numerically, it is more convenient to use expressions that monotonically increase the compressive contribution without becoming unbounded. In this work, we use an expression with the
Fig. 2 (a) Schematic of a fractured particulate monolayer. The fractured zone is represented as a smooth transition zone across which the damage field d takes values between zero and unity. (b) The microstructure of the elastic domain O is captured by an additional packing fraction field f representing the solid to fluid ratio of the medium.
Soft Matter Paper
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5835
form 1.0 + f (1.0 + tanh((f 0.5(fm fj))/l)), with f being the factor of amplification and l being the length scale controlling the width of the regularization.
For the contribution to c associated with fracture, we choose a modified version
Qðd;rdÞ ¼ GðdÞ 2
(7)
of the expression proposed by Miehe et al.,17 with G defined by
GðdÞ ¼ G0 gs gf þ 1 gs
gf
gf
d
(8)
where gf and gs are the surface tension of the liquid layer and the surfactant, respectively, G0 4 0 is a constant that represents the fracture toughness of the undamaged monolayer, and l4 0 is a constant proportional to the characteristic thickness of a layer between damaged and undamaged material. This modifica- tion accounts for a reduction in the fracture toughness with increasing surfactant concentration. Accordingly, it is convenient to introduce
Gr ¼ Gð1Þ ¼ G0 gs gf
(9)
as a measure of the reduced fracture toughness of the monolayer. Finally, the surfactant contribution to c accounts for the
interplay between the monolayer and the surfactant (Fig. 4) and is assumed to be of the form
Fðe; dÞ ¼ F0 gf gsð ÞðtreÞd2
2 ; (10)
2.2 Governing equations
Following Silva et al.19 (but neglecting inertia), the governing equations of the phase-field model consist of the macroforce balance
div @cðe; d;rdÞ
8< : (12)
where b Z 0 represents the kinetic modulus that controls the rate at which cracks can propagate through the monolayer. The alternative on the right-hand side of (12) embodies the requirement that, consistent with experimental observations, cracks that form in the monolayer never heal. This requirement takes the form of the constraint
: d Z 0, and a reactive micro-
force pr is needed to ensure satisfaction of that constraint. In particular, pr vanishes for
: d 4 0 and is determined by the
left-hand side of (12) for : d = 0.
The kinetic modulus b incorporates two effects. First, it captures the capacity of the surfactant to penetrate the parti- culate monolayer and thereby generate damage, which is influenced by f. As f increase, the higher density of particles impede surfactant spreading. Second, it captures the resistance of the underlying liquid layer to rearrangements of the parti- cles, a resistance which is directly related to the viscosity mf
of the liquid comprising that layer. For the foregoing reasons, we assume that b increases with f and mf in accord with the relation
b ¼ b0 1 fm
1 f ; b0 ¼
; (13)
where C 4 0 is a dimensionless constant and b0 represents the kinetic modulus for maximal packing f = fm.
Since the particular choice of the kinetic modulus b defined in (13) is positive, damage only increases when the left-hand side of (12) is positive. Accordingly, the evolution eqn (12) for d becomes
1 fðeÞ b0 1 fmð Þ div
@cðe; d;rdÞ @rd
hhi ¼ 0; h 0;
h; h4 0;
( (15)
denotes its Macaulay bracket. With reference to the right-hand side of the definition (6) of
the free-energy density c and using (7)–(10) and (13) with reference to the parameters listed in Table 1, we find the
Fig. 3 Variation in the jamming factor x with the packing fraction f. The jamming factor tends towards N as f surpasses the jamming threshold fj = 0.84 and approaches the maximum packing fraction fm = 0.9.
Fig. 4 Representation of the surfactant force driving the system. A collection of particles, initially in equilibrium at a constant surface tension gf, are reached by the surfactant. The gradient between its lower surface tension gs and the surrounding media generates a Marangoni flow which pushes the particulate monolayer outward.
Paper Soft Matter
governing eqn (11) and (14) can be written as
div ð1 dÞ2@WþðeÞ @e
2 I
¼ 0
1 fðeÞ b0 1 fmð Þ G0l 1 1 gs
gf
d
G0d
respectively, and where I denotes the (two-dimensional) identity tensor. We consider (16) and (17) subject to the vanishing displacement condition
u = 0 (18)
on the boundary qO of O and, letting n denote a unit normal on qO, the natural boundary condition
rdn = 0. (19)
f(x,0) = f0(x), d(x,0) = d0(x), (20)
with f0 representing the initial distribution of particles and d0 the initial damage.
2.3 Model characterization: dimensionless number v
Introducing characteristic measures L and T of length and time, we define dimensionless counterparts x* and t* of x and t by
x ¼ x
T : (21)
rd ¼ Lrd; Dd ¼ L2Dd; and @d
@t ¼ T _d (22)
of the gradient, Laplacian, and partial time-derivative, dimen- sionless parameters,
G ¼ Gr
E ; L ¼ L
; g ¼ gs gf ; and l ¼ l
L ;
and dimensionless strain-energy densities
W ðeÞ ¼ 1
2ð1þ nÞ ej j 2 þ n ð1 nÞ2 treð Þ2; (24)
we arrive at dimensionless counterparts
div ð1 dÞ2@Wþ ðeÞ
@e þ F
G
g l 1 1 gð Þdð ÞDd G
g l
þ 2ð1 dÞLWþðeÞ LFðtreÞd ¼ @d
@t
(26)
of the governing eqn (16) and (17). While we have not undertaken an exhaustive suite of simula-
tions spanning the complete range of parameter space associated with the six dimensionless variables in (23), the extensive testing that we have conducted indicates that the fracture pattern is largely governed by the dimensionless driving force F* and the dimensionless fracture toughness G* of the monolayer. In particular, the number of fractures in the final configuration appears to be dictated by the ratio
w ¼ Fð Þ2
G : (27)
In Section 4.2, we discuss the threshold levels of the dimen- sionless parameter w that delineate different regimes, ranging from no fractures at all to a configuration with multiple branches.
The expression (27) for w can be further simplified by invoking the approximate scaling of Young’s modulus with the surface tension as derived by Vella et al.8 (i.e., E p gf/%rp). The expression for w can therefore be rewritten to isolate the influence of the surface tensions, yielding
w / F0rp 2 gf gsð Þ2
G0gs : (28)
Daniels et al.22 proposed a similar relationship between the elastic and Marangoni energies to explain fracture patterns in agar gels. In that work, the number of crack branches was found to scale linearly with the difference in surface tension. As described in Section 4, our experimental observations and
Table 1 List of parameters in the model
Parameter Description
fj Jamming packing fraction of monolayer fm Maximum packing fraction of monolayer %rp Average particle radius gf Surface tension of liquid layer gs Surface tension of surfactant G0 Fracture toughness of undamaged monolayer F0 Surfactant force constant E Young’s modulus of monolayer n Poisson’s ratio of monolayer mf Dynamic viscosity of monolayer b0 Kinetic parameter
Soft Matter Paper
model-based simulations suggest the system considered here exhibits a much stronger sensitivity to surface tension differences.
3 Experiments 3.1 Materials
All videos from experimental work were acquired using a Phantom V641 high-speed camera equipped with an AF Nikkor 50 mm f/1.8 D lens. Videos were saved using PCC software provided by Phantom. A 5 W white LED board used to illuminate the mixture was supplied by VANCO (series #33342). Dispersion of silica microballoons (SIG MFG) was performed using the Active Motif Q120AM (120 W, 20 kHz) ultrasonicator equipped with a CL18 3.2 mm probe. Oleic acid (surfactant) was supplied by WAKO Chemicals and the acetone used as a cosolvent was supplied by Nacalai Tesque. The surfactant was dropped using Terumo 2.5 mL syringes (SS-02SZ) equipped with a 25 G (0.50 16 mm) needle. Surface tension measurements of mixtures were obtained using the KSV NIMA LB Small trough (33473).
3.2 Methodology
Fig. 5 shows a schematic of the experimental setup. A clean low- form cylindrical glass vessel 17 cm in diameter and 9 cm in height filled to a height of 2.5 cm with Milli-Q water (or mixed with a cosolvent) was placed atop an LED light panel. Micro- balloons were carefully weighed to either 0.150 g or 0.250 g, introduced at the air–water interface, and dispersed by means of ultrasonication at 25% power to form a particulate mono- layer. Using a 2.5 mL syringe fitted with a 25 G metal tip, a single drop of surfactant (oleic acid with or without solvent) was introduced to the approximate center of the vessel at time t = 0 s. The water immiscible surfactant was observed to spread and push the microballoon particles radially outwards resulting in compaction of the particulate monolayers. Observations were captured using a Phantom high-speed camera at different frame rates. Data analysis was performed by first converting
raw data files (.cine) to multipage image files (.tif) and these files were subsequently analyzed using ImageJ. In cases where the surfactant or the bulk water phase was mixed with a solvent to reduce surface tension, acetone was added carefully and aliquots of the liquids were quickly transferred to an LB trough for surface tension measurements.
4 Results and discussion
We discretize the model described in Section 2 by recasting the evolution equations in a variational form and approximating the displacement and surfactant damage fields with finite element basis functions. Temporal integration is performed with a second-order accurate, unconditionally stable backward Euler algorithm (BDF2) with adaptive time stepping. The coupled system of nonlinear equations are solved with the Newton method in a staggered manner (i.e., alternating the updates for the displacement and damage fields).
4.1 Pattern comparison
Initial-boundary value problems are constructed over circular domains representative of the petri dishes used in the experi- ments. The experiments have a stochastic aspect corresponding to the initial particle locations. This effect is modeled at the continuum scale by using initial packing fields with spatial variability. Specifically, we numerically construct initial packing fraction fields over the circular domain that have a mean value of f = 0.795 with a 6% variation, corresponding to uniform random distributions with minimum and maximum packing fractions of 0.75 and 0.84, respectively. The initial drop of surfactant is approximated by a sufficiently small (on the order of l), empty circular region (f0 = 0) concentric with the center of the dish. Table 2 provides the reference values of the parameters used for the simulations in this section. All simulations were per- formed using MOOSE (Multiphysics Object-Oriented Simulation Environment), a finite-element framework developed primarily by Idaho National Laboratory.23
A first qualitative assessment of the model is shown in Fig. 6. In this example, the parameter w is adjusted to reproduce a three-branch pattern, which is a commonly observed configu- ration in the experiments for our given set of parameters. The result indicates that the model captures the general fracture process and reproduces the final three-branch configuration. Additionally, the cracks are sensitive to the packing fraction structure, which gives rise to more subtle features. These features include the particular angle in which the three- branch pattern is generated, the kinking and bending of cracks when gradients in the particle density are encountered, and the bifurcation of a crack when it encounters a high density region directly in line with its propagation. In this work, stochastic features of the experiments are approximated through random variations in the initial packing fraction field. The specific details of the numerically generated fracture patterns thus stem from the initial conditions and are sensitive to the particular random seed used in the simulations. As a result, our
Fig. 5 Setup: a clean low-form cylindrical vessel containing distilled water was placed on an LED tablet. A camera suspended vertically above records particulate monolayer compaction dynamics at the air–water interface when surfactant was introduced using a 25 G steel needle.
Paper Soft Matter
5838 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
simulation are designed to capture only the generic features of the experiments. Nevertheless, the model reproduces the number of cracks observed in the experiments and, at a qualitative level, the features of crack kinking and the for- mation of secondary branches.
We now provide a detailed description of a representative simulation in which a three-branch fracture pattern arises. Fig. 7 shows the results from a simulation of this common pattern, obtained with the parameters given in Table 2 and w = 0.0125 (see Table 3 for the correspondent dimensionless quantities). The columns in Fig. 7 show the evolution of three fields, namely the surfactant damage, the tensile energy, and the packing fraction, along with a comparison with snapshots of a matching experiment at comparable moments in time. The distribution of d closely matches the fractured portion of the particle raft in the experiment. The tensile strain-energy density field W+, which is representative of the processes that underpin crack propagation, transitions from an initially circular tensile distribution due to the surfactant force to one with broken symmetry depending on the available energy in the system. The tensile strain-energy density drives the three- branches until they eventually arrest due to an increase in the compressive energy. The field f depicts the evolution of the
initial packing fraction. As the fractured region expands, the monolayer contracts (to approximately f = 0.87) and the mobility of the surfactant in the fractured region (f = 0.3) increases.
4.2 Phase diagram for fracture pattern characterization
We now study whether the dimensionless parameter w can be used to construct a phase diagram for the fracture patterns that are produced. In so doing, we eliminate the random aspect of the previous simulations and consider initial packing fractions that are spatially uniform with f0 = 0.80. As illustrated in Fig. 8, the results allow us to identify threshold levels of w that delineate regions of a phase diagram for the fracture patterns obtained through simulations. The threshold w = 0.005 defines a limit below which fracture is not observed. As w is increased,
Table 2 Dimensional variables8,10
Property Description Value Unit
rd Radius of Petri dish 8.5 102 m f0 Range of initial packing fraction (7.5–8.4) 101 None fm Packing fraction of jammed monolayer 9.0 101 None G0 Fracture toughness of undamaged monolayer 5.7 (104–102) N m1
F0 Surfactant force constant 1.0 (101–102) m1
gf Surface tension of liquid layer (water) 7.2 102 N m1
gs Surface tension of surfactant (oleic acid) 4.0 102 N m1
mf Range of dynamic viscosity of liquid layer 1.0 (103–101) Pa s %rp Mean particle radius 1.0 104 m E Young’s modulus of monolayer 3.178 103 Pa n Poisson’s ratio of monolayer 1
ffiffiffi 3 p
None
b0 Kinetic parameter 1.0 Pa s l Characteristic length scale 1.0 102 m
Fig. 6 Comparison between the result of a simulation (left) and an experimental image (right) for a three-branch fracture pattern. The model captures the main traits of the final configuration, as represented by the number of fractures, and more detailed features corresponding to varia- tions in the distribution of particles, such as secondary fracturing, crack bending, and branching. The particular orientation of the fracture pattern shown in the simulation result is arbitrary, as it is sensitive to the random seed used to construct the initial packing field.
Fig. 7 Numerical fields obtained for a simulation of a three-branch fracture pattern. From left to right, each column shows the initial (top) and final (bottom) states of surfactant damage d, tensile strain-energy density W+ and f fields. The rightmost column shows the reference experiment, obtained by adding oleic acid to a 0.25 g mass monolayer of particles laying on a 18% salt saturation water solution.
Table 3 Dimensionless parameters
L* 8.5 102
b*0 2.7 (108–106) g* 5.6 101
l* 1.5 101
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5839
the fracture pattern starts to develop as a single crack, and then into three-branch shapes with further increases. We note that the three-branch configuration is quite ubiquitous in nature, as the presence of acute angles and the compression generated by the three cracks stabilizes the system. For larger values of w, a fourth and a fifth branch are gradually realized. Larger values of w are required to produce larger number of crack branches, resulting in unstable fracture patterns that cannot be clearly delineated in well defined regions of our phase diagram.
To validate the zones and thresholds predicted by the simulations, we first calibrated our model to experiments using oleic acid spreading over a pure water layer. The resulting fracture pattern is shown in the bottom of the center column of Fig. 9, which presents images of the final configurations for three different experiments, along with simulation results for comparable values of w. Numerical simulations indicate that such a pattern can be obtained using w = 0.018, which we label as wref. As indicated by expression (28), changing the surface tension difference corresponds to modulating the value of w in the model. Accordingly, we conducted two additional experi- ments. In the first of these experiments, the surface tension of the liquid layer was reduced from 72 to 56 mN m1 (colsolvent). In the second, the surface tension of the surfactant was reduced from 40 to 28 mN m1 (acetone). These changes in surface tension correspond approximately to modifying factors of 1/3 and 2, respectively, on the reference wref = 0018. We denote the corresponding levels of w as wdown and wup.
As illustrated in Fig. 9, simulations based on the values of wup and wdown yield fracture patterns that are remarkably consistent with the corresponding experiments. Simulations based on doubling w from wref to wup result in an increase in the number of crack branches, from three to five. Conversely, simulations based on a reduction of w from wref to wdown result in a decrease in the number of crack branches, from three to two. In each case, a qualitative match with the corresponding experiment is clear. As w is increased, the model presents a gradual change between patterns, from the formation of a
single crack to a configuration with multiple branches. The model nicely reproduces the central damage zone and multiple crack branches observed in experiments for the range (104–101) of w in this study.
This result is particularly important because the conven- tional view has been that the packing fraction distribution is the primary variable governing the morphological features of the crack pattern. Our results suggest that, although the randomness of the initial packing fraction can underpin features such as crack deflections and arrest, the branching and general patterns observed in the experiments correspond to a much more fundamental principle related to the interplay between the surfactant force and the fracture toughness. In other words, the difference in surface tension between the surfactant and underlying liquid are central to determining the crack patterns that emerge.
4.3 Temporal scaling
We now investigate the extent to which the model captures the speed of crack propagation. We examine the overall time required to fracture the particulate monolayer as well as the temporal scaling of the process. The snapshots in Fig. 10 show two fracture patterns that evolve at very different rates, but for which the final configurations are visually indistinguishable. For comparison purposes, we show our results using a dimensional version of the model (i.e., viscosity in Pa s). Previous experimental observa- tions have indicated that, for the same initial packing fractions, the geometry of the fracture pattern is insensitive to the viscosity of the underlying liquid layer.10 In the experiments, the liquid viscosity can be adjusted, for example, by increasing the ratio of glycerol in the glycerol–water mix based liquid layer, raising the viscosity from 103 Pa s to 101 Pa s. This is consistent with work showing that the rate at which the process occurs decreases by the same order of magnitude.10 We can capture this effect by adjusting the mobility implicitly defined in (13), which is propor- tional to the viscosity mf of the liquid layer.
Fig. 8 Phase diagram for fracture pattern characterization. Threshold levels of the dimensionless number w = (F*)2/G* delineate various patterns. In the region below w = 0.005 (zone i), the initially fractured region does not grow with time. As w increases, the fracture patterns gradually change from the single branch configuration (zone ii) to the three (zone iii), four (zone iv) and five (zone v) branch configurations.
Fig. 9 Comparison between fracture patterns from experiments (top) and simulations (bottom). In all cases, the final configuration is shown. Simula- tion results were obtained for values of the dimensionless parameter w (wdown = 0.006, wref = 0018, wup = 0.035), achieved by modulating the surface tension of the liquid layer and the surfactant used in each of the experiments, consistent with (28).
Paper Soft Matter
5840 | Soft Matter, 2017, 13, 5832--5841 This journal is©The Royal Society of Chemistry 2017
Additionally, experimental measurements indicate that the fracture proceeds in two stages. Crack growth rates that scale with the 3/4 and higher powers of t are observed in the early stages of the fracture process, followed by a gradual flattening as the crack advances and eventually arrests.10 Arrest may occur in response to the increase in packing around the perimeter of the domain, which gradually makes it more difficult to displace the particles of the monolayer. As shown in Fig. 11, our model qualitatively reproduces this behavior with time, in terms of the normalized fractured area A* (namely the ratio of the area of the surface of the liquid layer that is exposed to the area of O). As argued by Bandi et al.,11 if the radially receding front of R(t) of a particle raft scales according to R(t) B t3/4, then the particle free area proceeds as R(t)2 B (t3/4)2 B t3/2. The growth rates in the early stages are slightly below t3/2, which suggests the presence of the monolayer slightly retards the expansion of the surfactant from what would be expected over a pure liquid surface.24 Finally, we recognize a shift between the results from experiment and simulation, with the initial expansion process
being noticeably faster in the experiment. This discrepancy at early times could be related to the rapid expansion in the drop of surfactant, a process that is not captured by the model. The subtleties of the physics at early times are likely sensitive to the volume and height of the drop of surfactant at the time of injection. Questions related to these subtleties are a subject for future work.
Following the rapid initial expansion, the dynamics slow down considerably once the particulate rafts jam together. As shown in the inset to Fig. 11, the late stage dynamics exhibit approximately logarithmically slow relaxation over a decade in time. This kind of dynamics, also known as ‘‘creep’’ or ‘‘aging’’, is observed in a variety of phenomena including granular compaction,25 electron glasses,26 polymer relaxation,27 and super- conductors,28 and is now considered a hallmark of non- exponential, slow relaxation in amorphous media.29 Although the observation times in our model and experiments do not extend over several decades, they are consistent with the expected behavior of slow relaxation in frustrated, amorphous media.
5 Conclusions
We have presented a phase-field model that captures the main features of the surfactant-driven fracture of closely-packed particulate monolayers. The numerical simulations and experi- mental results indicate that there is a competition between the spreading force of the surfactant and the fracture toughness of the monolayer that determines the number of fractures. The model gives rise to a dimensionless parameter that can be written in terms of the surface tensions, and which allows for a straightforward comparison with experimental conditions. Experiments were conducted to validate the model and delineate the regimes of fracture pattern as a function of that parameter.
Our model rests on a number of simplifying assumptions. In particular, we used the damage to the particulate monolayer as an indicator function for the surfactant concentration. As such, our model precludes investigations of the extent to which the surfactant can penetrate into the particulate monolayer network, ahead of (or behind) the crack front. We have also assumed that the strains exhibited by the monolayer are sufficiently small to justify modeling it as a linearly elastic solid. This assumption ceases to be reasonable if the surface tension difference is increased above some critical threshold. Moreover, its applic- ability is limited to systems in which the initial packing is low enough to allow the particulate layer to exhibit both loosely- packed and jammed behaviors as the surfactant spreads.11 In either case, transitions between linear and nonlinear response, or between fluid and solid behavior, are worthy of further study. Finally, the interplay between liquid and solid constituents in the monolayer may endow it with a viscoelastic behavior that could also be explored. As the rate of surfactant transport is a function of the Marangoni stress, modulating the surface tension difference might also be used to explore any rate depen- dence in the monolayer. Future work will focus on enhancing the
Fig. 10 Identical fracture patterns are produced at different speeds due to changes in the viscosity mf of the underlying liquid layer.
Fig. 11 Comparison between simulation and experiment for the temporal scaling for the normalized fractured area A* (the ratio of the portion of the area of the liquid layer that is exposed to the area of O), indicating two regimes of crack evolution. For reference, a t3/2 curve is provided (red dashed line). The inset presents a zoom of the latter stages using semi-log axes, revealing a slow logarithmic relaxation of the fracture pattern.
Soft Matter Paper
This journal is©The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 5832--5841 | 5841
model to incorporate these effects and to allow for detailed exploration of their consequences and significance.
Acknowledgements
Christian Peco, Yingjie Liu, and John E. Dolbow are grateful for the support of NSF grant CMMI-1537306, to Duke University. Wei Chen, M. M. Bandi, and Eliot Fried gratefully acknowledge support from the Okinawa Institute of Science and Technology Graduate University with subsidy funding from the Cabinet Office, Government of Japan. John E. Dolbow would also like to acknowledge support from the Okinawa Institute of Science and Technology during a sabbatical stay in late 2014, during which time this work began.
References
1 A. Varshney, A. Sane, S. Ghosh and S. Bhattacharya, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2012, 86, 31402.
2 S. U. Pickering, J. Chem. Soc., Trans., 1907, 91, 2001–2021. 3 B. P. Binks and P. D. I. Fletcher, Langmuir, 2001, 17,
4708–4710. 4 P. Cicuta, E. J. Stancik and G. G. Fuller, Phys. Rev. Lett., 2003,
90, 236101. 5 P. Aussillous and D. Quere, Nature, 2001, 411, 924–927. 6 N. Tsapis, D. Bennett, B. Jackson, D. A. Weitz and D. A. Edwards,
Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 12001–12005. 7 A. B. Subramaniam, M. Abkarian and H. A. Stone, Nat. Mater.,
2005, 4, 553–556. 8 D. Vella, P. Aussillous and L. Mahadevan, Europhys. Lett.,
2004, 68, 212–218. 9 L. Scriven and C. Sterling, Nature, 1960, 187, 186–188.
10 D. Vella, H. Y. Kim, P. Aussillous and L. Mahadevan, Phys. Rev. Lett., 2006, 96, 1–4.
11 M. M. Bandi, T. Tallinen and L. Mahadevan, Europhys. Lett., 2011, 96, 36008.
12 P. Cicuta and D. Vella, Phys. Rev. Lett., 2009, 102, 138302. 13 C. Planchette, E. Lorenceau and A.-L. Biance, Soft Matter,
2012, 8, 2444–2451. 14 C. W. MacMinn, E. R. Dufresne and J. S. Wettlaufer,
Phys. Rev. X, 2015, 5, 11020. 15 C. Spandagos, T. B. Goudoulas, P. F. Luckham and O. K.
Matar, Langmuir, 2012, 28, 7197–7211. 16 E. H. Yoffe, Philos. Mag., 1951, 42, 739–750. 17 C. Miehe, M. Hofacker and F. Welschinger, Comput. Methods
Appl. Mech. Eng., 2010, 199, 2765–2778. 18 M. J. Borden, C. V. Verhoosel, M. A. Scott, T. J. R. Hughes
and C. M. Landis, Comput. Methods Appl. Mech. Eng., 2012, 217-220, 77–95.
19 M. N. da Silva Jr., F. P. Duda and E. Fried, J. Mech. Phys. Solids, 2013, 61, 2178–2195.
20 C. Miehe, Comput. Struct., 1998, 66, 37–43. 21 C. S. O’Hern, L. E. Silbert, A. J. Liu and S. R. Nagel, Phys. Rev.
E: Stat., Nonlinear, Soft Matter Phys., 2003, 68, 11306. 22 K. E. Daniels, S. Mukhopadhyay, P. J. Houseworth and
R. P. Behringer, Phys. Rev. Lett., 2007, 99, 124501. 23 D. Gaston, C. Newman, G. Hansen and D. Lebrun-Grandie,
Nucl. Eng. Des., 2009, 239, 1768–1778. 24 O. E. Jensen, J. Fluid Mech., 1995, 293, 349–378. 25 E. Ben-Naim, J. B. Knight, E. R. Nowak, H. M. Jaeger and
S. R. Nagel, Phys. D, 1998, 123, 380–385. 26 A. Amir, Y. Oreg and Y. Imry, Phys. Rev. Lett., 2009,
103, 126403. 27 J. D. Ferry, Viscoelastic Properties of Polymers, John Wiley and
Sons, 1980. 28 A. Gurevich and H. Kupfer, Phys. Rev. B: Condens. Matter
Mater. Phys., 1993, 48, 6477–6487. 29 A. Amir, Y. Oreg and Y. Imry, Proc. Natl. Acad. Sci. U. S. A.,
2012, 109, 1850–1855.