Fluids Review Design

This section is meant to be a refresher on fluid mechanics. It will only cover the topics of fluids mechanics that will be used heavily in the course.

If you wish to review fluid mechanics in (much) more detail, please refer to this guide. If you wish to review from a legitimate textbook, you can find a pdf of good book by Frank White here.

Important Terms and Equations

Terms:

  1. Head
  2. Streamline
  3. Head loss
  4. Laminar
  5. Turbulent
  6. Moody Diagram
  7. Viscosity
  8. Driving head
  9. Vena Contracta/Coefficient of Contraction

Equations:

  1. Bernoulli equation
  2. Energy equation
  3. Darcy-Weisbach equation
  4. Reynolds number
  5. Swamee-Jain equation
  6. Hagen-Poiseuille equation
  7. Orifice equation

Introductory Concepts

Before diving in to the rest of the fluids review document, there are a few important concepts which will be the foundation for building your understanding of fluid mechanics. One must walk before they can run, and similarly, the basics of fluid mechanics must be understood before moving on to the more fun sections of this document.

Continuity Equation

Continuity is simply an application of mass balance to fluid mechanics. It states that the cross sectional area \(A\) that a fluid flows through multiplied by the fluid’s average flow velocity \(\bar v\) must equal the fluid’s flow rate \(Q\):

(4)\[Q = \bar v A\]

Note

The line above the \(v\) is called a ‘bar,’ and represents an average. Any variable can have a bar. In this case, we are adding the bar to velocity \(v\), turning it into average velocity \(\bar v\). This variable is pronounced ‘v bar.’

In CEE 4540, we deal primarily with flow through pipes. For a circular pipe, \(A = \pi r^2\). Substituting diameter in for radius, \(r = \frac{D}{2}\), we get \(A = \frac{\pi D^2}{4}\). You will often see this form of the continuity equation being used to relate the flow rate in a pipe to the fluid velocity and pipe diameter:

(5)\[Q = \bar v \frac{\pi D^2}{4}\]

The continuity equation is also useful when flow is going from one geometry to another. In this case, the flow in one geometry must be the same as the flow in the other, \(Q_1 = Q_2\), which yields the following equations:

(6)\[\bar v_1 A_1 = \bar v_2 A_2\]
(7)\[\bar v_1 \frac{\pi D_1^2}{4} = \bar v_2 \frac{\pi D_2^2}{4}\]
Such that:
\(Q =\) fluid flow rate
\(\bar v =\) fluid average velocity
\(A =\) pipe area
\(r =\) pipe radius
\(D =\) pipe diameter

An example of changing flow geometries is when the a change in pipe size occurs in a circular piping system, as is demonstrated below. The flow through \({\rm pipe} \, 1\) must be the same as the flow through \({\rm pipe} \, 2\).

internal figure

Fig. 3 Flow going from a small diameter pipe to a large one. The flow through each pipe must be the same.

Laminar and Turbulent Flow

Considering that this class deals with the flow of water through a water treatment plant, understanding the characteristics of the flow is very important. Thus, it is necessary to understand the most common characteristic of fluid flow: whether it is laminar or turbulent. Laminar flow is very smooth and highly ordered. Turbulent flow is chaotic, messy, and disordered. The best way to understand each flow and what it looks like is visually, like in this video or the wikipedia image below. Please ignore the part of the video after the image of the tap.

Laminar flow, turbulent flow, and the transition

Fig. 4 This is a beautiful example of the difference between ordered, smooth laminar and chaotic turbulent flow.

A numeric way to determine whether flow is laminar or turbulent is by finding the Reynolds number, \({\rm Re}\). The Reynolds number is a dimensionless parameter that compares inertia, represented by the average flow velocity \(\bar v\) times a length scale \(D\) to viscosity, represented by the kinematic viscosity \(\nu\). Click here for a brief video explanation of viscosity. If the Reynolds number is less than 2,100 the flow is considered laminar. If it is more than a certain value, it is considered turbulent.

(8)\[{\rm Re = \frac{inertia}{viscosity}} = \frac{\bar vD}{\nu}\]

