13 Dimensional Analysis

Introduction

New engineers quickly begin to recognize from even the most casual observations that the solution to engineering problems usually involves the consideration of many parameters as well as parameter dependencies and interdependencies. In other words, many quantities could affect the physical behavior of a given system, and different effects will be obtained when varying one parameter versus another.

In particular, the behavior of fluids is relatively difficult to describe because they deform, especially as they flow around objects (often referred to as “bodies”). The fluid also produces forces and moments on bodies, including the familiar lift and drag forces. 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 orientation to the flow. The forces will 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.

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 they form the basis of a tool called computational fluid dynamics (CFD). However, this is not a straightforward approach, not just because of the complexity of the equations that must be solved but also 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 a wind tunnel. Again, however, this would not be an easy task, partly because of the number of measurements needed.

At least in principle, either or both of these latter approaches will allow the effects of the fluid on the body to be determined and the reciprocal effects 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 non-dimensional groupings called similarity parameters. These parameters then allow  a significant simplification in understanding the effects of the original problem parameters. Dimensional analysis is not restricted to fluid problems and is helpful in many fields of engineering and science.

Learning Objectives

  • Review the system of SI (metric) and USC (British Imperial) units of measurement and their proper use.
  • Understand how to perform dimensional analysis, a method used to bring out the primary dependencies in an engineering problem in terms of a non-dimensional grouping called similarity parameters.
  • Better understand the significance of similarity parameters and how they can be used to simplify a problem and reduce the number of dependencies, i.e., making a problem easier to understand.
  • Develop and discuss some exemplar problems to become familiar with the method and practice of dimensional analysis and the use of similarity parameters.

All About Units

Before examining the process of dimensional analysis, it is helpful to review the subject of units and measurement systems. A unit is a standardized measurement of a parameter or quantity that has been adopted by consensus and defined by convention, e.g., the SI system. The expression of quantities that are measured in terms of standardized units allows their magnitude or numerical value to be evaluated and compared. Adopting a units system for the measurement of a broad range of physical parameters allows engineers and scientists across the world to communicate and work with each other using a common standard.

The value of a particular physical quantity is always given the product of a number and a unit (unless unitless), e.g., 12.4 meters per second (m/s) in SI units is equal to 40.68 feet per second (SI) in USC or British Imperial units, where the conversion factor from meters to feet is 3.281. This means the numerical value of a quantity depends on the system of units in which it is expressed. Engineers in the U.S. typically use both SI and USC units. However, there is a preference to use SI units because of their widespread adoption and use in scientific research and compliance with international standards. Conversion factors are available so that numerical values in one unit system can be expressed as equivalent values in the other.

SI Units

The Système International or International System of Units (SI), sometimes known as the metric system, has become the international standard for measurement units. The SI system is also referred to as the decimal system because it is based on powers of 10. The SI system comprises seven base units with many other derived units that are given special names and symbols, with a subset being shown in the table below. The SI system is based on the meter (m) (or the metre if outside the U.S.) as the base unit of length, the kilogram (kg) as the base unit of mass, the Newton (N) as the base unit of force, and the second (s) as the base unit of time. Notice that unit symbols are never italicized.
Definitions for common SI units.
Quantity Quantity symbol SI unit name Unit symbol
Length L meter m
Mass M kilogram kg
Time T second s
Temperature \theta Kelvin K or oK
Velocity V meter per second m/s
Acceleration a meter per second2 m/s2
Force F Newton N
Energy E Joule J
Power P Joule per second = Watt W
Pressure p Newtons per m2 = Pascal Pa
Viscosity \mu Pascal-second Pa s

USC Units

U.S. customary (USC) units is a system of measurement units most commonly used in the U.S. The USC system was developed from traditional British measurement units. The British measurement system was updated in 1824, which created the so-called Imperial system or British Imperial units of measurement, often known in the U.S. as “English” or “British Engineering” units. For engineering use, the USC or “British system” is based on the foot (ft) as the base unit of length, the slug (slug or sl) as the base unit of mass, the pound (lb) as the base unit of force, and the second (s) as the base unit of time, as reviewed in the table below.
Definitions for USC units. *Power is commonly measured in horsepower (hp) where 1 hp = 550 ft-lb/s.
Quantity Quantity symbol USC unit name Unit symbol
Length L foot ft
Mass M slug slug or sl
Time T second s
Temperature \theta Rankine R or oR
Velocity V feet per second ft/s
Acceleration a feet per second2 ft/s2
Force F pound lb
Energy E foot-pound ft-lb
Power* P foot-pounds per second ft-lb/s
Pressure p pounds per foot2 lb/ft2
Viscosity \mu slug/(ft s) or pound seconds per foot slug/(ft s) or lb s/ft

Notice that the name “slug” as a unit of mass originates from the concept of mass as an inertia or “sluggish” and has been referred to as the “engineer’s mass unit.” The unit of  a “slug” has also been associated with the mass of a lead projectile, but because a slug as a unit mass has a corresponding weight of 32.17 lb, then this is a more unlikely association.

