Homework 3 Solution - Purdue College of Engineering

3rd Homework Solution Manual
Problem 1 Part A: (4 points for derivation, 0.5 points for the answer and 0.5 points for the units)
The free body diagram for this problem is presented above. In order for the door to swing open the
moment created by the water pressure has to offset the weight moment. Figure above shows the
positive moment (M), x, y, orientations water height (H) and pressure vectors (P) but you can define
these as you want as long as you stay consistent with your definition.
Because the door in this example is rectangular the weight is applied at the centroid and the moment
can be calculated as:
๐‘€๐‘‘ = ๐‘Š๐‘ก โˆ—
๐ฟ
2
The moment generated by the water pressure is calculated bellow
๐ฟ
๐‘€๐‘ค = ๐น โˆ— ๐‘ฅ = ๐‘ƒ โˆ— ๐ด โˆ— ๐‘ฅ = โˆซ ๐‘ฅ โˆ— ๐‘ ๐‘‘๐ด = โˆซ ๐‘ ๐‘ค ๐‘ฅ ๐‘‘๐‘ฅ = ๐‘๐‘ค โˆซ ๐‘ฅ ๐‘‘๐‘ฅ
0
Setting the weight moment equal to water pressure moment would yield the following
๐‘=
๐‘Š๐‘ก
๐ฟ๐‘ค
Since the pressure in a static incompressible fluid equals
๐‘ƒ = ๐œŒ๐‘”๐ป
the height of the water needed to open the door can be computed as:
๐ป=
๐‘Š๐‘ก
=
๐ฟ ๐‘ค ๐œŒ ๐‘” (3๐‘“๐‘ก)(3๐‘“๐‘ก)
400๐‘™๐‘๐‘š
= 0.71๐‘“๐‘ก
๐‘ ๐‘™๐‘ข๐‘”๐‘ 
๐‘™๐‘๐‘š
(1.94
)
(32.174
)
๐‘ ๐‘™๐‘ข๐‘”
๐‘“๐‘ก 3
Part B: (4 points for derivation, 0.5 points for the answer and 0.5 points for the units)
This is similar to part A but now but because the hinge is replaced by a rail, the door can only move
up/down and the door will open when the force created by the water will offset the doorโ€™s weight.
๐ฟ
๐‘Š๐‘ก = ๐‘ƒ โˆ— ๐ด = โˆซ ๐‘ ๐‘ค ๐‘‘๐‘ฅ
0
Again, since the fluid is incompressible we can rewrite equation above as
๐ป=
๐‘Š๐‘ก
400๐‘™๐‘๐‘š
=
= 0.71๐‘“๐‘ก
๐ฟ ๐‘ค ๐œŒ ๐‘” (3๐‘“๐‘ก)(3๐‘“๐‘ก) (1.94 ๐‘ ๐‘™๐‘ข๐‘”๐‘ ) (32.174 ๐‘™๐‘๐‘š )
๐‘ ๐‘™๐‘ข๐‘”
๐‘“๐‘ก 3
Note: answers for parts A and B are identical only because the door was rectangular and horizontal.
Problem 2 Part A (2 for a, d, e, f and 1 point for b, c)
u ๏€ฝ F (r ) x means that the velocity is a product of the position vector and a scalar function F. To
multiply a vector by a scalar, multiply each component of a vector by that scalar.
The position vector is:
๐‘ฅโƒ— = ๐‘ฅ ๐‘– + ๐‘ฆ ๐‘— + ๐‘ง๐‘˜
And F * x becomes
๐น โˆ— ๐‘ฅ = ๐น๐‘ฅ ๐‘– + ๐น๐‘ฆ ๐‘— + ๐น๐‘ง ๐‘˜
Since:
๐‘Ÿ = โˆš๐‘ฅ 2 + ๐‘ฆ 2 + ๐‘ง 2
i,j and k components of velocity become:
๐‘ข = ๐น (โˆš๐‘ฅ 2 + ๐‘ฆ 2 + ๐‘ง 2 ) ๐‘ฅ
๐‘ฃ = ๐น (โˆš๐‘ฅ 2 + ๐‘ฆ 2 + ๐‘ง 2 ) ๐‘ฆ
๐‘ค = ๐น (โˆš๐‘ฅ 2 + ๐‘ฆ 2 + ๐‘ง 2 ) ๐‘ง
Part B
This vector field below was plotted using Maple. As can be seen, the fluid is moving away from the origin
at an increasing rate.
Part C
This vector field below was plotted using Maple. As can be seen, the fluid is moving toward the origin at
a decreasing rate.
Part D
The divergence of U is defined as:
โˆ‡โˆ—๐‘‰ =
๐‘‘๐‘ข ๐‘‘๐‘ฃ ๐‘‘๐‘ค
+
+
๐‘‘๐‘ฅ ๐‘‘๐‘ฆ ๐‘‘๐‘ง
Since F is constant
โˆ‡โˆ—๐‘‰ =
๐‘‘(๐น๐‘ฅ) ๐‘‘(๐น๐‘ฆ) ๐‘‘(๐น๐‘ง)
+
+
= ๐น + ๐น + ๐น = 3๐น
๐‘‘๐‘ฅ
๐‘‘๐‘ฆ
๐‘‘๐‘ง
Positive values value of velocity divergence means that the fluid is expanding. See the 2nd homework
help session notes for more information about the divergence.
Part E
Since F is defined as a function of r, you should use the spherical coordinates to solve this part and
define X as the following:
๐‘ฅโƒ— = ๐‘Ÿ ๐‘Ÿโƒ—
After you multiply x by F you should get
๐‘ข
โƒ—โƒ— = ๐ด ๐‘Ÿ โˆ’๐‘› โˆ— ๐‘Ÿ
The divergence of u in spherical coordinates is
Since we only have the radial component of u, the divergence of u will become
โˆ‡โˆ—๐‘ˆ =
(๐ด ๐‘Ÿ 3โˆ’๐‘› )
1
๐‘Ÿ 2 (๐ด ๐‘Ÿ1โˆ’๐‘› )
1
1
(
๐‘‘
=
(
๐‘‘
(
))
(
)) = 2 ( (3 โˆ’ ๐‘›)๐ด๐‘Ÿ 2โˆ’๐‘› ) = (3 โˆ’ ๐‘›)๐ด๐‘Ÿ โˆ’๐‘›
2
2
๐‘Ÿ
๐‘‘๐‘Ÿ
๐‘Ÿ
๐‘‘๐‘Ÿ
๐‘Ÿ
The fluid will be incompressible if n = 3 or A = 0
The fluid will be expanding if A > 0 and n <3 or A < 0 and n > 3
The fluid will be contracting if A > 0 and n >3 or A < 0 and n <3
Part F
The volumetric flow rate is
๐‘„ = โˆฎ ๐‘ข ๐‘‘๐ด
And the Gaussesโ€™ theorem states
โˆฎ ๐‘ข ๐‘‘๐ด = โˆญ โˆ‡๐‘ข ๐‘‘๐‘‰
Since the volume of a sphere is
4
๐‘‰ = ๐œ‹ ๐‘Ÿ3
3
And a derivative of that is
๐‘‘๐‘‰ = 4๐œ‹๐‘Ÿ 2 ๐‘‘๐‘Ÿ
We can calculate the flow rate as
๐‘…
๐‘…
๐‘…
๐‘„ = โˆซ โˆ‡๐‘ข 4๐œ‹๐‘Ÿ 2 ๐‘‘๐‘Ÿ = โˆซ โˆ‡๐‘ข 4๐œ‹๐‘Ÿ 2 ๐‘‘๐‘Ÿ = 4๐œ‹๐ด โˆซ (3 โˆ’ ๐‘›)๐‘Ÿ 2โˆ’๐‘› ๐‘‘๐‘Ÿ = 4๐œ‹๐ด๐‘… 3โˆ’๐‘›
0
0
0
Problem 3(Part A -> 3 pts for derivations, 1 pts for answer, Part B -> 5pts for derivation, 1 pt 6 points for
answer)
Schematic
3
1
2
Part A
To solve this problem simply apply the definition of Volume Flow Rate
10๐‘๐‘š3
๐‘„
12.73๐‘๐‘š
๐‘ 
๐‘„ =๐‘‰โˆ—๐ด โ†’ =๐‘‰ =
=
2
๐ด
0.25๐œ‹ ๐‘๐‘š
๐‘ 
Part B
Apply the conservation of volume (since density is constant) to the control volume shown in a schematic
above. This control volume has 3 surfaces. The first is a hemisphere, with a radius of 10cm. The second
is the circular hole and the third is a circle with a hole in the middle. Remember that if the fluid is leaked
through one of the surfaces of this control volume then some other fluid has to replace it through
another control volume. Since fluid cannot flow through the third control surface, there is a wall there;
the flow rate through the 1st control surface has to match the flow rate through the 2nd control surface.
๐‘„1 = ๐‘„2 = ๐‘‰2 โˆ— ๐ด2
Area 2 is the surface area of the hemisphere (2ฯ€r2) and V2 is the mean velocity flowing perpendicular to
the hemisphere. Therefore
10๐‘๐‘š3
๐‘„1
๐‘๐‘š
๐‘ 
๐‘‰2 =
=
= 0.016
2
๐ด2
2๐œ‹(10๐‘๐‘š)
๐‘ 
Problem 4(2 points for the plot, 4 points for correct integral setup, 1 point for flow rate answer, 2 points
for part c derivation, 1 point for correct answer on part c)
Schematic
Part A
Velocity Profile
12
Velocity (cm/s)
10
8
6
4
2
0
-2
-1
0
1
Distance from the center (cm)
Part B
๐‘Ÿ 2
๐‘„ = โˆซ( 10 ( 1 โˆ’ ( ) ) (2ฯ€r dr) = 20๐œ‹ ๐‘๐‘š3 /๐‘ 
๐‘…
Part C
๐‘‰=
๐‘„ 20๐œ‹ ๐‘๐‘š3 /๐‘ 
=
= 5 ๐‘๐‘š/๐‘ 
๐ด ๐œ‹(2๐‘๐‘š)^2
Note: Centerline velocity is twice the average velocity in a round pipes.
2
Problem 5(2 points for the plot, 4 points for correct integral setup, 1 point for flow rate answer, 2 points
for part c derivation, 1 point for correct answer on part c)
Part A
Velocity Profile
Velocity (cm/s)
12
10
8
6
4
2
0
0
0.5
1
1.5
Distance from the lower surface (cm)
Part B
2
๐‘„ = โˆซ 10๐‘ฆ (2 โˆ’ ๐‘ฆ)๐‘‘๐‘ฆ =
0
40
๐‘๐‘š2
= 13.33
3
๐‘ 
Part C
๐‘‰=
๐‘„ 13.3 20
๐‘๐‘š
=
=
= 6.67
๐ด
2
3
๐‘ 
Note: Q in this problem is volume flow rate per width of the channel.
2