A fluid of constant density rho enters a duct of width W and
Solution
a)
Consider a strip of height dy1 at distance y from the centre. Cross-section of the strip = W*dy1
For bulk mean velocity at entry, we have
(h1*W)*V1_mean = Integral [(W*dy1)*V] ............from y1 = -h1/2 to y1 = h1/2
(h1*W)*V1_mean = 2*Integral [(W*dy1)*V1* (1 - (2y1/h1)2)] ............from y1 = 0 to y1 = h1/2
h1*V1_mean = 2V1* Integral [dy1* (1 - (2y1/h1)2)] ............from y1 = 0 to y1 = h1/2
h1*V1_mean = 2V1* [y1 - (4/3)y13 / h12]............from y1 = 0 to y1 = h1/2
h1*V1_mean = 2V1* [(h1/2) - (4/3)(h1/2)3 / h12]
V1_mean = (2/3)V1
Similarly, we can get V2_mean = (2/3)V2
By mass conservation, rho*(h1*W)*V1_mean = rho*(h2*W)*V2_mean
h1*(2/3)V1 = h2*(2/3)V2
V2 = V1*(h1/h2)
b)
By momentum conservation,
F = Momentum in - Momentum out
= P1*(h1*W) + [rho*(h1*W)*V1_mean]*V1_mean - P2*(h2*W) - [rho*(h2*W)*V2_mean]*V2_mean
= rho*W*[h1*V1_mean2 - h2*V2_mean2] + (P1*h1*W - P2*h2*W)
= rho*W*[h1*(2V1/3)2 - h2*(2V2/3)2] + (P1*h1*W - P2*h2*W)
= rho*W*[h1*(2V1/3)2 - h2*(2V1*(h1/h2)/3)2] + (P1*h1*W - P2*h2*W)
= (4/9)*rho*W*h1*V1*[V1 - V1*(h1/h2)] + (P1*h1*W - P2*h2*W)
= (4/9)*rho*W*h1*V12*[1 - (h1/h2)] + (P1*h1*W - P2*h2*W)
