Open access peer-reviewed chapter - ONLINE FIRST

Viscosity Models Based on the Free Volume and Entropy Scaling Theories for Pure Hydrocarbons over a Wide Range of Temperatures and Pressures

By Hseen O. Baled and Isaac K. Gamwo

Submitted: January 8th 2019Reviewed: April 12th 2019Published: October 7th 2019

DOI: 10.5772/intechopen.86321

Downloaded: 26

Abstract

Viscosity is a critical fundamental property required in many applications in the chemical and oil industry. Direct measurements of this property are usually expensive and time-consuming. Therefore, reliable predictive methods are often employed to obtain the viscosity. In this work, two viscosity models based on the free-volume and entropy scaling theories are assessed and compared for pure hydrocarbons. The modeling results are compared to experimental data of 52 pure hydrocarbons including straight-chain alkanes, branched alkanes, cycloalkanes, and aromatics. This study considers viscosity data to extremely high-temperature and high-pressure (HTHP) conditions up to 573 K and 300 MPa. The results obtained with the free-volume theory viscosity in conjunction with the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state are characterized by an overall average absolute deviation (AAD%) of 3% from the experimental data. The overall AAD% obtained with the predictive entropy scaling method by Lötgering-Lin and Gross is 8%.

Keywords

  • high temperature
  • high pressure
  • hydrocarbons
  • modeling
  • viscosity

1. Introduction

Viscosity is a key property in many engineering disciplines, including chemical and petroleum engineering. For instance, viscosity influences the fluid flow through porous media and pipelines; hence, it is required for the design of pipelines and transport equipment as well as for the estimation of recoverable oil and flow rates in porous media or wellbores. Viscosity can be determined through experimental measurements. However, carrying out viscosity measurements at all conditions of interest is not only expensive and time-consuming but also may not be possible at extreme conditions such as those encountered in ultra-deep reservoirs including pressures up to 300 MPa and temperatures up to 573 K. Reliable prediction models provide an alternative approach to generating predicted and correlated viscosity data at conditions where experimental data are not readily available.

Unlike the viscosity of gases at low pressures which is well defined by the kinetic theory of the gases, the viscosity theory of liquids is still inadequately developed due to the complications caused by the intermolecular forces between the molecules [1]. Therefore, there is no widely accepted simple theoretical method for predicting liquid viscosities, and most estimation techniques used for viscosity prediction of liquids are of empirical or semiempirical nature. The empirical models are correlations based on experimental observation with no theoretical background, whereas semi-theoretical models have a fundamental basis but contain adjustable parameters determined by fitting the model to experimental data.

In the present study, two viscosity models based on the free-volume and entropy scaling theories are assessed and tested against viscosity data for pure hydrocarbons from different chemical families that are commonly found in crude oil at high-temperature and high-pressure (HTHP) conditions up to 573 K and 300 MPa. Pure components are well-suited for initial evaluation of the viscosity models because a viscosity model that has difficulty in correctly describing the viscosity of a single hydrocarbon is likely to fail when predicting multicomponent mixtures.

2. Free-volume theory (FVT)

The FVT model is based on the free-volume concept. The idea that the viscosity depends upon the free space was first introduced by Batschinski [2] about 100 years ago. The viscosity, η, can be expressed as a sum of two contributions given in Eq. (1):

η=η0+ΔηE1

where η0 is the dilute gas viscosity and the Δη term dominates for liquid viscosity. The dilute gas term η0 is determined from the kinetic gas theory at very low pressures. It should be noted, though, that for liquids and supercritical fluids, the dilute gas viscosity term η0 is negligibly small in comparison to the total viscosity η; hence, η0 can be neglected for such fluids. Doolittle [3] found that the viscosity of liquid n-alkanes can be represented by a simple function of the free space fraction, fv=vfv0=vv0v0:

η=AexpB/fvE2

where v0is the molecular volume of reference or hard-core volume, vis the specific molecular volume, B is characteristic of the free-volume overlap, and A is a material-specific constant. A viscosity model based on the relation between free volume, friction coefficient, and viscosity has been proposed by Allal and coauthors [4, 5]:

η=AexpB/fvE3

