6.1.5.1: Momentum for For Constant Pressure and Frictionless Flow
- Page ID
- 716
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Another important sub category of simplification deals with flow under approximation of the frictionless flow and uniform pressure. This kind of situations arise when friction (forces) is small compared to kinetic momentum change. Additionally, in these situations, flow is exposed to the atmosphere and thus (almost) uniform pressure surrounding the control volume. In this situation, the mass flow rate in and out are equal. Thus, equation (15) is further reduced to
\[
\label{mom:eq:mom:eq:egRTTststS}
\pmb{F} = \int_{out} \rho \pmb{U} \overbrace{\left( \pmb{U} \cdot \hat{n} \right)}^{U_{rn}} dA
- \int_{in} \rho \pmb{U} \overbrace{\left( \pmb{U} \cdot \hat{n} \right)}^{U_{rn}} dA
\]
\[
\label{mom:eq:egRTTststSU}
\pmb{F} = \dot{m} \overline{\pmb{U}_o} - \dot{m} \overline{\pmb{U}_i}
\]
The average velocity is related to the velocity profile by the following integral
\[
\label{mom:eq:UmaxUaveI}
{\overline{U}}^2 = \dfrac{1}{A}\int_A \left[U(r)\right]^2\,dA
\]
Equation (18) is applicable to any velocity profile and any geometrical shape.
Example 6.1
Calculate the average velocity for the given parabolic velocity profile for a circular pipe.
Solution 6.1
The velocity profile is
\[
\label{mom:eq:Uprofilec}
U\left(\dfrac{r}{R}\right) = U_{max} \left[ 1 - \left( \dfrac{r}{R} \right)^2 \right]
\]
\[
\label{mom:eq:UmaxUaveIC}
{\overline{U}}^2 = \dfrac{1}{2\,\pi\,R^2}\int_0^R \left[U(r)\right]^2\,2\,\pi\,r\, dr
\]
results in
\[
\label{mom:eq:Uave}
{\overline{U}}^2 = \left(U_{max}\right)^2\, \int_0^1\left(1-\bar{r}^2\right)^2 \bar{r} d\bar{r} = \dfrac{1}{6}
\left(U_{max}\right)^2
\]
Thus,
\begin{align*}
\overline{U} = \dfrac{U_{max}}{\sqrt{6}}
\end{align*}
Figure 6.2 Schematics of maximum angle for impinged by a jet
Example 6.2
A jet is impinging on a stationary surface by changing only the jet direction (see Figure 6.2). Neglect the friction, calculate the force and the angle which the support has to apply to keep the system in equilibrium. What is the angle for which maximum force will be created?
Solution 6.2
Equation (11) can be reduced, because it is a steady state, to
\[
\label{areaJet:gov}
\pmb{F} = \int_{out} \rho \pmb{U} \overbrace{\left( \pmb{U} \cdot \hat{n} \right)}^{U_{rn}} dA
- \int_{in} \rho \pmb{U} \overbrace{\left( \pmb{U} \cdot \hat{n} \right)}^{U_{rn}} dA
= \dot{m} \pmb{U_o} - \dot{m} \pmb{U_i}
\]
\begin{align*}
F_x = \dot{m}\,\left(\cos\theta\, {U_o} - {U_i} \right)
\end{align*}
or since $U_i = U_o$
\begin{align*}
F_x = \dot{m}\,{U_i}\,\left(\cos\theta - 1 \right)
\end{align*}
It can be observed that the maximum force, \(F_x\) occurs when \(\cos\theta=\pi\). It can be proved by setting \(dF_x/d\theta=0\) which yields \(\theta=0\) a minimum and the previous solution. Hence,
\begin{align*}
\left. F_x\right|_{max} = - 2\,\dot{m}\,{U_i}
\end{align*}
and the force in the \(y\) direction is
\begin{align*}
F_y = \dot{m} \, {U_i} \,\sin\theta
\end{align*}
the combined forces are
\begin{align*}
F_{total} = \sqrt{{F_x}^2 + {F_y}^2} = \dot{m}\,{U_i}
\sqrt{\left( \cos\theta - 1\right)^2 + \sin^2\theta}
\end{align*}
Which results in
\begin{align*}
F_{total} = \dot{m}\,{U_i} \sin \left( \theta/2\right)
\end{align*}
with the force angle of
\begin{align*}
\tan \phi = \pi - \dfrac{F_y}{F_x} = \dfrac{\pi}{2} - \dfrac{\theta}{2}
\end{align*}
For angle between \( 0 <\theta < \pi\) the maximum occur at \(\theta = \pi\) and the minimum at \(\theta \sim 0\). For small angle analysis is important in the calculations of flow around thin wings.
Example 6.3
Liquid flows through a symmetrical nozzle as shown in the Figure 6.3 with a mass
Fig. 6.3 Nozzle schematic for the discussion on the forces and for example .
flow rate of 0.01 [\(gk/sec\)]. The entrance pressure is 3[\(Bar\)] and the entrance velocity is 5 [\(m/sec\)]. The exit velocity is uniform but unknown. The exit pressure is 1[Bar]. The entrance area is 0.0005[\(m^2\)] and the exit area is 0.0001[\(cm^2\)]. What is the exit velocity? What is the force acting the nozzle? Assume that the density is constant \(\rho = 1000 [kg/m^3]\) and the volume in the nozzle is 0.0015 [\(m^3\)].
Solution 6.3
The chosen control volume is shown in Figure 6.3. First, the velocity has to be found. This situation is a steady state for constant density. Then
\begin{align*}
A_{1}\,U_{1} = A_{2}\,U_{2}
\end{align*}
and after rearrangement, the exit velocity is
\begin{align*}
U_{2} = \dfrac{A_{1} }{A_{2} } \, U_{1} = \dfrac{0.0005}{0.0001} \times 5 = 25 [m/sec]
\end{align*}
Equation (12) is applicable but should be transformed into the z direction which is
\begin{multline*}
\sum F_z + \int_{c.v.} \pmb{g}\cdot \hat{k} \,\rho\, dV + \int_{c.v.} \pmb{P}\cos\theta_z\, dA +
\int _{c.v.} \boldsymbol{\tau}_z \,dA = \\
\overbrace{\dfrac{t}{dt} \int_{c.v.} \rho\,\pmb{U}_z\,dV }^{=0} +
\int_{c.v.} \rho\,\pmb{U}_z\cdot\pmb{U}_{rn} dA
\end{multline*}
The control volume does not cross any solid body (or surface) there is no external forces. Hence,
\begin{multline*}
\overbrace{\sum F_z}^{=0} + \int_{c.v.} \pmb{g}\cdot \hat{k} \,\rho\, dV +
\overbrace{\int_{c.v.} \pmb{P}\cos\theta_z\, dA}^{\text{ liquid surface}} + \\
\overbrace{\overbrace{\int_{c.v.} \pmb{P}\cos\theta_z\, dA}^{\text{solid surface}} + \int _{c.v.} \boldsymbol{\tau}_z \,dA}^{\text{forces on the nozzle $F_{nozzle}$}} =
\int_{c.v.} \rho\,\pmb{U}_z\cdot\pmb{U}_{rn} dA
\end{multline*}
All the forces that act on the nozzle are combined as
\[
\label{mom:eq:egRTTcompbined}
\sum F_{nozzle} + \int_{c.v.} \pmb{g}\cdot \hat{k} \,\rho\, dV +
\int_{c.v.} \pmb{P}\cos\theta_z\, dA
= \int_{c.v.} \rho\,\pmb{U}_z\cdot\pmb{U}_{rn} dA
\]
\begin{align*}
\pmb{F}_b = - \int_{c.v.} \pmb{g}\cdot\hat{n}\,\rho \, dV = - g\,\rho V_{nozzle}
\end{align*}
Notice that in the results the gravity is not bold since only the magnitude is used. The part of the pressure which act on the nozzle in the \(z\) direction is
\begin{align*}
- \int_{c.v.} P dA = \int_{1} P dA - \int_{2} P dA
= \left.PA\right|_{1} - \left.PA\right|_{2}
\end{align*}
The last term in equation (23) is
\begin{align*}
\int_{c.v.} \rho\,\pmb{U}_z\cdot\pmb{U}_{rn} dA =
\int_{A_{2}} U_{2} \left( U_{2} \right) dA -
\int_{A_{1}} U_{1} \left( U_{1} \right) dA
\end{align*}
which results in
\begin{align*}
\int_{c.v.} \rho\,\pmb{U}_z\cdot\pmb{U}_{rn} dA =
\rho\, \left( {U_{2}}^2 A_{2} - {U_{1}}^2 A_{1} \right)
\end{align*}
Combining all transform equation (23) into
\[
\nonumber
F_z = - g\,\rho V_{nozzle} + \left.PA\right|_{2} - \left.PA\right|_{1} +
\rho\, \left( {U_{2}}^2 A_{2} - {U_{1}}^2 A_{1} \right)
\]
\begin{align*}
F_z = 9.8 \times 1000 \times
\end{align*}
Contributors and Attributions
Dr. Genick Bar-Meir. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or later or Potto license.