Dimensional Analysis (Partial Analysis)

Dimensional Analysis (Partial Analysis) DA is a mathematical method of considerable value to problems in science and engineering especially physics and fluid mechanics. All physical quantities can usually be expressed in terms of certain primary quantities which in mechanics are: Length (L), Mass (M), and Time (T). e.g. Force = Mass x Acceleration = Mass x Length/Time = M L T -

The principle of homogeneity of dimensions can be used to:. Check whether an equation has been correctly formed (apples = apples) ;. Establish the form of an equation relating a number of variables (dimensional analysis); and 3. Assist in the analysis of experimental results (experimental design). Checking Equations: Eg. : Show by D. A. that the equation: p + ½ ρ v + ρ g z = H

is a possible relationship between pressure p, velocity v, and height z for a frictionless flow along a streamline of a fluid of density ρ, and determine the dimension of H. Solution: p = Force/Area = m a/area = M L T - L - = M L - T - ρ = mass/volume = M L -3 v = length/time = L T - g = L T -, Z = L p = M L- T- ½ ρ v = M L -3 x L T - = M L - T - ρ g z = M L -3 x L T - x L = M L - T - 3

All terms on the l.h.s have the same dimensions and the equation is physically possible if H also has dimension of M L - T -. ½ is a pure number hence not amenable to D. A. Constant could have been another number and the equation is still valid. Therefore, the analysis is not complete. We need additional information to confirm that the constant is indeed ½. Eg.. Check that Einstein s famous equation is indeed dimensionally correct: E = mc 4

Methods of dimensional analysis: grouping of variables. Many methods are available but we will consider the main ones:. Rayleigh indicial method. Buckingham π 3. Hunsaker and Rightmire (Whittington s) 4. Matrix Method The methods above are given roughly in chronological order. Also from most tedious to use to almost automatic (if matrix inversion routine available). 5

Rayleigh s method: Example: The velocity of propagation of a pressure wave through a liquid can be expected to depend on the elasticity of the liquid represented by the bulk modulus K, and its mass density ρ. Establish by D. A. the form of the possible relationship. Assume: u = C K a ρ b U = velocity = L T -, ρ = M L -3, K = M L - T - L T - = M a L -a T -a x M b L -3b M: 0 = a + b L: = -a 3b T: - = - a 6

Therefore: a = ½, b = -½, and a possible equation is: u = C K ρ Rayleigh s method is not always so straightforward. Consider the situation of flow over a U-notched weir. Q = f(ρ, µ, H, g) [Q] = [C ρ a µ b H c g d ] [ ] => dimensions of Using the M, L, T system, [L 3 T - ] = [ML -3 ] a [M L - T - ] b [L] c [L T - ] d M: 0 = a + b () L: 3 = -3a b +c + d () T: - = - b d (3) 7

8 We have only 3 equations, but there are 4 unknowns. Need to express a, b, c, in terms of d. b = d a = -b = d - c = 3 + 3a + b d = + 3d Q = C ρ (d-) µ (-d) H (+3d) g (d) d g H H C = 3 µ ρ ρ µ = 3 µ ρ φ µ ρ g H H Q dimensionless groups. Please check. This method can be very tedious if there are more variables.

Buckingham s π theorem This is the most well-known method and is the basis for the matrix and other methods. If there are n variables in a problem and these variables contain m primary dimensions (e.g. M, L, T), the equation relating the variables will contain n-m dimensionless groups. Buckingham referred to these dimensionless groups as π, π,, π n-m, and the final equation obtained is: π = φ (π, π 3,, π n-m ). 9

E.g. consider the previous problem with 5 variables, hence there will be 5 3 = dimensionless groups. The dimensionless groups can be formed as follows: ) Select 3 repeating variables from the list of variables which together must contain M, L,T. eg. ρ (ML 3 ), H (L), g (LT - ) ) Combine repeaters with one other variable in turn, 3) Never pick the dependent variable as a repeater. Consider: Q = f(ρ, µ, H, g) 5 variables Pick repeaters: ρ, H, g 0

Combine with one other variable, say, Q. x y z ρ H g π = Q ( 3 ) x ( ) ( ) ( 3 ) y z ML L LT L = T M: x = 0 L: -3x + y + z + 3 = 0 T: -z = 0 Therefore: z = -½, y = -5/ 0 π = ρ H π = g / Q H 5/ 5/ g / Q

π = ρ H g µ x y z ( 3 ) x y ( ) ( ) ( ) z ML L LT ML = T M: x + = 0 L: -3x + y + z - = 0 T: -z = 0 x = -, z = -/, y = 3/ π π = 3/ = ρ H g ρh µ 3/ g / / µ That is: g / Q H 5/ = φ ρh µ 3/ g /

Ex. The head loss per unit length ( h/l) of pipe in turbulent flow through a smooth pipe depends on v, D, ρ, g, and µ. Use Buckingham s π theorem to determine the general form of the equation. F( h/l, v, D, ρ, µ, g) = 0; 6 variables, 3 dimensions, 3 π terms. That is: π = φ(π, π 3 ). Choose 3 repeating variables: v, D, and ρ π = v ρ µ x D x x ( ) x x ( ) ( 3) x LT L ML ML = T 3

M: z + = 0 L: x + y 3z = 0 T: -x = 0 Therefore: x = -, z = -, and y = -. π = v D ρ µ or µ ρvd or ρvd µ = Re h L Final solution: = φ{ Re, Fr } Where Re = Reynold s number, and Fr = Froude number = v gd 4

Hunsaker and Rightmire s Method (947) This method is easier to use and quicker than Buckingham s method. The method is similar the Buckingham s method by use of repeating variables, but express them in term of the variables themselves. e.g. D to represent Length, [L] = D v to represent Time, [T] = L/v = D/v r to represent Mass, [M] = ρl 3 = ρd 3 π : g = L T - = v D D - = v D - v, therefore π = gd 5

π : µ = M L - T - = Therefore: π = ρvd µ ρd D = 3 v Re = ρvd h π 3 =, as before. L Example: A spherical drop of liquid of diameter D oscillates under the influence of its surface tension. Investigate the frequency of oscillation f. F(f,σ, D, ρ) = 0 Ans: f = K 3 σ ρd 6

Matrix Method Best for problems with many variables. Can be solved using a matrix inversion routine. Consider a problem with 7 variables in 4 dimensions (M, L, T, θ) First form the dimensional matrix: M L T θ a b c d e f g A B 7

We need to transform the above matrix to: a b c d e f g a b c d 0 0 0 0 0 0 0 0 0 0 0 0 I = A - A D = A - B A - = adj A/ A Whatever operation was done to get the unit matrix on the left must be also done on the right to get D. 8

Eg. Consider the earlier problem of the head loss per unit length. F( h/l, v, g, D, ρ, µ) = 0 M ρ v D g µ h/l 0 0 0 0 L T -3 0-0 - 0 - - 0 A B Inverse of A can be obtained using Minitab, Matlab, or by hand. 9

0 = 3 0 0 0 0 A D = 0 0 0 0 3 0 0 0 0 D = 0 0 0 0 The three π terms are: L h vd v gd,, ρ µ The matrix method can obtain the π terms all in one go instead of one term at a time.

Comments about Dimensional Analysis Most important but most difficult problem in applying DA to any problem is the selection of the variables involved. There is no easy way of identifying the correct variables without specialized knowledge about the phenomenon being investigated. If you select too many variables, you get too many π terms and may require much additional experimentation to eliminate them. If important variables are omitted, then an incorrect result will be obtained, which may prove to be costly and difficult to ascertain.

Type of variables:. Geometry length, angles, diameter, or area.. Material properties ρ, µ, elasticity, etc. 3. External effects any variable that tends to produce a change in the system e.g. forces, pressures, velocities, gravity, etc. You need to keep the number of variables to a minimum, and that they are independent. E.g. D and Area need not be included together because one is derived from the other. Therefore, heavy thinking is required in variable selection similar to DOE.

Points to consider in the selection of variables:. Clearly define the problem. What is the main response variable of interest? That is, what is Y?. Consider the basic laws that govern the phenomenon. Even a crude theory may be useful. 3. Start the selection process by grouping the variables in the 3 broad classes: geometry, material properties, and external effects. 4. Consider other variables that may not fall into one of the three categories, e.g. time, temperature, colour, equipment, etc. 3

5. Be sure to include all quantities that enter the problem even though some of them may be held constant e.g. g. For D.A. it is the dimensions of the quantities that are important not specific values. 6. Make sure that all variables are independent look for relationships among subsets of the variables (same as DOE). Remember that after a dimensional analysis, you still need to carry out the experiment to relate the dimensionless groups. Hence DOE may be needed unless you have only or π terms. More on this later. 4

More on π terms Specific π terms obtained depend on the somewhat arbitrary selection of repeating variables. For example, if we choose: µ, D, g instead of ρ, D, v, we would end up with a different set of π terms. Both results would be correct, and both would lead to the same final equation for h/l, however, the function relating the different π terms would be different. 5

Hence, there is not a unique set of π terms which arises from a dimensional analysis. However the required number of π terms is fixed, and once a correct set is determined all other possible sets can be developed from this set by combination of products of the original set. e.g. if we have problem involving 3 π terms, π = φ(π, π 3 ) we can combine the π terms to give a new π term: ' a b π π 3 π = Then the relationship could be expressed as: ( ' π = φ π ) Or even: ( ' π = φ π ),π3,π3 6

All these would be correct; however, the required number of π terms cannot be reduced by this manipulation; only the form of the π terms is altered. e.g. π pd = ; π ρv = ρvd µ π = π π = pd ρv ρvd µ = pd µ v Which form of p terms is best? There is no simple answer. Best to keep it simple. Some p terms are well-known dimensionless numbers like Reynolds, Froude, Mach, Weber, Cauchy, Euler, etc. 7

Correlation of Experimental Data One of the most important uses of dimensional analysis is an aid in the efficient handling, interpretation, and correlation of experimental data. As noted earlier, DA cannot provide a complete answer and a specific relationship among the π terms cannot be determined without experimentation. The degree of difficulty obviously depends on the number of π terms. 8

Problems with: π term π = C where C = a constant π terms π = f(π ) simple regression problem > π terms π = f(π, π 3 ) multiple regression problem With more and more π terms, a DOE approach may be needed as in Problem 4 of assignment and may require the use of RSM if relationship is nonlinear. 9

Eg : The pressure drop per unit length, p/l for the flow of blood through a horizontal small diameter tube is a function of flow rate Q, diameter D, and the blood viscosity µ. For a series of test with D = mm and µ = 0.004 Ns/m, the following data were obtained for p measured over a length of 300 mm. Q (m 3 /s) 3.6 x 0-6 4.9 x 0-6 6.3 x 0-6 7.9 x 0-6 9.8 x 0-6 p (N/m ). x 0 4.5 x 0 4.9 x 0 4.4 x 0 4 3.0 x 0 4 Perform a DA for this problem and make use of the data to determine a general relationship between p and Q, one that is valid for other values of D, L and µ. Solution: 4 variables, F ( p/l, D, Q, µ) = 0, i.e. 4 3 = π term 30

From D.A. (try this yourself), we get: ( p / L) D µ Q 4 = constant Substituting the values used in the experiment, ( p / L) D µ Q 4 = p / 0.3(0.00) 0.004 Q 4 =.33 0 8 p Q. Using the data obtained from the experiment, ( p / L) D µ Q 4 = [ 40.6, 40.7, 40. 40.4 40.7] p µ Q 40 L D Average for constant C = 40.5, hence: =.5 4 3

Eg : A liquid flows with a velocity v through a hole in the side of a tank. Assume that v = f(h, g, ρ, σ). Where h is the height of water above the hole, ρ is the density of the fluid, and σ the surface tension. The density ρ is 000 kg/m 3, and s = 0.074 N/m. The data obtained by changing h and measuring v are: V (m/s) 3.3 4.43 5.4 6.5 7.00 h (m) 0.50.00.50.00.50 Plot the data using appropriate dimensional variables. Could any of the original variables have been omitted? Solution: 5 variables, F(v, h, g, ρ, σ) = 0, π terms. From dimensional analysis, v gh ρgh = φ σ 3

ρgh /σ 3.3 x 0 4 3.3 x 0 4 8.8 x 0 4 53.0 x 0 4 8.9 x 0 4 v/(gh) /.4.4.4.4.4 Plotting the data will show that v gh means that ρ and σ can be omitted. is independent of ρ gh, which σ Of course this is well-known if one were to apply the Bernoulli equation to solve the problem. 33

References: Chapter, 3, 4 notes from course website. Thomas Szirtes (998): Applied Dimensional Analysis and Modeling, McGraw Hill, 790 pages. Most books in physics and fluid mechanics. 34