Buried Pipe Heat Loss by the Soil Shape Factor: Conduction Spreading, Insulation Diameter, and Moisture Sensitivity

A buried pipe loses heat by conduction through the soil in curved spreading paths to the ground surface, and the geometry of that spreading is captured by a conduction shape factor, not a convection coefficient. Getting buried pipe heat loss right requires understanding why the soil resistance takes the form acosh(2z/D)/(2πk), which surface the soil actually contacts on an insulated line, and why soil moisture is the dominant uncertainty in nearly every buried-pipe thermal calculation.

Why a Buried Pipe Is Not an Above-Ground Insulation Problem

An above-ground insulated pipe sheds heat through its insulation and then across an outer surface film to the air. The outer resistance is a convection coefficient, driven by air velocity and pipe geometry. Bury the same pipe and the air film disappears entirely — the soil itself becomes the outer resistance. Heat leaves the buried surface and fans outward and upward through the ground to the surface, following curved conduction paths that expand as they go.

The resistance of that spreading path is the conduction shape factor for a horizontal cylinder buried in a semi-infinite medium, per Incropera, Bergman, Lavine and DeWitt (Fundamentals of Heat and Mass Transfer, Wiley) and Carslaw and Jaeger (Conduction of Heat in Solids, Oxford):

R_soil = acosh(2z / D_boundary) / (2π × k_soil)    [m·K/W per unit length]

A deeper pipe has more soil around it; the resistance is higher. A larger buried diameter gives the soil a shorter effective path to the surface; the resistance is lower. The shape factor bundles the three-dimensional spreading geometry into one expression.

The calculation divides the pipe-to-ground temperature difference by a chain of series resistances: soil, insulation if any, and the pipe wall for plastic. For a bare pipe the soil is the entire resistance. Two details govern every result: the soil sees the outside of the insulation on an insulated line, not the pipe itself, and soil conductivity swings nearly ten to one with moisture content.

This closes the heat-transfer pair with the Insulated Pipe Heat Loss Calculator and its companion article: that analysis uses an outer air film above ground; this one replaces the air film with soil shape-factor conduction below grade. Same series-resistance logic, different outer boundary.

Calculator Inputs: Installation, Pipe, Depth, Soil Moisture, Insulation

The calculator works through a linear sequence of inputs.

Installation type selects bare or insulated. Insulation fields appear only for insulated lines.

Unit system sets US (°F, in, ft, Btu/h·ft) or Metric (°C, mm, m, W/m). All outputs report in the chosen system.

Pipe and fluid temperature is the outer-surface temperature of the pipe. For bare steel or copper this is effectively the fluid temperature; for plastic or insulated pipe the model treats the entered value as the pipe outer-surface temperature.

Ground boundary temperature is the design ground-surface or design boundary value from the project basis. The calculator does not compute seasonal ground-temperature swings; use the design condition for the analysis.

Pipe size can be entered as a nominal size (ASME B36.10M OD lookup, NPS ½ to 24) or as a direct outer diameter. The B36.10M lookup covers Schedule 40 steel and standard copper tubing.

Burial depth has three input options: depth to pipe centre (direct z), cover over the pipe crown (z = cover + pipe radius), or cover over the insulation crown (z = cover + insulation radius). The shape factor formula needs the centre depth; the calculator converts from whichever basis the drawing uses.

Soil type and moisture sets the soil thermal conductivity from three defaults: dry (k ≈ 0.3 W/m·K, 0.17 Btu/h·ft·°F), moist (k ≈ 1.0 W/m·K, 0.58), and wet/saturated (k ≈ 2.5 W/m·K, 1.44), per ASHRAE Handbook Fundamentals and IGSHPA ground-source references. A conductivity override field accepts a geotechnical-measured value.

Insulation inputs (insulated lines) include material (polyurethane foam k ≈ 0.025 W/m·K; mineral/glass wool 0.040; cellular glass 0.050, per manufacturer data from Logstor, Perma-Pipe, and Rovanco), insulation thickness (which sets insulation OD = pipe OD + 2 × thickness), and an optional jacket OD for pre-insulated pipe systems per EN 13941.

