Scribd
Upload
57 views2 pages

Continuity Equation in Other Coordinate Systems: An Internet Book On Fluid Dynamics

The document discusses the continuity equation in various coordinate systems. It shows that in rectangular coordinates the general continuity equation is ∂ρ/∂t + ∇·(ρu) = 0. For planar flows the continuity equation simplifies to ∂ρ/∂t + ∂(ρu)/∂x + ∂(ρv)/∂y = 0. In polar coordinates the continuity equation is ∂ρ/∂t + (1/r)(∂(ρrur)/∂r) + (1/r)(∂(ρuθ)/∂θ) = 0. For cylindrical and spherical coordinates, the document provides the general continuity equations and examples of

Uploaded by

Krishna Jaiswal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
57 views2 pages

Continuity Equation in Other Coordinate Systems: An Internet Book On Fluid Dynamics

The document discusses the continuity equation in various coordinate systems. It shows that in rectangular coordinates the general continuity equation is ∂ρ/∂t + ∇·(ρu) = 0. For planar flows the continuity equation simplifies to ∂ρ/∂t + ∂(ρu)/∂x + ∂(ρv)/∂y = 0. In polar coordinates the continuity equation is ∂ρ/∂t + (1/r)(∂(ρrur)/∂r) + (1/r)(∂(ρuθ)/∂θ) = 0. For cylindrical and spherical coordinates, the document provides the general continuity equations and examples of

Uploaded by

Krishna Jaiswal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

An Internet Book on Fluid Dynamics

Continuity equation in other coordinate systems

We recall that in a rectangular Cartesian coordinate system the general continuity equation is
∂ρ ∂(ρu) ∂(ρv) ∂(ρw)
+ + + =0 (Bce1)
∂t ∂x ∂y ∂z
or in tensor notation
∂ρ ∂(ρuj )
+ =0 (Bce2)
∂t ∂xj
or in vector notation
∂ρ
+ ∇.(ρu) = 0 (Bce3)
∂t

If the flow is planar, the velocity and the derivatives in one direction (say the z-direction) become zero
and the continuity equation reduces to
∂ρ ∂(ρu) ∂(ρv)
+ + =0 (Bce4)
∂t ∂x ∂y
and if the flow is incompressible this is further reduced to
∂u ∂v
+ =0 (Bce5)
∂x ∂y
a form that is repeatedly used in this text.

In a planar flow such as this it is sometimes convenient to use a polar coordinate system (r, θ). Then
the continuity equation becomes
∂ρ 1 ∂(ρrur ) 1 ∂(ρuθ )
+ + =0 (Bce6)
∂t r ∂r r ∂θ
where ur , uθ are the velocities in the r and θ directions. If the fluid is incompressible this further reduces
to
∂(rur ) ∂uθ
+ =0 (Bce7)
∂r ∂θ
and, in the very simple circumstances in which the flow is purely radial so that uθ = 0, it follows that
∂(rur )
=0 (Bce8)
∂r
and therefore ur = C/r where C is a constant or only a function of time. This is the flow field for a simple
”line source or sink” at r = 0 which is perpendicular to the plane of r and θ. We will discuss these flows
elsewhere in this book. Another simple planar flow that is valuable to identify is one in which the vortical

flow in which ur is everywhere zero and the flow proceeds in a circular path. Then the continuity equation
for incompressible flow becomes
∂uθ
=0 (Bce9)
∂θ
and therefore uθ is a function only of r and t. There are several specific types of vortical motion that
are discussed elsewhere in this book. Reverting to the more general three-dimensional flow, the continuity

equation in cylindrical coordinates (r, θ, z) is


∂ρ 1 ∂(ρrur ) 1 ∂(ρuθ ) ∂(ρuz )
+ + + =0 (Bce10)
∂t r ∂r r ∂θ ∂z
where ur , uθ , uz are the velocities in the r, θ and z directions of the cylindrical coordinate system. A
particular subset of such flows is axisymmetric flow in which the derivatives in the θ direction are zero so
that the continuity equation becomes
∂ρ 1 ∂(ρrur ) ∂(ρuz )
+ + =0 (Bce11)
∂t r ∂r ∂z
where ur , uθ , uz are the velocities in the r, θ and z directions of the cylindrical coordinate system. For
incompressible axisymmetric flow it follows that the continuity equation becomes
1 ∂(rur ) ∂(uz )
+ =0 (Bce12)
r ∂r ∂z
Another subset is vortical flow about the z axis (with no flow in the z direction) in which case the continuity
equation becomes
∂ρ 1 ∂(ρrur ) 1 ∂(ρuθ )
+ + =0 (Bce13)
∂t r ∂r r ∂θ

Turning now to spherical coordinates (r, θ, φ) the continuity equation becomes


∂ρ 1 ∂(ρr2 ur ) 1 ∂(ρuθ sin θ) 1 ∂(ρuφ )
+ 2 + + =0 (Bce14)
∂t r ∂r r sin θ ∂θ r sin θ ∂φ
and in some particular flows in which this simplifies further. In purely radial flow such as that due to a
point source or a sink (as opposed to the line source or sink described above) it becomes
∂ρ 1 ∂(ρr2ur )
+ 2 =0 (Bce15)
∂t r ∂r
For incompressible flow this reduces to
∂(r2ur )
=0 (Bce16)
∂r
and therefore ur = C/r2 where C is only a function of time. Clearly this can be arrived at by arguing
that for an incompressible fluid the spherical surface area times the radial velocity must be the same at
all radii.

You might also like