laminar below Re ≈ 2000, turbulent above Re ≈ 4000
Mean velocity
V = Q/A, A = π·D²/4
Darcy–Weisbach
hf = f·(L/D)·V²/(2g)
Minor losses
hm = ΣK·V²/(2g)
Typical K values
entrance 0.5, exit 1.0, 90° bend ≈ 0.3
Pressure drop
Δp = ρ·g·hf
Two open reservoirs
Δz = (f·L/D + ΣK)·V²/(2g)
Free outlet
H = (1 + f·L/D + ΣK)·V²/(2g)
Open-channel flow (Manning)
Manning equation
V = (1/n)·R2/3·S1/2
Flow rate
Q = A·V
Hydraulic radius
R = A/P (P = wetted perimeter, free surface excluded)
Rectangle
A = b·y, P = b + 2y
Trapezoid (side slope z:1)
A = (b + z·y)·y, P = b + 2y·√(1 + z²)
Circular pipe, full or half full
R = D/4
Required slope
S = (n·Q/(A·R2/3))²
Pumps & pump power
Hydraulic power
Ph = ρ·g·Q·H = Q·Δp
Shaft power
Pshaft = Ph/ηpump
Electrical power
Pel = Pshaft/ηmotor
Pump head
H = Δz + Δp/(ρ·g) + hL
System losses
hL = (f·L/D + ΣK)·V²/(2g)
Time value of money
Future value
FV = PV·(1 + i)n
Present value
PV = FV/(1 + i)n
Several periods per year
i = r/m, N = m·n
Annuity, present value
PV = PMT·[1 − (1 + i)−n]/i
Annuity, future value
FV = PMT·[(1 + i)n − 1]/i
Method
Draw a timeline with every cash flow and its timing.
Match the rate to the period: a monthly rate for monthly payments.
Use FV = PV·(1 + i)ⁿ for single amounts and the annuity formulas for equal payments.
Check the sign and the order of magnitude of the answer.
NPV & IRR
Net present value
NPV = −C0 + Σ CFt/(1 + r)t
Equal cash flows
NPV = −C0 + CF·[1 − (1 + r)−n]/r
Internal rate of return
the r at which NPV = 0
Decision
accept when NPV > 0 (or IRR > required return)
Method
List all cash flows with their year; the investment at year 0 is negative.
Discount each cash flow: CFt/(1 + r)t.
Add them up to get the NPV. Positive means the project earns more than r.
For the IRR, find the rate that makes NPV zero by trial or interpolation.
Normal distribution
z-score
z = (x − μ)/σ
Below x
P(X < x) = Φ(z)
Above x
P(X > x) = 1 − Φ(z)
Between
Φ(z2) − Φ(z1)
Percentile
x = μ + z·σ
Method
Convert the value to a z-score: z = (x − μ)/σ.
Look up Φ(z), the area to the left of z.
For the area to the right use 1 − Φ(z); for an interval subtract the two left areas.
For a percentile, work backwards: find z for the area, then x = μ + z·σ.
2Practice problems
1.Beam support reactionsMedium
The simply supported beam has a pin at A and a roller at B. Find the vertical support reaction RB.
L = 8 mP1 = 30 kN at 4 mP2 = 30 kN at 7 m
RB = kN
2.Beam support reactionsHard
The simply supported beam has a pin at A and a roller at B. Find the vertical support reaction RB.
L = 8 mP = 30 kN at 7 mw = 6 kN/m from 0 to 7 m
RB = kN
3.Truss: method of jointsMedium
Use the method of joints to find the force FBC in member BC. Give tension as positive and compression as negative.
span = 6 mheight = 4 mB is 2 m from Avertical load = 45 kN
FBC = kN
4.Truss: method of jointsHard
Exactly one member of this truss carries no force under these loads. Which is the zero-force member? Click it in the diagram.
D: 20 kN ↓F: 30 kN →G: 15 kN ↓H: 15 kN ↓
1member AB
2member BC
3member CD
4member DE
5member FG
6member GH
7member AF
8member HE
9member BF
10member CG
11member DH
12member BG
13member DG
5.Shear force & bending momentMedium
Find the bending moment Mx at the section x = 2.5 m from support A. Use sagging moments positive and upward shear on the left face positive.
L = 5 mw = 5 kN/m from 0 to 5 mx = 2.5 m
Mx = kN·m
6.Shear force & bending momentHard
Find the bending moment Mx at the section x = 7 m from support A. Use sagging moments positive and upward shear on the left face positive.
L = 8 mP = 40 kN at 1 mw = 3 kN/m from 0 to 8 mx = 7 m
Mx = kN·m
7.Beam deflectionMedium
A cantilever of length L is fixed at A and carries a point load P at a distance a = 4.5 m from A. Find the deflection δB at the free end B.
P = 5 kNa = 4.5 mL = 5 mE = 200 GPa (steel)I = 84 × 10⁶ mm⁴
δB = mm
8.Beam deflectionHard
A simply supported beam carries a uniformly distributed load w. Its deflection may not exceed span/360. Find the smallest second moment of area Ireq that meets this limit.
w = 12 kN/mL = 5 mE = 70 GPa (aluminium)
Ireq = × 10⁶ mm⁴
9.Column bucklingMedium
A column of length L made of aluminium is fixed at both ends. Using Euler’s formula, find its critical buckling load Pcr.
L = 5 mE = 70 GPa (aluminium)b × h = 20 × 80 mm
Pcr = kN
10.Column bucklingHard
A column of length L made of steel is fixed at the base and free at the top. Find its slenderness and the critical buckling stress σcr from Euler’s formula.
L = 1.5 mE = 200 GPa (steel)d = 50 mm
σcr = MPa
11.Centroid & second moment of areaMedium
A circular tube has outside diameter D = 190 mm and inside diameter d = 170 mm. Find its radius of gyration r about the horizontal centroidal axis.
D = 190 mmd = 170 mm
r = mm
12.Centroid & second moment of areaHard
A welded I-section has a top flange 140 × 10 mm, a web 290 mm high and 12 mm thick, and a bottom flange 170 × 16 mm (width × thickness). Find its second moment of area Ix about the horizontal axis x through its centroid.
top flange: 140 × 10 mmweb: 12 × 290 mmbottom flange: 170 × 16 mm
Ix = × 10⁶ mm⁴
13.Effective stress in soilMedium
For the soil profile shown, find the pore water pressure u at point X.
A 5 m thick overconsolidated clay layer has Cc = 0.4, Cr = 0.03 and e₀ = 0.6. At its mid-depth the effective vertical stress is 105 kPa and the overconsolidation ratio is OCR = 2.5. A new building adds 40 kPa of vertical stress at mid-depth. Find the primary consolidation settlement S in mm.
A clay layer 3 m thick has a coefficient of consolidation cv = 0.8 m²/year. It lies on impervious rock, with sand above it, so water drains out only at the top. For 60% average consolidation the time factor is Tv = 0.287. How many years t does the layer take to reach 60% consolidation?
A vertical rectangular gate 1.5 m wide sits in the wall of a tank of water, with its top edge 2 m below the surface. Find the resultant hydrostatic force F on the gate.
A vertical rectangular gate 1 m wide sits in the wall of a tank of oil, with its top edge 2 m below the surface. Find the depth ycp of the centre of pressure below the surface.
d = 2 mH = 2 mb = 1 mρ = 880 kg/m³g = 9.81 m/s²
ycp = m
19.Pipe friction & head lossMedium
Water (ρ = 1000 kg/m³) flows at 2.5 m/s through a horizontal pipe, 1000 m long, with an inside diameter of 300 mm. Use the Darcy friction factor f = 0.022 and g = 9.81 m/s². Find the pressure drop Δp over the pipe in kPa.
Two open reservoirs are joined by a pipe 300 m long with an inside diameter of 200 mm and a Darcy friction factor f = 0.018. The water surface of the upper reservoir is 37 m above that of the lower one. The minor losses are a square-edged entrance (K = 0.5), 2 bends (K = 0.3 each), a valve (K = 2) and the exit into the lower reservoir (K = 1.0). Use g = 9.81 m/s². Find the flow rate Q in L/s.
A rectangular clean earth channel is 4 m wide. Water flows uniformly 0.8 m deep. Manning n = 0.022 and the bed slope is S = 0.002. Find the flow rate Q in m³/s.
b = 4 my = 0.8 mn = 0.022S = 0.002
Q = m³/s
22.Open-channel flow (Manning)Hard
A trapezoidal grass-lined earth channel with a bottom width of 5 m and side slopes of 2:1 (horizontal:vertical) must carry 42 m³/s in uniform flow at a depth of 2.3 m. Manning n = 0.03. Find the required bed slope S (in m/m).
Q = 42 m³/sb = 5 my = 2.3 mz = 2n = 0.03
S = m/m
23.Pumps & pump powerMedium
In a test, a pump delivers 40 L/s of water (ρ = 1000 kg/m³) against a head of 12 m while its shaft takes 5.7 kW. Use g = 9.81 m/s². Find the pump efficiency η in percent.
Q = 40 L/sH = 12 mPshaft = 5.7 kWg = 9.81 m/s²
η = %
24.Pumps & pump powerHard
A pump moves 15 L/s of water (ρ = 1000 kg/m³) from an open lower reservoir to an upper tank whose water surface is 7 m higher. The pipe is 200 m long with an inside diameter of 100 mm and a Darcy friction factor f = 0.018; the minor-loss coefficients add up to ΣK = 8 (entrance and exit included). Use g = 9.81 m/s². The pump efficiency is 75%. Find the shaft power Pshaft in kW.
You invest €20 000 at 7% per year with continuous compounding. Find the value FV after 12 years.
PV = €20 000r = 7%t = 12
FV = €
26.Time value of moneyHard
You want €75 000 in 13 years. You save the same amount at the end of every month in an account paying 4.2% per year, compounded monthly. Find the monthly deposit A.
FV = €75 0004.2% / 12N = 156
A = €
27.NPV & IRRMedium
A project costs €400 000 today. It generates €29 000 at the end of the first year, and the cash flow then grows by 2% per year forever. The cost of capital is 10%. Find the net present value NPV.
C0 = €400 000CF1 = €29 000g = 2%r = 10%
NPV = €
28.NPV & IRRHard
A project costs €10 000 and returns €7 500, €4 000, €5 500 at the end of years 1 to 3. Find its internal rate of return IRR in percent.
CF0 = −10 000CF1 = 7 500CF2 = 4 000CF3 = 5 500
IRR = %
29.Normal distributionMedium
X is normal with mean 250 g and standard deviation 4 g; Y is independent and normal with mean 500 g and standard deviation 9 g. Find the probability P that X + Y is more than 735.2 g. Round z to two decimals.
μX = 250, σX = 4μY = 500, σY = 9X + Y > 735.2 g
P =
30.Normal distributionHard
The exam scores are normally distributed with mean 65 points and standard deviation 10 points. Find the value x below which 25% of them lie (the 25th percentile).
μ = 65σ = 10P = 25%
x = points
3Solutions
1.Beam support reactions
Answer: RB = 41.25 kN
Static equilibrium: a body at rest has ΣFx = 0, ΣFy = 0 and ΣM = 0 about any point. Taking moments about one support removes that support’s reaction from the equation.
Free-body diagram
Replace the supports with reaction forces. A roller gives one vertical reaction; with only vertical loads the pin gives one vertical reaction too.
Equilibrium
Moments about A, counter-clockwise positive. Then vertical force balance.
ΣMA = 0: RB·8 = 30·4 + 30·7
RB = 330/8 = 41.25 kN
ΣFy = 0: RA = 60 − 41.25 = 18.75 kN
Answer
RB = 41.25 kN
2.Beam support reactions
Answer: RB = 44.63 kN
Static equilibrium: a body at rest has ΣFx = 0, ΣFy = 0 and ΣM = 0 about any point. Taking moments about one support removes that support’s reaction from the equation.
Free-body diagram
Replace the supports with reaction forces. A roller gives one vertical reaction; with only vertical loads the pin gives one vertical reaction too.
Equilibrium
Moments about A, counter-clockwise positive. Then vertical force balance.
UDL resultant: W = w·ℓ = 6·7 = 42 kN, acting at x̄ = 3.5 m
ΣMA = 0: RB·8 = 30·7 + 6·7·3.5
RB = 357/8 = 44.63 kN
ΣFy = 0: RA = 72 − 44.63 = 27.38 kN
Answer
RB = 44.63 kN
3.Truss: method of joints
Answer: FBC = −21.21 kN
Method of joints: every pin of a truss is in equilibrium, so ΣFx = 0 and ΣFy = 0 at each joint. Start at a joint with at most two unknown member forces. A negative result means the member is in compression.
Support reactions
Equilibrium of the whole truss. Moments about A, counter-clockwise positive.
ΣMA = 0: Cy·6 = 45·2 → Cy = 15 kN
ΣFy = 0: Ay = 45 − 15 = 30 kN
Joint C
The roller gives only a vertical reaction, so joint C has just two unknown member forces. Assume both in tension.
sin γ = 4/5.657 = 0.7071, cos γ = 4/5.657 = 0.7071
Zero-force members: at a joint, the forces must balance in every direction. If a joint has no load and no support, and two of its three members are in line, then the balance across that line involves only the third member, so its force is zero. If only two members meet at an unloaded joint and they are not in line, both are zero.
Joint C
BC and CD lie on one straight line and the joint carries no load. Balancing forces across that line leaves only CG.
ΣF⊥ = 0 → FCG = 0
Check with the method of joints
Solving all joints gives these member forces (T = tension, C = compression).
Method of sections: cut the beam where you want the internal forces and keep one side. The shear force V balances the vertical forces on that side, and the bending moment M balances their moments about the cut.
Support reactions
Equilibrium of the whole beam.
UDL resultant: W = w·ℓ = 5·5 = 25 kN, acting at x̄ = 2.5 m
ΣMA = 0: RB·5 = 5·5·2.5
RB = 62.5/5 = 12.5 kN
ΣFy = 0: RA = 25 − 12.5 = 12.5 kN
Cut at x = 2.5 m and keep the left part
Only forces left of the cut appear. A UDL portion acts as one force at the middle of that portion.
M = RA·x − (each load left of the cut × its distance to the cut)
Mx = 12.5·2.5 − 5·2.5·1.25 = 15.63 kN·m
Answer
Mx = 15.63 kN·m
6.Shear force & bending moment
Answer: Mx = 15.5 kN·m
Method of sections: cut the beam where you want the internal forces and keep one side. The shear force V balances the vertical forces on that side, and the bending moment M balances their moments about the cut.
Support reactions
Equilibrium of the whole beam.
UDL resultant: W = w·ℓ = 3·8 = 24 kN, acting at x̄ = 4 m
ΣMA = 0: RB·8 = 40·1 + 3·8·4
RB = 136/8 = 17 kN
ΣFy = 0: RA = 64 − 17 = 47 kN
Cut at x = 7 m and keep the left part
Only forces left of the cut appear. A UDL portion acts as one force at the middle of that portion.
M = RA·x − (each load left of the cut × its distance to the cut)
Mx = 47·7 − 40·(7 − 1) − 3·7·3.5 = 15.5 kN·m
Answer
Mx = 15.5 kN·m
7.Beam deflection
Answer: δB = 10.55 mm
Beam deflection: deflections come from integrating the curvature M/(EI) twice. Standard cases are tabulated, for example PL³/(3EI) for a cantilever with an end load and 5wL⁴/(384EI) for a simply supported beam under a full uniform load. Deflections from separate loads add (superposition).
Consistent units
Work in N and mm: P in N, w in N/mm (kN/m = N/mm), L in mm, E in N/mm² (GPa × 1000), I in mm⁴.
P = 5000 N, a = 4500 mm, L = 5000 mm
EI = 1.68 × 10¹³ N·mm²
Deflection at the free end
Deflection at the load plus the slope at the load times the remaining length.
δP = Pa³/(3EI) = 9.04 mm
θP = Pa²/(2EI) = 0.003013 rad
δB = δP + θP·(L − a) = Pa²(3L − a)/(6EI) = 10.55 mm
Answer
δB = 10.55 mm
8.Beam deflection
Answer: Ireq = 100.4 × 10⁶ mm⁴
Beam deflection: deflections come from integrating the curvature M/(EI) twice. Standard cases are tabulated, for example PL³/(3EI) for a cantilever with an end load and 5wL⁴/(384EI) for a simply supported beam under a full uniform load. Deflections from separate loads add (superposition).
Deflection limit
δallow = L/360 = 5000/360 = 13.89 mm
Required second moment of area
I = 5wL⁴/(384·E·δallow) = 5·12·5000⁴/(384·70000·13.89) = 1.004 × 10⁸ mm⁴
Ireq = 100.4 × 10⁶ mm⁴
Answer
Ireq = 100.4 × 10⁶ mm⁴
9.Column buckling
Answer: Pcr = 5.895 kN
Euler buckling: a slender column becomes unstable at Pcr = π²EI/(KL)², where KL is the effective length: K = 1 for pinned–pinned, 2 for fixed–free, 0.7 for fixed–pinned and 0.5 for fixed–fixed. It buckles about the axis with the smallest I. The formula holds for slender columns, whose critical stress stays below the yield stress.
Smallest second moment of area
A rectangle buckles about its weak axis, using the smaller dimension cubed.
Euler buckling: a slender column becomes unstable at Pcr = π²EI/(KL)², where KL is the effective length: K = 1 for pinned–pinned, 2 for fixed–free, 0.7 for fixed–pinned and 0.5 for fixed–fixed. It buckles about the axis with the smallest I. The formula holds for slender columns, whose critical stress stays below the yield stress.
Smallest second moment of area
A round section has the same I about every axis.
I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴
Radius of gyration and slenderness
A = 1963 mm²
r = √(I/A) = √(306796/1963) = 12.5 mm
KL = 2·1500 = 3000 mm
KL/r = 3000/12.5 = 240
Critical stress
σcr = π²E/(KL/r)² = π²·200000/240² = 34.27 MPa
Answer
σcr = 34.27 MPa
11.Centroid & second moment of area
Answer: r = 63.74 mm
Second moment of area: for a rectangle about its own centroidal axis I = bh³/12, and about its base bh³/3 (parallel-axis theorem I = Ic + A·d²). A hole is handled by subtracting its I, as long as both parts share the same axis. The elastic section modulus is S = I/c, where c is the distance to the extreme fibre. The radius of gyration r = √(I/A) is the distance at which the whole area would give the same I; it is used in column buckling.
Net area and second moment of area
A = π(D² − d²)/4 = 5655 mm²
I = π(D⁴ − d⁴)/64 = 2.297 × 10⁷ mm⁴
Radius of gyration
r = √(I/A) = √(2.297 × 10⁷/5655) = 63.74 mm
Answer
r = 63.74 mm
12.Centroid & second moment of area
Answer: Ix = 114 × 10⁶ mm⁴
Composite areas: split the section into rectangles. The centroid is the area-weighted average ȳ = ΣAᵢyᵢ/ΣAᵢ. The second moment of area about the centroidal axis adds each part’s own bh³/12 and its parallel-axis term A·d², where d is the distance from the part’s centroid to the overall centroid. Bending stress is σ = M·y/I.
Split into three rectangles
Measure every y from the bottom edge.
top flange: A = 1400 mm², y = 311 mm
web: A = 3480 mm², y = 161 mm
bottom flange: A = 2720 mm², y = 8 mm
Centroid
ȳ = (1400·311 + 3480·161 + 2720·8)/7600 = 133.9 mm
Parallel-axis theorem
top flange: 140·10³/12 + 1400·177.1² = 4.393 × 10⁷ mm⁴
Effective stress (Terzaghi): the soil skeleton carries σ′ = σ − u. The total stress σ is the weight of everything above the point, Σγ·h, using the dry unit weight above the water table and the saturated unit weight below it. Below the water table the pore water pressure is hydrostatic: u = γw·(depth below the water table).
Pore water pressure
Water pressure grows with depth below the water table.
u = γw·(z − zw) = 9.81·6 = 58.86 kPa
Answer
u = 58.86 kPa
14.Effective stress in soil
Answer: σ′ = 106.1 kPa
Effective stress (Terzaghi): the soil skeleton carries σ′ = σ − u. The total stress σ is the weight of everything above the point, Σγ·h, using the dry unit weight above the water table and the saturated unit weight below it. Below the water table the pore water pressure is hydrostatic: u = γw·(depth below the water table).
Total stress
Add the weight of each soil band above X. Above the water table use the dry unit weight, below it the saturated one.
σ = γd·3 + γsat,sand·2 + γsat,clay·4
σ = 16·3 + 20.5·2 + 19·4 = 165 kPa
Pore water pressure
Water pressure grows with depth below the water table.
u = γw·(z − zw) = 9.81·6 = 58.86 kPa
Effective stress
σ′ = σ − u = 165 − 58.86 = 106.1 kPa
Answer
σ′ = 106.1 kPa
15.Consolidation settlement
Answer: S = 13.14 mm
Primary consolidation of a clay layer of thickness H, using the effective stress at the middle of the layer: σ′f = σ′0 + Δσ. Normally consolidated clay: S = Cc·H/(1 + e0)·log₁₀(σ′f/σ′0). Overconsolidated clay recompresses along Cr up to the preconsolidation pressure σ′p and follows Cc only beyond it. Time: t = Tv·Hdr²/cv, where the drainage path Hdr is H/2 when water can leave at the top and bottom, and H when only one side drains.
Preconsolidation pressure
σ′p = OCR·σ′0 = 2.5·105 = 262.5 kPa
Compare with the preconsolidation pressure
σ′f = 105 + 40 = 145 kPa ≤ 262.5 kPa
Settlement
The clay stays on its recompression line.
S = Cr·H/(1 + e0)·log₁₀(σ′f/σ′0) = 0.03·5/(1 + 0.6)·log₁₀(145/105) = 13.14 mm
Answer
S = 13.14 mm
16.Consolidation settlement
Answer: t = 3.229 years
Primary consolidation of a clay layer of thickness H, using the effective stress at the middle of the layer: σ′f = σ′0 + Δσ. Normally consolidated clay: S = Cc·H/(1 + e0)·log₁₀(σ′f/σ′0). Overconsolidated clay recompresses along Cr up to the preconsolidation pressure σ′p and follows Cc only beyond it. Time: t = Tv·Hdr²/cv, where the drainage path Hdr is H/2 when water can leave at the top and bottom, and H when only one side drains.
Drainage path
Hdr = H = 3 m
Time
t = Tv·Hdr²/cv = 0.287·3²/0.8 = 3.229 years
Answer
t = 3.229 years
17.Hydrostatic force
Answer: F = 88.29 kN
Hydrostatic force on a plane surface: F = ρ·g·yc·A. It acts at the centre of pressure, ycp = yc + Ic/(ycA), which always lies below the centroid.
Centroid depth and area
yc = d + H/2 = 2 + 2/2 = 3 m
A = H·b = 2·1.5 = 3 m²
Resultant force
F = ρ·g·yc·A = 1000·9.81·3·3 = 88.29 kN
Answer
F = 88.29 kN
18.Hydrostatic force
Answer: ycp = 3.111 m
Hydrostatic force on a plane surface: F = ρ·g·yc·A. It acts at the centre of pressure, ycp = yc + Ic/(ycA), which always lies below the centroid.
Centroid depth and area
yc = d + H/2 = 2 + 2/2 = 3 m
A = H·b = 2·1 = 2 m²
Resultant force
F = ρ·g·yc·A = 880·9.81·3·2 = 51.8 kN
Centre of pressure
Pressure grows with depth, so the resultant acts below the centroid.
Ic = b·H³/12 = 1·2³/12 = 0.6667 m⁴
ycp = yc + Ic/(yc·A) = 3 + 0.6667/(3·2) = 3.111 m
Answer
ycp = 3.111 m
19.Pipe friction & head loss
Answer: Δp = 229.2 kPa
Pipe flow: the Reynolds number Re = V·D/ν = ρ·V·D/μ tells laminar flow (Re below about 2000) from turbulent flow (above about 4000). Friction along a pipe costs the Darcy–Weisbach head loss hf = f·(L/D)·V²/(2g). Fittings, valves, the entrance and the exit add minor losses hm = ΣK·V²/(2g). Between two open reservoir surfaces the energy equation gives Δz = hf + hm, with the exit loss K = 1 included.
Head loss
hf = f·(L/D)·V²/(2g) = 0.022·(1000/0.3)·2.5²/(2·9.81) = 23.36 m
Pressure drop
The pipe is horizontal and the velocity does not change, so the energy equation reduces to Δp = ρ·g·hf.
Pipe flow: the Reynolds number Re = V·D/ν = ρ·V·D/μ tells laminar flow (Re below about 2000) from turbulent flow (above about 4000). Friction along a pipe costs the Darcy–Weisbach head loss hf = f·(L/D)·V²/(2g). Fittings, valves, the entrance and the exit add minor losses hm = ΣK·V²/(2g). Between two open reservoir surfaces the energy equation gives Δz = hf + hm, with the exit loss K = 1 included.
Loss coefficients
f·L/D = 0.018·300/0.2 = 27
ΣK = 0.5 + 2·0.3 + 2 + 1.0 = 4.1
Energy equation
Both surfaces are open to the air and practically at rest, so all of Δz is lost in the pipe.
Manning’s equation (SI): for uniform flow in an open channel V = (1/n)·R2/3·S1/2 and Q = A·V, with n the roughness coefficient, S the bed slope (m/m) and R = A/P the hydraulic radius: flow area divided by wetted perimeter. Only the boundary in contact with the water counts in P, not the free surface. Rectangle: A = b·y, P = b + 2y. Trapezoid with side slopes z:1 (horizontal:vertical): A = (b + z·y)·y, P = b + 2y·√(1 + z²). A circular pipe flowing full or half full has R = D/4.
Geometry
A = b·y = 4·0.8 = 3.2 m²
P = b + 2y = 4 + 2·0.8 = 5.6 m
R = A/P = 0.5714 m
Velocity
V = (1/0.022)·0.57142/3·0.0021/2 = 1.4 m/s
Flow rate
Q = A·V = 3.2·1.4 = 4.479 m³/s
Answer
Q = 4.479 m³/s
22.Open-channel flow (Manning)
Answer: S = 0.001994 m/m
Manning’s equation (SI): for uniform flow in an open channel V = (1/n)·R2/3·S1/2 and Q = A·V, with n the roughness coefficient, S the bed slope (m/m) and R = A/P the hydraulic radius: flow area divided by wetted perimeter. Only the boundary in contact with the water counts in P, not the free surface. Rectangle: A = b·y, P = b + 2y. Trapezoid with side slopes z:1 (horizontal:vertical): A = (b + z·y)·y, P = b + 2y·√(1 + z²). A circular pipe flowing full or half full has R = D/4.
Geometry
A = (b + z·y)·y = (5 + 2·2.3)·2.3 = 22.08 m²
P = b + 2y·√(1 + z²) = 5 + 2·2.3·√(1 + 2²) = 15.29 m
Pump power: the power given to the water is the hydraulic power Ph = ρ·g·Q·H = Q·Δp. The shaft must deliver more because the pump loses some: Pshaft = Ph/ηpump, and the motor draws Pel = Pshaft/ηmotor. The pump head H follows from the energy equation: static lift + pressure head + all friction and minor losses.
Hydraulic power
Ph = ρ·g·Q·H = 1000·9.81·0.04·12 = 4.709 kW
Efficiency
η = Ph/Pshaft = 4.709/5.7 = 82.61%
Answer
η = 82.61%
24.Pumps & pump power
Answer: Pshaft = 2.978 kW
Pump power: the power given to the water is the hydraulic power Ph = ρ·g·Q·H = Q·Δp. The shaft must deliver more because the pump loses some: Pshaft = Ph/ηpump, and the motor draws Pel = Pshaft/ηmotor. The pump head H follows from the energy equation: static lift + pressure head + all friction and minor losses.
Time value of money: money today is worth more than the same amount later, because it can earn interest. With interest rate i per period, an amount grows as FV = PV·(1 + i)n. A series of equal payments (an annuity) is valued by adding up each payment, which gives the annuity formulas.
Exponent
r·t = 0.07·12 = 0.84
e0.84 = 2.316
Future value
FV = 20 000·2.316 = €46 327.34
Answer
FV = €46 327.34
26.Time value of money
Answer: A = 362.2 €
Time value of money: money today is worth more than the same amount later, because it can earn interest. With interest rate i per period, an amount grows as FV = PV·(1 + i)n. A series of equal payments (an annuity) is valued by adding up each payment, which gives the annuity formulas.
Net present value: discount every cash flow back to today with the required rate of return r, then add them up: NPV = Σ CFₜ/(1 + r)ᵗ. A project with NPV > 0 earns more than r and adds value. The internal rate of return (IRR) is the rate at which NPV = 0.
Value of the cash flows
PV = 29 000/(0.1 − 0.02) = 362 500.00
NPV
NPV = 362 500.00 − 400 000 = −€37 500.00
Answer
NPV = −€37 500.00
NPV is negative, so the project earns less than the required return.
28.NPV & IRR
Answer: IRR = 34.88 %
Net present value: discount every cash flow back to today with the required rate of return r, then add them up: NPV = Σ CFₜ/(1 + r)ᵗ. A project with NPV > 0 earns more than r and adds value. The internal rate of return (IRR) is the rate at which NPV = 0.
Bracket the IRR
NPV falls as the rate rises, so the IRR lies between a rate with positive NPV and one with negative NPV.
NPV(30%) = 639.51
NPV(35%) = −14.23
Narrow down
Linear interpolation gives a first estimate; repeating it (or a calculator's IRR function) gives the exact rate.
IRR ≈ 30 + 5·639.51/(639.51 − (−14.23)) = 34.89%
exact: NPV(34.88%) = 0 → IRR = 34.88%
Answer
IRR = 34.88%
29.Normal distribution
Answer: P = 0.9332
Sums and differences of normal variables: if X and Y are independent and normal, X + Y and X − Y are normal too. Means add or subtract: E(X ± Y) = μX ± μY. Variances always add, also for a difference: Var(X ± Y) = σX² + σY², so σ = √(σX² + σY²). Standard deviations never add directly. The total of n independent items with the same distribution has mean n·μ and standard deviation σ·√n. Then standardise as usual: z = (x − mean)/σ and read Φ(z).
Distribution of the sum
μ = 750 g, σ = √(4² + 9²) = 9.849 g
Standardise
z = (735.2 − 750)/9.849 = −1.5
Upper tail
P = 1 − Φ(−1.5) = 1 − 0.06681 = 0.9332
Answer
P = 0.9332
30.Normal distribution
Answer: x = 58.26 points
Normal distribution: convert a value x to a z-score, z = (x − μ)/σ, which counts standard deviations from the mean. A z-table (or Φ(z)) gives the area to the left of z, which is P(X < x). The area to the right is 1 − Φ(z), and the area between two values is the difference of their two areas.