Derivation of Continuity Equation in Cylindrical Coordinates
Derivation of Continuity Equation in Cylindrical Coordinates
CYLINDRICAL COORDINATES
V = vr r+ v + vz
z
where r , and z are the unit vectors in the r , and z direction
For the net flow through the control volume, we deal with it one face at
a time. Starting with the r faces, the net inflow is
m
r,in
= v r rd d z
r,out
=(v r +
r
r )(r+dr)ddz
So that the net flow in the r direction is Being 0,0(d2r), the last term in
this equation can be dropped so that the net flow on the r faces is
m
out
in
r
dr
)(r+dr)ddz
r
v rrd d z
=(
r
v r + r
r)drddz
The net flow in the direction is slightly easier to compute since the
areas of the inflow and outflow faces are the same. At the outset, the net
flow in the theta direction is
m
,in
= v d r d z
,out
=(v +
d) dr dz
out
- m
in
d)drdz
v d rd z
ddrdz
We now turn our attention to the z direction. The face area is that of
a sector of angle
1
2
Az =
( r + dr )
= r dr d +
1
2
1
2
d2r d
= r dr d
then, the inflow at the lower z face is
m
z,in
= v z r d r d
z,out
=(v z +
z
z dz) r dr d
out
= (v z +
- m
in
z dz) r dr d
v z r d r d
z
z
r dr ddz
r dr ddz
r dr ddz + ( v r +
r
r
r)drddz +
ddrdz +
z
z
dr ddz
By dividing ( r dr ddz )
=
1
r
1
r
r
v r + r
r r
r
1
r
1
r
z
z
z
z
=
z
z
1
r
r r
r
1
r