The use of “pounds-mass” (lbm) as a measure of mass is generally considered to be a source of confusion in the engineering and scientific fields. Therefore, the unit of “pounds-mass” is is best avoided in preference of using the “slug” as the base mass unit. While USC units are largely similar to British Imperial units, there are differences in regard to volumetric measurements.

What is a horsepower?

The output of an engine is often measured in “horsepower,” which is given the unit symbol “hp.” This unit is attributed to the Scottish engineer James Watt, whose work on steam engines was fundamental to the changes brought about by the Industrial Revolution. James Watt wanted to compare the power output of his new steam engines to that of horses to help market them. Watt did experiments and determined that a typical farm horse could, on average, steadily lift a 600 lb weight through a simple pulley system an average distance of 63.9 ft in 69 seconds. Therefore, the work done is

    \[ \mbox{Work} = F \, d = 600 \times 63.9 = 38,340~\mbox{ft-lb} \]

Power is the rate of doing work, so the work per unit of time is

    \[ \mbox{Power} = \frac{\mbox{Work}}{\mbox{Time}} = \frac{38,340}{69} = 555.65~\mbox{ft-lb s$^{-1}$} \]

James Watt settled on the result that 1 hp = 550 ft-lb s^{-1}; he was not too concerned about accuracy. All he wanted was a simple but representative quantitative measure of the power delivered by a horse relative to what his steam engines could produce and so better market his engines. The unit made sense to farmers and others using horses to move equipment, so for a steam engine of 10 hp, for example, the purchaser knew exactly what they were buying. The unit of  “horsepower” is still used today almost worldwide.

In 1975, the U.S. government declared the international SI system to be the “preferred” unit system but did not require it. Indeed, the USC system still prevails in the U.S. as the preferred measurement system, and the U.S. is one of the few countries in the world that does not exclusively use the SI system. In Britain, the use of mixed SI and Imperial units is also common, but not as much as in the U.S. The SI units system is preferred in all engineering and scientific work, and many situations it may be required by specification, international code, or regulation.  However, in today’s global community, engineers and scientists from all countries need to be conversant in the use of both SI and USC units.

NOTE: In aviation, distances are measured in terms of nautical miles, which is neither Imperial or SI, airspeed is measured in terms of knots (nautical miles per hour), and altitude is measured in feet, the use of such units being in accordance with ICAO aviation and aeronautical standards.

Working with units

Student engineers need to become comfortable with the use of both SI and USC units. This end, when presented with a problem in SI units, then the problem should be solved entirely in SI units. Likewise, when presented with a problem in USC units, then the problem should be solved entirely in USC units. Never convert back and forth to suit yourself because this is when mistakes happen. Sometimes, information (values) will be presented in mixed units, with some values in SI and others in USC. In this case, it is important to convert all of the information to either SI or USC, depending on what outcome is required. In either case, it is important that all information available s converted to base units before evaluations are conducted, in which case the final desired value or values will also be in base units. If unit conversions are required, then handy conversion factors can be found by access the NIST website from the links below. It is important that when using unit conversions, the full accuracy of the conversion factor should be used. Rounding of numerical values should only be done at the end of the calculations.

CGS Units

The centimeter–gram–second or CGS system of units is an increasingly obsolete variant of the metric system. CGS units are based on the centimeter as the base unit of length, the gram as the base unit of mass, and the second as the base unit of time. The CGS system has been largely supplanted by the use of SI units, but sometimes quantities measured in CGS units are still encountered in engineering and scientific practice. Converting between CGS and SI units is easy because the conversion factors are all powers of 10, as reviewed in the table blow, For example, 100 cm = 1 m and 1,000 g = 1 kg. The CGS unit of force is called the dyne, which is 1 g cm/s2, so a Newton (N) in SI units is equal to 100,000 dynes in CGS units.

Definitions of CGS units and conversion factors from CGS to SI units
Quantity Quantity symbol CGS name Unit symbol Unit definition SI units
Length L centimeter cm 1/100 of meter 10−2 m
Mass M gram g 1/1000 of kilogram 10−3 kg
Time T second s 1 second 1 s
Velocity V centimeter per second cm/s cm/s 10−2 m/s
Acceleration a gal (Galileo) Gal cm/s2 10−2 m/s2
Force F dyne dyn g cm/s2 10−5 N
Energy E erg erg g cm2/s2 10−7 J
Power P erg per second erg/s g cm2/s3 10−7 W
Pressure p barye Ba g/(cm s2) 10−1 Pa
Dynamic viscosity \mu poise P g/(cm s) 10−1 Pa s
Kinematic viscosity \nu Stokes St cm2/s 10−4 m2/s

When NASA lost a spacecraft because of a units mistake!

In September of 1999, the NASA Mars Climate Orbiter unexpectedly burned up and crashed into the surface of Mars after almost ten months of space travel. The project cost over $100M. It turned out that engineers had made a massive units blunder, muddling quantities between the SI and USC systems. NASA’s Jet Propulsion Laboratory (JPL) used the SI system, while Lockheed used the USC system. Because of this blunder, the navigation system flew the spacecraft too fast and too close to the atmosphere, where it burned up. Read more about this mistake here as well as some other engineering disasters that have occurred because of the muddling of units. Be sure not to make the same mistakes!