where the free-volume fraction fvwas defined by means of the fluctuation-dissipation theory as

fvRTE32E4

where R is the gas constant and T is the temperature. In this expression, E=E0+PMρ, where E0=αρis related to the energy barrier that the molecule has to overcome in order to diffuse, ρis the density, and Pis the pressure. The viscosity of the dense state is linked to the fluid microstructure using the friction coefficient, ζ, which is related to molecular mobility and to the diffusion of linear momentum:

η=ρNAζL2ME5

where NAis Avogadro’s number and Lis a characteristic molecular length parameter. By combining the relations between free volume, friction coefficient, and viscosity, the following expression was obtained for the viscosity:

η=ρlαρ+PMρ3RTMexpBαρ+PMρRT3/2E6

The unit for the viscosity is [Pa·s], when all other variables are in SI units. The term PMρis linked to the energy necessary to form vacant vacuums required for the diffusion of the molecules. lis the characteristic length parameter in Å. The unitless parameter B is characteristic of the free-volume overlap. The density appears explicitly in Eq. (6), and hence the values of the free-volume theory parameters are directly dependent on whether experimental or calculated densities are used.

The three pure component parameters l, α, and B are determined by fitting Eq. (1) to experimental viscosity data. The use of the FVT requires density information, either experimental or calculated values. In this study, FVT is used in conjunction with the hybrid group-contribution perturbed-chain statistical associating fluid theory equation of state (G-C PC-SAFT EoS) [6] since this EoS provides reliable density predictions over wide ranges of pressure and temperature. In this equation, the PC-SAFT parameters are determined using two different sets of group-contribution (G-C) parameters for two pressure ranges: low-to-moderate pressures (≲7 MPa) and high pressures (≳7 MPa).

3. Entropy scaling model by Lötgering-Lin and Gross (ES-LG)

The basic idea of this method is to relate the viscosity to the residual entropy. The residual entropy is defined as the difference between a real state value and ideal gas state value at the same temperature and density, sresρT=sρTsidρT. Lötgering-Lin and Gross [7] proposed a predictive entropy scaling method for viscosities using a group-contribution (G-C) method based on the group-contribution perturbed-chain polar statistical associating fluid theory equation of state (G-C PCP-SAFT EoS) [8, 9]. Lötgering-Lin and Gross linked the Chapman-Enskog viscosity to PCP-SAFT segments in terms of G-C parameters, with

ηCE,gc=516MmolkBT/mgcNAπσgc2Ωgc22E7

where Mmol is the molar mass, kB is the Boltzmann constant, and T is the absolute temperature, NAis Avogadro’s number, mis the segment number, σis the segment diameter, and Ω22is the reduced collision integral. The index gc indicates pure component parameters that are calculated with the group-contribution method based on G-C-PCP-SAFT EoS.

A reduced viscosity is then defined as

η=ηηCE,gcE8

where the pure component reduced viscosity, ηi, is empirically correlated as

lnηi=Ai+Biz+Ciz2+Diz3E9

with

z=sreskBmgc,iE10

The residual entropy, sres, is calculated from the G-C-PCP-SAFT EoS originally proposed by Vijande et al. [8] and reparametrized by Sauer et al. [9]:

sresρT=aresTρE11

where ares=Ares/Nis the specific Helmholtz energy given by Gross and Sadowski [10]. Nis the total number of molecules.

The viscosity parameters Aito Diof pure substances are obtained from parameters Aαto Dαof functional group α, respectively. The following empirical expressions are proposed by Lötgering-Lin and Gross [7] for mixing group-contribution parameters:

Ai=αnα,imασα3AαE12
Bi=αnα,imασα3Vtot,iγBαE13
Ci=αnα,iCαE14
Di=Dαnα,iE15

With

Vtot,i=αnα,imασα3E16

where nα,idenotes the number of functional groups of type αin the substance i. The exponent γand the parameter D are kept constant for all studied substances and are optimized for n-alkanes (D = −0.01245 and γ= 0.45) [7]. The group-contribution parameters Aα, Bα, and Cαof all groups αare given in [7].

4. Modeling results

