Calculates the friction factor for the current guess

Assignment Help Physics
Reference no: EM13278267

Question 1

Introduction

You are designing a pipe network system that transfers water from the upper pipe to the lower pipe. Note that Figure 1 is a plan view and the elevation is constant across all pipes. The static pressure difference between points A and D is designed to be PA - PD = 3 atm (1 atm = 1 standard atmospheric pressure = 101.3 kPa). It is necessary to ensure that the speed of the flow through every pipe is at least 2 m/s so that there is no sediment build-up. Determine if this is the case.

96_Designing a pipe network system.png

Theory

The change in pressure between two points along a streamline (a flow path) is modelled by the Bernoulli equation

230_Designing a pipe network system1.png

where P is the static pressure, Ρ the density, V the speed, g gravitational acceleration, h the elevation and ΔΡloss is the reduction in pressure due to any losses in the system. The most important loss (and the only one to be accounted for here) is caused by friction:

1108_Designing a pipe network system2.png

where f is the Darcy-Weisbach friction factor. Fluid flow is governed by the continuity equation (which is conservation of mass); for incompressible (constant-density) flow, this is:

Q = VA = Const

where Q is the volume flow rate (m3/s) and A is the cross-sectional area. Incompressible flow is a good assumption for liquids. A consequence of Eq. (3) for Eq. (1) is that if the cross-sectional area is constant for a given pipe, the flow speed at the start is equal to the speed at the end and can be defined based on the pipe ID rather than an end-point ID.

Pipe networks can be considered to be equivalent to electrical circuits in series and parallel, with pressure change equivalent to potential difference and volume flow rate equivalent to current.

However, the resistance cannot be treated as constant: it is a non-linear function of the flow speed. The rules of potential difference and current are still maintained however:

a) The pressure change (potential difference) across multiple branches in parallel is equal (e.g. PB - PC is the same regardless of whether pipe 2 or 3 is taken).

b) The sum of the volume flow rates (currents) entering a junction is equal to the sum of the volume flow rates exiting a junction (e.g. Q1 = Q2 + Q3).

Because of the non-linear nature of the system [flow rate is squared in Eq. (1) and f is a non-linear function of Q], iteration is required to determine the flow speeds in each pipe section. To do this, certain constraints can be applied based on the rules of potential difference and current:

1. The pressure change between points A and D is known and is equal to the sum of the pressure changes along all pipes in series that connect points A and D. This constraint should be used only once.

2. The pressure change across pipes in parallel is the same for each pipe.

3. The sum of volume flow rates entering junction B is equal to the sum of volume flow rates exiting B; similarly for C.

These constraints provide a set of equations that can be solved for the speeds in each pipe.

To calculate the friction factor, the common formula that is used is the Colebrook formula (Colebrook 1938-39):

113_Designing a pipe network system3.png

where e is the absolute roughness of the pipe wall and the Reynolds number is

with v (the Greek letter "nu") the kinematic viscosity. Equation (4) is only valid for turbulent pipe flow (Re > 2300), otherwise f = 64/Re. The Moody diagram (Figure 2) is a standard method of determining f (the left axis) for hand calculations; the range of possible values for f given the different blue lines that represent the range of plausible pipe roughnesses can be seen.

Requirements

For this assessment item, you must produce MATLAB code which:

1. Iterates until convergence is reached (the previous guess for the velocities is not significantly different from the current guess).

2. Calculates the friction factor for the current guess of velocities.

3. Reports the flow velocity in each pipe to the Command Window in addition to confirming that the speeds meet the required level

4. Validates the code by checking that the volume flow rates in the branches are consistent.

5. Verifies the code by ensuring that the friction factor calculation is valid.

6. Displays to the Command Window a brief discussion (fewer than 5 lines) stating the value you selected for the under-relaxation factor and why you selected that value.

7. Has appropriate comments throughout.

An important component of quality assurance is to test (verify) that each function works correctly:

8. Write a test program that supplies block of code with a known input and confirm (verify) that the output is correct. (Writing blocks of code as functions makes the code more transparent, more portable and easier to test!)


Attachment:- matlab.rar

Reference no: EM13278267

Questions Cloud

Calculate reactive power of motor and the net power factor : A 220V , 13.47KVA, 0.707 lagging power factor load is in parallel with a three phase wye- connected round rotor synchronous motor having a synchronous reactance of 1.27W. The motor draws 33KW at apower angle of 30 degree.
Explain a strong electrolyte and weak electrolyte : Classify each compound as a strong electrolyte, weak electrolyte, or non-electrolyte. HI HClO4 CH3OH HC2H3O2
How much work is done by the man : A man pushes a lawn mower on a level lawn with a force of 200N. If 40% of this force is directed downward, then how much work is done by the man
Derive the state diagram for fsm that has input and output : Derive the state diagram for an FSM that has an input w and an output z. The machine has to generate z = 1 when the previous four values of w were 1001 or 1111; otherwise, z = 0. Overlapping input patterns are allowed.
Calculates the friction factor for the current guess : Iterates until convergence is reached (the previous guess for the velocities is not significantly different from the current guess and calculates the friction factor for the current guess of velocities.
Define compound set has the highest and lowest boiling point : Which compound of each set has the highest and lowest boiling point? Decribe the factors that mitigate these values. A) CH4, CH3CH3, CH3CH2CH3, CH3CH2CH2CH3 B) CH3Cl, CH2Cl2, CHCl3, CCl4
Determine the direction of the velocity of the two vehicles : A 1100kg car traveling south at 24 m/s with respect to the ground collides with and attaches to a 2550kg delivery truck traveling east at 16 m/s. Determine the direction of the velocity of the two vehicles
What quantity of ice melts : A 4.40-g lead bullet at 32.0°C is fired at a speed of 230 m/s into a large block of ice at 0°C, What quantity of ice melts
Explain the monobromination of 4-ethylheptane : Consider the monobromination of 4-ethylheptane using Br2 and UV light. How many unique alkyl halide products could be formed during this reaction

Reviews

Write a Review

 

Physics Questions & Answers

  What is an angular displacement of the propeller

A flight attendant pulls her 72.0 N flight bag a distance of 282 m along a level airport floor at a constant speed. The force she exerts is 33.0 N at an angle of 52.0° above the horizontal.

  Illustrate where will the end of spring line up with marks

A spring with 58 hangs vertically next to a ruler. The end of the spring is next to the 11- mark on the ruler. If a 2.5- mass is now attached to the end of the spring, illustrate where will the end of the spring line up with the ruler marks?

  What is the smallest force f

A 21.8-kg box rests on a frictionless ramp with a 15.7° slope, what is the smallest force F the mover will have to exert to move the box up the ramp

  What is the apparent weight of a 94 kg astronaut

An insulated Thermos contains 210 cm3 of hot coffee at 88 °C. You put a 19.00 g ice cube at its melting point to cool the coffee. By how many degrees has the coffee cooled once the ice has melted and equilibrium is reached? Treat the coffee as though..

  How fast was the ball travelling

A 70.0 kg football player leaps straight up into the air (with no horizontal velocity) to catch a pass. If the ball was moving horizontally when it was caught, how fast was the ball travelling

  Determine the minimum area of the capacitor plates

A parallel-plate capacitor is constructed using a dielectric material whose dielectric constant is 3.90 and whose dielectric strength is 1.80 multiply.gif 108 V/m. Find the minimum area of the capacitor plates

  Find the torque on the dipole

An electric dipole consists of charges +e and -e separated by 0.76 nm. It is in an electric field of strength 2.50 10^4 N/C, What is the torque on the dipole when it is perpendicular to the field

  What value does this experiment yield for plancks constant

When a wavelength is 300nm, its stopping potential is 1.10V. When it is changed to 200nm, the stopping potential is 3.06V. What value does this experiment yield for Planck's constant

  Calculate the power emitted by the human body

The emissivity of the human skin is 97.0 percent. Use 35.0 °C for the skin temperature and approximate the human body by a rectangular block with a height of 1.58 m, Calculate the power emitted by the human body

  Determine the effective spring constant of the molecule

A molecule of DNA (deoxyribonucleic acid) is 2.28 µm long, Determine the effective spring constant of the molecule

  How many stairs must woman climb to work off her breakfast

A 65-kg woman cheats on her diet and eats a 612-Calorie (612 kcal) jelly doughnut for breakfast. If the human body is only 21% efficient in converting chemical energy to mechanical energy, how many stairs must the woman climb to work off her breakf..

  Find the magnitude and direction of the electric force on a

Particle A carries a charge of 5.0 uC, Particle B is 12.0cm, Find the magnitude and direction of the electric force on A due to B

Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd