Journal ArticleDOI

# Derivation of the Marrucci model from transient-network theory

01 Jan 1981-Journal of Non-newtonian Fluid Mechanics (Elsevier)-Vol. 8, pp 183-190

TL;DR: Many theories about the rheological behaviour of melts and concentrated solutions of high molecular weight polymers are based upon the transientnetwork model, originally developed by Green and Tobolsky, Lodge, and Yamamoto.

AbstractMany theories about the rheological behaviour of melts and concentrated solutions of high molecular weight polymers are based upon the transientnetwork model, originally developed by Green and Tobolsky [l], Lodge [2], and Yamamoto [ 31. Among these theories a model proposed by Marrucci and coworkers [4,5] turns out to be rather successful [6-111. The basic equations of this model, as formulated in ref. [ 51, are:

Topics:

### G Oi = cnoikT. (35)

• At first sight, this equation differs from (5), since the creation terms in both equations are not of the same form.
• It has been shown that, although introduced originally in a more or less empirical way, eqns.
• Finally it should be noted that the present derivation offers the possibility of detailed comparison of the Marrucci model with other transient-network theories based on different forms of the functions gi and hi in eqn. (23).

Did you find this useful? Give us your feedback

Content maybe subject to copyright    Report

183
Journal of Non-Newtonian Fluid Mechanics, 8 (1981) 183-190
Elsevier Scientific Publishing Company, Amsterdam - Printed in The Netherlands
Short Communication
DERIVATION OF THE MARRUCCI MODEL FROM
TRANSIENT-NETWORK THEORY
R. J. J.
JONGSCHAAP
Department of Applied Physics, Twente University of Technology, Enschede (The
Netherlands)
1. Introduction
Many theories about the rheological behaviour of melts and concentrated
solutions of high molecular weight polymers are based upon the transient-
network model, originally developed by Green and Tobolsky [l], Lodge [2],
and Yamamoto [ 31. Among these theories a model proposed by Marrucci
and coworkers [4,5] turns out to be rather successful [6-111.
The basic equations of this model, as formulated in ref. [ 51, are:
Gi = GoPi 9
(3)
Ai = xerx;.4 )
(4)
(5)
where (1) represents a spectral decomposition of the stress tensor, Xi and Gi
are the i-th relaxation time and elastic modulus, respectively (the same quan-
tities with zero subscript are the corresponding equilibrium values), 6/St is
the contravariant convective derivative (see eqn. (22) below), D is the rate-of-
strain tensor and where the dimensionless quantities xi (3~~ Q 1) can be
regarded as structural variables which describe how far the existing structure
is away from equilibrium. They are related to the degree of connectivity of
0377-0257/81/0000-0000/\$2.50 @ 1981 Elsevier Scientific Publishing Company

184
the macromolecular network. The constant Q is an adjustable parameter.
The model has proven to predict correctly the mechanical response of vari-
ous polymer melts in different stress and deformation histories. From a physi-
cal point of view it is appealing because of the explicit introduction of an
expression (eqn. (5)) that describes the change of structure of the system.
Despite these attractive features, the derivation of the model, as given in refs.
[4] and [5], is not entirely satisfactory. In fact, eqns. (l-6) are introduced
in a semi-empirical way: the form of the constitutive equation (2) is
suggested by a springdashpot model with a variable spring modulus. The con-
travariant convective derivative is chosen in accordance with results of the
network theory for concentrated systems [2,12] and of bead-spring models
for dilute polymer solutions [13]. The dependence (eqn. (3)) of Gi on xi is
suggested by the theory of rubber elasticity. The dependence (eqn. (4)) of Xi
on Xi is chosen in such a way that the zero shear viscosity becomes propor-
tional to c3R, where c is the number of macromolecules per unit volume.
Finally, the kinetic expression (5) is proposed on the basis of certain micro-
scopic arguments. An objection that may be raised against this procedure is
that the constitutive equation (2) and the kinetic equation (5) are not derived
from the transient-network model (or any other consistent microscopic
In this note it will be shown that it is possible to derive both the constitu-
tive and the kinetic equations of the Marrucci model from a balance law of
the segment-distribution function in the transient-network model. This makes
made in the Marrucci theory and allows a comparison of this theory with
other theories based upon the transient-network model.
2. Distribution functions, balance law
Consider an incompressible fluid with c macromolecules per unit volume.
Let Ni be the number per unit volume of segments, consisting of i freely
jointed rigid links of length I (segments of this type will be called i-segments),
and \ki(q, t) d3q the number of these segments that have their end-to-end vec-
tor q in the volume element d3q in the configuration space. It follows that
s
\ki(q, t) d3q = Ni *
(f-5)
The number of i-segments per molecule will be denoted by ni, so
ni E Ni/C .
For highly entangled systems ni >> 1; then the quantity
(7)
n=
C ni
i
equals the average number of entanglements per molecule. The equilibrium

185
value of ni will be denoted by noi.
Finally we introduce the structural vari-
ables [ 4,5]
Xi E Yti/Tloi ,
(9)
which describe how far the existing structure is away from equilibrium.
Besides the total distribution functions \ki, we also introduce the distribu-
tion functions
9i(q9 t) G tllNi) *it49 t, -
(10)
Assuming Gaussian statistics, in equilibrium we have:
It/p = (bi/x)3’2 eXp(--bifJ2) ,
(11)
where bi = 3/2i12 .
The time dependence of the segment-distribution functions is governed by a
balance equation of the following form:
a\ki
- = -a: m (\kii) + ki - hi\ki e
at
(12)
The first two terms are of the usual form of an equation of continuity in
q-space; the terms ki and hi\ki represent the formation and annihilation of
i-segments respectively. Equation (12) was first derived by Yamamoto [3 3
and can be shown to be consistent with the theory of Lodge [2,12]. An exten-
sive discussion of (12) can be found in two papers by Wiegel and de Bats [14],
u51.
The function 4 = \$4, t) in eqn. (12) represents the motion of the segment
vectors in the configuration space. If affine deformation of the network is
assumed,
cj=L.q,
(13)
where L is the macroscopic velocity gradient of the fluid flow. In this note we
will assume that eqn. (13) holds; the subsequent discussion, however, also
applies for a special type of non&fine deformation, namely if 4 = L - q,
where z = AL) is some effective velocity gradient of the network. In that case
one simply has to replace L by 1 in all our results. A theory of this type has
been formulated recently by Phan Thien and Tanner [16,17]. For our pur-
pose it is of interest to note that the assumption of non&fine deformation in
this case leads to a constitutive equation which is not of the contravariant
Maxwellian type, the form proposed in the Marrucci model.
The creation function ki and the destruction function hi will depend [3,15]
in general on the variables q and i. Usually certain specific forms of these
functions are assumed. The theory of Lodge [12], for instance, essentiaIly
corresponds to the assumption that
ki =giILiO(q) ,
(14)

186
where gi are constants and J/p is given by eqn. (ll), and that the hi are con-
stants. The form (14) of the creation functions is a consequence of the
assumption that the segments are created at a constant rate, and that they
have, at the instant of creation, the same distribution as free chains. The con-
stancy of the hi means that all i-segments have the same constant probability
per unit time of leaving the network
*. In more elaborate theories (see for
instance [16]) the form (14) is often used in combination with the assump-
tion that gi and hi are functions of the mean square extension of the segments,
i.e. gi =gi(<q’>) and hi = hi(<q2>).
We shah make no specific assumptions about the functions k,(q, t) and
h,(q, t) at this stage of the development. For convenience, however, ki will be
written in the form (l4), with gi =gi(q, t).
Substitution of (14) in (12) and integration over all configurations gives
the following differential equations for Ni :
dNi/dt = <gi>O - Ni <hi> 9
(15)
with
and
<hi> = \$-
i
Jhi9id3q = JhiJltd3q m
Substitution of (7) and (9) gives the following expression for the rate of
change of the structural variables Xi :
&i
dt = \$ <gi>O -xi <hi> s
(17)
(18)
If \ki is written as ‘Pi = N~I,!J~ and substituted in (12), one obtains, using (15):
atii
-=
at
-a:*(\$iQ)+i(g’\$P -<gi>O Jli)-(hi-<h>) \$1.
(19)
By multiplying this equation with the diadic qq and integrating over the con-
figuration space, the following expression for the averages
(20)
* In fact, in the theory of Lodge [ 121 an additional parameter, the complexity, on
which the quantities gi and hi may depend, is used. This, however, is not essential for
the present discussion.

187
can be obtained:
it <44>i = NI (<giqq>O - <gi>O <qq>i) - <hiqq> + <hi> <qq>i 3
with the contravariant convected derivative 6/6t defined as
(21)
(22)
for an arbitrary second-order tensor A.
3. Specific assumptions
From now on we assume that gi and hi do not depend explicitly on q but
only on average properties of the state of the fluid. It will be shown that for
deriving the Marrucci model this assumption is sufficient. The specific forms
of the functions gi and hi in that case will be discussed later; at present it is
sufficient to know that gi and hi do not depend on q. The basic equations
(15) and (21) then reduce to
dNi/dt = gi - Nihi
and
(23)
~t<qq>i = ~
i
(<qq>P -<qq>i) P
(24)
respectively .
4. Constitutive equation
The contribution of the i-segments to the stress tensor is given by the usual
expression
Z =Ni<fiq) 9
(25)
in which the entropic forcefi is related to the parameter bi of the equilibrium
distribution (11) in the following way:
fi = 2kTbiq s
(26)
Here k is the Boltzmann constant and T the absolute temperature. Since the
fluid is incompressible, the partial stresses (25) are determined up to an arbi-
trary isotropic pressure. So we may also use partial stresses Ti defined as
Ti = Ti -PiI 9
(27)
where pi are chosen in such a way that Ti = 0 at equilibrium. From (25) and
(26) it follows that
Ti = 2NikTbi(<qq>i -+<q2>fl) 9
(28)

##### Citations
More filters

Journal ArticleDOI
Abstract: A new constitutive equation for whole human blood is derived using ideas drawn from temporary polymer network theory to model the aggregation and disaggregation of erythrocytes in normal human blood at different shear rates. Each erythrocyte is represented by a dumbbell. The use of a linear spring law in the dumbbells leads to a multi-mode generalized Maxwell equation for the elastic stress and both the relaxation times and viscosities are functions of a time-dependent structure variable. An approximate constitutive equation is derived by choosing a single mode corresponding to the cell aggregate size where the largest number of cells are to be found. This size is identified in the case of steady flows. The model exhibits shear-thinning, viscoelasticity and thixotropy and these are clearly related to the microstructural properties of the fluid. Agreement with the experimental data of Bureau et al. [M. Bureau, J.C. Healy, D. Bourgoin, M. Joly, Rheological hysteresis of blood at low shear rate, Biorheology 17 (1980) 191–203] in the case of a simple triangular step shear rate flow is convincing.

133 citations

Journal ArticleDOI
Abstract: Structural rheological modelling of complex fluids developed in Part I of this series and applied to shear thickening systems (Parts II & III), is now used to improve such a modelling in the case of unsteady behaviour, that is, in the presence of thixotropy. The model is based on an explicit viscosity-structure relationship, η ( S ), between the viscosity and a structural variable S. Under unsteady conditions, characterized by a reduced shear, Γ ( t ), shear-induced structural change obeys a kinetic equation (through shear-dependent relaxation times). The general solution of this equation is a time-dependent function, S ( t ) ≡ S [ t, Γ ( t )]. Thixotropy is automatically modelled by introducing S [ t, Γ ( t )] into η ( S ) which leads directly to η ( t ) ≡ η [ t, Γ ( t )], without the need for any additional assumptions in the model. Moreover, whilst observation of linear elasticity requires small enough deformation i.e. no change in the structure, larger deformations cause structural buildup/breakdown, i.e. the presence of thixotropy, and hence leads to a special case of non-linear viscoelasticity that can be called “thixoelasticity”. Predictions of a modified Maxwell equation, obtained by using the above-defined η ( S ) and assuming G = G 0 S (where G 0 is the shear modulus in the resting state defined by S = 1) are discussed in the case of start-up and relaxation tests. Similarly modified Maxwell-Jeffreys and Burger equations are used to predict creep tests and hysteresis loops. Discussion of model predictions Maynly concerns (i) effects of varying model variables or/and applied shear rate conditions and (ii) comparison with some experimental data.

78 citations

Journal ArticleDOI
Abstract: In the presented model elements of polymer melt rheology and polymer crystallization kinetics are combined. In particular, the proneness of the melt to the special type of crystallization which is characteristic for shear treatment is supposed to emerge only gradually during shear flow. Following Avrami's early ideas on crystal growth, an induction time is introduced. In principle, the model can be applied to any flow and temperature history. The special case of isothermal flow at constant shear rate is covered in greater detail: A favorable comparison is made with experimental results, as published by Lagasse and Maxwell [10].

50 citations

Journal ArticleDOI
Abstract: A consitutive theory for polymetric liquids based on the transient network concept is developed, following Wiegel and Jongschaap. The Phan-Thien-Tanner equation is shown to follow from the general theory with two critical assumption, one of these is quasi-equilibrium of the internal structure, which preludes consistency of application in ”fast flows“. The Marrucci model can be made consistent with the general format with a small change in the kinetic equation that can be deduced from asymptotic behaviour and leaves the steady viscometric behavior unchanged. The simplest genaral formulation requires the linear viscoelastic spectrum and two parameters; the latter cannot be determined uniquely from steady viscometric flow data.

42 citations

Journal ArticleDOI
Abstract: A new numerical approach is presented to exactly solve the convection equation arising in network theories. The method is based on a direct stochastic interpretation of the convection equation. We show that with this approach models can be studied extensively which are not solvable analytically. It turns out that a conceptually simple approach to network theories predicts a qualitatively satisfying rheological behavior.

38 citations

##### References
More filters

Journal ArticleDOI
Abstract: A constitutive equation is derived from a Lodge—Yamamoto type of network theory for polymeric fluids. The network junctions are not assumed to move strictly as points of the continuum but allowed a certain “effective slip”. The rates of creation and destruction of junctions are assumed to depend on the instantaneous elastic energy of the network, or equivalently, the average extension of the network strand, in a simple manner. Agreement between model predictions and the I.U.P.A.C. data on L.D.P.E. is good.

979 citations

Journal ArticleDOI
Abstract: A nonlinear constitutive equation for polymer melts and concentrated solutions is derived from a Lodge‐Yamamoto type of network theory. The network junctions are postulated to move nonaffinely in a well‐defined manner. The functional form of the creation and destruction rates of junctions is assumed to depend on the average extension of the network strand and the absolute temperature in such a way to allow for the time‐temperature superposition principle. The theory shows good agreement with all data examined. The paper concludes with a strong flow problem (melt spinning). The results indicate the validity of the model in this flow regime.

449 citations

Journal ArticleDOI

378 citations

Journal ArticleDOI
Abstract: Many kinds of amorphous high polymer substances show the non-ideal rheological behaviours even in the “static” observations. In this paper, the general stress-strain-time relation and energy of dissipation of these systems are treated by adopting the modified rubber-like network model in which the junctions connecting the polymer chains break up and reform continually. The probability of the breakage of the junction, namely the breakage of the chain, per unit time is not a constant but a functions of both the number of segments which construct the chain and the end-to-end distance of the chain. It is assumed that the velocities of deformation of the system and of the change of the situation of network structure are extremely slower than those of the micro-Brownian movements of polymer segments. Therefore the statistical mechanical considerations of equilibrium states may be used to derive the stress-strain-time ( S – S – T )relation at each “instance” of macroscopic observations. From our general expressi...

260 citations

Journal ArticleDOI
Abstract: In this mainly expository paper, constitutive equations based on the network models ofYamamoto,Lodge, andKaye are re-derived in a common notation involving the use of base vectors embedded in the deforming macroscopic continuum. The derivations are thereby simplified in some respects and the differences of detail between the models are clarified. InLodges theory, the sub-network superposition assumption is replaced by alternative assumptions concerning the creation and loss of network segments, and the theory is extended to non-Gaussian networks.Kayes theory is extended to allow for the presence of entanglement junctions of different complexities.

209 citations