Buckingham π theorem

The Buckingham π theorem is a key theorem in dimensional analysis. The theorem loosely states that if we have a physically meaningful equation involving a certain number, "n", of physical variables, and these variables are expressible in terms of "k" independent fundamental physical quantities, then the original expression is equivalent to an equation involving a set of "p" = "n" − "k" dimensionless variables constructed from the original variables. This provides a method for computing sets of dimensionless parameters from the given variables, even if the form of the equation is still unknown. However, the choice of dimensionless parameters is not unique: Buckingham's theorem only provides a way of generating sets of dimensionless parameters, and will not choose the most 'physically meaningful'.

Statement

More formally, the number of dimensionless terms that can be formed, "p", is equal to the nullity of the dimensional matrix, and "k" is the rank. For the purposes of the experimenter, different systems which share the same description in terms of these dimensionless numbers are equivalent.

In mathematical terms, if we have a physically meaningful equation such as

:$f(q\_1,q\_2,ldots,q\_n)=0,!$

where the "q_i" are the "n" physical variables, and they are expressed in terms of "k" independent physical units, then the above equation can be restated as

:$F(pi\_1,pi\_2,ldots,pi\_p)=0,!$

where the π_i are dimensionless parameters constructed from the "q_i" by "p" = "n" − "k" equations of the form

:$pi\_i=q\_1^\{m\_1\},q\_2^\{m\_2\}cdots\; q\_n^\{m\_n\}$

where the exponents "m_i" are rational numbers (they can always be taken to be integers: just raise it to a power to clear denominators).

The use of the π_"i" as the dimensionless parameters was introduced by Edgar Buckingham in his original 1914 paper on the subject from which the theorem draws its name.

Significance

The Buckingham π theorem provides a method for computing sets of dimensionless parameters from the given variables, even if the form of the equation is still unknown. However, the choice of dimensionless parameters is not unique: Buckingham's theorem only provides a way of generating sets of dimensionless parameters, and will not choose the most 'physically meaningful'.

Two systems for which these parameters coincide are called "similar" (as with similar triangles, they differ only in scale); they are equivalent for the purposes of the equation, and the experimentalist who wants to determine the form of the equation can choose the most convenient one.

Proof

Outline

Begin by considering the space of fundamental and derived physical units as a vector space over the rational numbers, with the fundamental units as basis vectors, and with multiplication of physical units as the "vector addition" operation, and raising to powers as the "scalar multiplication" operation:represent a dimensional variable as the set of exponents needed for the fundamental units (with a power of zero if the particular fundamental unit is not present). For instance, the gravitational constant "g" has units of $ell/t^2=ell^1t^\{-2\}$ (distance over time squared), so it is represented as the vector $(1,-2)$ with respect to the basis of fundamental units (distance,time).

Making the physical units match across sets of physical equations can then be regarded as imposing linear constraints in the physical unit vector space.

Formal proof

(The examples make this clearer.)

Given a system of "n" dimensional variables (physical variables), in "k" (physical) dimensions, write the "dimensional matrix" "M", whose rows are the dimensions and whose columns are the variables: the "(i,j)"th entry is the power of the "i"th unit in the "j"th variable. The matrix can be interpreted as taking in a combination of the dimensional quantities and giving out the dimensions of this product. So:is the units of:$q\_1^\{m\_1\},q\_2^\{m\_2\}cdots\; q\_n^\{m\_n\}$

A dimensionless variable is a combination whose units are all zero (hence, dimensionless), which is equivalent to the kernel of this matrix; a dimensionless variable is a linear relation between units of dimensional variables.

By the rank-nullity theorem, a system of "n" vectors in "k" dimensions (where all dimensions are necessary) satisfies a "(p=n-k)"-dimensional space of relations. Any choice of basis will have "p" elements, which are the dimensionless variables.

The dimensionless variables can always be taken to be integer combinations of the dimensional variables (by clearing denominators). There is mathematically no natural choice of dimensionless variables; some choices of dimensionless variables are more physically meaningful, and these are what are ideally used.

Examples

Speed

This example is elementary, but demonstrates the general procedure.

Suppose a car is driving at 100 km/hour; how long does it take it to go 200 km?

This question has 2 fundamental physical units: time "t" and length "l", and 3 dimensional variables: distance "D", time taken "T", and velocity "V". Thus there is 3-2=1 dimensionless quantity.

