Surface Tension Interactive Calculator

Surface tension is the elastic tendency of liquid surfaces to minimize their surface area, arising from cohesive forces between liquid molecules. This interactive calculator enables engineers, scientists, and students to analyze capillary rise, contact angles, Young-Laplace pressure differences, and force balance in systems involving liquid interfaces. Applications span microfluidics, coating processes, ink-jet printing, flotation, and biological systems where interfacial phenomena dominate flow behavior.

📐 Browse all free engineering calculators

Visual Diagram: Surface Tension Phenomena

Interactive Surface Tension Calculator

Governing Equations

Theory & Practical Applications

Molecular Origin of Surface Tension

Surface tension arises from the imbalance of intermolecular forces at a liquid-air (or liquid-vapor) interface. Molecules in the bulk liquid experience cohesive forces equally in all directions, resulting in a net force of zero. However, molecules at the surface lack neighbors on the air side, creating a net inward cohesive force. This imbalance causes the interface to behave as if under tension, minimizing surface area to reduce the number of molecules at the energetically unfavorable interface. The energy required to increase surface area by one unit is the surface tension γ, typically expressed in N/m or equivalently J/m².

For water at 20°C, γ ≈ 0.0728 N/m. For mercury, γ ≈ 0.486 N/m, reflecting its much stronger cohesive metallic bonding. Organic solvents like ethanol have lower surface tensions (~0.022 N/m at 20°C) due to weaker intermolecular forces. Surface tension decreases with temperature as increased thermal energy disrupts intermolecular bonds, following approximately linear decline for most liquids until near the critical point where liquid-vapor distinction vanishes.

Capillary Action and Jurin's Law Derivation

Capillary rise in a narrow tube occurs when the upward surface tension force exceeds the weight of the liquid column. At equilibrium, the vertical component of surface tension force around the tube perimeter balances the gravitational force on the liquid column. The surface tension force acts along the contact line with magnitude F_γ = γ (2πr), and its vertical component is F_γ,vertical = γ (2πr) cos θ, where θ is the contact angle measured through the liquid.

The weight of the risen liquid column is W = ρ g π r² h. At equilibrium, γ (2πr) cos θ = ρ g π r² h. Solving for h yields Jurin's law: h = (2γ cos θ) / (ρ g r). This inverse dependence on radius explains why capillary rise is dramatic in very narrow tubes. A 0.5 mm diameter glass tube filled with water (θ ≈ 0°) exhibits approximately 30 mm of rise, whereas a 5 mm tube shows only 3 mm. This r⁻¹ scaling is critical in microfluidic design where channel dimensions approach tens of micrometers, producing capillary pressures exceeding 1000 Pa.

An important non-obvious limitation: Jurin's law assumes quasi-static equilibrium and neglects dynamic effects like viscous dissipation during the rise process. For rapid capillary filling, the Lucas-Washburn equation governs the time-dependent position: h²(t) = (γ r cos θ / 2η) t, where η is dynamic viscosity. This reveals that capillary rise velocity decreases as the column grows, asymptotically approaching the Jurin height. For a 100 μm radius tube with water (η = 1 mPa·s), the initial rise rate is approximately 10 mm/s, but slows to 1 mm/s after 5 mm of rise. In high-speed inkjet printing, these dynamic effects dominate over equilibrium predictions.

Young-Laplace Equation and Curved Interfaces

The Young-Laplace equation quantifies the pressure jump across a curved interface due to surface tension. For a surface with principal radii of curvature R₁ and R₂, the pressure difference is ΔP = γ (1/R₁ + 1/R₂). This equation is fundamental to understanding droplet formation, bubble stability, and emulsion behavior. A spherical water droplet of radius 1 μm in air experiences an internal overpressure of ΔP = 2γ/R ≈ 145,600 Pa (~1.44 bar), explaining why small droplets evaporate faster — the curvature-induced pressure increases vapor pressure according to the Kelvin equation.

For a cylindrical meniscus in a capillary tube, one radius of curvature is the tube radius r (in the radial direction), and the other is infinite along the tube axis. Thus, ΔP = γ/r. At the base of a capillary rise of height h, this pressure difference must equal the hydrostatic pressure of the liquid column: γ/r = ρ g h, which rearranges to Jurin's law when including the contact angle term. This connection reveals that capillary rise is fundamentally a pressure-balance phenomenon driven by interface curvature.

