17.11.2012 Views

Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian

Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian

Dimensional Analysis – Drag on a Sphere R. Shankar Subramanian

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<str<strong>on</strong>g>Dimensi<strong>on</strong>al</str<strong>on</strong>g> <str<strong>on</strong>g>Analysis</str<strong>on</strong>g> – <str<strong>on</strong>g>Drag</str<strong>on</strong>g> <strong>on</strong> a <strong>Sphere</strong><br />

R. <strong>Shankar</strong> <strong>Subramanian</strong><br />

In a typical experiment, we look for how a dependent parameter varies as we change a variety of<br />

independent parameters. By independent parameters, we mean those that we can vary<br />

c<strong>on</strong>veniently, and the dependent parameter is the <strong>on</strong>e whose behavior we seek to establish.<br />

<str<strong>on</strong>g>Dimensi<strong>on</strong>al</str<strong>on</strong>g> analysis permits us to organize the process by which we vary the independent<br />

parameters. In fact, it helps us identify the true dependent and independent parameters in a<br />

situati<strong>on</strong>; in the process, the number of parameters that we must c<strong>on</strong>sider is minimized. But<br />

dimensi<strong>on</strong>al analysis is not foolproof – we must be careful in c<strong>on</strong>sidering all possible parameters<br />

that can affect the dependent parameter. If we omit a crucial parameter in making the list of<br />

independent parameters, dimensi<strong>on</strong>al analysis cannot help us find it. We always need comm<strong>on</strong><br />

sense and physical intuiti<strong>on</strong> in selecting the lists of parameters in a problem.<br />

Here, by c<strong>on</strong>sidering a simple example, I’ll show you how to use dimensi<strong>on</strong>al analysis. We aim<br />

to develop an organized way of examining how the drag <strong>on</strong> a sphere settling through a fluid<br />

varies with the relevant parameters. The first step is to identify the parameters <strong>on</strong> which the drag<br />

F D is likely to depend. These are the diameter D of the sphere, the density ρ and the viscosity<br />

μ of the fluid, and the settling velocity of the sphere V .<br />

Step 1<br />

Make a list of parameters and identify their dimensi<strong>on</strong>s using the fundamental dimensi<strong>on</strong>s of<br />

Mass ( M ) Length ( L ) and Time ( T ) . From now <strong>on</strong>, we’ll use the term “variables” to<br />

designate these parameters, but recognize that they are parameters in the usual sense of that term.<br />

Variable Symbol Dimensi<strong>on</strong><br />

<strong>Sphere</strong> Diameter<br />

D L<br />

Fluid Density ρ<br />

M<br />

3<br />

L<br />

Settling Velocity<br />

V<br />

L<br />

T<br />

Fluid Viscosity μ<br />

M<br />

LT<br />

<str<strong>on</strong>g>Drag</str<strong>on</strong>g> D F ML<br />

2<br />

T<br />

1


Step 2<br />

Identify the repeat variables. There should be as many repeat variables as there are dimensi<strong>on</strong>s.<br />

These are variables that will, in principle, appear in every dimensi<strong>on</strong>less group that we form.<br />

The requirement is that it should not be possible to form a dimensi<strong>on</strong>less group in the form of a<br />

product of powers of these variables. The simplest way to meet this requirement is to select the<br />

repeat variables so that each has a unique dimensi<strong>on</strong> in it. Also, it is helpful to choose repeat<br />

variables that are as simple as possible in their dimensi<strong>on</strong>s. With this in mind, we choose the<br />

first three, namely, the diameter, density, and velocity as our repeat variables, because there are<br />

three dimensi<strong>on</strong>s in this problem. You can see that density has Mass ( M ) appearing uniquely in<br />

it, and velocity has time ( T ) uniquely appearing in it, and the diameter has neither. Therefore, it<br />

is not possible to use these three variables to form a dimensi<strong>on</strong>less group.<br />

Step 3<br />

The number of dimensi<strong>on</strong>less groups is always equal to the number of variables minus the<br />

number of repeat variables. Therefore, we can expect to form two dimensi<strong>on</strong>less groups in this<br />

problem. The group involving the drag will be the dependent dimensi<strong>on</strong>less group and that<br />

involving the viscosity will be the independent dimensi<strong>on</strong>less group. Each group is obtained by<br />

forming a product of each repeat variable to an unknown power and then multiplying by <strong>on</strong>e of<br />

the remaining variables. If we call these two dimensi<strong>on</strong>less groups Π 1 and Π2<br />

, then we might<br />

define Π as follows.<br />

Π =<br />

1<br />

1<br />

a b c<br />

μ D ρ V<br />

To find the values of ab , ,andcthat<br />

would make this group dimensi<strong>on</strong>less, first we write out the<br />

equati<strong>on</strong> in terms of the dimensi<strong>on</strong>s of each side.<br />

b c<br />

0 0 0 M a ⎛M ⎞ ⎛ L⎞<br />

1+ b −+ 1 a− 3b+ c −− 1 c<br />

= × × × =<br />

3<br />

M LT L ⎜ ⎟ ⎜ ⎟ M L T<br />

LT ⎝ L ⎠ ⎝T⎠ Now, match the powers of M , L, and T<br />

the exp<strong>on</strong>ents ab , , and cas<br />

follows.<br />

1+ b = 0<br />

− 1+ a− 3b+<br />

c=<br />

0<br />

−1− c = 0<br />

<strong>on</strong> both sides. This yields three algebraic equati<strong>on</strong>s for<br />

The soluti<strong>on</strong> of these equati<strong>on</strong>s is straightforward, and the result is a = b= c=−<br />

1.<br />

Therefore,<br />

μ<br />

Π 1 =<br />

DρV 2


You can verify that the group in the right side of the above result is indeed dimensi<strong>on</strong>less. Now,<br />

there is nothing in dimensi<strong>on</strong>al analysis that prevents us from using powers of this group as a<br />

dimensi<strong>on</strong>less group in which these four variables appear. It is c<strong>on</strong>venti<strong>on</strong>al to invert this group<br />

and call it the Reynolds number.<br />

So, let us invent a new dimensi<strong>on</strong>less group to replace Π 1 .<br />

Reynolds Number Re =<br />

D V ρ<br />

μ<br />

We shall have more to say about the Reynolds number and its physical significance later in the<br />

course. Next, let us find the dimensi<strong>on</strong>less group that includes the dependent variable F D .<br />

Π = F D ρ V<br />

2<br />

D<br />

a b c<br />

Again, we write out the equati<strong>on</strong> in terms of the dimensi<strong>on</strong>s of each side.<br />

b c<br />

0 0 0 ML a ⎛M⎞ ⎛L⎞ 1+ b 1+ a− 3b+ c −2−c = × × × =<br />

2 3<br />

M LT L ⎜ ⎟ ⎜ ⎟ M L T<br />

T ⎝ L ⎠ ⎝T ⎠<br />

As before, matching the powers of M , L,andT the exp<strong>on</strong>ents ab , ,andc.<br />

1+ b = 0<br />

1+ a− 3b+<br />

c=<br />

0<br />

−2− c = 0<br />

The soluti<strong>on</strong> is a=− 2, b=− 1, c=−2.<br />

Therefore,<br />

Π =<br />

F<br />

DV ρ<br />

2 2<br />

D<br />

2<br />

<strong>on</strong> both sides yields three algebraic equati<strong>on</strong>s for<br />

In practice, a dimensi<strong>on</strong>less drag known as the <str<strong>on</strong>g>Drag</str<strong>on</strong>g> Coefficient, C D , is defined as follows.<br />

C<br />

D<br />

=<br />

8FD DV<br />

2 2<br />

π ρ<br />

The reas<strong>on</strong> for the multiplicative factor ( 8/π ) is that the drag coefficient is defined as the drag<br />

divided by the product of the projected area of the sphere and the velocity head. Thus,<br />

3


C<br />

D<br />

FD 8FD<br />

= = 2 2<br />

π 2 1 2<br />

D × ρV<br />

π DV ρ<br />

4 2<br />

Therefore, dimensi<strong>on</strong>al analysis tells us that the drag coefficient is a universal functi<strong>on</strong> of<br />

Reynolds number, regardless of the choice of fluid, sphere diameter or the settling velocity.<br />

C = φ<br />

D<br />

( Re)<br />

The nature of the functi<strong>on</strong> φ has to be established by experiments or theory as appropriate. The<br />

universality means that when we plot the drag coefficient against the Reynolds number, we<br />

obtain a single curve, regardless of the choice of sphere diameter, fluid, and the velocity at which<br />

the sphere moves through the fluid. This is a remarkable result, because it permits us to use any<br />

fluid(s) of our choice to carry out the experiments to find the universal relati<strong>on</strong>ship between<br />

these two dimensi<strong>on</strong>less variables, and <strong>on</strong>ce determined, use it to infer the drag <strong>on</strong> any sphere<br />

settling through any fluid at any velocity. You can see how powerful a tool this is.<br />

Figure 7.16 (b) in Fluid Mechanics, Sixth Editi<strong>on</strong>, by F.M. White shows the universal curve that<br />

results for the drag <strong>on</strong> a sphere as a functi<strong>on</strong> of the Reynolds number.<br />

4

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!