1 / 29

J. Stegen + , J. Billen ° , M. Wilson ° , A.R.C. Baljon ° . A.V. Lyulin +

Structural origin of non-Newtonian rheology Computer simulations on a solution of telechelic associating polymers. J. Stegen + , J. Billen ° , M. Wilson ° , A.R.C. Baljon ° . A.V. Lyulin + + Eindhoven University of Technology (The Netherlands) ° San Diego State University (USA).

wade-carson
Download Presentation

J. Stegen + , J. Billen ° , M. Wilson ° , A.R.C. Baljon ° . A.V. Lyulin +

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Structural origin of non-Newtonian rheologyComputer simulations on a solution of telechelic associating polymers J. Stegen+, J. Billen°, M. Wilson °, A.R.C. Baljon °. A.V. Lyulin+ + Eindhoven University of Technology (The Netherlands) ° San Diego State University (USA)

  2. Introduction

  3. Polymeric gels Reversible junctions between end groups (telechelic associating polymers) Concentration Sol Gel Temperature

  4. Constitutive relation for gel • Stress • Shear rate • Viscosity • Constitutive relation for gel • Regime where stress decreases with increasing shear due to shear induced structure: • decrease in number of elastic junctions • increased orientation in shear direction stress shear rate

  5. Hybrid MD/MC simulation of a polymeric gel

  6. Molecular dynamics simulation Molecular dynamics: Grest-Kremer bead-spring model Equations of motion: (Langevin equation, coupling to heat bath through fluctuation dissipation theorem)

  7. U [e] Distance [s] Bead-spring model [K. Kremer and G. S. Krest. J. Chem. Phys 1990] Attraction beads in chain 1s Repulsion all beads

  8. Ubond U [e] Unobond Distance [s] Associating polymer [A. Baljon et al., J. Chem. Phys., 044907 2007] • Junction between end groups : LJ + FENE + Association energy

  9. P=1 form P<1 possible form Dynamics of associating polymer • Monte Carlo: attempt to form or destroy junction ? D U [e] Uassoc=-22 Distance [s]

  10. Simulation details • 1000 polymeric chains, 8 beads/chain • Units: s (length), e (energy & temperature), m (mass), t=s(m/e)1/2 (time); • Box size: (23.5 x 20.5 x 27.4) s3 with: • periodic boundary conditions in x,y direction. • Fixed walls in z-direction • Average volume density in system: 0.32 • NVT simulation

  11. Shearing the system Move wall with constant shear rate. Some chains grafted to wall to minimise wall slip (50 per wall) moving wall fixed wall

  12. Nomenclature Bead (8 per chain) • Chain bead (6 per chain, white/gray) • End group (2 per chain) • Dangler (blue) • Loop (orange) • Aggregate (red & orange) Network structure of 4 chains Single chain

  13. Structural properties in equilibrium

  14. Structural properties in mechanical equilibrium I

  15. Structural properties in mechanical equilibrium II

  16. Structural properties in mechanical equilibrium II

  17. Structural properties in mechanical equilibrium III T=1.0

  18. Structural properties in mechanical equilibrium III T=0.55

  19. Structural properties in mechanical equilibrium III T=0.35

  20. Structural properties in mechanical equilibrium IV: Conclusions • Aggregates increase in size with decreasing temperature • Gel network immobile, macroscopic lifetime • Spatial ordering of aggregates observed in gel phase • Boundary effects visible at all temperatures, induces structure and ordering at lower temperature

  21. Shear Banding

  22. Shear banding: theory Instable region in constitutive relation (striped) Stable configuration through two shear bands coexisting at a stress σ Lever rule: Plateau in shear-stress curve Difference in mesoscopicstructure between bands

  23. Shear banding: force and velocity profile Simulation details: T=0.35ε wall velocity 0.01 σ/τ shear rate 3.6*10-4τ-1 total wall displacement ~700 σ

  24. Shear banding: aggregate size distribution • More small and large aggregates in shear banding state • Large aggregates strong influence on velocity profile?

  25. xx zz xz Shear banding: orientation function Orientation in xx-direction, xz-direction and perpendicular to zz-direction: effects of applied shear on chains decrease No significant differences between shear bands

  26. Shear banding: spatial distribution High shear band very small (~5σ), too small to contain mesoscopic structure? Fluctuations in density of ~10% at bottom of high shear band. No stationary flow but hopping like behaviour of end groups at interface? Shear direction

  27. Conclusion • Shear bands in velocity profile observed. • High shear band too small to accommodate a mesoscopic structure different from the low shear band.No significant differences in structure observed between bands. • More large aggregates in a sheared system, these could be responsible for the observed shear banding. • Fluctuations in end-group density at interface, no steady flow. • Validity of lever rule has not been checked. Uncertain if observed shear banding corresponds to the shear banding observed in experiment.

  28. Other work… Jammed system at constant stress & fluctuation relation • Elastic behaviour visible • Two types of behaviour observed in time • Deviations from fluctuation relation observed

  29. Questions?

More Related