Industrial Applications Across Length Scales

In microfluidics and lab-on-a-chip devices, capillary forces dominate over gravity at the sub-millimeter scale. Passive capillary pumping eliminates the need for external pressure sources in diagnostic assays. A typical lateral flow immunoassay strip (e.g., rapid COVID-19 test) uses a porous nitrocellulose membrane with effective pore radii of 5-15 μm. Fluid wicks along the strip at 1-2 mm/s driven purely by capillarity. Design engineers tune pore size distributions and surface treatments to control flow rates, ensuring proper reagent mixing and reaction times. Non-ideal wetting (θ > 20°) can cause flow stalling, a common failure mode addressed by plasma treatment to reduce contact angles below 10°.

In inkjet printing, both thermal and piezoelectric printheads exploit surface tension to form droplets. After ejection, a liquid thread breaks up due to Rayleigh-Plateau instability — surface tension minimizes surface area by fragmenting the thread into droplets. Typical droplet diameters are 20-50 μm, with ejection velocities of 5-10 m/s. Surface tension must overcome viscous resistance and inertia within microseconds. Ink formulations balance γ (typically 0.025-0.035 N/m) with viscosity (2-10 mPa·s) to achieve reliable jetting without satellite droplet formation. Contact angle on paper substrates (~30-50°) governs dot spreading and print quality; overly hydrophilic papers cause feathering while hydrophobic coatings produce sharp features but poor adhesion.

In enhanced oil recovery (EOR), surfactants reduce interfacial tension between oil and water from ~30 mN/m to as low as 0.001 mN/m (ultra-low IFT). Capillary pressure traps residual oil in porous rock pores after primary production. The capillary number Ca = μv/γ quantifies viscous forces relative to capillary forces; at typical reservoir flow velocities (v ~ 1 m/day ≈ 10⁻⁵ m/s, μ ~ 1 mPa·s), Ca ~ 10⁻⁷, indicating capillary-dominated flow. Reducing γ by 1000× increases Ca to 10⁻⁴, mobilizing trapped oil ganglia. This surface tension reduction can recover an additional 10-20% of original oil in place, economically transforming marginal fields.

Worked Example: Capillary Rise in a Microfluidic Channel

Problem: A microfluidic diagnostic device uses rectangular glass channels coated with a hydrophilic layer to transport reagent solutions via capillary action. The channel has a hydraulic diameter of 180 μm, and the reagent is an aqueous buffer with density 1015 kg/m³ and surface tension 0.0685 N/m at 25°C. The contact angle of the buffer with the treated glass is 12°. Calculate: (a) the maximum capillary rise height, (b) the capillary pressure at the meniscus, (c) the time required to fill a 15 mm length of channel assuming Lucas-Washburn dynamics with viscosity 1.15 mPa·s, and (d) the required contact angle to achieve a 25 mm rise.

Solution:

(a) Maximum capillary rise: The hydraulic diameter D_h = 180 μm = 180 × 10⁻⁶ m serves as the effective diameter for circular capillary approximation. The hydraulic radius r_h = D_h/2 = 90 × 10⁻⁶ m. Contact angle θ = 12° = 0.2094 rad. Using Jurin's law:

h = (2γ cos θ) / (ρ g r_h)

h = [2(0.0685 N/m)(cos 12°)] / [(1015 kg/m³)(9.81 m/s²)(90 × 10⁻⁶ m)]

cos(12°) = 0.9781

Numerator: 2(0.0685)(0.9781) = 0.1340 N/m

Denominator: (1015)(9.81)(90 × 10⁻⁶) = 0.8961 kg/(m·s²) = 0.8961 N/m

h = 0.1340 / 0.8961 = 0.1496 m = 149.6 mm

This large theoretical rise indicates the channel can easily transport fluid vertically if needed, though real devices are typically horizontal.

(b) Capillary pressure at meniscus: For a cylindrical meniscus approximation:

ΔP = (2γ cos θ) / r_h

ΔP = [2(0.0685)(0.9781)] / (90 × 10⁻⁶)

ΔP = 0.1340 / (90 × 10⁻⁶) = 1489 Pa ≈ 1.49 kPa

This pressure drives flow along the channel. Alternatively, ΔP = ρ g h = (1015)(9.81)(0.1496) = 1489 Pa, confirming consistency.

(c) Filling time via Lucas-Washburn equation: The Lucas-Washburn equation for position as a function of time is:

h(t)² = (γ r_h cos θ / 2η) t

Solving for t when h = 15 mm = 0.015 m:

t = [2η h²] / (γ r_h cos θ)

η = 1.15 mPa·s = 1.15 × 10⁻³ Pa·s

t = [2(1.15 × 10⁻³)(0.015)²] / [(0.0685)(90 × 10⁻⁶)(0.9781)]

Numerator: 2(1.15 × 10⁻³)(2.25 × 10⁻⁴) = 5.175 × 10⁻⁷ Pa·s·m²

Denominator: (0.0685)(90 × 10⁻⁶)(0.9781) = 6.029 × 10⁻⁶ N·m/m² = 6.029 × 10⁻⁶ N

Wait, dimensional check: denominator should be in units that cancel with numerator to yield seconds. Correcting:

t = [2η h²] / (γ r_h cos θ) has units [Pa·s·m²] / [N/m·m·1] = [Pa·s·m²] / [N] = [N·s·m²/m²] / [N] = s ✓

t = 5.175 × 10⁻⁷ / 6.029 × 10⁻⁶ = 0.0858 s ≈ 86 ms

The channel fills rapidly in under 0.1 second, demonstrating the power of capillary action at this scale.

(d) Required contact angle for 25 mm rise: Rearranging Jurin's law to solve for θ:

cos θ = (ρ g r_h h) / (2γ)

h = 25 mm = 0.025 m

cos θ = [(1015)(9.81)(90 × 10⁻⁶)(0.025)] / [2(0.0685)]

Numerator: (1015)(9.81)(90 × 10⁻⁶)(0.025) = 0.02240 N/m

Denominator: 2(0.0685) = 0.1370 N/m

cos θ = 0.02240 / 0.1370 = 0.1635

θ = arccos(0.1635) = 80.59° ≈ 80.6°

Achieving only 25 mm rise would require significantly reducing surface wettability, approaching neutral wetting. This illustrates why microfluidic devices demand highly hydrophilic coatings (θ < 20°) to maintain strong capillary driving forces.

Edge Cases and Practical Limitations

Real systems deviate from idealized models in several critical ways. Contact angle hysteresis — the difference between advancing and receding contact angles on the same surface — typically ranges from 5° to 30° depending on surface roughness and chemical heterogeneity. During capillary rise, the advancing angle governs the initial ascent, but the receding angle determines whether the liquid retracts when external disturbances occur. Hysteresis can trap liquid columns at intermediate heights below the equilibrium prediction.

Gravity-viscosity competition: Jurin's law assumes equilibrium, but viscous resistance delays the approach to final height. For highly viscous fluids (η > 100 mPa·s), filling times can extend to minutes or hours in fine capillaries. Glycerin (η ≈ 1400 mPa·s at 20°C) in a 100 μm tube requires approximately 200 seconds to reach 90% of equilibrium height, compared to 0.2 seconds for water. This dramatically affects process cycle times in coating, impregnation, and wicking applications.

Evaporation effects: In open capillaries, evaporation from the meniscus competes with capillary rise, establishing a steady-state where evaporative mass loss equals capillary flow. For water in ambient air, evaporation rates of 10⁻⁷ kg/(m²·s) are typical. In a 100 μm tube, this corresponds to a descent rate of approximately 1.3 μm/s. Over hours, evaporation-induced flow can deposit solutes at the meniscus, creating "coffee ring" patterns exploited in deposition processes but detrimental to uniform coatings.

Dynamic contact angle effects: At high flow velocities, the contact angle deviates from its static value. The advancing contact angle increases with velocity, reducing capillary driving force. Empirical correlations like the Hoffman-Voinov-Tanner law quantify this: (θ_dynamic)³ ≈ (θ_static)³ + Ca, where Ca is the capillary number. For inkjet ejection at Ca ~ 0.1, contact angle increases can reach 10-20°, significantly impacting droplet impact dynamics and spreading on substrates.

Frequently Asked Questions