15 Dimensional Analysis

Introduction

Engineers quickly begin to recognize from even the most casual observations that understanding and solving real-world problems involves considering many parameters, as well as parameters with other dependencies and interdependencies. In other words, many quantities and factors could affect the physical behavior of a given system, and different effects will be obtained when one parameter is varied versus another. In particular, the behavior of air and other fluids is relatively complex to describe because they can be sheared and deformed, primarily as they flow around objects; notice that “objects” are often referred to as “bodies” in the field of fluid dynamics.

The fluid also produces forces and moments on such bodies, including the familiar lift and drag forces, as shown in the figure below. There is undoubtedly an intuitive expectation that what the flow does (e.g., where the fluid goes) and what effects it has on the body will be determined at least partially by the geometrical shape of the body, as well as its angular orientation to the flow. The forces can also be expected to depend on the fluid properties already discussed, such as flow velocity (airspeed), density, viscosity, temperature, pressure, speed of sound, etc. However, other parameters may also be involved.

A body in a flow will produce forces and moments on itself, which depend on several fluid parameters as well as its shape and orientation to the flow.

To solve this type of fluid dynamics problem, one approach would be to set down and use the equations of motion governing the behavior of the fluid and then try to solve these equations. These nonlinear partial differential equations are called the Navier-Stokes equations and form the basis of a class of methods called computational fluid dynamics (usually referred to as “CFD”). However, this is not a straightforward approach because of the complexity of the nonlinear partial differential equations that must be solved and because considerable proficiency is needed in the solution methods. Another approach would be to do experiments to find the relationships between the forces on a body and the fluid properties, e.g., in tests done in a wind tunnel. At least in principle, either or both of these latter approaches will allow the effects of the fluid on the body to be determined, as well as the reciprocal influences of the body on the fluid. However, even with CFD and the ability to do experiments, sequential variations in the values of the problem parameters to determine cause and effect will take a very long time.

Nevertheless, another method can be used to help determine the general effects of key physical quantities on aerodynamic behavior. This method is called dimensional analysis, which reduces the number of unknowns governing the understanding of the problem into sets of dimensionless or non-dimensional groupings called similarity parameters. These similarity parameters then allow a significant simplification in understanding the effects of the original problem parameters. Another advantage of dimensional analysis is the ability to help predict the behavior of a given system at a larger scale from results known at a smaller scale, an important principle used in many aspects of engineering design. Dimensional analysis is not restricted to fluid problems and is helpful in many fields of engineering and science.

Learning Objectives

  • Appreciate the principle of dimensional homogeneity and balancing the dimensions of an equation.
  • Understand how to perform dimensional analysis, a method used to bring out the primary dependencies in an engineering problem in terms of dimensionless groupings called similarity parameters.
  • Realize the significance of dimensionless similarity parameters and how they can simplify a problem and reduce the number of dependencies, i.e., making a problem more straightforward to comprehend.
  • Develop and discuss some exemplar problems to become familiar with the method and practice of dimensional analysis and the use of similarity parameters.

Physical Quantities

Recall that physical quantities are of two types: base quantities and derived quantities. Base quantities such as mass, length, time, and temperature can be readily measured. Derived quantities are things like velocity, pressure, energy, and power, all of which are derived from these base quantities. In dimensional analysis, only base quantities are used; those encountered in mechanics and engineering are listed in the table below. Other base quantities that may arise in problem-solving are electrical current (Amp), the intensity of light (candela), and concentration (mole). Notice also that unit symbols are generally not italicized. For example, the symbol for length is L, but the symbol for the base unit of length is L. Likewise, the symbol for temperature is T, but the symbol for the base unit of temperature is \theta.

List of key base quantities used in engineering mechanics.
Base quantity Symbol Base unit SI unit USC unit
Mass M M kilogram (kg) slug (or slugs or sl)
Length L L meter (m) foot (ft)
Time t T second (s) second (s)
Temperature T \theta Kelvin (^{\circ}K or K) Rankine (^{\circ}R or R)

Many of the parameters routinely used in engineering are derived parameters, in that their units are derived from combinations of base units, e.g., volume, density, force, work, energy, power, etc. Consider the following examples of derived parameters, noting that a parameter inside a set of square brackets [ ] means the dimensions of that parameter, i.e.,

(1)   \begin{equation*} \left[ \text{parameter} \right] \equiv \text{dimensions of that parameter} \end{equation*}

  • Velocity, V, is distance per unit time, i.e., L/T = L T^{-1}. Notice also that scientific notation must be used throughout dimensional analysis to avoid ambiguity in the final units of the derived parameter. Therefore, [V] = L T^{-1}.
  • Volume, {\cal{V}}, is the volume of something such as a cuboid, then its width (units of L) must be multiplied by its length (units of L) by its height (units of L). Therefore, the derived units for volume are L{^3} or [{\cal{V}}] = L{^3}.
  • The density, \varrho,  of some substance, is its mass (units of M) per unit volume (units of L{^3}), so density will have units of M L{^{-3}}, i.e., [\varrho] = M L{^{-3}}.
  • A force is equivalent to the product of a mass and its acceleration (Newton’s second law), so in terms of its base units, then [F] = M LT^{-2}.
  • Torque, given the symbol \tau or T, is a moment, so the product of a force times a distance or “arm.” Therefore, [\tau] \equiv [T]  = (M LT^{-2}) L = M L{^2} T^{-2}.

Obtaining the correct base units of any given parameter takes some practice, but it is essential for correctly using dimensional analysis. Again, the scientific notation should be used when expressing the units of a parameter to avoid ambiguity, e.g., always use the form M L T^{-2} and not M L / T^2.

Aerodynamic Force on a Wing

To introduce dimensional analysis methods, consider the forces produced on a wing placed in the test section of a wind tunnel, as shown in the photo below. The resultant aerodynamic force on the wing is R. For convenience, it can be assumed that this force acts at a center of pressure, which, by definition, is a point where there are no pitching moments. In 1687, Isaac Newton, in his Philosophiae Naturalis Principia Mathematica (Mathematical Principles of Natural Philosophy) made three statements say that the forces acting on two geometrically similar bodies which move in fluids with different densities are proportional to:

  1. The square of the velocity of the fluid.
  2. The square of the linear dimension of the body.
  3. The density of the fluid.

He further notes that the proportionality of the force to the square of the velocity is restricted to lower flow speeds because ballistic experiments, even in Newton’s time, suggested that these statements do not apply to flow velocities greater than the speed of sound. Under these conditions, the thermodynamics of the fluid flow will also become important.

A wing set at an angle to the flow will produce some resultant aerodynamic force, the magnitude and direction of which depends on the size of the wing and the free-stream parameters.

In general, for a subsonic flow, it would be expected that the resultant aerodynamic force R on the wing for a given orientation to the flow will depend on:

  1. The free-stream velocity, V_{\infty}.
  2. The free-stream density, \varrho_{\infty}.
  3. The size of the wing, as described by its planform area S, which for a rectangular wing is the product of its chord c and span b, i.e., S = b \, c, as shown in the figure below.
  4. The sonic velocity a_{\infty}, i.e., the potential for compressibility effects.
  5. The viscosity of the air, i.e., on the coefficient of viscosity \mu_{\infty}.
A wing set at an angle to the flow will produce an aerodynamic force, the magnitude and direction of which depends on the size of the wing and the free-stream parameters.

The subscript \infty or “infinity” means undisturbed conditions, i.e., the “free-stream” or reference values of these quantities that are constant well upstream of the wing. Defining fixed reference values is logical because otherwise, all of these values will change as the flow approaches the wing. Notice that the angle of attack of the wing, \alpha, is the angle between the chordline of the wing (i.e., a line running from the leading edge or nose to the trailing edge or tail of the wing) and the direction of the free-stream velocity V_{\infty}.

The relationship between R and the air properties may now be written in a general functional form as

(2)   \begin{equation*} R = \phi \left( \varrho_{\infty}, \, V_{\infty}, \, S, \, \mu_{\infty}, \, a_{\infty} \right) \end{equation*}

where R is referred to as the dependent variable and \varrho_{\infty}, V_{\infty}, \mu_{\infty} and a_{\infty} are called the independent variables. The functional dependence, i.e., the relationship(s) between the dependent and independent variables, denoted by \phi in this case, is still to be determined, which, more often than not, is determined from experiments.

Although this latter equation expresses a general relationship, its use in this form is mostly impractical for the direct measurement or calculation of R. To see why, consider how R could be measured in the wind tunnel by a student of aerodynamics. The procedure might start with the student placing a wing at a given angle of attack, \alpha, and then systematically varying each of the dependencies V_{\infty}, \varrho_{\infty}, \mu_{\infty} and a_{\infty} (i.e., a total of five parameters) in turn, and in each case measuring the corresponding values of R.

Varying the wind speed V_{\infty} can be controlled easily in the wind tunnel, and different wings could be tested to change the wing area S. The other parameters, however, are more challenging to vary directly or systematically. However, if it could be done, the total number of possible combinations of parameters required in the wind tunnel would be 5^5 = 3,125 (!). Then, that process would need to be repeated by the student for different values of the wing origination to the flow, \alpha. If ten angles were selected, over 30,000 runs would be needed. Eventually, the relationships required between R and all dependent quantities at each angle of attack would be determined.

Therefore, it takes little to figure out that this process of systematically varying each dependent parameter, in turn, and holding the other parameters constant could take the student a very long time in the wind tunnel. For example, if 5 parameters have to be varied sequentially, the total number of test cases can be determined by knowing the number of different values each parameter can take. Denote the number of values each parameter can take as follows: v_1 for the first parameter, v_2 for the second parameter, v_3 for the third parameter, v_4 for the fourth parameter, and v_5 for the fifth parameter. Then the total number of test cases is the product of the number of values each parameter can take, is

(3)   \begin{equation*} \text{Total test cases} = v_1 \times v_2 \times v_3 \times v_4 \times v_5 \end{equation*}

For example, if each parameter can take 5 different values, the total number of test cases would be

(4)   \begin{equation*} 5 \times 5 \times 5 \times 5 \times 5 = 5^5 = 3,125 \end{equation*}

Most would agree this is a lot of tests in the wind tunnel or runs on the computer!