Physical Quantities

Science and engineering involve observations of the physical world; engineers aim to express observations in the most general and the simplest type of formulation. The role of mathematics is essential here, and the resulting formulation will stem directly from the properties or quantities considered in the quantitative analysis. The allowed types of properties are called physical quantities.

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 mostly 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).

List of key base quantities used in engineering mechanics.
Base quantity Symbol SI unit USC unit
Mass M kilogram (kg) slug (or slugs or sl)
Length L meter (m) foot (ft)
Time T second (s) second (s)
Temperature \theta Kelvin (^{\circ}K or K) Rankine (^{\circ}R or R)
Confused about units? Here is a short video lesson on units (SI and USC) from Dr. Leishman’s “Math & Physics Hints and Tips” series.
Uncertain about the differences between mass and weight? Here is a short video lesson on the differences between mass and weight from Dr. Leishman’s “Math & Physics Hints and Tips” series.

Finding Base Dimensions

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. From now on, notice that the parameter inside a set of square brackets, i.e., “[parameter],” is used to mean the dimensions of that parameter. Notice also that scientific notation must be used throughout to avoid ambiguity in the final units of the derived parameter, e.g., use L \, T^{-1} and not L / T. Consider the following examples:

  • Volume, {\cal{V}}, is a derived quantity because to calculate 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.
  • Density, \varrho, is also a derived parameter because the density \varrho of something is its mass (M) per unit volume (L^3), so density will have units of M L^{-3}, i..e., [\varrho] = M L^{-3}.
  • Force, F, is a derived parameter, although in some cases, “force” could be carried through the problem as a base unit. A force is the product of a mass and acceleration (Newton’s second law), so in terms of its base units, then [F] = M \,  LT^{-2}.
  • Torque, \tau, is a derived parameter because a torque is a moment, so the product of a force times a distance or “arm.” Therefore, [\tau] = (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 never the form M L / T^2.

Need to brush up on the concepts of work, energy, and power? Here is a short video lesson on these concepts from Dr. Leishman’s “Math & Physics Hints and Tips” series.

Aerodynamic Force on a Wing

To introduce the subject and methods of dimensional analysis, consider the forces produced on a wing in the 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.

A wing set at an angle to the flow, in this case being tested in the ERAU wind tunnel, will produce some resultant aerodynamic force on the wing depending on the size of the wing and the free-stream parameters.

On a physical intuitive basis, then it would be expected that the resultant aerodynamic force R on the wing for a given orientation to the flow will depend on:

  1. Free-stream velocity, V_{\infty}.
  2. Free-stream density, \varrho_{\infty}.
  3. The viscosity of the air, i.e., on the coefficient of viscosity \mu_{\infty}.
  4. 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.
  5. The sonic velocity a_{\infty}, i.e., the potential for compressibility effects.

The subscript \infty or “infinity” is used to mean undisturbed conditions, i.e., the “free-stream” or reference values of these quantities that are constant well upstream of the wing. This approach makes sense because as the flow approaches the wing, these values are all likely to change. 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 nose to the tail of the wing) and the direction of the free-stream velocity V_{\infty}.

The resultant force can also be decomposed into a vertical component that acts perpendicular to the relative wind (called the lift) and a horizontal component parallel to the relative wind (called the drag).

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

(1)   \begin{equation*} R = f \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 between the variables, which in this case is denoted by “f” is still to be determined. Although this latter equation expresses a very 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 (somehow) 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} is easy because it 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, then the number of possible experimental runs 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. Eventually, the needed relationships between R and all dependent quantities at each angle of attack would be determined. Therefore, it does not take much to figure out that this process of systematically varying each of the dependent parameters, in turn, could take the student a very long time, even by doing parallel calculations using aerodynamic theory.

This approach would also result in a large quantity of data, 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 was used, then just a single column 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 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 challenging. 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.

Basis of Dimensional Analysis

Fortunately, the approach to the preceding problem can be simplified by employing dimensional analysis, which can be used to determine sets of non-dimensional (or dimensionless) groups of parameters that affect the problem of interest. Such an approach will also reduce the number of independent variables, so it becomes an essential technique for understanding many engineering problems, not just aerodynamics.

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

(2)   \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 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. Furthermore, to make sense then each parameter must also be expressed numerically in the same units system.

This previous equation can be made non-dimensional or dimensionless, or what is called dimensionally homogeneous by dividing through by one of the terms, i.e., the same result can be expressed as

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

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

Lord Rayleigh and Joseph Bertrand originally set down these 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) is a generalized formalization of the method of dimensions, which can now be examined.

The Buckingham Π Method

Let K be equal to the number of fundamental dimensions required to describe the physical variables. In engineering, most physical variables can be expressed a combinations of mass, length, time, and temperature, so in general K=4. In mechanics, then mass, length, and time are involved, so K=3. Let P_{1}, P_{2}, ...P_{N} represent a list of N physical variables in the general physical relationship

(4)   \begin{equation*} P_{1} = f \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 f represents a general functional relationship, as yet unknown. The previous equation can also be written is implicit form as

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

