Modeling Polymer Rheology: Constitutive Equations and Simulation Approaches

Modeling Polymer Rheology: Constitutive Equations and Simulation Approaches

Overview

Modeling polymer rheology links microscopic polymer physics to macroscopic flow behavior using constitutive equations and numerical simulation. Models predict stress, viscosity, normal stresses, and time-dependent responses (relaxation, creep, start-up) needed for processing, formulation, and product design.

Key constitutive equation families

• Newtonian and generalized Newtonian

  • Newtonian: constant viscosity, valid for low-molecular-weight fluids.
  • Generalized Newtonian (e.g., Carreau, Cross, Power-law): shear-rate–dependent viscosity, no elastic memory — good for steady shear flows.

• Linear viscoelastic models

  • Maxwell model (single or multi-mode): captures stress relaxation and linear viscoelastic spectra.
  • Kelvin–Voigt: better for creep-dominated responses.
  • Generalized linear models represented by relaxation spectrum G(t) or complex modulus G*(ω).

• Nonlinear viscoelastic / integral-type models

  • K-BKZ (Kaye–Bernstein–Kearsley–Zapas): integral constitutive model for large deformations with strain-history dependence.
  • Lodge rubber-like liquid: useful for moderately entangled polymer melts/solutions.

• Differential-type molecularly inspired models

  • Oldroyd-B and upper-convected Maxwell (UCM): simple viscoelastic models including flow kinematics; UCM predicts normal stresses but has unbounded extensional viscosity.
  • Giesekus: adds anisotropic drag to produce shear-thinning and finite extensional viscosity.
  • Phan–Thien–Tanner (PTT): models strain-dependent viscosity and finite extensibility.
  • FENE-P / FENE-CR: finitely extensible nonlinear elastic dumbbell models, useful for dilute solutions.

• Tube/entanglement-based models (for concentrated/entangled polymers)

  • Doi–Edwards tube model and its extensions (convective constraint release, chain stretch): capture reptation, stress relaxation, and nonlinear shear/thinning for high-molecular-weight melts.
  • Rolie-Poly: a simplified differential form combining reptation, stretch, and CCR — widely used in simulations of entangled polymers.

Important model features to choose

• Shear vs extensional behavior: Some models predict extensional viscosity poorly (UCM diverges); choose models with finite extensibility or CCR modifications for extensional flows.
• Memory and timescales: Multi-mode spectra or multiple relaxation times are needed to capture broad molecular weight distributions.
• Nonlinearities: Large-deformation flows require models that include strain-dependent relaxation, stretching, and orientation.
• Thermorheological complexity: Include temperature dependence (time–temperature superposition, WLF/Arrhenius) if processing spans temperatures.
• Constitutive stability: Avoid models that produce numerical instabilities in strong flows without appropriate stabilization.

Numerical simulation approaches

• Finite Element Method (FEM): Widely used for complex geometries (extrusion dies, molds). Requires stabilization techniques (DEVSS, SUPG) and often log-conformation methods to improve numerical stability at high Weissenberg numbers.
• Finite Volume Method (FVM): Popular in CFD codes for flows in process equipment; well-suited to conservation-law formulations.
• Spectral / spectral-element methods: High accuracy for smooth problems, limited to simpler geometries.
• Particle-based methods: Brownian dynamics, multi-particle collision dynamics, dissipative particle dynamics useful for mesoscale or microstructure-resolving simulations.
• Hybrid multiscale methods: Couple continuum constitutive models to molecular dynamics or kinetic simulations (e.g., CONNFFESSIT, heterogeneous multiscale modeling) for parameter estimation or bridging scales.

Parameterization and validation

• Rheometry data: Use steady shear, oscillatory shear (G’, G”), stress relaxation, creep, and extensional rheometry to fit model parameters.
• Multi-mode fitting: Represent relaxation spectra with discrete modes; fit molecular-weight-dependent timescales.
• Microstructural input: For molecular models, input chain length, entanglement molecular weight, and friction parameters from experiments or molecular simulations.
• Benchmarking: Validate against processing-relevant flows (capillary extrusion, start-up shear, step strains, contraction flows) and compare predicted pressure drops, die swell, and stress growth.

Practical tips for simulations

• Start with simpler models (generalized Newtonian or single-mode Maxwell) to get baseline behavior, then increase complexity.
• Use log-conformation transformation for high Weissenberg number problems.
• Apply mesh refinement near boundaries, steep gradients, and geometric singularities.
• Implement consistent boundary conditions for stress/extra-stress and account for slip if present.
• Calibrate on small-strain linear viscoelastic data before fitting nonlinear parameters.

Example workflow (concise)

1. Characterize rheology: steady shear, oscillatory, extensional tests.
2. Choose model family (e.g., Rolie-Poly for entangled melts).
3. Fit parameters (multi-mode spectrum, relaxation times, CCR factors).
4. Implement model in FEM/CFD solver with stabilization and log-conformation.
5. Run benchmark flows, refine mesh and parameters.
6. Validate against processing experiments and iterate.

Further reading (select topics)

• Tube models and Rolie-Poly for entangled polymers
• Log-conformation numerical stabilization
• FENE-P and dumbbell models for dilute solutions

If you want, I can: provide a ready-to-run Rolie-Poly implementation for your solver of choice, fit a multi-mode Maxwell spectrum from your rheometry data, or generate simulation setup suggestions for a specific processing geometry.