Misplaced Pages

Buckingham π theorem: Difference between revisions

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Browse history interactively← Previous editContent deleted Content addedVisualWikitext
Revision as of 01:31, 18 October 2006 edit18.38.6.40 (talk)No edit summary← Previous edit Latest revision as of 06:56, 9 October 2024 edit undoWikiCleanerBot (talk | contribs)Bots925,164 editsm v2.05b - Bot T20 CW#61 - Fix errors for CW project (Reference before punctuation - Title linked in text)Tag: WPCleaner 
(382 intermediate revisions by more than 100 users not shown)
Line 1: Line 1:
{{More citations needed|date=April 2022}}
The '''Buckingham π theorem''' is a key theorem in ]. 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. More accurately, the number of dimensionless terms that can be formed, ''p'', is equal to the nullity of the dimensional matrix. For the purposes of the experimenter, different systems which share the same description in terms of these ]s are equivalent.
{{Short description|Theorem in dimensional analysis}}
{{DISPLAYTITLE:Buckingham {{pi}} theorem}}
] circa 1886]]
In ], ], and ], the '''Buckingham {{pi}} theorem''' is a key ] in ]. It is a formalisation of ]. Loosely, the theorem states that if there is a physically meaningful equation involving a certain number ''n'' of physical variables, then the original equation can be rewritten in terms of a set of ''p''&nbsp;=&nbsp;''n''&nbsp;&minus;&nbsp;''k'' dimensionless parameters {{pi}}<sub>1</sub>, {{pi}}<sub>2</sub>, ..., {{pi}}<sub>''p''</sub> constructed from the original variables, where ''k'' is the number of physical dimensions involved; it is obtained as the ] of a particular ].


The theorem provides a method for computing sets of dimensionless parameters from the given variables, or ], even if the form of the equation is still unknown.
In mathematical terms, if we have a physically meaningful equation such as


The Buckingham {{pi}} theorem indicates that validity of the ] does not depend on a specific ] system. A statement of this theorem is that any physical law can be expressed as an ] involving only dimensionless combinations (ratios or products) of the variables linked by the law (for example, pressure and volume are linked by ] – they are ]). If the dimensionless combinations' values changed with the systems of units, then the equation would not be an identity, and the theorem would not hold.
:<math>f(q_1,q_2,\ldots,q_n)=0\,\!</math>


== History ==
where the ''q<sub>i</sub>''&nbsp; are the ''n''&nbsp; physical variables, and they are expressed in terms of ''k''&nbsp; independent physical units, then the above equation can be restated as


Although named for ], the {{pi}} theorem was first proved by the French mathematician ] in 1878.<ref>{{cite journal|last=Bertrand|first=J.|year=1878|title=Sur l'homogénéité dans les formules de physique|url=https://archive.org/details/comptesrendusheb86acad|journal=Comptes Rendus|volume=86|number=15|pages=916–920}}</ref> Bertrand considered only special cases of problems from electrodynamics and heat conduction, but his article contains, in distinct terms, all the basic ideas of the modern proof of the theorem and clearly indicates the theorem's utility for modelling physical phenomena. The technique of using the theorem ("the method of dimensions") became widely known due to the works of ]. The first application of the {{pi}} theorem ''in the general case''<ref group=note>When in applying the {{pi}}–theorem there arises an ''arbitrary function'' of dimensionless numbers.</ref> to the dependence of pressure drop in a pipe upon governing parameters probably dates back to 1892,<ref>{{cite journal|last=Rayleigh|year=1892|title=On the question of the stability of the flow of liquids|url=http://gidropraktikum.narod.ru/Rayleigh-1892.djvu|journal=Philosophical Magazine|volume=34|issue=206|pages=59–70|doi=10.1080/14786449208620167}}</ref> a heuristic proof with the use of series expansions, to 1894.<ref>{{cite book|last=Strutt|first=John William|year=1896|title=The Theory of Sound|url=https://archive.org/details/theorysound05raylgoog|publisher=Macmillan|volume=II|edition=2nd}}</ref>
:<math>F(\pi_1,\pi_2,\ldots,\pi_p)=0\,\!</math>


Formal generalization of the {{pi}} theorem for the case of arbitrarily many quantities was given first by {{ill|A. Vaschy|fr|Aimé Vaschy}} in 1892,<ref>Quotes from Vaschy's article with his statement of the pi–theorem can be found in: {{cite journal|last=Macagno|first=E. O.|year=1971|title=Historico-critical review of dimensional analysis| url=http://gidropraktikum.narod.ru/Macagno-1971.djvu|journal=Journal of the Franklin Institute|issue=6|volume=292|pages=391–402|doi=10.1016/0016-0032(71)90160-8}}</ref><ref>{{Cite journal |last=De A. Martins |first=Roberto |date=1981 |title=The origin of dimensional analysis |url=https://linkinghub.elsevier.com/retrieve/pii/0016003281904750 |journal=Journal of the Franklin Institute |language=en |volume=311 |issue=5 |pages=331–337 |doi=10.1016/0016-0032(81)90475-0}}</ref> then in 1911—apparently independently—by both A. Federman<ref>{{cite journal|last=Федерман|first=А.|year=1911|title=О некоторых общих методах интегрирования уравнений с частными производными первого порядка|url=http://gidropraktikum.narod.ru/Federman.djvu|journal=Известия Санкт-Петербургского политехнического института императора Петра Великого. Отдел техники, естествознания и математики|issue=1|volume=16|pages=97–155}} (Federman A., On some general methods of integration of first-order partial differential equations, Proceedings of the Saint-Petersburg polytechnic institute. Section of technics, natural science, and mathematics)</ref> and ],<ref>{{cite journal|last=Riabouchinsky|first=D.|year=1911|title=Мéthode des variables de dimension zéro et son application en aérodynamique|url=http://gidropraktikum.narod.ru/Riabouchinsky-Aerophile-1911.djvu|journal=L'Aérophile|pages=407–408}}</ref> and again in 1914 by Buckingham.{{sfn|Buckingham|1914}} It was Buckingham's article that introduced the use of the symbol "<math>\pi_i</math>" for the dimensionless variables (or parameters), and this is the source of the theorem's name.
where the &pi;<sub>i</sub> are dimensionless parameters constructed from the ''q<sub>i</sub>''&nbsp; by ''k'' = ''n'' &minus; ''p''&nbsp; equations of the form


== Statement ==
:<math>\pi_i=q_1^{m_1}\,q_2^{m_2}\ldots q_n^{m_n}</math>
More formally, the number <math>p</math> of dimensionless terms that can be formed is equal to the ] of the ], and <math>k</math> is the ]. For experimental purposes, different systems that share the same description in terms of these ]s are equivalent.


In mathematical terms, if we have a physically meaningful equation such as
where the exponents ''m<sub>i</sub>''&nbsp; are rational numbers. The use of the &pi;<sub>i</sub> as the dimensionless parameters was introduced by ] in his original 1914 paper on the subject from which the theorem draws its name.
<math display=block>f(q_1,q_2,\ldots,q_n)=0,</math>
where <math>q_1, \ldots, q_n</math> are any <math>n</math> physical variables, and there is a maximal dimensionally independent subset of size <math>k</math>,<ref group=note>A dimensionally independent set of variables is one for which the only exponents <math>q_1^{a_1} \, q_2^{a_2} \cdots q_k^{a_k}</math> yielding a dimensionless quantity are <math>a_1 = a_2 = \cdots = 0</math>. This is precisely the notion of ].</ref> then the above equation can be restated as
<math display=block>F(\pi_1,\pi_2,\ldots,\pi_p)=0,</math>
where <math>\pi_1, \ldots, \pi_p</math> are dimensionless parameters constructed from the <math>q_i</math> by <math>p = n - k</math> dimensionless equations — the so-called ''Pi groups'' — of the form
<math display=block>\pi_i=q_1^{a_1}\,q_2^{a_2} \cdots q_n^{a_n},</math>
where the exponents <math>a_i</math> are rational numbers. (They can always be taken to be integers by redefining <math>\pi_i</math> as being raised to a power that clears all denominators.) If there are <math>\ell</math> fundamental units in play, then <math>p \geq n - \ell</math>.


== Significance ==
Most importantly, the Buckingham &pi; 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'.
The Buckingham {{pi}} theorem provides a method for computing sets of dimensionless parameters from given variables, even if the form of the equation remains unknown. However, the choice of dimensionless parameters is not unique; Buckingham's theorem only provides a way of generating sets of dimensionless parameters and does not indicate the most "physically meaningful".


Two systems for which these parameters coincide are called ''similar'' (as with ], 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. Most importantly, Buckingham's theorem describes the relation between the number of variables and fundamental dimensions.
== Proving the &pi; theorem ==


== Proof ==
Proofs of the &pi; theorem often begin by considering the space of fundamental and derived physical units as a ], 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.