where g represents another functional relationship between the parameters yet to be determined. The Buckingham \Pi method then states that the above relationship 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. One key to the method, therefore, is how to go about the process of selecting the most suitable repeating variables.

To see how this method works, let P_{1}, P_{2}, ....P_{K} be the selected K repeating variables, i.e.,

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

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

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

where \phi is the functional relationship (still to be determined) but now it is expressed in terms of each of the \Pi products. In explicit form, this latter equation can also be written as

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

where \psi is some other functional relationship.

The upshot of the dimensional analysis process means that the functional relationship has been reduced from the initial N dimensional variables to N-K non-dimensional or dimensionless variables. The reduction the number of dependent parameters, therefore, is a significant simplification in the analysis of the original problem. For example, if originally N = 6 and K = 3, i.e., the dependent variable is a function of five other parameters as in the wing problem discussed above, the dimensional analysis will reduce the problem to three non-dimensional groupings. If N = 5 and K = 3, then the dimensional analysis will reduce the problem to two non-dimensional groupings in that one grouping becomes a sole function of the other grouping.

Π 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. Perhaps the drag of a sphere sounds boring compared to a wing, but this problem helps first to introduce the process of dimensional analysis without making the process too lengthy.

Furthermore, the aerodynamic characteristics of a sphere is particularly relevant to all sorts of ball games. The motion of balls through air is a complex fluid dynamic problem that is influenced by many factors, including the density, velocity, viscosity, the size of the ball, and its surface finish. For example, the aerodynamics of a golf ball affect the distance, spin, and overall trajectory of the ball. Ever wonder why a golf ball has dimples in its surface? In baseball, the aerodynamics of a baseball 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 non-dimensional groupings that will describe 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.

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

    \[ D = f_1\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. The functional dependence in this case is denoted by f_1

The functional dependence in implicit form is

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

where f_2 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 in the form y = f(x) versus 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 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 g is NOT the same as the function f. If y = x^2 +5, for example, then y = f(x), which is an explicit dependency of 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 starts off 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. For each dependency, then

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

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

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. Setting up the dimensional matrix in this case gives

    \[ \begin{array}{l|r|r|r|r|r|r} &   D  & \varrho_{\infty} & V_{\infty} & d &  \mu_{\infty} \\ \hline \mbox{Mass: } M  &   1  & 1 & 0 & 0 & 1 \\ \mbox{Length: } L  &  1  & -3 & 1 & 1 & -1 \\ \mbox{Time: } T  &  -2 & 0 & -1 & 0 & -1 \end{array} \]

One purpose of the dimensional matrix is to lay the problem out in a clear manner and test for 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 a 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 based 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 all of 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.

If these prior rules are followed, the dimensional analysis process will likely succeed. If not, then it will fail, often catastrophically. For complex problems, of which there are many in engineering, the choice of repeating variables may not be obvious but should still be selected to be quantities that are 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 have significant effects on 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

(9)   \begin{equation*} \begin{array}{l|r|r|r|} & \varrho_{\infty} & V_{\infty} & d \\ \hline M & 1 & 0 & 0 \\ L & -3 & 1 & 1 \\ T & 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

(10)   \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

(11)   \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 zero, i.e., the test is whether

(12)   \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*}

In this case, the determinant is 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} & = & g_{1} \left( \underbrace{\varrho_{\infty}, V_{\infty}, d}_{\rm{~Repeating}}  , \underbrace{D}_{\rm {Other}} \right)  \\ \Pi_{2} & = & g_{2} \left( \underbrace{ \varrho_{\infty}, V_{\infty}, d}_{\rm{~Repeating}}  , \underbrace{\mu_{\infty}}_{\rm {Other}} \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 non-dimensional.  In terms of the dimensions of the problem then

    \[ \left[ \Pi_{1} \right] = 1 = M^{0}L^{0}T^{0} = \left( ML^{-3} \right)^{\alpha} \left( LT^{-1} \right)^{\beta} \left( L \right)^{\gamma} \left( 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. Therefore, for \Pi_{1} to be dimensionless, then

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

so that

    \begin{eqnarray*} M: 0 & = & \alpha+1  \\ L: 0 & = & -3\alpha +\beta+ \gamma+1  \\ T: 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. Normally, 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 as 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} \]

Dealing with exponents or powers

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

    \[ 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} \]

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

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

When working with exponents, the rule to remember is that the base parameter has to be the same and the process of 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

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

Equating the powers gives

    \begin{eqnarray*} M: 0 & = & \alpha +1  \\ L: 0 & = & -3\alpha +\beta+ \gamma-1  \\ T: 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 non-dimensional grouping is still obtained, i.e.,

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

Recall from previously that this non-dimensional parameter is called the Reynolds number.

Final Result for Functional Dependency

As a result of the dimensional analysis on 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 \]

or in explicit form

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

The functional dependency to be measured is now C_{D} = \phi \left( Re \right). The force can be measured with a balance, and with measurements of the density and viscosity (via measurements of pressure and temperature) the 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 results are shown in the figure below for the drag coefficient on the sphere as measured in the wind tunnel as a function of Reynolds number. 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. This behavior is tied to the characteristics of the flow on the surface of the sphere, which change from being dominated by laminar flow events to turbulent flow events.

Variation of the drag coefficient of a sphere as a function of its 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 foregoing plot is that the drag coefficient of the sphere 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 useful to think about what might happen if surface roughness was added as fifth variable to the parameter list, i.e., in this case assume that

    \[ D = f_1 \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

    \[ f_2 \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} = f_{3} \left( Re , \frac{\epsilon}{d} \right) \]

where \epsilon/d is a non-dimensional roughness height.

Π Method Applied to a Wing

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

(13)   \begin{equation*} R = f_1 \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.,

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

where f_2 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.

Dimensions of the Variables

Now, the dimensions of the variables must be determined. Remember that the square brackets [-] around a variable mean “the dimensions of” that variable. So, for each dependency in this case then

    \begin{eqnarray*} \left[ R \right] & = & MLT^{-2} \\ \left[ \varrho_{\infty} \right] & = & ML^{-3} \\ \left[ V_{\infty} \right] & = & LT^{-1} \\ \left[ S \right] & = & L^2 \\ \left[ \mu_{\infty} \right] & = & ML^{-1}T^{-1} \\ \left[ a_{\infty} \right] & = & 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.,

(15)   \begin{equation*} \begin{array}{l|r|r|r|r|r|r} & R & \varrho_{\infty} & V_{\infty} & S & \mu_{\infty} & \alpha_{\infty} \\ \hline \mbox{Mass: } M & 1 & 1 & 0 & 0 & 1 & 0 \\ \mbox{Length: } L & 1 & -3 & 1 & 2 & -1 & 1 \\ \mbox{Time: } T & -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. In a practical sense, 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 be used to 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). For example, length and area scales cannot both be used as repeating variables, which 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 non-dimensional groupings. Even though non-dimensional 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} & = & g_{1} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~ Repeating}} , \underbrace{R}_{\mbox{~Other}} \right) \\ \Pi_{2} & = & g_{2} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~Repeating}} , \underbrace{\mu_{\infty}}_{\mbox{~Other}} \right) \\ \Pi_{3} & = & g _{3} \left( \underbrace{\varrho_{\infty}, V_{\infty}, S}_{\mbox{~Repeating}} , \underbrace{a_{\infty} }_{\mbox{~Other}} \right) \end{eqnarray*}

First Π Product

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

(16)   \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 non-dimensional with units of 1. In terms of the dimensions, then

(17)   \begin{equation*} \left[ \Pi_{1} \right] = M^{0}L^{0}T^{0} = \left( ML^{-3} \right)^{\alpha} \left( LT^{-1} \right)^{\beta} \left( 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

(18)   \begin{equation*} \left[ \Pi_{1} \right] = M^{0}L^{0}T^{0} = M^{(\alpha+1)} \, L^{(-3\alpha+\beta+2\gamma+1)} \, 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*} M: 0 & = & \alpha+1 \\ L: 0 & = & -3\alpha+\beta+2\gamma+1 \\ T: 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

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

or

(20)   \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

(21)   \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

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

so

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

and

    \begin{eqnarray*} M: 0 & = & \alpha+1 \\ L: 0 & = & -3\alpha+\beta+2\gamma-1 \\ T: 0 & = & -\beta-1 \end{eqnarray*}

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

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

Notice that this derived \Pi product is a non-dimensional parameter that will become infinite as the velocity tends to zero. This behavior is not very convenient in practice, but the \Pi product can simply be inverted to prevent that from happening. The outcomes is still a non-dimensional parameter, but now one that is more convenient, i.e., the \Pi product becomes

(25)   \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

(26)   \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) but generally in an aerodynamic problem a parameter describing a downstream distance should be used, e.g., wing chord c in this case.

This latter non-dimensional 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

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

and so in terms of dimensions, then

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

giving

    \begin{eqnarray*} M: 0 & = & \alpha \\ L: 0 & = & -3\alpha+\beta+2\gamma+1 \\ T: 0 & = & -\beta-1 \end{eqnarray*}

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

(29)   \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 a non-dimensional quantity so

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

This particular non-dimensional 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. The Mach number is an important quantity in aerodynamics that not only affects the magnitude of the forces on a body but also the compressible nature of the flow, in general.

Final Result for Functional Dependency

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

(31)   \begin{equation*} \phi \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

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

or in explicit form

(33)   \begin{equation*} C_{R} = f_{0} \left( Re, M_{\infty} \right) \end{equation*}

The previous result is a significant outcome because it means the original aerodynamic problem in the wind tunnel has been reduced from one that involves five independent variables (for which is practically impossible to determine the functional relationships to R) to one that now only has only 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 non-dimensional parameters Re (the Reynolds number) and M_{\infty} (the Mach number) are known as similarity parameters.

Both the Reynolds number and the Mach number are known to have substantial effects on the aerodynamic forces produced on a wing. In the figure below, the wing lift coefficient (non-dimensional lift) is plotted versus the angle of attack for a series of free-stream Mach numbers from 0.3 (subsonic) up to 0.75 (transonic). Notice that there are two primary effects. The first effect is that increasing the Mach number increases the lift coefficient for a given angle of attack, at least in the lower angle of attack regime. The second effect is that the maximum value of the attainable lift coefficient is reduced with increasing Mach number, which is a consequence of the flow detaching or separating from the wing and is called stall. In other words, the angle of attack at which stall occurs decreases with increasing Mach number. The effects of Reynolds number on the lift on a wing will be discussed later.