Pipe wall resistance is optional. For steel and copper it is negligible. For plastic (PVC, HDPE, PP) or thick-wall pipe it matters; the calculator adds the annular wall resistance when pipe ID and wall conductivity are entered.

Outputs include heat loss per unit length (primary), total heat for an optional entered run length, the resistance breakdown (soil/insulation/pipe wall shares), the governing resistance, bare-pipe baseline and percent reduction on insulated lines, and a pass/fail against an optional allowable heat loss limit.

The Soil Shape Factor: Conduction Spreading to the Surface

The soil resistance comes from a conduction shape factor for a horizontal cylinder buried in a semi-infinite medium. Per Incropera and Carslaw and Jaeger, that shape factor converts the three-dimensional spreading of heat through the ground into a single resistance:

R_soil = acosh(2z / D_boundary) / (2π × k_soil)    [m·K/W per unit length]

  z          = depth to pipe CENTRE [m]            (typical range: 0.3 – 3.0 m)
  D_boundary = diameter of the surface touching soil [m]
               (pipe OD for bare; insulation OD for insulated)
  k_soil     = soil thermal conductivity [W/m·K]   (0.3 – 2.5)
  acosh      = inverse hyperbolic cosine

The equivalent shape-factor form cross-checks the resistance form:

S = 2πL / acosh(2z / D_boundary)     [m, conduction shape factor for length L]
Q = S × k_soil × ΔT                  [W, total heat for length L]

Heat leaving a buried cylinder does not flow through a flat slab or a simple annulus; it spreads outward and curves up to the surface through an expanding volume of soil. The shape factor is the analytical solution to that geometry, per Incropera and Carslaw and Jaeger, packaging the three-dimensional spreading into a single resistance expression.

The ratio 2z/D sets the resistance. A pipe buried deep (large z/D) has a longer, more resistant soil path; shallow burial (small z/D) has less. Doubling depth raises resistance through acosh, not linearly.

Worked for the bare NPS 4 case (D_boundary 0.1 m = 3.94 in, z 0.5 m = 1.64 ft, k 0.9 W/m·K = 0.520 Btu/h·ft·°F):

2z/D = 2 × 0.5/0.1 = 10     (z/D = 5; deep burial, ln form valid — see next section)
R_soil = ln(4 × 0.5/0.1) / (2π × 0.9) = ln(20) / 5.655 = 2.996 / 5.655
       = 0.530 m·K/W  (0.309 h·ft·°F/Btu)

Using a flat-wall or convection formula for the soil path is wrong; the geometry is the buried-cylinder shape factor, not a slab.

acosh versus ln: The Deep-Burial Simplification and Its Limit

The exact shape factor uses the inverse hyperbolic cosine. For pipes buried deep relative to their diameter, it simplifies to a natural logarithm, and knowing when that simplification holds keeps shallow-burial results accurate.

Exact:              R_soil = acosh(2z/D) / (2πk)
Deep burial only:   R_soil ≈ ln(4z/D)   / (2πk)    [valid for z/D > 1.5]

When the pipe is deep enough (depth-to-centre more than about 1.5 times the buried diameter), acosh(2z/D) ≈ ln(4z/D). This log form is the one most textbooks reproduce. At shallow depths the pipe's proximity to the surface makes the geometry more sensitive to depth, and the log approximation drifts. Per Incropera: use the full acosh for shallow burial.

At z/D = 5 (deep):

acosh(10) = ln(10 + √(100−1)) = ln(19.950) = 2.993
ln(4 × 5) = ln(20)            = 2.996
Difference: 0.1%  — negligible; log form applies

At z/D = 0.75 (shallow):

acosh(1.5) = 0.9624
ln(3)      = 1.0986
Difference: 14%  — significant; full acosh required

The calculator uses the exact acosh throughout. Most buried distribution mains are at z/D ≥ 3 and the log form holds, but shallow buried pipe (thin cover) and shallow conduit at z/D < 1.5 require the full expression.

Resistances in Series: Soil, Insulation, and the Pipe Wall

A buried pipe's total thermal resistance is the soil path, any insulation layer, and the pipe wall in series. Heat loss per unit length is:

q = (T_pipe − T_ground) / R_total     [W/m]   positive = loss; negative = gain

R_total = R_soil + R_ins + R_pipe

R_soil = acosh(2z/D_boundary) / (2π k_soil)
R_ins  = ln(D_ins_OD / D_pipe_OD)   / (2π k_ins)    [0 if bare]
R_pipe = ln(D_pipe_OD / D_pipe_ID)  / (2π k_pipe)   [optional, plastic/thick-wall]

For a bare pipe the soil is the only resistance: R_total = R_soil, and heat loss scales almost directly with soil conductivity. Adding insulation introduces R_ins, which typically becomes 80 to 95 percent of R_total once the insulation thickness is substantial, reducing sensitivity to soil moisture (see Section 8). The soil resistance never disappears; it is always in the circuit.

The calculator identifies the governing resistance. Bare lines are soil-governed; insulated lines are typically insulation-governed. That label tells you the design lever: on a bare line, soil moisture dominates; on a heavily insulated line, insulation thickness and conductivity dominate.

For steel or copper, R_pipe = ln(OD/ID)/(2π × k_metal) is negligible (steel k ≈ 50 W/m·K, copper k ≈ 400 W/m·K). For HDPE (k ≈ 0.4 W/m·K) or PVC (k ≈ 0.16 W/m·K) on thick-wall pipe, the pipe-wall term is worth including.

The same series logic applies in the Insulated Pipe Heat Loss article (above-ground composite cylinder), except that article uses an outer air-film convection resistance where this one uses R_soil. Same framework, different outer boundary.

The Insulation Diameter Rule: What Surface the Soil Actually Sees

On an insulated buried pipe, the soil contacts the outside of the insulation jacket, not the pipe. The soil shape factor must therefore use the insulation outer diameter as D_boundary, not the pipe OD. Per Incropera: using the pipe diameter is the most common buried-pipe thermal error.

Bare pipe:      D_boundary = pipe OD
Insulated pipe: D_boundary = insulation OD   ← soil contacts this surface

The error direction: using the pipe OD (smaller value) instead of the insulation OD (larger) for a given depth gives a larger 2z/D ratio, a larger acosh, and an overstated R_soil. The calculation understates soil-path heat loss on an insulated line.

For the insulated example (D_pipe 0.1 m = 3.94 in, insulation 50 mm = 1.97 in foam on each side):

Insulation OD = 0.1 + 2 × 0.05 = 0.2 m (7.87 in)
D_boundary = 0.2 m (insulated) vs 0.1 m (bare)

Bare:      R_soil = ln(4×0.5/0.1)/(2π×0.9) = ln(20)/5.655 = 0.530 m·K/W
Insulated: R_soil = ln(4×0.5/0.2)/(2π×0.9) = ln(10)/5.655 = 2.303/5.655 = 0.407 m·K/W

The soil resistance drops from 0.530 to 0.407 m·K/W (0.309 to 0.237 h·ft·°F/Btu) when insulation enlarges D_boundary: the bigger buried surface has a shorter effective soil path to the surface. Using the pipe's 0.1 m would wrongly hold R_soil at 0.530.

The same rule governs the depth basis: when cover is measured to the insulation crown, z = cover + insulation radius, not pipe radius (see next section).

Soil Moisture: The Eight-to-One Conductivity Range

Soil thermal conductivity is the input that most often controls a buried-pipe result, and it spans nearly ten to one from dry to saturated soil, driven almost entirely by moisture content.

Dry soil:        k ≈ 0.3 W/m·K  (0.17 Btu/h·ft·°F)
Moist soil:      k ≈ 1.0 W/m·K  (0.58)
Wet/saturated:   k ≈ 2.5 W/m·K  (1.44)

Per ASHRAE Handbook Fundamentals and IGSHPA: water conducts heat far better than the air it displaces in soil pores. Wetting replaces insulating air with conducting water, raising k sharply. The dry-to-saturated range spans approximately 8.3 to one (2.5/0.3).

For a bare pipe, R_total = R_soil, which is proportional to 1/k_soil. Heat loss is therefore proportional to k_soil directly:

