Outline and Equation Sheet for M E 345

Outline and Equation Sheet for M E 345
Author: John M. Cimbala, Penn State University
Latest revision, 13 November 2014
Introduction
 Primary dimensions – mass, length, time, temperature, current, amount of light, and amount of matter.
 Significant digits – the rules for multiplication and division, and for addition and subtraction.
 Rounding off – round up if the least significant digit is odd, and truncate if the least significant digit is even.
Dimensional Analysis
 Law of dimensional homogeneity – Every additive term in an equation must have the same dimensions.
 The method of repeating variables – there are 6 steps:
1. List the parameters and count them, n.
2. List the primary dimensions of each parameter.
3. Set the reduction, j, as the number of primary dimensions in the problem. Then k  n  j . [Reduce j by
one if necessary.]
4. Choose j repeating variables.
5. Construct the k s, and manipulate as necessary.
6. Write the final functional relationship 1  function   2 ,  3 ,... k  and check your algebra.
 Dimensional analysis is often extremely useful in setting up and designing experiments.
Errors and Calibration
 Systematic errors (bias errors) – consistent, repeatable errors.
 Random errors (precision errors) – scatter in data, a lack of repeatability, unrepeatable, inconsistent errors.
 Accuracy – accuracy error is the measured value minus the true value.
 Precision – precision error is the reading minus the average of readings.
 Other errors – zero, linearity, sensitivity, resolution, hysteresis, instrument reapeatability, drift.
 Calibration – static (time not relevant) vs. dynamic (time is relevant) calibration.
1 n x  xtrue
 Mean bias error – defined as MBE   i
, and usually expressed as a percentage.
n i 1 xtrue
Basic Electronics Review



 1 
Unity conversion factors – 
 1 ,
 1 V/A 
 1F 
 1 H  A/s 
1 A s 

 1, 
  1 , and 
 1 .
 1 C/V 
 1V 
 1C 
dV
dI
Ohm’s law and basic equations – V  IR , I  C
and V  L
at any instant in time.
dt
dt
R2
Voltage divider – Vout  Vin
where Vout is measured across resistor R2.
R1  R2

Resistors in series – Rtotal  R1  R2  ...  RN and I  I1  I 2  ...  I N .
1
Resistors in parallel – Rtotal 
and I total  I1  I 2  ...  I N .
1
1
1

 ... 
R1 R2
RN

Capacitors in series – Ctotal 

Capacitors in parallel – Ctotal  C1  C2  ...  C N and I total  I1  I 2  ...  I N .


Inductors in series – Ltotal  L1  L2  ...  LN and I  I1  I 2  ...  I N .
1
Inductors in parallel – Ltotal 
and I total  I1  I 2  ...  I N .
1
1
1

 ... 
L1 L2
LN

Impedance – Z 

Vpeak
I peak
1
and I  I1  I 2  ...  I N .
1
1
1

 ... 
C1 C2
CN
: resistor Z  R , capacitor Z 
1
, inductor Z  i L .
iC
Basic Statistics
n


n
 di 2
 xi  x 

1 n
i 1
i 1

Definitions – sample mean x   xi , sample standard deviation S 
, sample
n i 1
n 1
n 1
variance (S2), sample median (half lower, half higher), sample mode (most probable value – one that occurs
most frequently), population (all values) vs. sample (a selected portion of the total population).
Excel – learn how to use Excel’s built-in statistics functions.
2
2

1 n  xi  xtrue 
Root mean square error – defined as RMSE 
 , and usually expressed as a percentage.

n i 1  xtrue 
Histograms and Probability Density Functions
 Histogram plots – frequency, bins or classes, bin width or class width, Sturgis and Rice rules.
 Normalized histograms – how to normalize the vertical scale and the horizontal scale.
 PDF – how to create a probability density function from a histogram.
 Expected value – same as population mean.
 Standard deviation – same as population standard deviation.
x
 Normalized PDF – transform from x to z using z 
and f  z    f  x  .





The Gaussian or normal PDF – how to predict probabilities for systems with purely random errors using the
2  
1
 z 
exp   2 d  . See table of A(z).
error function, A  z   erf 
where erf   


 0
2

 2
Other PDFs – lognormal, chi-squared distribution, student’s t, degrees of freedom, confidence level and
significance level.
S
How to estimate the population mean  from a sample using the student’s t PDF:   x  t / 2
with
n
confidence level 1 – ; use table of critical values for the student’s t PDF.
How to estimate the population standard deviation  from a sample using the 2 PDF:
S2
S2
with confidence level 1 – ; use table of critical values for the 2 PDF.
 n  1 2   2   n  1 2


 /2