For simplicity, it will be assumed that the space of fundamental and derived physical units forms a ] over the ]s, 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:
Making the physical units match across sets of physical equations can then be regarded as imposing linear constraints in the physical unit vector space.
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 ] <math>g</math> has units of <math>\mathsf{L} / \mathsf{T}^2 = \mathsf{L}^1 \mathsf{T}^{-2}</math> (length over time squared), so it is represented as the vector <math>(1, -2)</math> with respect to the basis of fundamental units (length, time). We could also require that exponents of the fundamental units be rational numbers and modify the proof accordingly, in which case the exponents in the pi groups can always be taken as rational numbers or even integers.


=== Rescaling units ===
The &pi;-theorem describes how every physically meaningful equation involving ''n'' variables can be equivalently rewritten as an equation of ''n'' &minus; ''k'' dimensionless parameters, where ''k'' is the number of fundamental units used. Furthermore, and most importantly, it provides a method for computing these dimensionless parameters from the given variables, even if the form of the equation is still unknown.
Suppose we have quantities <math>q_1, q_2, \dots, q_n</math>, where the units of <math>q_i</math> contain length raised to the power <math>c_i</math>. If we originally measure length in meters but later switch to centimeters, then the numerical value of <math>q_i</math> would be rescaled by a factor of <math>100^{c_i}</math>. Any physically meaningful law should be invariant under an arbitrary rescaling of every fundamental unit; this is the fact that the pi theorem hinges on.


=== Formal proof ===
Two systems for which these parameters coincide are called ''similar''; 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.
Given a system of <math>n</math> dimensional variables <math>q_1, \ldots, q_n</math> in <math>\ell</math> fundamental (basis) dimensions, the ''dimensional matrix'' is the <math>\ell \times n</math> matrix <math>M</math> whose <math>\ell</math> rows correspond to the fundamental dimensions and whose <math>n</math> columns are the dimensions of the variables: the <math>(i, j)</math>th entry (where <math>1 \leq i \leq \ell</math> and <math>1 \leq j \leq n</math>) is the power of the <math>i</math>th fundamental dimension in the <math>j</math>th variable.
The matrix can be interpreted as taking in a combination of the variable quantities and giving out the dimensions of the combination in terms of the fundamental dimensions. So the <math>\ell \times 1</math> (column) vector that results from the multiplication
<math display=block>M\begin{bmatrix}a_1\\ \vdots \\ a_n\end{bmatrix}</math>
consists of the units of
<math display=block>q_1^{a_1}\,q_2^{a_2}\cdots q_n^{a_n}</math>
in terms of the <math>\ell</math> fundamental independent (basis) units.<ref group=note name name=ProofDetailsOfMatrixWithExample />


If we rescale the <math>i</math>th fundamental unit by a factor of <math>\alpha_i</math>, then <math>q_j</math> gets rescaled by <math>\alpha_1^{-m_{1j}}\, \alpha_2^{-m_{2j}} \cdots \alpha_\ell^{-m_{\ell j}}</math>, where <math>m_{ij}</math> is the <math>(i, j)</math>th entry of the dimensional matrix. In order to convert this into a linear algebra problem, we take ] (the base is irrelevant), yielding <math display=block>\begin{bmatrix} \log{q_1} \\ \vdots \\ \log{q_n} \end{bmatrix} \mapsto \begin{bmatrix} \log{q_1} \\ \vdots \\ \log{q_n} \end{bmatrix} - M^\operatorname{T} \begin{bmatrix} \log{\alpha_1} \\ \vdots \\ \log{\alpha_\ell} \end{bmatrix},</math> which is an ] of <math>\mathbb{R}^\ell</math> on <math>\mathbb{R}^n</math>. We define a physical law to be an arbitrary function <math>f \colon (\mathbb{R}^+)^n \to \mathbb{R}</math> such that <math>(q_1, q_2, \dots, q_n)</math> is a permissible set of values for the physical system when <math>f(q_1, q_2, \dots, q_n) = 0</math>. We further require <math>f</math> to be invariant under this action. Hence it descends to a function <math>F \colon \mathbb{R}^n / \operatorname{im}{M^\operatorname{T}} \to \mathbb{R}</math>. All that remains is to exhibit an isomorphism between <math>\mathbb{R}^n/\operatorname{im}{M^\operatorname{T}}</math> and <math>\mathbb{R}^p</math>, the (log) space of pi groups <math>(\log{\pi_1}, \log{\pi_2}, \dots, \log{\pi_p})</math>.
The &pi;-theorem uses ]: the space of all possible physical units can be seen as a ] over the ]s if we represent a unit as the set of exponents needed for the fundamental units (with a power of zero if the particular fundamental unit is not present). Multiplication of physical units is then represented by vector addition within this vector space. The algorithm of the &pi;-theorem is essentially a ] carried out in this vector space.

We construct an <math>n \times p</math> matrix <math>K</math> whose columns are a basis for <math>\ker{M}</math>. It tells us how to embed <math>\mathbb{R}^p</math> into <math>\mathbb{R}^n</math> as the kernel of <math>M</math>. That is, we have an ]
:<math> 0 \to \mathbb{R}^p \xrightarrow{\ K \ } \mathbb{R}^n \xrightarrow{\ M\ } \mathbb{R}^\ell.</math>
Taking tranposes yields another exact sequence
:<math>\mathbb{R}^\ell \xrightarrow{\ M^\operatorname{T}\ } \mathbb{R}^n \xrightarrow{\ K^\operatorname{T}\ } \mathbb{R}^p \to 0.</math>

The ] produces the desired isomorphism, which sends the ] <math>v + M^\operatorname{T} \mathbb{R}^\ell</math> to <math>K^\operatorname{T} v</math>. This corresponds to rewriting the tuple <math>(\log q_1, \log q_2, \dots, \log q_n)</math> into the pi groups <math>(\log\pi_1, \log\pi_2, \dots, \log\pi_p)</math> coming from the columns of <math>K</math>.

The ] defines seven base units, which are the ], ], ], ], ], ] and ]. It is sometimes advantageous to introduce additional base units and techniques to refine the technique of dimensional analysis. (See ] and reference.<ref name="SCADE2006">{{cite book|last1=Schlick|first1=R.|last2=Le Sergent|first2=T.|title=Computer Safety, Reliability, and Security |chapter=Checking SCADE Models for Correct Usage of Physical Units |date=2006|publisher=Springer|location=Berlin|volume=4166|pages= 358–371|doi= 10.1007/11875567_27|series=Lecture Notes in Computer Science|isbn=978-3-540-45762-6}}</ref>)


== Examples == == Examples ==


=== The simple pendulum=== === Speed ===
This example is elementary but serves to demonstrate the procedure.


Suppose a car is driving at 100&nbsp;km/h; how long does it take to go 200&nbsp;km?
We wish to determine the period ''T''&nbsp; of small oscillations in a simple pendulum. It will be assumed that it is a function of the length ''L''&nbsp;, the mass ''M''&nbsp;, 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


This question considers <math>n = 3</math> dimensioned variables: distance <math>d,</math> time <math>t,</math> and speed <math>v,</math> and we are seeking some law of the form <math>t = \operatorname{Duration}(v, d).</math> Any two of these variables are dimensionally independent, but the three taken together are not. Thus there is <math>p = n - k = 3 - 2 = 1</math> dimensionless quantity.
:<math>f(T,M,L,g) = 0.\,</math>


The dimensional matrix is
There are only three fundamental physical units in this equation: mass, time, and length. Thus we need only 4−3=1 dimensionless parameter, denoted &pi;, and the model can be re-expressed as
<math display=block>M = \begin{bmatrix}
1 & 0 & \;\;\;1\\
0 & 1 & -1
\end{bmatrix}</math>
in which the rows correspond to the basis dimensions <math>L</math> and <math>T,</math> and the columns to the considered dimensions <math>L, T, \text{ and } V,</math> where the latter stands for the speed dimension. The elements of the matrix correspond to the powers to which the respective dimensions are to be raised. For instance, the third column <math>(1, -1),</math> states that <math>V = L^0 T^0 V^1,</math> represented by the column vector <math>\mathbf{v}=,</math> is expressible in terms of the basis dimensions as <math>V = L^1 T^{-1} = L/T,</math> since <math>M\mathbf{v} = .</math>


For a dimensionless constant <math>\pi=L^{a_1}T^{a_2}V^{a_3},</math> we are looking for vectors <math>\mathbf{a}=</math> such that the matrix-vector product <math>M \mathbf{a}</math> equals the zero vector <math>.</math> In linear algebra, the set of vectors with this property is known as the ] (or nullspace) of the dimensional matrix. In this particular case its kernel is one-dimensional. The dimensional matrix as written above is in ], so one can read off a non-zero kernel vector to within a multiplicative constant:
:<math>f(\pi) = 0\,</math>
<math display=block>\mathbf{a} = \begin{bmatrix} -1\\ \;\;\;1\\ \;\;\;1\\ \end{bmatrix}.</math>


