# 9.3: Nusselt's Technique

- Page ID
- 766

The Nusselt's method is a bit more labor intensive, in that the governing equations with the boundary and initial conditions are used to determine the dimensionless parameters. In this method, the boundary conditions together with the governing equations are taken into account as opposed to Buckingham's method. A common mistake is to ignore the boundary conditions or initial conditions. The parameters that results from this process are the dimensional parameters which control the problems. An example comparing the Buckingham's method with Nusselt's method is presented. In this method, the governing equations, initial condition and boundary conditions are normalized resulting in a creation of dimensionless parameters which govern the solution. It is recommended, when the reader is out in the real world to simply abandon Buckingham's method all together. This point can be illustrated by example of flow over inclined plane. For comparison reasons Buckingham's method presented and later the results are compared with the results from Nusselt's method.

Example 9.10

Utilize the Buckingham's method to analyze a two dimensional flow in incline plane. Assume that the flow infinitely long and thus flow can be analyzed per width which is a function of several parameters. The potential parameters are the angle of inclination, \(\theta\), liquid viscosity, \(\nu\), gravity, \(g\), the height of the liquid, \(h\), the density, \(\rho\), and liquid velocity, \(U\). Assume that the flow is not affected by the surface tension (liquid), \(\sigma\). You furthermore are to assume that the flow is stable. Develop the relationship between the flow to the other parameters.

Solution 9.10

Under the assumptions in the example presentation leads to following

\[
\label{dim:eq:buckinghamNusseltCompar}

\dot{m} = f\left( \theta, \nu, g, \rho, U\right)

\]

The number of basic units is three while the number of the parameters is six thus the difference is \(6-3=3\). Those groups (or the work on the groups creation) further can be reduced the because angle \(\theta\) is dimensionless. The units of parameters can be obtained in Table 9.3 the following table.

Parameter |
Units |
Parameter |
Units |

\(\nu\) |
\(L^2t^{-1}\) |
\(g\) |
\(L^{1}t{-2}\) |

\(\dot{m}\) |
\(M\,t^{-1}L^{-1}\) |
\(\theta\) |
none |

\(U\) |
\(L^{1}t^{-1}\) |
\(\rho\) |
\(M\,L^{3}\) |

The basic units are chosen as for the time, \(U\), for the mass, \(\rho\), and for the length \(g\). Utilizing the building blocks technique provides

\[
\label{inclidePlane:ini}

\overbrace{\dfrac{M}{t\,L} }^{\dot{m}}

= \left( \overbrace{\dfrac{M}{L^3}}^{\rho} \right)^a

\left( \overbrace{\dfrac{L}{t^2} }^{g} \right)^b

\left( \overbrace{\dfrac{L}{t} }^{U} \right)^c

\]

\[
\label{inclidePlane:gov0}

\left.

\begin{array}{rrl}

\text{Mass}, M & a =& 1 \\

\text{Length}, L & -3a + b +c =& -1 \\

\text{time}, t & -2b -c =& - 1

\end{array}

\right\} \Longrightarrow

\pi_1 = \dfrac{\dot{m} \,g}{\rho\,\,U^3}

\]

\[
\label{inclidePlane:ini1}

\overbrace{\dfrac{L^2}{t} }^{\nu}

= \left( \overbrace{\dfrac{M}{L^3}}^{\rho} \right)^a

\left( \overbrace{\dfrac{L}{t^2} }^{g} \right)^b

\left( \overbrace{\dfrac{L}{t} }^{U} \right)^c

\]
The equations obtained from equation (2) are

\[
\label{inclidePlane:gov}

\left.

\begin{array}{rrl}

\text{Mass}, M & a =& 0 \\

\text{Length}, L & -3a + b +c =& 2 \\

\text{time}, t & -2b -c =& - 1

\end{array}

\right\} \Longrightarrow

\pi_2 = \dfrac{\nu \,g}{U^3}

\]
Thus governing equation and adding the angle can be written as

\[
\label{inclidePlane:}

0 = f\left(\dfrac{\dot{m} \,g}{\rho\,\,U^3} , \dfrac{\nu \,g}{U^3} ,\theta\right)

\]
The conclusion from this analysis are that the number of controlling parameters totaled in three and that the initial conditions and boundaries are irrelevant.

