Open access peer-reviewed chapter

Energy Transfer and Dissipation in non-Newtonian Flows in non-Circular Tubes

By Mario F. Letelier

Submitted: June 4th 2014Reviewed: November 12th 2014Published: July 8th 2015

DOI: 10.5772/59907

Downloaded: 1242

1. Introduction

In general, physiological processes are basically isothermal and isobaric. Several mechanisms contribute to keep constant pressure in systems such as blood vessels, among them diameter distribution and peristaltic motion. In this, it is to be understood that pressure constancy refers to time-averaged or extreme values.

Many pathologies make the human body develop fever, which raises temporarily body temperature. It is common practice to use, in some cases, medicines that restore, albeit for a given period of time, normal temperature. Fever is not the only cause of body heat-transfer. This process also occurs with body adjustment to changes in ambient temperature, physical exertion, metabolism acceleration due to food ingestion, drug effects, and others. During these body changes, blood and other fluid vessels undergo heat-transfer processes that pose considerable difficulties to physical modeling. Some pertinent variables are fluid composition, vessel elasticity, peristaltic motion in some cases, unsteadiness, and complex vessel geometry.

Heat-transfer in tube-flow has been investigated for several decades. One main practical driving force for such studies is the need to understand and design efficient heat-transfer equipment, widely used in several industrial areas. Models of simple tube-flow heat-transfer problems are now textbook´s standard contents.

These include laminar Newtonian steady flow in round tubes with simple boundary conditions. More recently, as it is reviewed in the next section, Newtonian and non-Newtonian flows in tubes of geometry other then circular have been modeled and analyzed. One important finding as to this chapter´s objective is the effect of the interplay between tube geometry and non-linear viscoelasticity.

Such coupling leads to the development of secondary flows, or helicoidal flows, that increase the transversal transport capacity of the flow, a phenomenon that can be applied, or related, to heat-transfer enhancement, or to other transversal transport processes, such as cross-sectional transport of particles immersed in the fluid.

Especially in blood flow, cells and other components introduce viscoelasticity [1, 2] that becomes relevant in smaller blood vessels. The understanding of its effect on the flow characteristics may become very relevant in cases such as heart arteries, and when artificial implants affect the blood flow. Also blood plastic effects appear is smaller vessels due to the aggregation of red blood cells at low shear rates, which develop a yield stress to be overcome for the flow to ensue.

Drawing on the above results, in this chapter it is presented a summary of relatively new analytical findings that may be useful for the better understanding of heat-transfer and complex flow phenomena in vessels that share some characteristics with biological vessels, particularly when these work under abnormal conditions. Also it is analyzed the effect of geometry in energy dissipation. This chapter is closely related to [3].


2. Mathematical models

Non-linear viscoelasticity and fluid plasticity are related in the following to heat- transfer and energy dissipation processes. Next are presented the mathematical models to be considered together with some remarks as to the corresponding state-of-the art, and relevant developments with related references.

2.1. Viscoelastic flow

The physical model considered is a straight tube of arbitrary cross-sectional shape, in which a non-Newtonian fluid moves along the axial coordinate zimpelled by a pressure gradient, which can be a function of t=time. Secondary flows are induced when the necessary conditions operate. The flow is assumed laminar, incompressible and with constant properties. Considering Cartesian coordinates, the velocity field can be expressed as (u,v,w)in which the velocity components align with the (x,y,z)axes respectively. The temperature and velocity fields are, in general, dependent on x,y,z and t.

Figure 1.

Definition diagram

Two well-known problems are here relevant, i.e.,

  1. Flow with constant wall heat-flux, in which the wall temperature is assumed constant at any cross-section, with prescribed axial variation.

  2. The “Graetz problem”, in which the flow is assumed isothermal up to a point z=0, from which on a different constant temperature is applied at the wall, so that the fluid is progressively heated or cooled.

Experimental findings concerning heat-transfer characteristics of aqueous polymer solutions flowing in straight tubes point at considerable enhancement as compared to its Newtonian counterpart driven by the same conditions and in the same geometry. Specifically it is reported that heat-transfer results for viscoelastic aqueous polymer solutions are considerably higher in flows fully developed both hydrodynamically and thermally, as much as by an order of magnitude depending primarily on the constitutive elasticity of the fluid and to some extent on the boundary conditions, than those found for water in laminar flow in rectangular ducts [4, 5]. Heat-transfer phenomena in laminar flow of non-linear fluids has not been the subject of many investigations with the exception of round pipes, and the case of inelastic shear-thinning fluids in tubes of rectangular cross-section, in spite of the widespread use of some specific contours in industry such as flattened elliptical tubes.

This statement is true for all cross-sectional shapes for both steady and unsteady phenomena including quasi-periodic flows. Heat-transfer with viscoelastic fluids has been declared to be a new challenge in heat-transfer research in the early nineties [6], but progress has been limited since that time. The physics of the phenomenon has not been entirely clarified.

Highly enhanced heat-transfer to aqueous solutions of polyacrylamide and polyethylene of the order of 40–45% as compared to the case of pure water in flattened copper tubes was observed by Oliver [7] and later by Oliver and co-workers as early as 1969. Recent numerical investigations in rectangular cross-sections of Gao and Hartnett [8, 9], Naccache and Souza Mendes [10], Payvar [11] and Syrjala [12] establish the connection between the enhanced heat-transfer observed and the secondary flows induced by viscoelastic effects. The former researchers as well as Naccache and Souza Mendes predict for instance viscoelastic Nusselt numbers as high as three times their Newtonian counterparts. Gao and Hartnett [8, 9].report numerical results in rectangular contours which provide evidence that the stronger the secondary flow (as represented by the dimensionless second normal stress coefficient Ψ2) the higher the value of the heat-transfer (as represented by the Nusselt number Nu) regardless the combination of thermal boundary conditions on the four walls. Constant heat flux is imposed everywhere on the heated walls in their numerical experiments with the remaining walls being adiabatic. The combination of boundary conditions plays some role in the enhancement reported with the largest enhancement occurring when two opposing walls are heated. Despite these efforts heat-transfer characteristics of viscoelastic fluids in steady laminar flow in rectangular tubes remains very much an open question (quoted from Siginer and Letelier [13]).

Coelho et al [14] presented an analytical solution for the Graetz problem for the MPTT fluid for a circular tube, steady flow, including several computations for negative heat flux (flow cooling). Valko [15] published a solution of the Graetz problem using a power-law fluid model; he determined the influence of the Brinkman number for several flow conditions in circular tubes. Kin and Özisik [17] published work on transient laminar forced convection of a power-law fluid in ducts with sudden change in wall temperature. In these and related references it is reflected the actual state-of-the-art in this subject.

Concerning numerical analysis, there is only one available commercial package for viscoelastic fluid flow computations, POLYFLOW. However POLYFLOW cannot handle even relatively high Weissenberg number flows and although convergent, gives erroneous results of the order of 400% as compared to analytical test cases, Filali et al. (2012). In addition it cannot handle heat-transfer in steady viscoelastic fluid flow in tubes. Thus to study problems of this type a numerical algorithm has to be built from scratch and tested for stability (Hadamard-type). Numerical analysis is not considered in this chapter.

Two well-known models of non-linear viscoelastic fluids are next described:

Modified Phan-Thien-Tanner (MPTT)

2ηmD=1+ϵ0ληm0tr ττ+λVτ-φτ-τφTE1


τ+λτˇ+αλη0ττ=η0Dτˇ=- VTτ+τVE2

In the following, the MPTT model is applied to the flow field analysis. For the purposes of this presentation, there is some evidence [17] that both models lead to qualitatively similar results.

The applicable equation of motion for 3-D, steady and incompressible flow are, in cylindrical coordinates:


V=0 E3



or, in expanded form


In the above, u, v, ware the radial, tangential and axial velocity components, σis the stress matrix and Pis the piezometric pressure. Scale factors applied are a (base radius) for r, w0for the velocity components, and ηN w0/afor the stress components, in which ηN is the Newtonian viscosity.

The MPTT model of viscoelastic fluid is next expressed, in dimensional variables, through the following equations, i.e. [18].

σ=-PI+2ηND+τ E8

where Dis the rate of deformation tensor and τis the non-Newtonian component of the shear stress.


φ=VT- ξDE11

then, a more compact form of (10), is

τ=Vτ-φτ- τφTE12

In this ηmis a viscosity, ε0 , ξare material parameters, and λis the relaxation time.

The function f, as defined in the MPTT model can be simplified for small values of ε0, yielding


where from

2ηmD=1+ε0ληm0tr ττ+λVτ-φτ-τφTE14

For this constitutive equation the viscosity is defined as


in which


and where ηm0, and mare additional material parameters. If m=1, then (15) reduces to


Coming back to dimensionless variables and parameters (14) becomes

21+2ξ2-ξWi2trD2D=1+ε0 Wi trττ+Wi τE18

in which Wiis the Weissenberg number, or dimensionless relaxation time, defined as


By combining equations (18), (8) and (5-7), it is obtained the set of working equations in full, i.e.


in which the functions Fr , Fθand Fzare the viscoelastic forcing functions that determine the transversal flow. For developed flow, all derivatives with respect to z, excepted for the pressure, must be put equal to zero.

The explicit expressions for the extra forcing terms are:

Fθ=τθ=τrθr+1rτθθθ+2rτrθr E24


For steady flow of an incompressible fluid with constant properties, the energy equation, in terms of the temperature T can be written as


Where ρ0is the density, Cpis the specific heat at constant pressure, ϕis the dissipation function, and kis the thermal conductivity coefficient. Further, it is assumed negligible dissipation and constant k. Under these assumptions (26) becomes


2.2. Viscoplastic flow

Viscoplastic fluids are fluids that exhibit yield stress, which must be overcome before the material develops deformation. As already mentioned plasticity appears in small blood vessels due to red cell aggregation.

Also, many industrial fluids exhibit yield stress, and are found in areas such as mining (slurries), food (pastes), construction (concrete and mud), cosmetics, etc. Flow and heat-transfer description in tubes and other configurations is compounded by their geometry, which brings in non-linear constitutive equations, except for very simple shapes. Some recent research in this field include flow around a cylinder by Tokpavi et al. [19] particle sedimentation, Yu and Wachs [20], flow in an eccentric annular tube, Wachs [21], and Walton and Bittleston [22],a general analysis for flow in non-circular ducts, Letelier and Siginer [23] and a preliminary analysis of the velocity field in non-circular pipes, Letelier, et al. [24]. Other pertinent references can be found in the previous ones.

The constitutive characteristics of viscoplastic flow determine complex structures of velocity and shear fields in tube flow when the tube cross-sectional contour differs from circular. This is true even for the case of the Bingham model of viscoplastic fluid, which is one of the simplest mathematical expressions for this kind of fluids.

Plasticity implies existence of fluid yield stress, which may induce both plug zones and stagnant zones within the tube cross-section, when in there the tube geometry determines zones of shear stress below the value of the yield stress. Such solid regions of the flow make it necessary to apply greater pressure gradients in order to keep a given rate of flow which, in turn, leads to greater energy dissipation and heating inside the fluid.

In the following Bingham´s model of fluid is used. Simple plastic flows in straight non-circular tubes do not develop secondary flows. The flow is, therefore, parallel when the motion is laminar, and only the axial component of the velocity exists. Under these conditions, the momentum equation, in terms of shear stress components, is


The constitutive relations are


in which


is an invariant related to the rate of deformation matrix. Also


where τyis the dimensional yield stress. The parameter Nis a dimensionless yield stress that greatly influence the flow characteristics.

The momentum and constitutive equations can be written in more compact form by using natural coordinates, i.e.


In the above, nis a coordinate normal to isovels and ρ=radius of curvature of isovels

The dissipation function ϕcan be used as an index of energy dissipation. For parallel flow ϕis given by


3. Analysis of secondary flows and their effect on heat-transfer

The equations of motion and energy are next solved by using a double perturbation method, as follows. First, and for all boundary conditions, the tube cross-section boundary is defined by


In which Gis here labeled as a “shape factor”, that can describe a wide array of shapes according to the given values of the parameters nand ε[1]. The shape perturbation parameter is ε, which can take values in between 0 and a limiting values for the curve (36) staying closed. The parameter nmust be given integer values in order to get regular shapes.

Next all velocity components are expanded in series in terms of the Weissenberg number, i.e.

u=Wi u1+Wi2u2+...v=Wiv1+Wi2v2+... w=w0+Wi w1+Wi2w2+...E37

Complementarily, all velocity component, at any order, are defined as

V=G f0+εf1+ε2f2+...E38

where Vis a generic velocity component, at any order in Wi, and f0,f1,...are functions specific for every V,to be determined by solving the momentum equations. Similarly, all other dependent variables are expressed in series, i.e.

τ=τ0+Wiτ1+Wi2τ2+... P=P0+WiP1+Wi2P2+...E39

More details can be found in references [13, 25]. The velocity field is first found by substituting (37-39) in (20-25) and related equations. For the axial velocity, the following expressions can be determined, i.e.

p=-14Pz   E41

Further, albeit more involved, exact solutions can be found for w3and higher order terms in Wi. In these results, since the parameter εhas a maximums value of 0,3849 for n=3, which decreases as nincreases, the series in εshown in brackets in (38) were developed up to (ε) yielding accurate results.

Computations show that up to 0Wi2the transversal velocity field is zero. Viscoelastic forcing terms in (20-22) are non-zero from 0 Wi3upwards, which implies u1=u2=v1=v2=0.

At third order inWi, the viscoelastic forcing terms become non-zero. These are


If the velocity components at third order inWiare expressed in terms of a stream-function Ψ3, i.e.


then it is found




The solution of (46) is


Where from


Plots of (47) are shown in figures 2, 3, 4 and 5

Figure 2.

Characteristics plots of transversal streamlines forn=3andε=0.3849

Figure 3.

Characteristics plots of transversal streamlines forn=3andε=0.33

Figure 4.

Characteristics plots of transversal streamlines forn= 4 andε=0.25

Figure 5.

Characteristics plots of transversal streamlines forn=4andε=0.2

It is to be noted that the vortical shape does not change with the slip parameter ξ, the pressure coefficient por the Weissenberg number. As these parameters change, the strength of the vortices change. Also, other parameters being equal, the strength of the vortices increases significantly with ε, and thus transversal transport capacity. In all these cases the Reynolds number is Re= 180. Similar results are found for n= 4, 5... in which the number of vortices is 2n.

In Figure 6 is shown the internal distribution of the normal axial shear stress for the case of ε=0.3849

Figure 6.

Characteristics plots of normal axial shear stress forn=3, ε=0.3849, Wi=0.3 andξ=0.2

Curves are similar for all values of parameters, but the values of the normal axial shear stress vary withWi, being of the order of 70% greater than the Newtonian counterpart for Re=180and Wi=0,3.

The temperature field can be computed once the velocity field is known, through the energy equation (27). To this end, the temperature Tis expressed as

T=T0+Wi T1+Wi2T2+...  E49

In the following the temperature field and heat-transfer are computed for the case in which there is a constant heat flux through de tube wall, so that the temperature difference between de wall and the average temperature, i.e., Tw-Taremains constant and also Ta/zis constant. This is problem 1 described in the section Mathematical Models.

The above leads to


in which


and Pris the Prandtl number, i.e.


Applying the condition that the temperature, at all orders in Wibe zero at the tube contour, i.e., that the fluid temperature is equal to Twat the wall, then it is found


Higher order terms of Tcan be found in similar way, but the expressions become very involved and are omitted here.

Heat exchange is computed through the Nusselt number Nu, i, e.

Nu=Dnhk  E54


h=kdTdndp-p-Tw-Ta  E55
Dn=4 Sp-=hydraulic diameter  E56
S=r dr dθ E57

In the above p-is de contour perimeter and Sis de cross-sectional area. A plot of Nuin terms of Wiis shown in Figure 7 for different values of the Reynolds number Re,defined as

Re=ρ0waDhηm0 E58

in which wais the axial average velocity.

Numerical values related to Figure 7 are given in the following table, where they can be compared with the Newtonian case

Typical plots of temperature distribution appears in Figures 8 and 9.

Figure 7.

Please add caption

Table 1.

Nusselt number values

Figure 8.

Characteristics plots of isothermal curves forn= 3 andε=0.3849

Figure 9.

Characteristics plots of isothermal curves forn=3andε=0.2

These results, as initially presented and discussed in [13, 25, 26] show that the Nusselt number, i.e., the heat-transfer between fluid and wall, increases as the viscoelastic parameter Wiincreases, with an asymptotic trend, for a given value of the Reynolds number. As Reincreases, also Nuincreases for a given value of Wi.

The relevance of these findings for biological flows may be related to vessel deformation due to wall elasticity and peristaltic motion.

The preceding analysis shows that a very small deviation of the cross-section contour from the circular geometry, as represented by the value of the parameter ε, may induce secondary flows, and so increase the transversal transport capacity of biological flows. It seems to be an open research area the study of physiological and morphological changes that, in this context, may induce pathologies that accelerate heat transfer inside the human body. Such changes may create conditions that improve or worsen transport processes that should lead to restore normal physiological states.

Previous results associated to the flow of viscoelastic fluid in channels of axially-varying cross-section [27] show also that axial change of geometry, as found in biological vessels, also augment transport capacity.

4. Plastic flow

The velocity field and related variables for the flow of a Bingham fluid in non-circular tubes is next found using a similar technique. The axial velocity is defined as

w =w0+εw1+ε2w2+..E59

The constitutive equations (29, 30) can be expanded in series around the parameter εby using (58), where from [24]


From the above the governing equation for w1are found, i.e.

r r-N22w1r2+1rw1r+2w1θ2=0 E62

which can be solved by separation of variables.

The dissipation function, in terms of the normal axial shear stress, is given in equation (35) and can be computed, as an indicator of energy dissipation due to friction.

Typical results for the velocity field and dissipation function are shown in the following figures.

Figure 10.

Plots of isovels forn=3,  ε=0.3849andN=0.2

Figure 11.

Plots of isovels forn=3, ε=0.3849andN=0.5

Figure 12.

Plots of isovels forn=3, ε=0.3849andN=1.0

Figure 13.

Plots of the dissipation function (35) forn=3, ε=0.3849andN=0.5

One significant feature of plastic flows are plug (or solid) and stagnant zones, which, as a general rule, tend to decrease the rate of flow for a given pressure gradient, or require that more energy by supplied in order to keep constant the rate of flow.

Figure 13 shows that energy dissipation is greatly affected by the tube geometry. The maximum energy dissipation in the case of an equilateral triangle appears in the middle of the sides, where the axial shear stress is maximum. Also the shear stress is zero at the corners, which is a factor that determines stagnant areas close to such corners.

The plug zones are defined through the following joint conditions

dwdn=0   ;    w=const E63

which can be met up to certain value of Naccording to the cross-section shape.

Abnormal vessel geometries may occur in biological flow arising from many sources. Especially in the case of stenosis, or geometry change due to solid deposition in artery walls, plastic effects may lead to blood clotting in corners, when the arteries are small.

5. Concluding remarks

Some analytical results concerning the effects of fluid elasticity and plasticity, coupled with tube cross-sectional geometry variation have been presented in this chapter.

The analysis is unified through some concepts that make it possible to explore, in a rather general fashion, the mechanisms of transversal transport that arise with the coupling of viscoelasticity and non-circular shape.

The circular shape is the most energy-efficient shape when only longitudinal mass transfer is considered. However, this not the case when transversal motion becomes important, as in heat-exchange processes, or in cases when also particle distribution is relevant.

General results indicate that viscoelasticity tends to increase transversal transport, which can be demonstrated by analytical means. These in particular, show that starting at third order in Wi, secondary flows appear when the fluid exhibits non-linear viscoelasticity, as prescribed by the fluid models of Giesekus and Phan-Tien-Tanner. Such secondary, or transversal motions, when coupled with the temperature, through the energy equation, determine a temperature field that improve heat-transfer between the fluid and the tube wall.

Plastic effects, on another perspective, do not induce transversal motion if the cross-section is not circular. Rather plasticity increases energy consumption when maintaining a given rate of flow is a priority. Shapes that include sharp corners lead to stagnant zones of fluids in the vicinity of those corners, thus decreasing the flow. The analytical method herein applied can be used for determining plug-zones of limiting values of the yield parameter Nfor which such zones exist properly.

The application of these findings to biological flows have been commented in previous sections of the chapter.



The author acknowledges the financial support of the National Fund for Scientific and Technological Development-FONDECYT, and of the Department of Scientific and Technological Research-DICYT at the University of Santiago, Chile.

© 2015 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

Mario F. Letelier (July 8th 2015). Energy Transfer and Dissipation in non-Newtonian Flows in non-Circular Tubes, Advances in Bioengineering, Pier Andrea Serra, IntechOpen, DOI: 10.5772/59907. Available from:

chapter statistics

1242total 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

Quantifying Stress in Crabs and Humans using Modified DFA

By Toru Yazawa

Related Book

First chapter

Enzyme-based Electrochemical Biosensors

By Zhiwei Zhao and Helong Jiang

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