Furthermore, it would likely be challenging to keep the other parameters constant while the effects of one parameter are being measured. This approach would also result in a large quantity of data values, from which cross-plotting and perhaps curve-fitting could be used to find the needed functional relationships as to how V_{\infty}, \varrho_{\infty}, etc., would each affect R. If just one dependent variable were used, then two columns of numbers would be obtained, and the functional relationship would be easy to establish. However, a book of tables would be needed for four or more dependent variables or the equivalent of many pages of charts and graphs. Hence, it soon becomes apparent to the student that managing an engineering problem with three, four, or more independent variables is exceptionally challenging.

Basis of Dimensional Analysis

It is clear that systematically varying each independent parameter to establish the “cause and effect” on the dependent parameter would be impractical from a measurement perspective, nor would it be attractive from a computational (modeling) approach. Fortunately, the approach to the preceding problem can be simplified by employing dimensional analysis, which can be used to determine sets of dimensionless groups of parameters that affect the problem of interest. Such an approach will also be shown to reduce the number of independent variables, thereby becoming an essential technique for understanding many engineering problems, not just aerodynamics.

Dimensional Homogeneity

Dimensional analysis is based on the fact that if an equation of parameters such as

(5)   \begin{equation*} P_1 + P_2 + P_3 = P_4 \end{equation*}

expresses some physical relation, then the parameters P_1, P_2, P_3, and P_4 must all have the same dimensions and so be dimensionally homogeneous for the equation to be mathematically valid, e.g., P_1, P_2, P_3, and P_4 could all be values of pressures or velocity components. Equation 5 can also be made dimensionless by dividing through by one of the terms, i.e., the same result can be expressed as

(6)   \begin{equation*} \frac{P_1}{P_4} + \frac{P_2}{P_4} + \frac{P_3}{P_4} = \frac{P_4}{P_4} = 1 \end{equation*}

Another advantage of this latter dimensionless form is that all quantities are now expressed as ratios, so the system of units is irrelevant.

Complete and Incomplete Equations

An equation arising from a physical problem is said to be complete if stated so that the physical dimensions of each variable involved can be explicitly recognized. Consider a falling body accelerating to some velocity, V, under gravity, g, from an initial velocity, U, then

(7)   \begin{equation*} V^2 = U^2 + 2 \, g \, s \end{equation*}

where s is vertical distance, then the form of this equation is said to be complete. However, if the equation is written as

(8)   \begin{equation*} V^2 = U^2 + 19.62 s \end{equation*}

then it is incomplete because the dimensions of the second term in the equation cannot be recognized, i.e., it is impossible to positively conclude that 19.62 equals 2g. However, that would certainly be a good guess with SI units (g = 9.81 m/s2).

Recall that a parameter inside a set of square brackets [ ] means the dimensions of that parameter. If the equation is complete, then it is possible to conclude that in terms of dimensions, then

(9)   \begin{equation*} \left[ V^2  \right] = \left[ U^2 \right] = \left[ 2 \, g \, s \right] \end{equation*}

Inserting the base units for each parameter gives

(10)   \begin{equation*} \left( \frac{\rm L}{ \rm T} \right)^2  = \left( \frac{ \rm L}{ \rm T} \right)^2 = \left( \frac{ \rm L}{ \rm T^2} \right) \rm L  = \left( \frac{ \rm L}{ \rm T} \right)^2 = \left( \frac{ \rm L}{ \rm T} \right)^2 \end{equation*}

and the equation is, therefore, dimensionally homogeneous.

Furthermore, to make sense of any equation, each parameter must be expressed numerically in the same units system, i.e., either entirely in SI units or entirely in USC units, in each case with the same base units. However, notice that Eq. 7 can also be written in non-dimensional form as

(11)   \begin{equation*} \left( \frac{V}{U} \right)^2 - \frac{ 2 \, g \, s}{U^2} = 1 \end{equation*}

so that the system of units (SI or USC) does not matter. When undertaking derivations and developing equations, it is often desirable to work in non-dimensional forms.

Satisfying Dimensional Homogeneity

Consider the kinetic energy, KE, of a mass, M, traveling at velocity, V, which is given by

(12)   \begin{equation*} KE = \frac{1}{2} m V^2 \end{equation*}

Dimensional homogeneity is required in this equation, i.e., it is dimensionally balanced. On the left-hand side, the units are energy, which is the ability to do work so

(13)   \begin{equation*} \left[ KE \right] \equiv \mbox{Force times distance} = (\rm M L T^{-2} ) L = \rm M L^2 T^{-2} \end{equation*}

On the right-hand side, then

(14)   \begin{equation*} \left[ \frac{1}{2} m V^2\right] \equiv \mbox{Mass times velocity squared} = \rm M (L T^{-1})^2 = M L^2 T^{-2} \end{equation*}

Therefore, the dimensions on the left side of this equation balance those on the right side, so the equation is dimensionally homogeneous.

Consider another problem in which dimensional homogeneity must be satisfied. To describe the laminar flow of a fluid between two moving plates, the shear stress, \tau, is some function of the viscosity, \mu, and the velocity difference \Delta V between the two plates separated by a distance h. The relationship can be expressed as

(15)   \begin{equation*} \tau = \mu^{\alpha} \, \Delta V^{\beta} \, h^{\gamma} \end{equation*}

By satisfying dimensional homogeneity, it is possible to determine the values of \alpha, \beta, and \gamma. In terms of the dimensions, then

(16)   \begin{equation*} \left[ \tau \right] = \left[ \mu^{\alpha} \right] \left[ \Delta V^{\beta} \right] \left[ h^{\gamma} \right] \end{equation*}

and so for each parameter, then

(17)   \begin{equation*} \left[ \tau \right] = \rm M \rm L^{-1} \rm T^{-2} \quad\quad \left[ \mu \right] = \rm M \rm  L^{-1} \rm T^{-1} \quad\quad \left[ \Delta V \right] = \rm L \rm T^{-1} \quad\quad \left[ h \right] = \rm  L \end{equation*}

Therefore,

(18)   \begin{equation*} \rm M \rm L^{-1} \rm T^{-2}   = \left( \rm M \rm L^{-1} \rm T^{-1} \right)^{\alpha} \, \left( \rm L \rm T^{-1} \right)^\beta \, \left( \rm L \right)^\gamma \end{equation*}

To obtain dimensional homogeneity, then

(19)   \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 1 & = & \alpha \\ \text{\small Exponents~of~L:} \,\,  -1 & = & -\alpha + \beta + \gamma\\ \text{\small Exponents~of~T:} \,\, -2 & = & -\alpha -\beta \end{eqnarray*}

so \alpha = 1, \beta = 1, and \gamma = -1. Inserting the values gives the final relationship as

(20)   \begin{equation*} \tau = \mu^{\alpha} \, \Delta V^{\beta} \, h^{\gamma} = K \, \mu \left( \frac{\Delta V}{h} \right) \end{equation*}

showing that the shear stress is proportional to the velocity gradient, with the constant of proportionality being the value of viscosity, \mu. This is called Newton’s law of viscosity.

The Buckingham Π Method

These ideas of dimensional analysis and dimensional homogeneity are fundamental and intrinsic to the Buckingham \Pi method. Lord Rayleigh and Joseph Bertrand initially set down these preceding simple, but profoundly important, ideas in what is called the “Method of Dimensions.” However, they are now formally embodied in what is known as the “Buckingham \Pi theorem” after Edgar Buckingham, who developed it into a more structured form suitable for general use. The Buckingham \Pi Theorem (usually just called the Buckingham \Pi method or just the \Pi method) is a generalized formalization of the method of dimensions, which can now be examined. The original publication by Edgar Buckingham, a scientist at NIST, is available, which is worth reading in conjunction with the following exposition of his method.

Let P_{1}, \, P_{2}, .....P_{N}, represent a list of N physical variables in the general physical relationship

(21)   \begin{equation*} P_{1} = \phi \left( P_{2}, .........P_{N} \right) \end{equation*}

where P_{1} is the dependent parameter and P_{2}, .........P_{N} are the independent parameters. The function \phi represents a general functional relationship, as yet unknown. Notice that the previous equation can also be written in an implicit form as

(22)   \begin{equation*} \psi \left( P_{1}, \, P_{2}, .........P_{N} \right) = 0 \end{equation*}

where \psi represents another functional relationship between the parameters yet to be determined.

Let K equal the number of fundamental dimensions required to describe the physical variables. In engineering, most physical variables can be expressed as a combination of mass, length, time, and temperature, so generally, K=4. In mechanics, then mass, length, and time are involved, so K=3.  The Buckingham \Pi method then states that the relationship in Eq. 22 may be expressed as a relationship of fewer (N-K) dimensionless products called \Pi products, where each \Pi product is a dimensionless product of a selected set of K repeating variables plus one other variable. Therefore, one key to applying the \Pi method is selecting the most suitable set of repeating variables.

To see how the \Pi method works, let P_{2}, \, P_{3}, ....P_{K+1} be the selected K repeating variables. The repeating variables are not unique in most cases, but their selection is not arbitrary and must be selected according to some specific rules. The primary and most important rule is that the selected repeating variables cannot include the dependent variable, P_1, which can only appear in one of the \Pi groups. Therefore,

(23)   \begin{eqnarray*} \Pi_{1} & = & \phi_{1} \left( \underbrace{ P_{2}, \, P_{3}......P_{K+1}}_{\mbox{~\scriptsize Repeating}} ,\underbrace{P_{1}}_{\mbox{~\scriptsize Dependent}} \right) \\[6pt] \Pi_{2} & = & \phi_{2} \left( \underbrace{ P_{2}, \, P_{3}, .....P_{K+1}}_{\mbox{~\scriptsize Repeating}} ,\underbrace{P_{K+2}}_{\mbox{~\scriptsize Other}} \right) \\[6pt] & & \bigcdot  \\[6pt] & & \bigcdot  \\[6pt] & & \bigcdot  \\[6pt] \Pi_{N-K} & = & \phi_{N-K} \left( \underbrace{ P_{2}, \, P_{3}, .......P_{K+1}}_{\mbox{~\scriptsize Repeating}}, \underbrace{P_{N}}_{\mbox{~\scriptsize Remaining}} \right) \end{eqnarray*}

where P_N is the last remaining variable.

