Open access peer-reviewed chapter

On the Thermodynamic Consistency of a Two Micro-Structured Thixotropic Constitutive Model

By Hilbeth P. Azikri de Deus and Mikhail Itskov

Submitted: November 22nd 2017Reviewed: February 26th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.75987

Downloaded: 517


The time-dependent rheological behavior of the thixotropic fluids is presented in various industrial fields (cosmetics, food, oil, etc.). Usually, a couple of equations define constitutive model for thixotropic substances: a constitutive equation based on linear viscoelastic models and a rate equation (an equation related to the micro-structural evolution of the substance). Many constitutive models do not take into account the micro-structural dependence of the shear modulus and viscosity in the dynamic principles from which are developed. The modified Jeffreys model (considering only one single micro-structure type) does not show this incoherence in its formulation. In this chapter, a constitutive model for thixotropic fluids, based on modified Jeffreys model, is presented with the addition of one more micro-structure type, besides of comments on some possible generalizations. The rheological coherence of this constitutive model and thermodynamic consistency are analyzed too. This model takes into account a simple isothermal laminar shear flows, and the micro-structures dynamics are relate to Brownian motion and de Gennes Reptation model via the Smoluchowski™s coagulation theory.


  • thixotropic fluids
  • thermodynamic consistency
  • modified Jeffreys model

1. Introduction

The thixotropic substances are considered structured fluids and have a rheological behavior (time-dependent) guided by their structural nature [15]. In terms of models (constitutive systems), this behavior is represented by a couple of time-dependent equations. These equations connect the micro-structural characteristics to the rheological behavior. In many works, this system of equations is based on a qualitative way, without any formal justification on the physical principles. In this chapter, a different approach is presented based on some well-established physical principles (Smoluchowski’s coagulation theory [610], Brownian motion [69] and de Gennes reptation model [8, 11, 12]).

The rheological properties evolution, in many thixotropy models [2, 5, 1317], is presented, generically, in terms of a couple of equations, as follows:

τ=τλγ̇,constitutive relationship basedonlinear viscoelastic models;E1
λ̇=Gγ̇λγ̇,relation associated to the microstructural evolution,E2

where τis the shear stress (from Cauchy stress measure), γ̇is the shear rate (infinitesimal strain measure), λis the structural parameter (a positive scalar quantity associated to the structural level of substance) and Gγ̇is a functional form. In large number of works [5, 1821], Eq. (1) is represented by a linear viscoelastic model/mechanism (Maxwell, Jeffreys, Kelvin-Voigt, etc.) or combinations of that (in parallel or in series). However, it exists a problem with these approaches, they do not consider, the time/micro-structural dependence of the shear modulus (G) and viscosity (ημand ην) in physical principles from which the constitutive equation system is obtained. Differently, the approach presented in [8] (the modified Jeffreys model) does not present this mistake, considering time/micro-structural evolution in the constitutive system development, and in this way, showing more coherence to represent thixotropic substances behavior.

Eqs. (1) and (2) connect micro-structural characteristic and viscoelasticity nature of the thixotropy. In the specialized literature, a considerable number of micro-structure types [2229] can be found, although only one single micro-structure type is considered in the development of the modified Jeffreys model. In this way, it is natural to imagine that it can exist more than one micro-structure type in a same thixotropic substance. In this chapter, a new perspective in terms of constitutive approach for thixotropic substances, with apparent-yield-stress nature and based on the modified Jeffreys model [8], is presented in terms of two different types of micro-structure (i.e., two structural parameters), and some of their properties are investigated in a rheological and thermodynamic sense. In the presented approach, each micro-structure is associated to a specific mechanism considered in the model (Maxwell-like element in parallel with a pure viscous Newtonian-like element), that is, one micro-structure type is associated to the Maxwellian element and another one to the pure viscous element. In this way, it is important to stand out that each of them (mechanisms) can react, under the shear loads, in distinct (but integrated) forms one from the other.

Aiming to extend the modified Jeffreys model for two distinguished micro-structures, in the next sections, the constitutive model is formally presented in terms of equations, the analysis of its thermodynamics consistence is done, points related to the characterization of the transition region are discussed and some illustrative numerical simulations of rheological tests are presented as well.


2. Main ideas on two types of micro-structures approach

The approach proposed in this work is based on a Jeffreys model [8] for the thixotropic system representation. In this case, a sketch of the model (Maxwell element in parallel with viscous element) is presented in Figure 1.

Figure 1.

Sketch of the two micro-structured thixotropic constitutive model.

It is important to point out the presence of two micro-structure types (λμand λν) that are intimately related to the behavior of the shear modulus G=Gλν, and viscosity coefficients ην=ηνλνand ημ=ημλμ. In this way, under a certain shear load (τ), for the Maxwell element (.ν), one has


where γ̇νvis the viscous strain rate and γ̇νeis the elastic strain rate, respectively, and the total strain rate for the Maxwell element is γ̇ν=γ̇νv+γ̇νe. After some algebraic manipulation, it follows:


From the viscous element (.μ),


and remembering that γ̇=γ̇ν=γ̇μand τ=τν+τμτ̇=τ̇ν+τ̇μ, one has


It is important to stand out that the time variation of the shear modulus and the viscosity coefficient are taken into account, which is in agreement with the expected physical behavior of thixotropic fluids [8, 30, 31].

It is important to keep in mind the fact that the parameters λμ(0λμ1) and λν(0λν1) are used to characterize the structural level of the substance. The structural parameter closer to 1implies highly building-up state of micro-structure, and when it is closer to 0, the micro-structure is close to a fully broken-down state, in respective mechanism (.νor .μ). In this sense, it is presented the two micro-structures version of constitutive model based on the work [8].


In this case, there are the following positive parameters: , kν, kμ, βν, βμ, tν, tμ, Go, m, ςν, ςμ, Kν, Kμ, ηo, α1, α2and α3. Some physical details related to these parameters can be found in [8, 32].

3. On the thermodynamic consistence

In this section, the thermodynamic consistence of constitutive model is analyzed in terms of the entropy-producing processes [33] associated with the concept of natural configuration [8, 3436]. The main objective is to verify the consistence of the constitutive equation to the Clausius-Duhem inequality.

Figure 2 presents a sketch of the configurations spaces used in this approach for the body B, where k̂oBis the reference configuration, k̂tBis the current/actual configuration, and k̂ntBis a family of natural configurations. In isothermal process, one is supposed that k̂ntB=k̂ntBλμλν, where λμand λνare functions of the flow history. In this sense, it is assumed that the family of natural configurations can be parametrized by the structural level of the substance (λμ,λν). In this way, it has F¯¯k̂o(gradient of the motion from reference configuration to current configuration), F¯¯k̂nt(gradient of the motion from natural configuration to current configuration) and F¯¯k̂ok̂nt(gradient of the motion from reference configuration to natural configuration). In this way, denoting for each ξi-th relaxation mechanism.


and the extensive thermodynamic Helmholtz potential (Ψ) can be represented as


where follows from chain rule Ψ¯¯,ξiB¯¯k̂nt=Ψ,Iξii+IξiiΨ,IIξiiI¯¯Ψ,IIξiiξiB¯¯k̂nt+IIIξiiΨ,IIIξiiξiB¯¯k̂ntT, with Iξii=trξiB¯¯k̂nt, IIξii=12trξiB¯¯k̂nt2trξiB¯¯k̂nt2and IIIξi=detξiB¯¯k̂nt. Without loss of generality it is chosen k̂ntsuch that ξiF¯¯k̂nt=ξiV¯¯k̂nt(ξiC¯¯k̂nt=ξiB¯¯k̂nt), that is, for an appropriate rotated natural configuration k̂nt, it follows [8]:


Figure 2.

Configuration spaces.

where D¯¯is the symmetric part of velocity gradient (L¯¯). Returning to Eq. (20), one has


