Showing posts with label dimensional analysis. Show all posts
Showing posts with label dimensional analysis. Show all posts

Saturday, March 21, 2020

Dimensional Analysis Revisited


The use of dimensions or units in physics may be used to check the homogeneity of physical equations, or obtain a useful empirical equation from basic measurements, observations.

Example:

Consider a volume V of a liquid flowing per second through a pipe under steady pressure as shown below:
Fig. 1: The fluid flow through a uniform pipe


It is reasonable to assume, in the first instance, that V (volume flow per unit time) depends upon:

a)     The coefficient of viscosity, h, of  the liquid
b)     The radius, r, of the pipe
c)      DP/L the pressure gradient arising from the flo

Then we may write: V = k hx ry (DP/L)z


Where k denotes a constant, and x, y and z are indices whose values are to be found.

We know:

V has dimensions L3 T - 1

h has dimensions M L-1 T -1

r has dimension L

DP has dimensions M L-1 T - 2

So, DP/L has dimensions: M L-2 T - 2

Then, equating dimensions:

L3 T – 1  = [M L-1 T -1]x [L]y [M L-2 T - 2]z

We next equate indices for M, L, and T on both sides:

For M:  0 = x + z

For L: z = -x + y – 2z

For T: -1 = -x – 2z

Next, solve for the indices:

a)     x = -z
b)     -1 = -(-z)- 2z = -z or z = 1
c)      -1 = -x -2z
d)     3 = 1 + y -2 or y = 4

Finally:  V = k h- 1 r 4 (DP/L)1

Or: V =  k DP r 4 /h L


We note here that k cannot be obtained by the method of dimensions, but a fuller analysis would reveal k = p/8, so:

V = p DP r 4 /8 h L    (Poiselle’s formula)

Problems:

1)      Use the method of dimensions to obtain a formula for the force experienced by a sphere of radius r moving at velocity v through a fluid of viscosity h.  (Hint: Set out the force as F = k hx v y r )

2) Repeat the exercise above to obtain the period (T) for a pendulum’s swing.


  1. Errors & Uncertainties in Measurements

In general, there are two types of errors with which we will deal in physics:

1) Systematic: an error repeated due to a fault in the instrument of the procedure.

2)     Random: an error which arises due to intrusion of some random factor, or the necessity to estimate and the limits of instrumental sensitivity.

An illustration of a systematic error is illustrated below for an ammeter.  We observe that in this case, the meter is incorrectly calibrated since the real zero is displaced from the value read by the ammeter, by and amount 0.2 A. Thus when the meter reads 1.2 A it is actually 1.0 A.


Fig. 2: Improperly calibrated ammeter.

In this case, X = 1.0 A represents the “true value” for the measurement, and the reading (x) is registered as:

x = X + e

Then:  + e  =  x – X = 1.2A – 1.0A =  + 0.2A

The nature of a random error is illustrated in Fig. 3.:



Fig. 3: A random error arising from multiple deviating readings.

In this case, reading of the ammeter yields a series of slightly different readings at different times, say: x1, x2, x3, and x4. In this case, x denotes the “true value” and we have:

x1 = x – c1

x2 = x – c2

x3 = x – c3

where c1, c2, and c3 are “correction factors” of the observed value to the true value, x. In such an instance, we generally reduce the effects of these random different values by taking the average value x by:

x = (x1 + x2 + x3 +…..xn)/ n

Which should be closer to the true value x, than any one of the other observed values.

In general the associated error will be:

e = (x – x )

So the fractional error is:  f = e/ x

And the percentage error is: p = (e/x) x 100%

Usually: f < < 1 and p < < 100%

Example:

As illustrated below, a measurement of 2.5 cm is made using a scale rule graduated in millimeters (mm or 0.1 cm) What would the possible error (p.e.) be in this case, and how would we write the length of the object?



Fig. 4.: Example illustrating possible error of measurement for an object

The possible error in this case is the equivalent of the limit of accuracy for the scale or 0.1 cm. Then:

p.e. =  + 0.1 cm

The length of the object is written: L = 2.5 + 0.1 cm

The percentage possible error (p.p.e.) is:

(+ 0.1 cm/ 2.5 cm) x 100% =  + 4%

Example (2):

For the same measurement a vernier caliper may find:

L = 2.46 + 0.01 cm

Since its accuracy is a factor 10 enhanced over the scale.

The possible error is then: p.e. = + 0.01 cm

And:  the percentage possible error is:

p.p.e. = (+ 0.01 cm/ 2.46 cm) x 100% =  + 0.4%


 3. Combining Errors of Measurements

It is possible and indeed, frequently necessary, to combine errors. We consider several cases below:

Case (1): The error sum Q = a + b

Then: The total p.e. in Q = p.e. in a + p.e. in b

(Also, for absolute errors: DQ = Da + D b)

Let: a = 5.1 + 0.1 cm  and b = 3.2 + 0.1 cm 

Then Q = (5.1 + 0.1 cm ) + (3.2 + 0.1 cm ) = 8.3 + 0.2 cm 

In addition:

Q max = (5.2 + 3.3) cm = 8.5 cm,   and

Q min = (5.0 + 3.1) cm = 8.1 cm

Case (2): The error difference Q = a – b

If Q = a – b then the same rule applies, in terms of addition of possible errors. Also: Da =  D b  = + e)

Case (3) Product or quotients and fractional errors.

If a quantity x is measured and recorded as x’ units, then let: (x’ – x) denote the error in x or x’ = x + e. But we may also write:

x’ = x(1 + f)

where f = e/x is the fractional error in x. Note also that: (x’/ x) = 1 + f.

For a sum or difference: Q = a + b = (a + e1) + (b + e2)

And: f = (e1 + e2)/ (a + b)

Also:

f1 = e1/ a and f2 = e2/b

Therefore: f = f1 + f2

Or f = (af1/ a) + bf2/b)

Error in a product, viz, Q = ab

Then: Q = a (1 + f1) x b (1 + f2)

Or: Q = ab(1 + f1 + f2)

Error in a quotient, viz. Q = a/b:

Then: Q = a(1 + f1)/ b(1 + f2)

Q = a/b [ (1 + f1)(1 – f2 + …) = a/b [1 + f1 – f2]


Problems for the ambitious reader :

1) Two lengths ℓ1 and ℓ2 are measured as (10.0 + 0.1 cm ) and (9.0+ 0.1 cm ).  Find the p.e. for each length. Express the sum (ℓ1 + ℓ2) with errors-uncertainties. Thence or otherwise obtain the fractional error derived from the sum.

2) The length of a sheet of paper is measured as: ℓ = 11.0 + 0.1 cm, and its width as w =  8.5 + 0.1 cm. Find the Area A = ℓ x w in terms of the associated possible errors. Find the fractional error associated with the area.

3) A glass thermometer is scaled in tenths of a degree. If a student reads the lab temperature as 30.5 o C, express the possible error in the measurement.

Thursday, March 19, 2015

An Introduction To Dimensional Analysis

In the post before last I cited a loony question received (at 'All Experts') from a clueless individual who sincerely believed her "pure light speed math" held the keys to the universe. She presented a link to a Canadian website 'Einstein now', e.g. http://members.shaw.ca/warmbeach/INDEX3.htm

Which confirmed that all the units were wacky and the math actually improper. Moreover, none of the equations offered were complete which she defended in a subsequent comment. Of course, when a site presents incomplete equations we know we've ventured into crank territory, or should. And cranks often try to gain respect by pressing legit sites to bestow some gravitas - but we never do. What we do is call them "rubbish".

The young lady never would have been trapped in such bunkum had she possessed even rudimentary skills at dimensional analysis (the basis for all good physics) and also had an acquaintance with proper units as well as significant figures.  In this post we examine some examples of dimensional analysis which can show whether an equation that purports to show a force, velocity or current, really does. This ought to be in the "kit bag" of any person who claims to do physics, "pure light speed" (sic) or not.

Units and dimensional analysis lay the groundwork for physics, since above all physics is a science of measurement and obtaining correct values while recognizing the associated errors is essential. We begin with the base  and derived units for the S.-I. system given in the link below:

http://physics.nist.gov/cuu/Units/units.html

The units are all metric which is used because it has the greatest functionality, convenience as well as rationality, as opposed to English units.  Let's look at a few basic examples (the readers should confirm familiarity with the assorted unites used from the table first):

1) Current as function of drift velocity:

By definition  : I = nevA

where v is the drift velocity, n is the density of electrons per unit volume, A is the area of the conductor, and e, the electronic charge.

We know v will be in m/s (meters per second) and assume here: 8.4 x 10 -4  m/s

e = 1.6 x 10 -19 C

Area A is: 2 x 10  -5 m 2


The number density is: n =  7.4 x 10 28 /m3

If then all the above factors are multiplied together they ought to give the current in amperes (A) or Coulombs (charge) per second.

We write:

I  = n e v A =

( 7.4 x 10 28 /m3)(1.6 x 10 -19 C) ( 8.4 x 10 -4  m/s) (2 x 10  -5 m )

Note when multiplying the factors together one must also attend to units and how they combine on multiplication.  Paying attention to the units factors alone we will see:

(m3) / (m3)  C/ s

Since the cubic meters cancel out that leaves only C/s or coulombs per second which we note from the table is A (amperes).

Thus, I  is indeed confirmed as a current.

2) The wavelength of the photon emitted from an atom, going from higher energy E2 to lower E1

We know this should be in m.


The equation to test is:  l =   hc/ (E2 – E1) 

Where h = 6.626069 x 10- 34 J-s   (The Planck constant)

c = 3 x 10 8 m/s

E2 - E1 =   16.4 x 10-19 J 

Then:

l =    (6.626069 x 10- 34 J-s)(3 x 10 8 m/s)/ 16.4 x 10-19 J 

Check first the unit product - that it is really in m:

(J-s) (m/s) / J    =   (J/ J ) (s/s ) m  = m

so the end unit is m, making the dimensional character correct.


3) The magnetic Reynolds number is supposed to be dimensionless

 It is defined:   Âm  =  L VA / h


The units for  magnetic diffusivity,  h   are m2 /s

 V is in m/s

L is in m

Check the end units, which ought to be NO units:

(m)    (m/s)/  (m2 /s)

So:   (m2 /s) /   (m2 /s)  =   1

 We see m2 cancels with  m2  and s with s leaving a dimensionless result.

There is also more detailed dimensional analysis which focuses on treating the dimensions of the equation as opposed to units. Thus, this approach may be used to check the homogeneity of physical equations, or obtain a useful empirical equation from basic measurements.


Consider a volume V of a liquid flowing per second through a pipe under steady pressure as shown below:




It is reasonable to assume, in the first instance, to assume that V (volume flow per unit time) depends upon:


a)     The coefficient of viscosity, h, of  the liquid


b)     The radius, r, of the pipe


c)      DP/L the pressure gradient arising from the flow


Then we may write as a first approximation to the actual equation: V = k hx ry (DP/L)z

Where k denotes a constant, and x, y and z are indices whose values are to be found.

We know:

V has dimensions L3 T - 1



h has dimensions M L-1 T -1



r has dimension L



DP has dimensions M L-1 T - 2



So, DP/L has dimensions: M L-2 T - 2



Then, equating dimensions:



L3 T – 1  = [M L-1 T -1]x [L]y [M L-2 T - 2]z


We next equate indices for M, L, and T on both sides:


For M:  0 = x + z

For L: z = -x + y – 2z

For T: -1 = -x – 2z

Next, solve for each of the indices:
a)     x = -z
b)     -1 = -(-z)- 2z = -z or z = 1
c)      -1 = -x -2z
d)     3 = 1 + y -2 or y = 4

Finally:  V = k h- 1 r 4 (DP/L)1

Or: V =  k DP r 4 /h L

We note here that k cannot be obtained by the method of dimensions, but a fuller analysis would reveal k = p/8, so:

V = p DP r 4 /8 h L    (Poiselle’s formula)

 Problems:
1)     Use the method of dimensions to obtain a formula for the force experienced by a sphere of radius r moving at velocity v through a fluid of viscosity h.  (Hint: Set out the force as             F = k hx v y r z  )

2) Repeat the exercise above to obtain the equation  for the period (T) of a pendulum’s swing. (Hint: The pendulum's period T should depend on its length L and the acceleration of gravity, g.)