Notice that there are now N - K new functional relationships, i.e., \phi_{1}, \phi_{2}….\phi_{N-K}, but they are expressed in terms of the dimensionless \Pi products. Therefore, it can be written that

(24)   \begin{equation*} \phi_1 \left( \Pi_1, \, \Pi_2, .., \Pi_{N-K} \right) = 0 \end{equation*}

where \phi_1 is the functional relationship (still to be determined), it is now expressed in terms of each of the \Pi products. The functional relationship is often found from experiments, as suggested in the figure below. The process of dimensional analysis does not allow theoretical relationships to be developed, i.e., equations. Notice that in explicit form, this latter equation can also be written as

(25)   \begin{equation*} \Pi_1 = \phi_2 \left( \Pi_2, .., \Pi_{N-K} \right) \end{equation*}

where \phi_2 is some other functional relationship.

 

The functional relationships are usually determined using experiments, but computational methods can also be used. The functional relationships will be formed from discrete data.

The upshot of the dimensional analysis process means that the functional relationship has been reduced from the initial N dimensional variables to N-K dimensionless variables. Therefore, decreasing the number of dependent parameters significantly simplifies analyzing the original problem. For example, suppose originally N = 6 and K = 3, i.e., the dependent variable is a function of five other parameters as in the aforementioned wing problem. In that case, the dimensional analysis will reduce the problem to three dimensionless groupings, i.e., one dimensionless grouping becomes a function of two other dimensionless groupings. If N = 5 and K = 3, then the dimensional analysis will reduce the problem to two dimensionless groupings in that one grouping becomes a sole function of the other grouping. The simplifications, therefore, afforded by dimensional analysis become clear.

Π Method Applied to the Drag on a Sphere

Consider a sphere to be tested in a low-speed wind tunnel to measure its drag, as shown in the photograph below. The drag of a sphere sounds relatively uninteresting compared to a wing, but this problem helps introduce the dimensional analysis process without making the process too lengthy. Furthermore, the aerodynamic characteristics of a sphere are particularly relevant to all sorts of ball games, such as tennis, baseball, soccer, and many others. The motion of balls through the air is a complex fluid dynamic problem that is influenced by many factors, including the flow density and viscosity, the velocity and size of the ball, and its surface finish. For example, the aerodynamics of a golf ball affects the distance, spin, and overall trajectory of the ball. Ever wonder why a golf ball has dimples on its surface? In baseball, the aerodynamics of the ball can influence the speed and trajectory of a pitch.

A sphere set up in a wind tunnel to measure its drag. The sphere is mounted on a force balance.

Functional Dependency

The objective is to find the dimensionless groupings describing this aerodynamic problem. On a physical intuitive basis, the aerodynamic drag force, D, on the sphere would be expected to depend on:

  1. Free-stream velocity, V_{\infty}.
  2. Free-stream density, \varrho_{\infty}.
  3. Free-stream viscosity, i.e., the coefficient of viscosity \mu_{\infty}.
  4. The size of the sphere, as described by its diameter d.

This list may be incomplete, and other parameters may affect the drag. On the one hand, omitting a relevant parameter may restrict the validity of the result or make it wrong. On the other hand, including unnecessary parameters may make the final result more complicated than it should be.

The relationship between D and the air properties may be written in a general functional form as

    \[ D = \phi\left( \varrho_{\infty}, V_{\infty}, \mu_{\infty}, d \right) \]

The drag D is the dependent variable. \varrho_{\infty}, V_{\infty}, \mu_{\infty} and d are the independent variables. In this case, the functional dependence is denoted by \phi.

The functional dependence in implicit form is

    \[ \psi \left( \varrho_{\infty}, \, V_{\infty}, \, \mu_{\infty}, \, d,  \, D \right) = 0 \]

where \psi is some other function. In this case, N = 5, K = 3, N - K = 2, so there are two \Pi products to be determined.

Understanding explicit versus implicit functional forms

What is the difference between writing an equation y = f(x) and the form g(x,y) = 0? In mathematics, the dependent function, say y, for example, can be written in terms of other parameters, e.g., y = f(x) where {\scriptstyle{f}} is some function, or more generally as y = f(x_1, x_2, \! ... x_N), which are called explicit functional dependencies. However, the relationship can also be written in implicit form, i.e., y = f(x) is equivalent to g (x, y) = 0. Notice that the function {\scriptstyle{g}} is NOT the same as the function {\scriptstyle{f}}. If y = x^2 +5, for example, then y = f(x), which is an explicit dependency of {\scriptstyle{y}} on {{x}}. Alternatively, it can be written that y - x^2 - 5 = 0, which is in the form g(x,y) = 0, which is the implicit form of the same relationship. Likewise, y = f(x_1, x_2, \! ...  x_N) as g(x_1, x_2, \! ...  x_N, y) = 0. In dimensional analysis, the process classically starts by writing the functional dependency of the problem parameters in implicit form.

Dimensions of the Variables

Now, the dimensions of the variables in this problem must be determined. Remember that the square brackets [-] around a variable means “the dimensions of” that variable. For each dependency, then

    \begin{eqnarray*} \left[ D \right]  & =  & \rm M L T^{-2} \\ \left[ \varrho_{\infty} \right]  & = & \rm ML^{-3} \\ \left[ V_{\infty} \right] & = & \rm LT^{-1} \\ \left[ \mu_{\infty} \right] & = & \rm ML^{-1} \rm T^{-1} \\ \left[ d \right] & = & \rm L \end{eqnarray*}

Recall that the units of force are obtained from mass times acceleration, i.e., MLT^{-2}.

Dimensional Matrix

The solution to the dimensional analysis problem proceeds by setting up a dimensional matrix with all the parameters in their base units. Setting up the dimensional matrix, in this case, gives

    \[ \begin{array}{l|r|r|r|r|r|r} \mbox{\small Base~unit} &   D  & \varrho_{\infty} & V_{\infty} & d &  \mu_{\infty} \\ \hline \text{\small Mass, M}  \quad & 1  & 1 & 0 & 0 & 1 \\ \mbox{\small Length, L}  \quad &  1  & -3 & 1 & 1 & -1 \\ \mbox{\small Time, T}  \quad &-2 & 0 & -1 & 0 & -1 \end{array} \]

One purpose of the dimensional matrix is to lay the problem out clearly and test for the linear independence of the repeating variables in terms of the chosen primary dimensions.

Repeating Variables

Although the selection of the repeating variables can be viewed as somewhat arbitrary, in general, they must always be selected based on four essential rules:

  1. The repeating variables must all have a primary dependency on the dependent variable. In this case, the variables \varrho_{\infty}, V_{\infty} and d can all be expected to have strong effects on the dependent variable, D.
  2. The dependent variable, D, cannot be used as a repeating variable, so it must appear only in one group by default. No exceptions.
  3. The repeating variables must include all of the physical base dimensions of the problem, i.e., they must collectively include mass and/or length and/or time, as appropriate, for the particular problem in question. In this case, \varrho_{\infty}, V_{\infty}, and d collectively include the base dimensions of mass, length, and time.
  4. The repeating variables must be linearly independent of one another, i.e., repeating variables must not be chosen that will have the same units or products of the same units, e.g., repeating variables with two different length scales or two different reference velocities cannot be used. In this case, \varrho_{\infty}, V_{\infty}, and d are all linearly independent of each other.

The dimensional analysis process will likely succeed if these prior rules are followed. If not, then it will fail, often catastrophically. For complex problems, many of which are in engineering, the choice of repeating variables may not be obvious. However, they should still be selected as quantities most likely to have substantial effects on the dependent variable. For example, the velocity and density of the flow in an aerodynamic problem are excellent examples of two quantities that will significantly affect the aerodynamic forces.

Formal Check for Linear Independence

To formally check for linear independence (although, in most cases, this will be obvious), the determinant of the sub-matrix formed by the selected repeating variables must be determined, which in this case is

(26)   \begin{equation*} \begin{array}{l|r|r|r|} \mbox{\small Base~unit} & \varrho_{\infty} & V_{\infty} & d \\ \hline \text{\small Mass, M} \quad & 1 & 0 & 0 \\ \text{\small Length, L} \quad & -3 & 1 & 1 \\ \text{\small Time, T} \quad & 0 & -1 & 0 \end{array} \end{equation*}

If c_1, c_2, and c_3 are non-zero constants in the following equation for linear independence, then

(27)   \begin{equation*} c_1 \left[ \begin{array}{r} 1 \\ -3 \\ 0 \end{array} \right] + c_2 \left[ \begin{array}{r} 0 \\ 1 \\ -1 \end{array} \right] + c_3 \left[ \begin{array}{r} 0 \\ 1 \\ 0 \end{array} \right] \ne 0 \end{equation*}

Reassembling this in matrix form gives

(28)   \begin{equation*} \left[ \begin{array}{rrr} 1 & 0 & 0 \\ -3 & 1 & 1 \\ 0 & -1 & 0 \end{array} \right] \left[ \begin{array}{r} c_1 \\ c_2 \\ c_3 \end{array} \right] \ne 0 \end{equation*}

For this matrix equation to have a nontrivial solution, then the determinant must be non-zero, i.e., the test in this case is whether

(29)   \begin{equation*} \det \left[ \begin{array}{rrr} 1 & 0 & 0 \\ -3 & 1 & 1 \\ 0 & -1 & 0 \end{array} \right] \equiv \left| \begin{array}{rrr} 1 & 0 & 0 \\ -3 & 1 & 1 \\ 0 & -1 & 0 \end{array} \right| \ne 0 \, ? \end{equation*}

Therefore, the determinant is

    \[ 1 \bigg( (1 \times 0 - (-1 \times 1) \bigg) -  0 \bigg( (-3 \times 0) - (0 \times 1) \bigg) + 0 \bigg( (-3 \times -1 - (0 \times 1) \bigg) = 1 - 0 + 0 = 1 \]

Therefore, it has now been formally verified that the selected repeating variables \varrho_{\infty}, V_{\infty} and d are indeed linearly independent of each other. In most cases, however, the formal test is not needed because if the variables are chosen such that they do not have the same units or products of the same units, then they will be linearly independent.

Development of the Π Products

Following the Buckingham \Pi method, then the \Pi products in this case will be

    \begin{eqnarray*} \Pi_{1} & = & \phi_{1} \left( \underbrace{\varrho_{\infty}, \, V_{\infty}, \, d}_{\mbox{~\scriptsize Repeating}}  , \underbrace{D}_{\mbox{~\scriptsize Dependent}} \right)  \\[8pt] \Pi_{2} & = & \phi_{2} \left( \underbrace{\varrho_{\infty}, \, V_{\infty}, \, d}_{\mbox{~\scriptsize Repeating}}  , \underbrace{\mu_{\infty}}_{\mbox{~\scriptsize Remaining}} \right) \end{eqnarray*}

First Π Product

For \Pi_{1} then

    \[ \Pi_{1} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (d)^{\gamma} \, D \]

The values of the coefficients \alpha, \beta, and \gamma must now be obtained to make the equation dimensionally homogeneous and dimensionless.  In terms of the dimensions of the parameters, then

    \[ \left[ \Pi_{1} \right] = 1 = \rm M^{0}L^{0}T^{0} = \left( \rm ML^{-3} \right)^{\alpha} \left( LT^{-1} \right)^{\beta} \left( \rm L \right)^{\gamma} \left( \rm MLT^{-2} \right) \]

Because the left-hand side of the above equation is dimensionless, it means that the powers or exponents of M, L, and T must add to zero to make the equation dimensionally homogeneous and dimensionless. Therefore, for \Pi_{1} to be dimensionless, then

    \[ \rm M^{0}L^{0}T^{0} = \rm M^{(\alpha+1)} \ \rm L^{(-3\alpha +\beta + \gamma+1)} \ \rm T^{(-\beta-2)} \]

so that

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha+1  \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha +\beta+ \gamma+1  \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta-2 \end{eqnarray*}

By inspection, \alpha = -1, \beta = -2, and \gamma = -2. Therefore, the first \Pi product is

    \[ \Pi_{1} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (d)^{\gamma} D = \varrho_{\infty}^{-1} V_{\infty}^{-2} \, d^{-2} \, D \]

or

    \[ \Pi_{1} = \frac{D}{\varrho_{\infty} V_{\infty}^{2} \, d^2} \]

which is a force coefficient. Usually, aerodynamic force coefficients are defined in terms of the dynamic pressure, i.e., (1/2)\varrho_{\infty} V_{\infty}^2. Also, in this case, a good reference is to use projected area \pi \, d^2/4. Therefore, the force coefficient would be written as

    \[ \Pi_{1} = C_D = \frac{D}{ \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} (\pi \, d^2/4)}  = \frac{D}{ \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} A} \]

It is always good to check the final units of the grouping to be sure it is dimensionless, i.e., it has units of 1. So, in this case

    \[ \left[ C_D \right]  = \frac{\rm M L T^{-2}}{ (\rm M L^{-3}) (\rm L^2 T^{-2}) (L^{2})} = \frac{\rm M L T^{-2}}{\rm M L T^{-2}}  =  1 \]

which confirms that the grouping is dimensionless.

Dealing with exponents or powers

In dimensional analysis, solving for the dimensionless groupings involves adding the values of the unknown exponents for each base dimension. Exponents are also called the power of the quantity because they represent the number of times a quantity is multiplied by itself. Here are some examples:

    \[\small 2^2 \, 2^3 = 2^{2+3} = 2^5 \]

    \[ x^a \, x^b = x^{a+b} \]

    \[ x^a \, y^b \, y^c = x^a \, y^{b+c} \]

    \[ \rm M^{\alpha} \, \rm M^{\beta + 1} = \rm M^{(\alpha + \beta + 1)} \]

    \[ \rm M^{\alpha} \, \rm L^{\gamma} \, \rm M^{\beta -2 } \, L^{-1} = \rm M^{(\alpha + \beta - 2)} \, \rm L^{\gamma -1} \]

When working with exponents, remember that the base parameter has to be the same, and the addition is performed on the exponent.

Second Π Product

For \Pi_{2} then

    \[ \Pi_{2} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (d)^{\gamma} \, \mu_{\infty} \]

so, in terms of dimensions, then

    \[ \rm M^{0} L^{0} T^{0} = \left(\rm ML^{-3} \right)^{\alpha} \left( \rm LT^{-1} \right)^{\beta} \left(\rm L \right)^\gamma \left( \rm ML^{-1}T^{-1} \right) \]

Equating the powers gives

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha +1  \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha +\beta+ \gamma-1  \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta-1 \end{eqnarray*}

Therefore, in this case \alpha = -1, \beta = -1, and \gamma = -1.

    \[ \Pi_{2} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (S)^{\gamma} \mu_{\infty} = \varrho_{\infty}^{-1} V_{\infty}^{-1} \, d^{-1} \, \mu_{\infty} \]

i.e.,

    \[ \Pi_{2} = \frac{\mu_{\infty}}{\varrho_{\infty} V_{\infty} d} \]

Or by inverting the quantity, a dimensionless grouping is still obtained, i.e.,

    \[ \Pi_{2} = \frac{\varrho_{\infty} V_{\infty} d}{\mu_{\infty}} \]

Remember from before that this dimensionless parameter is called the Reynolds number. Again, a check should be made to be sure that the grouping is dimensionless, so in this case

    \[ \left[ Re \right]  = \frac{(\rm M L^{-3}) ( \rm L T^{-1}) (\rm L) }{ \rm M L^{-1} T^{-1}} = \frac{\rm M L^{-1} \rm T^{-1} }{ \rm M L^{-1} T^{-1}} =  1 \]

which indeed confirms that the grouping is dimensionless.

Final Result for Functional Dependency

As a result of the dimensional analysis of the sphere, then

    \[ \phi_1 \left( \frac {D}{\frac{1}{2} \varrho_{\infty} V_{\infty}^{2} A}, \frac{\varrho_{\infty} V_{\infty} d}{\mu_{\infty}} \right) = 0 \]

or

    \[ \phi_1 \left( C_{D}, Re \right) = 0 \]

In explicit form, the relationship between the drag coefficient and the Reynolds number is

    \[ C_{D} = \phi_2 \left( Re \right) \]

Benefits of Dimensional Analysis

The benefits of dimensional analysis have now become much clearer. In this latter problem of finding the effects of different parameters on the drag of a sphere, the question for the analyst comes down to finding the explicit relationship between drag and four parameters versus finding the relationship between a force coefficient and the Reynolds number, as summarized in the figure below. Most analysts would side with the latter approach, not only because it involves less work but because the results from other experiments and at other Reynolds numbers, such as in other wind tunnels, can be directly compared on the same non-dimensional basis.

The benefits of dimensional analysis and the use of dimensionless parameters soon become apparent when engineers embark on experiments.

Dimensional analysis is also essential for generalizing data and establishing causal relationships, including for design purposes. In general, dimensional analysis provides a systematic approach for consistency checking, unit conversion, estimation, modeling, experimental design, equation verification, and problem-solving across various engineering and scientific disciplines.

Doing the Experiment

The functional dependency to be measured is now C_{D} = \phi_2 \left( Re \right). The drag force on the sphere, D, can be measured with a balance and with measurements of the flow velocity V_{\infty}, as well as density, \varrho_{\infty}, and viscosity \mu_{\infty}, (via measurements of pressure and temperature), calculations are made of:

1. The drag coefficient, C_D  = \displaystyle{\frac {D}{\frac{1}{2} \varrho_{\infty} V_{\infty}^{2} A}}

2. The Reynolds number, Re = \displaystyle{\frac{\varrho_{\infty} V_{\infty} d}{\mu_{\infty}}}

The figure below shows the results for the sphere’s drag coefficient as measured in the wind tunnel as a function of Reynolds number. The data collection can proceed rapidly because the values of the groupings are being plotted rather than the values of the individual parameters. Notice the rapid decrease in the drag coefficient between Re = 3 \times 10^5 and Re = 4 \times 10^5, which suggests that something very profound is happening to the flow about the sphere. Indeed, this behavior is tied to the flow characteristics on the sphere’s surface, called the boundary layer, and how the flow separates behind the sphere and creates a wake. In this dimensionless form, these results can also be directly compared to measurements made, for example, in another wind tunnel or added to a pool of data covering other ranges of Reynolds numbers.

Drag coefficient of a sphere as a function of Reynolds number based on diameter. Notice the different relationships depending on the surface roughness.

Effects of surface roughness?

Notice also from the results given in the preceding plot that the sphere’s drag coefficient depends on whether the sphere has a smooth or rough surface. Indeed, this is a very interesting outcome as it affects the drag of a golf ball. Therefore, it is helpful to think about what might happen if surface roughness was added as the fifth variable to the parameter list, i.e., in this case, assume that

    \[ D = \phi \left( \varrho_{\infty}, V_{\infty}, \mu_{\infty}, d , \epsilon \right) \]

where \epsilon is a linear measure of the roughness height. The new functional dependence of D in the implicit form is

    \[ \psi\left( \varrho_{\infty}, \, V_{\infty}, \, \mu_{\infty}, \, d, \, \epsilon, \, D \right) = 0 \]

In this case, N = 6, K = 3, N - K = 3, so there are three \Pi products. The first two \Pi products have already been determined, and it is easily shown that now

    \[ C_{D} = \phi_3 \left( Re , \frac{\epsilon}{d} \right) \]

where \epsilon/d is a dimensionless roughness height. Notice that the effects of roughness decrease the drag of a sphere in the range where golf balls fly, which allows them to go further.

Π Method Applied to a Wing

Returning to the wing problem initially discussed, the explicit dependency of the resultant aerodynamic force R (the dependent quantity) in terms of the independent variables is

(30)   \begin{equation*} R = \phi \left( \varrho_{\infty}, V_{\infty}, S, \mu_{\infty}, a_{\infty} \right) \end{equation*}

The functional dependence can also be written in implicit form, i.e.,

(31)   \begin{equation*} \psi \left( \varrho_{\infty}, \, V_{\infty}, \, S, \, \mu_{\infty}, \, a_{\infty}, \, R \right) = 0 \end{equation*}

where \psi is some functional dependency to be determined. Following the \Pi method then N = 6, K = 3, N - K = 3, and so three \Pi products must be found in this case.

