Publications

Results 201–225 of 265

Search results

Jump to search filters

Substructured multibody molecular dynamics

Crozier, Paul C.; Grest, Gary S.; Ismail, Ahmed I.; Lehoucq, Richard B.; Plimpton, Steven J.; Stevens, Mark J.

We have enhanced our parallel molecular dynamics (MD) simulation software LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, lammps.sandia.gov) to include many new features for accelerated simulation including articulated rigid body dynamics via coupling to the Rensselaer Polytechnic Institute code POEMS (Parallelizable Open-source Efficient Multibody Software). We use new features of the LAMMPS software package to investigate rhodopsin photoisomerization, and water model surface tension and capillary waves at the vapor-liquid interface. Finally, we motivate the recipes of MD for practitioners and researchers in numerical analysis and computational mechanics.

More Details

Report on ASC project degradation of organic materials

Lo, Chi S.; Thompson, Aidan P.; Grest, Gary S.

Using molecular dynamics simulations, a constitutive model for the chemical aging of polymer networks was developed. This model incorporates the effects on the stress from the chemical crosslinks and the physical entanglements. The independent network hypothesis has been modified to account for the stress transfer between networks due to crosslinking and scission in strained states. This model was implemented in the finite element code Adagio and validated through comparison with experiment. Stress relaxation data was used to deduce crosslinking history and the resulting history was used to predict permanent set. The permanent set predictions agree quantitatively with experiment.

More Details

MD simulations of chemically reacting networks : analysis of permanent set

Curro, John G.; Budzien, Joanne L.; Grest, Gary S.; Thompson, Aidan P.

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.

More Details

Constitutive models for rubber networks undergoing simultaneous crosslinking and scission

Budzien, Joanne L.; Lo, Chi S.; Curro, John G.; Thompson, Aidan P.; Grest, Gary S.

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.

More Details

Elucidating the mysteries of wetting

Brooks, Carlton F.; Emerson, John A.; Grest, Gary S.; Grillet, Anne M.; Sackinger, Philip A.; Ash, Benjamin J.; Webb, Edmund B.; Gorby, Allen D.; Bourdon, Christopher B.

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.

More Details

Surface wetting of liquid nanodroplets: Droplet-size effects

Physical Review Letters

Heine, David R.; Grest, Gary S.; Webb, Edmund B.

The spreading of liquid nanodroplets of different initial radii R0 is studied using molecular dynamics simulation. Results for two distinct systems, Pb on Cu(111), which is nonwetting, and a coarse-grained polymer model, which wets the surface, are presented for Pb droplets ranging in size from ∼55000 to 220000 atoms and polymer droplets ranging in size from ∼200000 to 780000 monomers. In both cases, a precursor foot precedes the spreading of the main droplet. This precursor foot spreads as rf2(t)=2Defft with an effective diffusion constant that exhibits a droplet-size dependence Deff∼R01/2. The radius of the main droplet rb(t)∼R04/5 is in agreement with kinetic models for the cylindrical geometry studied. © 2005 The American Physical Society.

More Details
Results 201–225 of 265
Results 201–225 of 265