Then, from the Clausius-Duhem inequality follows the specific rate of dissipation Ξ


where τ¯¯is the Cauchy stress tensor.

This approach makes possible the analysis of the constitutive law and the nonnegativity of the entropy production in the context of irreversible thermodynamic processes. In this sense, it is assumed that total rate dissipation Ξcan be divided into three parts, each one associated to a specific process, as follows:


where Ξk̂tis the dissipation rate due to changes in k̂tB), Ξλis the dissipation rate due to changes in micro-structure and Ξk̂ntis the dissipation rate due to changes in k̂ntB.

Following similar steps to the work [8], one has




Thus, one has for each ξi-relaxation mechanism.


where =̇L¯¯L¯¯Tis the “Oldroyd time derivative”. Considering the two relaxation mechanisms, with rν=rμ=1ξντ¯¯=τ¯¯νand ξμτ¯¯=τ¯¯μ, where ηξν=ην, Gξν=G, ηξμ=ημand Gξμ+, one has.


Therefore, as τ¯¯=τ¯¯ν+τ¯¯μone has


Note that the isochoric motion constraint was not required in this analysis. In this way, associating each relaxation mechanism to a micro-structure, it is important to point out that depending on the nature (level of complexity) of the analyzed thixotropic substance, the reasoning presented here can be extrapolated for more than two relaxation mechanisms (i.e., more than two micro-structure types) aiming a consistent/coherent (in thermodynamic and rheological senses) approach of the model.

4. On the transition region

The objective of this topic is to analyze the transition/yielding region criteria under a theoretical and formal point of view. In this sense, it is important to stand out that the constitutive equation


can be rewritten in a more appropriate form as follows:




It can be seen that the form of Eq. (36) is quite similar to the standard Jeffreys model. The quantities Π, Γ̇, η, θ1and θ2can be interpreted as a dimensionless stress, dimensionless strain, dimensionless apparent viscosity, relaxation time and retardation time, respectively. Note that since


Eq. (36) corresponds to the standard Jeffreys model. In this way, it is clear the relation between the nonlinear viscoelastic region (associated to the thixotropic effect) and λ̇. Many works propose that an indicative for the beginning of the transition/yielding region is a region where occurs the modification in the behavior of the thixotropic substance, from the linear viscoelastic region to nonlinear viscoelastic region [3741] associating an specific strain value for that. However, this specific value is not constant [4145].

In this sense, taking into account the mapping in terms of variables in the actual configurations (x¯=xig¯i) and the respective displacements (u¯=uig¯i), it follows:


and in this context, one has the Jacobian (J˜)








where δklijand δrstijkare the generalized Kronecker delta [46, 47]. Thus, the relationship between dV(differential of volume in the reference configuration) and dv(differential of volume in the actual configuration) can be described in the following way:


Taking into account that the yield region can be characterized as a transition region. This region, in fact, can be treated as a singularity region. As the continuity axiom holds a suitable asymptotic behavior [4850], the concept of transition phenomenon has an asymptotic nature, and can be treated by limiting approaches. Physically, the transition region (transition/yielding) can be noted when a substance loses part of their original intrinsic properties, and from this point a new set of properties are raised in this substance, resulting in a new behavior. For example, in an elasto-plastic transition behavior, the transition state is related to a strain ellipsoid degeneracy, in a geometric sense, which can be transformed to an infinite cylinder, point or a pair of planes [51].

Considering the abovementioned lines, it can be imagined the nonlinear viscoelastic region as a direct mapping image from the linear viscoelastic region. Thus, the Jacobian (J˜), of this mapping, should translate the singularity on the transition region vicinity as an asymptotic extremely large deformation of the original microscopic element. In other words, taking into account Eq. (47), it can be seen that for a finite initial volume dV, with J˜0(singularity on the neighbor of the transition region) implicates in dv, and in this way, it can write for the transition region the following condition:


or in an asymptotic sense on the transition region (t.r.)