Example 3 — bare 100 mm pipe, same geometry, 70°C (126°F) difference:
Dry   (k 0.3):  R = ln(20)/(2π×0.3) = 2.996/1.885 = 1.590 m·K/W;  q = 70/1.590 =  44 W/m  (45.8 Btu/h·ft)
Moist (k 0.9):  R = 0.530 m·K/W;                                     q = 70/0.530 = 132 W/m  (137 Btu/h·ft)
Wet   (k 2.5):  R = ln(20)/(2π×2.5) = 2.996/15.708 = 0.191 m·K/W;  q = 70/0.191 = 367 W/m  (381 Btu/h·ft)
Dry-to-wet: 44 → 367 W/m, 8.3× range (equals k_wet/k_dry = 2.5/0.3)

The wet pipe loses 8.3 times the dry because soil is the only resistance on a bare line. When moisture is uncertain, bracket the design with dry and wet runs. Adding insulation reduces this sensitivity substantially: once insulation governs 80 to 95 percent of R_total, the soil contribution is a minor fraction and the moisture swing matters far less.

Use the soil moisture condition for the site and season; geotechnical data is the most reliable source. The calculator shows the conductivity value with its moisture label.

Depth to the Pipe Centre: Converting Cover to Centreline

The shape factor formula uses depth to the pipe centreline. Drawings typically show cover to the top of the pipe or insulation. Three input options convert automatically:

Basis 1 (centre directly):           z = entered value
Basis 2 (cover to pipe crown):       z = cover + pipe radius        (bare pipe)
Basis 3 (cover to insulation crown): z = cover + insulation radius  (insulated)

For cover to insulation crown, add the insulation radius, not the pipe radius. The formula uses depth to the centre of the buried circle, and for an insulated line that circle is the insulation jacket. Adding the pipe radius in this case places the centreline too high and understates z.

Example (insulated, insulation OD 200 mm = 7.87 in, radius 100 mm = 3.94 in, cover over insulation crown 450 mm = 17.7 in):

z = 0.450 + 0.100 = 0.550 m (1.80 ft) to pipe centre
Using 0.450 m directly: z too shallow by one insulation radius
Adding pipe radius (50 mm) instead: z = 0.500 m, still 50 mm low

On shallow pipe (small z), one missed radius shifts z/D noticeably and changes R_soil. For deeply buried mains at 1.5 m or more, the error is a small fraction of z and affects R_soil modestly. The calculator converts from whichever cover basis matches the project drawing.

Bare Pipe Worked Example: 100 mm Main at 132 Watts per Metre

Scenario: Bare buried steel main, 100 mm OD (NPS 4, 3.94 in), centre 0.5 m (1.64 ft) deep, moist soil k 0.9 W/m·K (0.520 Btu/h·ft·°F), pipe surface 80°C (176°F), ground 10°C (50°F), 30 m (98.4 ft) run.

Step 1. Depth ratio and formula form:

D_boundary = pipe OD = 0.1 m (bare)
z = 0.5 m;   z/D = 0.5/0.1 = 5 → deep (z/D > 1.5), ln form applies

Step 2. Soil resistance:

R_soil = ln(4z/D) / (2π k_soil) = ln(4×0.5/0.1) / (2π×0.9)
       = ln(20) / 5.655 = 2.996 / 5.655 = 0.530 m·K/W  (0.309 h·ft·°F/Btu)

Step 3. Total resistance:

R_total = R_soil = 0.530 m·K/W   (bare, soil governs 100%)

Step 4. Heat loss per metre:

q = (80 − 10) / 0.530 = 70 / 0.530 = 132 W/m  (137 Btu/h·ft)

Step 5. Total over run:

Q = 132 × 30 = 3,960 W  (13,510 Btu/h)

Step 6. Moisture sensitivity (Example 3):

Dry   (k 0.3):   q =  44 W/m  (45.8 Btu/h·ft)
Moist (k 0.9):   q = 132 W/m  (137 Btu/h·ft)  ← base case
Wet   (k 2.5):   q = 367 W/m  (381 Btu/h·ft)
Dry-to-wet range: 8.3×

Step 7. Why the eight-to-one range: soil is the only resistance, so heat loss tracks k_soil directly. The 8.3:1 moisture range in k produces an 8.3:1 range in heat loss.

Step 8. Depth effect: deeper burial raises R_soil through acosh/ln, reducing loss, but the dependence is logarithmic. Doubling z from 0.5 to 1.0 m changes R_soil from 0.530 to ln(40)/5.655 = 3.689/5.655 = 0.652 m·K/W, a 23% increase in resistance and 19% reduction in heat loss, not a 50% reduction. Diminishing returns apply; trench depth is governed by frost depth and clearances, not thermal optimisation.

Step 9. Result: bare 100 mm main, moist soil, 70°C (126°F) difference: 132 W/m (137 Btu/h·ft), 3,960 W (13,510 Btu/h) for 30 m. Moisture spans 44 to 367 W/m. Soil governs entirely. Add insulation to cut loss and reduce moisture sensitivity.

Insulated and Chilled Worked Examples: 89 Percent Reduction and Heat Gain

Insulated example (Example 2): Same 100 mm pipe, 50 mm polyurethane foam (k 0.025 W/m·K, 0.0145 Btu/h·ft·°F), same conditions (80°C pipe, 10°C ground, k_soil 0.9 W/m·K).

Step 1. Insulation OD becomes the new D_boundary:

Insulation OD = 0.1 + 2×0.05 = 0.2 m  (7.87 in)   ← D_boundary for soil

Step 2. Insulation resistance (annular logarithm):

R_ins = ln(0.2/0.1) / (2π×0.025) = ln(2) / 0.1571 = 0.693 / 0.1571 = 4.415 m·K/W  (2.571 h·ft·°F/Btu)

Step 3. Soil resistance (using new D_boundary = 0.2 m):

R_soil = ln(4×0.5/0.2) / (2π×0.9) = ln(10) / 5.655 = 2.303 / 5.655 = 0.407 m·K/W  (0.237 h·ft·°F/Btu)

Step 4. Total and breakdown:

R_total = 4.415 + 0.407 = 4.822 m·K/W  (2.808 h·ft·°F/Btu)
Insulation share: 4.415/4.822 = 91.6% ≈ 92%
Soil share: 8%

Step 5. Heat loss:

q = 70 / 4.822 = 14.5 W/m  (15.1 Btu/h·ft)
Q = 14.5 × 30 = 435 W  (1,484 Btu/h)

Step 6. Reduction from bare:

Bare: 132 W/m → Insulated: 14.5 W/m
Reduction: (132 − 14.5) / 132 = 89%

Step 7. Why R_soil dropped: R_soil fell from 0.530 m·K/W (bare, D = 0.1 m) to 0.407 m·K/W (insulated, D = 0.2 m) because the larger buried insulation surface has a shorter effective soil path to the surface. The insulation-diameter rule from Section 7 directly reduces R_soil.

Chilled example (Example 4): Same insulated pipe, chilled water 5°C (41°F), ground 25°C (77°F).

Step 8. Temperature difference and direction:

ΔT = 5 − 25 = −20°C (−36°F)  → heat GAIN (ground warmer than pipe)
q = 20 / 4.822 = 4.1 W/m  (4.3 Btu/h·ft) heat gain

Step 9. Same resistances, reversed direction: the insulated resistances are identical; only the sign of ΔT changes. The ground delivers heat to the chilled pipe. Condensation on the outer jacket and vapour barrier adequacy are separate checks outside this calculator.

Step 10. Results: insulated hot line 14.5 W/m (15.1 Btu/h·ft), 89% reduction from bare, insulation governs at 92%. Chilled line 4.1 W/m (4.3 Btu/h·ft) heat gain, same resistance chain reversed. For chilled lines in warm saturated soil, heat gain can be considerably higher; bracket using the wet-soil conductivity.

Steady-State Heat Loss Is Not a Freeze Prediction

Steady-state buried-pipe heat loss and freeze protection address related but distinct questions. This calculator gives the rate the pipe sheds heat at fixed boundary temperatures. Whether the pipe freezes and when is a transient question — the steady-state loss is one input, not the answer.

Steady-state heat loss: rate of heat exchange at fixed T_pipe and T_ground [W/m]
Freeze prediction: does water reach 0°C (32°F), and when? [transient — hours to days]

A high-loss pipe may never freeze if water flows continuously, because incoming flow continuously replenishes thermal energy. A modest-loss pipe on a stagnant dead leg may freeze in hours. Freezing depends on heat loss AND fluid heat capacity, flow rate, inlet temperature, and stagnation time; the steady-state loss is one input, not the answer.

What freeze protection adds: trace heating sized to replace the per-metre steady-state loss, flow maintenance to keep warm water circulating, and insulation to reduce the rate. A transient time-to-freeze analysis then confirms whether these measures hold for the design stagnation period.

Per heat-transfer practice and ASHRAE district heating guidance: use the steady-state heat loss as the input to freeze-protection sizing and trace-heat specification, not as a freeze guarantee. Freezing is a transient energy balance involving heat capacity, flow, and time.

Application Boundaries: Adjacent Pipes, Groundwater, Layered Soil, Casings

The calculator models a single isolated horizontal pipe in uniform semi-infinite soil at steady state. Several scenarios fall outside that scope:

Adjacent pipes. One pipe only. A supply-and-return pair or shared-trench bundle interact thermally; each pipe warms the soil around the other, reducing the effective temperature gradient for both. Multi-pipe buried analysis by superposition or numerical method is needed for close spacing.

Groundwater flow. Assumes still soil conduction. Flowing groundwater carries heat away by advection, raising effective heat loss above the conduction-only result. Significant in high-water-table or flowing-aquifer sites.

Layered soil and surface cover. Treats soil as uniform and semi-infinite. Pavement, snow cover, surface insulation boards, and soil layers with differing conductivities change both the effective boundary temperature and the resistance path. Numerical methods or adjusted conductivity/boundary are needed for layered profiles.

Seasonal and transient conditions. Steady-state at a fixed ground temperature. Seasonal ground-temperature swings (significant for shallow burial), startup transients, and freeze/thaw cycles require transient analysis.

Duct, conduit, casing, or tunnel. Assumes the buried surface contacts soil directly. A pipe inside a casing or conduit has an annular air space that adds resistance not captured here; a pipe in a utility tunnel or valve pit is an entirely different geometry.

Fluid temperature drop along run. Assumes uniform pipe temperature. Long hot-water distribution runs lose heat continuously, lowering fluid temperature along the way. Segment the run or apply an energy balance for accuracy on runs longer than a few hundred metres.

Condensation on chilled lines. Gives heat gain, not condensation risk or vapour-barrier adequacy. A chilled buried line with an inadequate vapour barrier can allow moisture intrusion into the insulation, degrading it toward ambient soil conductivity.

Per Incropera and ASHRAE: single isolated direct-buried pipe, uniform soil, steady-state conduction is the appropriate model for preliminary design and energy accounting. Adjacent pipes, groundwater, layered soil, surface cover, seasonal transients, casings and tunnels, fluid temperature drop, freeze prediction, and condensation require separate qualified analysis and a project-specific design basis.

Buried Pipe Heat Loss Calculator

Buried pipe heat loss by the soil conduction shape factor: computes the soil resistance from the buried-cylinder shape factor acosh(2z/D)/2πk, adds insulation and pipe-wall resistances in series, and divides the pipe-to-ground temperature difference by the total. The soil sees the insulation outer diameter on an insulated line, and soil conductivity is set by moisture condition (dry to saturated spans nearly eight to one). Reports heat loss per unit length, the resistance breakdown, and the governing resistance. Bare and insulated, hot and chilled, SI and US units.

FAQ

Related Calculators

• Insulated Pipe Heat Loss Calculator: The above-ground counterpart, using an outer air film instead of soil shape-factor conduction for insulated pipe in air (article).
• Condensate Return Line Sizing Calculator: Flash-steam condensate sizing for buried district-heating steam mains (article).
• Pipe Expansion Loop Sizing Calculator: Guided-cantilever thermal loops for buried and above-ground hot piping (article).
• Hazen-Williams Pipe Flow Calculator: Pressurized water pipe flow per AWWA M22 for buried water mains (article).
• Expansion Tank Sizing Calculator: Thermal expansion control per IPC Section 607.3 (article).
• Pump Power Calculator: Hydraulic, brake, and motor power for buried-main distribution pumps (article).
• Hanger Load Calculator: Support load per MSS SP-58 for above-ground sections of the run (article).