Numerical simulation methods for the Rouse model in flow

Michael P. Howard, Scott T. Milner

Research output: Contribution to journalArticlepeer-review

2 Scopus citations


Simulation of the Rouse model in flow underlies a great variety of numerical investigations of polymer dynamics, in both entangled melts and solutions and in dilute solution. Typically a simple explicit stochastic Euler method is used to evolve the Rouse model. Here we compare this approach to an operator splitting method, which splits the evolution operator into stochastic linear and deterministic nonlinear parts and takes advantage of an analytical solution for the linear Rouse model in terms of the noise history. We show that this splitting method has second-order weak convergence, whereas the Euler method has only first-order weak convergence. Furthermore, the splitting method is unconditionally stable, in contrast to the limited stability range of the Euler method. Similar splitting methods are applicable to a broad class of problems in stochastic dynamics in which noise competes with ordering and flow to determine steady-state order parameter structures.

Original languageEnglish (US)
Article number051804
JournalPhysical Review E - Statistical, Nonlinear, and Soft Matter Physics
Issue number5
StatePublished - Nov 28 2011

All Science Journal Classification (ASJC) codes

  • Statistical and Nonlinear Physics
  • Statistics and Probability
  • Condensed Matter Physics


Dive into the research topics of 'Numerical simulation methods for the Rouse model in flow'. Together they form a unique fingerprint.

Cite this