The units of the dimensional quantities are

:$D\; =\; ell,\; T\; =\; t,\; V\; =\; ell/t.$

The dimensional matrix is:

:

(The rows correspond to the dimensions $ell$, and "t", and the columns to the dimensional variables "D", "T", "V". For instance, the 3rd column, (1, −1), states that the "V" (velocity) variable has units of $ell^1\; t^\{-1\}\; =\; ell/t$.)

This is in reduced row echelon form, so one can read off that the
kernel is generated by

:

(Were it not already reduced, one could perform Gauss-Jordan elimination on the dimensional matrix.)

Thus:(or some power thereof).

In dimensions:

:$pi=(ell)^\{-1\}(t)^1(ell/t)^1\; =\; 1$

is dimensionless.

Admittedly, here the relation is simply $D=VT$ (or rather $D-VT=0$), so $VT/D=1$ is dimensionless, and the dimensionless equation ($f(pi)=0$) is

:$VT/D\; -\; 1\; =\; 0$

which can be solved for time to

:$T\; =\; frac\{D\}\{V\}.$

However, the above dimensional analysis does not require any physical understanding, and is useful in less familiar situations.

The simple pendulum

We wish to determine the period "T" of small oscillations in a simple pendulum. It will be assumed that it is a function of the length "L" , the mass "M" , and the acceleration due to gravity on the surface of the Earth "g", which has units of length divided by time squared. The model is of the form

:$f(T,M,L,g)\; =\; 0.,$

