1. Introduction
The rheological behaviors of gelled waxy crude oils during yielding process is a very relevant problem for theoretical rheology, as it exhibits the consequences of drastic structural modifications at a microscopic level, leading the system from a solid-like behavior to a thixotropic fluid behavior. Furthermore, it is also very useful because it may help engineers make smart decisions on the operation strategies for pipelines. For instance, before it enters the pipeline, waxy crude oil is often heated to a relatively high temperature to keep it from gelation. While during pipeline shutdown, heat supplement no longer enters the pipe and drastic temperature drop could happen. As a result, waxy components once dissolved in waxy crude may precipitate from the oil when the temperature falls below its wax appearance temperature (WAT). If the shutdown lasts for a long time, irreversible gelation could happen in the pipe since a more stable structure has developed due to high fraction of wax precipitates. Quantitative description of the viscoelastic behavior of a waxy crude gel is helpful in estimating the minimum pressure and time needed to restart the flow after a period of shutdown [
1,
2,
3,
4]. Such practical relevance and consequent economic implications of yielding behavior add to its theoretical interest and experimental challenges, making it a topic of great significance. In this work, we aim to model this nonlinear viscoelastic behavior of waxy crude oil as the yielding of a fractal gel.
The yielding behavior of waxy crude oil has many complexities [
5,
6,
7,
8,
9,
10,
11]. The flow curves in
Figure 1 (whose experimental details will be fully discussed in
Section 3) are typical ones of yielding. To quantitatively describe them, two forms of thixotropy models have been developed, both of which consist of an evolution equation for stress and for the structure, characterized by an abstract structure parameter
, varying from 1 (pure elastic behavior) to 0 (pure viscous behavior). Depending on the rheology of the stress equation, the two categories are often named as the viscoplastic thixotropic models and the viscoelastic thixotropic models.
The viscoplastic thixotropic models widely adopted in the oil industry are developed on the basis of a pseudo-plastic model, such as the Bingham model or the Herschel-Buckley model [
13,
14,
15,
16]. These models consider the waxy gel as a rigid body before yielding and as a viscous fluid with decreasing viscosity with time immediately after a critical stress is loaded. The viscoelastic models, or the elasto-viscoplastic models, incorporate the description of viscoelasticity during pre-yielding period by using stress equation [
17,
18,
19,
20]. In most of these models, the yield stress term of a steady viscoplastic model is simply replaced with elastic stress. Since the structure-dependent yield stress well reflects the qualitative features of thixotropy during yielding, the predictive ability of these two kinds of models is usually found to be quite satisfactory [
21].
The shortcoming of the viscoplastic thixotropic models is obvious. They fail to give a rheological description in the pre-yielding region. Moreover, the presence of the yield stress even after yielding in those models is theoretically problematic. Further problems of the existing viscoelastic-thixotropic models for crude oils are: (1) Their stress equations are not developed on mechanical basis, and the steady-state shear stress-shear rate relationship could not be recovered by setting the time derivatives in the evolution equation to zero; (2) Too many parameters with no specific physical meaning have to be determined through complex numerical fitting.
Alternatively, the models based on a linear viscoelastic stress equation developed from classic mechanical models, such as Maxwell model or Jeffrey model, are able to describe the pre-yielding rheological behavior, and thus enhance their theoretical interest [
22,
23,
24,
25,
26,
27], among which the model proposed by de Souza Mendes and R. L. Thompson [
27] has been modified and applied to describe the rheological behavior of waxy crude oil by Van Der Geest et al. [
28]. However, quantitative description of the gel microstructure evolution in the early stage of yielding does not contribute to any of these models, since most them are not developed on the basis of actual material, so that experimental aspects could not be taken advantage of in understanding the yielding mechanics. As a result, all of them use an ambiguous structure parameter
to describe the structure, which seems very difficult to validate. Moreover, using a linear viscoelastic stress model (such as Maxwell or Kelvin-Voigt model), the effects of the thixotropy caused by drastic structural changes during yielding process on the mechanical properties, i.e., the dependence of the elastic modulus
G and the viscosity
on time
t, is not accounted for. However, since thixotropy is caused by structural evolution of the gel, it should have been taken into consideration as a main character when modeling the yielding of structured gels such as waxy crude oil itself. Mendes et al. [
29] have chosen a different way of modelling the rheological behavior of waxy crude oil by involving both flow and temperature history dependence, based on new findings about the structure destruction characteristics from the experimental results of waxy crude oil. Since waxy crude oils are complicated mixtures with special characteristics, their rheological modelling needs further studies.
In this paper, for the evolution equation of shear stress, we employ a nonlinear viscoelastic model proposed by Marrucci and used to model the viscoelasticity of structured polymers [
30,
31,
32,
33,
34]. We propose an explicit identification for the structure parameter
, which is often used in an abstract way in other thixotropic models [
16,
17,
18,
19,
22,
23,
24,
25,
26,
27], and incorporate the structure change by relating
to the several material coefficients in the stress equation. In order to do so, we consider the gel as a fractal structure with a characteristic fractal dimension
D. Based on a model for the yielding of weakly flocculated suspensions [
35,
36], we are able to estimate the viscosity and the discussion of the yielding point from a microscopic point of view, in terms of a characteristic radius
. Another particular aspect we consider here is the time variation of shear rate
(and therefore of shear strain
). Instead of assuming
as an imposed constant along the whole process, we find the inertial transient behavior from 0 s
to the final imposed steady-state value of
with an overshoot in between. We use expressions closer to the experimental observations and analyze their influence on the results. This noteworthy influence occurs in the yielding experiments when the imposed shear rate is driven by a controlled-stress rheometer, which is not rare in References [
7,
8,
20,
29,
37], considering the fact that controlled-stress rheometers are more common than controlled-strain ones.
The article is organized as follows—in
Section 2, we discuss the evolution of the structure during yielding from the microscopic perspective, based on which the theoretical justification of the model and its viscoelastic description for both elastic modulus
G and viscosity
is proposed. In
Section 3, experimental protocols and material properties used to obtain the rheological parameters of the model so as to validate the model are presented. In
Section 4 and
Section 5, we summarize detailed computing method of the model, and then present the results under different hypotheses and comment on the merit and weakness of each one. In
Section 6, the concluding remarks are given and the direction of future investigation is discussed.
2. Theory
Our aim is to express the shear stress
as a function of strain
or of time
t in the yielding region. First, we use a fractal model proposed by Snabre and Mills [
35,
36] to interpret the structural evolution of the waxy crude gel corresponding to its mechanical behavior in yielding, and to analyze the critical strain of yielding, which will be used in the following modelization. For the macroscopic description of the evolution of
containing both elastic and viscous contributions, we use a thixotropy modelization similar to that in earlier works [
23,
25,
32,
38].
2.1. Microscopic Description
In this section, we use a blob model to characterize the flocculated clusters in the waxy crude gel. Under this assumption the critical transition conditions between two successive stages and microstructural changes during each stage will be analysed, among which the most important transition condition—the critical strain of fracture yielding—will be estimated theoretically.
2.1.1. Fractal Model of the Gel
Pioneer studies in References [
8,
39,
40] show that rheological properties, for instance, the viscosity and elastic modulus of waxy crude oil, may be well described by fractal models. Also, waxy crude gels and fractal systems have been found to have many morphological features as well as gelation mechanisms in common [
9]. Therefore, from experimental point of view, waxy crude oil shares many basic features with fractal models. Here, we consider waxy crude gel as a fractal cluster network constituted by gel clusters of characteristic size
accumulating primary spherical particles of smaller size
a (precipitated wax crystals, in our case), and characterized by fractal dimension
D, meaning that the number
N of primary particles in each cluster depends on
a, and
as
Furthermore, following Snabre and Mill’s work [
35,
36], we assume that the clusters in waxy crude oil may be conceptually considered as spherical blobs of radius
, filling the volume of the system with a volume fraction of
and a maximum packing volume fraction of
(
for random packing of rigid spheres). The sketch of the configuration of the blobs in different stage of yielding in
Figure 2 is to help understand the yielding of the gel from the microscopic prospective. As long as it is in the gel state, the blob radius is assumed to stay
, and the volume fraction occupied by the blobs
will always be
. After the blobs are broken and the system enters the sol state,
will decrease with
following the time evolution equation which may be obtained from the kinetic equation of the structure parameter
(see
Section 2.2.2). Because of the fractal assumption Equation (
1), both
and
may be related to the volume fraction of the primary particles
as [
8]
Thus, the relationship between
and
, both before and after the breaking of the blobs, could be described by
2.1.2. Microstructure Evolution and Characteristic Transition Points
As from the macroscopic point of view, the mechanisms during different stages of yielding—elastic deformation region, creeping stage, and fracture stage—are significantly different, we put our focus on interpreting what happens to the microstructure in each stage, and discuss the rheological behavior separately as well. From the flow curves in
Figure 1a, four distinct regions can be identified:
The initial overlapping part of the two curves (corresponding to different nominal values of the shear rate )—from the beginning up to point A—can be interpreted as an elastic deformation of the gel, with an infinite or extremely high viscosity;
The region between point A and peak-stress point B can be interpreted as to a viscoelastic creep, where elastic modulus decreases and viscosity becomes finite, but elastic behavior still dominates;
In the region beyond the peak point B up to the inflexion point C, viscous flow becomes increasingly dominant and the microstructure of the initial gel breaks in a steep way, leading to a reduction in both elastic modulus and viscosity;
The region between point C and the steady-state point D corresponds to thixotropic behavior of the the remaining structure, mostly a viscous behavior with the viscosity decrease slowing down until a steady-state is achieved.
Upon the blob assumption of the microstructure of the gel proposed in
Section 2.1.1, we may interpret the critical points on the flow curves in
Figure 1a as follows: (1) point A as the elastic-limit point where bonds between gel particles inside the blobs begin to break (when the spherical blobs are increasingly turned into ellipsoidal blobs); (2) point B as the structure fractural yielding point (corresponding to the strain value
) where the blobs become able to move relatively to each other because of the removing of the geometric constrictions to their motion so that the gel network collapses as a whole and the blobs begin to break; (3) point C (corresponding to the strain value
) as the inflexion point, where breaking process begins to slow down.
From the beginning up to point A, the gel deforms elastically and one may imagine this as a relative deformation of the blobs with the density of bonds between them staying constant (
Figure 2a). From point A to Point B, further deformation breaks internal bonds and stretches out the blobs from spheres into ellipsoids in the direction of shear. The breaking of the bonds will lead to a reduction of the elastic effects, while the orientation of the blobs will make more space for them to move, thus generating viscous effect. But as long as
, the bond density still stays high enough to hold all those primary particles inside the blobs where they originally belong to (
Figure 2b). And since the volume fraction of the blobs does not change, it will remain in the gel phase. For
, that is, from point B to C, the bonds break with an increasingly reduced proportion and there are small particles leaving their original blob. Thus, the blobs become small enough to freely move with respect to each other, and the system enters the sol phase (
Figure 2c). For
that is, from point C to D, the bonds between particles in individual blobs become weaker due to larger strain and more intense shear, and thus the radius of the blobs will decrease with time (
Figure 2d).
2.1.3. Estimation of the Yield Strain
While the thixotropy models of waxy crude oil have always given a satisfactory description of the fourth stage (from C to D), our focus will, alternatively, be the first three stages, that is, from the beginning to C, which are closely related to the transition from gel to sol state.
To estimate the value of
, let us consider four of the spherical blobs in contact in the compact gel, each one of radius
, two in the upper line and the other two in the lower line (see
Figure 3a). The separation
l between the upper and lower horizontal planes in tangent contact with the spheres is
. In the configuration of the
Figure 3a, the upper sphere cannot move with respect to the lower ones, because of geometric restrictions. In
Figure 3b, two deformed ellipsoidal blobs may stand vertically one over the other, and thus relative motion is enabled. In such critical configuration, the critical radius
of the inscribed circle of any blob ellipse, whose semi-minor axis is equal to
, may be given by
.
A direct estimation of critical strain is given by the tangent of the angle
formed by the lines AC and BC in
Figure 3a. In order that the upper sphere may slide over the lower one (with some vertical contraction of both of them), the contact point A of both spheres must arrive at position B. The corresponding tangent is given by
, which provides an alternative estimation for the critical strain. In fact, the actual yield point is expected be somewhat smaller than 0.27 if one takes into account that the spheres will be compressed in the vertical direction as they are sheared one with respect to the other.
Alternatively, we may estimate the yield strain in a more accurate way. The deformation is a combined effect of shear and contraction. We assume a horizontal displacement
along the upper surface due to simple shear, plus a vertical contraction. In order to start relative motion, the blob must be tilted with an angle
, whose cosine may be expressed as
, with
half the distance between the two points (B’ and C’ in
Figure 3b) on the elliptical blob in contact with neighboring blobs, and
equal to its critical vertical projection. The vertical projection
(as obtained in the former paragraph), and
may be taken
as an approximation, in such a way that
, and
. Actually, the distance between B’ and C’ is slightly lower than
, so the actual yield strain must be lower than 0.40. We may estimate this strain more precisely with the aid of analytic geometry, and the results is shown to be 0.25 (see
Appendix A).
In our analysis we will take
(see Table 3). The experimental results of yield strain 0.2 at 31
C to 1.0 at 35
C [
41] are indeed in this range.
2.2. The Model
In this section, the theoretical concept of the development of the model will be given in detail. Based on a viscoelastic model for the description of mechanics, a kinetic equation for the structure and two material functions for the rheological parameters of the gel at a certain structured level are also proposed.In the viscoelastic model, two material parameters appear—the elastic modulus G and viscosity , which depend on the structure of the material, and on and . To close the description, an evolution equation of the structure and its relation to G and are needed. Thus, two evolution equations, for the stress and for the structure parameter are needed, complemented with functions for , appearing in the stress equation, and for the building-up and the destruction terms of structure and , appearing in the evolution equation for , respectively. For the sake of higher efficiency, we will slightly modify the general scheme of modeling the kinetics for the structure by expressing it in terms of .
2.2.1. Mechanical Description: The Nonlinear Viscoelastic Model
In previous models, the evolution of shear stress
is often described by linear viscoelastic models like Maxwell-like or Jeffrey-like models, which do not account for the time-dependent behavior of
G. In our model, instead, the shear stress
will be described by the nonlinear viscoelastic model
Equation (
4) could be seen as the nonlinear version of Maxwell model. It has the same mechanical analogical model as Maxwell model, only with the time dependence, or the strain dependence of the material parameters considered. And it reduces to the Maxwell model upon constant
and
G. This kind of model was first used by Marrucci et al. [
30,
31,
32,
33,
34] for the description of structured suspensions viscoelasticity, but so far has not been applied to waxy crude oil, whose elastic modulus
G dependence on time
t or on shear strain
should not be ignored as well due to the dramatic structural changes during its yielding. The shear modulus
and the viscosity
will be given in
Section 2.2.3 as Equations (
14) and (
15) respectively. We will show that Equation (
4) leads to better results than using the linear Maxwell equation in the discussion
Section 5.2.
2.2.2. Kinetic Equation and Physical Interpretation of the Structure Parameter
The second aspect to be considered is the structure evolution of the system, which in its turn influences
G and
. Since the microstructure is very complicated, it is usually simplified as a scalar parameter
which goes from
for the completely structured gel to
for the completely unstructured fluid [
16,
17,
18,
19,
22,
23,
24,
25,
26]. Of course, this is a very radical simplification, but it may cause ambiguity in modelling and validation. So in this section, we will first make our definition for the structure parameter
on the basis of the features of the microstructure of the gel.
a. Physical interpretation of
In searching for an explicit physical interpretation for the structure parameter
, let us recall the microstructure of the gel, which has been discussed in
Section 2.1.1. Since
is the main feature that describes the gel structure, we define our structure parameter
as
In the system of a fractal gel, while the microstructure is often characterized by size
, the material parameters such as the elastic modulus and viscosity are mostly affected by the the concentration of gel clusters. Therefore, an expression for
in terms of the volume fraction of the mentioned blobs
, rather than in terms of their radii
is needed. According to the relation between
and
described by Equation (
3), the structure parameter
may also be interpreted as
Note that for the fully gel phase,
and
; and to the completely unstructured fluid
, and
, as usually required to recover the steady-state viscosity relation. Note, finally, that a yet more detailed analysis could assume a lamellar structure of the gel instead of an isotropic fractal structure, as observed in microscopic analyses [
42], but we will leave these details for future work. Furthermore, the fractal dimension
D itself could change along the yielding process, but here we kept it constant.
b. Kinetic equation for
The evolution of the structure parameter
consists of a build-up term
describing structure formation and a destruction term
describing structural degradation. Both terms depend on the shear conditions, and thus may be expressed in terms of stress [
25,
26,
27], of the shear strain [
37], or of the dissipation rate [
18,
23]. In the literature, this equation is often expressed as
with
,
the times characterizing the growth and decay of
, and
,
characteristic exponents.
describes the structural evolution dependence on the instantaneous shear condition, which could be assumed, for instance, to be a univariate dependence on
as References [
17,
19,
37] did in their work.
Departing from the assumption that
is decreasing from the very beginning of loading (
) made in other models [
25,
26,
27], we assume that
as long as the gel structure maintains intact, and only begins to decrease after the critical strain
for yielding (which has been discussed in
Section 2.1.2). Thus,
in Equation (
7) could be
in which
is the Heaviside step function with a zero value for
and unit value for
,
is the shear rate above which most of the blobs have become primary particles, and
n is an exponent. If we take
for simplicity, we obtain the steady state value of
In the recent work of Mendes et al. [
29], the microstructure evolution was found to depend more directly on shear strain, rather than on shear rate. Moreover, as is discussed in
Section 2.1.2 , the yielding point may also be interpreted in terms of shear strain, therefore, instead of expressing the changing rate of structure in terms of
t like Equation (
7), here we choose to express the evolution of
directly as a function of shear strain
. They also proved that the extent of structure destruction depends on the most severe shear condition that the oil gel has experienced during its yielding process, which in our case would be the steady-state shear rate. Thus, we consider a same steady-state structure
as given by Equation (
9), and propose the evolution of structure parameter following relaxation equation of
as
where
is a value of the shear strain characterizing the change of
from the initial value to the steady-state value of
. Of course, the time-dependence rate of
may be obtained as in Equation (
10), by applying the chain rule
. This is especially easy if the dependence of shear rate on time
is explicitly given, for instance, if a constant, stepwise or oscillatory shear rate is applied on the system. However, in the particular situation under examination, where the shear rate depends on time in a more complicated way (see
Section 4.1), it would take some serious effort to obtain the solution of Equation (
7) as a function
t, so we focus our attention on Equation (
10) instead. Below, once
is given, we may obtain
. This is a further original aspect of the present work, which allows to consider the
curve in a more direct way. Anyway, after considering
, in
Section 6, we will also be able to express the time evolution of the shear stress
.
With the physical interpretation of
in Equation (
5), a shear-strain dependence of the evolution of
, describing how the blob breaking in the region from point B to C in
Figure 1a, may be obtained. Since the evolution equation of
as a function of
has been given by Equation (
10), thus by solving it, one obtains the transition behavior of
between
and
as
in which
is the initial value of the blob radius corresponding to the radius of the maximum packing situation (gel state), and
the steady-state radius. According to Equation (
9),
is also to depend on the imposed
as
in which
is phenomenological parameter related to the relaxation time of both blob break and particle reunion,
is the critical strain above which most blobs have been breaking into primary particles, and
q is an empirical exponent. This expression reflects the obvious fact that the characteristic blob radius will decrease with increasing
, because primary particles will be stripped off more easily from the clusters. Thus, in the initial breaking region B ∼ C,
decreases with
, from its initial value
to its steady-state value
. Practically, all the parameters in Equation (
12) could be obtained by fitting the steady-state viscosity at several
, which are related to
through relation Equation (
15). Technically, these parameter could also be obtained through microscopic experiments. As for our analysis, it turns out
,
and
fits the results well.
2.2.3. Material Functions and Rheological Properties
The other two functions needed to describe the process are
G and
appearing in Equation (
4). The second one is dependent on
, because viscosity of suspensions is known to strongly depend on the radius of particles. In contrast, no such strong dependence should be expected for
G, because in the gel phase
, and the only effect on
G is that internal bonds between the blobs are increasingly broken with higher
in the gel phase when the blobs deform.
a. Elastic modulus
We assume the elastic modulus
G is proportional to the number density of full chains of bonds times the elastic constant of each bond, and a certain number of such bonds break at a given shear strain. According to this assumption, the nonlinear elastic modulus
of the system is given by expression
with
the equilibrium value of elastic modulus,
the characteristic shear strain, and
p the phenomenological exponent. Equation (
13) of
G exhibits a steep drop at
after a plateau equal to
for
, and followed by a very small value for
, with
p controlling the steepness of the reduction of
around
. This form of
is analogous to that used in References [
16,
43], with the difference that they used elastic part of shear strain
instead of the directly measured strain
.
depends on the conditions of gelation, for instance, whether it was in a quiescent state or under shear, and how fast the cooling rate was [
8]. It could be tested with rheometer by applying a small amplitude oscillatory shear to the waxy gel sample after it rests at the test temperature for a long time. Comparing the behavior of
to that of our situation, and since a sudden reduction in
G reasonably corresponds to the collapse of the gel network,
will be taken close to the yield strain
.
As the imposed shear rate is as low as 0.1 or 1 s
, it is only reasonable to consider the initial behavior as purely elastic. Therefore, the elastic modulus
G within the strain range
could be calculated by
. The calculated elastic modulus shown in
Figure 4 is seen to start falling from a very high value at the very beginning of the deformation. This could arise from a pre-collapse of gel network at a rather small strain
, which is hardly noticed because of its short-term existence. Thus, one may imagine, that besides the gel network which we have seen to break around the yield strain
, with an equilibrium elastic modulus
which can be obtained from rheometric test, there is another more brittle and rigid network in the gel, with smaller characteristic strain of breaking
and higher initial elastic modulus
. Such double network structure would have an elastic modulus equal to the summation of moduli of both networks, each of which is given in the form of Equation (
13)
In practical terms, the critical shear strain
would imply the strain for which the elastic modulus of the network
i drops with a steepness characterized by
(related to the fractal dimension of each network) due to the breaking of the corresponding type of bonds. Therefore, for small strains
, one would have a plateau of modulus
G of approximately
, and between
and
, there might be another plateau of approximately
, whose width and position depend on the separation between
and
and on steepness of the two drop behaviors, that is, on the values of exponents
. This behavior may be caused by asphaltene or resin macromolecules which coexist with wax crystals. Its influence will be further explored and justified on experimental grounds in
Section 5.1.
b. Viscosity
To estimate the viscous effect beyond point A, we use the viscosity model for suspensions with viscoelastic aggregates proposed by Snabre and Mills [
35], which relates
to the volume fraction of aggregates
and deformation tensor
as
where
is the solvent viscosity,
the maximum packing volume fraction (64% for random close packing of rigid spheres),
is the
component of the deformation tensor
, and
a parameter characterizing the viscous effect generated during creeping. For
, we take a simple assumption that it is proportional to
until the structural failure of the network, beyond which it becomes a saturated constant. Equation (
15) is developed on the basis of a viscosity model for suspensions with rigidly flocculated blobs [
35], by taking into account the effect of the elastic deformation of the blobs. Instead of being constant as in the model of the rigid-blob suspension, the effective maximum packing volume fraction
in Equation (
15) is assumed to increase with deformation as
. This correction makes the volume fraction of the blobs
smaller than
so that a high but finite viscosity is expected before fracture yielding. It reflects the physics discussed in
Section 2.1.1—since the blobs are elongated and become more orientated, bigger room is made for their movement. Thus, the internal flow through the network becomes less congested, leading to a higher fluidity, but still an inconspicuous viscous effect before the blobs are broken.
6. Conclusions
We have presented a simplified theoretical model for the yielding of a waxy crude oil by assuming the gel as a fractal structure. The predicted results of the model show a very close form in comparison with experimental curves. It gives a quantitative description of the isothermal nonlinear viscoelastic behavior of waxy crude gels formed under a given constant cooling rate, and may be used to estimate the minimum pressure and time needed to restart the flow after a period of shutdown. It has sufficient simplicity and wide flexibility to explore the influence of different physical factors, such as fractal dimension, shear rate history or the initial elastic modulus. Moreover, since we have not used any specific details of a waxy crude oil in the modeling process, the main framework of our model is sufficiently general to be also applicable to the yielding of other fractal gels.
Our model is based on evolution equations for stress and for structural parameter , and on a modelization for and . It starts from a definite physical model for the structure parameter instead of from more abstract parameterizations. Compared to previous works, five novel factors especially considered in our work have consequent advantages, namely: (1) the use of Marrucci’s nonlinear viscoelastic model in descriptions for the stress evolution; (2) The use of an explicit microscopic interpretation for the structure parameter in terms of the radius of the gel clusters, which allows us to estimate a range of the yield strain with simple arguments; (3) the material function for shear modulus G inspired by the double-network model; (4) the application of well-established sophisticated viscosity model to the description of the viscosity of a transient process; and (5) the modification the actual imposed shear rate history which deviates from the setting one due to device inertia.
Though the main framework of the model is very general, we have focused our attention on a concrete application. Further studies should be the following ones:
(1) An extension to more values of nominal shear rate should be examined. Furthermore, applications to other typical rheometric situations such as oscillating shear and small-step increase of strain should also be carried out.
(2) Physical grounds for the inertia effect of device as well as the double-network structure of elastic modulus still need further exploration. The empirical function Equation (21) for the actual shear rate history and the material function Equation (18) for the elastic modulus should be improved accordingly.
(3) Based on the presented model developed for the yielding region and applied to a fixed temperature history and composition situation, a comprehensive description of the rheological behavior may be extended to. First, the composition and temperature dependence of the material functions should be studied to compare the predictions of the model in different circumstances. Secondly, the structural evolution needs to be understood microscopically, leading to a more detailed evolution equation on the basis of a rate function incorporating the build-up and destruction of structure.