Adding another dependency to the problem

Think about what might happen if span b was added as a seventh variable to the parameter list in this previous problem, i.e., in this case assume that

    \[ g \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 first three have already been determined so it should be fairly obvious that the forth non-dimensional grouping will be

    \[ \Pi_{4} = \frac{b}{c} = A\!R \]

which is a span to chord ratio and is known as the aspect ratio of the wing, A\!R. Notice in the figure below, the profound influence of aspect ratio on the lift coefficient. Wings with high aspect ratio are more efficient in that for a given angle of attack and all other parameters held constant then a wing with high aspect will produce more lift and less drag.
Diagram of the effects of wing aspect ratio on the lift coefficient of a finite wing.

Other Worked Examples

The best way of learning the method of dimensional analysis is by doing examples. After working through two or three more examples, the student 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 useful in all branches of engineering and the physical sciences.

Worked Example #1 – Drag on a Tiny Particle

A tiny spherical particle of diameter d (assume in the sub-micron size range) is falling freely vertically at velocity V in air, as shown in the figure below. Find the similarity parameter that governs the drag on the particle.Illustration of a tiny spherical particle falling freely in the vertical direction.

Functional Dependence

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

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

where \mu is the viscosity of the air. The drag can be written is explicit functional form as

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

where “f” is some function to be determined. In implicit form then

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

where “f_1” 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} & d & V & \mu & D \\ \hline \mbox{Mass} ~M: & 0 & 0 & 1 & 1 \\ \mbox{Length} ~L: & 1 & 1 & -1 & 1 \\ \mbox{Time} ~T: & 0 & -1 & -1 & -2 \end{array} \]

\Pi 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] = M^{0}L^{0}T^{0} = (L)^{\alpha} (L T^{-1})^{\beta} (M L^{-1} T^{-1})^{\gamma} M L T^{-2} \]

For the \Pi group to be dimensionless, then

    \begin{eqnarray*} M: 0 & = & \gamma + 1 \\ L: 0 & = & \alpha + \beta - \gamma + 1 \\ T: 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

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

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

which indeed confirms that the grouping is non-dimensional.

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 non-dimensional groupings. The accepted solution to this dilemma is to reduce the number of repeating variables by one and create an extra \Pi grouping.

 

Worked Example #2 – The Bobbing Booster

A solid rocket booster parachutes back to Earth and falls into the sea. The booster bobs around in the sea in an upright position at a frequency of \omega. Using dimensional analysis, show that the non-dimensional frequency of this motion, k, is given by

    \[ k = \omega \sqrt{\frac{d}{g}} \]

It may be assumed that for this problem the frequency of motion depends on the diameter of the booster d, the mass of the booster, m, the density of the water \varrho, and acceleration under gravity, g.

Functional Dependency

By using the information provided then the frequency of motion can be written in explicit functional form as

    \[ \omega = f(D, m, \varrho, g) \]

where “f” is some function to be determined, or in implicit form as

    \[ f_1 (d, m, \varrho, g, \omega) = 0 \]

where “f_1” is some other function. Therefore, N = 5, K = 3 and so there are two \Pi products.

Repeating Variables

The repeating variables must collectively include all the units of mass M, length L and time T, so the best choice here is d, m and g. Notice that if \varrho were to be chosen instead of g as a repeating variable, then collectively the repeating variables would not include time T and so the Buckingham \Pi method will fail in this case. The dependent variable, \omega, cannot be used as a repeating variable.

\Pi Products

Therefore, the two \Pi products are

    \[ \Pi_1 = f_2(d, m, g, \omega) \]

and

    \[ \Pi_2 = f_3(d, m, g, \varrho) \]

Dimensional Matrix

The dimensional matrix is

    \[ \begin{array}{l|r|r|r|r|r} &   d  & m & g & \omega & \varrho \\ \hline \mbox{Mass} ~M: &    0 & 1 & 0 & 0 & 1 \\ \mbox{Length} ~L: &    1 & 0 & 1 & 0 & -3 \\ \mbox{Time} ~T:  &  0 & 0 & -2 & -1 & 0 \end{array} \]

Derivation of the \Pi Products