If the dimensional matrix were not already reduced, one could perform ] on the dimensional matrix to more easily determine the kernel. It follows that the dimensionless constant, replacing the dimensions by the corresponding dimensioned variables, may be written:
where &pi; is given by
<math display=block>\pi = d^{-1}t^1v^1 = tv/d.</math>


Since the kernel is only defined to within a multiplicative constant, the above dimensionless constant raised to any arbitrary power yields another (equivalent) dimensionless constant.
:{|
|-
| <math>\pi\,\!</math> || <math>=(T)^{m_1}(M)^{m_2}(L)^{m_3}(g)^{m_4}\,</math>
|-
| || <math> =(T)^{m_1}(M)^{m_2}(L)^{m_3}(L/T^2)^{m_4}\,</math>
|}


Dimensional analysis has thus provided a general equation relating the three physical variables:
for some values of m<sub>1</sub>&hellip;m<sub>4</sub>. With a little thought or experimentation it can be found that only
<math display=block>F(\pi)=0,</math>
or, letting <math>C</math> denote a ] of function <math>F,</math>
<math display=block>\pi=C,</math>
which can be written in the desired form (which recall was <math>t = \operatorname{Duration}(v, d)</math>) as
<math display=block>t = C\frac{d}{v}.</math>


The actual relationship between the three variables is simply <math>d = vt.</math> In other words, in this case <math>F</math> has one physically relevant root, and it is unity. The fact that only a single value of <math>C</math> will do and that it is equal to 1 is not revealed by the technique of dimensional analysis.
:{|
|-
| <math>\pi\,\!</math> || <math>=(T)^2(M)^0(L)^{-1}(L/T^2)^1\,</math>
|-
| || <math> =gT^2/L\,</math>
|}


]
(or some power thereof) satisfies this requirement. Note that if m<sub>2</sub> were non-zero there would be no way to cancel the M value&mdash;therefore m<sub>2</sub> ''must'' be zero.
Dimensional analysis has allowed us to conclude that the period of the pendulum is not a function of its mass.


===The simple pendulum===
The model can now be expressed as
We wish to determine the period <math>T</math> of ]. It will be assumed that it is a function of the length <math>L,</math> the mass <math>M,</math> and the ] on the surface of the Earth <math>g,</math> which has dimensions of length divided by time squared. The model is of the form
<math display=block>f(T,M,L,g) = 0.</math>


(Note that it is written as a relation, not as a function: <math>T</math> is not written here as a function of <math>M, L, \text{ and } g.</math>)
:<math>f(gT^2/L) = 0.\,</math>


Period, mass, and length are dimensionally independent, but acceleration can be expressed in terms of time and length, which means the four variables taken together are not dimensionally independent. Thus we need only <math>p = n - k = 4 - 3 = 1</math> dimensionless parameter, denoted by <math>\pi,</math> and the model can be re-expressed as
Assuming the zeroes of ''f''&nbsp; are discrete, we can say
<math display=block>F(\pi) = 0,</math>
''gT<sup>2</sup>/L = K<sub>n</sub>''&nbsp; where ''K<sub>n</sub>''&nbsp; is the nth zero. If there is only one zero, then ''gT<sup>2</sup>/L=K''&nbsp;. 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&pi;<sup>2</sup>''&nbsp;.
where <math>\pi</math> is given by
<math display=block>\pi = T^{a_1}M^{a_2}L^{a_3}g^{a_4}</math>
for some values of <math>a_1, a_2, a_3, a_4.</math>


The dimensions of the dimensional quantities are:
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.
<math display=block>T = t, M = m, L = \ell, g = \ell/t^2.</math>


The dimensional matrix is:
===The Atomic bomb===
<math display=block>\mathbf{M} = \begin{bmatrix}
1 & 0 & 0 & -2\\
0 & 1 & 0 & 0\\
0 & 0 & 1 & 1
\end{bmatrix}.</math>


(The rows correspond to the dimensions <math>t, m,</math> and <math>\ell,</math> and the columns to the dimensional variables <math>T, M, L, \text{ and } g.</math> For instance, the 4th column, <math>(-2, 0, 1),</math> states that the <math>g</math> variable has dimensions of <math>t^{-2}m^0 \ell^1.</math>)
In 1941, Sir ] used dimensional analysis to estimate the energy released in an ] explosion (Taylor, 1950a,b). The first atomic bomb was detonated near ], ] on July 16, 1945. In 1947, movies of the explosion were declassified, allowing
Sir 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.


We are looking for a kernel vector <math>a = \left</math> such that the matrix product of <math>\mathbf{M}</math> on <math>a</math> yields the zero vector <math>.</math> The dimensional matrix as written above is in reduced row echelon form, so one can read off a kernel vector within a multiplicative constant:
Taylor supposed that the description of the process was adequately described by five
<math display=block>a = \begin{bmatrix}2\\ 0 \\ -1 \\ 1\end{bmatrix}.</math>
physical quantities, the time ''t''&nbsp; since the detonation, the energy ''E''&nbsp;
which is released at a single point in space at detonation, the radius ''R''&nbsp; of the shock wave at time ''t''&nbsp;, the atmospheric pressure ''p''&nbsp; and the ambient
density &rho;. 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


Were it not already reduced, one could perform ] on the dimensional matrix to more easily determine the kernel. It follows that the dimensionless constant may be written:
:<math>\pi_0=R\,\left(\frac{\rho}{Et^2}\right)^{1/5}</math>
<math display=block>\begin{align}
and
\pi &= T^2M^0L^{-1}g^1\\
:<math>\pi_1=p\,\left(\frac{t^6}{E^2\rho^3}\right)^{1/5}.</math>
&= gT^2/L
The process can now be described by an equation of the form
:<math>f(\pi_0,\pi_1)=0,\,</math> \end{align}.</math>
In fundamental terms:
or, equivalently
:<math>R=\left(\frac{Et^2}{\rho}\right)^{1/5}g(\pi_1),</math> <math display=block>\pi = (t)^2 (m)^0 (\ell)^{-1} \left(\ell/t^2\right)^1 = 1,</math>
which is dimensionless. Since the kernel is only defined to within a multiplicative constant, if the above dimensionless constant is raised to any arbitrary power, it will yield another equivalent dimensionless constant.
where ''g''(&pi;<sub>1</sub>) is some function of &pi;<sub>1</sub>. 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 &pi;<sub>1</sub> 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


In this example, three of the four dimensional quantities are fundamental units, so the last (which is <math>g</math>) must be a combination of the previous.
:<math>R\approx\left(\frac{Et^2}{\rho}\right)^{1/5}.</math>
Note that if <math>a_2</math> (the coefficient of <math>M</math>) had been non-zero then there would be no way to cancel the <math>M</math> value; therefore <math>a_2</math> {{em|must}} be zero. Dimensional analysis has allowed us to conclude that the period of the pendulum is not a function of its mass <math>M.</math> (In the 3D space of powers of mass, time, and distance, we can say that the vector for mass is linearly independent from the vectors for the three other variables. Up to a scaling factor, <math>\vec g + 2 \vec T - \vec L</math> is the only nontrivial way to construct a vector of a dimensionless parameter.)


The model can now be expressed as:
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)
<math display=block>F\left(gT^2/L\right) = 0.</math>

Then this implies that <math>gT^2/L = C_i</math> for some zero <math>C_i</math> of the function <math>F.</math> If there is only one zero, call it <math>C,</math> then <math>gT^2/L = C.</math> 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 <math>C = 4\pi^2.</math>

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 as the ].

=== Electric power ===
To demonstrate the application of the {{pi}} theorem, consider the ] consumption of a ] with a given shape.
The power, ''P'', in dimensions , is a function of the ], ''ρ'' , and the ] of the fluid to be stirred, ''μ'' , as well as the size of the stirrer given by its ], ''D'' , and the ] of the stirrer, ''n'' . Therefore, we have a total of ''n'' = 5 variables representing our example. Those ''n'' = 5 variables are built up from ''k'' = 3 independent dimensions, e.g., length: L (] units: ]), time: T (]), and mass: M (]).

According to the {{pi}}-theorem, the ''n'' = 5 variables can be reduced by the ''k'' = 3 dimensions to form ''p'' = ''n'' − ''k'' = 5 − 3 = 2 independent dimensionless numbers. Usually, these quantities are chosen as <math display=inline>\mathrm{Re} = {\frac{\rho n D^2}{\mu}}</math>, commonly named the ] which describes the fluid flow regime, and <math display=inline>N_\mathrm{p} = \frac{P}{\rho n^3 D^5}</math>, the ], which is the dimensionless description of the stirrer.

Note that the two dimensionless quantities are not unique and depend on which of the ''n'' = 5 variables are chosen as the ''k'' = 3 dimensionally independent basis variables, which, in this example, appear in both dimensionless quantities. The Reynolds number and power number fall from the above analysis if <math display=inline>\rho</math>, ''n'', and ''D'' are chosen to be the basis variables. If, instead, <math display=inline>\mu</math>, ''n'', and ''D'' are selected, the Reynolds number is recovered while the second dimensionless quantity becomes <math display=inline>N_\mathrm{Rep} = \frac{P}{\mu D^3 n^2}</math>. We note that <math display=inline>N_\mathrm{Rep}</math> is the product of the Reynolds number and the power number.

