Dimensionless Numbers in Chemical Engineering: Transport Phenomena Guide

1. Introduction to Dimensional Analysis in Chemical Engineering

Dimensional analysis and dimensionless numbers form the universal bedrock of chemical engineering and transport phenomena. Whether analyzing turbulent momentum transfer in high-pressure oil and gas pipelines, convective heat transfer inside shell-and-tube exchangers, or species mass transfer across packed distillation column packings, dimensionless numbers allow process engineers and researchers to compress complex multi-variable differential equations into elegant, physically interpretable scaling ratios.

According to the Buckingham π Theorem, any physically meaningful equation involving $n$ physical variables expressible in terms of $k$ independent fundamental dimensions can be reformulated as an equivalent relationship among $p = n - k$ dimensionless parameters (π-groups):

$$\Pi_1 = \Phi(\Pi_2, \Pi_3, \dots, \Pi_p)$$

In chemical engineering practice, these dimensionless ratios represent the relative magnitude of competing physical mechanisms???such as inertial forces versus viscous dissipation, convective heat transfer versus molecular conduction, or bulk fluid motion versus molecular diffusion. Understanding these balances is critical for plant scale-up, pilot plant prototyping, hydraulic line sizing, and reactor design.

2. Core Dimensionless Numbers: Mathematical Formulation & Physical Meaning

2.1 Reynolds Number ($Re$) : Fluid Dynamics

The Reynolds number represents the ratio of inertial forces to viscous shear forces within a moving fluid element:

$$Re = \frac{\text{Inertial Forces}}{\text{Viscous Forces}} = \frac{\rho v D}{\mu} = \frac{v D}{\nu}$$

Where:

• $\rho$ = Fluid density ($\text{kg/m}^3$)
• $v$ = Characteristic fluid velocity ($\text{m/s}$)
• $D$ = Characteristic length or pipe inside diameter ($\text{m}$)
• $\mu$ = Dynamic fluid viscosity ($\text{Pa}\cdot\text{s}$ or $\text{kg}/(\text{m}\cdot\text{s})$)
• $\nu = \mu / \rho$ = Kinematic viscosity or momentum diffusivity ($\text{m}^2/\text{s}$)

2.2 Prandtl Number ($Pr$) : Thermal Boundary Layers

Named after Ludwig Prandtl, the Prandtl number is a pure thermophysical property ratio of the fluid (independent of flow geometry or velocity). It quantifies the ratio of momentum diffusivity to thermal diffusivity:

$$Pr = \frac{\text{Kinematic Viscosity (Momentum Diffusivity)}}{\text{Thermal Diffusivity}} = \frac{\nu}{\alpha} = \frac{\mu / \rho}{k / (\rho C_p)} = \frac{C_p \mu}{k}$$

Where:

• $C_p$ = Specific heat capacity at constant pressure ($\text{J}/(\text{kg}\cdot\text{K})$)
• $\mu$ = Dynamic viscosity ($\text{Pa}\cdot\text{s}$)
• $k$ = Thermal conductivity ($\text{W}/(\text{m}\cdot\text{K})$)
• $\alpha = k / (\rho C_p)$ = Thermal diffusivity ($\text{m}^2/\text{s}$)

The magnitude of $Pr$ dictates the relative growth rate of the hydrodynamic boundary layer ($\delta$) versus the thermal boundary layer ($\delta_t$):

$$\frac{\delta}{\delta_t} \approx Pr^n \quad (n \approx 1/3 \text{ for laminar flow})$$

• Liquid Metals ($Pr \ll 1$, e.g., Liquid Sodium $Pr \approx 0.005$): Heat conducts exceptionally fast relative to viscous momentum. The thermal boundary layer extends far beyond the hydrodynamic boundary layer, ideal for fast-breeder nuclear reactors.
• Gases ($Pr \approx 0.7 - 0.8$, e.g., Air, Steam, Flue Gas): Momentum and thermal boundary layers grow at almost identical rates.
• Water ($Pr \approx 2 - 10$ at ambient to elevated temperatures): Hydrodynamic boundary layer is moderately thicker than thermal boundary layer.
• Heavy Lubricating Oils ($Pr \approx 100 - 10,000+$): Viscous momentum spreads rapidly throughout the channel, while heat conducts sluggishly through the core.

2.3 Schmidt Number ($Sc$) : Mass Boundary Layers

The Schmidt number is the mass transfer analogue of the Prandtl number. It represents the ratio of momentum diffusivity to molecular mass diffusivity:

$$Sc = \frac{\text{Kinematic Viscosity (Momentum Diffusivity)}}{\text{Mass Diffusivity}} = \frac{\nu}{\mathcal{D}_{AB}} = \frac{\mu}{\rho \mathcal{D}_{AB}}$$

Where $\mathcal{D}_{AB}$ is the binary molecular diffusion coefficient of solute $A$ diffusing through solvent $B$ ($\text{m}^2/\text{s}$). The Schmidt number controls the relative thickness of the velocity boundary layer ($\delta$) to the concentration boundary layer ($\delta_c$):

$$\frac{\delta}{\delta_c} \approx Sc^{1/3}$$

• Gas Systems ($Sc \approx 0.2 - 2.0$): Molecular mass diffusion is rapid; concentration boundary layers are comparable in thickness to momentum layers.
• Liquid Systems ($Sc \approx 500 - 2,000+$): Molecular diffusion in liquids is roughly 10,000 times slower than in gases ($\mathcal{D}_{AB} \sim 10^{-9} \text{ m}^2/\text{s}$). The concentration boundary layer is microscopically thin, creating steep concentration gradients near phase interfaces.

2.4 Lewis Number ($Le$) : Heat vs. Mass Diffusivity

The Lewis number compares thermal diffusivity directly to mass diffusivity:

$$Le = \frac{\alpha}{\mathcal{D}_{AB}} = \frac{Sc}{Pr} = \frac{k}{\rho C_p \mathcal{D}_{AB}}$$

In wet cooling tower design and air conditioning psychrometrics (water evaporating into unsaturated air), $Le \approx 1.0$. This thermodynamic coincidence forms the basis of the Lewis Relation, allowing direct linkage between convective heat transfer coefficients ($h$) and convective mass transfer coefficients ($k_c$):

$$h \approx k_c \rho C_p$$

2.5 Nusselt Number ($Nu$) & Sherwood Number ($Sh$)

In thermal and mass separation processes, surface transfer is characterized by:

• Nusselt Number ($Nu$): Ratio of convective heat transfer to pure conductive heat transfer across a fluid layer of characteristic thickness $L$: $$Nu = \frac{h L}{k}$$
• Sherwood Number ($Sh$): Ratio of convective mass transfer to pure molecular mass diffusion: $$Sh = \frac{k_c L}{\mathcal{D}_{AB}}$$

For fully developed turbulent flow inside smooth circular pipes, the classic Dittus-Boelter correlation and Gnielinski correlation link these numbers directly to $Re$ and $Pr$:

$$Nu = 0.023 \cdot Re^{0.8} \cdot Pr^n \quad (n = 0.4 \text{ for heating, } n = 0.3 \text{ for cooling})$$

2.6 Péclet Number ($Pe$) : Advection vs. Diffusion

The Péclet number measures the relative importance of bulk advective transport to molecular diffusion:

$$Pe_{thermal} = \frac{\text{Heat Advected}}{\text{Heat Conducted}} = \frac{v L}{\alpha} = Re \cdot Pr$$ $$Pe_{mass} = \frac{\text{Mass Advected}}{\text{Mass Diffused}} = \frac{v L}{\mathcal{D}_{AB}} = Re \cdot Sc$$

In continuous chemical reactors (tubular PFRs and packed beds), high Péclet numbers ($Pe > 100$) indicate ideal plug flow where axial dispersion is negligible, whereas low Péclet numbers ($Pe < 1$) indicate that axial back-mixing dominates.

3. Summary Comparison Table: Transport Phenomena Analogies

4. Step-by-Step Practical Worked Example

Design Scenario: A cooling water line in a petrochemical complex circulates treated cooling water at $T = 20^\circ ext{C}$ through a DN100 Schedule 40 carbon steel pipe (Inside Diameter $D = 102.3 ext{ mm} = 0.1023 ext{ m}$). The design liquid velocity is $v = 1.8 ext{ m/s}$.

Water Thermophysical Properties @ 20??C:

• Density $ ho = 998.2 ext{ kg/m}^3$
• Dynamic Viscosity $\mu = 1.002 ext{ cP} = 1.002 imes 10^{-3} ext{ Pa}\cdot ext{s}$
• Specific Heat $C_p = 4.182 ext{ kJ/(kg}\cdot ext{K)} = 4,182 ext{ J/(kg}\cdot ext{K)}$
• Thermal Conductivity $k = 0.598 ext{ W/(m}\cdot ext{K)}$

Step 1: Calculate Reynolds Number ($Re$)

$$Re = \frac{\rho v D}{\mu} = \frac{(998.2\text{ kg/m}^3)(1.8\text{ m/s})(0.1023\text{ m})}{1.002 \times 10^{-3}\text{ Pa}\cdot\text{s}} = \frac{183.81}{0.001002} = 183,443$$

Result: $Re = 1.83 \times 10^5 \gg 4,000$. Flow is fully turbulent.

Step 2: Calculate Prandtl Number ($Pr$)

$$Pr = \frac{C_p \mu}{k} = \frac{(4,182\text{ J/(kg}\cdot\text{K)})(1.002 \times 10^{-3}\text{ Pa}\cdot\text{s})}{0.598\text{ W/(m}\cdot\text{K)}} = \frac{4.1904}{0.598} = 7.01$$

Result: $Pr = 7.01$. The hydrodynamic boundary layer is approximately $Pr^{1/3} = 7.01^{0.333} \approx 1.91$ times thicker than the thermal boundary layer.

Step 3: Calculate Thermal Péclet Number ($Pe$)

$$Pe = Re \cdot Pr = (183,443)(7.01) = 1,285,935 \approx 1.29 \times 10^6$$

Result: Advective convective heat transport dominates axial conduction by over a million-fold.

Step 4: Estimate Inside Heat Transfer Coefficient ($h_i$) using Dittus-Boelter

$$Nu = 0.023 \cdot Re^{0.8} \cdot Pr^{0.4} = 0.023 \cdot (183,443)^{0.8} \cdot (7.01)^{0.4}$$ $$Nu = 0.023 \cdot (16,076) \cdot (2.176) \approx 804.7$$ $$h_i = \frac{Nu \cdot k}{D} = \frac{(804.7)(0.598\text{ W/m}\cdot\text{K})}{0.1023\text{ m}} = 4,704\text{ W/(m}^2\cdot\text{K)}$$

This computed value provides the exact film heat transfer coefficient needed for thermal rating in shell-and-tube heat exchangers.

5. Common Engineering Pitfalls & Scale-Up Rules of Thumb