A small note, it is well established that the combination of angle gravity or effective body force is significant to the results. Hence, this analysis misses, at the very least, the issue of the combination of the angle gravity. Nusselt's analysis requires that the governing equations along with the boundary and initial conditions to be written. While the analytical solution for this situation exist, the parameters that effect the problem are the focus of this discussion. In Chapter 8, the Navier–Stokes equations were developed. These equations along with the energy, mass or the chemical species of the system, and second laws governed almost all cases in thermo–fluid mechanics. This author is not aware of a compelling reason that this fact should be used in this chapter. The two dimensional NS equation can obtained from equation (??) as

\[
\label{dim:eq:twoDNSx}

\begin{array} {rll}

\rho \left(\dfrac{\partial U_x}{\partial t} +

\right. & U_x\dfrac{\partial U_x}{\partial x} +

U_y \dfrac{\partial U_x}{\partial y} + \left. U_z

\dfrac{\partial U_x}{\partial z}\right) = &\\

&-\dfrac{\partial P}{\partial x} + \mu \left(\dfrac{\partial^2 U_x}{\partial x^2} +

\dfrac{\partial^2 U_x}{\partial y^2} + \dfrac{\partial^2 U_x}{\partial z^2}\right) +

\rho g\,\sin\theta

\end{array}

\]

\[
\label{dim:eq:twoDNSy}

\begin{array} {rll}

\rho \left(\dfrac{\partial U_y}{\partial t} +

\right. & U_x\dfrac{\partial U_y}{\partial x} +

U_y \dfrac{\partial U_y}{\partial y} + \left. U_z

\dfrac{\partial U_y}{\partial z}\right) = &\\

&-\dfrac{\partial P}{\partial x} + \mu \left(\dfrac{\partial^2 U_y}{\partial x^2} +

\dfrac{\partial^2 U_y}{\partial y^2} + \dfrac{\partial^2 U_y}{\partial z^2}\right) +

\rho g\,\sin\theta

\end{array}

\]
With boundary conditions

\[
\label{dim:eq:twoDNSbcBotton}

\begin{array}{l}

U_x (y=0) = U_{0x} f(x) \\

\dfrac{\partial U_x}{\partial x} (y=h) = \tau_0 f(x)

\end{array}

\]
The value \(U_0x\) and \(\tau_0\) are the characteristic and maximum values of the velocity or the shear stress, respectively. and the initial condition of

\[
\label{dim:eq:twoDNSic}

U_x (x=0) = U_{0y}\, f(y)

\]
where \( U_{0y}\) is characteristic initial velocity. These sets of equations (7-??) need to be converted to dimensionless equations. It can be noticed that the boundary and initial conditions are provided in a special form were the representative velocity multiply a function. Any function can be presented by this form.

In the process of transforming the equations into a dimensionless form associated with

some intelligent guess work. However, no assumption is made or required about whether or not the velocity, in the \(y\) direction. The only exception is that the \(y\) component of the velocity vanished on the boundary. No assumption is required about the acceleration or the pressure gradient etc.

The boundary conditions have typical velocities which can be used. The velocity is selected according to the situation or the needed velocity. For example, if the effect of the initial condition is under investigation than the characteristic of that velocity should be used. Otherwise the velocity at the bottom should be used. In that case, the boundary conditions are

\[
\label{dim:eq:twoDNScbles}

\begin{array}{l}

\dfrac{ U_x (y=0)}{ U_{0x}} = f(x) \\

\mu \dfrac{\partial U_x}{\partial x} (y=h) = \tau_0 \,g(x)

\end{array}

\]
Now it is very convenient to define several new variables:

\[
\label{dim:eq:twoDNSdefVer1}

\begin{array}{rlcrl}

\overline{U} = & \dfrac{ U_x (\overline{x} ) }{ U_{0x}}\\

where:\

\overline{x} = & \dfrac{x}{h} &\qquad& \overline{y} &= \dfrac{y}{h} \\

\end{array}

\]
The length \(h\) is chosen as the characteristic length since no other length is provided. It can be noticed that because the units consistency, the characteristic length can be used for "normalization'' (see Example 9.11). Using these definitions the boundary and initial conditions becomes

\[
\label{dim:eq:twoDNScbles1}

\qquad

\begin{array}{l}

\dfrac{ \overline{U_x} (\overline{y}=0)}{ U_{0x}} = f^{'}(\overline{x}) \\

\dfrac{h \, \mu}{U_{0x}}\,

\dfrac{\partial \overline{U_x}}{\partial \overline{x}} (\overline{y}=1) = \tau_0\, g^{'}(\overline{x})

\end{array}

\]
It commonly suggested to arrange the second part of equation (13) as

\[
\label{dim:eq:twoDNScbles2}

\dfrac{\partial \overline{U_x}}{\partial \overline{x}} (\overline{y}=1) =

\dfrac{\tau_0\,U_{0x}}{h \, \mu}\,

g^{'}(\overline{x})

\]
Where new dimensionless parameter, the shear stress number is defined as

\[
\label{dim:eq:tauDef}

\overline{\tau_0} = \dfrac{\tau_0\,U_{0x}}{h \, \mu}

\]
With the new definition equation (14) transformed into

\[
\label{dim:eq:twoDNScbles3}

\dfrac{\partial \overline{U_x}}{\partial \overline{x}} (\overline{y}=1) =

\overline{\tau_0} \, g^{'}(\overline{x})

\]

Example 9.11

Non–dimensionalize the following boundary condition. What are the units of the coefficient in front of the variables, \(x\). What are relationship of the typical velocity, \(U_0\) to \(U_{max}\)?

\[
\label{twoDNSbc:bc}

U_x (y = h) = U_0 \left( a\,x^2 + b\,xp(x) \right)

\]

Solution 9.11

The coefficients \(a\) and \(b\) multiply different terms and therefore must have different units. The results must be unitless thus \(a\)

\[
\label{twoDNSbc:aDef1}

L^0 = a \, \overbrace{ {L^2} }^{x^2}

\Longrightarrow a = \left[ \dfrac{1}{L^2} \right]

\]

\[
\label{twoDNSbc:bDef1}

{L}^{0} = b \, xp\left( h\,\dfrac{x}{h}\right) =

b \, xp\left( h \right)

\, xp\left( \dfrac{x}{h}\right) =

b \, xp\left( h \right)

\, xp\left( \overline{x}\right)

\]
Hence defining

\[
\label{twoDNSbc:bDef12}

\overline{b} = \dfrac{1}{xp{h}}

\]
With the new coefficients for both terms and noticing that \(y=h\longrightarrow \overline{y} =1\) now can be written as

\[
\label{twoDNSbc:}

\dfrac{ U_x (\overline{y} =1)}{U_{0}} = \overbrace{a\,h^2}^{\overline{a}}\,x^2 +

\overbrace{b \, xp\left( h \right) }^{\overline{b}}

\, xp\left( \overline{x}\right)

= \overline{a}\,\overline{x}^2 + \overline{b} xp{\overline{x}}

\]
Where \(\overline{a}\) and \(\overline{b}\) are the transformed coefficients in the dimensionless presentation.

After the boundary conditions the initial condition can undergo the non–dimensional process. The initial condition (10) utilizing the previous definitions transformed into

\[
\label{dim:eq:twoDNSicDless}

\dfrac{U_x(\overline{x}=0)}{U_{0x}} = \dfrac{U_{0y}}{U_{0x}}

f(\overline{y})

\]

\[
\label{dim:eq:twoDNSdefDevT}

\dfrac{\partial U_x}{\partial t} =

\dfrac{\partial \overbrace{\overline{U_x}}^{\dfrac{U_x}{U_{0x}} } U_{0x}}

{\partial \underbrace{\overline{t}}_{ \dfrac{t\,U_{0x}}{ h } }\dfrac{h}{U_{0x}} }

=

\dfrac

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[15]/span[1], line 1, column 5
```

\dfrac{\partial \overbrace{\overline{U_x} }^{\dfrac{U_x}{U_{0x}}} U_{0x}}

{\partial \underbrace{\overline{x} }_{\dfrac{x}{h}} h }

=

\dfrac

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[15]/span[2], line 1, column 5
```

\] The second derivative of velocity looks like

\[
\label{dim:eq:twoDNSdef2Dir}

\dfrac{\partial^2 U_x}{\partial x^2} =

\dfrac{\partial}{\partial \left( \overline{x} h \right) }

\dfrac{\partial \left( \overline{U_x} U_{0x} \right) }{\partial \left( \overline{x} h \right) } =

\dfrac{U_{0x}}{h^2}

\dfrac{\partial^2 \overline{U_x} }{\partial \overline{x}^2}

\]
The last term is the gravity \(g\) which is left for the later stage. Substituting all terms and dividing by density, \(\rho\) result in

\[
\label{dim:eq:twoDNSxx}

\begin{array} {rll}

\dfrac

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[17]/span[1], line 1, column 5
```

\overline{U_y} \dfrac{\partial \overline{U_x}}{\partial \overline{y}} +

\left. \overline{U_z} \dfrac{\partial \overline{U_x}}{\partial \overline{z}}\right) = &\\

&- \dfrac{P_0-P_\infty}{h\,\rho}

\dfrac{\partial \overline{P}}{\partial \overline{x}} +

\dfrac

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[17]/span[2], line 1, column 5
```

\end{array}

\] Dividing equation (27) by \({U_{0x}}^2/ {h}\) yields

\[
\label{dim:eq:twoDNSx1i}

\begin{array} {rll}

\left(\dfrac{\partial \overline{U_x} }{\partial \overline{t} } +

\right. & \overline{U_x}\dfrac{\partial \overline{U_x} }{\partial \overline{x}} +

\overline{U_y} \dfrac{\partial \overline{U_x}}{\partial \overline{y}} +

\left. \overline{U_z} \dfrac{\partial \overline{U_x}}{\partial \overline{z}}\right) = &\\

&- \dfrac{P_0-P_\infty}

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[18]/span[1], line 1, column 5
```

\dfrac{\mu}{U_{0x}\,h\,\rho} \left(\dfrac{\partial^2 U_x}{\partial x^2} +

\dfrac{\partial^2 U_x}{\partial y^2} + \dfrac{\partial^2 U_x}{\partial z^2}\right) +

\dfrac{ g \,h }

```
Callstack:
at (Bookshelves/Civil_Engineering/Book:_Fluid_Mechanics_(Bar-Meir)/09:_Dimensional_Analysis/9.3:_Nusselt's_Technique), /content/body/p[18]/span[2], line 1, column 5
```

\] The third kind of boundary condition is the mix condition. This category includes combination of the function with its derivative. For example a typical heat balance at liquid solid interface reads

\[
\label{dim:eq:mixBCex}

h(T_0 - T) = - k \dfrac{\partial T}{ \partial x}

\]
This kind of boundary condition, since derivative of constant is zero, translated to

\[
\label{dim:eq:mixBCexDim1}

h\, \cancel{(T_0 -T_{max})}\, \left( \dfrac{T_0 - T }{ T_0 -T_{max}} \right) =

- \dfrac{ k \,\cancel{\left( T_0 -T_{max} \right)} }{ \ell} \,

\dfrac{ - \partial \left( \dfrac{ T - T_0} { T_0 -T_{max} } \right) }

{ \partial \left(\dfrac{x}{\ell} \right) }

\]
or

\[
\label{dim:eq:mixBCexDim}

\left( \dfrac{T_0 - T }{ T_0 -T_{max}} \right) =

\dfrac{ k }{h\, \ell}

\dfrac{\partial \left( \dfrac{ T - T_0} { T_0 -T_{max} } \right) }

{ \partial \left(\dfrac{x}{\ell} \right) }

\Longrightarrow

\Theta = \dfrac{1}{Nu} \dfrac{\partial \Theta}{\partial \overline{x}}

\]
temperature are defined as

\[
\label{dim:eq:NuDef}

Nu = \dfrac{h\,\ell}{k} && \Theta = \dfrac{ T - T_0} { T_0 -T_{max} }

\]
and \(T_{max}\) is the maximum or reference temperature of the system. The last category is dealing with some non–linear conditions of the function with its derivative. For example,

\[
\label{dim:eq:surfaceTensionBCi}

\Delta P \approx {\sigma}\left({ \dfrac{1}{r_1} + \dfrac{1}{r_2} } \right) =

\dfrac{\sigma}{ r_1}\, \dfrac{r_1+ r_2}{r_2}

\]
Where \(r_1\) and \(r_2\) are the typical principal radii of the free surface curvature, and, \(\sigma\), is the surface tension between the gas (or liquid) and the other phase. The surface geometry (or the radii) is determined by several factors which include the liquid movement instabilities etc chapters of the problem at hand. This boundary condition (38) can be rearranged to be

\[
\label{dim:eq:surfaceTensionBC}

\dfrac{\Delta P \, r_1}{\sigma} \approx \dfrac{r_1+ r_2}{r_2} \Longrightarrow

Av \approx \dfrac{r_1+ r_2}{r_2}

\]
The Avi number represents the geometrical characteristics combined with the material properties. The boundary condition (39) can be transferred into

\[
\label{dim:eq:surfaceTensionBCdim}

\dfrac{ \Delta P\,r_1}{\sigma} = Av

\]
Where \(\Delta P\) is the pressure difference between the two phases (normally between the liquid and gas phase). One of advantage of Nusselt's method is the object-Oriented nature which allows one to add additional dimensionless parameters for addition "degree of freedom.'' It is common assumption, to initially assume, that liquid is incompressible. If greater accuracy is needed than this assumption is removed. In that case, a new dimensionless parameters is introduced as the ratio of the density to a reference density as

\[
\label{dim:eq:refereceRho}

\overline{\rho} = \dfrac{\rho}{\rho_0}

\]
In case of ideal gas model with isentropic flow this assumption becomes

\[
\label{dim:eq:gasDensity}

\bar{\rho} = {\rho \over \rho_0} = \left(\dfrac{ P_{0}}{ P } \right)^{\dfrac{1}{n}}

\]
The power \(n\) depends on the gas properties.

## Characteristics Values

Normally, the characteristics values are determined by physical values e.g. The diameter of cylinder as a typical length. There are several situations where the characteristic length, velocity, for example, are determined by the physical properties of the fluid(s). The characteristic velocity can determined from \(U_0 =\sqrt{2P_{0} / \rho}\). The characteristic length can be determined from ratio of \(\ell = \Delta P/\sigma\).

Example 9.12

One idea of renewable energy is to use and to utilize the high concentration of brine water such as in the Salt Lake and the Salt Sea (in Israel). This process requires analysis the mass transfer process. The governing equation is non-linear and this example provides opportunity to study nondimensionalizing of this kind of equation. The conversion of the species yields a governing nonlinear equation (??) for such process is

\[
\label{highMass:gov}

U_0 \dfrac{\partial C_A}{\partial x} =

\dfrac{\partial }{ \partial y }

\dfrac {D_{AB}}{ \left( 1 - X_A\right) } \dfrac{\partial C_A}{ \partial y }

\]

\[
\label{highMass:BCy}

\dfrac{\partial C_A}{\partial y} \left(y=\infty\right) = 0

\]

\[
\label{highMass:BCy2i}

C_A (y=0) = C_e

\]
Where \(C_e\) is the equilibrium concentration. The initial condition is

\[
\label{highMass:BCy2}

C_A (x=0) = C_0

\]
Select dimensionless parameters so that the governing equation and boundary and initial condition can be presented in a dimensionless form. There is no need to discuss the physical significance of the problem.

Solution 9.12

This governing equation requires to work with dimension associated with mass transfer and chemical reactions, the "mole.'' However, the units should not cause confusion or fear since it appear on both sides of the governing equation. Hence, this unit will be canceled. Now the units are compared to make sure that diffusion coefficient is kept the units on both sides the same. From units point of view, equation (43) can be written (when the concentration is simply ignored) as

\[
\label{highMass:govDim}

\overbrace{\dfrac{L}{t}} ^{U}

\overbrace{\dfrac{\cancel{C}}{L} }^{\dfrac{\partial C}{\partial x}} =

\overbrace{\dfrac{1}{L}}^{\dfrac{\partial }{\partial y}}

\overbrace{\dfrac{D_{AB}}{1}}^{\dfrac{D_{AB} }{ \left(1 - X \right)}}

\overbrace{\dfrac{\cancel{C}}{L} }^{\dfrac{\partial C}{\partial y}}

\]

\[
\label{highMass:govGov}

\dfrac{1}{t} = \dfrac{1}{L^2} D_{AB} \Longrightarrow D_{AB} = \dfrac{L^2}{t}

\]
Hence the units of diffusion coefficient are typically given by \(\left[m^2/sec\right]\) (it also can be observed that based on Fick's laws of diffusion it has the same units). The potential of possibilities of dimensionless parameter is large. Typically, dimensionless parameters are presented as ratio of two quantities. In addition to that, in heat and mass transfer (also in pressure driven flow etc.) the relative or reference to certain point has to accounted for. The boundary and initial conditions here provides the potential of the "driving force'' for the mass flow or mass transfer. Hence, the potential definition is

\[
\label{highMass:phi}

\Phi = \dfrac{C_A - C_0 }{C_e - C_0}

\]
With almost "standard'' transformation

\[
\label{highMass:coordinates}

\overline{x} = \dfrac{x}{\ell} &\qquad & \overline{y} = \dfrac{y}{\ell}

\]
Hence the derivative of \(\Phi\) with respect to time is

\[
\label{highMass:PhiD}

\dfrac{\partial \Phi}{\partial \overline{x}} =

\dfrac{\partial \dfrac{C_A - C_0 }{C_e - C_0} }{\partial \dfrac{x}{\ell} }

= \dfrac{\ell}{C_e - C_0}

\dfrac{\partial \left(C_A - \cancelto{0}{C_0} \,\,\,\,\,\right)} {\partial {x} }

= \dfrac{\ell}{C_e - C_0}

\dfrac{\partial C_A } {\partial {x} }

\]

In general a derivative with respect to \(\overline{x}\) or \(\overline{y}\) leave yields multiplication of \(\ell\). Hence, equation (43) transformed into

\[
\label{highMass:totalDim}

\begin{array}{rl}

U_0\dfrac{\cancel{\left(C_e - C_0\right)}} {\ell}\dfrac{\partial \Phi}{\partial \overline{x}} &=

\dfrac{1}{\ell}\, \dfrac{\partial }{ \partial \overline{y} }

\dfrac {D_{AB}}{ \left( 1 - X_A\right) }

\dfrac{\cancel{\left(C_e - C_0\right)}} {\ell}\dfrac{\partial \Phi}{\partial \overline{y}} \

\mbox{\Huge $\displaystyle\leadsto$}

\dfrac{U_0} {\ell}\dfrac{\partial \Phi}{\partial \overline{x}} &=

\dfrac{1}{\ell^2}\, \dfrac{\partial }{ \partial \overline{y} }

\dfrac {D_{AB}}{ \left( 1 - X_A\right) }

\dfrac{\partial \Phi}{\partial \overline{y}}

\end{array}

\]
Equation (52) like non–dimensionalized and proper version. However, the term \(X_A\), while is dimensionless, is not proper. Yet, \(X_A\) is a function of \(\Phi\) because it contains \(C_A\). Hence, this term, \(X_A\) has to be converted or presented by \(\Phi\). Using the definition of \(X_A\) it can be written as

\[
\label{highMass:X}

X_A = \dfrac{C_A} {C} = \left( C_e- C_0 \right)\dfrac{C_A - C_0}{ C_e- C_0} \dfrac{1}{C}

\]
Thus the transformation in equation (53) another unexpected dimensionless parameter as

\[
\label{highMass:Xf}

X_A = \Phi \, \dfrac{C_e- C_0}{C}

\]
Thus number, \(\dfrac{C_e- C_0}{C}\) was not expected and it represent ratio of the driving force to the height of the concentration which was not possible to attend by Buckingham's method.

## Contributors and Attributions

Dr. Genick Bar-Meir. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or later or Potto license.