The resultant force on a wing of a given size and orientation to the flow will be dependent on several parameters, including the velocity, density, and viscosity of the flow. This wing has an aspect ratio of 4.

Dimensions of the Variables

Now, the dimensions of the variables must be determined. For each dependency in this case, then

    \begin{eqnarray*} \left[ R \right] & = & \rm  MLT^{-2} \\ \left[ \varrho_{\infty} \right] & = & \rm ML^{-3} \\ \left[ V_{\infty} \right] & = &\rm LT^{-1} \\ \left[ S \right] & = & \rm L^2 \\ \left[ \mu_{\infty} \right] & = & \rm ML^{-1}T^{-1} \\ \left[ a_{\infty} \right] & = & \rm LT^{-1} \end{eqnarray*}

Notice that the units of force are obtained from the product of mass and acceleration, i.e., [F] = ( M ) ( LT^{-2}). It is critically important that each of the variables is assigned the correct base units. If not, the Buckingham \Pi method will fail.

Forming the Dimensional Matrix

The solution to the dimensional analysis problem proceeds by setting up a dimensional matrix with all of the parameters being given in terms of their base units, i.e.,

(32)   \begin{equation*} \begin{array}{l|r|r|r|r|r|r} \mbox{\small Base~unit} & R & \varrho_{\infty} & V_{\infty} & S & \mu_{\infty} & \alpha_{\infty} \\ \hline \text{\small Mass, M}  \quad & 1 & 1 & 0 & 0 & 1 & 0 \\ \text{\small Length, L}  \quad & 1 & -3 & 1 & 2 & -1 & 1 \\ \text{\small Time, T} \quad & -2 & 0 & -1 & 0 & -1 & -1 \end{array} \end{equation*}

Choice of Repeating Variables

In this problem, a good choice of the repeating variables is \varrho_{\infty}, V_{\infty}, and S, which is a common choice for aerodynamic problems. Remember not to choose the dependent variable as a repeating variable. Remember also that the repeating variables must be linearly independent, i.e., the units of any one repeating variable cannot be linearly obtained from the units of any other. Choosing repeating variables with different dimensions or products of dimensions is necessary. For example, two repeating variables comprising a length and an area or two repeating variables with units, cannot be used, or the method will fail. Furthermore, the repeating variables must collectively include all of the base dimensions of the problem. In this case, the choice of the repeating variables \varrho_{\infty}, V_{\infty}, and S satisfy all of these requirements.

Recap: Rules for choosing the repeating variables

The formal determination that the repeating variables are linearly independent takes some work, although it is just arithmetic. However, in some cases, this outcome can be verified by inspection. Some basic rules can help choose the repeating variables, so the formal test for linear independence can be foregone. These rules are:

  • NEVER choose the dependent variable as a repeating variable! If this rule is ignored, then the Buckingham \Pi method will fail catastrophically.
  • Always choose repeating variables that will, or are likely to, have a primary effect on the dependent variable. For aerodynamic problems, these are usually the density of the flow, the speed of the flow, and some area or length scale.
  • Do NOT choose repeating variables that have the same dimensions, so don’t choose two repeating variables, such as flow velocity and the speed of sound (which both have units of length/unit time), or two length variables, such as a wing chord and a wing span (both have units of length). Furthermore, two variables with units of length and area, respectively, cannot be used as repeating variables because they are also not linearly independent.
  • Remember that the choice of repeating variables may not, in most cases, be unique! The choice is not arbitrary, but a different choice of repeating variables will inevitably lead to a different set of dimensionless groupings. Even though dimensionless groupings are obtained, the process can be conducted again with a different set of repeating variables.

Development of Π Products

Continuing to follow the Buckingham \Pi procedure, then the \Pi products in this problem are

    \begin{eqnarray*} \Pi_{1} & = & \phi_{1} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~\scriptsize Repeating}} , \underbrace{R}_{\mbox{\scriptsize ~Dependent}} \right) \\[8pt] \Pi_{2} & = & \phi_{2} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~\scriptsize Repeating}} , \underbrace{\mu_{\infty}}_{\mbox{~\scriptsize Other}} \right) \\[8pt] \Pi_{3} & = & \phi_{3} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~Repeating}} , \underbrace{a_{\infty} }_{\mbox{~\scriptsize Remaining}}\right) \end{eqnarray*}

First Π Product

Take the first \Pi product, \Pi_{1}, then

(33)   \begin{equation*} \Pi_{1} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (S)^{\gamma} \, R \end{equation*}

where the coefficients \alpha, \beta, and \gamma are unknowns, and their values must be obtained appropriately to make the equation dimensionally homogeneous, i.e., finding their values to make the equation dimensionless with units of 1. In terms of the dimensions, then

(34)   \begin{equation*} \left[ \Pi_{1} \right] = \rm  M^{0}L^{0}T^{0} = \left( \rm ML^{-3} \right)^{\alpha} \left(\rm LT^{-1} \right)^{\beta} \left( \rm L^2 \right)^{\gamma} \left( MLT^{-2} \right) \end{equation*}

where \left[x\right] is used to mean “the dimensions of variable x.” Collecting the terms gives

(35)   \begin{equation*} \left[ \Pi_{1} \right] = \rm M^{0}L^{0}T^{0} = \rm M^{(\alpha+1)} \, \rm L^{(-3\alpha+\beta+2\gamma+1)} \, \rm T^{(-\beta-2)} \end{equation*}

Because the left-hand side of the above equation is dimensionless, it means that the powers or exponents of M, L, and T must add to zero to make the equation dimensionally homogeneous, i.e.,

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha+1 \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha+\beta+2\gamma+1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta-2 \end{eqnarray*}

Therefore, in this case by inspection then \alpha = -1, \beta = -2, and \gamma = -1 and so the first \Pi product is

(36)   \begin{equation*} \Pi_{1} = \varrho_{\infty}^{-1} V_{\infty}^{-2} S^{-1} R \end{equation*}

or

(37)   \begin{equation*} \Pi_{1} = \frac{R}{\varrho V_{\infty}^{2} S} \end{equation*}

Hence, \Pi_{1} is a form of force coefficient. Aerodynamic force coefficients are nearly always defined in terms of the dynamic pressure, i.e., (1/2)\varrho_{\infty} V_{\infty}^2, which is referred to as q_{\infty}, so more usually the force coefficient would be written as

(38)   \begin{equation*} \Pi_{1} = C_R = \frac{R}{ \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} S} = \frac{R}{ q_{\infty} S} \end{equation*}

Notice that multiplying a similarity parameter by a constant value does not change the units of that parameter.Paralleling the above approach used to find \Pi_{1}, the two remaining \Pi products, namely \Pi_{2} and \Pi_{3}, can now be derived.

Second Π Product

For \Pi_2 then

(39)   \begin{equation*} \Pi_{2} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (S)^{\gamma} \mu_{\infty} \end{equation*}

so

(40)   \begin{equation*} \left[ \Pi_{2} \right] = \rm M^{0} L^{0} T^{0} = \left(\rm ML^{-3} \right)^{\alpha} \left( \rm LT^{-1} \right)^{\beta} \left( \rm L^2 \right)^\gamma \left( \rm ML^{-1}T^{-1} \right) \end{equation*}

and

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha+1 \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha+\beta+2\gamma-1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta-1 \end{eqnarray*}

Therefore, in this case \alpha = -1, \beta = -1, and \gamma = -1/2 and so this \Pi product becomes

(41)   \begin{equation*} \Pi_{2} = \varrho_{\infty}^{-1} V_{\infty}^{-1} S^{-1/2} \mu_{\infty} \end{equation*}

This derived \Pi product is a dimensionless parameter that will become infinite as the velocity tends to zero. This behavior is inconvenient, but the \Pi product can be inverted to prevent that. The outcomes are still a dimensionless parameter, but now one that is more convenient, i.e., the \Pi product becomes

(42)   \begin{equation*} \Pi_{2} = \frac{\varrho_{\infty} V_{\infty} S^{1/2}}{\mu_{\infty}} \end{equation*}

Because, in this case, there is a square root of an area, it can be replaced by a linear dimension, such as the chord of the wing c, i.e., the area of a rectangular wing will be its chord c times the span of the wing b or S = c \, b, so now

(43)   \begin{equation*} \Pi_{2} = \frac{\varrho_{\infty} V_{\infty} \, c}{\mu_{\infty}} \end{equation*}

The choice of the appropriate length scale is somewhat arbitrary (e.g., the wing span could have been used). Still, generally, in an aerodynamic problem, a parameter describing a downstream distance should be used, e.g., wing chord c in this case.

This latter dimensionless parameter is known as the Reynolds number, which is denoted by the symbol Re. The Reynolds number is one of the most significant parameters in aerodynamics because it governs the relative magnitude of viscous effects to inertia effects in the flow.

Third Π Product

Similarly, for the third \Pi product, then

(44)   \begin{equation*} \Pi_{3} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (S)^{\gamma} \, a_{\infty} \end{equation*}

and so, in terms of dimensions, then

(45)   \begin{equation*} \left[ \Pi_{3} \right] = \rm M^{0}L^{0}T^{0} = \left( \rm ML^{-3} \right)^{\alpha} \left( \rm LT^{-1} \right) ^{\beta} \left( \rm L^2 \right)^{\gamma} \left(\rm LT^{-1} \right) \end{equation*}

giving

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha+\beta+2\gamma+1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta-1 \end{eqnarray*}

Therefore, in this case \alpha = 0, \beta = -1, and \gamma = 0, i.e.,

(46)   \begin{equation*} \Pi_{3} = \varrho_{\infty}^{0} V_{\infty}^{-1} S^{0} a_{\infty} = \frac{a_{\infty}}{V_{\infty}} \end{equation*}

Again, the \Pi product can be inverted for convenience (to put V_{\infty} on the numerator) and still have adimensionless quantity so

(47)   \begin{equation*} \Pi_{3} = \frac{V_{\infty}}{a_{\infty}} \end{equation*}

This particular dimensionless parameter is called the Mach number (the free-stream Mach number in this case), which is the ratio of the free-stream velocity to the free-stream speed of sound. Mach number is an essential quantity in aerodynamics that affects the magnitude of the forces on a body and the general compressible nature of the flow.

Final Result for Functional Dependency

Therefore, as a result of the dimensional analysis of the wing problem, then