===Other examples===

An example of dimensional analysis can be found for the case of the mechanics of a thin, solid and parallel-sided rotating disc. There are five variables involved which reduce to two non-dimensional groups. The relationship between these can be determined by numerical experiment using, for example, the finite element method.<ref>{{cite web|last1=Ramsay|first1=Angus|title=Dimensional Analysis and Numerical Experiments for a Rotating Disc|url=http://www.ramsay-maunder.co.uk/knowledge-base/technical-notes/dimensional-analysis--numerical-experiments-for-a-rotating-disc/|website=Ramsay Maunder Associates|access-date=15 April 2017}}</ref>

The theorem has also been used in fields other than physics, for instance in ].<ref>{{cite journal|last1=Blondeau|first1=J.|title=The influence of field size, goal size and number of players on the average number of goals scored per game in variants of football and hockey: the Pi-theorem applied to team sports|journal=Journal of Quantitative Analysis in Sports|date=2020|volume=17|issue=2|pages=145–154|doi=10.1515/jqas-2020-0009|s2cid=224929098|url=https://doi.org/10.1515/jqas-2020-0009}}</ref>


==See also== ==See also==
{{Portal|Mathematics|Physics}}
* ]
* ] * ]
* ]
* ]
* ] * ]
* ]

== External links ==
* http://www.math.ntnu.no/~hanche/notes/buckingham/
* http://www.isd.uni-stuttgart.de/~rudolph/
* http://scienceworld.wolfram.com/physics/BuckinghamsPiTheorem.html


== References == == References ==


=== Notes ===
* {{cite journal
{{reflist|group=note|refs=
| last=Buckingham
<ref group=note name name=ProofDetailsOfMatrixWithExample>If these basis units are <math>b_1, \ldots, b_\ell</math> and if the units of <math>q_j = m_{1j} b_1 + \cdots + m_{\ell j} b_\ell</math> for every <math>1 \leq j \leq n</math>, then
| first=E.
<math display="block">M = \begin{bmatrix}
| authorlink=Edgar Buckingham
m_{11} & \cdots & m_{1j} & \cdots & m_{1n} \\
| title=On physically similar systems; illustrations of the use of dimensional equations
\vdots & & \vdots & & \vdots \\
| journal=Phys. Rev.
m_{\ell 1} & \cdots & m_{\ell j} & \cdots & m_{\ell n} \\
| year=1914
\end{bmatrix}</math>
| volume= 4
so that, for instance, the units of <math>q_1</math> in terms of these basis units are
| pages=345-376
<math display="block">
| url=http://prola.aps.org/abstract/PR/v4/i4/p345_1}}
M \left(\left^{\operatorname{T}}\right) =
\begin{bmatrix}
m_{11} \\
\vdots \\
m_{\ell 1} \\
\end{bmatrix}.</math>
For a concrete example, suppose that the <math>\ell = 2</math> fundamental units are ]s <math>b_1 = m</math> and ]s <math>b_2 = s,</math> and that there are <math>n = 3</math> dimensional variables: <math>q_1 = m/s^2, q_2 = 1/m, q_3 = s / m.</math> By definition of vector addition and scalar multiplication of units,
<math display=block>
q_1 = ms^{-2} = 1 m + (-2) s,\quad q_2 = m^{-1} = (-1) m + 0 s,\quad\text{and}\quad q_3 = m^{-1} s = (-1) m + 1 s,
</math>
so that
<math display="block">M = \begin{bmatrix}
m_{11} & m_{12} & m_{13} \\
m_{21} & m_{22} & m_{23} \\
\end{bmatrix}
= \begin{bmatrix}
1 & -1 & -1 \\
-2 & 0 & 1 \\
\end{bmatrix}.</math>
By definition, the dimensionless variables are those whose units are <math>m^0 s^0,</math> which are exactly the vectors in
<math display="block">\ker M = \operatorname{span} \left\{ ^{\operatorname{T}} \right\} = \left\{ \left(q_1 - q_2 + 2 q_3\right)^s : s \in \Q \right\}.</math>
This can be verified by a direct computation:
<math display=block>q_1 - q_2 + 2 q_3 = \left(m s^{-2}\right)^1 + \left(m^{-1}\right)^{-1} + \left(s m^{-1}\right)^2 = m^1 s^{-2} + m^1 + m^{-2} s^2 = m^{1+1+(-2)} s^{-2 + 0 + 2} = m^0 s^0, </math>
which is indeed dimensionless.
Consequently, if some physical law states that <math>q_1, q_2, q_3</math> are necessarily related by a (presumably unknown) equation of the form <math>f\left(q_1, q_2, q_3\right) = 0</math> for some (unknown) function <math>f</math> with <math>\operatorname{domain}(f) \subseteq \R^3</math> (that is, the tuple <math>\left(q_1, q_2, q_3\right)</math> is necessarily a zero of <math>f</math>), then there exists some (also unknown) function <math>F : \R^1 \to \R</math> that depends on only <math>p = 3 - 2 = 1</math> variable, the dimensionless variable <math>\pi_1 := q_1 - q_2 + 2 q_3 = q_1 q_3^2 / q_2</math> (or any non-zero rational power <math>\hat{\pi}_1 := \pi_1^s</math> of <math>\pi_1,</math> where <math>0 \neq s \in \Q</math>), such that <math>F\left(\pi_1\right) = 0</math> holds (if <math>\hat{\pi}_1 := \pi_1^s</math> is used instead of <math>\pi_1</math> then <math>F</math> can be replaced with <math>\hat{F}(x) := F\left(x^{1/s}\right)</math> and once again <math>\hat{F}\left(\hat{\pi}_1\right) = 0</math> holds).
Thus in terms of the original variables, <math>F\left(q_1 q_3^2 / q_2\right) = 0</math> must hold (alternatively, if using <math>\hat{\pi}_1 := \pi_1^{1/2} = \sqrt{\pi_1}</math> for instance, then <math>\hat{F}\left(\sqrt{q_1 q_3^2 / q_2}\right) = 0</math> must hold).
In other words, the Buckingham {{pi}} theorem implies that <math>q_1 q_3^2 / q_2 \in F^{-1}(0),</math> so that if it happens to be the case that this <math>F</math> has exactly one zero, call it <math>C,</math> then the equation <math>q_1 q_3^2 / q_2 = C</math> will necessarily hold (the theorem does not give information about what the exact value of the constant <math>C</math> will be, nor does it guarantee that <math>F</math> has exactly one zero).
</ref>
}}


===Citations===
* {{cite journal
{{reflist}}
|last=Buckingham
|first=E.
|authorlink=Edgar Buckingham
|year=1915
|title=The principle of similitude
|journal=Nature
|volume=96
|pages=396-397}}


=== Bibliography ===
* {{cite journal
{{refbegin}}
|last=Buckingham
* {{cite book |author-link=Garrett Birkhoff |first=Garrett |last=Birkhoff |chapter=4. Modeling and Dimensional Analysis §63 Pi Theorem |chapter-url={{GBurl|Vj7WCgAAQBAJ|p=93}} |title=Hydrodynamics |publisher=Princeton University Press |orig-year=1960 |date=2015 |isbn=978-1-4008-7777-5 |pages=93– |edition=2nd}}
|first=E.
* {{cite web
|authorlink=Edgar Buckingham
| last=Hanche-Olsen |first=Harald | year=2004
|year=1915
| url=http://www.math.ntnu.no/~hanche/notes/buckingham/buckingham-a4.pdf
|title=Model experiments and the forms of empirical equations
| title=Buckingham's pi-theorem | publisher=NTNU | access-date=April 9, 2007
|journal=Trans. A.S.M.E
}}
|volume=37
*{{cite book
|pages=263-296}}

* {{ cite book
| last = Hart | last = Hart
| first = George W. | first = George W. |author-link=George W. Hart
| year = March 1, 1995 | date = 1995
| title = Multidimensional Analysis: Algebras and Systems for Science and Engineering | title = Multidimensional Analysis: Algebras and Systems for Science and Engineering
| publisher = Springer-Verlag | publisher = Springer
| isbn = 978-0-387-94417-3
| id = ISBN 0387944176
| url = http://www.georgehart.com/research/multanal.html}} | url = http://www.georgehart.com/research/multanal.html
}}
*{{cite book

* {{ cite book
| last = Kline | last = Kline
| first = Stephen J. | first = Stephen J.
| year = 1986 |chapter-url={{GBurl|GHvnCAAAQBAJ|p=8}} |chapter=2. Dimensional Analysis and the Pi Theorem Units and Dimensions
| year = 1986
| title = Similitude and Approximation Theory | title = Similitude and Approximation Theory
| publisher = Springer-Verlag, New York | publisher = Springer |url=https://books.google.com/books?id=GHvnCAAAQBAJ&pg=PR1
| isbn = 978-0-387-16518-9 |pages=8–35
| id = ISBN 0387165185 }}
}}
*{{cite news
| last = Hartke
| first = Jan-David
| year = 2019
| title = On Buckingham's Π-theorem
| arxiv = 1912.08744
}}
*{{cite book
| last = Wan
| first = Frederic Y.M.
| year = 1989
| title = Mathematical Models and their Analysis
| publisher = Harper & Row
| isbn = 978-0-06-046902-3
}}
*{{cite web
| last=Vignaux |first=G.A. | year=1991
| url=http://www.mcs.vuw.ac.nz/~vignaux/docs/maxent.pdf
| title=Dimensional analysis in data modelling | publisher= Victoria University of Wellington | access-date=December 15, 2005
}}
*{{cite web |first=Mike |last=Sheppard |date=2008 |url=http://www.msu.edu/~sheppa28/constants/constants.html |title=Systematic Search for Expressions of Dimensionless Constants using the NIST database of Physical Constants |archive-url=https://web.archive.org/web/20120928062000/http://www.mit.edu/~mi22295/constants/constants.html |archive-date=2012-09-28 }}
*{{cite book
| last = Gibbings
| first = J.C.
| year = 2011
| title = Dimensional Analysis
| publisher = Springer
| isbn = 978-1-84996-316-9
}}
*{{cite book
| last = Schuring
| first = Dieterich J.
| year = 1977
| title = Scale models in engineering : fundamentals and applications
| publisher = Pergamon Press, Oxford
| isbn = 978-0080208602
}}
{{refend}}


=== Original sources ===
{{refbegin}}
*{{cite journal
| last=Vaschy | first=A.
| year=1892
| title=Sur les lois de similitude en physique
| journal=]
| volume= 19 | pages=25–28
| url=https://babel.hathitrust.org/cgi/pt?id=uiug.30112098055970&view=1up&seq=31
}}
*{{cite journal
|last=Riabouchinsky | first=D.
|year=1911
|title=Мéthode des variables de dimension zéro et son application en aérodynamique
|journal=]
|pages=407–408
}}
* {{cite journal * {{cite journal
| first = Sir G. | last=Buckingham | first=E.
| author-link=Edgar Buckingham
| last = Taylor
| year=1914
| authorlink = Geoffrey Ingram Taylor
| title=On physically similar systems; illustrations of the use of dimensional equations
| year = 1950
| journal=]
| title = The Formation of a Blast Wave by a Very Intense Explosion. I. Theoretical Discussion
| volume=4 | issue=4 | pages=345–376
| journal = Proc. Roy. Soc. A
| doi=10.1103/PhysRev.4.345|url=https://babel.hathitrust.org/cgi/pt?id=uc1.31210014450082&view=1up&seq=905
| volume = 201
|bibcode = 1914PhRv....4..345B | hdl=10338.dmlcz/101743
| pages = 159-174
| hdl-access=free
}}
*{{cite journal
|last=Buckingham |first=E.
|author-link=Edgar Buckingham
|year=1915
|title=The principle of similitude
|journal=]
|volume=96 |pages=396–397
|doi=10.1038/096396d0
|bibcode = 1915Natur..96..396B
|issue=2406|s2cid=3956628
|url=https://zenodo.org/record/1429603
}}
*{{cite journal
|last=Buckingham |first=E.
|author-link=Edgar Buckingham
|year=1915
|title=Model experiments and the forms of empirical equations
|journal=]
|volume=37 |pages=–296
|doi=10.1115/1.4059754
|url=https://archive.org/details/transactionsof37amer
}}
*{{cite journal
| first = Sir G. | last = Taylor
| author-link = G. I. Taylor
| year = 1950
| title = The Formation of a Blast Wave by a Very Intense Explosion. I. Theoretical Discussion
| journal = ]
| volume = 201 | pages = 159–174
| doi = 10.1098/rspa.1950.0049
|bibcode = 1950RSPSA.201..159T
| issue = 1065 | s2cid = 54070514
}} }}
*{{cite journal

| first = Sir G. | last = Taylor
* {{cite journal
| first = Sir G. | author-link = G. I. Taylor
| last = Taylor
| authorlink = Geoffrey Ingram Taylor
| year = 1950 | year = 1950
| title = The Formation of a Blast Wave by a Very Intense Explosion. II. The Atomic Explosion of 1945 | title = The Formation of a Blast Wave by a Very Intense Explosion. II. The Atomic Explosion of 1945
| journal = Proc. Roy. Soc. A | journal = ]
| volume = 201 | volume = 201 | pages = 175–186
| doi = 10.1098/rspa.1950.0050
| pages = 175-186
|bibcode = 1950RSPSA.201..175T
| issue = 1065 | doi-access = free
}} }}
{{refend}}


==External links==
* {{ 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 }}


{{DEFAULTSORT:Buckingham Pi Theorem}}
* {{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 | accessdate=December 15 | accessyear=2005 }}



''To do: continue proof sketch...''
] ]
] ]

]
]
]
]
]

Latest revision as of 06:56, 9 October 2024

This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Buckingham π theorem" – news · newspapers · books · scholar · JSTOR (April 2022) (Learn how and when to remove this message)
Theorem in dimensional analysis
Edgar Buckingham circa 1886

In engineering, applied mathematics, and physics, the Buckingham π theorem is a key theorem in dimensional analysis. It is a formalisation of Rayleigh's method of dimensional analysis. Loosely, the theorem states that if there is a physically meaningful equation involving a certain number n of physical variables, then the original equation can be rewritten in terms of a set of p = n − k dimensionless parameters π1, π2, ..., πp constructed from the original variables, where k is the number of physical dimensions involved; it is obtained as the rank of a particular matrix.

The theorem provides a method for computing sets of dimensionless parameters from the given variables, or nondimensionalization, even if the form of the equation is still unknown.

The Buckingham π theorem indicates that validity of the laws of physics does not depend on a specific unit system. A statement of this theorem is that any physical law can be expressed as an identity involving only dimensionless combinations (ratios or products) of the variables linked by the law (for example, pressure and volume are linked by Boyle's law – they are inversely proportional). If the dimensionless combinations' values changed with the systems of units, then the equation would not be an identity, and the theorem would not hold.

History

Although named for Edgar Buckingham, the π theorem was first proved by the French mathematician Joseph Bertrand in 1878. Bertrand considered only special cases of problems from electrodynamics and heat conduction, but his article contains, in distinct terms, all the basic ideas of the modern proof of the theorem and clearly indicates the theorem's utility for modelling physical phenomena. The technique of using the theorem ("the method of dimensions") became widely known due to the works of Rayleigh. The first application of the π theorem in the general case to the dependence of pressure drop in a pipe upon governing parameters probably dates back to 1892, a heuristic proof with the use of series expansions, to 1894.

Formal generalization of the π theorem for the case of arbitrarily many quantities was given first by A. Vaschy [fr] in 1892, then in 1911—apparently independently—by both A. Federman and D. Riabouchinsky, and again in 1914 by Buckingham. It was Buckingham's article that introduced the use of the symbol " π i {\displaystyle \pi _{i}} " for the dimensionless variables (or parameters), and this is the source of the theorem's name.

Statement

More formally, the number p {\displaystyle p} of dimensionless terms that can be formed is equal to the nullity of the dimensional matrix, and k {\displaystyle k} is the rank. For experimental purposes, different systems that 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 , , q n ) = 0 , {\displaystyle f(q_{1},q_{2},\ldots ,q_{n})=0,} where q 1 , , q n {\displaystyle q_{1},\ldots ,q_{n}} are any n {\displaystyle n} physical variables, and there is a maximal dimensionally independent subset of size k {\displaystyle k} , then the above equation can be restated as F ( π 1 , π 2 , , π p ) = 0 , {\displaystyle F(\pi _{1},\pi _{2},\ldots ,\pi _{p})=0,} where π 1 , , π p {\displaystyle \pi _{1},\ldots ,\pi _{p}} are dimensionless parameters constructed from the q i {\displaystyle q_{i}} by p = n k {\displaystyle p=n-k} dimensionless equations — the so-called Pi groups — of the form π i = q 1 a 1 q 2 a 2 q n a n , {\displaystyle \pi _{i}=q_{1}^{a_{1}}\,q_{2}^{a_{2}}\cdots q_{n}^{a_{n}},} where the exponents a i {\displaystyle a_{i}} are rational numbers. (They can always be taken to be integers by redefining π i {\displaystyle \pi _{i}} as being raised to a power that clears all denominators.) If there are {\displaystyle \ell } fundamental units in play, then p n {\displaystyle p\geq n-\ell } .

Significance

The Buckingham π theorem provides a method for computing sets of dimensionless parameters from given variables, even if the form of the equation remains unknown. However, the choice of dimensionless parameters is not unique; Buckingham's theorem only provides a way of generating sets of dimensionless parameters and does not indicate 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. Most importantly, Buckingham's theorem describes the relation between the number of variables and fundamental dimensions.

Proof

For simplicity, it will be assumed that the space of fundamental and derived physical units forms a vector space over the real 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 standard gravity g {\displaystyle g} has units of L / T 2 = L 1 T 2 {\displaystyle {\mathsf {L}}/{\mathsf {T}}^{2}={\mathsf {L}}^{1}{\mathsf {T}}^{-2}} (length over time squared), so it is represented as the vector ( 1 , 2 ) {\displaystyle (1,-2)} with respect to the basis of fundamental units (length, time). We could also require that exponents of the fundamental units be rational numbers and modify the proof accordingly, in which case the exponents in the pi groups can always be taken as rational numbers or even integers.

Rescaling units

Suppose we have quantities q 1 , q 2 , , q n {\displaystyle q_{1},q_{2},\dots ,q_{n}} , where the units of q i {\displaystyle q_{i}} contain length raised to the power c i {\displaystyle c_{i}} . If we originally measure length in meters but later switch to centimeters, then the numerical value of q i {\displaystyle q_{i}} would be rescaled by a factor of 100 c i {\displaystyle 100^{c_{i}}} . Any physically meaningful law should be invariant under an arbitrary rescaling of every fundamental unit; this is the fact that the pi theorem hinges on.

Formal proof

Given a system of n {\displaystyle n} dimensional variables q 1 , , q n {\displaystyle q_{1},\ldots ,q_{n}} in {\displaystyle \ell } fundamental (basis) dimensions, the dimensional matrix is the × n {\displaystyle \ell \times n} matrix M {\displaystyle M} whose {\displaystyle \ell } rows correspond to the fundamental dimensions and whose n {\displaystyle n} columns are the dimensions of the variables: the ( i , j ) {\displaystyle (i,j)} th entry (where 1 i {\displaystyle 1\leq i\leq \ell } and 1 j n {\displaystyle 1\leq j\leq n} ) is the power of the i {\displaystyle i} th fundamental dimension in the j {\displaystyle j} th variable. The matrix can be interpreted as taking in a combination of the variable quantities and giving out the dimensions of the combination in terms of the fundamental dimensions. So the × 1 {\displaystyle \ell \times 1} (column) vector that results from the multiplication M [ a 1 a n ] {\displaystyle M{\begin{bmatrix}a_{1}\\\vdots \\a_{n}\end{bmatrix}}} consists of the units of q 1 a 1 q 2 a 2 q n a n {\displaystyle q_{1}^{a_{1}}\,q_{2}^{a_{2}}\cdots q_{n}^{a_{n}}} in terms of the {\displaystyle \ell } fundamental independent (basis) units.

If we rescale the i {\displaystyle i} th fundamental unit by a factor of α i {\displaystyle \alpha _{i}} , then q j {\displaystyle q_{j}} gets rescaled by α 1 m 1 j α 2 m 2 j α m j {\displaystyle \alpha _{1}^{-m_{1j}}\,\alpha _{2}^{-m_{2j}}\cdots \alpha _{\ell }^{-m_{\ell j}}} , where m i j {\displaystyle m_{ij}} is the ( i , j ) {\displaystyle (i,j)} th entry of the dimensional matrix. In order to convert this into a linear algebra problem, we take logarithms (the base is irrelevant), yielding [ log q 1 log q n ] [ log q 1 log q n ] M T [ log α 1 log α ] , {\displaystyle {\begin{bmatrix}\log {q_{1}}\\\vdots \\\log {q_{n}}\end{bmatrix}}\mapsto {\begin{bmatrix}\log {q_{1}}\\\vdots \\\log {q_{n}}\end{bmatrix}}-M^{\operatorname {T} }{\begin{bmatrix}\log {\alpha _{1}}\\\vdots \\\log {\alpha _{\ell }}\end{bmatrix}},} which is an action of R {\displaystyle \mathbb {R} ^{\ell }} on R n {\displaystyle \mathbb {R} ^{n}} . We define a physical law to be an arbitrary function f : ( R + ) n R {\displaystyle f\colon (\mathbb {R} ^{+})^{n}\to \mathbb {R} } such that ( q 1 , q 2 , , q n ) {\displaystyle (q_{1},q_{2},\dots ,q_{n})} is a permissible set of values for the physical system when f ( q 1 , q 2 , , q n ) = 0 {\displaystyle f(q_{1},q_{2},\dots ,q_{n})=0} . We further require f {\displaystyle f} to be invariant under this action. Hence it descends to a function F : R n / im M T R {\displaystyle F\colon \mathbb {R} ^{n}/\operatorname {im} {M^{\operatorname {T} }}\to \mathbb {R} } . All that remains is to exhibit an isomorphism between R n / im M T {\displaystyle \mathbb {R} ^{n}/\operatorname {im} {M^{\operatorname {T} }}} and R p {\displaystyle \mathbb {R} ^{p}} , the (log) space of pi groups ( log π 1 , log π 2 , , log π p ) {\displaystyle (\log {\pi _{1}},\log {\pi _{2}},\dots ,\log {\pi _{p}})} .

We construct an n × p {\displaystyle n\times p} matrix K {\displaystyle K} whose columns are a basis for ker M {\displaystyle \ker {M}} . It tells us how to embed R p {\displaystyle \mathbb {R} ^{p}} into R n {\displaystyle \mathbb {R} ^{n}} as the kernel of M {\displaystyle M} . That is, we have an exact sequence

0 R p   K   R n   M   R . {\displaystyle 0\to \mathbb {R} ^{p}\xrightarrow {\ K\ } \mathbb {R} ^{n}\xrightarrow {\ M\ } \mathbb {R} ^{\ell }.}

Taking tranposes yields another exact sequence

R   M T   R n   K T   R p 0. {\displaystyle \mathbb {R} ^{\ell }\xrightarrow {\ M^{\operatorname {T} }\ } \mathbb {R} ^{n}\xrightarrow {\ K^{\operatorname {T} }\ } \mathbb {R} ^{p}\to 0.}

The first isomorphism theorem produces the desired isomorphism, which sends the coset v + M T R {\displaystyle v+M^{\operatorname {T} }\mathbb {R} ^{\ell }} to K T v {\displaystyle K^{\operatorname {T} }v} . This corresponds to rewriting the tuple ( log q 1 , log q 2 , , log q n ) {\displaystyle (\log q_{1},\log q_{2},\dots ,\log q_{n})} into the pi groups ( log π 1 , log π 2 , , log π p ) {\displaystyle (\log \pi _{1},\log \pi _{2},\dots ,\log \pi _{p})} coming from the columns of K {\displaystyle K} .

The International System of Units defines seven base units, which are the ampere, kelvin, second, metre, kilogram, candela and mole. It is sometimes advantageous to introduce additional base units and techniques to refine the technique of dimensional analysis. (See orientational analysis and reference.)

Examples

Speed

This example is elementary but serves to demonstrate the procedure.

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

This question considers n = 3 {\displaystyle n=3} dimensioned variables: distance d , {\displaystyle d,} time t , {\displaystyle t,} and speed v , {\displaystyle v,} and we are seeking some law of the form t = Duration ( v , d ) . {\displaystyle t=\operatorname {Duration} (v,d).} Any two of these variables are dimensionally independent, but the three taken together are not. Thus there is p = n k = 3 2 = 1 {\displaystyle p=n-k=3-2=1} dimensionless quantity.

The dimensional matrix is M = [ 1 0 1 0 1 1 ] {\displaystyle M={\begin{bmatrix}1&0&\;\;\;1\\0&1&-1\end{bmatrix}}} in which the rows correspond to the basis dimensions L {\displaystyle L} and T , {\displaystyle T,} and the columns to the considered dimensions L , T ,  and  V , {\displaystyle L,T,{\text{ and }}V,} where the latter stands for the speed dimension. The elements of the matrix correspond to the powers to which the respective dimensions are to be raised. For instance, the third column ( 1 , 1 ) , {\displaystyle (1,-1),} states that V = L 0 T 0 V 1 , {\displaystyle V=L^{0}T^{0}V^{1},} represented by the column vector v = [ 0 , 0 , 1 ] , {\displaystyle \mathbf {v} =,} is expressible in terms of the basis dimensions as V = L 1 T 1 = L / T , {\displaystyle V=L^{1}T^{-1}=L/T,} since M v = [ 1 , 1 ] . {\displaystyle M\mathbf {v} =.}

For a dimensionless constant π = L a 1 T a 2 V a 3 , {\displaystyle \pi =L^{a_{1}}T^{a_{2}}V^{a_{3}},} we are looking for vectors a = [ a 1 , a 2 , a 3 ] {\displaystyle \mathbf {a} =} such that the matrix-vector product M a {\displaystyle M\mathbf {a} } equals the zero vector [ 0 , 0 ] . {\displaystyle .} In linear algebra, the set of vectors with this property is known as the kernel (or nullspace) of the dimensional matrix. In this particular case its kernel is one-dimensional. The dimensional matrix as written above is in reduced row echelon form, so one can read off a non-zero kernel vector to within a multiplicative constant: a = [ 1 1 1 ] . {\displaystyle \mathbf {a} ={\begin{bmatrix}-1\\\;\;\;1\\\;\;\;1\\\end{bmatrix}}.}

If the dimensional matrix were not already reduced, one could perform Gauss–Jordan elimination on the dimensional matrix to more easily determine the kernel. It follows that the dimensionless constant, replacing the dimensions by the corresponding dimensioned variables, may be written: π = d 1 t 1 v 1 = t v / d . {\displaystyle \pi =d^{-1}t^{1}v^{1}=tv/d.}

Since the kernel is only defined to within a multiplicative constant, the above dimensionless constant raised to any arbitrary power yields another (equivalent) dimensionless constant.

Dimensional analysis has thus provided a general equation relating the three physical variables: F ( π ) = 0 , {\displaystyle F(\pi )=0,} or, letting C {\displaystyle C} denote a zero of function F , {\displaystyle F,} π = C , {\displaystyle \pi =C,} which can be written in the desired form (which recall was t = Duration ( v , d ) {\displaystyle t=\operatorname {Duration} (v,d)} ) as t = C d v . {\displaystyle t=C{\frac {d}{v}}.}

The actual relationship between the three variables is simply d = v t . {\displaystyle d=vt.} In other words, in this case F {\displaystyle F} has one physically relevant root, and it is unity. The fact that only a single value of C {\displaystyle C} will do and that it is equal to 1 is not revealed by the technique of dimensional analysis.

The simple pendulum

We wish to determine the period T {\displaystyle T} of small oscillations in a simple pendulum. It will be assumed that it is a function of the length L , {\displaystyle L,} the mass M , {\displaystyle M,} and the acceleration due to gravity on the surface of the Earth g , {\displaystyle g,} which has dimensions of length divided by time squared. The model is of the form f ( T , M , L , g ) = 0. {\displaystyle f(T,M,L,g)=0.}

(Note that it is written as a relation, not as a function: T {\displaystyle T} is not written here as a function of M , L ,  and  g . {\displaystyle M,L,{\text{ and }}g.} )

Period, mass, and length are dimensionally independent, but acceleration can be expressed in terms of time and length, which means the four variables taken together are not dimensionally independent. Thus we need only p = n k = 4 3 = 1 {\displaystyle p=n-k=4-3=1} dimensionless parameter, denoted by π , {\displaystyle \pi ,} and the model can be re-expressed as F ( π ) = 0 , {\displaystyle F(\pi )=0,} where π {\displaystyle \pi } is given by π = T a 1 M a 2 L a 3 g a 4 {\displaystyle \pi =T^{a_{1}}M^{a_{2}}L^{a_{3}}g^{a_{4}}} for some values of a 1 , a 2 , a 3 , a 4 . {\displaystyle a_{1},a_{2},a_{3},a_{4}.}

The dimensions of the dimensional quantities are: T = t , M = m , L = , g = / t 2 . {\displaystyle T=t,M=m,L=\ell ,g=\ell /t^{2}.}

The dimensional matrix is: M = [ 1 0 0 2 0 1 0 0 0 0 1 1 ] . {\displaystyle \mathbf {M} ={\begin{bmatrix}1&0&0&-2\\0&1&0&0\\0&0&1&1\end{bmatrix}}.}

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

We are looking for a kernel vector a = [ a 1 , a 2 , a 3 , a 4 ] {\displaystyle a=\left} such that the matrix product of M {\displaystyle \mathbf {M} } on a {\displaystyle a} yields the zero vector [ 0 , 0 , 0 ] . {\displaystyle .} The dimensional matrix as written above is in reduced row echelon form, so one can read off a kernel vector within a multiplicative constant: a = [ 2 0 1 1 ] . {\displaystyle a={\begin{bmatrix}2\\0\\-1\\1\end{bmatrix}}.}

Were it not already reduced, one could perform Gauss–Jordan elimination on the dimensional matrix to more easily determine the kernel. It follows that the dimensionless constant may be written: π = T 2 M 0 L 1 g 1 = g T 2 / L . {\displaystyle {\begin{aligned}\pi &=T^{2}M^{0}L^{-1}g^{1}\\&=gT^{2}/L\end{aligned}}.} In fundamental terms: π = ( t ) 2 ( m ) 0 ( ) 1 ( / t 2 ) 1 = 1 , {\displaystyle \pi =(t)^{2}(m)^{0}(\ell )^{-1}\left(\ell /t^{2}\right)^{1}=1,} which is dimensionless. Since the kernel is only defined to within a multiplicative constant, if the above dimensionless constant is raised to any arbitrary power, it will yield another equivalent dimensionless constant.

In this example, three of the four dimensional quantities are fundamental units, so the last (which is g {\displaystyle g} ) must be a combination of the previous. Note that if a 2 {\displaystyle a_{2}} (the coefficient of M {\displaystyle M} ) had been non-zero then there would be no way to cancel the M {\displaystyle M} value; therefore a 2 {\displaystyle a_{2}} must be zero. Dimensional analysis has allowed us to conclude that the period of the pendulum is not a function of its mass M . {\displaystyle M.} (In the 3D space of powers of mass, time, and distance, we can say that the vector for mass is linearly independent from the vectors for the three other variables. Up to a scaling factor, g + 2 T L {\displaystyle {\vec {g}}+2{\vec {T}}-{\vec {L}}} is the only nontrivial way to construct a vector of a dimensionless parameter.)

The model can now be expressed as: F ( g T 2 / L ) = 0. {\displaystyle F\left(gT^{2}/L\right)=0.}

Then this implies that g T 2 / L = C i {\displaystyle gT^{2}/L=C_{i}} for some zero C i {\displaystyle C_{i}} of the function F . {\displaystyle F.} If there is only one zero, call it C , {\displaystyle C,} then g T 2 / L = C . {\displaystyle gT^{2}/L=C.} 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 C = 4 π 2 . {\displaystyle C=4\pi ^{2}.}

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 as the angle approaches zero.

Electric power

To demonstrate the application of the π theorem, consider the power consumption of a stirrer with a given shape. The power, P, in dimensions , is a function of the density, ρ , and the viscosity of the fluid to be stirred, μ , as well as the size of the stirrer given by its diameter, D , and the angular speed of the stirrer, n . Therefore, we have a total of n = 5 variables representing our example. Those n = 5 variables are built up from k = 3 independent dimensions, e.g., length: L (SI units: m), time: T (s), and mass: M (kg).

According to the π-theorem, the n = 5 variables can be reduced by the k = 3 dimensions to form p = nk = 5 − 3 = 2 independent dimensionless numbers. Usually, these quantities are chosen as R e = ρ n D 2 μ {\textstyle \mathrm {Re} ={\frac {\rho nD^{2}}{\mu }}} , commonly named the Reynolds number which describes the fluid flow regime, and N p = P ρ n 3 D 5 {\textstyle N_{\mathrm {p} }={\frac {P}{\rho n^{3}D^{5}}}} , the power number, which is the dimensionless description of the stirrer.

Note that the two dimensionless quantities are not unique and depend on which of the n = 5 variables are chosen as the k = 3 dimensionally independent basis variables, which, in this example, appear in both dimensionless quantities. The Reynolds number and power number fall from the above analysis if ρ {\textstyle \rho } , n, and D are chosen to be the basis variables. If, instead, μ {\textstyle \mu } , n, and D are selected, the Reynolds number is recovered while the second dimensionless quantity becomes N R e p = P μ D 3 n 2 {\textstyle N_{\mathrm {Rep} }={\frac {P}{\mu D^{3}n^{2}}}} . We note that N R e p {\textstyle N_{\mathrm {Rep} }} is the product of the Reynolds number and the power number.

Other examples

An example of dimensional analysis can be found for the case of the mechanics of a thin, solid and parallel-sided rotating disc. There are five variables involved which reduce to two non-dimensional groups. The relationship between these can be determined by numerical experiment using, for example, the finite element method.

The theorem has also been used in fields other than physics, for instance in sports science.

See also

References

Notes

  1. When in applying the π–theorem there arises an arbitrary function of dimensionless numbers.
  2. A dimensionally independent set of variables is one for which the only exponents q 1 a 1 q 2 a 2 q k a k {\displaystyle q_{1}^{a_{1}}\,q_{2}^{a_{2}}\cdots q_{k}^{a_{k}}} yielding a dimensionless quantity are a 1 = a 2 = = 0 {\displaystyle a_{1}=a_{2}=\cdots =0} . This is precisely the notion of linear independence.
  3. If these basis units are b 1 , , b {\displaystyle b_{1},\ldots ,b_{\ell }} and if the units of q j = m 1 j b 1 + + m j b {\displaystyle q_{j}=m_{1j}b_{1}+\cdots +m_{\ell j}b_{\ell }} for every 1 j n {\displaystyle 1\leq j\leq n} , then M = [ m 11 m 1 j m 1 n m 1 m j m n ] {\displaystyle M={\begin{bmatrix}m_{11}&\cdots &m_{1j}&\cdots &m_{1n}\\\vdots &&\vdots &&\vdots \\m_{\ell 1}&\cdots &m_{\ell j}&\cdots &m_{\ell n}\\\end{bmatrix}}} so that, for instance, the units of q 1 {\displaystyle q_{1}} in terms of these basis units are M ( [ 1   0     0 ] T ) = [ m 11 m 1 ] . {\displaystyle M\left(\left^{\operatorname {T} }\right)={\begin{bmatrix}m_{11}\\\vdots \\m_{\ell 1}\\\end{bmatrix}}.} For a concrete example, suppose that the = 2 {\displaystyle \ell =2} fundamental units are meters b 1 = m {\displaystyle b_{1}=m} and seconds b 2 = s , {\displaystyle b_{2}=s,} and that there are n = 3 {\displaystyle n=3} dimensional variables: q 1 = m / s 2 , q 2 = 1 / m , q 3 = s / m . {\displaystyle q_{1}=m/s^{2},q_{2}=1/m,q_{3}=s/m.} By definition of vector addition and scalar multiplication of units, q 1 = m s 2 = 1 m + ( 2 ) s , q 2 = m 1 = ( 1 ) m + 0 s , and q 3 = m 1 s = ( 1 ) m + 1 s , {\displaystyle q_{1}=ms^{-2}=1m+(-2)s,\quad q_{2}=m^{-1}=(-1)m+0s,\quad {\text{and}}\quad q_{3}=m^{-1}s=(-1)m+1s,} so that M = [ m 11 m 12 m 13 m 21 m 22 m 23 ] = [ 1 1 1 2 0 1 ] . {\displaystyle M={\begin{bmatrix}m_{11}&m_{12}&m_{13}\\m_{21}&m_{22}&m_{23}\\\end{bmatrix}}={\begin{bmatrix}1&-1&-1\\-2&0&1\\\end{bmatrix}}.} By definition, the dimensionless variables are those whose units are m 0 s 0 , {\displaystyle m^{0}s^{0},} which are exactly the vectors in ker M = span { [ 1 , 1 , 2 ] T } = { ( q 1 q 2 + 2 q 3 ) s : s Q } . {\displaystyle \ker M=\operatorname {span} \left\{^{\operatorname {T} }\right\}=\left\{\left(q_{1}-q_{2}+2q_{3}\right)^{s}:s\in \mathbb {Q} \right\}.} This can be verified by a direct computation: q 1 q 2 + 2 q 3 = ( m s 2 ) 1 + ( m 1 ) 1 + ( s m 1 ) 2 = m 1 s 2 + m 1 + m 2 s 2 = m 1 + 1 + ( 2 ) s 2 + 0 + 2 = m 0 s 0 , {\displaystyle q_{1}-q_{2}+2q_{3}=\left(ms^{-2}\right)^{1}+\left(m^{-1}\right)^{-1}+\left(sm^{-1}\right)^{2}=m^{1}s^{-2}+m^{1}+m^{-2}s^{2}=m^{1+1+(-2)}s^{-2+0+2}=m^{0}s^{0},} which is indeed dimensionless. Consequently, if some physical law states that q 1 , q 2 , q 3 {\displaystyle q_{1},q_{2},q_{3}} are necessarily related by a (presumably unknown) equation of the form f ( q 1 , q 2 , q 3 ) = 0 {\displaystyle f\left(q_{1},q_{2},q_{3}\right)=0} for some (unknown) function f {\displaystyle f} with domain ( f ) R 3 {\displaystyle \operatorname {domain} (f)\subseteq \mathbb {R} ^{3}} (that is, the tuple ( q 1 , q 2 , q 3 ) {\displaystyle \left(q_{1},q_{2},q_{3}\right)} is necessarily a zero of f {\displaystyle f} ), then there exists some (also unknown) function F : R 1 R {\displaystyle F:\mathbb {R} ^{1}\to \mathbb {R} } that depends on only p = 3 2 = 1 {\displaystyle p=3-2=1} variable, the dimensionless variable π 1 := q 1 q 2 + 2 q 3 = q 1 q 3 2 / q 2 {\displaystyle \pi _{1}:=q_{1}-q_{2}+2q_{3}=q_{1}q_{3}^{2}/q_{2}} (or any non-zero rational power π ^ 1 := π 1 s {\displaystyle {\hat {\pi }}_{1}:=\pi _{1}^{s}} of π 1 , {\displaystyle \pi _{1},} where 0 s Q {\displaystyle 0\neq s\in \mathbb {Q} } ), such that F ( π 1 ) = 0 {\displaystyle F\left(\pi _{1}\right)=0} holds (if π ^ 1 := π 1 s {\displaystyle {\hat {\pi }}_{1}:=\pi _{1}^{s}} is used instead of π 1 {\displaystyle \pi _{1}} then F {\displaystyle F} can be replaced with F ^ ( x ) := F ( x 1 / s ) {\displaystyle {\hat {F}}(x):=F\left(x^{1/s}\right)} and once again F ^ ( π ^ 1 ) = 0 {\displaystyle {\hat {F}}\left({\hat {\pi }}_{1}\right)=0} holds). Thus in terms of the original variables, F ( q 1 q 3 2 / q 2 ) = 0 {\displaystyle F\left(q_{1}q_{3}^{2}/q_{2}\right)=0} must hold (alternatively, if using π ^ 1 := π 1 1 / 2 = π 1 {\displaystyle {\hat {\pi }}_{1}:=\pi _{1}^{1/2}={\sqrt {\pi _{1}}}} for instance, then F ^ ( q 1 q 3 2 / q 2 ) = 0 {\displaystyle {\hat {F}}\left({\sqrt {q_{1}q_{3}^{2}/q_{2}}}\right)=0} must hold). In other words, the Buckingham π theorem implies that q 1 q 3 2 / q 2 F 1 ( 0 ) , {\displaystyle q_{1}q_{3}^{2}/q_{2}\in F^{-1}(0),} so that if it happens to be the case that this F {\displaystyle F} has exactly one zero, call it C , {\displaystyle C,} then the equation q 1 q 3 2 / q 2 = C {\displaystyle q_{1}q_{3}^{2}/q_{2}=C} will necessarily hold (the theorem does not give information about what the exact value of the constant C {\displaystyle C} will be, nor does it guarantee that F {\displaystyle F} has exactly one zero).

Citations

  1. Bertrand, J. (1878). "Sur l'homogénéité dans les formules de physique". Comptes Rendus. 86 (15): 916–920.
  2. Rayleigh (1892). "On the question of the stability of the flow of liquids". Philosophical Magazine. 34 (206): 59–70. doi:10.1080/14786449208620167.
  3. Strutt, John William (1896). The Theory of Sound. Vol. II (2nd ed.). Macmillan.
  4. Quotes from Vaschy's article with his statement of the pi–theorem can be found in: Macagno, E. O. (1971). "Historico-critical review of dimensional analysis". Journal of the Franklin Institute. 292 (6): 391–402. doi:10.1016/0016-0032(71)90160-8.
  5. De A. Martins, Roberto (1981). "The origin of dimensional analysis". Journal of the Franklin Institute. 311 (5): 331–337. doi:10.1016/0016-0032(81)90475-0.
  6. Федерман, А. (1911). "О некоторых общих методах интегрирования уравнений с частными производными первого порядка". Известия Санкт-Петербургского политехнического института императора Петра Великого. Отдел техники, естествознания и математики. 16 (1): 97–155. (Federman A., On some general methods of integration of first-order partial differential equations, Proceedings of the Saint-Petersburg polytechnic institute. Section of technics, natural science, and mathematics)
  7. Riabouchinsky, D. (1911). "Мéthode des variables de dimension zéro et son application en aérodynamique". L'Aérophile: 407–408.
  8. Buckingham 1914.
  9. Schlick, R.; Le Sergent, T. (2006). "Checking SCADE Models for Correct Usage of Physical Units". Computer Safety, Reliability, and Security. Lecture Notes in Computer Science. Vol. 4166. Berlin: Springer. pp. 358–371. doi:10.1007/11875567_27. ISBN 978-3-540-45762-6.
  10. Ramsay, Angus. "Dimensional Analysis and Numerical Experiments for a Rotating Disc". Ramsay Maunder Associates. Retrieved 15 April 2017.
  11. Blondeau, J. (2020). "The influence of field size, goal size and number of players on the average number of goals scored per game in variants of football and hockey: the Pi-theorem applied to team sports". Journal of Quantitative Analysis in Sports. 17 (2): 145–154. doi:10.1515/jqas-2020-0009. S2CID 224929098.

Bibliography

Original sources

External links

Categories: