differential
differential
Assuming the continuum approximation is valid, we then consider the flux of mass into the differential
control volume and equate it with the rate of increase of mass inside the control volume. Consider first
the flux of mass through the two sides perpendicular to the x−axis. We will define u as the velocity in
the x−direction at the center of lefthand of these two sides. Similarly we define the density of the fluid
at the center of lefthand side as ρ. Then the flux of mass into the control volume through that lefthand
side is given by ρudy dz. Then it follows that by Taylor’s series (neglecting all terms of order (dx)2 and
higher which can be shown to have no contribution to the result) the flux of mass out of the control volume
through the righthand side is given by [ρu + {∂(ρu)/∂x}dx]dy dz. Combining the fluxes through these
two sides perpendicular to the x-direction, it follows that the net flux of mass out of the control volume
through the sides perpendicular to the x-direction is {∂(ρu)/∂x}dx dy dz. The fluxes through the other
two pairs of sides follow from a similar construction so that the net flux of mass out of the control volume
through all of its sides becomes
∂(ρu) ∂(ρv) ∂(ρw) ∂(ρuj )
+ + dx dy dz = dx dy dz (Bcd1)
∂x ∂y ∂z ∂xj
By conservation of mass this must be equal to minus the rate of increase of mass inside the control volume.
Since the mass inside the control volume is ρ dx dy dz the differential form of the continuity equation
becomes
∂(ρu) ∂(ρv) ∂(ρw) ∂(ρuj ) ∂(ρ dx dy dz)
+ + dx dy dz = dx dy dz = − (Bcd2)
∂x ∂y ∂z ∂xj ∂t
where we must use the Eulerian time derivative since the control volume is defined as an Eulerian volume.
Re-arranging and cancelling the differential form of the continuity equation becomes
∂ρ ∂(ρu) ∂(ρv) ∂(ρw)
+ + + =0 (Bcd3)
∂t ∂x ∂y ∂z
or in tensor notation
∂ρ ∂(ρuj )
+ =0 (Bcd4)
∂t ∂xj
or in vector notation
∂ρ
+ ∇.(ρu) = 0 (Bcd5)
∂t
These are the differential forms of the continuity equation in a rectangular Cartesian coordinate system.
There are also many problems in which it is much more convenient to use an alternate coordinate system
such as a polar coordinate system, a cylindrical coordinate system or a spherical coordinate system. We
detail the forms of the continuity equation in these alternate coordinate systems on another page.