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 = 5 mP1 = 40 kN at 1 mP2 = 50 kN at 2 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 = 6 mP = 25 kN at 3.5 mw = 3 kN/m from 0 to 5 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 = 2 mB is 5 m from Avertical load = 25 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: 15 kN ↓H: 10 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 = 1.5 m from support A. Use sagging moments positive and upward shear on the left face positive.
L = 8 mw = 6 kN/m from 0 to 8 mx = 1.5 m
Mx = kN·m
6.Shear force & bending momentHard
The simply supported beam carries a uniformly distributed load w over its whole span and a point load P. Find the maximum bending moment Mmax in the beam.
L = 8 mP = 60 kN at 6 mw = 2 kN/m from 0 to 8 m
Mmax = kN·m
7.Beam deflectionMedium
A simply supported beam carries a uniformly distributed load w over its whole span. Find the maximum deflection δmax, at midspan.
w = 10 kN/mL = 3 mE = 11 GPa (timber)I = 120 × 10⁶ mm⁴
δmax = mm
8.Beam deflectionHard
A propped cantilever is fixed at A and simply supported (propped) at B. It carries a uniformly distributed load w over its whole length. EI is constant. Find the reaction RB at the prop.
L = 6 mw = 5 kN/m
RB = kN
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)d = 50 mm
Pcr = kN
10.Column bucklingHard
A column of length L made of steel is pinned at both ends. Find its slenderness and the critical buckling stress σcr from Euler’s formula.
L = 5 mE = 200 GPa (steel)d = 60 mm
σcr = MPa
11.Centroid & second moment of areaMedium
A channel section lies with its open side up: a base plate 160 mm wide with two upright legs, 170 mm high overall. Base and legs are all 14 mm thick. Find its second moment of area Ix about the horizontal axis x through its centroid.
b = 160 mmh = 170 mmt = 14 mm
Ix = × 10⁶ mm⁴
12.Centroid & second moment of areaHard
A 100 mm × 150 mm plate has a rectangular cut-out 40 mm wide and 40 mm high, centred on the vertical centre line with its centre 97.5 mm above the bottom. Find the second moment of area Ix about the horizontal axis through the centroid.
b = 100 mmh = 150 mmcut-out: 40 × 40 mmhole centre: y = 97.5 mm
Ix = × 10⁶ mm⁴
13.Effective stress in soilMedium
For the soil profile shown, find the pore water pressure u at point X.
A 7 m thick overconsolidated clay layer has Cc = 0.5, Cr = 0.02 and e₀ = 0.7. At its mid-depth the effective vertical stress is 120 kPa and the overconsolidation ratio is OCR = 1.5. A new building adds 20 kPa of vertical stress at mid-depth. Find the primary consolidation settlement S in mm.
In a laboratory oedometer test, a clay specimen 25 mm thick, drained at the top and bottom, reaches 50% consolidation in 15 min. In the field the same clay forms a layer 4.5 m thick. It lies between two sand layers, so water drains out at the top and the bottom. How many days t does the field layer take to reach 50% consolidation?
A vertical rectangular gate 3 m wide sits in the wall of a tank of water, with its top edge 2.5 m below the surface. Find the resultant hydrostatic force F on the gate.
A vertical rectangular gate 2 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 = 3 mb = 2 mρ = 880 kg/m³g = 9.81 m/s²
ycp = m
19.Pipe friction & head lossMedium
Water flows at 25 L/s through a pipe 800 m long with an inside diameter of 150 mm. Use the Darcy friction factor f = 0.022 and g = 9.81 m/s². Find the friction head loss hf.
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.02. The minor losses are a square-edged entrance (K = 0.5), 2 bends (K = 0.3 each), a valve (K = 10) and the exit into the lower reservoir (K = 1.0). Use g = 9.81 m/s². What difference in water level Δz between the reservoirs gives a flow rate of 40 L/s?
A rectangular gravel channel is 3.5 m wide. Water flows uniformly 1.5 m deep. Manning n = 0.025 and the bed slope is S = 0.003. Find the flow rate Q in m³/s.
b = 3.5 my = 1.5 mn = 0.025S = 0.003
Q = m³/s
22.Open-channel flow (Manning)Hard
A trapezoidal gravel channel with a bottom width of 6 m and side slopes of 1:1 (horizontal:vertical) must carry 21 m³/s in uniform flow at a depth of 1.6 m. Manning n = 0.025. Find the required bed slope S (in m/m).
Q = 21 m³/sb = 6 my = 1.6 mz = 1n = 0.025
S = m/m
23.Pumps & pump powerMedium
A pump with an efficiency of 85% delivers 20 L/s of water (ρ = 1000 kg/m³) against a head of 78 m. Use g = 9.81 m/s². Find the shaft power Pshaft needed, in kW.
Q = 20 L/sH = 78 mη = 85%g = 9.81 m/s²
Pshaft = kW
24.Pumps & pump powerHard
A pump moves 60 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 300 m long with an inside diameter of 200 mm and a Darcy friction factor f = 0.022; the minor-loss coefficients add up to ΣK = 4 (entrance and exit included). Use g = 9.81 m/s². The pump efficiency is 75%. Find the shaft power Pshaft in kW.
You invest €8 000 at a nominal rate of 2% per year, compounded monthly. Find the value FV of the investment after 10 years.
PV = €8 0002% / 12n = 10
FV = €
26.Time value of moneyHard
You borrow €15 000 at 4.2% per year, compounded monthly, and repay €350 at the end of every month. How many months does it take until the loan is repaid (the last payment may be smaller)? Give the number of payments N.
P = €15 000A = €3504.2% / 12
N = payments
27.NPV & IRRMedium
A machine costs €50 000 today. It generates €10 000 at the end of each year for 5 years and can then be sold for €5 000. The cost of capital is 10%. Find the net present value NPV.
C0 = €50 000CF = €10 000S = €5 000n = 5r = 10%
NPV = €
28.NPV & IRRHard
A project costs €40 000 and returns €12 500, €17 000, €22 000, €12 500 at the end of years 1 to 4. Cash flows can be reinvested at the cost of capital of 8%. Find the modified internal rate of return MIRR in percent.
The masses of bags of cement are independent and normal with mean 25 kg and standard deviation 0.5 kg. Find the probability P that the total mass of 6 of them is below 148.2 kg. Round z to two decimals.
μ = 25 kgσ = 0.5 kgn = 6< 148.2 kg
P =
30.Normal distributionHard
The heights of adult men are normally distributed with mean 178 cm and standard deviation 7 cm. Find the value x below which 25% of them lie (the 25th percentile).
μ = 178σ = 7P = 25%
x = cm
3Solutions
1.Beam support reactions
Answer: RB = 28 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·5 = 40·1 + 50·2
RB = 140/5 = 28 kN
ΣFy = 0: RA = 90 − 28 = 62 kN
Answer
RB = 28 kN
2.Beam support reactions
Answer: RB = 20.83 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·ℓ = 3·5 = 15 kN, acting at x̄ = 2.5 m
ΣMA = 0: RB·6 = 25·3.5 + 3·5·2.5
RB = 125/6 = 20.83 kN
ΣFy = 0: RA = 40 − 20.83 = 19.17 kN
Answer
RB = 20.83 kN
3.Truss: method of joints
Answer: FBC = −23.29 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 = 25·5 → Cy = 20.83 kN
ΣFy = 0: Ay = 25 − 20.83 = 4.167 kN
Joint C
The roller gives only a vertical reaction, so joint C has just two unknown member forces. Assume both in tension.
sin γ = 2/2.236 = 0.8944, cos γ = 1/2.236 = 0.4472
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·ℓ = 6·8 = 48 kN, acting at x̄ = 4 m
ΣMA = 0: RB·8 = 6·8·4
RB = 192/8 = 24 kN
ΣFy = 0: RA = 48 − 24 = 24 kN
Cut at x = 1.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 = 24·1.5 − 6·1.5·0.75 = 29.25 kN·m
Answer
Mx = 29.25 kN·m
6.Shear force & bending moment
Answer: Mmax = 102 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. Maximum moment: since dM/dx = V, the bending moment has its peak where the shear force passes through zero, either at a point load or inside a distributed load.
Support reactions
Equilibrium of the whole beam.
UDL resultant: W = w·ℓ = 2·8 = 16 kN, acting at x̄ = 4 m
ΣMA = 0: RB·8 = 60·6 + 2·8·4
RB = 424/8 = 53 kN
ΣFy = 0: RA = 76 − 53 = 23 kN
Where is V = 0?
V(6⁻) = RA − w·6 = 11 kN
V(6⁺) = 11 − 60 = −49 kN
The shear changes sign at the load, so the maximum is under the load.
Maximum moment
Mmax = RA·6 − w·6²/2 = 23·6 − 2·6²/2 = 102 kN·m
Answer
Mmax = 102 kN·m
It occurs at x = 6 m from A.
7.Beam deflection
Answer: δmax = 7.99 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⁴.
w = 10 N/mm, L = 3000 mm
EI = 1.32 × 10¹² N·mm²
Deflection
δmax = 5wL⁴/(384EI) = 5·10·3000⁴/(384·1.32 × 10¹²) = 7.99 mm
Answer
δmax = 7.99 mm
8.Beam deflection
Answer: RB = 11.25 kN
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). For a statically indeterminate beam, remove a redundant support and require the deflection there to be zero (compatibility).
Remove the prop
Without the prop the beam is a cantilever. Its free end would deflect under the load.
δload = wL⁴/(8EI)
Compatibility at the prop
The prop force RB alone lifts the free end by RBL³/(3EI). The prop holds B at zero deflection, so the two must cancel. EI cancels out.
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.
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 = π·60⁴/64 = 636173 mm⁴
Radius of gyration and slenderness
A = 2827 mm²
r = √(I/A) = √(636173/2827) = 15 mm
KL = 1·5000 = 5000 mm
KL/r = 5000/15 = 333.3
Critical stress
σcr = π²E/(KL/r)² = π²·200000/333.3² = 17.77 MPa
Answer
σcr = 17.77 MPa
11.Centroid & second moment of area
Answer: Ix = 19.59 × 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 rectangles
Take the base at full width and the two legs above it (height h − t), so nothing is counted twice.
base: A1 = 160·14 = 2240 mm², y1 = 7 mm
two legs: A2 = 2·14·156 = 4368 mm², y2 = 14 + 156/2 = 92 mm
Sections with holes: treat a hole as a negative area. The centroid is ȳ = (A₁y₁ − A₂y₂)/(A₁ − A₂), with A₁, y₁ for the full plate and A₂, y₂ for the hole. For the second moment of area about the centroidal axis, use the parallel-axis theorem for each part and subtract the hole: I = (Ic1 + A₁d₁²) − (Ic2 + A₂d₂²), where d is the distance from each part's own centroid to the axis through ȳ. A rectangle has Ic = bh³/12 and a circle Ic = πd⁴/64.
Centroid
ȳ = (15000·75 − 1600·97.5)/(15000 − 1600) = 72.31 mm
Full plate about ȳ
100·150³/12 + 15000·2.687² = 28.23 × 10⁶ mm⁴
Hole about ȳ
40·40³/12 + 1600·25.19² = 1.228 × 10⁶ mm⁴
Subtract
Ix = 27 × 10⁶ mm⁴
Answer
Ix = 27 × 10⁶ mm⁴
13.Effective stress in soil
Answer: u = 9.81 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).
Pore water pressure
Water pressure grows with depth below the water table.
u = γw·(z − zw) = 9.81·1 = 9.81 kPa
Answer
u = 9.81 kPa
14.Effective stress in soil
Answer: u = 49.05 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).
Pore water pressure
Water pressure grows with depth below the water table.
u = γw·(z − zw) = 9.81·5 = 49.05 kPa
Answer
u = 49.05 kPa
15.Consolidation settlement
Answer: S = 5.513 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 = 1.5·120 = 180 kPa
Compare with the preconsolidation pressure
σ′f = 120 + 20 = 140 kPa ≤ 180 kPa
Settlement
The clay stays on its recompression line.
S = Cr·H/(1 + e0)·log₁₀(σ′f/σ′0) = 0.02·7/(1 + 0.7)·log₁₀(140/120) = 5.513 mm
Answer
S = 5.513 mm
16.Consolidation settlement
Answer: t = 337.5 days
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 paths
Hdr,lab = 25/2 = 12.5 mm = 0.0125 m
Hdr,field = 4.5/2 = 2.25 m
Scale the time
t = Tv·Hdr²/cv with the same T_v and cv in the lab and the field.
tfield = 15·(2.25/0.0125)² = 486000 min
tfield = 486000/1440 = 337.5 days
Answer
t = 337.5 days
17.Hydrostatic force
Answer: F = 275.9 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.5 + 2.5/2 = 3.75 m
A = H·b = 2.5·3 = 7.5 m²
Resultant force
F = ρ·g·yc·A = 1000·9.81·3.75·7.5 = 275.9 kN
Answer
F = 275.9 kN
18.Hydrostatic force
Answer: ycp = 3.714 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 + 3/2 = 3.5 m
A = H·b = 3·2 = 6 m²
Resultant force
F = ρ·g·yc·A = 880·9.81·3.5·6 = 181.3 kN
Centre of pressure
Pressure grows with depth, so the resultant acts below the centroid.
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.
Mean velocity
A = π·0.15²/4 = 0.01767 m²
V = Q/A = 0.025/0.01767 = 1.415 m/s
Velocity head
V²/(2g) = 1.415²/(2·9.81) = 0.102 m
Darcy–Weisbach
hf = f·(L/D)·V²/(2g) = 0.022·(800/0.15)·1.415²/(2·9.81) = 11.97 m
Answer
hf = 11.97 m
20.Pipe friction & head loss
Answer: Δz = 3.479 m
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.
Velocity
V = Q/A = 0.04/0.03142 = 1.273 m/s
V²/(2g) = 0.08263 m
Loss coefficients
f·L/D = 0.02·300/0.2 = 30
ΣK = 0.5 + 2·0.3 + 10 + 1.0 = 12.1
Total head loss
hf = 30·0.08263 = 2.479 m
hm = 12.1·0.08263 = 0.9998 m
Δz = hf + hm = 3.479 m
Answer
Δz = 3.479 m
21.Open-channel flow (Manning)
Answer: Q = 9.976 m³/s
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 = 3.5·1.5 = 5.25 m²
P = b + 2y = 3.5 + 2·1.5 = 6.5 m
R = A/P = 0.8077 m
Velocity
V = (1/0.025)·0.80772/3·0.0031/2 = 1.9 m/s
Flow rate
Q = A·V = 5.25·1.9 = 9.976 m³/s
Answer
Q = 9.976 m³/s
22.Open-channel flow (Manning)
Answer: S = 0.001538 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 = (6 + 1·1.6)·1.6 = 12.16 m²
P = b + 2y·√(1 + z²) = 6 + 2·1.6·√(1 + 1²) = 10.53 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.02·78 = 15.3 kW
Shaft power
Pshaft = Ph/η = 15.3/0.85 = 18 kW
Answer
Pshaft = 18 kW
24.Pumps & pump power
Answer: Pshaft = 10.89 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.
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.
Rearrange the annuity formula
i = 0.0035
(1 + i)−N = 1 − 15 000·0.0035/350 = 0.85
Solve with logarithms
N = −ln(0.85)/ln(1.004) = 46.52
Round up to whole payments
46 full payments leave a small balance, which the last, smaller payment clears.
N = 47
Answer
N = 47 payments
27.NPV & IRR
Answer: NPV = −8988 €
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.
Present value of the yearly cash flows
10 000·3.791 = 37 907.87
Present value of the salvage value
5 000/1.611 = 3 104.61
NPV
NPV = −50 000 + 37 907.87 + 3 104.61 = −€8 987.53
Answer
NPV = −€8 987.53
28.NPV & IRR
Answer: MIRR = 15.76 %
Modified IRR: compound every inflow forward to the end of the project at the reinvestment rate (usually the cost of capital) to get the terminal value TV. MIRR is the rate that turns the initial investment into TV: MIRR = (TV/C₀)^(1/n) − 1.
Terminal value of the inflows
12 500·1.083 = 15 746.40
17 000·1.082 = 19 828.80
22 000·1.081 = 23 760.00
12 500·1.080 = 12 500.00
TV = 71 835.20
MIRR
MIRR = (71 835.20/40 000)1/4 − 1 = 15.76%
Answer
MIRR = 15.76%
29.Normal distribution
Answer: P = 0.07078
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 total
μT = 150 kg, σT = 0.5·√6 = 1.225 kg
Standardise
z = (148.2 − 150)/1.225 = −1.47
Probability
P = Φ(−1.47) = 0.07078
Answer
P = 0.07078
30.Normal distribution
Answer: x = 173.3 cm
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.