(Note that it is written as a relation, not as a function: "T" isn't here written as a function of "M", "L", and "g".)

There are 3 fundamental physical units in this equation: time "t", mass "m", and length "l", and 4 dimensional variables, "T", "M", "L", and "g". Thus we need only 4−3=1 dimensionless parameter, denoted π, and the model can be re-expressed as

:$f(pi)\; =\; 0$

where π is given by

:$pi\; =T^\{m\_1\}M^\{m\_2\}L^\{m\_3\}g^\{m\_4\}$for some values of "m"₁, ..., "m"₄.

The units of the dimensional quantities are:

:$T\; =\; t,\; M\; =\; m,\; L\; =\; ell,\; g\; =\; ell/t^2.$

The dimensional matrix is:

:

(The rows correspond to the dimensions "t", "m", and "l", and the columns to the dimensional variables "T", "M", "L" and "g". For instance, the 4th column, (−2, 0, 1), states that the "g" variable has units of $t^\{-2\}m^0\; ell^1$.)

This is in reduced row echelon form, so one can read off that the
kernel is generated by

:

(Were it not already reduced, one could perform Gauss-Jordan elimination on the dimensional matrix.)

Thus:(or some power thereof).

In dimensions:

:$pi=(t)^2(m)^0(ell)^\{-1\}(ell/t^2)^1\; =\; 1$

is dimensionless.

This example is easy because 3 of the dimensional quantities are fundamental units, so the last ("g") is a combination of the previous.

Note that if m₂ were non-zero there would be no way to cancel the M value—therefore m₂ "must" be zero.Dimensional analysis has allowed us to conclude that the period of the pendulum is not a function of its mass.

The model can now be expressed as

:$f(gT^2/L)\; =\; 0.$

Solving

The above is how far dimensional analysis takes us; solving the equation requires further analysis or experimentation.

Assuming the zeroes of "f" are discrete, we can say"gT"²/"L" = "K"_"n" where "K_n" is the "n"th zero. If there is only one zero, then "gT"²/"L" = "K" . It requires more physical insight or an experiment to show that there is indeed only one zero and that the constant is in fact given by "K" = 4π² .

Generalizations

For large oscillations of a pendulum, the analysis is complicated by an additional dimensionless parameter, the maximum swing angle. The above analysis is a good approximation in the limit that this angle is zero.

The atomic bomb

In 1941, Sir Geoffrey I. Taylor used dimensional analysis to estimate the energy released in an atomic bomb explosion (Taylor, 1950a,b). The first atomic bomb was detonated near Alamogordo, New Mexico on July 16, 1945. In 1947, movies of the explosion were declassified, allowingSir Geoffrey to complete the analysis and estimate the energy released in the explosion, even though the energy release was still classified. The actual energy released was later declassified and its value was remarkably close to Taylor's estimate.

Taylor supposed that the description of the process was adequately described by fivephysical quantities: the time "t" since the detonation, the energy "E" which is released at a single point in space at detonation, the radius "R" of the shock wave at time "t" , the atmospheric pressure "p" and the ambient density ρ. There are only three fundamental physical units in this equation: mass, time, and length. Thus we need only 5 − 3 = 2 dimensionless parameters, which can be found to be

:$pi\_0=R,left(frac\{\; ho\}\{Et^2\}\; ight)^\{1/5\}$

and

:$pi\_1=p,left(frac\{t^6\}\{E^2\; ho^3\}\; ight)^\{1/5\}.$

The process can now be described by an equation of the form

:$f(pi\_0,pi\_1)=0,,$

or, equivalently

:$R=left(frac\{Et^2\}\{\; ho\}\; ight)^\{1/5\}g(pi\_1),$

where "g"(π₁) is some function of π₁. The energy in the explosion is expected to be huge, so that for times of the order of a second after the explosion, we can estimate π₁ to be approximately zero, and experiments using light explosives can be conducted to determine that "g"(0) is on the order of unity so that

:$Rapproxleft(frac\{Et^2\}\{\; ho\}\; ight)^\{1/5\}.$

This is Taylor's equation which, once he knew the radius of the explosion as a function of the time, allowed him to calculate the energy of the explosion. (Wan, 1989)

ee also

* Dimensional analysis
* Dimensionless number
* Similitude (model)

References

Exposition

* cite web
author=Hanche-Olsen, Harald | year=2004
url=http://www.math.ntnu.no/~hanche/notes/buckingham/buckingham-a4.pdf
title=Buckingham's pi-theorem | format=PDF
publisher=NTNU | accessmonthday=April 9 | accessyear=2007

* cite book
last = Hart
first = George W.
date = March 1, 1995
title = Multidimensional Analysis: Algebras and Systems for Science and Engineering
publisher = Springer-Verlag
id = ISBN 0387944176
url = http://www.georgehart.com/research/multanal.html

* cite book
last = Kline
first = Stephen J.
year = 1986
title = Similitude and Approximation Theory
publisher = Springer-Verlag, New York
id = ISBN 0387165185

* cite book
last = Wan
first = Frederic Y.M.
year = 1989
title = Mathematical Models and their Analysis
publisher = Harper & Row Publishers, New York
id = ISBN 0060469021

* cite web
author=Vignaux, G.A. | year=1991
url=http://www.mcs.vuw.ac.nz/~vignaux/docs/maxent.pdf
title=Dimensional analysis in data modelling | format=PDF
publisher= Victoria University of Wellington | accessmonthday=December 15 | accessyear=2005

* Mike Sheppard, 2007 [http://www.msu.edu/~sheppa28/constants/constants.html Systematic Search for Expressions of Dimensionless Constants using the NIST database of Physical Constants]

Original sources

* cite journal
last=Buckingham
first=E.
authorlink=Edgar Buckingham
title=On physically similar systems; illustrations of the use of dimensional equations
journal=Phys. Rev.
year=1914
volume= 4
pages=345–376
url=http://prola.aps.org/abstract/PR/v4/i4/p345_1
doi=10.1103/PhysRev.4.345

* cite journal
last=Buckingham
first=E.
authorlink=Edgar Buckingham
year=1915
title=The principle of similitude
journal=Nature
volume=96
pages=396–397
doi=10.1038/096396d0

* cite journal
last=Buckingham
first=E.
authorlink=Edgar Buckingham
year=1915
title=Model experiments and the forms of empirical equations
journal=Trans. A.S.M.E
volume=37
pages=263–296

* cite journal
first = Sir G.
last = Taylor
authorlink = Geoffrey Ingram Taylor
year = 1950
title = The Formation of a Blast Wave by a Very Intense Explosion. I. Theoretical Discussion
journal = Proc. Roy. Soc. A
volume = 201
pages = 159–174
doi = 10.1098/rspa.1950.0049

* cite journal
first = Sir G.
last = Taylor
authorlink = Geoffrey Ingram Taylor
year = 1950
title = The Formation of a Blast Wave by a Very Intense Explosion. II. The Atomic Explosion of 1945
journal = Proc. Roy. Soc. A
volume = 201
pages = 175–186
doi = 10.1098/rspa.1950.0050

External links