26.11.2012 Views

Non-Newtonian turbulence: viscoelastic fluids and binary mixtures.

Non-Newtonian turbulence: viscoelastic fluids and binary mixtures.

Non-Newtonian turbulence: viscoelastic fluids and binary mixtures.

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2.1. Polymer dynamics in <strong>fluids</strong> 37<br />

2.1.3 Dumbbell model<br />

A model which is very often used for its simplicity is the dumbbell model [41],<br />

where the complex structure of the polymer molecule is replaced by a couple of<br />

beads of negligible mass, connected by a spring. Such a spring has the same<br />

properties of an entire freely jointed chain (see sec. 2.1.1) of Hook modulus<br />

H = 3kBT<br />

2Rmaxb . The evolution of the dumbbell end-to-end vector R = x2 − x1<br />

is determined by different contributions: the hydrodynamic drag force acting on<br />

the polymer, thermal noise, <strong>and</strong> the elastic recalling force of the spring. In a homogeneous<br />

flow, the equation of motion for R then is:<br />

˙R = − H<br />

R + Bξ (2.14)<br />

β<br />

where β is the friction coefficient <strong>and</strong> ξ is a zero mean Brownian process with<br />

correlation 〈ξi(t)ξj(t ′ )〉 = δijδ(t − t ′ ). The relaxation time is clearly introduced<br />

as:<br />

τ ≡ β<br />

H<br />

(2.15)<br />

The dependence of the constant B on the gyration radius R0 <strong>and</strong> the relaxation<br />

time τ can be derived as follows. The formal solution of eq. (2.14) allows to<br />

estimate the long time behaviour of the square mean elongation as 〈R2 〉 ∼ τB2<br />

. At 2<br />

is comparable<br />

equilibrium, the elastic energy of the spring U ∼ H〈R2 〉eq ∼ HR2 0<br />

to the thermal energy kBT , which implies B = (2R2 0 /τ)1/2 .<br />

U ( x1)<br />

R(t)<br />

U (x2 )<br />

Figure 2.4: Sketch of the dumbbell in a velocity field u(x).<br />

In a non-homogeneous flow, polymers can also get stretched, because of different<br />

velocities of the two beads (see fig. 2.4). Therefore, a term ˙ R = u(x2, t) −<br />

u(x1, t) must be added to the evolution equation (2.14). Since the flow is smooth<br />

at the scale of the polymer, we can approximate this stretching term with the ve-<br />

37

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!