1 / 2
The central limit theorem (CLT) – standard error of the mean = standard deviation of sample means =
S

S x  overall . As n gets large, we write the CLT in terms of population standard deviations:  x 
.
n
n
Correlation and Regression
i n

Linear correlation coefficient – rxy 
 x
i
i 1
in
 x
i 1
i
 x  yi  y 
 x
2
in
 y
i 1
i
 y
. By definition, rxy must always lie between
2
1 and 1, i.e., 1  rxy  1 . To determine if there is a trend in the data, use the table of critical values of rxy.

Regression analysis – least-squares fits (linear, polynomial, and multiple variables), standard error of
estimate. Equations are provided in the lecture notes, but we typically use Excel to do the analysis.
Standard error of estimate – also called standard error: S y , x 
i n
in
i n
i 1
i 1
i 1
 yi 2  b yi  a xi yi
n2
.
Outlier Points
 Outliers in a set of data points – use the modified Thompson tau technique: compare maximum absolute
value of the deviation with Thompson  times sample standard deviation ( S) to determine outliers:
If  >  S, reject the outlier.
in

 y
i
 Yi 
2
, df  n  ( m  1) .
df
Examine the standardized residual ei / Sy,x to determine outliers. There are two criteria for a data pair to be
called an outlier:
o
ei / S y , x  2
Outliers in a set of data pairs – for polynomial fit of order m, S y , x 
o
i 1
The standardized residual is inconsistent with its neighbors (as judged by a plot of ei / Sy,x vs. x).
Experimental Uncertainty Analysis
 Maximum uncertainty – assumes all errors have the same sign and add up (unlikely). Normally, we do not
use the maximum uncertainty. Typically, we use instead the RSS uncertainty, as described below:
 RSS uncertainty or expected uncertainty – assumes some errors are positive and some are negative (more
likely). Use uR ,RSS 

2
 R 
uR u R ,RSS
aN
a1
a2


 u xi
 or, if R  C  x1  x2 ...  xN ,

R
R
xi 
i 1 
i N
Combining elemental uncertainties using RSS uncertainty analysis – use u x 
iK
u
i 1
2
 u xi 
 ai
 .

xi 
i 1 
iN
2
i
,where K is the
number of elemental uncertainties, and ui are the individual elemental uncertainties.
Experimental Design
 Choosing a test matrix – one parameter at a time (full factorial) vs. Taguchi’s experimental design arrays
P
(fractional factorial) test matrices. N  L (L = levels, P = parameters) for a full factorial experiment.
 Level averages – E.g., average over all the runs where parameter a is at level 1.
 Optimum Taguchi design array – Meets two criteria:
o Each level of each parameter appears the same number of times in the array.
o Repetitions of parameter-level combinations are minimized as much as possible.
 Taguchi orthogonal arrays – know how to spot a poorly designed experimental design array.
 Response surface methodology (RSM) – goal is to efficiently hunt for the optimum values of parameters a, b,
c, etc. such that response y is maximized or minimized.
 a  amid value 
 Coded variables – transform the physical variables into coded variables, x1  2 
 ,
 a
range


a  amax
amid value  min
, arange  amax  amin (with similar equations for the other variables). This is necessary to
2
properly use the RSM technique.
 y y y 
,
,
 Direction of steepest ascent – vector 
 determined by regression analysis on coded variables.
 x1 x2 x3 
For known x1, x2  x1

y
x2
y
y
and x3  x1
x3
x1
Marching – Transform back to physical variables, a 
ascent until y is no longer improving.
y
.
x1
arange
2
x1 , then march in the direction of steepest
Hypothesis Testing
 Null hypothesis – a theory that is being considered or tested. Typically, the null hypothesis represents
“nothing is happening.”
 There are two parts to the null hypothesis:
o The critical value (extreme value), o.
o The “sides” or “tails” of the null hypothesis – select the least likely scenario ( = 0,  < 0, or  > 0).
 Alternative hypothesis – also called the research hypothesis – the complement of the null hypothesis.
x  0
 t-statistic – calculate for the extreme value (critical value) of the null hypothesis, t 
.
S/ n
 p-value – p-value is the probability of wrongly rejecting the null hypothesis, if it is in fact true. In Excel, use