There is a transition between laminar and turbulent flow which is not yet well understood. To simplify this phenomenon and make it possible to code for laminar or turbulent flow, we assume that the transition occurs at \(\rm{Re} = 2100\). The flow regime is assumed to be laminar below this value and turbulent above it. This variable is coded into aide_design as pc.RE_TRANSITION_PIPE. We will neglect transitional flow.

Fluid can flow through very many different geometries like a pipe, a rectangular channel, or any other shape. To account for this, the characteristic length scale is quantified as the hydraulic diameter, which can be applied to any geometry. For circular pipes, which are the most common geometry you’ll encounter in this class, the hydraulic diameter is simply the pipe diameter.

Here are other commonly used forms of the Reynolds number equation. They are the same as the one above, just with the substitutions \(Q = \bar v \frac{\pi D^2}{4}\) and \(\nu = \frac{\mu}{\rho}\)

(9)\[{\rm{Re}} = \frac{\bar vD}{\nu} = \frac{4Q}{\pi D\nu} = \frac{\rho \bar vD}{\mu}\]
Such that:
\(Q\) = fluid flow rate in pipe
\(D\) = pipe diameter
\(\bar v\) = fluid velocity
\(\nu\) = fluid kinematic viscosity
\(\mu\) = fluid dynamic viscosity

See also

Function in aide_design: pc.re_pipe(FlowRate, Diam, Nu) Returns the Reynolds number in a circular pipe. Functions for finding the Reynolds number through other conduits and geometries can also be found in physchem.py within aide_design.

Note

Laminar and turbulent flow are described as two different flow regimes. When there is a characteristic of flow and different categories of the characteristic, each category is referred to as a flow regime. For example, the Reynolds number describes a flow characteristic, and its categories, referred to as flow regimes, are laminar or turbulent.

Streamlines and Control Volumes

Both streamlines and control volumes are tools to compare different sections of a system. For this class, this system will always be hydraulic.

Imagine water flowing through a pipe. A streamline is the path that a particle would take if it could be placed in the fluid without changing the original flow of the fluid. A more technical definition is “a line which is everywhere parallel to the local velocity vector.” Computational tools, dyes (in water), or smoke (in air) can be used to visualize streamlines.

A control volume is just an imaginary 3-dimensional shape in space. Its boundaries can be placed anywhere by the person applying the control volume, and once set the boundaries remain fixed in space over time. These boundaries are usually chosen to compare two relevant surfaces to each other. The entirety of a control volume is usually not shown, as it is often unnecessary. This is shown in the following image:

Control volume simplification

Fig. 5 While the image on the left indicates a complete control volume, control volumes are usually shortened to only include the relevant surfaces, in which the control volume intersects the fluid. This is shown in the image on the right.

Important

Many images will be used over the course of this class to show hydraulic systems. A standardized system of lines will be used throughout them all to distinguish reference elevations from control volumes from streamlines. This system is described in the image below.

Image control volumes

Fig. 6 A convention for figure control volume and streamlines will be very helpful throughout this course.

The Bernoulli and Energy Equations

As explained in CEE 3310 with more details than most of you wanted to know, the Bernoulli and energy equations are incredibly useful in understanding the transfer of the fluid’s energy throughout a streamline or through a control volume. The Bernoulli equation applies to two different points along one streamline, whereas the energy equation applies across a control volume. The energy of a fluid has three forms: pressure, potential (deriving from elevation), and kinetic (deriving from velocity).

The Bernoulli Equation

These three forms of energy expressed above make up the Bernoulli equation:

(10)\[\frac{p_1}{\rho g} + {z_1} + \frac{v_1^2}{2g} = \frac{p_2}{\rho g} + {z_2} + \frac{v_2^2}{2g}\]
Such that:
\(p\) = pressure
\(\rho\) = fluid density
\(g\) = acceleration due to gravity, in aide_design as pc.gravity
\(z\) = elevation relative to a reference
\(v\) = fluid velocity

Notice that each term in this form of the Bernoulli equation has units of \([L]\), even though the terms represent the energy of water, which has units of \(\frac{[M] \cdot [L]^2}{[T]^2}\). When energy of water is described in units of length, the term used is called head.

