Review on the Brownian Dynamics Simulation of BeadRodSpring Models Encountered in Computational Rheology
Type
Articles dans des revues avec comité de lectureDate
2012Journal
Archives of Computational Methods in EngineeringAbstract
Kinetic theory is a mathematical framework intended to relate directly the most relevant characteristics of the molecular structure to the rheological behavior of the bulk system. In other words, kinetic theory is a microtomacro approach for solving the flow of complex fluids that circumvents the use of closure relations and offers a better physical description of the phenomena involved in the flow processes. Cornerstone models in kinetic theory employ beads, rods and springs for mimicking the molecular structure of the complex fluid. The generalized beadrodspring chain includes the most basic models in kinetic theory: the freely jointed beadspring chain and the freelyjointed beadrod chain. Configuration of simple coarsegrained models can be represented by an equivalent FokkerPlanck (FP) diffusion equation, which describes the evolution of the configuration distribution function in the physical and configurational spaces. FP equation can be a complex mathematical object, given its multidimensionality, and solving it explicitly can become a difficult task. Even more, in some cases, obtaining an equivalent FP equation is not possible given the complexity of the coarsegrained molecular model. Brownian dynamics can be employed as an alternative extensive numerical method for approaching the configuration distribution function of a given kinetictheory model that avoid obtaining and/or resolving explicitly an equivalent FP equation. The validity of this discrete approach is based on the mathematical equivalence between a continuous diffusion equation and a stochastic differential equation as demonstrated by Itô in the 1940s. This paper presents a review of the fundamental issues in the BD simulation of the linear viscoelastic behavior of beadrodspring coarse grained models in dilute solution. In the first part of this work, the BD numerical technique is introduced. An overview of the mathematical framework of the BD and a review of the scope of applications are presented. Subsequently, the links between the rheology of complex fluids, the kinetic theory and the BD technique are established at the light of the stochastic nature of the beadrodspring models. Finally, the pertinence of the present stateoftheart review is explained in terms of the increasing interest for the stochastic microtomacro approaches for solving complex fluids problems. In the second part of this paper, a detailed description of the BD algorithm used for simulating a smallamplitude oscillatory deformation test is given. Dynamic properties are employed throughout this work to characterise the linear viscoelastic behavior of beadrodspring models in dilute solution. In the third and fourth part of this article, an extensive discussion about the main issues of a BD simulation in linear viscoelasticity of diluted suspensions is tackled at the light of the classical multibeadspring chain model and the multibeadrod chain model, respectively. Kinematic formulations, integration schemes and expressions to calculate the stress tensor are revised for several classical models: Rouse and Zimm theories in the case of multibeadspring chains, and Kramers chain and semiflexible filaments in the case of multibeadrod chains. The implemented BD technique is, on the one hand, validated in front of the analytical or exact numerical solutions known of the equivalent FP equations for those classic kinetic theory models; and, on the other hand, is controlset thanks to the analysis of the main numerical issues involved in a BD simulation. Finally, the review paper is closed by some concluding remarks.
Files in this item
Related items
Showing items related by title, author, creator and subject.

Shearstrain step response in linear regime of dilute suspensions of naturally bent carbon nanotubes CRUZ, Camilo; ILLOUL, Lounès; CHINESTA, Francisco; REGNIER, Gilles (Wiley, 2012)Impressive enhancements of the storage modulus have been documented when low volume fractions of single wall carbon nanotubes (SWNTs) are added to a Newtonian solvent for obtaining dilute suspensions. The intrinsic bending ...

DIANI, Julie; FREDY, Carole; GILORMINI, Pierre; MERCKEL, Yannick; REGNIER, Gilles; ROUSSEAU, Ingrid (Elsevier, 2011)A torsion device was designed and built for testing the shape fixity and shape recovery of shape memory polymers at large deformation. A simple thermal chamber was used to regulate the thermal environment during testing ...

GILORMINI, Pierre; CHEVALIER, Luc; REGNIER, Gilles (Wiley, 2010)In order to simulate the thermoforming of a transparency, constitutive equations are proposed for the nonlinear viscoelastic behaviour of poly(methyl methacrylate) near glass transition temperature, which include large ...

TEYSSEDRE, Hubert; GILORMINI, Pierre; LANDIS, Stefan; REGNIER, Gilles (Elsevier, 2013)Nanoimprint lithography is an efficient way to reproduce nanostructures down to 20 nanometers in submicrometer polymeric films. To optimize this process, simulation using a Newtonian behavior is a cheap and efficient way ...

WALCZAK, M; CIESIELSKI, W; GALESKI, A; POTRZEBOWSKI, M.J.; REGNIER, Gilles; HILTNER, A; BAER, Eric (Wiley, 2012)Multilayered film polycarbonate/polystyrene (PC/PS) comprising 257 layers with total thickness 125 m was made by coextrusion process. The nominal thickness of PC layers was 680 nm, and the nominal thickness of PS layers ...