Question 5 of 6: Buckingham Pi analysis of a sudden contraction
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Exams, May 2016 —
07-Mec-B6 Advanced Fluid Mechanics. Three hours, OPEN BOOK, any
approved Sharp or Casio calculator permitted. Six questions are printed; any five
of them constitute a complete paper, each carries an equal 20 marks, and the item weights
are shown in the left margin. No aid sheet is bound into this paper — the open-book rule is the candidate’s table source, so every compressible-flow ratio below is quoted in closed form rather than read off a chart. All six questions are solved here.
Reference texts. Solutions follow the conventions of the texts the
EGBC syllabus recommends for this subject:
F. M. White, Fluid Mechanics, 8th ed. — potential-flow building blocks
(§4.4, §8.2–8.3), laminar boundary layers and the momentum integral
(§7.2–7.4), dimensional analysis (§5.2–5.4), compressible duct
flow (§9.5–9.7).
F. M. White, Viscous Fluid Flow, 3rd ed. — exact solutions of the
Navier–Stokes equations in cylindrical coordinates (§3.2–3.4) and
integral boundary-layer methods (§4.5).
J. D. Anderson, Modern Compressible Flow, 3rd ed. — quasi-one-dimensional
nozzle flow, normal shocks and Fanno flow (Ch. 3, Ch. 5).
P. K. Kundu, I. M. Cohen & D. R. Dowling, Fluid Mechanics, 6th ed. —
complex potential, images and the Blasius/Lagally force theorems (Ch. 6); boundary
layers (Ch. 9).
B. R. Munson et al., Fundamentals of Fluid Mechanics, 8th ed. —
Buckingham Pi method and the repeating-variable procedure (Ch. 7).
Question 5: Buckingham Pi analysis of a sudden contraction (20 marks)
Given. The functional statement
Δp = f(D1, D2, U,
ρ, μ), where U is the velocity in the larger pipe, with
D1, U and μ nominated as the repeating variables.
Find. A complete, independent set of dimensionless groups, and a reasoned
explanation of why the velocity in the smaller pipe must not be added to the variable list.
The sudden contraction. The pressure drop is measured between a tap upstream in the large pipe and one downstream in the small pipe, beyond the vena contracta where the jet has refilled the bore.
Approach. Count variables and independent dimensions to fix the
number of groups, confirm the nominated repeating set is dimensionally independent, then
solve one small linear system per non-repeating variable.
Count. The variable list is
Δp, D1, D2, U, ρ, μ, so
n = 6. In the MLT system their dimensions are
$$[\Delta p]=ML^{-1}T^{-2},\quad [D_{1}]=[D_{2}]=L,\quad [U]=LT^{-1},
\quad [\rho]=ML^{-3},\quad [\mu]=ML^{-1}T^{-1},$$
and the rank of the dimensional matrix is k = 3. Buckingham’s theorem
therefore predicts
$$n-k=6-3=3\ \text{dimensionless groups}.$$
Confirm the nominated repeating variables are admissible. The set
{D1, U, μ} must be dimensionally independent and must not by
itself form a dimensionless group. Mass appears only in μ, time only in U and μ,
and length in all three, and no product
D1aUbμc
is dimensionless unless a = b = c = 0. The set is valid, and note
in passing that it deliberately excludes the dependent variable Δp, as the
method requires.
First group: the dependent variable. Seek
Π1 = Δp D1aUbμc.
Balancing exponents,
$$M:\ 1+c=0,\qquad T:\ -2-b-c=0,\qquad L:\ -1+a+b-c=0,$$
which gives c = −1, b = −1, a = 1, so
$$\boxed{\;\Pi_{1}=\frac{\Delta p\,D_{1}}{\mu U}\;}$$
a viscous-scaled pressure drop — the natural form when viscosity is a repeating
variable.
Second group: the geometry. With
Π2 = D2D1aUbμc,
the mass balance gives c = 0, the time balance b = 0, and the length balance
a = −1:
$$\boxed{\;\Pi_{2}=\frac{D_{2}}{D_{1}}\;}$$
the contraction ratio, which is the only geometric parameter the problem contains.
Third group: the remaining property. With
Π3 = ρD1aUbμc,
$$M:\ 1+c=0,\qquad T:\ -b-c=0,\qquad L:\ -3+a+b-c=0,$$
giving c = −1, b = 1, a = 1, so
$$\boxed{\;\Pi_{3}=\frac{\rho U D_{1}}{\mu}=\mathrm{Re}_{D_{1}}\;}$$
the Reynolds number based on the upstream pipe.
State the result. The dimensionless statement of the problem is
$$\frac{\Delta p\,D_{1}}{\mu U}=\phi\!\left(\frac{D_{2}}{D_{1}},\ \mathrm{Re}_{D_{1}}\right),$$
a reduction from a function of five variables to a function of two. Any independent
recombination is equally valid, and the one used in practice divides out the viscosity:
$$\frac{\Pi_{1}}{\Pi_{3}}=\frac{\Delta p}{\rho U^{2}}=\mathrm{Eu}
=\phi'\!\left(\frac{D_{2}}{D_{1}},\ \mathrm{Re}_{D_{1}}\right),$$
the Euler number, which is preferred because at high Reynolds number the loss coefficient
becomes almost independent of Re and collapses onto a single curve in
D2/D1. A worked instance: with Δp = 4.2 kPa,
D1 = 0.20 m, D2 = 0.10 m, U = 3 m/s and water,
Π1 = 2.79 × 105, Π2 = 0.5,
Π3 = 5.98 × 105 and Eu = 0.468.
Part (b) — why the small-pipe velocity must be excluded. The
velocity in the smaller pipe is not an independent quantity. Conservation of mass for an
incompressible fluid fixes it completely from variables already on the list:
$$U_{2}=U\left(\frac{D_{1}}{D_{2}}\right)^{2},\qquad\text{i.e.}\qquad
\frac{U_{2}}{U}=\Pi_{2}^{-2}.$$
Buckingham’s theorem requires the variable list to be independent; a quantity
that is an algebraic function of others is redundant.
What actually goes wrong if you add it. Adding U2
raises n to 7 while the rank stays at k = 3, so the method dutifully
produces a fourth group,
U2D22/(UD12).
But that group is identically equal to unity for every physically realisable flow: it is not a
free parameter, it is a constraint the fluid has already satisfied. Carrying it forward is
misleading in two ways. It suggests the problem has one more degree of freedom than it does,
so an experimenter might waste a whole test matrix trying to vary it independently; and it
breaks the uniqueness of the correlation, because the same physical state can now be labelled
by infinitely many nominally different parameter sets. The rule generalises: include the
independent variables that define the problem, and let the conservation laws supply everything
they determine.
Part
Group
Form
Name
(a)
Π1
Δp D1 / (μU)
viscous-scaled pressure drop
(a)
Π2
D2 / D1
contraction ratio
(a)
Π3
ρU D1 / μ
Reynolds number
(a)
Result
Π1 = φ(Π2, Π3)
equivalently Eu = Δp/ρU2 = φ′(Π2, Re)
(b)
Reason
U2 = U(D1/D2)2
not independent — fixed by continuity; yields a group identically equal to 1