15.1 Fluid Density, Specific Gravity, and Viscosity Correlations

15.1 Fluid Density, Specific Gravity, and Viscosity Correlations

Accurate evaluation of fluid density and viscosity is essential for solving core chemical engineering problems on the NCEES PE Chemical Exam, including line sizing, pump and compressor power calculations, frictional head loss determination via the Moody diagram, and boundary layer transport analysis. Chemical process streams span pure solvents, cryogenic liquefied gases, heavy crude oils, non-Newtonian polymer melts, and high-pressure supercritical gases.

---

1. Fluid Density and Specific Gravity Formulations

Mass Density, Specific Volume, and Specific Gravity

Mass density ($\rho$) is defined as mass per unit volume ($\text{kg/m}^3$ or $\text{lbm/ft}^3$). Its reciprocal is specific volume ($v = 1/\rho$, with units $\text{m}^3/\text{kg}$ or $\text{ft}^3/\text{lbm}$). On a molar basis, molar density is $\rho_m = \rho / M$ and molar volume is $V = M / \rho$, where $M$ is the molecular weight.

Specific gravity ($SG$) is a dimensionless ratio comparing the density of a target fluid to the density of a standard reference fluid at specified reference temperatures:

$$SG = \frac{\rho_{\text{fluid}}(T)}{\rho_{\text{ref}}(T_{\text{ref}})}$$

• Liquids (Standard Reference): Pure liquid water at its maximum density point ($4.0^\circ\text{C}$ or $39.2^\circ\text{F}$ and $1.0\text{ atm}$): $$\rho_{\text{water}}(4^\circ\text{C}) = 1000.0\text{ kg/m}^3 = 1.0000\text{ g/cm}^3 = 62.428\text{ lbm/ft}^3 = 8.3454\text{ lb/gal}$$
• Petroleum Industry Standard ($SG_{60/60^\circ\text{F}}$): In petroleum refining and custody transfer, specific gravity is evaluated at $60.0^\circ\text{F}$ ($15.56^\circ\text{C}$) relative to water at $60.0^\circ\text{F}$: $$\rho_{\text{water}}(60^\circ\text{F}) = 999.01\text{ kg/m}^3 = 62.366\text{ lbm/ft}^3 = 8.3371\text{ lb/gal}$$
• Gases (Standard Reference): Dry air at $0^\circ\text{C}$ ($32^\circ\text{F}$) or $60^\circ\text{F}$ and $1.0\text{ atm}$ ($M_{\text{air}} = 28.964\text{ kg/kmol}$ or $\text{lbm/lbmol}$): $$SG_{\text{gas}} = \frac{M_{\text{gas}}}{M_{\text{air}}}$$

API Gravity ($^\circ\text{API}$) and Baumé Scales ($^\circ\text{Bé}$)

In petroleum refining, the American Petroleum Institute defined the API Gravity scale to characterize crude fractions:

$$^\circ\text{API} = \frac{141.5}{SG_{60/60^\circ\text{F}}} - 131.5 \iff SG_{60/60^\circ\text{F}} = \frac{141.5}{^\circ\text{API} + 131.5}$$

Notice the strict inverse mathematical relationship: as fluid density increases, API gravity decreases.

• Water ($SG = 1.000$) corresponds exactly to $10.0^\circ\text{API}$.
• Light crudes ($> 31.1^\circ\text{API}$) have $SG < 0.870$ (e.g., gasoline $\approx 60^\circ\text{API}$, kerosene $\approx 42^\circ\text{API}$).
• Heavy crudes ($< 22.3^\circ\text{API}$) have $SG > 0.920$.
• Bitumen and extra-heavy oils have $^\circ\text{API} < 10.0^\circ$, meaning $SG > 1.000$, causing them to sink in water.

For inorganic aqueous process solutions (e.g., sulfuric acid, caustic soda, brine), the Baumé scale ($^\circ\text{Bé}$) is utilized:

$$\text{For liquids heavier than water } (SG > 1.0): \quad ^\circ\text{Bé} = 145 - \frac{145}{SG_{60/60^\circ\text{F}}} \iff SG = \frac{145}{145 - ^\circ\text{Bé}}$$ $$\text{For liquids lighter than water } (SG < 1.0): \quad ^\circ\text{Bé} = \frac{140}{SG_{60/60^\circ\text{F}}} - 130 \iff SG = \frac{140}{130 + ^\circ\text{Bé}}$$

Saturated Liquid Density: The Rackett Equation

