Differential relations for a fluid particle
edit
We are interested in the distribution of field properties at each point in space. Therefore, we analyze an infinitesimal region of a flow by applying the RTT to an infinitesimal control volume, or , to a infinitesimal fluid system.
Conservation of mass
edit
The differential control volume dV and the mass flux through its surfaces
The conservation of mass according to RTT
∂
∂
t
∫
C
V
ρ
d
V
+
∫
C
S
ρ
U
→
⋅
n
→
d
A
=
0
{\displaystyle \displaystyle {\frac {\partial }{\partial t}}\int _{CV}\rho \ dV+\int _{CS}\rho \;{\vec {U}}\cdot {\vec {n}}\ dA=0}
or in tensor form
∂
∂
t
∫
C
V
ρ
d
V
+
∫
C
S
ρ
U
i
n
i
d
A
=
0
{\displaystyle \displaystyle {\frac {\partial }{\partial t}}\int _{CV}\rho \ dV+\int _{CS}\rho \;U_{i}n_{i}\ dA=0}
The differential volume is selected to be so small that density
(
ρ
)
{\displaystyle \displaystyle (\rho )}
can be accepted to be uniform within this volume. Thus the first integral in 1 is:
∂
∂
t
∫
C
V
ρ
d
V
≈
∂
ρ
∂
t
d
x
1
d
x
2
d
x
3
=
∂
ρ
∂
t
d
V
{\displaystyle \displaystyle {\frac {\partial }{\partial t}}\int _{CV}\rho \ dV\approx {\frac {\partial \rho }{\partial t}}dx_{1}dx_{2}dx_{3}={\frac {\partial \rho }{\partial t}}dV}
The flux term (second integral term) in the equation of conservation of mass can be analyzed in groups:
∫
C
S
ρ
U
i
n
i
d
A
=
∫
C
S
x
1
ρ
U
i
n
i
d
A
+
∫
C
S
x
2
ρ
U
i
n
i
d
A
+
∫
C
S
x
3
ρ
U
i
n
i
d
A
{\displaystyle \displaystyle \int _{CS}\rho U_{i}n_{i}\ dA=\int _{CS\,x_{1}}\rho U_{i}n_{i}\;dA+\int _{CS\,x_{2}}\rho U_{i}n_{i}\;dA+\int _{CS\,x_{3}}\rho U_{i}n_{i}\;dA}
Let's look to the surfaces perpendicular to
x
1
−
a
x
i
s
{\displaystyle \displaystyle x_{1}-axis}
∫
C
S
x
1
ρ
U
i
n
i
d
A
=
−
ρ
U
1
d
x
2
d
x
3
+
[
ρ
U
1
+
∂
(
ρ
U
1
)
∂
x
1
d
x
1
]
d
x
2
d
x
3
{\displaystyle \displaystyle \int _{CS\,x_{1}}\rho U_{i}n_{i}\ dA=-\rho U_{1}dx_{2}dx_{3}+\left[\rho U_{1}+{\frac {\partial \left(\rho U_{1}\right)}{\partial x_{1}}}dx_{1}\right]dx_{2}dx_{3}}
=
∂
(
ρ
U
1
)
∂
x
1
d
x
1
d
x
2
d
x
3
=
∂
(
ρ
U
1
)
∂
x
1
d
V
{\displaystyle \displaystyle ={\frac {\partial \left(\rho U_{1}\right)}{\partial x_{1}}}dx_{1}dx_{2}dx_{3}={\frac {\partial \left(\rho U_{1}\right)}{\partial x_{1}}}dV}
Similarly, the flux integrals through surfaces perpendicular to
x
2
−
a
x
i
s
{\displaystyle \displaystyle x_{2}-axis}
and
x
3
−
axis
{\displaystyle \displaystyle x_{3}-{\textrm {axis}}}
are
∫
C
S
x
2
ρ
U
i
n
i
d
A
=
∂
(
ρ
U
2
)
∂
x
2
d
V
,
{\displaystyle \displaystyle \int _{CS\,x_{2}}\rho U_{i}n_{i}\ dA={\frac {\partial \left(\rho U_{2}\right)}{\partial x_{2}}}dV,}
∫
C
S
x
3
ρ
U
i
n
i
d
A
=
∂
(
ρ
U
3
)
∂
x
3
d
V
.
{\displaystyle \displaystyle \int _{CS\,x_{3}}\rho U_{i}n_{i}\ dA={\frac {\partial \left(\rho U_{3}\right)}{\partial x_{3}}}dV.}
Hence the flux integral reads;
∫
C
S
ρ
U
i
n
i
d
A
=
∂
(
ρ
U
i
)
∂
x
i
d
V
{\displaystyle \displaystyle \int _{CS}\rho U_{i}n_{i}\ dA={\frac {\partial \left(\rho U_{i}\right)}{\partial x_{i}}}dV}
The conservation of mass equation becomes:
∂
ρ
∂
t
d
V
+
∂
(
ρ
U
i
)
∂
x
i
d
V
=
0
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}dV+{\frac {\partial \left(\rho U_{i}\right)}{\partial x_{i}}}dV=0}
Droping the
d
V
{\displaystyle \displaystyle dV}
, we reach to the final form of the conservation of mass:
∂
ρ
∂
t
+
∂
(
ρ
U
i
)
∂
x
i
=
0
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}+{\frac {\partial \left(\rho U_{i}\right)}{\partial x_{i}}}=0}
This equation is also called continuity equation . It can be written in vector form as:
∂
ρ
∂
t
+
∇
⋅
(
ρ
U
→
)
=
0
where
∇
=
∂
∂
x
1
e
1
→
+
∂
∂
x
2
e
2
→
+
∂
∂
x
3
e
3
→
gradient operator
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}+\nabla \cdot \left(\rho {\vec {U}}\right)=0\ \ {\text{where}}\ \ \nabla ={\frac {\partial }{\partial x_{1}}}{\vec {e_{1}}}+{\frac {\partial }{\partial x_{2}}}{\vec {e_{2}}}+{\frac {\partial }{\partial x_{3}}}{\vec {e_{3}}}\ \ {\text{gradient operator}}}
For a steady flow, continuity equation becomes:
∂
(
ρ
U
i
)
∂
x
i
=
0
{\displaystyle \displaystyle {\frac {\partial \left(\rho U_{i}\right)}{\partial x_{i}}}=0}
.
For incompresible flow, i.e.
ρ
=
constant
{\displaystyle \displaystyle \rho ={\text{constant}}}
:
∂
U
i
∂
x
i
=
0
i.e.
∂
U
1
∂
x
1
+
∂
U
2
∂
x
2
+
∂
U
3
∂
x
3
=
0
{\displaystyle \displaystyle {\frac {\partial U_{i}}{\partial x_{i}}}=0\ \ {\text{i.e.}}\ \ {\frac {\partial U_{1}}{\partial x_{1}}}+{\frac {\partial U_{2}}{\partial x_{2}}}+{\frac {\partial U_{3}}{\partial x_{3}}}=0}
For a two dimensional, steady and incompressible flow in
x
1
x
2
{\displaystyle \displaystyle x_{1}x_{2}}
plane given by:
U
1
=
A
x
1
{\displaystyle \displaystyle \displaystyle U_{1}=Ax_{1}}
Find how many possible
U
2
{\displaystyle \displaystyle U_{2}}
can exist.
For incompressible steady flow:
ρ
∂
U
i
∂
x
i
=
0
{\displaystyle \displaystyle \rho {\frac {\partial U_{i}}{\partial x_{i}}}=0}
in two diemnsions
∂
U
1
∂
x
1
+
∂
U
2
∂
x
2
=
0
→
∂
U
2
∂
x
2
=
−
∂
U
1
∂
x
1
{\displaystyle \displaystyle {\frac {\partial U_{1}}{\partial x_{1}}}+{\frac {\partial U_{2}}{\partial x_{2}}}=0\ \ \rightarrow \ \ {\frac {\partial U_{2}}{\partial x_{2}}}=-{\frac {\partial U_{1}}{\partial x_{1}}}}
Thus,
∂
U
2
∂
x
2
=
−
A
{\displaystyle \displaystyle {\frac {\partial U_{2}}{\partial x_{2}}}=-A}
This is an expression for the rate of change of
U
2
{\displaystyle \displaystyle U_{2}}
velocity while keeping
x
1
{\displaystyle \displaystyle x_{1}}
constant. Therefore the integral of this equation reads
U
2
=
−
A
x
2
+
f
(
x
1
)
{\displaystyle \displaystyle \displaystyle U_{2}=-Ax_{2}+f(x_{1})}
Thus, any function
f
(
x
1
)
{\displaystyle \displaystyle f(x_{1})}
is allowable.
Compressible and unsteady flow inside a piston
Consider one-dimensional flow in the piston. The piston suddenly moves with the velocity
V
p
{\displaystyle \displaystyle V_{p}}
. Assume uniform
ρ
(
t
)
{\displaystyle \displaystyle \rho (t)}
in the piston and a linear change of velocity
U
1
{\displaystyle \displaystyle U_{1}}
such that
U
1
=
0
{\displaystyle \displaystyle U_{1}=0}
at the bottom (
x
1
=
0
{\displaystyle \displaystyle x_{1}=0}
) and
U
1
=
V
p
{\displaystyle \displaystyle U_{1}=V_{p}}
on the piston (
x
1
=
L
{\displaystyle \displaystyle x_{1}=L}
), i.e.
U
1
=
x
1
L
V
p
{\displaystyle \displaystyle U_{1}={\frac {x_{1}}{L}}V_{p}}
Obtain a function for the density as a function of time.
The conservation of mass equation is:
∂
ρ
∂
t
+
∂
(
ρ
U
i
)
∂
x
i
=
0
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}+{\frac {\partial \left(\rho U_{i}\right)}{\partial x_{i}}}=0}
For one-dimensional flow and uniform
ρ
{\displaystyle \displaystyle \rho }
, this equaiton simplifies to
∂
ρ
∂
t
+
ρ
∂
U
1
∂
x
1
=
0
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}+\rho {\frac {\partial U_{1}}{\partial x_{1}}}=0}
∂
ρ
∂
t
=
−
ρ
∂
U
1
∂
x
1
=
−
ρ
V
p
L
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}=-\rho {\frac {\partial U_{1}}{\partial x_{1}}}=-\rho {\frac {V_{p}}{L}}}
L
=
L
0
+
V
p
t
{\displaystyle \displaystyle \displaystyle L=L_{0}+V_{p}t}
∂
ρ
∂
t
=
d
ρ
d
t
=
−
ρ
V
p
L
0
+
V
p
t
{\displaystyle \displaystyle {\frac {\partial \rho }{\partial t}}={\frac {d\rho }{dt}}=-\rho {\frac {V_{p}}{L_{0}+V_{p}t}}}
∫
ρ
0
ρ
d
ρ
ρ
=
∫
0
t
−
V
p
L
0
+
V
p
t
d
t
{\displaystyle \displaystyle \int _{\rho _{0}}^{\rho }{\frac {d\rho }{\rho }}=\int _{0}^{t}-{\frac {V_{p}}{L_{0}+V_{p}t}}dt}
l
n
(
ρ
ρ
0
)
=
l
n
(
L
0
L
0
+
V
p
t
)
{\displaystyle \displaystyle ln\left({\frac {\rho }{\rho _{0}}}\right)=ln\left({\frac {L_{0}}{L_{0}+V_{p}t}}\right)}
ρ
(
t
)
=
ρ
0
(
L
0
L
0
+
V
p
t
)
{\displaystyle \displaystyle \rho (t)=\rho _{0}\left({\frac {L_{0}}{L_{0}+V_{p}t}}\right)}
The same problem can be solved by using the integral approach with a deforming control volume.
The differential equation of linear momentum
edit
The differential control volume dV and the flux of
ρ
U
i
{\displaystyle \displaystyle \rho \;U_{i}}
(momentum per unit volume in i -direction) through the surfaces perpendicular to
x
j
{\displaystyle \displaystyle x_{j}}
axis
The integral equation for the momentum conservation is
∑
F
i
=
∂
∂
t
∫
C
V
ρ
U
i
d
V
+
∫
C
S
ρ
U
i
U
j
n
j
d
A
{\displaystyle \displaystyle \sum F_{i}={\frac {\partial }{\partial t}}\int _{CV}\rho \;U_{i}\;dV+\int _{CS}\rho \;U_{i}\;U_{j}\;n_{j}dA}
For the first integral we assume
ρ
{\displaystyle \displaystyle \rho }
and
U
i
{\displaystyle \displaystyle U_{i}}
are uniform within dV, and dV is so small that:
∂
∂
t
∫
C
V
ρ
U
i
d
V
≈
∂
∂
t
(
ρ
U
i
)
d
x
1
d
x
2
d
x
3
{\displaystyle \displaystyle {\frac {\partial }{\partial t}}\int _{CV}\rho \;U_{i}\;dV\approx \ {\frac {\partial }{\partial t}}(\rho \;U_{i})dx_{1}dx_{2}dx_{3}}
Analyze the flux of the
ρ
U
i
{\displaystyle \displaystyle \rho U_{i}}
momentum terms through the faces perpendicular to the axis:
∫
C
S
ρ
U
i
U
j
n
j
d
A
=
∫
C
S
x
1
ρ
U
i
U
j
n
j
d
A
+
∫
C
S
x
2
ρ
U
i
U
j
n
j
d
A
+
∫
C
S
x
3
ρ
U
i
U
j
n
j
d
A
{\displaystyle \displaystyle \int _{CS}\rho \;U_{i}\;U_{j}\;n_{j}\;dA=\int _{CS\,x_{1}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA\ +\int _{CS\,x_{2}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA\ +\int _{CS\,x_{3}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA}
First consider the flux of
ρ
U
i
{\displaystyle \displaystyle \rho \;U_{i}}
(momentum per unit volume in i -direction) through the surfaces perpendicular to
x
1
{\displaystyle \displaystyle x_{1}}
axis:
∫
C
S
x
1
ρ
U
i
U
j
n
j
d
A
=
−
ρ
U
i
U
1
d
x
2
d
x
3
+
[
ρ
U
i
U
1
d
x
2
d
x
3
+
∂
∂
x
1
(
ρ
U
i
U
1
)
d
x
1
]
d
x
2
d
x
3
{\displaystyle \displaystyle \int _{CS\,x_{1}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA=-\rho \;U_{i}\;U_{1}\;dx_{2}dx_{3}\;+\left[\rho \ U_{i}\ U_{1}\ dx_{2}\ dx_{3}+{\frac {\partial }{\partial x_{1}}}(\rho \ U_{i}\ U_{1})\ dx_{1}\right]dx_{2}\ dx_{3}}
=
∂
∂
x
1
(
ρ
U
i
U
1
)
d
V
{\displaystyle \displaystyle ={\frac {\partial }{\partial x_{1}}}(\rho \;U_{i}\;U_{1})dV}
Similarly, the momentum flux through the surfaces in other directions read
∫
C
S
x
2
ρ
U
i
U
j
n
j
d
A
=
∂
∂
x
2
(
ρ
U
i
U
2
)
d
V
{\displaystyle \displaystyle \int _{CS\,x_{2}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA={\frac {\partial }{\partial x_{2}}}(\rho \;U_{i}\;U_{2})dV}
,
∫
C
S
x
3
ρ
U
i
U
j
n
j
d
A
=
∂
∂
x
3
(
ρ
U
i
U
3
)
d
V
{\displaystyle \displaystyle \int _{CS\,x_{3}}\rho \;U_{i}\;U_{j}\;n_{j}\;dA={\frac {\partial }{\partial x_{3}}}(\rho \;U_{i}\;U_{3})dV}
.
Rearranging the equation for
∑
F
i
{\displaystyle \displaystyle \sum F_{i}}
we obtain:
∑
F
i
=
[
∂
∂
t
(
ρ
U
i
)
+
∂
∂
x
j
(
ρ
U
i
U
j
)
]
d
V
{\displaystyle \displaystyle \sum F_{i}=\left[{\frac {\partial }{\partial t}}(\rho \ U_{i})+{\frac {\partial }{\partial x_{j}}}(\rho \ U_{i}\ U_{j})\right]dV}
We can simplify further:
∑
F
i
=
[
U
i
∂
ρ
∂
t
+
ρ
∂
U
i
∂
t
+
U
i
∂
∂
x
j
(
ρ
U
j
)
+
ρ
U
j
∂
U
i
∂
x
j
]
d
V
{\displaystyle \displaystyle \sum F_{i}=\left[U_{i}{\frac {\partial \rho }{\partial t}}+\rho {\frac {\partial U_{i}}{\partial t}}+U_{i}{\frac {\partial }{\partial x_{j}}}(\rho \;U_{j})+\rho \;U_{j}\;{\frac {\partial U_{i}}{\partial x_{j}}}\right]dV}
∑
F
i
=
U
i
[
∂
ρ
∂
t
+
∂
∂
x
j
(
ρ
U
j
)
]
⏟
c
o
n
t
i
n
u
i
t
y
e
q
u
a
t
i
o
n
=
0
d
V
+
ρ
[
∂
U
i
∂
t
+
U
j
∂
U
i
∂
x
j
]
⏟
D
U
i
D
t
;
s
u
b
t
a
n
t
i
a
l
d
e
r
i
v
a
t
i
v
e
o
f
U
i
d
V
{\displaystyle \displaystyle \sum F_{i}=U_{i}\underbrace {\left[{\frac {\partial \rho }{\partial t}}\;+{\frac {\partial }{\partial x_{j}}}(\rho \;U_{j})\right]} _{continuity\ equation=0}dV\;+\rho \underbrace {\left[{\frac {\partial U_{i}}{\partial t}}\;+U_{j}{\frac {\partial U_{i}}{\partial x_{j}}}\right]} _{{\frac {DU_{i}}{Dt}};subtantial\ derivative\ of\ U_{i}}dV}
Hence
∑
F
i
=
ρ
D
U
i
D
t
d
V
{\displaystyle \displaystyle \sum F_{i}=\rho {\frac {DU_{i}}{Dt}}dV}
Let's look to the forces on the exposed on the diffrential control volume:
d
P
i
d
t
=
∑
F
i
=
d
F
b
o
d
y
i
+
d
F
s
u
r
f
a
c
e
i
=
ρ
d
U
i
d
t
d
V
{\displaystyle \displaystyle {\frac {dP_{i}}{dt}}=\sum F_{i}=dF_{body\;i}+dF_{surface\;i}=\rho {\frac {dU_{i}}{dt}}dV}
Here, only gravitational force is considered as a body force. Thus,
d
F
b
o
d
y
i
=
ρ
d
V
g
i
{\displaystyle \displaystyle dF_{body\;i}=\rho \;dV\;g_{i}}
Differential surface forces
Surface forces are the stresses acting on the control surfaces.
F
s
{\displaystyle \displaystyle F_{s}}
can be resolved into three components.
d
F
n
{\displaystyle \displaystyle dF_{n}}
is normal to dA.
d
F
t
{\displaystyle \displaystyle dF_{t}}
are tangent to dA:
σ
n
=
lim
d
A
→
0
d
F
n
d
A
{\displaystyle \displaystyle \sigma _{n}=\lim _{dA\rightarrow 0}{\frac {dF_{n}}{dA}}}
σ
t
=
lim
d
A
→
0
d
F
t
d
A
{\displaystyle \displaystyle \sigma _{t}=\lim _{dA\rightarrow 0}{\frac {dF_{t}}{dA}}}
σ
n
{\displaystyle \displaystyle \sigma _{n}}
is a normal stress whereas
σ
t
{\displaystyle \displaystyle \sigma _{t}}
is a shear stress. The shear stresses are also designated by
τ
{\displaystyle \displaystyle \tau }
.
Stresses on the surface of differential control volume
Thus, the surface forces are due to stresses on the surfaces of the control surface.
σ
i
j
=
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
{\displaystyle \displaystyle \sigma _{ij}=\left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{matrix}}\right]}
The positive stress directions
We define the positive direction for the stress as the positive coordinate direction on the surfaces (e.g. on ABCD) for which the outwards normal is in the positive coordinate direction . If the outward normal represents the negative direction (A'B'C'D'), then the stresses are considered positive if directed in the negative coordinate directions.
The stresses on the surface
(
σ
i
j
)
{\displaystyle \displaystyle (\sigma _{ij})}
are the sum of pressure plus the viscous stresses which arise from motion with velocity gradients:
σ
i
j
=
[
−
P
0
0
0
−
P
0
0
0
−
P
]
+
[
τ
11
τ
12
τ
13
τ
21
τ
22
τ
23
τ
31
τ
32
τ
33
]
=
[
−
P
+
τ
11
τ
12
τ
13
τ
21
−
P
+
τ
22
τ
23
τ
31
τ
32
−
P
+
τ
33
]
{\displaystyle \displaystyle \sigma _{ij}=\left[{\begin{matrix}-P&0&0\\0&-P&0\\0&0&-P\end{matrix}}\right]+\left[{\begin{matrix}\tau _{11}&\tau _{12}&\tau _{13}\\\tau _{21}&\tau _{22}&\tau _{23}\\\tau _{31}&\tau _{32}&\tau _{33}\end{matrix}}\right]=\left[{\begin{matrix}-P+\tau _{11}&\tau _{12}&\tau _{13}\\\tau _{21}&-P+\tau _{22}&\tau _{23}\\\tau _{31}&\tau _{32}&-P+\tau _{33}\end{matrix}}\right]}
p
{\displaystyle \displaystyle p}
has a minus sign since the force due to pressure acts opposite to the surface normal.
Stresses on the surface of the differential control volume in the x1 direction
Let us look to the differential surface force in the
x
1
{\displaystyle \displaystyle x_{1}}
direction:
d
F
s
u
r
f
a
c
e
1
=
∂
σ
11
∂
x
1
d
x
1
d
x
2
d
x
3
+
∂
σ
21
∂
x
2
d
x
1
d
x
2
d
x
3
+
∂
σ
31
∂
x
3
d
x
1
d
x
2
d
x
3
{\displaystyle \displaystyle dF_{surface\ 1}={\frac {\partial \sigma _{11}}{\partial x_{1}}}\;dx_{1}dx_{2}dx_{3}+{\frac {\partial \sigma _{21}}{\partial x_{2}}}\;dx_{1}dx_{2}dx_{3}+{\frac {\partial \sigma _{31}}{\partial x_{3}}}\;dx_{1}dx_{2}dx_{3}}
Noting that
d
V
=
d
x
1
d
x
2
d
x
3
{\displaystyle \displaystyle dV=dx_{1}dx_{2}dx_{3}}
and
σ
i
j
=
−
p
δ
i
j
+
τ
i
j
{\displaystyle \displaystyle \sigma _{ij}=-p\delta _{ij}+\tau _{ij}}
(4),
d
F
s
u
r
f
a
c
e
1
=
(
−
∂
P
∂
x
1
+
∂
τ
11
∂
x
1
+
∂
τ
21
∂
x
2
+
∂
τ
31
∂
x
3
)
d
V
{\displaystyle \displaystyle dF_{surface\ 1}=\left(-{\frac {\partial P}{\partial x_{1}}}+{\frac {\partial \tau _{11}}{\partial x_{1}}}+{\frac {\partial \tau _{21}}{\partial x_{2}}}+{\frac {\partial \tau _{31}}{\partial x_{3}}}\right)dV}
Thus in tensor form the differential surface forces in
i
{\displaystyle \displaystyle i}
'th direction can be written as
d
F
s
u
r
f
a
c
e
i
=
(
−
∂
P
∂
x
i
+
∂
τ
j
i
∂
x
j
)
d
V
{\displaystyle \displaystyle dF_{surface\ i}=\left(-{\frac {\partial P}{\partial x_{i}}}+{\frac {\partial \tau _{ji}}{\partial x_{j}}}\right)dV}
Note that
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
is a symmetric tensor, i.e.
τ
j
i
=
τ
i
j
{\displaystyle \displaystyle \tau _{ji}=\tau _{ij}}
Hence, the diffential surface forces reads:
d
F
s
u
r
f
a
c
e
i
=
(
−
∂
P
∂
x
i
+
∂
τ
i
j
∂
x
j
)
d
V
{\displaystyle \displaystyle dF_{surface\ i}=\left(-{\frac {\partial P}{\partial x_{i}}}+{\frac {\partial \tau _{ij}}{\partial x_{j}}}\right)dV}
Inserting
d
F
b
o
d
y
i
{\displaystyle \displaystyle dF_{body\ i}}
and
d
F
s
u
r
f
a
c
e
i
{\displaystyle \displaystyle dF_{surface\ i}}
into (2),
ρ
D
U
i
D
t
d
V
=
ρ
g
i
d
V
+
(
−
∂
P
∂
x
i
+
τ
i
j
∂
x
j
)
d
V
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}dV=\rho \ g_{i}\ dV+\left(-{\frac {\partial P}{\partial x_{i}}}+{\frac {\tau _{ij}}{\partial x_{j}}}\right)\;dV}
and canceling
d
V
{\displaystyle \displaystyle dV}
we obtain
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
+
∂
τ
i
j
∂
x
j
.
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho \;g_{i}-{\frac {\partial P}{\partial x_{i}}}+{\frac {\partial \tau _{ij}}{\partial x_{j}}}.}
Expanding the substatial derivative at the left hand side,
ρ
∂
U
i
∂
t
+
ρ
U
j
∂
U
i
∂
x
j
=
ρ
g
i
−
∂
P
∂
x
i
+
∂
τ
i
j
∂
x
j
{\displaystyle \displaystyle \rho {\frac {\partial U_{i}}{\partial t}}+\rho U_{j}{\frac {\partial U_{i}}{\partial x_{j}}}=\rho g_{i}-{\frac {\partial P}{\partial x_{i}}}+{\frac {\partial \tau _{ij}}{\partial x_{j}}}}
We obtain the the most general form of momentum equation which is valid for any fluid (Newtonian, Non-newtonian, Compressible, etc.). It is non-linear due to the
2
n
d
{\displaystyle \displaystyle 2^{nd}}
term at the LHS. Efect of Newtonian and Non-newtonian properties appears in the formulation of the viscous stresses
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
.
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
will introduce also non-linearity when the fluid is non-Newtonian.
It should be noted that these formulations are based on stress conception which was thought to exist in fluids in motion. However it is known that
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
can be expressed as momentum transfer per unit area and time. Thus it can be considered as molecular momentum transport term. Derivations based on this concept requires a molecular approach (which is lengthy). The students should be aware that
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
causes momentum transport when there is a gradient of velocity.
Linear momentum equation for Newtonian Fluid: "Navier-Stokes Equation"
edit
For a Newtonian fluid, the viscous stresses are defined as:
τ
i
j
=
μ
[
∂
U
i
∂
x
j
+
∂
U
j
∂
x
i
]
−
2
3
δ
i
j
μ
∂
U
k
∂
x
k
{\displaystyle \displaystyle \tau _{ij}=\mu \left[{\frac {\partial U_{i}}{\partial x_{j}}}+{\frac {\partial U_{j}}{\partial x_{i}}}\right]-{\frac {2}{3}}\delta _{ij}\mu {\frac {\partial U_{k}}{\partial x_{k}}}}
Note that derivation of this relation is beyond the scaope of this course.
Thus, the momentum equation (6) becomes
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
+
∂
∂
x
j
[
μ
(
∂
U
i
∂
x
j
+
∂
U
j
∂
x
i
)
−
2
3
μ
δ
i
j
∂
U
k
∂
x
k
]
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho g_{i}-{\frac {\partial P}{\partial x_{i}}}+{\frac {\partial }{\partial x_{j}}}\left[\mu \left({\frac {\partial U_{i}}{\partial x_{j}}}+{\frac {\partial U_{j}}{\partial x_{i}}}\right)-{\frac {2}{3}}\mu \ \delta _{ij}\ {\frac {\partial U_{k}}{\partial x_{k}}}\right]}
For a flow with constant viscosity (
μ
=
constant
{\displaystyle \displaystyle \mu ={\textrm {constant}}}
):
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
+
μ
∂
∂
x
j
[
∂
U
i
∂
x
j
+
∂
U
j
∂
x
i
−
2
3
δ
i
j
∂
U
k
∂
x
k
]
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho g_{i}-{\frac {\partial P}{\partial x_{i}}}+\mu {\frac {\partial }{\partial x_{j}}}\left[{\frac {\partial U_{i}}{\partial x_{j}}}+{\frac {\partial U_{j}}{\partial x_{i}}}-{\frac {2}{3}}\delta _{ij}\ {\frac {\partial U_{k}}{\partial x_{k}}}\right]}
since,
∂
2
U
i
∂
x
j
∂
x
i
=
∂
2
U
i
∂
x
i
∂
x
j
=
∂
∂
x
j
(
∂
U
i
∂
x
i
)
=
∂
∂
x
i
(
∂
U
i
∂
x
j
)
{\displaystyle \displaystyle {\frac {\partial ^{2}U_{i}}{\partial x_{j}\partial x_{i}}}={\frac {\partial ^{2}U_{i}}{\partial x_{i}\partial x_{j}}}={\frac {\partial }{\partial x_{j}}}\left({\frac {\partial U_{i}}{\partial x_{i}}}\right)={\frac {\partial }{\partial x_{i}}}\left({\frac {\partial U_{i}}{\partial x_{j}}}\right)}
then,
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
+
μ
∂
2
U
i
∂
x
j
x
j
+
μ
∂
∂
x
i
(
∂
U
j
∂
x
j
)
−
2
3
δ
i
j
μ
∂
∂
x
j
(
∂
U
k
∂
x
k
)
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho \ g_{i}-{\frac {\partial P}{\partial x_{i}}}+\mu {\frac {\partial ^{2}U_{i}}{\partial x_{j}x_{j}}}+\mu {\frac {\partial }{\partial x_{i}}}\left({\frac {\partial U_{j}}{\partial x_{j}}}\right)-{\frac {2}{3}}\delta _{ij}\mu {\frac {\partial }{\partial x_{j}}}\left({\frac {\partial U_{k}}{\partial x_{k}}}\right)}
For an incompressible flow
∂
U
k
∂
x
k
=
∂
U
j
∂
x
j
=
0
{\displaystyle \displaystyle {\frac {\partial U_{k}}{\partial x_{k}}}={\frac {\partial U_{j}}{\partial x_{j}}}=0}
hence assuming that the viscosity is constant, it can be easily shown that the momentum equation reduces to
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
+
μ
∂
2
U
i
∂
x
j
2
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho \ g_{i}-{\frac {\partial P}{\partial x_{i}}}+\mu {\frac {\partial ^{2}U_{i}}{\partial x_{j}^{2}}}}
Euler's equation: Inviscid flow
edit
When the velocity gradients in the flow is negligible and/or the Reynolds number takes very high values, the viscous stresses can be neglected:
τ
i
j
=
0
{\displaystyle \displaystyle \displaystyle \tau _{ij}=0}
Since, the viscous stresses are proportional to viscosity:
τ
i
j
∝
μ
{\displaystyle \displaystyle \tau _{ij}\propto \mu }
for flows, where
τ
i
j
{\displaystyle \displaystyle \tau _{ij}}
is neglected, the flow is called frictionless or inviscid, although there is a finite viscosity of the flow. Accordingly, the linear momentum equation reduces to
ρ
D
U
i
D
t
=
ρ
g
i
−
∂
P
∂
x
i
.
{\displaystyle \displaystyle \rho {\frac {DU_{i}}{Dt}}=\rho g_{i}-{\frac {\partial P}{\partial x_{i}}}.}
Euler's equation in streamline coordinates
edit
Differential control volume along streamline coordinates and the forces on it for a inviscid flow
Differential control volume along streamline coordinates and the forces on it for a inviscid flow
Euler's equation take a special form along and normal to a streamline with which one can see the dependency between the pressure, velocity and curvature of the streamline.
To obtain Euler's equation in s-direction, apply Newton's second law in s-direction in the absence of viscous forces.
ρ
d
V
[
∂
U
s
∂
t
+
U
s
∂
U
s
∂
s
]
=
−
∂
P
∂
s
d
V
−
ρ
g
s
i
n
β
d
V
{\displaystyle \displaystyle \rho dV\left[{\frac {\partial U_{s}}{\partial t}}+U_{s}{\frac {\partial U_{s}}{\partial s}}\right]=-{\frac {\partial P}{\partial s}}dV-\rho \ gsin\beta dV}
Omitting
d
V
{\displaystyle \displaystyle dV}
would deliver
D
U
s
D
t
=
−
1
ρ
∂
P
∂
s
−
g
s
i
n
β
{\displaystyle \displaystyle {\frac {DU_{s}}{Dt}}=-{\frac {1}{\rho }}{\frac {\partial P}{\partial s}}-gsin\beta }
Since
s
i
n
β
≈
d
x
2
d
s
=
∂
x
2
∂
s
{\displaystyle \displaystyle sin\beta \approx {\frac {dx_{2}}{ds}}={\frac {\partial x_{2}}{\partial s}}}
then the Equler's equation along a streamline reads
D
U
s
D
t
=
−
1
ρ
∂
P
∂
s
−
g
∂
x
2
∂
s
(
8
)
{\displaystyle \displaystyle {\frac {DU_{s}}{Dt}}=-{\frac {1}{\rho }}{\frac {\partial P}{\partial s}}-g{\frac {\partial x_{2}}{\partial s}}\ \ \ (8)}
For a steady flow and by neglecting body forces,
1
ρ
∂
P
∂
s
=
−
U
s
∂
U
s
∂
s
{\displaystyle \displaystyle {\frac {1}{\rho }}{\frac {\partial P}{\partial s}}=-U_{s}{\frac {\partial U_{s}}{\partial s}}}
it can be seen that decrease in velocity means increase in pressure as is indicated by the Bernoulli equation.
To obtain Euler's equation in n direction, apply Newton's second law in the absence of viscous forces and for a steady flow.
ρ
d
V
D
U
n
D
t
=
(
P
−
∂
P
∂
n
d
n
2
)
d
s
d
x
3
−
(
P
+
∂
P
∂
n
d
n
2
)
d
s
d
x
3
−
ρ
g
d
V
c
o
s
β
{\displaystyle \displaystyle \rho dV{\frac {DU_{n}}{Dt}}=\left(P-{\frac {\partial P}{\partial n}}{\frac {dn}{2}}\right)ds\ dx_{3}-\left(P+{\frac {\partial P}{\partial n}}{\frac {dn}{2}}\right)ds\ dx_{3}-\rho gdVcos\beta }
ρ
D
U
n
D
t
=
(
−
∂
P
∂
n
−
ρ
g
c
o
s
β
)
{\displaystyle \displaystyle \rho \ {\frac {DU_{n}}{Dt}}=\left(-{\frac {\partial P}{\partial n}}-\rho \ g\ cos\beta \right)}
Since,
c
o
s
β
≈
d
x
2
d
n
=
∂
x
2
∂
n
{\displaystyle \displaystyle cos\beta \approx {\frac {dx_{2}}{dn}}={\frac {\partial x_{2}}{\partial n}}}
Then,
D
U
n
D
t
=
−
1
ρ
∂
P
∂
n
−
g
∂
x
2
∂
n
{\displaystyle \displaystyle {\frac {DU_{n}}{Dt}}=-{\frac {1}{\rho }}{\frac {\partial P}{\partial n}}-g\ {\frac {\partial x_{2}}{\partial n}}}
For a steady flow, the normal acceleration of the fluid is towards the center of curvature of the streamline:
D
U
n
D
t
=
−
U
s
2
R
{\displaystyle \displaystyle {\frac {DU_{n}}{Dt}}=-{\frac {U_{s}^{2}}{R}}}
Hence,
U
s
2
R
=
1
ρ
∂
P
∂
n
+
g
∂
x
2
∂
n
{\displaystyle \displaystyle {\frac {U_{s}^{2}}{R}}={\frac {1}{\rho }}{\frac {\partial P}{\partial n}}+g{\frac {\partial x_{2}}{\partial n}}}
For an unsteady flow,
D
U
n
D
t
=
−
U
s
2
R
+
∂
U
n
∂
t
{\displaystyle \displaystyle {\frac {DU_{n}}{Dt}}=-{\frac {U_{s}^{2}}{R}}+{\frac {\partial U_{n}}{\partial t}}}
For steady flow neglecting body forces, the Euler's equation normal to the streamline is
1
ρ
∂
P
∂
n
=
U
s
2
R
{\displaystyle \displaystyle {\frac {1}{\rho }}{\frac {\partial P}{\partial n}}={\frac {U_{s}^{2}}{R}}}
which indicates that pressure increases in a direction outwards from the center of the curvature of the streamlines. In other words, pressure drops towards the center of curvature, which, consequently creates a potential difference in terms of pressure and forces the fluid to change its direction. For a straight streamline
R
→
∞
{\displaystyle \displaystyle R\rightarrow \infty }
, there is no pressure variation normal to the streamline.
Bernoulli equation: Integration of Euler's equation along a streamline for a steady flow
edit
For a steady flow, Euler's equation along a streamline reads,
1
ρ
∂
P
∂
s
+
g
∂
x
2
∂
s
=
−
U
s
∂
U
s
∂
s
{\displaystyle \displaystyle {\frac {1}{\rho }}{\frac {\partial P}{\partial s}}+g{\frac {\partial x_{2}}{\partial s}}=-U_{s}{\frac {\partial U_{s}}{\partial s}}}
If a fluid particle moves a distance ds, along a streamline, since every variable becomes a function of
s
{\displaystyle \displaystyle s}
:
∂
P
∂
s
d
s
=
d
P
:Change in pressure along s,
{\displaystyle \displaystyle {\frac {\partial P}{\partial s}}ds=dP{\text{:Change in pressure along s,}}}
∂
x
2
∂
s
d
s
=
d
x
2
:Change in elevation along s,
{\displaystyle \displaystyle {\frac {\partial x_{2}}{\partial s}}ds=dx_{2}\ {\text{:Change in elevation along s,}}}
∂
U
s
∂
s
d
s
=
d
U
s
:Change in velocity along s.
{\displaystyle \displaystyle {\frac {\partial U_{s}}{\partial s}}ds=dU_{s}\ {\text{:Change in velocity along s.}}}
Integration of the Euler equation between two locations, 1 and 2, along
s
{\displaystyle \displaystyle s}
reads
∫
1
2
(
1
ρ
∂
P
∂
s
+
g
∂
x
2
∂
s
+
U
s
∂
U
s
∂
s
)
d
s
=
0
{\displaystyle \displaystyle \int _{1}^{2}\left({{\frac {1}{\rho }}{\frac {\partial P}{\partial s}}+g{\frac {\partial x_{2}}{\partial s}}+U_{s}{\frac {\partial U_{s}}{\partial s}}}\right)ds=0}
For incompressible flow
ρ
=
constant
{\displaystyle \displaystyle \rho ={\textrm {constant}}}
and after changing the notation as:
U
=
U
s
{\displaystyle \displaystyle U=U_{s}}
and
z
=
x
2
{\displaystyle \displaystyle z=x_{2}}
, the integration results in
P
2
−
P
1
ρ
+
g
(
z
2
−
z
1
)
+
U
2
2
−
U
1
2
2
=
0
\ along\ s
{\displaystyle \displaystyle {\frac {P_{2}-P_{1}}{\rho }}+g(z_{2}-z_{1})+{\frac {U_{2}^{2}-U_{1}^{2}}{2}}=0{\textrm {\ along\ s}}}
or in its most beloved form:
P
1
ρ
+
g
z
1
+
U
1
2
2
=
P
2
ρ
+
g
z
2
+
U
2
2
2
{\displaystyle \displaystyle {\frac {P_{1}}{\rho }}+gz_{1}+{\frac {U_{1}^{2}}{2}}={\frac {P_{2}}{\rho }}+gz_{2}+{\frac {U_{2}^{2}}{2}}}
In other words along a streamline:
P
ρ
+
g
z
+
U
2
2
=
constant
{\displaystyle \displaystyle {\frac {P}{\rho }}+gz+{\frac {U^{2}}{2}}={\textrm {constant}}}
Note that due to the assumptions made during the derivation, the following restrictions applies to this equation: The flow should be steady, incompressible, frictionless and the equation is valid only along a streamline.
The common forms of Bernoulli equation are as follows:
Energy form (per unit mass)
U
2
2
⏟
kinetic energy
+
P
ρ
⏟
pressure energy
+
g
z
⏟
potential energy
=
ζ
{\displaystyle \displaystyle \underbrace {\frac {U^{2}}{2}} _{\text{kinetic energy}}+\underbrace {\frac {P}{\rho }} _{\text{pressure energy}}+\underbrace {gz} _{\text{potential energy}}=\zeta }
Pressure form
ρ
U
2
2
⏟
dynamic pressure
+
P
⏟
static pressure
+
ρ
g
z
⏟
geodesic pressure
=
K
{\displaystyle \displaystyle \underbrace {\rho {\frac {U^{2}}{2}}} _{\text{dynamic pressure}}+\underbrace {P} _{\text{static pressure}}+\underbrace {\rho gz} _{\text{geodesic pressure}}=K}
Head form
U
2
2
g
⏟
velocity head
+
P
ρ
g
⏟
pressure head
+
z
⏟
geodesic head
=
k
{\displaystyle \displaystyle \underbrace {\frac {U^{2}}{2g}} _{\text{velocity head}}+\underbrace {\frac {P}{\rho g}} _{\text{pressure head}}+\underbrace {z} _{\text{geodesic head}}=k}
Static, stagnation and dynamic pressures
edit
How do we measure pressure? When the streamlines are parallel to the wall we can use pressure taps.
If the measured location is far from the wall, static pressure measurements can be made by a static pressure probe.
The stagnation pressure is the value obtained when a flowing fluid is decelerated to zero velocity by a frictionless flow process. The Stagnation pressure can be calculated as follows:
P
ρ
+
U
2
2
=
constant
{\displaystyle \displaystyle {\frac {P}{\rho }}+{\frac {U^{2}}{2}}=\;{\textrm {constant}}}
P
0
ρ
+
U
0
2
2
=
P
1
ρ
+
U
1
2
2
{\displaystyle \displaystyle {\frac {P_{0}}{\rho }}+{\frac {U_{0}^{2}}{2}}={\frac {P_{1}}{\rho }}+{\frac {U_{1}^{2}}{2}}}
when
U
0
2
=
0
{\displaystyle \displaystyle U_{0}^{2}=0}
P
0
=
P
1
+
ρ
U
1
2
2
,
{\displaystyle \displaystyle P_{0}=P_{1}+\rho {\frac {U_{1}^{2}}{2}},}
where the last term is the dynamic pressure.
Pitot tube used in racing car
Figure: a)Jet is impinging to the wall and stagnating at the point of impingment b)Schematics of a of a pressure tap in a channel c)Static pressure probe
If we know the pressure difference
P
0
−
P
1
{\displaystyle \displaystyle P_{0}-P_{1}}
, we can calculate the
U
1
{\displaystyle \displaystyle U_{1}}
velocity.
U
1
=
2
(
P
0
−
P
1
)
ρ
{\displaystyle \displaystyle U_{1}={\sqrt {\frac {2\left(P_{0}-P_{1}\right)}{\rho }}}}
The stagnation pressure is measured in the laboratory using a probe that faces directly upstream flow.
Such a probe is called a stagnation pressure probe or Pitot tube . Thus, using a pressure tap and a Pitot tube one can measure local velocity:
P
0
=
P
A
+
ρ
U
A
2
{\displaystyle \displaystyle P_{0}=P_{A}+\rho {\frac {U_{A}}{2}}}
P
0
−
P
A
=
(
P
A
+
ρ
U
A
2
2
)
−
(
P
A
)
=
ρ
U
A
2
2
{\displaystyle \displaystyle P_{0}-P_{A}=\left(P_{A}+\rho {\frac {U_{A}^{2}}{2}}\right)-(P_{A})=\rho {\frac {U_{A}^{2}}{2}}}
Thus, measuring
P
0
−
P
A
{\displaystyle \displaystyle P_{0}-P_{A}}
one can determine
U
A
{\displaystyle \displaystyle U_{A}}
.
Pitot-static tube for velocity measurement
However, in the absence of a wall with well defined location, the velocity can be measured by a Pitot-static tube . The pressure is measured at B and C; assuming
P
B
=
P
C
{\displaystyle \displaystyle P_{B}=P_{C}}
.
Hence,
U
B
=
2
(
P
0
B
−
P
B
)
ρ
.
{\displaystyle \displaystyle U_{B}={\sqrt {\frac {2(P_{0_{B}}-P_{B})}{\rho }}}.}
Unsteady Bernoulli equation
edit
The Euler's equation along a streamline is:
−
1
ρ
∂
P
∂
s
−
g
∂
z
∂
s
=
D
U
s
D
t
=
U
s
∂
U
s
∂
s
+
∂
U
s
∂
t
{\displaystyle \displaystyle -{\frac {1}{\rho }}{\frac {\partial P}{\partial s}}-g{\frac {\partial z}{\partial s}}={\frac {DU_{s}}{Dt}}=U_{s}{\frac {\partial U_{s}}{\partial s}}+{\frac {\partial U_{s}}{\partial t}}}
along ds,
−
1
ρ
∂
P
∂
s
d
s
−
g
∂
z
∂
s
d
s
=
U
s
∂
U
s
∂
s
d
s
+
∂
U
s
∂
t
d
s
{\displaystyle \displaystyle -{\frac {1}{\rho }}{\frac {\partial P}{\partial s}}ds-g{\frac {\partial z}{\partial s}}ds=U_{s}{\frac {\partial U_{s}}{\partial s}}ds+{\frac {\partial U_{s}}{\partial t}}ds}
hence,
−
1
ρ
d
P
−
g
d
z
=
U
s
d
U
s
+
∂
U
s
∂
t
d
s
.
{\displaystyle \displaystyle -{\frac {1}{\rho }}dP-gdz=U_{s}dU_{s}+{\frac {\partial U_{s}}{\partial t}}ds.}
Integration between two points along a streamline is:
−
∫
1
2
d
P
ρ
−
∫
1
2
g
d
z
=
∫
1
2
U
s
d
U
s
+
∫
1
2
∂
U
s
∂
t
d
s
{\displaystyle \displaystyle -\int _{1}^{2}{\frac {dP}{\rho }}-\int _{1}^{2}{g}dz=\int _{1}^{2}{U_{s}dU_{s}}+\int _{1}^{2}{\frac {\partial U_{s}}{\partial t}}ds}
For incompressible flow,
ρ
=
constant
{\displaystyle \displaystyle \rho ={\textrm {constant}}}
, thus the integral reads
P
1
−
P
2
ρ
+
g
(
z
1
−
z
2
)
=
U
2
2
−
U
1
2
2
+
∫
1
2
∂
U
s
∂
t
d
s
{\displaystyle \displaystyle {\frac {P_{1}-P_{2}}{\rho }}+g(z_{1}-z_{2})={\frac {U_{2}^{2}-U_{1}^{2}}{2}}+\int _{1}^{2}{{\frac {\partial U_{s}}{\partial t}}ds}}
The unsetady Bernoulli equation involves the integration of the time gradient of the velocity between two points.:
P
1
ρ
+
g
z
1
+
U
1
2
2
=
P
2
ρ
+
g
z
2
+
U
2
2
2
+
∫
1
2
∂
U
s
∂
t
d
s
{\displaystyle \displaystyle {\frac {P_{1}}{\rho }}+gz_{1}+{\frac {U_{1}^{2}}{2}}={\frac {P_{2}}{\rho }}+gz_{2}+{\frac {U_{2}^{2}}{2}}+\int _{1}^{2}{{\frac {\partial U_{s}}{\partial t}}ds}}