The first \Pi product is

    \[ \Pi_{1} = (d)^{\alpha} (m)^{\beta} (g)^{\gamma} \omega \]

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

    \[ \left[ \Pi_{1} \right] = M^{0}L^{0}T^{0} = \left( L \right)^{\alpha} \left( M \right)^{\beta} \left( L T^{-2} \right)^{\gamma} \left( 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., in this case

    \begin{eqnarray*} M: 0 & = & \beta  \\ L: 0 & = & \alpha + \gamma  \\ T: 0 & = & -2\gamma - 1 \end{eqnarray*}

Therefore, \alpha = 1/2, \beta = 0, and \gamma = -1/2, so the \Pi_1 product is

    \[ \Pi_1 = (d)^{1/2} (m)^{0} (g)^{-1/2} \omega = \omega \sqrt{\frac{d}{g}} \]

For the second \Pi product then

    \[ \Pi_{2} = (d)^{\alpha} (m)^{\beta} (g)^{\gamma} \varrho \]

In terms of the dimensions then

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

For \Pi_{2} to be dimensionless, then

    \begin{eqnarray*} M: 0 & = & \beta + 1 \\ L: 0 & = & \alpha + \gamma - 3  \\ T: 0 & = & -2\gamma \end{eqnarray*}

Therefore, \beta = -1, \gamma = 0, and \alpha = 3, so the \Pi_2 product is

    \[ \Pi_2 = (d)^{3} (m)^{-1} (g)^{0} \varrho = \frac{\varrho d^3}{m} \]

which is a buoyancy similarity parameter, i.e., the ratio of the mass \varrho d^3 (or weight \varrho d^3  g) of the water displaced to the mass m (or weight m g) of the rocket booster.

Final Non-Dimensional Functionality

Finally, in non-dimensional form then

(34)   \begin{equation*} \Pi_{1} = f_0 (\Pi_2) \end{equation*}

or that

(35)   \begin{equation*} \omega \sqrt{\frac{d}{g}} = f_0 \left( \frac{\varrho d^3}{m}\right) \end{equation*}

CAUTION: Pound force versus pound mass? 

In engineering, we generally always use units of “slugs” (or sometimes just “slug”) for mass in the USC system and never “pounds-mass.” Using “slugs” for mass and “pound” for force prevents the muddling of “pounds-force” lb_f and “pounds mass” lb_m, which, if misused, can result in catastrophic errors in engineering calculations. A pound-force is defined as the product of a pound mass and acceleration under gravity, i.e., 1 lb = 1 lb_m x g, i.e., 1 pound force = 1 pound mass times acceleration under gravity or about 32.17 ft s^{-2}. Therefore, a 1 lb force will accelerate a 32.17 lb mass at 1 ft s^{-2}. However, a 1 lb force will also accelerate 1 slug of mass at 1 ft s^{-2}. This means that a “slug” has a mass equivalent to 32.17 “pounds mass.” So, when we are using “pounds mass,” there is always a multiplicative factor of “g” or “32.17 ft s^{-2}” floating around, and most of the time, it causes nothing but problems in calculations for students. However, the multiplicative factor is “1” when using slugs for the above reason. So, we see immediately why the preference is to keep everything simple and in good order by using “slugs” and never “pounds mass,” and so remove the extra “g” or 32.17 from our engineering calculations.

 

Worked Example #3 – Drag on a Banner

The AIAA Design Build & Fly (DBF) team need to determine the factors that influence the aerodynamic drag on a rectangular banner being towed behind their airplane. The size of the banner is determined by its length, l, and height, h. Use the Buckingham \Pi method to determine the non-dimensional groupings that will govern this problem. You may also assume that the problem is governed by the airspeed of the airplane, as well as the density and viscosity of the air. Assume further that the  banner remains flat and does not flutter in the flow behind the airplane.

Functional Relationship

The relationship between the drag on the banner D and the air properties can be written in the general functional form as

    \[ D = f \left( \varrho_{\infty}, V_{\infty}, \mu_{\infty}, h, l \right) \]

where the size of the banner is represented by its length, l, and height, h. Remember that D is called the dependent variable and \varrho_{\infty}, V_{\infty}, \mu_{\infty}, h and l are are the independent variables. The functional dependence of D in implicit form is

    \[ g \left( \varrho_{\infty}, V_{\infty}, \mu_{\infty}, h, l, D \right) = 0 \]

Counting up the variables gives N = 6 and K = 3 because there are 3 fundamental dimensions in this problem. Therefore, N - K = 3 and so there will be three \Pi products.

For each variable the dimensions are

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

Dimensional Matrix

The dimensional matrix is

    \[ \begin{array}{l|r|r|r|r|r|r} &   D  & \varrho_{\infty} & V_{\infty} & \mu_{\infty} & h &  l  \\ \hline \mbox{Mass: } M  &   1  & 1 & 0 & 1 & 0 & 0 \\ \mbox{Length: } L  &  1  & -3 & 1 &  -1  & 1 & 1 \\ \mbox{Time: } T  &  -2 & 0 & -1 &  -1 & 0 & 0 \end{array} \]

Repeating Variables

Choose \varrho_{\infty}, V_{\infty} and l as the repeating variables, which will all have primary effects on the drag of the banner. These variables also collectively include all the fundamental dimensions of this problem and are linearly independent of each other.

\Pi Products

Following the Buckingham \Pi method, then the three \Pi products are:

    \begin{eqnarray*} \Pi_{1} & = & g_{1} \left( \varrho_{\infty}, V_{\infty}, l, D \right)  \\ \Pi_{2} & = & g_{2} \left( \varrho_{\infty}, V_{\infty}, l, \mu_{\infty} \right) \\ \Pi_{3} & = & g_{3} \left( \varrho_{\infty}, V_{\infty}, l, h \right) \end{eqnarray*}

Derivation of the \Pi Products

For \Pi_{1}:

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

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

    \[ \left[ \Pi_{1} \right] = 1 = M^{0}L^{0}T^{0} = \left( ML^{-3} \right)^{\alpha} \left( LT^{-1} \right)^{\beta} \left( L \right)^{\gamma} \left( MLT^{-2} \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*} M: 0 & = & \alpha+1  \\ L: 0 & = & -3 \alpha +\beta+\gamma+1  \\ T: 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} (l)^{\gamma} D = \varrho_{\infty}^{-1} V_{\infty}^{-2} l^{-2} \, D \]

or

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

i.e., a form of drag coefficient. The aerodynamic force coefficients are usually defined in terms of the dynamic pressure, i.e., (1/2)\varrho_{\infty} V_{\infty}^2 so that more conventionally the force coefficient would be written as

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

It would also be legitimate to write the drag coefficient as

    \[ C_D = \frac{D}{ \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} \, l \, h } = \frac{D}{ \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} \, A } \]