For saturated pure liquids, the Rackett equation provides high-precision estimates of molar volume from the triple point up to reduced temperatures of $T_r \approx 0.95$:

$$V_{\text{sat}} = \frac{R T_c}{P_c} Z_{RA}^{\left[ 1 + (1 - T_r)^{2/7} \right]}$$

Where:

• $V_{\text{sat}}$ = saturated liquid molar volume ($\text{m}^3/\text{kmol}$ or $\text{cm}^3/\text{mol}$).
• $T_c, P_c$ = critical temperature ($\text{K}$) and critical pressure ($\text{kPa}$ or $\text{bar}$).
• $T_r = T / T_c$ = reduced temperature.
• $Z_{RA}$ = pure-component Rackett compressibility parameter (tabulated in the NCEES PE Chemical Reference Handbook; if unavailable, approximate with critical compressibility factor $Z_c = P_c V_c / (R T_c)$).
• $R$ = universal gas constant ($8.3145\text{ kPa}\cdot\text{m}^3/(\text{kmol}\cdot\text{K})$).

The saturated liquid mass density is calculated directly from the molar volume:

$$\rho_{\text{sat}} = \frac{M}{V_{\text{sat}}}$$

Ideal and Real Gas Density

Gas density is calculated using the equation of state:

$$\rho_{\text{gas}} = \frac{P M}{Z R T}$$

Where $Z$ is the compressibility factor ($Z = 1.0$ for ideal gases). For real gases at elevated pressures, $Z$ is obtained from generalized compressibility charts using reduced properties ($T_r = T/T_c, P_r = P/P_c$) or cubic equations of state (Peng-Robinson, SRK).

---

2. Viscosity Fundamentals: Dynamic vs. Kinematic

Newton's Law of Viscosity

Viscosity quantifies internal resistance to shear-induced deformation. In laminar parallel flow, shear stress ($\tau$) is directly proportional to the velocity gradient (shear rate $\dot{\gamma}$):

$$\tau_{yx} = -\mu \frac{dv_x}{dy} = \mu \dot{\gamma}$$

Where:

• $\tau_{yx}$ = shear stress parallel to the $x$-plane ($\text{N/m}^2$ or $\text{Pa}$, $\text{lbf/ft}^2$).
• $\dot{\gamma} = |dv_x/dy|$ = shear rate ($\text{s}^{-1}$).
• $\mu$ = dynamic (absolute) viscosity.

          Moving Plate (Area A, Velocity V) ------> Force F
          ===================================================
          Fluid Layer: Velocity Gradient dv/dy (Shear Rate)
          ---------------------------------------------------
          Stationary Plate (v = 0)
          ===================================================
          Shear Stress: tau = F / A = mu * (dv/dy)

Dynamic Viscosity Units and Conversions

• SI Unit: $\text{Pa}\cdot\text{s} = \text{kg}/(\text{m}\cdot\text{s}) = \text{N}\cdot\text{s/m}^2$.
• CGS Unit: $\text{Poise (P)} = \text{g}/(\text{cm}\cdot\text{s}) = 0.10\text{ Pa}\cdot\text{s}$. The industrial sub-unit is the centipoise (cP): $$1.0\text{ cP} = 10^{-2}\text{ P} = 10^{-3}\text{ Pa}\cdot\text{s} = 1.0\text{ mPa}\cdot\text{s}$$ (Convenient reference: Liquid water at $20.0^\circ\text{C}$ has $\mu = 1.002\text{ cP} \approx 1.0\text{ cP}$).
• US Customary Units: $\text{lbm}/(\text{ft}\cdot\text{s})$, $\text{lbm}/(\text{ft}\cdot\text{hr})$, and $\text{lbf}\cdot\text{s/ft}^2$: $$1.0\text{ cP} = 6.7197 \times 10^{-4}\text{ lbm}/(\text{ft}\cdot\text{s}) = 2.4191\text{ lbm}/(\text{ft}\cdot\text{hr}) = 2.0885 \times 10^{-5}\text{ lbf}\cdot\text{s/ft}^2$$

Kinematic Viscosity ($\nu$)

Kinematic viscosity is the ratio of dynamic viscosity to fluid density, representing the momentum diffusivity of the fluid:

$$\nu \equiv \frac{\mu}{\rho}$$