(48)   \begin{equation*} \phi_1 \left( \frac {R}{\frac{1}{2} \varrho_{\infty} V_{\infty}^{2} S}, \frac{\varrho_{\infty} V_{\infty} \, c }{\mu_{\infty}}, \frac{V_{\infty}}{a_{\infty}} \right) = 0 \end{equation*}

or

(49)   \begin{equation*} \phi _1 \! \left( C_{R}, Re, M_{\infty} \right) = 0 \end{equation*}

or in explicit form

(50)   \begin{equation*} C_{R} = \phi_3 \! \left( Re, M_{\infty} \right) \end{equation*}

In dimensional form, then

(51)   \begin{equation*} R =  \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} S \, \phi_{3} \! \left(Re, M_{\infty} \right) \end{equation*}

which is sometimes called Rayleigh’s equation.

The previous result is a significant outcome because it means the original aerodynamic problem in the wind tunnel has been reduced from one involving five independent variables (for which it is practically impossible to determine the functional relationships to R) to one that now only has two variables! Therefore, with a little effort, it has been shown that to determine the force R on a given wing at a given angle of attack \alpha, then only the dependence of Re and M_{\infty} on R need to be established. The dimensionless parameters Re (the Reynolds number) and M_{\infty} (the Mach number) are known as similarity parameters.

The plot below shows the actual results from the wind tunnel test to measure the resultant force coefficient, C_R, as a function of the wing’s angle of attack, \alpha, to the free stream flow. The measurements were performed at two different values of the Reynolds number, which in this case was obtained by varying the wind speed. Notice the significant effect of a higher Reynolds number, which increases the angle of attack and the maximum force coefficient before the onset of stall. Look carefully at the plot below. Why does the resultant force vary depending on whether the angle of attack is increasing or decreasing? In general, the values of the Reynolds and Mach numbers will substantially affect the aerodynamic forces produced on an airfoil section and on an entire wing, which is discussed in more detail in the following chapters.

The resultant force coefficient on a wing as a function of its angle of attack at two values of the Reynolds number based on the wing chord.

Why cannot both a and V be chosen as repeating variables?

The answer is that they have the same units and are not linearly independent. This can be mathematically verified by finding the determinant of the dimensional matrix in this case, i.e.

    \[ \begin{array}{l|r|r|} \mbox{\small Base~unit} & a_{\infty} & V_{\infty} \\ \hline \text{\small Length, L} \quad & 1 & 1 \\ \text{\small Time, T} \quad & -1 & -1 \end{array} \]

The determinant of the matrix gives

    \[ \det \left[ \begin{array}{rr} 1 & 1 \\ -1 & -1 \end{array} \right] = 1 \times (- 1 ) - (-1) \times 1 = 0 \]

Because the determinant is zero, then a_{\infty} and V_{\infty} are not linearly independent, and so both variables together cannot be chosen as repeating variables. As a final step, it is easy to prove that powers of parameters with the same dimensions are also not linearly independent.

Adding Wing Span to the Problem

In a final step, consider what would happen if wing span b was added as a seventh variable to the parameter list in this previous problem, i.e., in this case, assume that

    \[ \psi \left( \varrho_{\infty}, V_{\infty}, \mu_{\infty}, a_{\infty}, c, b  \right) = 0 \]

Therefore, N = 7, K = 3, N - K = 4, and so four \Pi products will result from the dimensional analysis. The first three have already been determined. It does not take much work to show that the fourth dimensionless grouping will be

    \[ \Pi_{4} = \frac{b}{c} = AR \]

which is a span-to-chord ratio known as the aspect ratio of the wing, A\!R. Therefore, in functional form, the resultant force can be written as

(52)   \begin{equation*} C_{R} =\psi \, \left( Re, M_{\infty}, AR \right) \end{equation*}

The outcome is that by measuring the aerodynamic characteristics of a wing as a function of its aspect ratio, free-stream Mach number, and Reynolds number, the wing’s behavior can be uniquely defined for any combination of flight conditions. Notice that the aspect ratio measures the slenderness of the wing. Wings that are long and narrow will have a higher aspect ratio. The results in the figure below show the effects on the resultant force on the wing by decreasing its aspect ratio from 4 to 3. This particular wing used in the wind tunnel has replaceable tip extensions to allow changes in aspect ratio and the addition of different tip shapes, including winglets. Notice that the higher aspect ratio gives a higher force for a given angle of attack. Furthermore, the drag is reduced with a higher aspect ratio wing.

Increasing the aspect ratio of a wing not only increases the lift for a given angle of attack but also decreases its drag.

Measurements, as well as theory, show that the value of C_R at a given Reynolds and Mach number is inversely proportional to the aspect ratio, i.e.,

(53)   \begin{equation*} C_R \quad \propto \quad \dfrac{C_R (A\!R = \infty)}{A\!R} \end{equation*}

so the higher the aspect ratio, the greater the force coefficient. The theory also gives the result that the drag (induced drag) is given by

(54)   \begin{equation*} C_{D_{i}} \quad \propto \quad \dfrac{(C_R (A\!R = \infty))^2}{A\!R} \end{equation*}

Therefore, wings with a high aspect ratio are more efficient because for a given angle of attack, and all other parameters such as Mach number and Reynolds number being held constant, a wing with a higher aspect will produce more lift and less drag. The additional results shown in the figure below confirm the influence of the aspect ratio on the resultant force coefficient.

Practicing Exemplars

The best way of learning the method of dimensional analysis is by example. Practice first by checking and balancing equations for dimensional homogeneity to ensure the terms have the correct dimensions and units. Then, practice problems using the Buckingham \Pi method to find dimensionless groupings. After working through two or three more examples, one will quickly realize the power of dimensional analysis and how this technique can help in engineering problem-solving. Naturally, the method is not restricted to aerodynamics, and the process of dimensional analysis is helpful in all branches of engineering and the physical sciences.

Check Your Understanding #1 – Balancing equations for dimensional homogeneity

The kinetic energy, KE of a rotating propeller, depends on its moment of inertia, {I}, and its angular velocity, \Omega. The relationship can be expressed as

    \[ KE = C \, I^{\alpha} \, \Omega^{\beta} \]

where C is a dimensionless constant. Write down the base dimensions of each of the parameters. Then, by satisfying dimensional homogeneity, determine the values of \alpha and \beta. Show all steps in your analysis. Hint: Moment of inertia has dimensions of mass times length squared. Show ALL steps in your analysis.

Show solution/hide solution.

We are given that

    \[ KE = C \, I^{\alpha} \, \Omega^{\beta} \]

where C is a dimensionless constant. In terms of the dimensions

    \[ \left[ KE\right] = \left[ I^{\alpha} \right] \left[ \Omega^{\beta} \right] \]

In terms of dimensions for each parameter, then

    \[ \left[ KE \right] = \rm M L^2 T^{-2} \quad\quad \left[ I \right] = \rm M L^{2} \quad\quad \left[ \Omega \right] = \rm T^{-1} \]

so

    \[ \rm M L^2 T^{-2} = \left(\rm M L^{2} \right)^\alpha \, \left( \rm T^{-1}\right)^\beta \]

For dimensional homogeneity, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad  1 & = & \alpha \\ \text{\small Exponents~of~L:} \quad  2 & = & 2 \alpha \\ \text{\small Exponents~of~T:} \,\, -2 & = & -\beta \end{eqnarray*}

so \alpha = 1 and \beta = 2. Inserting the values gives the relationship as

    \[ KE = C \, I^{1} \, \Omega^{2} = C \, I \, \Omega^{2} \]

Therefore, the final answer is

    \[ KE = C \, I \, \Omega^{2} \]

Now, we can check the final units to confirm that the equation is dimensionally balanced. In this case,

    \[ \left[ KE \right]  = \rm M L^2 T^{-2} \]

and

    \[ \left[ I \Omega^2 \right]  = \left( \rm M L^{2} \right) \rm T^{-2} = \rm M L^2 T^{-2} \]

which confirms that the equation KE = C \, I \, \Omega^{2} is indeed dimensionally balanced.

Check Your Understanding #2 – Drag on an aerosol particle

As shown in the figure below, a spherical aerosol particle of diameter d (assumed in the sub-micron size range) falls freely vertically at velocity V in air. Assume that the fluid’s viscosity is much more important than its density, so density effects can be ignored. Find the similarity parameter that governs the drag on the particle.

Illustration of a tiny spherical particle falling freely in the vertical direction.

Show solution/hide solution.

Functional Dependence

The aerodynamic drag on the particle D can be written in functional form as

    \[ D = \phi ( d, \, V, \, \mu ) \]

where \mu is the viscosity of the air and \phi is some function to be determined. In implicit form, then

    \[ \psi (d, \, V, \, \mu, \, D ) = 0 \]

where \psi is some other function. In this case, N = 4 and K = 3 (by inspection, all of mass, length, and time are involved in this problem), so there are N - K = 1, i.e., one \Pi product.

Dimensional Matrix

Setting down the dimensional matrix, then

    \[ \begin{array}{l|r|r|r|r|r} \mbox{\small Base~unit} & d & V & \mu & D \\ \hline \text{\small Mass, M}  \quad & 0 & 0 & 1 & 1 \\ \text{\small Length, L}  \quad & 1 & 1 & -1 & 1 \\ \text{\small Time, T}  \quad & 0 & -1 & -1 & -2 \end{array} \]

Π product

The dependent variable D in this case) cannot be a repeating variable, and so for the only \Pi product, then

    \[ \Pi = (d)^{\alpha} (V)^{\beta} (\mu)^{\gamma} \, D \]

and in terms of dimensions, then

    \[ \left[ \Pi \right] = \rm  M^{0}L^{0}T^{0} =\rm (L)^{\alpha} \rm (L T^{-1})^{\beta} \rm (M L^{-1} T^{-1})^{\gamma} \rm M L T^{-2} \]

For the \Pi group to be dimensionless, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \gamma + 1 \\ \text{\small Exponents~of~L:} \quad 0 & = & \alpha + \beta - \gamma + 1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\beta - \gamma - 2 \end{eqnarray*}

Therefore, \gamma = -1, and \beta = -1 and \alpha = -1 so the \Pi product is

    \[ \Pi = (d)^{-1}(V)^{-1} (\mu)^{-1} D = \frac{D}{d V \mu}  =  C_D \]

