Chapter 10

Conservation of Mass: The Continuity Equation

10.1 Fluid Flow in the Bulk

10.1.1 Field Variables

As we will see, the properties of a fluid flow are described by the field variables. The field variables are the physical properties of the fluids that are usually required in order to adequately describe the behavior of the fluid. These variables are either scalar or vector quantities. The scalar field variables of interest are:

 pressure p

 density ρ

 temperature T

The only vector variable that needs to be considered is the velocity

v=(vxvyvz)

si2_e

with v=dsdtsi3_e and a=dvdtsi4_e. However, this vector has three components (along each axis of the chosen reference system). So in total, these amount to 6 variables that need to be found.

A total of six equations is required to find the 6 field variables. These are provided by the following fundamental laws of conservation

 Conservation of mass (continuity equation) - one equation

 Conservation of momentum (Navier-Stokes equation) - three equations (one for each coordinate axis)

 Conservation of energy - one equation

We will need a sixth equation to calculate all field variables. The problem becomes a little more complicated as the conservation of energy will introduce the internal energy U as an additional field variable. Usually, we are not interested in the internal energy. We therefore need a thermodynamic equation of state that links the internal energy U to the desired field variables temperature T , pressure p, density ρ, and velocity v. We will derive this in section 12.8. Introducing this equation allows us to reconvert the internal energy U into our field variables. However, we still need a sixth equation.

If the fluid in question is compressible (i.e., if it is a gas), this sixth equation is the ideal gas equation (see Eq. 6.4). If the fluid in question is a liquid (incompressible), we do not need a sixth equation, as the density ρ of liquids is assumed to be constant. This leaves only five field variables for which we have a sufficient number of equations.

10.1.2 Common Simplifications

Fluidmechanics is all about solving very complex partial differential equations. In many cases, allowing a certain term to be neglected can significantly simplify finding a solution often by only sacrificing a very small degree of precision. This is why simplifications should be made as often as feasible. Several simplifications are often used.

Fully Developed.A fully developed flow was a flow profile that does not change as one moves along a given axis. Consider the example of a fluid flowing out of a big reservoir into a small pipe. Right at the beginning of the pipe, the flow profile will not be fully developed yet, as there are effects from the abrupt reduction in diameter which change the flow. As the liquid moves along the pipe, the flow will eventually reach a certain form that will not change as long as the ambiance (pipe diameter, etc.) does not change. In this state, the flow is fully developed and there will be no changes to the flow profile along the axis variable x of the pipe. In this case, all (partial) derivatives with respect to x will be zero ddx=x=0si5_e.

Stationary. A stationary flow is a flow that does not change over time. As an example, consider the flow of a water film along an inclined plane. The established flow profile is driven by gravity. As gravity is constant during the experiment, the system will eventually reach a state in which the liquid is flowing down the plane at a constant velocity. It will not change over time. In this case, all (partial) derivatives with respect to time are zero: ddt=t=0si6_e

Incompressible. As discussed in section 9.2, section 9.3.3, and section 9.3.4, liquids and gases have significantly different densities. A very important result from this is the fact that liquids are generally considered to be incompressible, whereas gases are compressible. This is a pretty good approximation of the fluid physics. Therefore, if discussing liquid flows, one can often assume the fluid to be incompressible. In these cases, all (partial) derivatives of the density are zero: dρd=ρ=0si7_e

10.1.3 Frames of Reference

Throughout the derivation of the fundamental equations of fluid mechanics we will work in the infinitesimal control volume dV = dx dy dz. There are two principal ways of placing these frames of reference in the flow field (see Fig. 10.1).

f10-01-9781455731411
Fig. 10.1 a) Eulerian frame of reference: fixed in space. b) Lagrangian frames of reference: moving in space.

Eulerian Frame of Reference. The Eulerian frame of reference uses control volumes which are fixed in space (see Fig. 10.1a). For the mathematical derivation of the fundamental equations of fluid mechanics, this is the simpler approach. In the Eulerian frame of reference, the direction of the flow inside of the control volume can change direction frequently. The control volume is located at the position s0 which is identical to s1, etc. The position vector therefore is not a function of time.