There are two important distinctions to keep in mind when using head to talk about energy. First is that head is dependent on the density of the fluid under consideration. Take mercury, for example, which is around 13.6 times more dense than water. 1 meter of mercury head is therefore equivalent to around 13.6 meters of water head. Second is that head is independent of the amount of fluid being considered, as long as all the fluid is the same density. Thus, raising 1 liter of water up by one meter and raising 100 liters of water up by one meter are both equivalent to giving the water 1 meter of water head, even though it requires 100 times more energy to raise the hundred liters than to raise the single liter. Since we are concerned mainly with water in this class, we will refer to ‘water head’ simply as ‘head’.

Going back to the Bernoulli equation, the \(\frac{p}{\rho g}\) term is called the pressure head, \(z\) the elevation head, and \(\frac{v^2}{2g}\) the velocity head. The following diagram shows these various forms of head via a 1 meter deep bucket (left) and a jet of water shooting out of the ground (right).

Different forms of head

Fig. 7 The three forms of hydraulic head.

Assumption in using the Bernoulli equation

Though there are many assumptions needed to confirm that the Bernoulli equation can be used, the main one for the purpose of this class is that energy is not gained or lost throughout the streamline being considered. If we consider more precise fluid mechanics terminology, then “friction by viscous forces must be negligible.” What this means is that the fluid along the streamline being considered is not losing energy to viscosity. Energy can only be transferred between its three forms if this equation is to be used, it can’t be gained or lost.

Example problems

Here is a simple worksheet with very straightforward example problems using the Bernoulli equation. Note that the solutions use the pressure-form of the Bernoulli equation. This just means that every term in the equation is multiplied by \(\rho g\), so the pressure term is just \(P\). The form of the equation does not affect the solution to the problem it helps solved.

The Energy Equation

The assumption necessary to use the Bernoulli equation, which is stated above, represents the key difference between the Bernoulli equation and the energy equation for the purpose of this class. The energy equation accounts for the (L)oss of energy from both the fluid flowing, \(h_L\), and any other energy drain, like the charging of a (T)urbine, \(h_T\). It also accounts for any energy inputs into the system, \(h_P\), which is usually caused by a (P)ump within the control volume.

(11)\[\frac{p_{1}}{\rho g} + z_{1} + \alpha_{1} \frac{\bar v_{1}^2}{2g} + h_P = \frac{p_{2}}{\rho g} + z_{2} + {\alpha_{2}} \frac{\bar v_{2}^2}{2g} + h_T + h_L\]

You’ll also notice the \(\alpha\) term attached to the velocity head. This is a correction factor for kinetic energy, and will be neglected in this class. In the Bernoulli equation, the velocity of the streamline of water is considered, \(v\). The energy equation, however compares control surfaces instead of streamlines, and the velocities across a control surface many not all be the same. Hence, \(\bar v\) is used to represent the average velocity. Since AguaClara does not use pumps nor turbines, \(h_P = h_T = 0\). With these simplifications, the energy equation can be written as follows:

(12)\[\frac{p_{1}}{\rho g} + z_{1} + \frac{\bar v_{1}^2}{2g} = \frac{p_{2}}{\rho g} + z_{2} + \frac{\bar v_{2}^2}{2g} + h_L\]

This is the form of the energy equation that you will see over and over again in CEE 4540. To summarize, the main difference between the Bernoulli equation and the energy equation for the purposes of this class is energy loss. The energy equation accounts for the fluid’s loss of energy over time while the Bernoulli equation does not. So how can the fluid lose energy?

Headloss

Head(L)oss, \(h_L\) is a term that is ubiquitous in both this class and fluid mechanics in general. Its definition is exactly as it sounds: it refers to the loss of energy of a fluid as it flows through space. There are two components to head loss: major losses caused by pipe-fluid (f)riction, \(h_{\rm{f}}\), and minor losses caused by fluid-fluid friction resulting from flow (e)xpansions, \(h_e\), such that \(h_L = h_{\rm{f}} + h_e\).

Major Losses

These losses are the result of friction between the fluid and the surface over which the fluid is flowing. A force acting parallel to a surface is referred to as shear. It can therefore be said that major losses are the result of shear between the fluid and the surface it’s flowing over. To help in understanding major losses, consider the following example: imagine, as you have so often in physics class, pushing a large box across the ground. Friction is what resists your efforts to push the box. The farther you push the box, the more energy you expend pushing against friction. The same is true for water moving through a pipe, where water is analogous to the box you want to move, the pipe is similar to the floor that provides the friction, and the major losses of the water through the pipe is analogous to the energy you expend by pushing the box.

In this class, we will be dealing primarily with major losses in circular pipes, as opposed to channels or pipes with other geometries. Fortunately for us, Henry Darcy and Julius Weisbach came up with a handy equation to determine the major losses in a circular pipe under both laminar and turbulent flow conditions. Their equation is logically but unoriginally named the Darcy-Weisbach equation and is shown below:

(13)\[h_{\rm{f}} \, = \, {\rm{f}} \frac{L}{D} \frac{\bar v^2}{2g}\]

Substituting the continuity equation \(Q = \bar vA\) in the form of \(\bar v^2 = \frac{16Q^2}{\pi^2 D^4}\) gives another, equivalent form of Darcy-Weisbach which uses flow, \(Q\), instead of velocity, \(\bar v\):

(14)\[h_{\rm{f}} \, = \,{\rm{f}} \frac{8}{g \pi^2} \frac{LQ^2}{D^5}\]
Such that:
\(h_{\rm{f}}\) = major loss, \([L]\)
\(\rm{f}\) = Darcy friction factor, dimensionless
\(L\) = pipe length, \([L]\)
\(Q\) = pipe flow rate, \(\frac{[L]^3}{[T]}\)
\(D\) = pipe diameter, \([L]\)

See also

Function in aide_design: pc.headloss_fric(FlowRate, Diam, Length, Nu, PipeRough) Returns only major losses. Works for both laminar and turbulent flow.

Darcy-Weisbach is wonderful because it applies to both laminar and turbulent flow regimes and contains relatively easy to measure variables. The one exception is the Darcy friction factor, \(\rm{f}\). This parameter is an approximation for the magnitude of friction between the pipe walls and the fluid, and its value changes depending on the whether or not the flow is laminar or turbulent, and varies with the Reynolds number in both flow regimes.

For laminar flow, the friction factor can be determined from the following equation:

(15)\[{\rm{f}} = \frac{64}{\rm{Re}}\]

For turbulent flow, the friction factor is more difficult to determine. In this class, we will use the Swamee-Jain equation:

(16)\[{\rm{f}} = \frac{0.25} {\left[ \log \left( \frac{\epsilon }{3.7D} + \frac{5.74}{{\rm Re}^{0.9}} \right) \right]^2}\]
Such that:
\(\epsilon\) = pipe roughness, \([L]\)
\(D\) = pipe diameter, \([L]\)

See also

Function in aide_design: pc.fric(FlowRate, Diam, Nu, PipeRough) Returns \(\rm{f}\) for laminar or turbulent flow. For laminar flow, use ‘0’ for the PipeRough input.

The simplicity of the equation for \(\rm{f}\) during laminar flow allows for substitutions to create a very useful, simplified equation for major losses during laminar flow. This simplification combines the Darcy-Weisbach equation, the equation for the Darcy friction factor during laminar flow, and the Reynold’s number formula:

(17)\[h_{\rm{f}} \, = \,{\rm{f}} \frac{8}{g \pi^2} \frac{LQ^2}{D^5}\]
(18)\[{\rm{f}} = \frac{64}{\rm{Re}}\]
(19)\[{\rm{Re}}=\frac{4Q}{\pi D\nu}\]

To form the Hagen-Poiseuille equation for major losses during laminar flow, and only during laminar flow:

(20)\[h_{\rm{f}} = \frac{128\mu L Q}{\rho g\pi D^4}\]
(21)\[h_{\rm{f}} = \frac{32\nu L\bar v}{ g D^2}\]

The significance of this equation lies in its relationship between \(h_{\rm{f}}\) and \(Q\). Hagen-Poiseuille shows that the terms are directly proportional (\(h_{\rm{f}} \propto Q\)) during laminar flow, while Darcy-Weisbach shows that \(h_{\rm{f}}\) grows with the square of \(Q\) during turbulent flow (\(h_{\rm{f}} \propto Q^2\)). As you will soon see, minor losses, \(h_e\), will grow with the square of \(Q\) in both laminar and turbulent flow. This has implications that will be discussed later, in the flow control section.

In 1944, Lewis Ferry Moody plotted a ridiculous amount of experimental data, gathered by many people, on the Darcy-Weisbach friction factor to create what we now call the Moody diagram. This diagram has the friction factor \(\rm{f}\) on the left-hand y-axis, relative pipe roughness \(\frac{\epsilon}{D}\) on the right-hand y-axis, and Reynolds number \(\rm{Re}\) on the x-axis. The Moody diagram is an alternative to computational methods for finding \(\rm{f}\).

Moody diagram

Fig. 8 This is the famous and famously useful Moody diagram.

Minor Losses

Unfortunately, there is no simple ‘pushing a box across the ground’ example to explain minor losses. So instead, consider a hydraulic jump. In the video, you can see lots of turbulence and eddies in the transition region between the fast, shallow flow and the slow, deep flow. The high amount of mixing of the water in the transition region of the hydraulic jump results in significant friction between water and water (recall that the measure of a fluid’s resistance to internal, fluid-fluid friction is called viscosity). This turbulent, eddy-induced, fluid-fluid friction results in minor losses, much like fluid-pipe friction results in major losses.

As is the case in a hydraulic jump, a flow expansion (from shallow flow to deep flow) creates the turbulent eddies that result in minor losses. This will be a recurring theme in throughout the course: minor losses are caused by flow expansions. Imagine a pipe fitting that connects a small diameter pipe to a large diameter one, as shown in the image below. The flow must expand to fill up the entire large diameter pipe. This expansion creates turbulent eddies near the union between the small and large pipes, and these eddies cause minor losses. You may already know the equation for minor losses, but understanding where it comes from is very important for effective AguaClara plant design. For this reason, you are strongly recommended to read through the full derivation, in Fluids Review Derivations.

There are three forms of the minor loss equation that you will see in this class:

(22)\[{\rm{ \mathbf{First \, form:} }} \,\,\, h_e = \frac{\left( \bar v_{in} - \bar v_{out} \right)^2}{2g}\]
(23)\[{\rm{ \mathbf{Second \, form:} }} \,\,\, h_e = \frac{\bar v_{in}^2}{2g}{\left( {1 - \frac{A_{in}}{A_{out}}} \right)^2} = \,\,\, \frac{\bar v_{in}^2}{2g} \mathbf{K_e^{'}}\]
(24)\[{\rm{ \mathbf{Third \, form:} }} \,\,\, h_e = \frac{\bar v_{out}^2}{2g}{\left( {\frac{A_{out}}{A_{in}}} -1 \right)^2} = \,\,\,\, \frac{\bar v_{out}^2}{2g} \mathbf{K_e}\]
Such that:
\(K_e^{'}, \,\, K_e\) = minor loss coefficients, dimensionless
Function in aide_design:
pc.headloss_exp_general(Vel, KMinor) Returns \(h_e\). Can be either the second or third form due to user input of both velocity and minor loss coefficient. It is up to the user to use consistent \(\bar v\) and \(K_e\).
pc.headloss_exp(FlowRate, Diam, KMinor) Returns \(h_e\). Uses third form, \(K_e\).

Note: You will often see \(K_e^{'}\) and \(K_e\) used without the \(e\) subscript, they will appear as \(K^{'}\) and \(K\).

The \(in\) and \(out\) subscripts in each of the three forms refer to the diagram that was used for the derivation:

The second and third forms are the ones which you are probably most familiar with. The distinction between them, however, is critical. First, consider the magnitudes of \(A_{in}\) and \(A_{out}\). \(A_{in}\) can never be larger than \(A_{out}\), because the flow is expanding. When flow expands, the cross-sectional area it flows through must increase. As a result, both \(\frac{A_{out}}{A_{in}} > 1\) and \(\frac{A_{in}}{A_{out}} < 1\) must always be true. This means that \(K^{'}\) can never be greater than 1, while \(K\) technically has no upper limit.

If you have taken CEE 3310, you have seen tables of minor loss coefficients like this one, and they almost all have coefficients greater than 1. This implies that these tables use the third form of the minor loss equation as we have defined it, where the velocity is \(\bar v_{out}\). There is a good reason for using the third form over the second one: \(\bar v_{out}\) is far easier to determine than \(\bar v_{in}\). Consider flow through a pipe elbow, as shown in the image below.

In order to find \(\bar v_{out}\), we first need to know which point is \(out\) and which point is \(in\). A simple way to distinguish the two points is that \(in\) occurs when the flow is most contracted, and \(out\) occurs when the flow has fully expanded after that maximal contraction. Going on these guidelines, point ‘B’ above would be \(in\), since it represents the most contracted flow in the elbow-pipe system. Therefore point ‘C’ would be \(out\), as it is the point where the flow has fully expanded after its compression in ‘B.’

\(\bar v_{out}\) is easy to determine because it is the velocity of the fluid as it flows through the entire area of the pipe. Thus, \(\bar v_{out}\) can be found with the continuity equation, since the flow through the pipe and its diameter are easy to measure, \(\bar v_{out} = \frac{4 Q}{\pi D^2}\). On the other hand, \(\bar v_{in}\) is difficult to find, as the area of the contracted flow is dependent on the exact geometry of the elbow. This is why the third form of the minor loss equation, as we have defined it, is the most common.

Head Loss = Elevation Difference Trick

This trick, also called the ‘control volume trick,’ or more colloquially, the ‘head loss trick,’ is incredibly useful for simplifying hydraulic systems and is used all the time in this class.

Consider the following image, which was taken from the Flow Control and Measurement powerpoint.

In systems like this, where an elevation difference is causing the flow of water, the elevation difference is called the driving head. In the system above, the driving head is the elevation difference between the water level and the end of the tubing. Usually driving head is written as \(\Delta z\) or \(\Delta h\), though above it is labelled as \(h_L\).

This image is violating the energy equation by saying that the elevation difference between the water in the tank and the end of the tube is \(h_L\). It implies that all of the driving head, \(\Delta z\), is lost to head loss and therefore that no water is flowing out of the tubing, which is not true. Let’s apply the energy equation between the two red points. Pressures are atmospheric at both points and the velocity of water at the top of tank is negligible.

(25)\[\rlap{\Bigg/}\frac{p_{1}}{\rho g} + z_{1} + \rlap{\Bigg/}\frac{\bar v_{1}^2}{2g} = \rlap{\Bigg/}\frac{p_{2}}{\rho g} + z_{2} + \frac{\bar v_{2}^2}{2g} + h_L\]

We now get:

(26)\[\Delta z = \frac{\bar v_2^2}{2g} + h_L\]

This contradicts the image above, which says that \(\Delta z = h_L\) and neglects \(\frac{\bar v_2^2}{2g}\). The image above is correct, however, if you apply the head loss trick. The trick incorporates the \(\frac{\bar v_2^2}{2g}\) term into the \(h_L\) term as a minor loss. See the math below:

(27)\[\Delta z = \frac{\bar v_2^2}{2g} + h_e + h_f\]
(28)\[\Delta z = \frac{\bar v_2^2}{2g} + \left( \sum K \right) \frac{\bar v_2^2}{2g} + h_f\]
(29)\[\Delta z = \left( 1 + \sum K \right) \frac{\bar v_2^2}{2g} + h_f\]

This last step incorporated the kinetic energy term of the energy equation, \(\frac{\bar v_2^2}{2g}\), into the minor loss equation by saying that its \(K\) is 1. From here, we reverse our steps to get \(\Delta z = h_L\)

(30)\[\Delta z = h_e + h_f\]
(31)\[\Delta z = h_L\]

By applying the head loss trick, you are considering the entire flow of water out of a control volume as lost energy. This is just an algebraic trick, the only thing to remember when applying this trick is that \(\sum K\) will always be at least 1, even if there are no ‘real’ minor losses in the system.

The Orifice Equation

This equation is one that you’ll see again and again throughout this class. Understanding it now will be invaluable, as future concepts will use and build on this equation.

Vena Contracta

Before describing the equation, we must first understand the concept of a vena contracta. Refer once more to this image of flow through a pipe elbow.

The flow contracts as the fluid moves from point ‘A’ to point ‘B.’ This happens because the fluid can’t make a sharp turn at the corner of the elbow. Instead, the streamline closest to the sharp turn makes a slow, gradual change in direction, as shown in the image. As a result of this gradual turn, the cross-sectional area the fluid is flowing through at point ‘B’ is less than the cross-sectional area it flows through at points ‘A’ and ‘C’. Written as an equation, \(A_{csB} < A_{csA} = A_{csC}\), where the \(_{csA}\) stands for ‘control surface \(A\)’ subscript

The term ‘vena contracta’ describes the phenomenon of contracting flow due to streamlines being unable to make sharp turns. \(\Pi_{vc}\) is a ratio between the flow area at the vena contracta, \(A_{csB}\), which is when the flow is maximally contracted, and the flow area before the contraction, \(A_{csA}\). In the image above, the equation for the vena contracta coefficient would be:

(32)\[\Pi_{vc} = \frac{A_{csB}}{A_{csA}}\]

Note that what this class calls \(\Pi_{vc}\) is often referred to as a ‘Coefficient of Contraction,’ \(C_c\), in other engineering courses and settings. When the most extreme turn a streamline must make is 90°, the value of the vena contracta coefficient is close to 0.62. This parameter is in aide_design as pc.RATIO_VC_ORIFICE. The vena contracta coefficient value is a function of the flow geometry.

A vena contracta coefficient is not a minor loss coefficient. Though the equations for the two both involve contracted and non-contracted areas, these coefficients are not the same. Refer to the flow through a pipe elbow image above. The minor loss coefficient equation uses the areas of points ‘B’ and ‘C,’ while the vena contracta coefficient uses the areas of points ‘A’ and ‘B.’ Additionally, the equations to calculate the coefficients themselves are not the same. Confusing the two coefficients is common mistake that this paragraph will hopefully help you to avoid.

Origin

The orifice equation is derived from the Bernoulli equation as applied to the red points in the following image:

At point A, the pressure is atmospheric and the instantaneous velocity is negligible as the water level in the bucket drops slowly. At point B, the pressure is also atmospheric. We define the difference in elevations between the two points, \(z_A - z_B\), to be \(\Delta h\). With these simplifications (\(p_A = \bar v_A = p_B = 0\)) and assumptions (\(z_A - z_B = \Delta h\)), the Bernoulli equation becomes:

(33)\[\Delta h = \frac{\bar v_B^2}{2g}\]

Substituting the continuity equation \(Q = \bar v A\) in the form of \(\bar v_B^2 = \frac{Q^2}{A_{vc}^2}\), the vena contracta coefficient in the form of \(A_{vc} = \Pi_{vc} A_{or}\) yields:

(34)\[\Delta h = \frac{Q^2}{2g \Pi_{vc}^2 A_{or}^2}\]

Which, rearranged to solve for Q gives The Orifice Equation:

(35)\[ Q = \Pi_{vc} A_{or} \sqrt{2g\Delta h}\]
Such that:
\(\Pi_{vc}\) = 0.62 = vena contracta coefficient, in aide_design as pc.RATIO_VC_ORIFICE
\(A_{or}\) = orifice area- NOT contracted flow area
\(\Delta h\) = elevation difference between orifice and water level
Equations in aide_design:
pc.flow_orifice(Diam, Height, RatioVCOrifice) Returns flow through a horizontal orifice.
pc.flow_orifice_vert(Diam, Height, RatioVCOrifice) Returns flow through a vertical orifice. The height parameter refers to height above the center of the orifice.

There are two configurations for an orifice in the wall of a reservoir of water, horizontal and vertical, as the image above shows. The orifice equation shown in the previous section is for a horizontal orifice, but for a vertical orifice the equation requires integration to return the correct flow. You will explore this in the Flow Control and Measurement Design Challenge.

Section Summary

  1. Bernoulli vs energy equations: The Bernoulli equation assumes that energy is conserved throughout a streamline or control volume. The Energy equation assumes that there is energy loss, or head loss \(h_L\). This head loss is composed of major losses, \(h_{\rm{f}}\), and minor losses, \(h_e\).

Bernoulli equation:

(36)\[\frac{p_1}{\rho g} + {z_1} + \frac{\bar v_1^2}{2g} = \frac{p_2}{\rho g} + {z_2} + \frac{\bar v_2^2}{2g}\]

Energy equation, simplified to remove pumps, turbines, and \(\alpha\) factors:

(37)\[\frac{p_{1}}{\rho g} + z_{1} + \frac{\bar v_{1}^2}{2g} = \frac{p_{2}}{\rho g} + z_{2} + \frac{\bar v_{2}^2}{2g} + h_L\]
  1. Major losses: Defined as the energy loss due to shear between the walls of the pipe/flow conduit and the fluid. The Darcy-Weisbach equation is used to find major losses in both laminar and turbulent flow regimes. The equation for finding the Darcy friction factor, \(\rm{f}\), changes depending on whether the flow is laminar or turbulent. The Moody diagram is a common graphical method for finding \(\rm{f}\). During laminar flow, the Hagen-Poiseuille equation, which is just a combination of Darcy-Weisbach, Reynolds number, and \({\rm{f}} = \frac{64}{\rm{Re}}\), can be used
Darcy-Weisbach equation:

(38)\[h_{\rm{f}} = {\rm{f}} \frac{L}{D} \frac{\bar v^2}{2g}\]

For water treatment plant design we tend to use plant flow rate, \(Q\), as our master variable and thus we have.

(39)\[h_{\rm{f}} = {\rm{f}} \frac{8}{g \pi^2} \frac{LQ^2}{D^5}\]

\(\rm{f}\) for laminar flow:

(40)\[{\rm{f}} = \frac{64}{\rm{Re}} = \frac{16 \pi D \nu}{Q} = \frac{64 \nu}{\bar v D}\]

\(\rm{f}\) for turbulent flow:

(41)\[{\rm{f}} = \frac{0.25} {\left[ \log \left( \frac{\epsilon }{3.7D} + \frac{5.74}{{\rm Re}^{0.9}} \right) \right]^2}\]

Hagen-Poiseuille equation for laminar flow:

(42)\[h_{\rm{f}} = \frac{32\mu L \bar v}{\rho gD^2} = \frac{128\mu Q}{\rho g\pi D^4}\]
  1. Minor losses: Defined as the energy loss due to the generation of turbulent eddies when flow expands. Once more: minor losses are caused by flow expansions. There are three forms of the minor loss equation, two of which look the same but use different coefficients (\(K^{'}\) vs \(K\)) and velocities (\(\bar v_{in}\) vs \(\bar v_{out}\)). Make sure the coefficient you select is consistent with the velocity you use.

First form:

(43)\[h_e = \frac{\left( \bar v_{in} - \bar v_{out} \right)^2}{2g}\]

Second form:

(44)\[h_e = \frac{\bar v_{in}^2}{2g}{\left( {1 - \frac{A_{in}}{A_{out}}} \right)^2} = \,\,\, \frac{\bar v_{in}^2}{2g} \mathbf{K^{'}}\]

Third and most common form:

(45)\[h_e = \frac{\bar v_{out}^2}{2g}{\left( {\frac{A_{out}}{A_{in}}} -1 \right)^2} = \,\,\,\, \frac{\bar v_{out}^2}{2g} \mathbf{K}\]
  1. Major and minor losses vary with flow: While it is generally important to know how increasing or decreasing flow will affect head loss, it is even more important for this class to understand exactly how flow will affect head loss. As the table below shows, head loss will always be proportional to flow squared during turbulent flow. During laminar flow, however, the exponent on \(Q\) will be between 1 and 2 depending on the proportion of major to minor losses.
Head loss scales with: Major Losses Minor Losses
Laminar \(Q\) \(Q^2\)
Turbulent \(Q^2\) \(Q^2\)
  1. The head loss trick, also called the control volume trick, can be used to incorporate the ‘kinetic energy out’ term of the energy equation, \(\frac{\bar v_2^2}{2g}\), into head loss as a minor loss with \(K = 1\), so the minor loss equation becomes \(\left( 1 + \sum K \right) \frac{\bar v^2}{2g}\). This is used to be able to say that \(\Delta z = h_L\) and makes many equation simplifications possible in the future.
  2. Orifice equation and vena contractas: The orifice equation is used to determine the flow out of an orifice given the elevation of water above the orifice. This equation introduces the concept of a vena contracta, which describes flow contraction due to the inability of streamlines to make sharp turns. The equation shows that the flow out of an orifice is proportional to the square root of the driving head, \(Q \propto \sqrt{\Delta h}\). Depending on the orientation of the orifice, vertical (like a hole in the side of a bucket) or horizontal (like a hole in the bottom of a bucket), a different equation in aide_design should be used.

The Orifice Equation:

(46)\[Q = \Pi_{vc} A_{or} \sqrt{2g\Delta h}\]