which is a drag coefficient C_d applicable, in this case, to what is known as a Stokes flow.

Check

It is always good to check the final units of the grouping to be sure it is dimensionless, i.e., it has units of 1. So, in this case

    \[ \left[ C_D \right]  = \frac{\rm M L T^{-2}}{(\rm L) (\rm L T^{-1}) (\rm M L^{-1} T^{-1})} = \frac{\rm M L T^{-2}}{\rm M L T^{-2}}  =  1 \]

which confirms that the grouping is dimensionless.

What happens when things go wrong?

Sometimes, when using the Buckingham \Pi method, it is impossible to find a set of repeating variables that are linearly independent or a set that includes all of the fundamental dimensions of the problem.  In some other cases, the solution to the problem of finding the powers of the parameters becomes indeterminate, i.e., we cannot uniquely solve for the dimensionless groupings. The accepted solution to this dilemma is to reduce the number of repeating variables by one and so create an extra \Pi grouping. If this approach is used, then the dimensional analysis process should now work, at least mathematically. However, the groupings may not be unique, and some experimentation may be required to eliminate different repeating variables to find acceptable groupings that make physical sense.

Check Your Understanding #3 – Equilibrium of a soap bubble

The equilibrium state of a soap bubble of diameter D depends on a differential pressure balance e, \Delta p,  between its internal pressure relative to the ambient pressure as well as its surface tension, \sigma. Find the non-dimensional grouping that expresses this relationship.

Show solution/hide solution.

Functional Dependency

The pressure difference can be written in functional form as

    \[ \Delta p = \phi (D, \, \sigma ) \]

or in implicit form as

    \[ \psi (D, \, \sigma, \, \Delta p )= 0 \]

Therefore, N = 3, K = 3 and so there are no \Pi products! This issue is resolved by following the rule, setting K = 2, and selecting D and \sigma as the repeating variables (there is no other choice). Therefore, there is one \Pi grouping, i.e.,

    \[ \Pi = \phi_1(D, \, \sigma, \, \Delta p ) \]

Dimensional Matrix

The dimensional matrix is

    \[ \begin{array}{l|r|r|r|r|r} \mbox{\small Base~unit}   & D  & \sigma & \Delta p  \\ \hline \text{\small Mass, M}  \quad  & 0 & 1 & 1 \\ \text{\small Length, L}  \quad  & 1 & 0 & -1 \\ \text{\small Time, T}  \quad & 0 & 0 & -2 \end{array} \]

Derivation of the Π Product

The \Pi product is

    \[ \Pi = (D)^{\alpha} (\sigma)^{\beta} \Delta p \]

where the values of the coefficients \alpha and \beta must be obtained to make the equation dimensionally homogeneous and dimensionless. In terms of the dimensions of the parameters, then

    \[ \left[ \Pi \right] = \rm M^{0}L^{0}T^{0} = \left( \rm L \right)^{\alpha} \left( \rm M T^{-2} \right)^{\beta}  \left(\rm M L^{-1} T^{-2}\right) \]

For \Pi to be dimensionless, then the powers or exponents of M, L, and T must add to zero, i.e., in this case

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \beta + 1 \\ \text{\small Exponents~of~L:} \quad 0 & = & \alpha - 1  \\ \text{\small Exponents~of~T:} \quad 0 & = & -2\beta - 2 \end{eqnarray*}

Therefore, \alpha = 1, \beta = -1, so the \Pi product is

    \[ \Pi = \dfrac{\Delta p \, D}{\sigma} \]

Check

Check the final units of the grouping to be sure it is dimensionless, i.e., it has units of 1. So, in this case

    \[ \left[ \dfrac{\Delta p \, D}{\sigma} \right]  = \dfrac{ (M L^{-1} T^{-2}) L }{M T^{-2} } = 1 \]

Check Your Understanding #4 – Orbital decay of a satellite in low Earth orbit

The orbital decay of a satellite in low earth orbit (LEO), given by \overbigdot{R_s} (= dR_s/dt), is known to be a function of the orbital radius, R_s, the density of the air, \varrho,  in the upper atmosphere, the size of the satellite in terms of a reference area, A_{\rm ref}, the aerodynamic drag on the satellite, D_s, the mass of the satellite, m_s, and acceleration under gravity, g. Find the dimensionless groupings that will govern this problem.

Show solution/hide solution.

Functional Dependency

By using the information provided then, the orbital decay in LEO can be written in an explicit functional form as

(55)   \begin{equation*} \overbigdot{R_s} = \psi_1 \left( \varrho, \, A_{\rm ref}, \, D, \,m_s, \, g, \, R_s, \,V_s \right) \end{equation*}

where \Psi_1 is some function, or the information can be written in implicit form as

    \[ \psi_2 \left( \overbigdot{R_s} , \varrho, \, A_{\rm ref}, \, D, \, m_s, \, g, \, R_s, \, V_s \right) = 0 \]

where \psi_2 is some other function. Therefore, N = 8, K = 3, and so there are five \Pi products to find in this case. This is a little more complex problem in dimensional analysis, but the principles are the same.

Repeating Variables

The repeating variables must collectively include the units of mass {M}, length {L}, and time {T}, so one choice here is \varrho, {R_s}, and {V_s},  all of which have a primary dependence on the orbital decay and are also linearly independent. Alternatively, \varrho, A_{\rm ref}, and {V_s} could be used. But if g were to be chosen instead of {\varrho} as a repeating variable, then collectively, the repeating variables would not include mass, {M}; the Buckingham \Pi method will fail in this case. Notice also that it would not be appropriate to use the satellite’s mass, m_s, as a repeating variable because the orbital characteristics, such as velocity and period, are independent of the mass. The dependent variable, \overbigdot{R_s}, cannot be used as a repeating variable.

Π Products

Therefore, the five \Pi products are formed from

    \[ \Pi_1 = \phi_1(\varrho, R_s,  V_s, \overbigdot{R_s} ) \]

    \[ \Pi_2 = \phi_2(\varrho, R_s,  V_s, D ) \]

    \[ \Pi_3 = \phi_2(\varrho, R_s,  V_s, m_s ) \]

    \[ \Pi_4 = \phi_2(\varrho, R_s,  V_s, g ) \]

    \[ \Pi_5 = \phi_2(\varrho, R_s,  V_s, A_{\rm ref} ) \]

Dimensional Matrix

Forming the dimensional matrix from the base units of each of the parameters gives

    \[ \begin{array}{l|r|r|r|r|r|r|r|r} \mbox{\small Base~unit} & \overbigdot{R_s} & \varrho & A_{\rm ref} & D & m_s & g & R_s & V_s \\ \hline \text{\small Mass, M} & 0 & 1 & 0 & 1 & 1 & 0 & 0 & 0 \\ \text{\small Length, L} & 1 & -3 & 2 & 1 & 0 & 1 & 1 & 1 \\ \text{\small Time, T} & -1 & 0 & 0 & -2 & 0 & -2 & 0 & -1 \end{array} \]

Derivation of the Π Products

The first \Pi product is formed from

    \[ \Pi_{1} = (\varrho)^{\alpha} (R_s)^{\beta} (V_s)^{\gamma} \, \overbigdot{R_s} \]

where the values of the coefficients \alpha, \beta, and \gamma must be obtained to make the equation dimensionally homogeneous and dimensionless. In terms of the dimensions of the parameters, then

    \[ \left[ \Pi_{1} \right] = \rm M^{0}L^{0}T^{0} = \left( \rm M L^{-3} \right)^{\alpha} \left( \rm L \right)^{\beta} \left( \rm L T^{-1} \right)^{\gamma} \left(\rm L T^{-1}\right) \]

For \Pi_{1} to be dimensionless, then the powers or exponents of M, L, and T must add to zero, i.e.,

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha + \beta + \gamma + 1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\gamma - 1 \end{eqnarray*}

Therefore, in this case \alpha = 0, \gamma = -1, and \beta = 0, so the \Pi_1 product is

    \[ \Pi_1 = (\varrho)^{0} (R_s)^{0} (V_s)^{-1} \, \overbigdot{R_s} = \dfrac{\overbigdot{R_s}}{V_s} \]

which is a dimensionless decay rate that is a fraction of the orbital velocity.

The second \Pi product involves the drag force, D, i.e.,

    \[ \Pi_2 = (\varrho)^{\alpha} (R_s)^{\beta} (V_s)^{\gamma} \, D \]

In terms of the dimensions, then

    \[ \left[ \Pi_2 \right] = \rm M^{0}L^{0}T^{0} = \left( \rm M L^{-3} \right)^{\alpha} \left( \rm L \right)^{\beta} \left( \rm L T^{-1} \right)^{\gamma} \left(\rm M L T^{-2}\right) \]

For \Pi_2 to be dimensionless, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha + 1 \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha + \beta + \gamma + 1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\gamma - 2 \end{eqnarray*}

Solving these equations gives \alpha = -1, \gamma = -2, and \beta = 2, so the \Pi_2 product is

    \[ \Pi_2 = (\varrho)^{-1} (R_s)^{2} (V_s)^{-2} \, D = \frac{D}{\varrho \, R_s^2 \, V_s^2} \]

which is a form of drag coefficient. It is legitimate to replace R_s^{2} by A_{\rm ref} and multiply the denominator by 1/2 to obtain the conventional form of drag coefficient, i.e.,

    \[ \Pi_2 = \dfrac{D}{\dfrac{1}{2} \varrho \, A_{\rm ref} \, V_s^2} \]

The third \Pi product involves the satellite mass, m_s, i.e.,

    \[ \Pi_3 = (\varrho)^{\alpha} (R_s)^{\beta} (V_s)^{\gamma} \, m_s \]

In terms of the dimensions, then

    \[ \left[ \Pi_3 \right] = \rm M^{0}L^{0}T^{0} = \left( \rm M L^{-3} \right)^{\alpha} \left( \rm L \right)^{\beta} \left( \rm L T^{-1} \right)^{\gamma} \, \rm M \]

For \Pi_3 to be dimensionless, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha + 1 \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha + \beta + \gamma \\ \text{\small Exponents~of~T:} \quad 0 & = & -\gamma \end{eqnarray*}