TDIST(t, df, tails) (p-value is the area under the tail(s) of the t-PDF). See also the tables of p-values for a
given value of df. Note: The values in the table are for one-tail hypothesis tests; multiply by 2 for two tails.
 A small p-value means that the null hypothesis is unlikely to be true – we reject the null hypothesis if the pvalue is less than some level, typically 0.05 (5%) to standard engineering confidence level.
 Paired samples hypothesis testing when n is the same for the two samples – use  = xB – xA as the variable,
  0
set the null hypothesis as 0    0 , and set the value of the critical t-statistic to t 
.
S / n
 Two sample hypothesis testing when n is not the same for the two samples – set the null hypothesis as
2


 S A2 S B 2 





nA
nB 
x A  xB



 A   B , set the critical t-statistic to t 
, and use df  NINT 
2
2
2 2
2 2 
SA
S
 1  S A   1  SB  
 B
 nA  1  n A  nB  1  nB  
nA
nB
(Welch’s equation) as a kind of weighted value of the degrees of freedom.
Digital Data Acquisition
 Binary to decimal and vice-versa – how to convert: sum columns or successive division, respectively.
Vmax  Vmin  quantization error   1 resolution   1 Vmax  Vmin 
 A/D systems – resolution 
,
.
2N
2
2
2N
 Discrete sampling – t  1/ f s ,   2 f , watch out for clipping and aliasing.

Ailasing – If f s  2 f , then there is no aliasing (Nyquist criterion), if
If f s 
 f
2
f , then f a   a
 f
3
 folding
2
f  f s  2 f , then f a  f s  f .
3

f
 f folding where f folding  s , and need to use the folding diagram.
2


 f 
OR, use the general equation f perceived  f  f s  NINT   , where NINT is the “nearest integer”.
 fs 
Signal reconstruction – the Cardinal series: reconstruct a digital signal back into an analog signal.

Fourier series – f  t   c0   an sin  n0t    bn cos  n0 t  where fundamental frequency is f 0  1/ T and


n 1
0  2 f 0 . Coefficients are c0 





n 1
T
T
T
1
2
2
f  t  dt , an   f  t  sin  n0t  dt , and bn   f  t  cos  n0 t  dt .

T0
T0
T0
Harmonic amplitude plots – plot of amplitude of each coefficient a1 , a2 , etc. versus harmonic number n.
fs N
 f where
2 2
N = total number of discrete data points taken, T = total sampling time, t  T / N = time between data
points, f s  1/ t  N / T = sampling frequency, and f  1/ T = frequency resolution.
Leakage – appears when the discrete data acquisition does not stop at exactly the same phase as it started.
Windowing – reduces but does not totally eliminate leakage.
Anti-aliasing filter – use a low-pass filter to remove high frequencies to avoid aliasing.
Fourier transforms, DFTs and FFTs – discrete version of Fourier series; f max  f folding 
Filters
Vout
 f
cutoff
1
1


   arctan 
, G
,
2
Vin
2
2 RC
 f cutoff
1   f f cutoff 

.

Vout
cutoff
 f
1
1


, G
,   arctan  cutoff
2
Vin
2
2 RC
 f
1   f cutoff f 

.


1st-order low-pass filter – f cutoff 

1st-order high-pass filter – f cutoff 

Higher-order filters – for a Butterworth low-pass filter of order n, G 

Gain expressed in decibels – for any filter or amplifier, GdB  20log10 G  20log10
Vout
Vin

1
1   f f cutoff 
Vout
Vin
2n
.
.
Operational Amplifiers (op-amps)
 Open-loop gain – Vo  g V p  Vn  .

Closed-loop configuration (with feedback loop) – V p  Vn .

Example circuits – buffer: Vout  Vin , inverting amplifier: Vout  
R2
Vin  GVin , inverter: Vout  Vin ,
R1
inverting summer: Vout   V1  V2  , noninverting amplifier: Vout  1   R2 R1  Vin  GVin .



Other circuits – active filters, clipping circuits (high and low voltage clipping), etc.
g
Common-mode rejection ratio (CMRR) – CMRR  20 log10
, and GCM is the common-mode gain.
GCM
Gain-bandwidth product (GBP) – op-amp acts like low-pass filter at high frequencies:
o Noninverting op-amp amplifier: GBP = GBPnoninverting  f c Gtheory, noninverting  constant [GBP supplied by
manufacturer’s specs].
 R2
GBPnoninverting , GBPinverting  f c Gtheory, inverting  constant .
R1  R2
Loading – input and output loading can cause voltage change due to internal resistances. Ri = input resistance
or input impedance.
o Inverting op-amp amplifier: GBPinverting 

