Update
Update
c
ADDIS ABABA SCIENCE AND
TECHNOLOGY UNIVERSITY
Department of Electromechanical
Industrial Automation and process control Project
Title – Water Treatment
Section – D
Name Id.no
1 Sena Mekonin ETS1466/14
2 Sem Wogayehu ETS 1462/14
3 Sinishawu Nigatu ETS 1481/14
4 Yonathan kefyalew ETS1742/14
Introduction
Water treatment is a crucial process for supplying safe and clean water for drinking,
industrial use, and environmental protection. The process involves several steps to
remove impurities, neutralize harmful chemicals, and improve water quality. Water
treatment facilities employ automation and control systems to ensure maximum efficiency,
reduce human intervention, and comply with regulatory standards.
Term Description
Setpoint Desired value for water quality parameters
Actuating Error Difference between setpoint and measured value
Disturbance External factors affecting the system (e.g., chlorine
addition)
To derive the transfer function for water treatment system, we will analyse each mass balance
principle.
The flow dynamics of each tank can be modelled as first-order
Step 1: Define the Transfer Function for a Single Tank
𝑄𝑖𝑛(𝑠)
H(s) =
1
𝐴𝑠 + 𝑅
1
𝐻12 (𝑠) =
(𝜏1 𝑠 + 1) (𝜏2 𝑠 + 1)
1
𝐻12 (𝑠) =
(5𝑠 + 1) (7𝑠 + 1)
1
𝐻12 (𝑠) = 35𝑠2 +12𝑠+ 1
Tank 2 & Tank 3 (Interactive with Chlorine Disturbance)
For interactive tanks, the output of Tank 2 is the input to Tank 3:
1
𝐻3 (𝑠) =
𝜏3 𝑠 + 1
1 𝐷(𝑠)
𝐻23 (𝑠) = + , 𝑏𝑢𝑡 𝑄𝑖𝑛2 (𝑠) = (7𝑠 + 1)
(𝜏2 𝑠 + 1)(𝜏3 𝑠 + 1) 𝑄𝑖𝑛2 (𝑠)
1 𝐷(𝑠)
𝐻23 (𝑠) = +
(7𝑠 + 1)(6𝑠 + 1) (7𝑠 + 1)
1 𝐷(𝑠)
𝐻23 (𝑠) = +
42𝑠 2 + 13𝑠 + 1 (7𝑠 + 1)
D(s) – Disturbance due to chlorine.
Tank 3 & Tank 4 (Non-Interactive)
1
𝐻4 (𝑠) =
𝜏4 𝑠 + 1
1
( )
𝐻34 𝑠 =
(𝜏3 𝑠 + 1) (𝜏4 𝑠 + 1)
1
𝐻34 (𝑠) =
(6𝑠 + 1) (4𝑠 + 1)
1
𝐻34 (𝑠) =
24𝑠 2 + 10𝑠 + 1
1
𝐻5 (𝑠) =
𝜏5 𝑠 + 1
𝐾𝑝
𝐻45 (𝑠) =
(𝜏4 𝑠 + 1) (𝜏5 𝑠 + 1)
2
𝐻45 (𝑠) =
(4𝑠 + 1) (3𝑠 + 1)
2
=
12𝑠 2 + 7𝑠+ 1
1 1 2
= ∗ ∗
(5𝑠 + 1) (7𝑠 + 1) (6𝑠 + 1)(4𝑠 + 1) (4𝑠 + 1) (3𝑠 + 1)
Transfer Function of with disturbance (Tank 3→4, 4→5)
𝐷(𝑠) 1 2
= ∗ ∗
7𝑠 + 1 4𝑠 + 1 3𝑠 + 1
𝐾𝑝𝑠
𝐻𝑡𝑜𝑡𝑎𝑙 (𝑠) =
(𝜏1 𝑠 + 1)(𝜏2 𝑠 + 1)(𝜏3 𝑠 + 1) (𝜏4 𝑠 + 1) (𝜏5 𝑠 + 1)
𝐾𝑝 ∗ 𝐷(𝑠)
+
(𝜏2 𝑠 + 1)(𝜏4 𝑠 + 1) (𝜏5 𝑠 + 1)
2 2𝐷(𝑠)
𝐻𝑡𝑜𝑡𝑎𝑙 (𝑠) = +
(5𝑠+ 1)(7𝑠+ 1) (6𝑠+ 1) (4𝑠+ 1)(3𝑠+ 1) (7𝑠+ 1)(4𝑠+1)(3𝑠+1)