In the macroscopic world, there is a simple analogy. If standing next to a river, keep staring at one specific location only. When standing on a bridge, one can concentrate on only one position. The flow field will, for the most part, not change and the view may become slightly fatiguing. Eventually, a small turbulence may be introduced into the flow field, e.g., if someone is throwing a stone into the water. Then the image changes as waves cross into the observer’s field of view. This is a good approximation of an Eulerian frame of reference.

Lagrangian Frame of Reference. The Lagrangian frame of reference moves with the flow field (see Fig. 10.1b). Returning to our analogy, the Lagrangian frame of reference of the flow field is the view the observer will have if sitting in a boat that is taken along by the moving water. The position of the observer, i.e., the frame of reference constantly changes. At the beginning for t0, the control volume will be at the position s0 at which the flow field has the velocity v0. At t1, the control volume will have traveled to s1 at which the flow field now has the velocity v1. The position vector s is changing and therefore a function of time. Connecting the individual position vectors over time allows reconstructing the particle trajectory.

Particle Trajectories. The particle trajectories are a good way of visualizing the flow field. For this, small particles are seeded into the flow field at different initial positions sA,0, sB,0, sC,0, etc. The flow field is then monitored with a camera with very long exposure settings. The image recorded will highlight the movement of the particles in the form of the trajectory lines as the particles move through the flow field. Of course, the same is also true for control volumes being moved in the flow field.

Conversion From Eulerian to Lagrangian Frames of Reference. Eulerian and Lagrangian frames of reference can be converted into each other, respectively. In the Eulerian frame of reference, position, velocity, and acceleration are only a function of time as the position vector is not changing

position:s(t)velocity:v(t)=ds(t)dtacceleration:a(t)=dv(t)dt=d2s(t)dt2

si8_e

In the Lagrangian frame of reference, the position vector is a function of the previous position, i.e., s1 depends on s0, etc. This can be expressed as

s(s(t),t)

si9_e

Therefore, for finding velocity and acceleration, a total differential (see section 3.1.5.2) has to be used. We find

position:s(s(t),t)

si10_e  (Eq. 10.1)

velocity:v(s(t),t)=ds(s(t),t)dt=st+sxdxt+sydydt+szdzdt=st+(v)s

si11_e  (Eq. 10.2)

acceleration:a(v(t),t)=dv(v(t),t)dt=vt+vxdxt+vydydt+szdzdtvt+(v·)v

si12_e  (Eq. 10.3)

Just as a short note, the expansion of the total integral of v (s (t) , t) is carried out as

dv(s(t),t)=stdt+sxdx+sydy+szdz

si13_e

Dividing this equation by dt results in the above total integral.

Interpretation. As we can see, the Lagrangian terms for speed and velocity contain the derivatives with respect to time, which we also find in the Eulerian frame of reference. However, they are partial derivatives stsi14_e and vtsi15_e and they are complemented, each with a second term. These terms can be thought of as the convective contribution of the moving fluid field. If the fluid was not moving, the terms (v)ssi16_e and (v)vsi17_e would disappear. The higher the flow speed of the fluid, the less important the terms originating from the particle movement will become. As a reference: If a boat was floating in a river, the first term in the velocity equation would be the Eulerian speed contribution, i.e., the movement of the boat due to the action of its engine. The second term would be due to the movement of the water in the river. If the engine is stopped, only this term will determine the overall movement speed of the boat.

10.1.4 Flow Tubes

Fig. 10.2b already highlights an important detail of some flow fields. As stated, the particle trajectories can be seen by observing the movement path of the particle as it is transported through the flow field. In flow fields with no turbulence, these distinct trajectories do not cross. If this is the case, the flow field can be effectively split up into individual flow tubes, each of which is a one-dimensional control volume with the position variable s being the only variable of the control volume. s will be winding along the axis of the flow tube that can be curved in space. However, the fluid mechanical problem to solve is effectively one-dimensional. One important implication of this is that there is no mass transport into or out of the flow tube. In numerous cases, using flow tubes instead of describing the two-dimensional fluid mechanical problem is significantly easier. We will discuss the fluid mechanics in these control volumes in section 14.