Stress, Strain, and Strain Gages
F
L
, a 
,  a  E a .
A
L

Axial stress and axial strain –  a 

Transverse strain & Poisson’s ratio –   Poisson's ratio  



t
 w t
D
 ; t 
, t 
for round rods.
a
w
t
D
1
1
Stress-strain relationship on a surface; principal axes –  x   x  y  ,  y   y  x  .
E
E
L  R
R/R
 S  a , where S is the strain gage factor, defined as S 
Wire resistance – R 
,
.
R
A
a
Wheatstone bridge – Vo  Vs
R3 R1  R4 R2
, bridge is balanced when R3 R1  R4 R2 . Approximate
 R2  R3  R1  R4 
relationship when resistances change:


R2,initial R3,initial   R1
Vo
 R3
 R2
 R4 





.
Vs  R2,initial  R3,initial 2  R1,initial R2,initial R3,initial R4,initial 
4 Vo 1
n
or Vo   a SVs , where n = 1 for a quarter bridge. n = 2 for a
n Vs S
4
half bridge, and n = 4 for a full bridge.
4 Vo  Vo,reference  1
n
Unbalanced bridge –  a 
or Vo  Vo,reference   a SVs .
n
Vs
S
4
Quarter, half, and full bridge –  a 
Dynamic System Response

General governing equation (ODE): an

First-order system: 
o
o

dny
d n 1 y
d2y
dy


...

 a1  a0 y  bx .
a
a
n 1
2
n
n 1
2
dt
dt
dt
dt
dy
 y  Kx , K  b / a0 ,   a1 / a0 .
dt
t
t
y  yf
y  yi
 e .
 nondimensionalized output = 1  e  ,   t  
For step-function input,
yi  y f
y f  yi


 t /
For ramp function input, ymeasured  KA t   1  e    , time lag   , error  ymeasured  yideal   KA .


Second-order system:
a
a1
1 d 2 y 2 dy

 y  Kx , K  b / a0 , n  0 ,  
. For step-function
a2
n 2 dt 2 n dt
2 a0 a2
input:
o
Underdamped,  < 1,
 1
y  yi
sin n t 1   2  
 1  e nt 
2
y f  yi
 1  

Undamped natural frequency and ringing frequency: f n 
n
1

2 2

 ,   sin 
1

1  2 .
a0
2
, and d  n 1   or
a2
fd  fn 1   2 .
o
Critically damped,  = 1,
o
Overdamped,  > 1,
y  yi
 1  e nt 1  n t  , and there is no phase shift.
y f  yi


y  yi

 1  e nt  cosh n t  2  1 
sinh n t  2  1  .
y f  yi


 2 1




Temperature Measurement
 Mechanical devices – liquid-in-glass thermometers, bimetallic strips, pressure thermometers.
 Thermojunctive devices – thermocouples, sensing junctions, reference junctions, type K, T, etc.
 Thermocouple laws – notation: V1 2  V1 R  V2  R .
o Law of intermediate metals: A third (intermediate) metal wire can be inserted in series with one of the
wires without changing the voltage reading (provided that the two new junctions are at the same
temperature).
o Law of intermediate temperatures: If identical thermocouples measure the temperature difference
between T1 and T2, and the temperature difference between T2 and T3, then the sum of the corresponding
voltages V1–2 + V2–3 must equal the voltage V1–3 generated by an identical thermocouple measuring the
temperature difference between T1 and T3, i.e., V13  V1 2  V2 3 for any 3 temperatures T1 , T2 , and T3 .
o Statement of the law of additive voltages: For a given set of 3 thermocouple wires, A, B, and C, all
measuring the same temperature difference T1  T2, the voltage measured by wires A and C must equal
the sum of the voltage measured by wires A and B and the voltage measured by wires B and C, i.e.,
V1 2,A&C  V1 2,A&B  V1 2,B&C .
 Thermopile – Several thermocouples in series.
 Thermoresistive devices – RTDs (metal, R increases with increasing T), thermistors (semiconductor, R
decreases with increasing T).
1/ 4




W
TH   assumed  Tind .
Radiative devices – e.g., infrared pyrometer: E   T ,   5.669  10
2 4 ,
mK
  actual 
Tables – how to use thermocouple tables and RTD tables, and interpolation. Reference temperature is 0oC for
all temperature tables.
4
8
Measurement of Mechanical Quantities
 Position – mechanical devices; interferometer; potentiometer; linear variable displacement transducer
a t
(LVDT); ultrasonic transducer: pulse-echo x 
, through transmission x  at ; capacitance sensor:
2
C2
C  K  0 A / d ,  0  8.854  1012
; laser displacement meter.
N  m2
rotation   radian  rotation  60 s 

 Angular velocity and rpm –  N rpm
  


 ; contacting tachometer,
min  
s  2 radian  min 

P
noncontacting tachometer: N rpm 
. For stroboscopic tachometers, watch out for aliasing.
nteeth

2
N rpmT ; dynamometers: prony brake, cradled DC motor, eddy current.
Torque and power – Wshaft  T 
60
Pressure Measurement
 Types of pressure – absolute, gage: Pgage  Pabs  Patm , vacuum: Pvac  Patm  Pabs , Pvac   Pgage .

Hydrostatics – Pbelow  Pabove   g z ; applications include:
o
Liquid manometers (U-tube manometers): typically, P    m    gh if measurement locations are at the
same elevation on either side of the manometer.
McCleod gage – a special type of liquid manometer used to measure very low pressures.
F
A
o Hydraulic jacks – the ideal mechanical advantage is 2  2 .
F1 A1
Mechanical pressure gages – work by mechanical means, without any electricity or liquid columns:
o Bourdon tube – flattened hollow tube that is coiled, bent, or twisted.
o Deadweight tester – weights applied to a piston of known area, thus P  F / Ae .
Electronic pressure transducers – diaphragm type, with displacement measured by:
o Strain gage – a strain gage is mounted on the diaphragm itself, sensing strain in the diaphragm.
o Capacitance – the diaphragm is mounted close to a fixed parallel plate, capacitance is measured.
o LVDT – diaphragm is attached to the core of a linear variable displacement transducer.
o Optical – various optical techniques used to measure the degree of diaphragm deformation.
Piezoelectric pressure transducers – compression causes a voltage that can be calibrated with pressure.
o



Velocity Measurement
 Linear velocity transducer (LVT) – Similar in principle to LVDT, but measure velocity, not displacement.
2V cos 
 Doppler radar velocimeter – Doppler frequency shift f D 
.

 Displacement sensors – V (t )  dx (t ) / dt , noise is amplified by differentiation.
t

Acceleration sensors – V (t )  V0   a (t ) dt , noise is attenuated by integration.


Velocity of fluids – Lagrangian methods follow fluid particle moving with the flow: V (t )  dx (t ) / dt ;
Eulerian methods measure velocity field with a probe sitting in the flow.
f
Laser Doppler velocimeter (LDV) – particles in fringe pattern of focal volume, V  fs 
.
2sin( / 2)
Particle image velocimetry (PIV) – velocity vectors calculated from moving particles in the plane illuminated
by a laser sheet, V  s / t .
Hot-wire and hot-film anemometer – wire heated to constant temperature; King’s law E 2  a  bV n .

Pitot-static probe – Pitot formula V 



t0
2  P1  P2 
.

Other – rotating velocimeter (cup anemometer, turbine or vane anemometer), electromagnetic velocimeter.
Volume Flow Rate Measurement
  V . Notation: V = volume, V = velocity, V = volume flow rate.
 Mass flow rate vs. volume flow rate – m
 End-line measurement – V  V / t vs. in-line measurement – various types:

2  P1  P2 
Obstruction – orifice, flow nozzle, and Venturi:   d / D , V  Ao Cd
, where d = orifice
 1   4 
diameter, D = pipe diameter, Ao = orifice area, V1= average speed through the pipe, V2= average speed
through the orifice (or throat), P1 – P2 = measured pressure drop from upstream to downstream, Cd =
V D V D
discharge coefficient, and Reynolds number is based on the pipe not the orifice, i.e., Re  1  1 :







o
For a standard orifice flowmeter, Cd  0.5959  0.0312  2.1  0.184 8.0  91.71 2.5 / Re0.75
o
For a standard flow nozzle flowmeter, Cd  0.9975  6.53 0.50 / Re0.50
o
For a standard Venturi flowmeter, Cd  0.98
Laminar flow element – pressure drop calibrated against volume flow rate, typically linear calibration
because flow through the small tubes is laminar even though flow through the pipe may be turbulent.
Positive displacement – measure how many “trapped” or “captured” volumes pass through per unit time.
Turbine and paddlewheel – spin a turbine or paddlewheel connected to a shaft, measure rpm and calibrate.
Roatmeters or floatmeters or variable-area flowmeters – float hovers when forces balance.
Miscellaneous flowmeters – ultrasonic (transit time and Doppler shift), electromagnetic, vortex, others.