Skip to content

Modified Berggren Formula

Reviewed by licensed engineers

Checked against its sources by licensed professional engineers on the previous site, where reviewers weren’t recorded by name.

How articles are reviewed

The modified Berggren formula was developed in the early 1950s to address the shortcomings of the Stefan formula. The modified Berggren formula assumes that the soil is a semi-infinite mass with uniform properties and existing initially at a uniform temperature (TiT_{i}). It is further assumed that surface temperature is suddenly changed from TiT_{i} to TsT_{s} (below freezing). The modified Berggren formula is simply the Stefan formula corrected for the effects of temperature changes in the soil mass:

x=λ48kavg nFILx = \lambda \sqrt{\frac{48 k_{avg}\, n\mathit{FI}}{L}}
Table 1.
where:xequalsdepth of freeze or thaw, (ft))
λequalsdimensionless coefficient which takes into consideration the effect of temperature changes in the soil mass (i.e., a fudge factor). Corrects the Stefan formula for the neglected effects of volumetric heats (accounts for "sensible heat" changes)
k_(avg)equalsthermal conductivity of soil, average of frozen and unfrozen (BTU/hr • ft • °F)
nequalsconversion factor for air freezing (or thawing) index to surface freezing (or thawing) index
FIequalsair freezing index (°F • days)
TIequalsair thawing index (°F • days)
Lequalslatent heat (BTU/ft³)

Determination of λ

l can be determined by chart (Figure 1) based on inputs of a (thermal ratio) and m (fusion parameter).

Table 2.
λequalsf (FI (or TI), mean annual air or ground temperature, thermal properties of soil)
equalsf(m,a) and can be read from Figure 1
μequalsfusion parameter=(∣Tf−Ts∣)(CL)= \left( \left| T_f - T_s \right| \right) \left( \dfrac{C}{L} \right)
Cequalsaverage volumetric heat capacity of a soil (BTU/ft³ • °F)
Lequalslatent heat (BTU/ft³)
equalssurface freezing (or thawing) index, nFI (or nTI) divided by length of freezing (or thawing) season. Represents temperature differential between average surface temperature and 32 °F taken over the entire freeze (or thaw) season.
equalsnFId or nTId\dfrac{n\mathit{FI}}{d} \text{ or } \dfrac{n\mathit{TI}}{d}
dequalslength of freezing or thawing duration. For example, if the winter freezing season is December through February, then the duration of freezing (d) equals about 90 days.
T_fequals32 °F
T_sequalsaverage surface temperature for the freezing (or thawing) period
αequalsthermal ratio=∣Tˉ−TfTf−Ts∣= \left| \dfrac{\bar{T} - T_f}{T_f - T_s} \right|
Tequalsaverage annual air or ground temperature
|T - T_f|equalsrepresents the amount that the mean annual temperature exceeds (or is less than) the freezing point of the soil moisture (assumed to be 32 °F).
l coefficient in the Modified Berggren Formula.
Figure 1. l coefficient in the Modified Berggren Formula.
Pictorial representation of variables involved with l.
Figure 2. Pictorial representation of variables involved with l.

Modified Berggren Formula Example

Determine the depth of frost penetration into a homogeneous sandy silt for the following conditions:

Given:

  • Mean annual temperature = 48 °F
  • Surface freezing index = nFI = 750 °F • days
  • Duration of freezing season = d = 100 days
  • Soil properties:
  • Dry density = gdg_{d}= 100 lb/ft³
  • Water content = w = 15% ''

Solution

  • Calculate soil thermal properties

Volumetric latent heat of fusion: L = (144 BTU/lb) (100 lb/ft³) (15/100) = 2160 BTU/ft³

Average volumetric specific heat: CavgC_{avg} = 100 (lb/ft³) (0.15) = 28.2 BTU/ft³• °F

Average thermal conductivity: kfk_{f} @ 0.80 BTU/hr • ft • °F kuk_{u} @ 0.72 BTU/hr • ft • °F

Therefore, kavgk_{avg} = 0.76 BTU/hr • ft • °F

  • Calculate λ (recall λ = f(α, μ))

α=∣Tˉ−TfTf−Ts∣=∣Tˉ−Tf∣(d)nFI=∣48−32∣(100)750=2.13\alpha = \left| \frac{\bar{T} - T_f}{T_f - T_s} \right| = \frac{\left| \bar{T} - T_f \right| (d)}{n\mathit{FI}} = \frac{\left| 48 - 32 \right| (100)}{750} = 2.13

μ=(∣Tf−Ts∣)(CavgL)=nFId(CavgL)=(750100)(28.22160)=0.10\mu = \left( \left| T_f - T_s \right| \right) \left( \frac{C_{avg}}{L} \right) = \frac{n\mathit{FI}}{d} \left( \frac{C_{avg}}{L} \right) = \left( \frac{750}{100} \right) \left( \frac{28.2}{2160} \right) = 0.10

From the chart, λ= 0.74

  • Calculate depth of freezing

x=λ48kavg nFIL=(0.74)48 (0.76) (750)2160=2.7 ft. (0.8 m)x = \lambda \sqrt{\frac{48 k_{avg}\, n\mathrm{FI}}{\mathrm{L}}} = (0.74) \sqrt{\frac{48\,(0.76)\,(750)}{2160}} = 2.7\ \mathrm{ft.}\ (0.8\ \mathrm{m})

Calculator

Frost depth (Stefan and modified Berggren)

Stefan · modified Berggren (Aldrich and Paynter)
Climate
Soil

Average of frozen and unfrozen values (Kersten charts).

Read from α and μ, shown in the working. Leave blank for Stefan only.

Example values: the Modified Berggren Formula article’s sandy silt.

Teaching this? Figures and worked examples are free to use with credit. Course reading lists

Cite this page

APA

Pavement Interactive. (2018, September 19). Modified berggren formula. HeadLight. Retrieved [the date you read it], from https://pavementinteractive.org/reference-desk/design/design-parameters/modified-berggren-formula/

ASCE

Pavement Interactive. (2018). “Modified Berggren Formula.” HeadLight, <https://pavementinteractive.org/reference-desk/design/design-parameters/modified-berggren-formula/> ([the date you read it]).

BibTeX

@misc{pavint-modified-berggren-formula,
  author = {{Pavement Interactive}},
  title = {Modified Berggren Formula},
  year = {2018},
  publisher = {HeadLight},
  howpublished = {\url{https://pavementinteractive.org/reference-desk/design/design-parameters/modified-berggren-formula/}},
  note = {Accessed [date]}
}
This address will not change. The citation records the date you read it.
Report an error

Reports go to the editors, who check them against the cited sources. Corrections are noted at the top of the article.

Only if you want a reply. It is not published.