The Independent Network Model (INM) has proven to be a useful tool for understanding the development of permanent set in strained elastomers. Our previous work showed the applicability of the INM to our simulations of polymer systems crosslinking in strained states. This study looks at the INM applied to theoretical models incorporating entanglement effects, including Flory's constrained junction model and more recent tube models. The effect of entanglements has been treated as a separate network formed at gelation, with additional curing treated as traditional phantom contributions. Theoretical predictions are compared with large-scale molecular dynamics simulations.
Constitutive models for chemically reacting networks are formulated based on a generalization of the independent network hypothesis. These models account for the coupling between chemical reaction and strain histories, and have been tested by comparison with microscopic molecular dynamics simulations. An essential feature of these models is the introduction of stress transfer functions that describe the interdependence between crosslinks formed and broken at various strains. Efforts are underway to implement these constitutive models into the finite element code Adagio. Preliminary results are shown that illustrate the effects of changing crosslinking and scission rates and history.
Nearly every manufacturing and many technologies central to Sandia's business involve physical processes controlled by interfacial wetting. Interfacial forces, e.g. conjoining/disjoining pressure, electrostatics, and capillary condensation, are ubiquitous and can surpass and even dominate bulk inertial or viscous effects on a continuum level. Moreover, the statics and dynamics of three-phase contact lines exhibit a wide range of complex behavior, such as contact angle hysteresis due to surface roughness, surface reaction, or compositional heterogeneities. These thermodynamically and kinetically driven interactions are essential to the development of new materials and processes. A detailed understanding was developed for the factors controlling wettability in multicomponent systems from computational modeling tools, and experimental diagnostics for systems, and processes dominated by interfacial effects. Wettability probed by dynamic advancing and receding contact angle measurements, ellipsometry, and direct determination of the capillary and disjoining forces. Molecular scale experiments determined the relationships between the fundamental interactions between molecular species and with the substrate. Atomistic simulations studied the equilibrium concentration profiles near the solid and vapor interfaces and tested the basic assumptions used in the continuum approaches. These simulations provide guidance in developing constitutive equations, which more accurately take into account the effects of surface induced phase separation and concentration gradients near the three-phase contact line. The development of these accurate models for dynamic multicomponent wetting allows improvement in science based engineering of manufacturing processes previously developed through costly trial and error by varying material formulation and geometry modification.
Slow, dense granular flows often exhibit thin, localized regions of particle motion, called shear bands, separating largely solid-like regions. Recent experiments using a split-bottom Couette cell found that the width of the shear zone grew as the pack height increased and the azimuthal velocities when rescaled fall on a universal curve regardless of the particle properties. Here we present large-scale Discrete Element simulations of a similar system for packs of varying height up to 180,000 monodisperse spheres. The onset and evolution of granular shear flow is investigated as a function of height. We find a transition in the nature of the shear as a characteristic height is exceeded. Below this height there is a central quasi-solid core; above this height we observe the onset of additional axial shear associated with a torsional failure mode of the inner core. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height while the axial width remains narrow and fixed.
The evolution of granular shear flow is investigated as a function of height in a split-bottom Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations, we find a transition in the nature of the shear as a characteristic height H* is exceeded. Below H* there is a central stationary core; above H* we observe the onset of additional axial shear associated with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height, while the axial width remains narrow and fixed.
Molecular-dynamics simulations are used to sample the single-chain form factor of labelled sub-chains in model polymer networks under elongational strain. We observe very similar results for randomly cross-linked and for randomly end-linked networks with the same average strand length and see no indication of lozenge-like scattering patterns reported for some experimental systems. Our data analysis shows that a recent variant of the tube model quantitatively describes scattering in the Guinier regime as well as the macroscopic elastic properties. The observed failure of the theory outside the Guinier regime is shown to be due to non-Gaussian pair-distance distributions.
Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state flow rate is reduced, the force distribution P(f) changes only slightly, while there is a large change in the impulse distribution P(i). In both cases, the distributions show an increase in small forces or impulses as the systems approach jamming, the opposite of that seen in previous Lennard-Jones simulations. This occurs dynamically as well for a hopper that transitions from a flowing to a jammed state over time. The final jammed P(f) is quite distinct from a poured packing P(f) in the same geometry. The change in P(i) is a much stronger indicator of the approach to jamming. The formation of a peak or plateau in P(f) at the average force is not a general feature of the approach to jamming.
In this study, we perform molecular dynamics simulations of adhesive contact and friction between alkylsilane Si(OH){sub 3}(CX{sub 2}){sub 10}CX{sub 3} and alkoxylsilane Si(OH){sub 2}(CX{sub 2}){sub 10}CX{sub 3} (where X = H or F) self-assembled monolayers (SAMs) on an amorphous silica substrate. The alkylsilane SAMs are primarily hydrogen-bonded or physisorbed to the surface. The alkoxylsilane SAMs are covalently bonded or chemisorbed to the surface. Previously, we studied the chemisorbed systems. In this work, we study the physisorbed systems and compare the tribological properties with the chemisorbed systems. Furthermore, we examine how water at the interface of the SAMs and substrate affects the tribological properties of the physisorbed systems. When less than a third of a monolayer is present, very little difference in the microscopic friction coefficient {mu} or shear stresses is observed. For increasing amounts of water, the values of {mu} and the shear stresses decrease; this effect is somewhat more pronounced for fluorocarbon alkylsilane SAMs than for the hydrocarbon SAMs. The observed decrease in friction is a consequence of a slip plane that occurs in the water as the amount of water is increased. We studied the frictional behavior using relative shear velocities ranging from v = 2 cm/s to 2 m/s. Similar to previously reported results for alkoxylsilane SAMs, the values of the measured stress and {mu} for the alkylsilane SAM systems decrease monotonically with v.
We present extensive simulations modeling the casting of multiblock polymer films by evaporation. The domain structure of the resulting film is strongly affected by varying the relative stiffness of the coblocks. The morphology changes from a bicontinuous lamellar phase when both blocks are flexible to a small-scale phase-separated phase with isolated domains as the stiffness of one of the blocks increases. As the relative stiffness of the blocks changes, the rate of evaporation, interfacial width, and morphology of the system changes. The findings can be used to tailor membrane morphology of interest to fuel-cell applications where the morphology is important for proton conduction.
Large-scale three dimensional Discrete Element simulations of granular flow in a modified split-bottom Couette cell for packs of up to 180,000 mono-disperse spheres are presented and compared with experiments. We find that the velocity profiles collapse onto a universal curve not only at the surface but also in the bulk of the pack until slip between layers becomes significant. In agreement with experiment, we find similar relations between the cell geometry and parameters involved in rescaling the velocities at the surface and in the bulk. Likewise, a change in the shape of the shear zone is observed as predicted for tall packs once the center of the shear zone is correctly defined; although the transition does not appear to be first order. Finally, the effect of cohesion is considered as a means to test the theoretical predictions.
Similar to entangled ropes, polymer chains cannot slide through each other. These topological constraints, the so-called entanglements, dominate the viscoelastic behavior of high-molecular-weight polymeric liquids. Tube models of polymer dynamics and rheology are based on the idea that entanglements confine a chain to small fluctuations around a primitive path which follows the coarse-grained chain contour. To establish the microscopic foundation for these highly successful phenomenological models, we have recently introduced a method for identifying the primitive path mesh that characterizes the microscopic topological state of computer-generated conformations of long-chain polymer melts and solutions. Here we give a more detailed account of the algorithm and discuss several key aspects of the analysis that are pertinent for its successful use in analyzing the topology of the polymer configurations. We also present a slight modification of the algorithm that preserves the previously neglected self-entanglements and allows us to distinguish between local self-knots and entanglements between distant sections of the same chain. Our results indicate that the latter make a negligible contribution to the tube and that the contour length between local self-knots, N{sub 1k} is significantly larger than the entanglement length N{sub e}.