where the banner area A = l \, h is used rather than l^2. Ultimately, how  C_D is defined is just matter of convenience and/or consistency with established convention.

For \Pi_{2}:

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

so

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

and

    \begin{eqnarray*} M: 0 & = & \alpha+1  \\ L: 0 & = & -3\alpha + \beta + \gamma-1  \\ T: 0 & = & -\beta-1 \end{eqnarray*}

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

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

or

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

Inverting the grouping gives

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

which in the latter case is a Reynolds number based on the banner length. Notice that the grouping can be inverted if needed, such as to follow established convention or just because it is otherwise convenient to do so.

For \Pi_{3}:

    \[ \Pi_{3} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (l)^{\gamma} \, h \]

so

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

giving

    \begin{eqnarray*} M: 0 & = & \alpha \\ L: 0 & = & -3\alpha + \beta+\gamma+1  \\ T: 0 & = & -\beta \end{eqnarray*}

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

    \[ \Pi_{3} = (\varrho_{\infty})^{\alpha} (V_{\infty})^{\beta} (l)^{\gamma} \, h = \varrho_{\infty}^{0} V_{\infty}^{0} l^{-1} \, h = \frac{h}{l} \]

so

    \[ \Pi_{3} =  \frac{h}{l} \]

or again inverting this grouping (for convenience) giving

    \[ \Pi_{3} =  \frac{l}{h} = A \! R \]

which is a length to height ratio or what would be called an aspect ratio A\!R.

Final Non-Dimensional Functional Dependency

As a result of our dimensional analysis then

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

or

    \[ \phi \left( C_{D}, Re, A\!R\right) = 0 \]

Finally, in explicit form the drag coefficient can be written as a function of Reynolds number based on banner length and the aspect ratio of the banner, i.e.,

    \[ C_{D} = f_{0} \left( Re, A\!R \right) \]

Note: Try this problem again using \varrho_{\infty}, \mu_{\infty} and h as the repeating variables. What happens to the groupings?

Other Dimensionless Quantities

Some physical variables are naturally dimensionless by virtue of 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, is often used, which is a dimensionless number because it is defined as the ratio of the density of a given liquid relative to the density of water.

Finally, it should be remembered that all angles are dimensionless. In engineering analysis, angles that are measured in radians are always used, which is defined as a ratio of the arc length to the radius of a circle. 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).

Conversion Factors

A ratio that converts one unit of measurement into another unit of measurement is called a conversion factor. Engineers often use conversion factors from common measurement units to engineering units, an example being converting from the temperature in Celsius to engineering temperature units of Rankine.

Often pressure may be measured in units of “pounds per square inch” or lb/in^2 or psi, but this is not an engineering (base) unit. It needs to be converted to “pounds per square foot” or lb/ft^2 or psf, i.e., 1 psf = 1/144 psf where the conversion factor, in this case, is 1/144 the factor 144 being 12 \times 12 where 12 inches equals one foot. Other units of pressure measurement include a “bar” where 1 bar = 100 kPa = 100,000 Pa with Pa (1 Pascal = 1 Newton per square meter or 1 Pa = 1 N/m^2) being the base engineering unit of pressure measurement. A Pascal is a small amount of pressure, so the use of hectoPascals (hPa), which is 100 Pa, kilo-Pa (kPa) or even Mega-Pa (MPa) is not uncommon. So, to obtain pressure units of Pa from measurements of pressure in bars, then 1 Pa = 1/100,000 bars, so the conversion factor, in this case, is 1/100,000 or 10^{-5}.

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 non-dimensional groupings called similarity parameters. The strength of dimensional analysis as a tool usually becomes apparent when it is 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 non-dimensional groupings helps understand what parameters are essential as they will affect the physical behavior and what parameters are not.

The dimensional analysis must be used in a systematic and structured approach, which is 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 adhered to 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 non-dimensional (unitless) parameters that should not only help in the understanding of 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 non-dimensional 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) as well as 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 the hull of a ship, in part, depends on the Froude number. The Froude number is defined as F\!r = 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 being dragged behind a ship. You propose to do this by doing small scale experiments in a water tunnel. Determine the non-dimensional 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 in the methods of dimensional analysis, navigate to the following online resources: