## Control Volume (Integral) Technique

**:**

__Introduction__**Techniques for solving flow problems**

There are three different techniques for solving these types of problems:

**Control Volume Analysis**- This is just like in thermodynamics class
- One can calculate the
(total power output, total heat transfer, etc.)__gross properties__ - With this however, we do not care about the details inside the control volume (In other words we can treat the control volume as a "black box.")
**Differential Analysis**- In this technique, one solves differential equations of motion everywhere (i.e. The Navier-Stokes equation).
- Here we solve for
__all__of the details in the flow **Dimensional Analysis and Experiment**- This method is used when methods one or two are not possible.
- One uses wind tunnels, models, etc. to employ this method.

__Introduction to systems and control volumes__

**Consider a**

__System__:__system__of a fluid:

- Recall from thermo class, that a
is defined as a volume of mass of fixed identity.__system__ - Let m = the mass of the system

Let= the velocity of the system__V__

Let= the acceleration of the system__a__ - Now we can write three basic conservation laws which apply to this system. Note: These conservation laws apply directly to a system.
states that the mass of a system is constant.__Conservation of mass__

This can be written as the following equation:

In this equation m = the mass of the system.which is a restatement of__Conservation of linear momentum__

Newton's Second Law.- In equation form this is written as:

Where m__V__= the linear momentum of the system. - Note that this is the same as Newton's second law, it is just written a little differently.
- Using the above equation, we can obtain a form of the equation in terms of acceleration:

- This equation is obtained through expansion using the product rule.
- We know that ,

and

by the conservation of mass found above. - Using these ideas, we can then see that
__Conservation of Energy__- For this, use the First Law of Thermodynamics in rate form to obtain the following equation:
- Where E = the total energy of the system. In the above equation

is the rate of change of system energy. - is the rate of heat added
the system*to* - is the rate of work done
the system . Because work is done by the system, the negative sign is in the equation for the first law of thermodynamics.*by* - Now, these conservation laws must always hold
.__for a system__ __Conservation of Angular Momentum__

We will not have time to study this, but see the text for details.

- One can think of the system approach as the
, which if recalled is the description where we follow the individual chunks of the fluid.__Lagrangian description__ - However, the trouble with this is that we prefer to use the
, where we define a control volume with fluid moving through it.__Eulerian Description__ - In other words, if we define some control volume,
, the system will flow__fixed in space__it.__through__ - Our goal here is to try to write the conservation laws above in a form applicable to a
(Eulerian description) rather than to a__control volume__(Lagrangian description).__system__ - There is a way to do this, and it is called the
__Reynolds Transport Theorem__(R.T.T.) - Now our goal is to write all three of these conservation laws in terms of a control volume. But, to save us some work, let's not derive it three times, but just once for a general property B, and
__then__use our result for any of the three conservation laws: - In other words, let's write R.T.T. for some arbitrary property B.

Let B = some arbitrary property (vector or scalar) and per unit mass.

This way, we can later substitute B for mass, momentum, energy, etc. - We are then able to make the following substitutions for :
- For
, let B = m, and .*conservation of mass* - For
, , and .*conservation of momentum* - For
, B = E, and .*conservation of energy* - SEE the text for the derivation.

__Control Volume__:

__REYNOLDS TRANSPORT THEOREM (R.T.T.)__

- This is for some property B, and is also for a
control volume.__fixed__

Deforming the control volume is more complicated; see the text for details. - The equation for Reynolds Transport Theorem is written as the following:

(This is equation 3.12 in the text.)For a fixed control volume, the__OR__of differentiation and integration doesn't matter. So alternatively,*order*

(This is equation 3.17 in the text. ) **Comments**on the above equation:**The total rate of change of B following the system.**This is the term to which the conservation laws directly apply.

- The time rate of change of B within the control volume. This is due to unsteadiness.
- The net flux of B out of the control surface.

Due to, B changes because system moves to a new part of the flow field, where conditions are different. A circle was put around the integral to emphasize that this is an integral over the__convection__control surface. (This is not done in the text however).*entire*

- Sometimes and are used to
that these are*emphasize**volume*

andintegrals. (Again this is not done in the text).*area*

## Conservation of Mass using Control Volumes

__Conservation of Mass__

In the Reynolds transport theorem for conservation of mass, let __B = m__, i.e. B_{sys} = m_{sys} = mass of the system,

For our system, we know that

So, __Reynold's Transport Theorem__ (R.T.T.) becomes

OR,

- Let's analyze the last term in the equation:
- n = outward normal unit vector (perpendicular to surface area element dA)
**V**= velocity vector, which can be in any arbitrary direction- Recall, from definition of dot product:
- Outlets:
**V**is__outward__(mass leaving C.V.) - Inlets:
**V**is__inward__(mass entering C.V.) - Now, it turns out that
= mass flow rate outward through a surface

so, if we integrate over the entire surface,

and the conservation of mass equation becomes:

*this is the most useful form. - e.g.
__Given__: A rigid tank of volume V with p = p_{0 }at t = 0. Air is pumped in at constant mass flow rate isothermally.

Find: p in tank as a function of time

__Solution:__first, draw a C.V. inside the entire tank.

Now use our conservation of mass equation.

assume density is equal everywhere in the tank, and only varies with time. There are no outlets and the only term remaining is the mass term at the inlet: This is a differential equation we must solve by separating the variables and integrating from t=0 to some arbitrary t.

Finally, use the ideal gas law to get the pressure. Thus,

(NOTE: pressure varies linearly with t)**Special Case**:For a steady flow, nothing is a function of time, the d/dt term in the conservation of mass equation goes away.__STEADY FLOW__

in other words,__"What goes in must come out"__If more is going in than out, mass will be accumulating with time inside the control volume, and it would not be steady state!! ( and you need to include the unsteady volume term)

__( density = constant)__

**Incompressible Steady Flow**- Since the density is constant,

take density outside of the integral

So the integral term is equal to Q, the volume flow rate

i.e. where mass flow rate = density X volume flow rate

So the conservation of mass equation becomes:

"volume flow rate it" = "volume flow rate out"

(for incompressible flow)**Another simplification: 1-D inlets & outlets**Consider an outlet:

Where:- 1-D outlet means velocity is parallel to
**n**, the unit normal vector. **V**= constant across the outlet.- Density = constant across the outlet

_{out}=**V**A. Therefore, at an__inlet__

and we usually say:

….since the negative sign is accounted for in the conservation of mass equation,

i.e.

- 1-D outlet means velocity is parallel to
**Average Velocity**

If an inlet or exit is

__not__1-D, we can still use but**V**_{AV}must be the__average velocity__- e.g.: At an outlet, for incompressible flow ( density = constant)

Let us define:but for most problems

**V**is parallel to**n**……

An equivalent 1-D outlet will have the same mass flow rate as the actual outlet.

- e.g.: At an outlet, for incompressible flow ( density = constant)
3.14 in text:**EXAMPLE**

A container of water,

__Given:__**V**_{AV}, Q_{3}, D_{1}, D_{2}and h equal constants.

V_{1}= 3 m/s; Q_{3}= 0.01 m^{3}/s; h = constant; D_{1}= 0.05m; D_{2}= 0.07m__Find:__Average exit velocity V_{2}__Solution:__Use conservation of mass. First draw the C.V. shown.

Since it is steady and incompressible,

Using the above equations, solve for V

_{2}:

## Conservation of Momentum using Control Volumes

**Conservation of Linear Momentum** Recall the conservation of linear momentum law for a system:

In order to convert this for use in a control volume, use RTT with B = m**V**, beta = **V **

we get:

NOTE: Recall that at any instant of time t, the system & CV occupy the SAME physical space.

So, the forces of the system are the same at the forces of the control volume at a given instant.

- For a
**fixed control volume**we have the following equation:

This is a vector equation so it has three components. - First, let us consider the component in the
. We will drop the cv subscript since it is understood. The conservation of linear momentum equation becomes:__X-direction__ - Next, let us consider the component in the
. The conservation of linear momentum applied to the y-direction becomes:__Y-direction__ - Similarly, the conservation of momentum could be applied to the
.__Z-direction__ - It is now time for a few simplifications for the right hand side of the CLM equation.
- First, let us assume that we have one dimensional inlets and outlets. This implies that our velocity vector
__V__is parallel to our normal to the surface vector__n__. We also assume that the velocity is constant across the inlet or outlet surface.

Assuming density is constant, we can rewrite the last term in our**CLM**equation:

This holds true because__V__dot__n__= V for outlets and -V for inlets. We can make a further simplification if we notice that the definition for mass flow rate is (density)(velocity)(area). Therefore our equation for the flux term becomes:

- If the flow is
**steady**we can drop the (d/dt) term. Most of the problems we deal with having steady flow also have one dimensional inlets and outlets with constant density. Our equation for conservation of linear momentum now becomes:

Notice that this is a vector equation. Therefore we can break this up into three components. Since the velocity vector = (u,v,w) and the force vector = (F_{x},F_{y},F_{z}), our equation can be rewritten into three equations:

- Now lets return to the left side of the CLM equation (the Force term). This term represents the sum of all the forces acting on the control volume. There are several types of forces that can act on our control volume.

__Body Forces__:

Gravity is a body force that acts in the negative z-direction.__Surface Forces:__These forces include:- Pressure-contributes only normal forces.
- Viscous-viscous and frictional forces contribute normal and shear forces.
- Other-forces due to control volume cutting through bolts or struts.
- Let us revisit the
**pressure force**. Recall that P is always defined as__positive inward__and__n__is always defined as__positive outward__.. These two directions are opposite. Therefore it turns out that

You can use either**absolute or gage**pressure in this equation as long as you are consistent everywhere! To prove this, consider an arbitrary control volume with a certain pressure field. Suppose there is a high pressure at two locations and atmospheric pressure everywhere else. Realizing the definition of gage pressure, Pg = Pabsolute - Patmospheric, we can subtract atmospheric pressure from everything. We realize that the sum of the forces due to pressure is the same in both cases. The**net**pressure force__will not change__. This holds true because any shape in a uniform pressure field has zero net pressure force.

- Now let us return to the
**viscous forces**.

This force is difficult to calculate because it has a normal and a tangential component acting in an arbitrary direction. In order to solve for this term you would have to integrate the shear stress over the entire control surface. Fortunately we usually do not have to integrate for this term. We can often choose a "wise" control volume to eliminate this difficult term. This is done by picking a control volume so that no part of the control surface has viscous forces acting on it. We could also keep the viscous term as are unknown term and solve for it by computing our other terms. - Finally let us consider
**other surface forces**. These forces include whatever is left. In other words, any force acting on the control volume that is not accounted for. Some examples of these types of forces include struts, bolts, cables, ropes in tension, forces holding a control volume in place, etc. - The final form of the momentum equation for steady flow with one dimensional inlets and outlets is:
__An incompressible jet has an exit pressure equal to the ambient pressure__. The ambient pressure is usually atmospheric. Therefore, when you draw your control volume, it is wise to slice through the jet at the exit plane. This will cause the pressure to be ambient everywhere.

__The pressure is approximately ambient in any slice through an incompressible jet__. As you move down the stream, the pressure is still equal to the ambient pressure.__The mass flow rate and momentum flow rate equal zero across a streamline__. Therefore it is often advantageous to pick a control surface so as to run along a streamline.

When you do this__V__dot__n__is equal to zero because they are perpendicular. This will cause the__flux terms__(mass, momentum, and energy) to be zero along the streamlines.__Mass flow rate and momentum flow rate equal zero along a solid wall__. This is true since a solid wall is always a streamline. No mass can go through a solid wall. Note that the pressure and viscous force terms may not be zero along streamlines, however sometimes it is "wise" to pick a control volume along a solid wall.

Notice that the __V__ dot __n__ term is a scalar, not a vector.

*Here are some rules when using the momentum equation.*

Example Problems

**Given**: A water jet of velocity V_{j} and thickness D_{j} impinges on a turning block which is held in place by a force F_{x}, as shown in the sketch.

As the water leaves the block, the round jet flattens out and slows down due to friction along the wall. The water turns a full 180 degrees and flattens into a rectangle shape of thickness D_{j}/6 and width 10D_{j} in cross section. The flow is steady.

a) **Find**: V_{e}, the exit velocity of the jet.

**Solution**:

- First, pick a control volume which cuts through the inlet and exit, but does not include the turning block.
- Use conservation of mass for steady 1-D flow:

b) **Find**: F_{x}, the force required to hold block in place.

**Solution**:

- Choose a new control volume because friction along the wall would be very hard to evaluate with the previous C.V. Pick a C.V. which cuts through the inlet and outlet and that includes the whole turning block. The wisest choice of control volume is one which cuts through where the force F
_{x}acts on the turning block, since F_{x}is our unknown. - Write momentum equation for steady 1-D inlet/outlet in the x-direction:

Examine each term on the left: (no grav in x-dir), (p = p_{a}everywhere), (wise choice of C.V), and - Thus, the momentum equation reduces to
- Now, recall, at an inlet or exit. Here at the inlet,

- So, solve for :

Or, finally, , which is our final answer. - Note: This is not merely an academic exercise. A testing facility at NASA Langley Research Center uses this principle to propel a test rig along a track to test airplane landing gear, etc. The experimental facility, called the Aircraft Landing Dynamics Facility, uses a high pressure water jet that hits a turning bucket on the test cart, much like in the above problem. The cart is propelled from zero to 250 miles per hour in two seconds flat!

Momentum Flux Correction Factor

- Again, just like in the conservation of mass equation, even if we don't really have one-dimensional (uniform) inlets and outlets, we would still like to use the simplified version of the conservation of momentum equation. I.e. we would like to use V
_{av}instead of V in the equation.

- The equivalent profile and actual profile thus have identical mass flow rates, .
- What about momentum flow rate, or momentum flux, MF? Do actual and equivalent profiles have the same momentum flux?
! It turns out, after integration, that:*No* - So, we can't just substitute V
_{av}for V in the momentum equation or we will get the wrong answer. So, instead, let's introduce a momentum flux correction factor, :

Define:

It turns out that: ,

and, . - Examples:

1) 1-D uniform flow, by definition

2) fully developed laminar pipe flow,

3) fully developed turbulent pipe flow, - Now, we can use the one-dimensional form of the momentum equation, but with these momentum flux correction factors thrown in:

For a fixed control volume with steady flow,

, and - Fortunately, most problems in real life applications are
__turbulent,__not laminar, and we can usually neglect the momentum flux correction factors since is close to 1.0.

Most Useful Form of the Momentum Equation

- For steady flow with a fixed control volume, the most useful form of the momentum equation is thus:

, where momentum flux correction factor, and .

Example Problem

**Given**: Consider incompressible flow in the entrance of a circular tube.

The inlet flow is uniform u_{1} = U_{0}. The flow at section 2 is developed laminar pipe flow.U_{0}, p_{1}, p_{2}, R, and , are also known. At section 2, .

**Find**: Total friction force on the fluid from 1 to 2.

**Solution**:

- First draw a C.V. Again cut through the inlet and outlet 1 and 2. But now, we want to cut along the walls since F
_{friction}is our unknown. - Draw all the forces acting on the C.V. in the x-direction.
- Simplifications: Assume steady, incompressible, uniform one-dimensional inlet, but not a uniform one-dimensional outlet.
- So, let's use the modified momentum equation for non-1-D inlets/outlets in the x-direction:
- Now look at each term:

(no grav. in x-dir)

(this is our unknown)

(no struts, bolts, etc.)

(inlet is 1-D so )

: This term is the hardest. First of all, what is u_{av}? I.e. what is the equivalent 1-D velocity? Well, since the area at the outlet is the same as that at the inlet, u_{av}at the outlet has to equal U_{0}. - What is at the outlet? By definition, . You can plug the equation in and integrate [Try it on your own]. You will get (the value for laminar pipe flow).
- So, the outlet momentum flux term becomes
- Finally, then sum everything and solve for the unknown:

- Notice, if is ignored, i.e. (set ), this last term would cancel completely, and the answer would be wrong.

**More Example Problems**

*Problem # 3.58 in the text:*

__Given__: Cart with water jet, deflector, as shown

Known in this problem are the jet area A, the average velocity V_{av}, the jet deflection angle, , and the momentum. flux correction factor of the jet, . Also, frictionless wheels are assumed.

__Find__: Tension in cable at time t=0.

__Solution__:

- As always, the first step in any control volume problem is to pick and draw a control volume. The forces and coordinate system have been labeled on the sketch. Here it is appropriate to slice through the jet and slice through the unknown force T (tension in cable), as shown in the sketch.
__Simplifications__:

Incompressible? Yes (water is the fluid, which is approximately an incompressible liquid)

Steady? - No, not really, since the water level in the tank is falling. But all we are asked for is T at time t=0, so think of this as a__pseudo-steady__problem in order to find T. Note: The falling water level will have no effect on the x-momentum equation anyway, since its velocity is vertical.- Now, use the x-component of the momentum equation in it steady form:

There is only one outlet, at which

.

The quantity u_{av}in the momentum flux outlet term in the above equation needs to be considered carefully. This is not V_{av}, but rather is the x-component of V_{av}. From a little trig one can see that

Thus the outlet term on the right hand side of the momentum equation is

and the final form of the x-momentum equation is (solving for T): - Plug in the numbers last: The density of water at room temperature is 998. kg/m
^{3}, the jet velocity, V_{av}is 8 m/s, the cross-sectional area A is that of a circle of diameter 0.04 m, the jet angle is given as 60^{o}, and the jet momentum flux correction factor is 1.0 for a uniform jet. This yields:

Or T = 40.1 N. __Question__: At what angle, is the tension a maximum?__Answer__: When the jet deflection angle is zero, i.e. the jet simply exits horizontally into the atmosphere. This is clearly seen in the above equation since the cosine of zero is unity.

*Example A water-mounted fire pump (an old exam problem.)*

__Given:__ A pump is anchored to the ground as shown, with V_{j}=35.0 m/s and d_{j}=3.00 cm. Assume the jet has a fully developed turbulent pipe flow profile at its exit.

picture wat_pump.gif

__Find:__ Horizontal force required to hold platform in place.

__Solution:__

- The first step in any control volume problem is to pick and draw a control volume. As shown, the control volume should slice through jet exit and through the bolts or whatever is holding the platform solidly to the ground. The inlet to the control volume is most easily taken at the surface of the water, where the pressure and velocity are known (p = p
_{a}and V = 0 at the surface). Note: When this problem was given as an exam question, lots of students took their control volume inlet at the pipe inlet. This makes the problem more difficult. Remember the first rule about selecting a control volume - Be Lazy! - Apply the x-momentum equation:

The gravity term on the left hand side is zero because gravity does not act in the x direction. Likewise, the pressure term is zero because everywhere on the control surface above the water, the pressure is atmospheric (including the portion of the control surface that slices through the jet). Below the water, the pressure is hydrostatic, and whatever pressure force is exerted on the left side of the control surface is exactly balanced by that on the right side of the control surface. In the control volume selected, there is no net viscous force acting on the control volume. Thus the first three terms on the left hand side of the x-momentum equation are zero. The only "other" force acting on the control surface is the force of the ground acting on the platform, as shown in the sketch. The direction of this force is assumed to act to the right; if this is wrong, the result will be negative.

On the right hand side (the momentum flux terms), there is no x-component of velocity at any inlet, so the last term on the right is zero. There is only one outlet, with u_{av}equal to the x-component of the velocity vector of the jet. Thus, the x-momentum equation reduces to

This is the answer in variable form. - Lastly, plug in the numbers, noting that for a fully developed turbulent pipe flow the momentum flux correction factor is around 1.02:

or finally,__the horizontal force required to hold platform in place is 799. N__.

*Problem # 3.51 in the text:*

__Given:__ A turbine wheel, powered by a water jet, as shown in the sketch (at time t = 0)

The turbine is spinning at a constant rotational speed.

__Assumptions:__

- 1-D inlet and outlet
- neglect friction in the turning bucket

(a) __Find:__ The force of the turning bucket on the turbine wheel at this instant of time.

__Solution:__

- Pick a control volume. Here we will pick a
__moving control volume__, which makes things a little tricky. When the control volume is moving, one must use the__relative velocity__, i.e.__relative to the moving C.V.__ - At the inlet, then, the relative velocity at the inlet is obtained by subtracting the velocity of the C.V. from the absolute velocity at the inlet, i.e.

where Vj is the magnitude of the absolute velocity at the inlet, and its direction is in the positive x-direction. The second term on the right represents the velocity of the control volume in an absolute reference frame, which is subtracted as shown. The x-component of this relative inlet velocity is simply the magnitude of the above velocity vector, since the vector acts only in the positive x-direction, i.e.

- Now apply conservation of mass to this moving control volume. Note: Here,
*relative*velocities must be used, since they represent the velocities actually entering and leaving the control volume as it moves along.

Note that the areas on either side of the equation are identical by the original assumption of negligible viscous effects. Also the incompressible assumption causes the density to drop out. At time t = 0, when the turning bucket is on the top of the turbine wheel as shown, the velocity of the inlet is in the positive x-direction, while that of the outlet is in the negative x-direction. Thus, the x-component of the relative outlet velocity is simply the negative of that of the relative inlet velocity, i.e.

In other words, the turning bucket changes the direction of the water jet by 180 degrees, but it does not change the magnitude of the jet velocity. (If friction along the walls of the turning vane was taken into account, the jet velocity magnitude would change as well.) __Now apply conservation of x-momentum__for our moving C.V.

Notice that the*relative*x-velocity components are used in the momentum flux terms on the right hand side since this is a moving control volume. The gravity term on the left hand side is zero because gravity does not act in the x direction. Likewise, the pressure term is zero because everywhere on the control surface, the pressure is atmospheric (including the portions of the control surface that slice through the jets). In the control volume selected, there is no net viscous force acting on the control volume. (This would be true even if viscous forces along the turning vane were not ignored, because the control surface does not pass along the turning vane wall.) Thus the first three terms on the left hand side of the x-momentum equation are zero. The only "other" force acting on the control surface is the force of the turbine wheel acting on the turning bucket, which is equal and opposite to the force of the turning bucket acting on the turbine wheel, as shown in the sketch. The direction of the force of the bucket on the wheel is assumed to act to the right; if this is wrong, the result will be negative.

On the right hand side (the momentum flux terms), there is only one inlet and one outlet. At both the inlet and the outlet, the x-component of the relative velocity vector of the jet is known from the conservation of mass analysis above. Since friction is being neglected, it is assumed that the momentum flux correction factors are unity. Thus, the x-momentum equation reduces to

This is the force of the turning bucket on the turbine wheel.

(b) __Find:__ The power, P, delivered to the wheel at this instant of time (t = 0).

__Solution:__

- Power on a rotating wheel is defined as the torque times the angular velocity of the wheel. The torque due to this turning bucket is the force just derived times the radius fo the wheel, R. Thus,

(c) __Find:__ The angular velocity which provides the maximum power to the wheel.

__Solution:__

- To find the maximum power, take the derivative of power with respect to angular velocity.

- The values of angular velocity at which this derivative is zero then represent conditions where the power is either a minimum or a maximum. There are two such roots, i.e.

- The first of these represents the case where the absolute jet speed and the speed of the turning bucket are the same. Under such conditions, no power is transferred to the turbine wheel at all. (This can be verified easily, i.e. P = 0 when the first root is plugged into the result of part (b) above.) In other words, the first root is a minimum. The second root is therefore the desired one; i.e. the angular velocity which supplies the maximum power to the turbine wheel is that in which the turbine wheel spins three times slower at its rim than the jet speed, i.e.

## 0 comments:

Post a Comment