Solving these, gives \alpha = -1, \gamma = 0, and \beta = -3, so the \Pi_3 product is

    \[ \Pi_3 = (\varrho)^{-1} (R_s)^{-3} (V_s)^{0} \, m_s = \frac{m_s}{\varrho \, R_s^3} \]

or inverting the grouping because a low value \varrho will make the parameter large, then

    \[ \Pi_3 = \frac{\varrho \, R_s^3}{m_s} \]

which is a form of buoyancy ratio. The value of this parameter in space is likely to be so small, even in a low Earth orbit, that it can be justifiably neglected.

The fourth \Pi product involves the gravitational acceleration, g, i.e.,

    \[ \Pi_4 = (\varrho)^{\alpha} (R_s)^{\beta} (V_s)^{\gamma} \, g \]

In terms of the dimensions, then

    \[ \left[ \Pi_4 \right] = \rm M^{0}L^{0}T^{0} = \left( \rm M L^{-3} \right)^{\alpha} \left( \rm L \right)^{\beta} \left( \rm L T^{-1} \right)^{\gamma} \left(\rm L T^{-2}\right) \]

For \Pi_4 to be dimensionless, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha + \beta + \gamma + 1 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\gamma - 2 \end{eqnarray*}

Solving these, gives \alpha = 0, \gamma = -2, and \beta = 1, so the \Pi_4 product is

    \[ \Pi_4 = (\varrho)^{0} (R_s)^{1} (V_s)^{-2} \, g = \frac{g \, R_s}{V_s^2} \]

which is a ratio of potential energy to kinetic energy.[1]

The fifth \Pi product involves the reference area, A_{\rm ref}, i.e.,

    \[ \Pi_5 = (\varrho)^{\alpha} (R_s)^{\beta} (V_s)^{\gamma} \, A_{\rm ref} \]

In terms of the dimensions, then

    \[ \left[ \Pi_5 \right] = \rm M^{0}L^{0}T^{0} = \left( \rm M L^{-3} \right)^{\alpha} \left( \rm L \right)^{\beta} \left( \rm L T^{-1} \right)^{\gamma} \left(\rm L^2\right) \]

For \Pi_5 to be dimensionless, then

    \begin{eqnarray*} \text{\small Exponents~of~M:} \quad 0 & = & \alpha \\ \text{\small Exponents~of~L:} \quad 0 & = & -3\alpha + \beta + \gamma + 2 \\ \text{\small Exponents~of~T:} \quad 0 & = & -\gamma \end{eqnarray*}

Solving these equations, gives \alpha = 0, \gamma = 0, and \beta = -2. Therefore, the \Pi_5 product is

    \[ \Pi_5 = (\varrho)^{0} (R_s)^{-2} (V_s)^{0} \, A_{\rm ref} = \frac{A_{\rm ref}}{R_s^2} \]

which is an area ratio, although of little significance in this problem.

Final Dimensionless Form

In summary, the five \Pi products are:

    \[ \Pi_1 = \frac{\overbigdot{R_s}}{V_s}, \quad \Pi_2 = \frac{D}{\dfrac{1}{2}\varrho \, A_{\rm ref} \, V_s^2}, \quad \Pi_3 = \frac{\varrho \, R_s^3}{m_s}, \quad \Pi_4 = \frac{g \, R_s}{V_s^2}, \quad \Pi_5 = \frac{A_{\rm ref}}{R_s^2} \]

The final relationship can be written as a function of these dimensionless groups as

    \[ \phi_1\left( \frac{\overbigdot{R_s}}{V_s}, \ \dfrac{D}{\dfrac{1}{2}\varrho \, A_{\rm ref} V_s^2}, \ \frac{\varrho \, R_s^3}{m_s}, \ \frac{g \, R_s}{V_s^2}, \ \frac{A_{\rm ref}}{R_s^2} \right) = 0 \]

Removing the irrelevant or redundant terms in writing in explicit form gives

    \[ \frac{\overbigdot{R_s}}{V_s} = \psi_2 \left( \dfrac{D}{\dfrac{1}{2}\varrho \, A_{\rm ref} \, V_s^2}, \ \frac{g R_s}{V_s^2} \right) \]

Therefore, the dimensionless radial decay of a satellite in LEO is a function of the aerodynamic drag coefficient of the satellite (C_D = 2.0 is often assumed), and the ratio of the potential to kinetic energy (usually about 1).

Other Dimensionless Quantities

In engineering practice, it is common to work with dimensionless quantities or to normalize dimensional quantities by appropriate reference values to make them dimensionless. This approach simplifies mathematical expressions and allows the use of standard mathematical functions without introducing units.

Some physical variables are naturally dimensionless by their definition, i.e., as ratios of the same dimensional quantities. One example often encountered in thermodynamics and aerodynamics is the ratio of specific heats, \gamma = C_P / C_{\cal{V}}. Another example of a naturally dimensionless parameter used in structural analysis is strain, which is defined as the ratio of the change in length of a structural component under loading, \epsilon, compared to its initial length, L, i.e., the strain is \epsilon / L. Finally, in hydrostatics, the specific gravity, SG, which is a dimensionless number defined as the ratio of the density of a given liquid relative to the density of water, is often used.

It should be remembered that all angles are dimensionless. In engineering analysis, angles measured in radians are always used, defined as a ratio of the arc length to the radius of a circle. The dimensions of a radian are 1, i.e., it is unitless. While angles are often measured or reported in degrees, units of degrees must be converted to radians for most engineering calculations. A complete circle subtends an angle of 2\pi radians. which is equivalent to 360 degrees or 360^{\circ}, so the angular conversion is 1 radian = 180/\pi degrees. The same should also be done with angular rates, which should always be converted to radians per second (rad/s) no matter what units the rates are originally given in, e.g., rates are often measured in revolutions per minute (rpm). The dimensions of angular rate are T^{-1}.

When dealing with trigonometric and exponential functions, the arguments must be dimensionless; using them with dimensional values is incorrect. Consider the following examples:

  • In trigonometry, angles are measured in radians, which is a dimensionless angle. Therefore, all trigonometric expressions such as \sin (x) and \tan(y) are valid only when x and y are in units of radians. For example, if x has units of length (typical in engineering), then the equation z = \sin (x) is invalid. However, if the lengths are non-dimensionalized by a fixed reference length such as L, then z/L = \sin (x/L) is a valid equation.
  • The exponential function \exp(t), which occurs frequently in engineering, is also defined for dimensionless arguments. If the variable t has units of time, then the equation y = \exp(-t)  is invalid. However, if \omega is a frequency (units of radians per second), then y = \exp(-\omega \, t)  is a valid equation. Likewise, if \tau is a reference time (units of seconds), then y = \exp(- t / \tau)  is a valid equation.

The original dimensional analysis paper by Edgar Buckingham.

The original dimensional analysis paper by Edgar Buckingham dates to 1914. Notice his comments on dealing with the square root of parameters (“dispense with fractional exponents”), the choice of repeating variables, and the inversion of dimensionless groupings. Later in the paper, he talks about the “screw propeller” problem and dynamic similarity, which he recognizes as a powerful tool for scaling measurements made on models and relating them to full scale. In 1914, he commented, “This subject seems to be less familiar to physicists than it deserves to be…” Today, we take advantage of dimensional analysis and the principles of dynamic similarity in much of what we do as engineers. It sometimes pays to step back and remember how we got this far.

Summary & Closure

Dimensional analysis can be a powerful tool for engineering analysis because it reduces and simplifies the size and scope of the problem in terms of dimensionless or non-dimensional groupings called similarity parameters. The strength of dimensional analysis as a tool usually becomes apparent when applied to increasingly more complicated engineering problems, such as where the effects of all the dependencies cannot be established explicitly and/or the governing equations for the solution of the problem can become formidably complicated. For example, it is common for engineering problems to have five, six, seven, or even more dependencies. However, sorting out the problem parameters into dimensionless groupings helps to understand what parameters are essential as they will affect the physical behavior and what parameters are not.

Dimensional analysis must be used in a systematic and structured approach, formally embodied in the Buckingham \Pi method. A critical first step is to carefully lay out the problem in terms of the dependent and independent quantities, which must all be assigned the correct base units. While the choice of the repeating quantities is somewhat arbitrary, specific rules must be followed for the method to work. For example, it is critical that the dependent quantity not be used as a repeating variable. The repeating variables must also have substantial effects on the dependent quantity and be linearly independent of each other, i.e., two repeating parameters with the same dimensions cannot be used, such as two reference length scales or two reference velocities. If the basic rules of the Buckingham \Pi method are followed, then the outcomes will be groups of dimensionless parameters that should not only help understand the problem but will simplify its scope.

5-Question Self-Assessment Quickquiz

For Further Thought or Discussion

  • You are tasked with explaining to a high-school student studying physics the advantages of using dimensional analysis and dimensionless quantities. What two primary arguments would you make?
  • The force on a body in the wind tunnel is known to depend on flow velocity (airspeed) and the static pressure and temperature in the flow. What other parameters may be involved? Hint: Do not suggest parameters such as viscosity or the speed of sound.
  • Use dimensional analysis to show that the drag of a ship’s hull, in part, depends on the Froude number. The Froude number is defined as Fr = V_{\infty}/\sqrt{L g}, where L is the length of the ship and V_{\infty} is the speed of the ship through the water.
  • You are tasked with determining the hydrodynamic drag force D on a fully submerged, very long, neutrally buoyant cable dragged behind a ship. You propose to do this by doing small-scale experiments in a water tunnel. Determine the dimensionless groupings that will affect this problem.
  • Show that the pressure drop caused by a restrictor in a flow can be expressed in terms of the Euler number Eu, where the pressure drop is expressed in terms of an energy loss in the flow per unit of kinetic energy per volume of the flow.

Other Useful Online Resources

To dive deeper into the methods of dimensional analysis, navigate to the following online resources:

 


  1. This outcome will be apparent by multiply the numerator and the denominator by the satellite's mass.

License

Icon for the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License

Introduction to Aerospace Flight Vehicles Copyright © 2022–2024 by J. Gordon Leishman is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, except where otherwise noted.

Digital Object Identifier (DOI)

https://doi.org/https://doi.org/10.15394/eaglepub.2022.1066.n12