The two viscosity methods, FVT and ES-LG, are tested on a database consisting of 52 hydrocarbons (21 normal alkanes, 13 branched alkanes, 4 cycloalkanes, 14 aromatics) typically present in most of the crude oils from ambient conditions to extremely high-temperature and high-pressure (HTHP) conditions up to 573 K and 300 MPa. The temperature and pressure ranges considered in this study are given for each substance in Table 1. The performance of each model is assessed by the following statistical measures:

CompoundRanges of conditionsReferenceFVTES-LG
T/KP/MPaAAD/%MD/%Bias/%AAD/%MD/%Bias/%
Straight-chain alkanes (normal alkanes)
CH4298–5730.1–300[11, 12, 13]333−2154713
C2H6298–5730.1–70[11, 14]251735−6
C3H8298–5000.1–100[14, 15]5211620−2
n-C4H10298–5730.1–69[11]42025212
n-C5H12298–5730.1–252[16, 17]3181615−3
n-C6H14298–5730.1–300[16, 18]418−2519−1
n-C7H16298–5730.1–100[16, 19]2604174
n-C8H18298–5230.1–242[18, 20, 21]260272
n-C9H20298–4730.1–300[16, 19, 21]2904104
n-C10H22298–5730.1–300[16, 17, 21, 22, 23, 24]2604142
n-C11H24303–3230.1–62[19]0.11014−1
n-C12H26298–5730.1–300[16, 22, 23, 24, 25, 26, 27]31304172
n-C13H28303–3530.1–100[28]14148−4
n-C14H30313–3930.69–60[29]28−1512−2
n-C15H32310–4080.1–320[25]2704110
n-C16H34298–5340.1–273[30, 31]39010241
n-C17H36323–5730.1–0.1[3]26−2816−8
n-C18H38326–5340.1–280[25, 31]31209288
n-C19H40333–523[32]3804102
n-C20H42326–5341.38–243[31, 33]413012409
n-C32H66373–4580.3–0.3[34]120242
Branched alkanes
Isobutane300–5110.1–55[35, 36]3112313−1
Isopentane303–5730.098–196[16, 37]622−31123−7
Neopentane311–4440.7–55[38]723−31660−13
2-Methylpentane298–5500.1–300[16]1905161
3-Methylpentane3130.1–147[39]1418108
2,2-Dimethylbutane3130.1–147[39]28−12737−27
2,3-Dimethylbutane3130.1–147[39]0.310711−7
3-Ethylpentane3130.1–147[39]23−2172017
2,4-Dimethylpentane31301.147[39]261232923
2,2,4-Trimethylpentane298–5230.1–300[20, 40, 41, 42]2805131
2,3,4-Trimethylpentane298–4530.1–195[41]3110715−7
Squalane303–4731–202[43, 44]633−2264923
2,2,4,4,6,8,8-Heptamethylnonane298–4530.1–195[41, 45, 46]749−6751−5
Cycloalkanes
Cyclopentane298–3530.1–300[47, 48]516−226−1
Cyclohexane298–3930.1–100[49, 50]2704141
Methylcyclohexane298–3430.1–300[51, 52]160410−4
Ethylcyclohexane300–5301–50[53]312−1723−5
Aromatics
Benzene298–3730.1–300[54]130513−4
Toluene298–3730.1–299[55]260392
Ethylbenzene298–4530.1–195[41]241164016
Butylbenzene313–3730.1–100[56]2904114
Hexylbenzene313–3730.1–100[56]130595
Octylbenzene313–3730.1–100[56]150493
1,2-Diphenylethane353–4530.1–195[41]1201217−12
m-Xylene298–4730.1–199[20]260410−2
o-Xylene298–3480.1–110[57]1302530−25
p-Xylene298–3480.1–110[57]27−168−6
Naphthalene375–4540.1–101[58]83311122−11
1-Methylnaphthalene298–4730.1–200[20]1135−73260−32
Tetralin298–4480.1–201[20]520−1295328
Phenanthrene396–5730.1–101[58]1341−61729−17
Overall AAD%38

Table 1.

Performance of FVT and ES-LG models for pure hydrocarbons over wide ranges of pressure and temperature (entries are rounded to nearest whole number).

Absolute Average DeviationAAD=100Ni=1Nηi,calηi,expηi,expE17
Maximum DeviationMD=100·maxηi,calηi,expηi,expE18
Bias=100Ni=1Nηi,calηi,expηi,expE19

where N is the total number of data points, ηi,calrepresents the calculated viscosity value, and ηi,expis the experimental data point obtained from the literature. The absolute average deviation AAD is a measure of how close the calculated values are to the experimental data, while the bias indicates how well the calculated values are distributed around the literature data. Low values of the bias imply that the deviations are evenly distributed about zero. A positive bias indicates overestimation of the calculated viscosity, whereas a negative value indicates underestimation. These statistical measures of the ability of each of the selected seven models to reproduce viscosity values at HTHP conditions for each of the 52 pure compounds are given in Table 1.

The overall AAD obtained with the FVT model in conjunction with the hybrid G-C PC-SAFT EoS is 3%. The three adjustable parameters (l, α, B) required in this method are obtained by fitting the FVT predictions for each pure compound to the corresponding literature data. These optimized parameters yield reliable viscosity values over the whole ranges of temperatures to 573 K and pressures to 300 MPa.

The results obtained with the entropy scaling method by Lötgering-Lin and Gross are generally in good agreement with experimental data with an overall AAD of 8%. The AADs obtained for n-alkanes, branched alkanes, cycloalkanes, and aromatics are within 1–15, 3–27, 2–7, and 3–32%, respectively. These results are impressive for a fully predictive model that requires only the input of the molecular mass and the number of functional groups in each molecule. Unfortunately, this model cannot differentiate between isomers, such as 2-methylpentane and 3-methylpentane, and xylene isomers. In addition, this model has not yet been extended to binary, ternary, and multicomponent mixtures, such as crude oils.

For comparison purposes, Figure 1(a–d) shows the performance of the two studied viscosity methods, FVT and ES-LG, for four pure compounds representative of straight-chain alkanes (n-hexane), branched alkanes (2,2,4-trimethylpentane), cycloalkanes (methylcyclohexane), and aromatics (toluene).

Figure 1.

Viscosity predictions obtained with FVT and ES-LG viscosity models compared with literature data (EXP) for (a) n-hexane, (b) 2,2,4-trimethylpentane, (c) methylcyclohexane, and (d) toluene.

5. Conclusions

This work provides an assessment of the capabilities of two viscosity methods based on the free-volume and entropy scaling theories to model the viscosity of pure hydrocarbons over wide ranges of temperatures and pressures. The performance of the two studied viscosity models is discussed and evaluated by comparison to experimental viscosity data of 52 pure hydrocarbons from four different chemical families, namely, straight-chain alkanes, branched alkanes, cycloalkanes, and aromatics, at ambient and extremely high-temperature and high-pressure (HTHP) conditions up to 573 K and 300 MPa. The viscosity of pure components is required in most mixture models as an input parameter, and hence accurate and reliable model for pure compounds, particularly under high-pressure conditions, is a prerequisite for the mixture viscosity to be accurately estimated using the mixture model. The predictive entropy scaling method proposed by Lötgering-Lin and Gross (ES-LG model) predicts the viscosity with an overall absolute average deviation of about 8%, and the predictions are reasonable for most engineering and industrial applications given that the accuracy of most experimental viscosity data is within 1–5%. The free-volume theory (FVT) viscosity model provides very satisfactory results with an overall AAD of 3%. However, it is important to note that unlike the entropy scaling method, the free-volume theory is not a predictive model and requires that sufficient experimental viscosity data are readily available over the temperature and pressure ranges of interest to determine the fluid specific parameters.

Acknowledgments

The Open Access Publishing fee for this chapter was fully paid by the University Library System, University of Pittsburgh, USA.

Download

chapter PDF

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

Hseen O. Baled and Isaac K. Gamwo (October 7th 2019). Viscosity Models Based on the Free Volume and Entropy Scaling Theories for Pure Hydrocarbons over a Wide Range of Temperatures and Pressures [Online First], IntechOpen, DOI: 10.5772/intechopen.86321. Available from:

chapter statistics

26total chapter downloads

More statistics for editors and authors

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

Access personal reporting

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