• SI Unit: $\text{m}^2/\text{s}$.
• CGS Unit: $\text{Stokes (St)} = \text{cm}^2/\text{s} = 10^{-4}\text{ m}^2/\text{s}$. The standard sub-unit is the centistokes (cSt): $$1.0\text{ cSt} = 10^{-2}\text{ St} = 10^{-6}\text{ m}^2/\text{s} = 1.0\text{ mm}^2/\text{s}$$ (Reference: Pure water at $20.0^\circ\text{C}$ has $\nu = 1.002 / 0.9982 \approx 1.004\text{ cSt}$).
• US Customary Unit: $\text{ft}^2/\text{s}$: $$1.0\text{ cSt} = 1.0764 \times 10^{-5}\text{ ft}^2/\text{s}$$
• Saybolt Universal Seconds (SUS): In older petroleum specifications, efflux time from a calibrated Saybolt viscometer is reported in SUS ($t$ in seconds at $100^\circ\text{F}$ or $210^\circ\text{F}$): $$\nu\text{ (cSt)} = 0.226 t - \frac{195}{t} \quad (\text{for } 32 \le t \le 100\text{ SUS})$$ $$\nu\text{ (cSt)} = 0.220 t - \frac{135}{t} \approx 0.216 t \quad (\text{for } t > 100\text{ SUS})$$

---

3. Rheological Classifications: Newtonian vs. Non-Newtonian Fluids

Fluids are classified by how their shear stress responds to imposed shear rate:

   Shear Stress (tau)
     ^
     |        / Dilatant (Shear-Thickening, n > 1)
     |       /   
     |      /      / Newtonian (n = 1, Constant mu)
     |     /      / 
     |    /      /        / Pseudoplastic (Shear-Thinning, n < 1)
     |   /      /        / 
     |  /      /        /
   tau0 +-----+--------+---------- Bingham Plastic (Yield Stress tau0)
     | /     /        /
     |/     /        /
     +-----------------------------------> Shear Rate (gamma_dot)

1. Newtonian Fluids

Newtonian fluids have a constant dynamic viscosity independent of shear rate or duration of shear. Examples include water, air, light hydrocarbons, simple alcohols, and benzene.

2. Ostwald-de Waele Power Law Model

For non-Newtonian fluids lacking yield stress, the power law model relates shear stress to shear rate:

$$\tau = K \dot{\gamma}^n$$

Where:

• $K$ = flow consistency index ($\text{Pa}\cdot\text{s}^n$ or $\text{lbf}\cdot\text{s}^n/\text{ft}^2$).
• $n$ = flow behavior index (dimensionless).

The apparent viscosity ($\mu_{\text{app}}$) is the effective viscosity at a specific local shear rate:

$$\mu_{\text{app}} = \frac{\tau}{\dot{\gamma}} = K \dot{\gamma}^{n-1}$$

• Pseudoplastic (Shear-Thinning, $n < 1$): Apparent viscosity decreases as shear rate increases. Molecular entanglements uncoil and align parallel to the flow streamlines. Examples: polymer melts, latex paint, fermentation broths, blood, paper pulp.
• Dilatant (Shear-Thickening, $n > 1$): Apparent viscosity increases as shear rate increases. High shear disrupts the packing of concentrated suspensions, forcing fluid into void spaces and increasing solid-solid interparticle friction. Examples: concentrated cornstarch suspensions (oobleck), wet beach sand, high-solids titanium dioxide slurries.

3. Bingham Plastics

Bingham plastics behave as rigid solids until an external yield stress ($\tau_0$) is exceeded, after which they flow with a constant plastic viscosity ($\mu_p$ or $\mu_\infty$):

$$\tau = \tau_0 + \mu_p \dot{\gamma} \quad (\text{for } |\tau| > \tau_0)$$ $$\dot{\gamma} = 0 \quad (\text{for } |\tau| \le \tau_0)$$

Examples include drilling muds, sewage sludge, toothpastes, greases, and mayonnaise.

4. Time-Dependent Rheology

• Thixotropic: Apparent viscosity decreases reversibly over time under constant, sustained shear rate (e.g., non-drip paints, drilling muds after standing).
• Rheopectic: Apparent viscosity increases reversibly over time under constant, sustained shear rate (e.g., gypsum slurries, bentonite clay pastes).

---

4. Temperature and Pressure Dependencies of Viscosity

Liquid Viscosity vs. Temperature: The Andrade / Arrhenius Equation

In liquids, molecules are closely packed in mutual potential energy wells. Momentum transfer occurs via short-range attractive intermolecular forces and molecular cage jumps into adjacent voids (free volume). Rising temperature increases thermal expansion, expanding intermolecular distances and providing thermal activation energy to overcome attractive cages. Consequently, liquid viscosity decreases exponentially with increasing temperature:

$$\ln(\mu_L) = A + \frac{B}{T} \iff \mu_L(T) = A' \exp\left( \frac{E_a}{R T} \right)$$

Where $T$ is absolute temperature ($\text{K}$ or $^\circ\text{R}$) and $E_a$ is the activation energy for viscous flow. For broad temperature spans, the three-parameter Vogel-Fulcher-Tammann (VFT) correlation is applied:

$$\ln(\mu_L) = A + \frac{B}{T - T_0}$$

Gas Viscosity vs. Temperature: Kinetic Theory & Sutherland's Law

In low-pressure gases, molecules are separated by distances much greater than their molecular diameters. Intermolecular attractive forces are negligible. Momentum transport occurs exclusively by the physical translation of gas molecules carrying momentum across parallel flow streamlines through collisions. Kinetic theory shows that the mean molecular thermal velocity is:

$$\bar{v} = \sqrt{\frac{8 R T}{\pi M}} \propto \sqrt{T}$$

Because hotter gas molecules travel faster and collide more frequently, gas viscosity increases with increasing temperature:

$$\mu_G \propto \sqrt{M T}$$

According to Chapman-Enskog kinetic theory for dilute gases:

$$\mu_G = 2.6693 \times 10^{-5} \frac{\sqrt{M T}}{\sigma^2 \Omega_v} \quad (\text{poise})$$

Where $\sigma$ is the collision diameter ($\text{Å}$) and $\Omega_v$ is the dimensionless collision integral (which varies weakly with $k_B T / \epsilon$).

For engineering calculations across wide temperature ranges, Sutherland's formula incorporates intermolecular attractions via an empirical Sutherland constant ($S$ in $\text{K}$):

$$\mu_G(T) = \mu_0 \left( \frac{T}{T_0} \right)^{3/2} \frac{T_0 + S}{T + S}$$

Pressure Dependence of Viscosity

• Liquids: Because liquids are virtually incompressible, viscosity is nearly independent of pressure from vacuum up to $20-50\text{ bar}$. At extreme pressures ($> 100-1000\text{ bar}$), free volume is compressed, causing viscosity to increase exponentially.
• Gases: At low to moderate pressures ($P < 10\text{ bar}$), gas viscosity is independent of pressure (as predicted by Maxwell, because gas density increases linearly while molecular mean free path decreases proportionally, cancelling out). At high pressures ($P_r > 1$), dense gases and supercritical fluids experience molecular clustering, causing viscosity to increase sharply toward liquid-like magnitudes.

---

5. Summary Comparison Table: Density and Viscosity Principles

---

6. Comprehensive Worked Numerical Example

Problem Statement

A chemical facility processes an organic process stream and a drilling mud slurry. Solve the following three engineering evaluations:

1. Crude Fraction Characterization: A petroleum condensate has an API gravity of $32.0^\circ\text{API}$ at $60.0^\circ\text{F}$. Determine its specific gravity $SG_{60/60^\circ\text{F}}$, mass density $\rho$ in $\text{lbm/ft}^3$, and mass flow rate $\dot{m}$ in $\text{lbm/hr}$ for a pipeline transporting $500.0\text{ gpm}$.
2. Rackett Saturated Liquid Density: Saturated liquid isobutane ($M = 58.12\text{ kg/kmol}$) is stored at $300.0\text{ K}$. Given $T_c = 407.8\text{ K}$, $P_c = 3640\text{ kPa}$, and $Z_{RA} = 0.2829$, calculate the saturated molar volume $V_{\text{sat}}$ in $\text{m}^3/\text{kmol}$ and mass density $\rho_{\text{sat}}$ in $\text{kg/m}^3$.
3. Power-Law Slurry Flow: A shear-thinning mineral slurry flows through a $D = 0.100\text{ m}$ ID pipe at an average velocity of $V_{\text{avg}} = 1.50\text{ m/s}$. The slurry obeys the Ostwald-de Waele model with $K = 0.800\text{ Pa}\cdot\text{s}^n$ and $n = 0.600$. Determine the wall shear rate $\dot{\gamma}_w$, wall shear stress $\tau_w$, and apparent viscosity $\mu_{\text{app},w}$ at the pipe wall in $\text{cP}$.

---

Step 1: Petroleum Condensate Characterization

From the API gravity definition:

$$SG_{60/60^\circ\text{F}} = \frac{141.5}{^\circ\text{API} + 131.5} = \frac{141.5}{32.0 + 131.5} = \frac{141.5}{163.5} = \mathbf{0.86544}$$

Mass density at $60^\circ\text{F}$:

$$\rho = SG \times \rho_{\text{water}}(60^\circ\text{F}) = 0.86544 \times 62.366\text{ lbm/ft}^3 = \mathbf{53.974\text{ lbm/ft}^3}$$

Volumetric flow rate conversion ($1\text{ ft}^3 = 7.4805\text{ gal}$):

$$Q = 500.0\text{ gpm} \times \frac{1\text{ ft}^3/\text{s}}{448.83\text{ gpm}} = 1.1140\text{ ft}^3/\text{s} = 66.840\text{ ft}^3/\text{min}$$ $$\dot{m} = Q \cdot \rho = (1.1140\text{ ft}^3/\text{s}) \times (53.974\text{ lbm/ft}^3) \times 3600\text{ s/hr} = \mathbf{2.1645 \times 10^5\text{ lbm/hr}}$$

---

Step 2: Saturated Isobutane via Rackett Equation

Compute reduced temperature:

$$T_r = \frac{T}{T_c} = \frac{300.0\text{ K}}{407.8\text{ K}} = 0.73565$$ $$1 - T_r = 1 - 0.73565 = 0.26435$$

Evaluate the exponent bracket:

$$(1 - T_r)^{2/7} = (0.26435)^{0.285714} = 0.68377$$ $$\text{Exponent} = 1 + (1 - T_r)^{2/7} = 1 + 0.68377 = 1.68377$$

Evaluate $Z_{RA}^{\text{Exponent}}$:

$$Z_{RA}^{1.68377} = (0.2829)^{1.68377} = \exp(1.68377 \times \ln(0.2829)) = \exp(1.68377 \times (-1.26267)) = \exp(-2.12604) = 0.11931$$

Calculate critical molar volume ratio:

$$\frac{R T_c}{P_c} = \frac{(8.3145\text{ kPa}\cdot\text{m}^3/(\text{kmol}\cdot\text{K})) \times (407.8\text{ K})}{3640\text{ kPa}} = \frac{3390.65}{3640} = 0.93150\text{ m}^3/\text{kmol}$$

Compute saturated molar volume and mass density:

$$V_{\text{sat}} = 0.93150 \times 0.11931 = \mathbf{0.11114\text{ m}^3/\text{kmol}} = 111.14\text{ cm}^3/\text{mol}$$ $$\rho_{\text{sat}} = \frac{M}{V_{\text{sat}}} = \frac{58.12\text{ kg/kmol}}{0.11114\text{ m}^3/\text{kmol}} = \mathbf{522.9\text{ kg/m}^3}$$

---

Step 3: Non-Newtonian Slurry Pipe Hydraulics

For fully developed laminar pipe flow of a power-law fluid, the wall shear rate is given by the Rabinowitsch-Mooney equation:

$$\dot{\gamma}_w = \left( \frac{3n + 1}{4n} \right) \left( \frac{8 V_{\text{avg}}}{D} \right)$$ $$\frac{3n + 1}{4n} = \frac{3(0.600) + 1}{4(0.600)} = \frac{2.800}{2.400} = 1.1667$$ $$\frac{8 V_{\text{avg}}}{D} = \frac{8 \times 1.50\text{ m/s}}{0.100\text{ m}} = 120.0\text{ s}^{-1}$$ $$\dot{\gamma}_w = 1.1667 \times 120.0\text{ s}^{-1} = \mathbf{140.0\text{ s}^{-1}}$$

Compute wall shear stress:

$$\tau_w = K (\dot{\gamma}_w)^n = 0.800 \times (140.0)^{0.600} = 0.800 \times 19.394 = \mathbf{15.515\text{ Pa}}$$

Compute apparent viscosity at the pipe wall:

$$\mu_{\text{app},w} = \frac{\tau_w}{\dot{\gamma}_w} = \frac{15.515\text{ Pa}}{140.0\text{ s}^{-1}} = 0.11082\text{ Pa}\cdot\text{s} = \mathbf{110.8\text{ cP}}$$

(Notice that without shear, the fluid would be much thicker; at $\dot{\gamma} = 1.0\text{ s}^{-1}$, $\mu_{\text{app}} = 0.800\text{ Pa}\cdot\text{s} = 800\text{ cP}$. Shearing reduces viscosity by over $86\%$).

---

7. Critical PE Exam Traps & Pitfalls