<< Chapter < Page Chapter >> Page >

Conservation of mass and the equation of continuity

Although the idea of mass is not a kinematical one, it is convenient to introduce it here and to obtain the continuity equation. Let ρ ( x , t ) be the mass per unit volume of a homogeneous fluid at position x and time t . Then the mass of any finite material volume V ( t ) is

m = v ( t ) ρ ( x , t ) d V .

If V is a material volume, that is, if it is composed of the same particles, and there are no sources or sinks in the medium we take it as a principle that the mass does not change. By inserting = ρ in Reynolds' transport theorem we have

d m d t = v ( t ) D ρ D t + ρ ( v ) d V = v ( t ) ρ t + ( ρ v ) d V = 0

This is true for an arbitrary material volume and hence the integrand itself must vanish everywhere. It follows that

D ρ D t + ρ ( v ) = ρ t + ( ρ v ) = 0

which is the equation of continuity.

A fluid for which the density ρ is constant is called incompressible . In this case the equation of continuity becomes

v = 0 incompressible flow

and the motion is isochoric or the velocity field solenoidal.

Combining the equation of continuity with Reynolds' transport theorem for a function = ρ F we have

d d t v ( t ) ρ F d V = v ( t ) D D t ( ρ F ) + ρ F ( v ) d V = v ( t ) ρ D F D t + F ( D ρ D t + ρ v ) d V = v ( t ) ρ D F D t d V

This equation is useful for deriving the conservation equation of a quantity that is expressed as specific to a unit of mass, e.g., specific internal energy and species mass fraction.

Deformation and rate of strain

The motion of fluids differs from that of rigid bodies in the deformation or strain that occurs with motion. Material coordinates give a reference frame for this deformation or strain.

Consider two nearby points P and Q with material coordinates ξ and ξ + d ξ . At time t they are to be found at x ( ξ , t ) and x ( ξ + d ξ , t ) . Now

x i ( ξ + d ξ , t ) = x i ( ξ , t ) + x i ξ j d ξ j + O ( d 2 )

where O ( d 2 ) represents terms of order d ξ 2 and higher which will be neglected from this point onward. Thus the small displacement vector d ξ has now become

d x = x ( ξ + d ξ , t ) - x ( ξ , t )

where

d x i = x i ξ j d ξ j .

It is clear from the quotient rule (since d ξ is arbitrary) that the nine quantities x i / ξ j are the components of a tensor. It may be called the displacement gradient tensor and is basic to the theories of elasticity and rheology. For fluid motion, its material derivative is of more direct application and we will concentrate on this.

If v = D x / D t is the velocity, the relative velocity of two particles ξ and ξ + d ξ has components

d v i = v i ξ k d ξ k = D D t x i ξ k d ξ k

However, by inverting the above relation, we have

d v i = v i ξ k ξ k x j d x j = v i x j d x j

expressing the relative velocity in terms of current position. Again it is evident that the ( v i / x j ) are components of a tensor, the velocity gradient tensor , for which we need to obtain a sound physical feeling.

We first observe that if the motion is a rigid body translation with a constant velocity u ,

x = ξ + u t

and the velocity gradient tensor vanishes identically. Secondly, the velocity gradient tensor can be written as the sum of symmetric and antisymmetric parts,

v i x j 1 2 v i x j + v j x i + 1 2 v i x j - v j x i = e i j + Ω i j or v = e + Ω

We have seen that a relative velocity d v i related to the relative position d x j by an antisymmetric tensor Ω i j , i.e., d v i = Ω i j d x j , represents a rigid body rotation with angular velocity ω = - v e c Ω . In this case

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Transport phenomena. OpenStax CNX. May 24, 2010 Download for free at http://cnx.org/content/col11205/1.1
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Transport phenomena' conversation and receive update notifications?

Ask