f10-02-9781455731411
Fig. 10.2 a) Trajectories of particles, seeded at different initial positions sA,0,sB,0,sC,0si1_e etc. as they are transported in the flow field. b) This behavior gives rise to the concept of flow tubes that can be used as one-dimensional control volumes in laminar flow fields.

10.2 Continuity Equation

10.2.1 Mass and Volume Flow

We will briefly introduce two equations that we will need for converting mass to volume flows. Mass flow is referred to as the change of mass over time dmdt=m˙si18_e. In analogy, the volume flow is the change of volume over time dVdt=V˙si19_e. Of course, mass and volume are related by the density

ρ=mV

si20_e  (Eq. 10.4)

=m˙V˙

si21_e  (Eq. 10.5)

10.2.2 Geometry, Flow Velocity, and Flowrate

In a fluidic system, the volume can be calculated as

V=Al

si22_e  (Eq. 10.6)

if the cross-section A and the length of the wetted segment l of the system are known. From Eq. 10.6 the flow velocity can be calculated as the derivative with respect to time

dVdt=Adldt=Av=V˙

si23_e  (Eq. 10.7)

from which we can derive

m˙=ρV˙=ρAv

si24_e  (Eq. 10.8)

We will use Eq. 10.7 often for calculating mass flows in and out of control volumes.

10.2.3 Flow rate

In many textbooks on microfluidics, you will find the symbol Q instead of dVdtsi25_e. It is used to refer to the flowrate, i.e., the volume of fluid flow through a system in a given amount of time.

10.3 Integral Representation of the Flowrate

Until now, we have studied cases in which the velocity v is constant over a given cross-section area A. However, as we will see, in most microfluidic systems, the velocity profile is a function of the position on the cross-section and therefore the flowrate must be calculated as the integral over the cross-section according to

dVdt=AvdA

si26_e  (Eq. 10.9)

m˙incompr.=ρAvdA

si27_e  (Eq. 10.10)

m˙compr.=ρAvdA

si28_e  (Eq. 10.11)

where we have also given the expression of the mass flow, which is simply the flowrate multiplied by the density if the density is constant (see Eq. 10.10). If this is not the case, the density is moved into the integral, resulting in Eq. 10.11.

10.4 Mass Balance

We will now derive the continuity equation that is one of the most important equations in fluid mechanics. It describes the conservation of mass, which we will study at the infinitesimal control volume dV = dx dy dz (see Fig. 10.3). Mass is flowing into and out of the control volume along the axes via the surface planes of the control volume. Per definition, mass flowing into the control volume has positive sign, whereas mass flowing out of the control volume has negative sign. Using Eq. 10.8, the sum of all in- and outflowing mass will change the total mass of the control volume over time

f10-03-9781455731411
Fig. 10.3 Schematic of the in- and outflux of mass in an infinitesimal control volume.

mt=(ρdV)t=(ρdxdydz)t=ρtdxdydz

si29_e  (Eq. 10.12)

The balance of mass flowing in and out of the control volume is

inx:(ρvx(ρvx+(ρvx)x)dx)dydz=(ρvx)xdxdydz

si30_e  (Eq. 10.13)

iny:(ρvy(ρvy+(ρvy)y)dy)dxdz=(ρvy)ydxdydz

si31_e  (Eq. 10.14)

inz:(ρvz(ρvz+(ρvz)z)dz)dxdy=(ρvz)zdxdydz

si32_e  (Eq. 10.15)

Here we have used the partial derivatives along the respective axes which account for changes in the respective variables over the lengths dx, dy, or dz of the control volume (see section 3.1.5.2). The sum of these mass flows must be equal to the change of mass in the control volume over time. So using Eq. 10.12 we derive

ptdxdydz=ρvxxdxdydzρvyydxdydzρvzzdxdydzρt=ρvxxρvyyρvzzρt+(ρv)=0