It is important to stand out that this condition was stated in terms of the strain invariants, which can be rewritten in terms of the stress invariants or in terms of energy, by the constitutive relationships. Note the difficulty involved in expressing the strain tensor invariants in terms of the stress tensor invariants, due to the constitutive model’s form.

Remark 1. If it is taken the time rate ofEq. (48),one has


Supposing the case that the strain rate is obtained as a response to known stress load, returning to Eq. (36), it follows:






Thus, one has




Note that


and in this way, from Eq. (53), it follows:


or in other words


where L¯¯=v¯(“Eulerian gradient”), and where v¯is the velocity field. Consequently,


since λμand λνare taken as constants (pre-transition/yield region). In this way, the criteria Eq. (48) and/or Eq. (50) can be appropriately analyzed for a stress load excitation.

5. Numerical results

In this section, some illustrative numerical results are presented. Two types of rheological tests are analyzed, the constant shear rate test (CR) and the constant shear stress test (CS). In both cases it is necessary to define consistent initial conditions in relation to physical and mathematical principles. The finite difference method was implemented in MATLAB, with the following set of parameters: =1, η0=0.08Pa.s, α1=5, α2=α3=0.5, kνςν=103J/m3, kμςμ=102J/m3, βν=5, βμ=3, m=106, Kν=0.001Kgm1K1, Kμ=0.01Kgm1K1, G0=6.5Pa, tνςν=104s1and tμςμ=103s1. It is important to point out that for both tests it was used Δt=103s. It is also important to comment that time-steps Δt=104s,105sand 106swere analyzed, but modifications in responses were not noted. The tolerance used for the Newton’s method was 106. Numerical results for real thixotropic substances (numerical/experimental comparative responses, regularization methods associated to nonlinear identification parameter problem, etc.) can be found in [52].

5.1. Constant shear rate test

In this section, it is considered the constant shear rate test, that is taken into account load conditions as γ̇t=Htγ̇0, where Htrepresents the standard Heaviside function and γ̇0it is a positive real constant. It is reasonable to think that in the begin of the test, the micro-structure is in a fully structured state (λμ0=λν0=1).

Figure 3 shows an interesting aspect, for the loads γ̇0=102s1and 101s1the behavior is close to a purely viscoelastic response. These behaviors are in agreement with Figures 4 and 5, for the same loads. Note that these two shear rates correspond to a low level of modification in λμand λν. In these both cases, the stress of steady state is the maximum stress reached.

Figure 3.


Figure 4.


Figure 5.


It is important to point out that for γ̇0=1s1the thixotropic effect can ever be noted in stress response (Figure 3). The loads γ̇0=10s1and 100s1show a typical thixotropic behavior (Figure 3). These responses are related to considerable modifications in the structural nature of the substance, as can be seen in Figures 4 and 5. In these both cases, the maximum stress occurs before the steady state. It can be proved, in a theoretical/mathematical sense, that the point where the maximum breakdown rate occurs is before the stress peak. This fact is in agreement with expected behavior of thixotropic substance.

Other interesting result is related to the energy behavior of the thixotropic substance (Figures 6 and 7). Note the presence of a specific slope change in the region where higher rates of decrease on the micro-structures (λμ,λμ) levels occurred. The same behavior can be observed on experimental results. In fact, this can be explained, in a theoretical sense, via the following relationship


that comes from the rate Eqs. (11) and (12).

Figure 6.


Figure 7.


5.2. Constant stress test

This chapter presents some points and results for the constant shear stress test. In this sense, it is considered loads as τt=Htτ0, where Htrepresents the standard Heaviside function and τ0is a positive real constant. It is assumed that in the begin of the test, the micro-structure is in a fully structured state (λμ0=λν0=1).

It can be seen in Figure 8, for the loads τ0=10, 20, 30, 40, 50and 100Pa, a typical behavior for thixotropic substances under low level of constant stress loads. In this test, it can be seen the “Avalanche effect,” and their relation with the micro-structural level (λμ,λν) of the substance (Figures 9 and 10). In these cases, it can be observed the strain rate decrease, related to an increase of the construction parcel (cpμ) and decrease of the destruction parcel (dpμ) of λμ, as can be seen in Figures 11 and 12 respectively, where


Figure 8.


Figure 9.


Figure 10.


Figure 11.


Figure 12.


For the Maxwell element (.ν), in the same region (the strain rate decreases), cpνis close to null value and dpνis increasing, but this increase level is not sufficient to stop the decreasing process of λ̇. A posterior increase of γ̇is due to an abrupt decrease of linkages number (λμ,λν) before the steady state.

It is important to note the relationship between the structural nature and the behavior of the substance (Figures 9 and 10). Small changes in the values of λμand λνwere detected for stress loads 10and 20Pa, in relation to the others. In these two cases γ̇presents an nonincreasing behavior.

Other interesting result is presented in Figure 13, where it can be seen the energy behavior along the test. Note a changing on slope of the energy lines, tending to horizontal slope, in the γ̇decreasing region and another one related to the γ̇increasing region.

Figure 13.



6. Concluding remarks

The main objective of this chapter is to investigate the constitutive model for thixotropic fluids based on [8] related to the existence of two different types of micro-structure, and their consistence in some rheological tests. The model and analysis presented here are based on well-established physics principles. In this sense, it is important to stand out some nonstandard points of the approach presented in this work for thixotropic modeling in respect to the others models. In this sense, it is important to stand out some points:

  • the shear modulus (Gλν), the viscosity coefficients (ημλμand ηνλν) and their dependences from the two different micro-structure types are considered in the development of the constitutive model (Eqs. (9)(14));

  • the set of the rate equations (Eqs. (11)(12)) are related to well-established physical principles as the “reptation” model and the “Smoluchowski” theory of coagulation. It is important to point out that in the “Smoluchowski” theory of coagulation the effect of the Brownian motion is clearly taken into account;

  • the thermodynamic consistence of the model was analyzed (Section 3);

  • a theoretical criterium for the transition region based on the strain gradient mapping degeneration was discussed and exploited (Section 4), however it is important to stand out the necessity of more discussions and analysis on the relationships (Eqs. (48) and (56)) taking into account some additional characteristics (as temperature, …) and their consequences. These characteristics shall be considered in future work goals;

  • the illustrative numerical examples attest the capability of the model to predict the expected behavior of real thixotropic substances under some typical rheological tests (Section 5).

It is also important to comment that the developed ideas and presented in this work can be easily extended to approaches including more than two types of micro-structure (related to each specific relaxation mechanism) that are taken into account in the constitutive equation system. It is clear that the complexity level of the considered substance determines the necessity of these incorporations. In this context, it can be noted the versatility of the approach exposed here, presenting some interesting perspectives on thixotropic modeling that can be more explored in future works.


τshear stress
γshear strain
Gshear modulus
ηdynamic viscosity
λstructural parameter
Gγ̇functional associated to structural parameter time evolution
δr…i…the generalized Kronecker delta
·total time derivative
··double total time derivative
Oldroyd time derivative
νassociated to the Maxwellian mechanism
μassociated to the pure viscous mechanism
¯first-order tensor
¯¯second-order tensor
inner product
:double inner product

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Hilbeth P. Azikri de Deus and Mikhail Itskov (November 5th 2018). On the Thermodynamic Consistency of a Two Micro-Structured Thixotropic Constitutive Model, Selected Problems of Contemporary Thermomechanics, Jerzy Winczek, IntechOpen, DOI: 10.5772/intechopen.75987. Available from:

chapter statistics

517total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Practical Methods for Online Calculation of Thermoelastic Stresses in Steam Turbine Components

By Mariusz Banaszkiewicz and Janusz Badur

Related Book

The Development and Application of Microwave Heating

Edited by Wenbin Cao

First chapter

Microwave Applications in Thermal Food Processing

By Mohamed S. Shaheen, Khaled F. El-Massry, Ahmed H. El-Ghorab and Faqir M. Anjum

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us