si33_e  (Eq. 10.16)

Eq. 10.16 is referred to as the continuity equation for compressible fluids. If we consider incompressible fluids, ρ is a constant and therefore ρt=0si34_e and (ρv)=vsi35_e. The equation then simplifies to

v=0

si36_e  (Eq. 10.17)

10.5 Derivation using Gauss’s Theorem

The continuity equation can also be derived using Gauss’s theorem (see section 7.2.1). In this case, we use ρv)si50_e (which is the directional mass flow in vector form) as the field in Eq. 7.10. We then have to consider the total mass transported across the boundaries of our control volume. Considering Fig. 7.2a as an arbitrary control volume, the mass flow across the boundary can be calculated as

(ρv)dV

si37_e  (Eq. 10.18)

The integral is assumed to be negative, as the flow is directed out of the control volume. The change in mass over time in the control volume is given by Eq. 10.12. Using Eq. 10.18 this results in

(ρv)dV=ρtdVρt+(ρv)=0

si38_e

which is identical to Eq. 10.16. In simple terms, the continuity equation describes that if mass is transported into the system via its boundaries, the density of the system must increase in order to store this mass. On the other hand, if mass is transported out of the system via its boundaries, the density of the system must decrease. If the liquid is incompressible, the density of the system cannot change. Therefore, if liquid is flowing into the system across one side of the boundary, the same amount of mass must be flowing out of the system across a different side of the boundary. If the total boundary surface of the system is constant, Eq. 10.8 allows simplifying this statement to “The sum of all influx and outflux into and out of the system must sum up to zero.”

10.6 Combined Convection and Diffusion: The Convection-Diffusion Equation

In many practical fluid mechanical problems, diffusion is often neglected due to the slow mass transport over long length scales. Mass transport is therefore solely based on the in- and outflow of mass into and out of the control volumes. This fluid movement is generally referred to as convection. However, in many microfluidic applications, the length scales are reasonable small and it may be necessary to account for both diffusion and convection.

We have already derived Fick’s second law of diffusion (see Eq. 6.90), which is the most general PDE for describing diffusion. It is very straightforward to add this equation to the continuity equation we just derived. The most general case of the continuity equation is Eq. 10.16. It is a homogeneous PDE because we assume no change of mass within the control volume (thus the conservation of mass). However, if we add diffusion, we expect a change of mass over time due to diffusive transport into and out of the control volume. Thus we add a right-hand side to Eq. 10.16 resulting in

ρt+(ρv)=t(masspervolumetransportedinto/outofthecontrolvolumebydiffusion)

si39_e

If we look at Eq. 6.90 we can see that this is equivalent to the left-hand side of the equation. Thus we find

ρt+(ρv)=ρt|bydiffusion=DΔρ

si40_e  (Eq. 10.19)

We can write this equation in this way because the mass concentration ρ effectively is nothing else then the density ρ, i.e., mass m per volume V . However, diffusion involves always more than one species. As an example, consider the movements of ions in a buffer which contains two types of ions (cations and anions). Here we have to apply Eq. 10.19 for both species, which is why we usually write it as a function of the mass concentration ρi of the ith species in the mixture. This results in

ρit+(ρiv)=DΔρi

si41_e  (Eq. 10.20)

which is the most commonly found version of this equation. It balances the change in concentration of the ith species in the mixture over time (first term of the left-hand side) as a consequence of diffusion (right-hand side) and convection (second term on the left-hand side). Eq. 10.20 is generally referred to as the convection-diffusion equation as it combined the two mass transport mechanisms.

10.7 Summary

In this section we have derived the continuity equation which is the first of the conservative equations that form the basis of fluid mechanics. As we have seen, the continuity equation is derived from the fact that mass must be conserved throughout a flow field. We have derived two versions of the equation: for compressible and for incompressible flows. As we have seen, we need a total of six equations to solve for all required field variables. One of these equations is the continuity equation. In the following section we will derive the Navier-Stokes equation which will give rise to three additional equations.

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset