Next Article in Journal
A Probabilistic Conductor Size Selection Framework for Active Distribution Networks
Previous Article in Journal
Progress in Development of Nanostructured Manganese Oxide as Catalyst for Oxygen Reduction and Evolution Reaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Feedback Control Strategies in a Plant-Wide Wastewater Treatment Plant for Simultaneous Evaluation of Economics, Energy Usage, and Removal of Nutrients

by
Abdul Gaffar Sheik
1,
Eagalapati Tejaswini
1,
Murali Mohan Seepana
1,
Seshagiri Rao Ambati
1,*,
Montse Meneses
2 and
Ramon Vilanova
2,*
1
Department of Chemical Engineering, National Institute of Technology, Warangal 506 004, Telangana, India
2
Department of Telecommunication and System Engineering, Universitat Autonoma de Barcelona, 08193 Barcelona, Spain
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(19), 6386; https://doi.org/10.3390/en14196386
Submission received: 16 August 2021 / Revised: 20 September 2021 / Accepted: 28 September 2021 / Published: 6 October 2021

Abstract

:
Simultaneous removal of nitrogen and phosphorous is a recommended practice while treating wastewater. In the present study, control strategies based on proportional-integral (PI), model predictive control (MPC), and fuzzy logic are developed and implemented on a plant-wide wastewater treatment plant. Four combinations of control frameworks are developed in order to reduce the operational cost and improve the effluent quality. As a working platform, a Benchmark simulation model (BSM2-P) is used. A default control framework with PI controllers is used to control nitrate and dissolved oxygen (DO) by manipulating the internal recycle and oxygen mass transfer coefficient (KLa). Hierarchical control topology is proposed in which a lower-level control framework with PI controllers is implemented to DO in the sixth reactor by regulating the KLa of the fifth, sixth, and seventh reactors, and fuzzy and MPC are used at the supervisory level. This supervisory level considers the ammonia in the last aerobic reactor as a feedback signal to alter the DO set-points. PI-fuzzy showed improved effluent quality by 21.1%, total phosphorus removal rate by 33.3% with an increase of operational cost, and a slight increase in the production rates of greenhouse gases. In all the control design frameworks, a trade-off is observed between operational cost and effluent quality.

1. Introduction

The increase in the global population and urbanization has emanated in an increase in the use of water and therefore in the production of contaminated water (wastewater). This causes the effect of the water cycle and generates disorder in natural functioning. The intensification of the demand for clean water with deficient water resources has eventually resulted in a growing interest towards resource recovery during the treatment of wastewater. It is perceived that valuable resources like clean water, energy, and nutrients can be recovered from wastewater [1]. This leads to the progression of wastewater treatment plants (WWTP) into a water resource recovery facility (WRRF) [2]. Modelling and control of WWTPs generated a lot of interest among the wastewater community to increase flexibility and reduce the costs while operating WWTPs. WWTP’s are intricate because of chemical and biological interactions in between the processes, the peculiar nature of microbes, a disorder of concentrations, and dynamic rates of flow [3,4]. Among all the benchmark simulation models, BSM1-P is a platform with a defined plant layout, bioprocess models (activated sludge models (ASM)), influent loads, sensors, and actuators. Hence, BSM1-P platform allows unbiased comparison of carbon, nitrogen, and phosphorus removal in activated sludge control strategies [5,6]. Based on the models of WWTPs, plant-wide models take the attention of researchers for a long time run as the whole plant is understood by considering all process interactions [7,8,9]. Different plant-wide models are studied in the literature, which includes sludge control approaches, biogas production in primary settler, handling of the anaerobic digester, and phosphorus modeling with interactions of sulfur and iron cycles [10,11,12,13,14,15,16,17]. In the literature, many control applications like fuzzy logic controller (FLC), model predictive controller (MPC), proportional-integral (PI), and ammonia-based aeration control (ABAC) with different hierarchical combinations of PI, MPC, and fuzzy were studied, and it is observed that there is a trade-off between operational cost and effluent quality [18,19,20,21]. Maheswari et al. (2020) designed the nested control loop on three-stage biological treatment for ammonia changes and they observed that Effluent Quality Index (EQI) is improved with higher operational costs [22]. Shiek et al. (2021) implemented an ammonia-based aeration control (ABAC) with four different combinations of controllers like PI-MPC, MPC-MPC, PI-fuzzy, and MPC-fuzzy, which resulted in a tradeoff between Operational Cost Index (OCI) and EQI. In their study, the ammonia removal rate was improved by 18% in the case of MPC-MPC but P removal was not affected much [23].
Based on the benchmark simulation model on the plant-wide model (BSM2), different control applications like PI, artificial neural network (ANN), and sludge-based strategies, hierarchal control approaches are reported [24,25]. It was reported that the application of FLC improves EQI by 8% and reduces the aeration energy consumption up to 13%. Revollar et al. (2020) implemented a cascade strategy by using PI control with DO and ammonia on the BSM2 platform, and they noticed a 9% improvement on operational cost [26]. Barbu et al. (2017) implemented controllers for DO, ammonia, and TSS (total suspended solids) on the BSM2-G (greenhouse gas emission) platform, and as a result, an improved EQI of 11% and reduced OCI of 7.7% were achieved [27]. Supervisory control approaches on BSM2 are developed with different combinations of PI, MPC, and fuzzy, and it was found that PI-MPC showed a reduced OCI of 6.7%. Similarly, PI-fuzzy provided an improved EQI of 9.2% [28]. Even though many researchers designed control approaches based on BSM2, controllers design based on BSM2-P is limited.
The benchmark simulation model (BSM2-P) is an integrated version of BSM1-P that includes both water and sludge treatment process units [29] and BSM2. In the BSM2-P framework, PI-based controlling DO in the sixth reactor and cascade PI in the control of ammonia and total suspended solids were developed and it was found that OCI and EQI are in the trade-off. Closed-loop control showed an improved result when compared with open-loop operation [30]. Sludge management strategies like bio-solids beneficiation facility (BBF) are studied to improve solubility, sludge dewaterability, and handle high sludge loads with change in the microbial population [31].
To the best of the authors’ knowledge, advanced hierarchical control strategies are not implemented on the BSM2-P plant level layout to improve effluent quality and operational costs. On this note, the present work focuses on implementing advanced hierarchical control strategies based on MPC, fuzzy controllers; and the plant performance is evaluated by developing different control approaches. A lower-level control strategy is developed for DO control in the sixth reactor by manipulating the KLa of the fifth, sixth, and seventh reactors in the biological treatment process. PI controller is used at the lower level, whereas fuzzy and MPC are used at the supervisory level. The supervisory level is based on the ABAC to determine the dissolved oxygen set-point. EQI; OCI; percentage of nutrient removals; and production rates of methane, carbon dioxide, and hydrogen are studied while evaluating the control strategies.

2. Materials and Methods

This section provides an explanation of the working regime of BSM2-P with details of the performance/economy assessment and indices. Model creation and simulation are performed using MATLAB/SIMULINK.

2.1. Influent Data

Model-based influent load generation as elucidated in [2,32,33] is used to generate dynamic influent load data to execute the performance of plant-wide scenarios of the wastewater treatment plants. The daily average dynamic mass flow rates are provided in Table 1. More information about the handling of influent generation is illustrated in [15,34,35]. Dynamic simulations are performed for 609 days with steady-state simulation for 300 days. S:COD is the ratio of added sulphate [2]. The last 1-year data are used for performance assessment of the plant. State variables of ASM2d, units with notations, and average influent data are reported in the Supplementary data Table S1.

2.2. Model Scenario

The plant-wide model of BSM2-P has a resemblance to BSM2 plant [9] but a modification is carried out in the activated sludge unit (ASU). In ASU, two anaerobic reactors are added ahead of anoxic and aerobic reactors (A2/O) to enhance phosphorus removal and to improve PAO’s with a competitive dominance over other nitrogenous bacteria. The plant-wide model of BSM2-P consists of ASU, primary (PSU) and secondary (SSU) sedimentation units, Thickener (THK), anaerobic digestion (ADU) unit, storage (SU), and dewatering (DU) unit, with internal and external recycles. Figure 1 depicts the plant-wide model of BSM2-P and Table 2 represents each process unit of WWTP of BSM2-P with their working function and physical configurations. The internal recycle is maintained at three times the influent data, and sludge flow is at the rate of 600 m3/d.

2.3. Plant-Wide Assessment Criteria

2.3.1. Effluent Quality Index

The Effluent Quality Index (EQI) defines the amount of effluent to surface waters averaged over the assessment time interval related to the weighting factors of discharge loads of composition. The EQI (kg pollutants units/D) [2,7] is integrated with additional phosphorus process reactions, and the pollution composition with individual variables is illustrated below:
EQI = 1 T · 1000 t start t end ( θ TSS TSS ef t + θ COD COD ef t + θ NKJ NKJ ef t + θ NO S NO ef t + θ BOD 5 BOD 5 ef t + θ P org P org ef t + θ P inorg P inorg ef t )   Q ef t   dt  
COD ef = S F ef + S A ef + S I ef + X I ef + X S ef + X B , H ef + X PAO ef + X PHA ef + X B , A ef + i COD S Fe II S Fe II ef + i COD S IS S IS ef + i COD X S O X S ef O + X SRB ef  
NKJ ef = S NH ef + i N S F S F ef + i N S I S I ef + i N X I X I ef + i N X S X S ef + i N BM X B , H ef + X PAO ef + X B , A ef + X SRB ef
P org ef = X PP ef + i P S F S F ef + i P S I S I ef + i P X I X I ef + i P X S X S ef + i P BM X B , H ef + X PAO ef + X B , A ef + X SRB ef
BOD 5 ef = 0.25 S F ef + S A ef + 1 f S I   X S ef + 1 f X IH X B , H ef + 1 f X IP X PAO ef + X PHA ef + 1 f X IA X B , A ef + X SRB ef
P inorg ef = S PO 4 ef
S NO ef = S NO 3
TSS ef = X TSS
The subscript ‘ef’ indicates the effluent discharge. θ i signifies the weighting factors of different pollutants to convert into basic pollution units, which are tabulated in Table 3. i COD i denotes the COD compounds, i N i denotes the nitrogen compounds, i P i denotes the phosphorus compounds, T signifies the total assessment time interval (364 days), and   Q ef denotes the discharge flow rate (m3/d). The corresponding conversion factors for fi are reported in [2,5,32,41]. All the concentrations are addressed in g/m3 units.

2.3.2. Operational Cost Index (OCI)

Operational cost is a weighted summation of costs associated with the production of sludge (SP) (kg ss/d), methane (PM) (kg.CH4/d), pumping (PE), aeration (AE), mixing (ME), and heating (HE) energies (KWh/d); internal and external recycles are provided (m3/d). All individual components are addressed in [2,7]. Thus, the OCI is estimated as
OCI = AE + PE + z PS · SP + ME z PM · PM + max ( 0 , HE 7 PM )
where z i denotes the weighting factors; z PS is 3 and   z PM is 6.
Aeration, pumping, and mixing energies are addressed in Equations (10), (11) and (13). Here, aeration power is needed to aerate bioreactors; pumping is used to alter the flow rate from one end to another end and for internal, external flow patterns [2,7].
The aeration energy (AE) is described as (kWh/d):
AE = S O sat 1800 · T t o t f i = 1 7 V i · K L a i t   dt
where KLai signifies the oxygen mass transfer coefficient, Vi notifies the volume of the reactors and oxygen saturation coefficient. T is the length of evaluation time (364 days).
The pumping energy (PE) is defined as (kWh/d):
PE = 1 T t o t f 0.004 · Q int + 0.008 · Q r + 0.050 · Q w + 0.075 · Q PU + 0.060   Q TU + 0.004 · Q DO dt  
where Qintr is the internal recycle (m3/d), Qexr is the external recycle (m3/d), Qw is the waste flow (m3/d), QPU is the primary clarifier underflow, QTU is the thickener underflow, and QDO is the dewatering overflow.
The sludge production (SP) is expressed as (kg/d):
SP = 1 T · 1000 · X TSS t f X TSS t o + t o t f TSS X t Q X t   dt   S  
X TSS t =   X TSS ,   ASU t + X TSS ,   SSU t + X TSS ,   PSU t + X TSS ,   ADU t + X TSS ,   SU t   S  
with X TSS , X t =   TSS X t · V X .
Where QX (t) is the sludge flow and TSSX is the total amount of solids in the sludge flow stream (after dewatering in BSM2-P). XTSS is elucidated as the sum of TSS mass present in an individual process unit. The subscripts refer to the concern process units.
The mixing energy (ME) is defined as (kWh/d):
The mixing is highly necessary to avoid the biomass settling in the non-aerated and aerated reactors like all ASU tanks; anaerobic digester and the mixing energy (ME) is defined as (kWh/d):
ME = ME ASU + ME ADU
where
ME ASU = 24 T t o t f i 1 i = 7 if   K L a i < 20 d 1 0.005 · V i if   K L a i 20 d 1 0 · dt
ME ADU = 24 · 0.005 · V ADU
where, Vi is the ith tank volume (m3) and 0.005 kW/m3 is the mixing power consumption factor in ASU. VADU is the volume of liquid in ADU and the mixing power consumption factor 0.005 kW/m3.
Methane production (kg CH4/d) is defined as:
The average methane production per day value is defined by using Equation (14).
PM = P atm · 16 T · R · T OT t o t f 1 P tg t · P g , CH 4 t · Q g t · dt
where, PgCH4 (bar) is the partial pressure of methane gas produced in the headspace; R denotes the universal gas law constant i.e., 8.3145.10−2 bar m3 kmol−1 k−1; TOT represents the operating temperature of the digester (308.15 K); Ptg is the total gas pressure in the headspace; Patm is atmospheric pressure (1.013 bar); and Qg is the gas flow rate of produced gas.
Net heating energy is described as:
HE net = max 0 ,   HE 7 · PM
where HE is the amount of energy required to get the anaerobic digester up to operating temperature, as shown in the below Equation (16):
HE = 24 86400 · T t o t f P H 2 O · C H 2 O · T OT T adu , i t · Q ad t · dt
T ADU , i = T PSU · Q PSU t + T THK t · Q THK t Q ADU t
Here, Q ADU t = Q PSU t + Q THK t .
Where, PH2O is the density of water (1000 kg/m3), CH2O is the specific heat capacity of water (4.186 KJ kg−1 0C−1). Tad,i is the temperature of ADU influent, and TOP is the optimal temperature of ADU. Qad is the flow rate to the ADU (m3/d).
Following stringent regulations is a top priority for wastewater treatment plants. The legal constraints to be followed are the same as BSM1, i.e., TP is less than 2 gP/m3; TN is less than 18 gN/m3; biological oxygen demand (BOD5) is less than 10 g/m3; COD is less than 100 gCOD/m3; TSS is less than 30 g/m3; and SNH is less than 4 gN/m3.

3. Control Approaches

3.1. Design of Proportional-Integral (PI) Controller

The PI controllers can be framed using a wide variety of techniques accessible in the literature. In the present work, the Skogested internal model control (SIMC) method is used to design the PI controllers [44]. The first-order plus time order delay (FOPTD) model as prescribed in Equation (17) is identified for the design of the PI controllers for each loop.
G s = K P e θ ds τ p   S + 1
where, K P denotes the process gain, θ d denotes the delay, and τp signifies the time constant of the system. For more clarification, the method of approach is depicted in the flow diagram in the Supplementary data Figure S1. The identification method and designed controllers are elucidated distinctly. Figure S2 depicts the PI controller in the wastewater flow diagram.

3.2. Design of Model Predictive Controller (MPC)

MPC is a well-known advanced controller in which the objective is to predict the future move of the variable of significance and by using optimization, the corresponding control moves are determined systematically. At each control time interval, MPC predicts the nature of the variables of outputs over the prediction horizon (P). The control move sequence is computed over the control horizon (M) based on output predictions [23,45]. A quadratic cost function (Q) is selected for the operation of MPC, which is defined in the mode of Equation (18).
Q = i = 1 P | | Γ L ( L A + l / A r ( A + l ) ) | | 2 + i = 1 M | | Γ Δ v ( Δ v ( A + l 1 ) ) | | 2
Here, Γ L and Γ Δ v denotes the input and output rate weights respectively. L A + l / A is the controlled variable output at the future time interval of A + l , where the model is predicted at a particular time of A. The identification method and designed controllers are elucidated distinctly. Figure S3 depicts the MPC controller in the wastewater flow diagram. The third order state-space model of the system as shown in Equations (19) and (20) is as below:
x A + 1 = Gx A + Nu A
y A = Sx A + Zu A
Here x A signifies the state of vector, and G, N, S, and Z denote the state-space matrices.

3.3. Design of Fuzzy-Logic Controller (FLC)

Fuzzy logic is a widely used tool in a variety of control applications. FLC’s have been used in all stages of wastewater treatment processing. According to the literature, the most advanced control and processing units in WWTP are solved using fuzzy control or rule (FLC). This is achieved by employing fuzzy rules that are identical to those used in human inference design. FLC is used on the WWTP in this study [19,23]. Figure S4 depicts the FLC controller in the wastewater flow diagram.

4. Results and Discussion

Implementation of control strategies on plant wide-models: The results after implementing all the four control strategies are discussed here.

4.1. PI Control (2 Loops)

In ASU, SO,7, in the last aerobic tank (tank7), and SNO,4, in the second reactor of anoxic (tank 4), are controlled. The corresponding input variables are mass transfer coefficient (KLa) and internal recycle (Qintr). The set points are chosen according to the requirements of the WWTP. Usually, the value of SO,7 in the oxic reactor is retained in the range of 1.5 to 4 gO2/m3, and the practiced value is 2 gO2/m3. Similarly, the advisable working values for the nitrate levels in the anoxic tank are in the range of 1–3 gN/m3 and the practiced value is 1 gN/m3 [21]. The corresponding control relevant models are developed using the attained open-loop data for each SNO,4 and SO,7 control loops. For the values of 88,000 and 73, the concentrations of NO and DO are reported as 2 gO2/m3 and 1 gN/m3 respectively. In the seventh and fourth reactors, a random input signal with a ± 10% variation in these values is given and the attained output data for NO and DO are collected. The NO and DO output data are now used to build state-space models and also first-order plus time order delay (FOPTD) models using the method of prediction error minimization. Based on these models, each loop is equipped with PI controllers that are designed using SIMC [44]. The respective obtained FOPTD model parameters are:
  • NO loop: KP = 0.000026144, τp = 0.012515, and θd = 0.000875.
  • DO loop: KP = 0.04538, τi = 0.010085, and θd = 0.
Based on these models, PI controllers are designed using the SIMC method and are obtained for NO and DO loops like Kc =35748.16, Ti= 0.01215, Kc =11.015, and Ti = 0.010085 respectively. With these controllers, the plant-wide model is simulated and the corresponding simulation results are given in Section 4.4. The resultant tracking performance of SO and SNO is depicted in Figure 2A,B. The default control approach associate with two control loops is shown in Figure 3A.

4.2. PI Control (One Loop)

In this approach, a close-loop control framework contains a single PI controller as shown in Figure 3B. According to this scheme, SO in the sixth tank is maintained at a set-point of 2 mgO2/l by regulating the KLa6. In addition, the oxygen mass transfer coefficients entering tank 5 and tank 7 are multiplied by a factor of 1 and 0.5, respectively, to have improved manipulation and thus efficient aeration in these tanks [30].

4.3. Ammonia-Based Aeration Control (ABAC) Approach

In this approach, cascade configuration, as shown in Figure 3C, is used, in which MPC/Fuzzy controllers for ammonia (SNH,6) control are used by manipulating SO,6 set point in the aeration tank 6. Here, SO,6 in aeration tank 6 is controlled by regulating the airflow rates of reactors 5, 6, and 7 as implemented in strategy 2. PI controller is designed at the lower level with combinations of fuzzy and MPC at the supervisory level. PI-MPC and PI-fuzzy combinations are proposed as overall control structures. The corresponding list of all four different control approaches is listed in Table 4.

4.3.1. PI-MPC Control Strategy

Here, PI controller is implemented for BSM2-P at the lower level and MPC is implemented in the supervisory layer. To design MPC in the supervisory layer, system identification of the model in the supervisory layer is carried out. In which, the SO set-point is varied by ± 10% around the operating point, and the resulting SNH concentration values are collected. The prediction error method is used to identify a third order state-space model based on this data set [46]. The identified state-space model for the supervisory layer is expressed below Equation (21):
A = 0.7231 0.1351 0.03826 0.3957 0.2845 0.01298 0.0068 0.0477 0.1704       B = 0.09427 0.75 1.994 C = 1.306 0.074 0.01624   D = 0
Based on this model, MPC is designed with a sampling time of control as 0.05 days (72 min); prediction and control horizons of 10 and 2, respectively; and rate of change of regulated variable as 0.1. The corresponding results are given in Figure 4 in which the computed values of SO by supervisory MPC controller and its tracking by the lower-layer controller for 245–252 days are shown. The performance evaluation was done in the period of 245–609 days. Figure 4 depicts that a good supervisory set-point tracking is achieved by using the PI-MPC controller design framework in the sixth bioreactor. The resultant average concentrations of nutrient removal, energy usages, and greenhouse gas emissions, and performance of the plant with cost assessment are reported in Section 4.4 and compared with the other three control frameworks.

4.3.2. PI-Fuzzy Control Strategy

In this case, the fuzzy controller manipulates the SO set-points at a supervisory layer to minimize ammonia peaks. The membership functions (MF) of SO,6 and SNH,6 are considered in the range of 0–4 mg O2/l and 0.1–20 mg N/l, respectively. The MF’s for both input and output variables are in a Gaussian bell-shaped curve, which is divided into three linguistic variables, “high,” “low,” and “medium,” as shown in Figure 5A,B. A total of three rules are considered according to the SO control loop requirement [23,28].
  • If SNH level is “low”, then SO level is “low”
  • If SNH level is “high”, then SO level is “high”
  • If SNH level is “medium”, then SO level is “medium”
The corresponding linguistic variables for “high,” “low,” and “medium,” are given in Table 5. With these controllers, closed-loop simulations are carried out and the results are given in Figure 5C. It can be observed that a good supervisory set-point tracking is achieved by using the PI-fuzzy controller design framework. The effluent quality and operational cost details along with other quantitative results are given in Section 4.4.

4.4. Comparison of Four Control Design Frameworks on BSM2-P

Nitrification oxidizes ammonium to nitrate, and denitrification reduces nitrate to nitrogen gas. Then a high DO improves nitrification, but an excess of nitrate perhaps is not fully denitrified in the anoxic reactors due to a lack of COD. Moreover, phosphorous removal is largely influenced by dissolved oxygen, which is directly proportional to the formation of orthophosphates. The effluent concentrations of ammonia, TN, TP, and TSS are compared for all four control frameworks with their corresponding pollutant limit value. These are depicted in Figure 6A–D.
Based on the simulation results of all four control designs [PI controllers (two loops), PI controller (one loop), PI–MPC at supervisory level, and PI–fuzzy at a supervisory level], the corresponding average values of effluent concentrations are given in Table 6. From Table 6, comparing with only lower-level controllers, the ammonia, TP, TSS, and TN removal concentrations are improved. For ammonia, TP, TSS, and TN, the removal efficiency is improved by 36%, 33.6%, 1.02%, and 11.3% in PI-MPC, PI-fuzzy, PI (one loop), and PI (two loops) respectively.
Figure 6E depicts the average values of energy usages like aeration, pumping, mixing, heating, consumed energies (kwh/d), and sludge production rate (kg ss/d). From the graph, it was observed that energy consumption was high in the case of PI-MPC and low in the case of the PI-fuzzy controller. The sludge production rate and heat energy are high with the PI-fuzzy controller. In comparison, the fuzzy controller shows less aeration energy to obtain better nutrient removal with a slight increase in the cost. In the literature, also it was showed that the fuzzy control is favorable for better removal of phosphorous. In comparison, PI (one loop) and PI-fuzzy showed an improved EQI of 21.1% with a 0.52% increase in OCI. In the comparison of all four control strategies, there is a trade-off between OCI and EQI. Overall, in comparison with the PI (one loop) controller, ammonia removal is improved with PI-MPC, and TP is improved with PI-fuzzy strategies. In both cases, PI-MPC and PI-fuzzy showed improved EQI with an increase of OCI. The corresponding greenhouse gases like methane, hydrogen, and carbon dioxides of averages production rates are also reported in Table 6. It is noticed that PI-fuzzy shows high methane, hydrogen, carbon dioxide, and gas flow production rates of 28.7%, 4.87%, 6.8%, and 3.2% in comparison with PI (one loop). PI-fuzzy and PI-MPC showed good improvement in minimizing the percentage of violations in TP and ammonia. The percentage of violations is also reported in Table 6. PI-fuzzy showed good removal efficiency in the phosphorus. Moreover, PI-fuzzy showed lower OCI when compared with PI-MPC.

4.5. Summary of Previous and Present Studies on BSM2-P

For the biological seven-reactor A2O process, different control frameworks have been designed [47,48]. Optimal set-points are computed by higher-level control loops [49,50]. In [30] used BSM2-P as a simulation platform. PI controllers are designed to control DO by regulating the oxygen mass transfer coefficient with additional aeration-based ammonia controller and TSS control by regulating the external recycle. The control strategy showed significant improvement in both effluent quality and operating costs. The designed controllers showed improved EQI of 31% with decreased OCI of 6.9% when compared with open loop. As far as pollutant concentration is concerned, TN and TP are improved by 17% and 42.1% respectively.
In the present work, an ammonia-based aeration controller at the supervisory level is designed. By using two different control combinations, PI-MPC and PI-fuzzy, the performance is compared with [30] PI-fuzzy showed improved EQI of 13.5% with an increase of 13.6% OCI. Phosphorus and TN removal are improved by 29.7% and 5.4%, respectively, with the PI-fuzzy control framework and produce better EQI. PI-fuzzy shows a high production rate of methane when compared to [30].

5. Conclusions

In this study, MPC and fuzzy were designed at the supervisory level, and PI was designed for lower-level control for BSM2-P in ASM2d as an activated sludge model. A total of four control frameworks were implemented to evaluate and test the plant performance, concentrations, as well as effluent quality. The resultant performance indices were compared with the PI strategy. In each control application case, there was a trade-off between EQI and OCI. In comparison with PI (one loop), PI-fuzzy showed an improved EQI of 21.1% with a 0.52% increase in OCI. Of all the compared outcomes, PI-fuzzy showed better EQI and increased OCI. When comparing all four control strategies, it was reported that average effluent pollutant concentrations like BOD5, COD, TN, ammonia, and TSS attained the regulatory limits, except for phosphorus. Optimized ammonia removal was noticed in PI-MPC, whereas better optimized phosphorous removal was noticed in PI-fuzzy. PI-fuzzy showed high production rates of greenhouse gas emissions and low consumption of aeration energy. The percentage of violations of total phosphorus was less in the case of PI-fuzzy.

Supplementary Materials

The following are available online athttps://www.mdpi.com/article/10.3390/en14196386/s1, Figure S1: Systematic procedure for model identification and control implementation, Figure S2: A simple PI controller in WWTP, Figure S3: Model Predictive Controller scheme in WWTP, Figure S4: A classical Fuzzy controller in WWTP, Table S1: State variables of ASM2d, units with notations, and average influent data.

Author Contributions

Conceptualization, A.G.S., R.V. and S.R.A.; methodology, A.G.S., E.T., M.M., M.M.S., S.R.A.; software, A.G.S.; validation, A.G.S., M.M., R.V., S.R.A. and E.T.; formal analysis, A.G.S., M.M., M.M.S., R.V., S.R.A.; investigation, A.G.S., M.M., R.V., S.R.A. and E.T.; resources, A.G.S., M.M., R.V., S.R.A. and M.M.S.; data curation, A.G.S. and E.T.; writing—original draft preparation, A.G.S., E.T.; writing—review and editing, M.M., R.V., M.M.S., S.R.A. and E.T.; supervision, M.M.S., R.V., S.R.A.; project administration, M.M.S., R.V. and S.R.A.; funding acquisition, M.M., R.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Catalan Government under Projects 2017 SGR 1202. The APC was funded by the Spanish Government under Project PID2019-105434RBC33 co-funded with European Regional Development Funds of the European Union.

Data Availability Statement

Data is available upon request from the authors.

Acknowledgments

Xavier Flores-Alsina, Senior Researcher, Denmark Technical University, Denmark, is acknowledged for sending the BSM2-P Matlab/Simulink codes.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AEAeration Energy (Kwh/d)
ASM1Activated Sludge Model No.1
ASM2Activated Sludge Model No.2
ASM2d Activated Sludge Model No.2d
ASM3Activated Sludge Model No.3
BOD5Biological Oxygen Demand
CODChemical Oxygen Demand
DODissolved Oxygen
EQIEffluent Quality Index
IQIInfluent Quality Index
KProportional gain
KLaOxygen transfer coefficient
TNTotal Nitrogen
TKNTotal Kjeldahl Nitrogen
TSSTotal suspended solids
NONitrate
PPhosphorus
PEPumping Energy (kWh/d)
KUtPollutant load corresponding to component
QoInfluent flow rate (m3/d)
QintrInternal recycle flow rate (m3/d)
QrReturn sludge flow rate (m3/d)
QwWaste sludge flow rate (m3/d)
SAFermentation products (g COD/m3)
SFReadily biodegradable organic substrate
SHCOAlkalinity of the waste water (HCO3/m3)
SIInert soluble organic material (g COD/m3)
SNHAmmonium and ammonia nitrogen (g N/m3)
SNONitrate and nitrite nitrogen (g N/m3)
SN2Dinitrogen (g N/m3)
SPO4Inorganic soluble phosphate (g P/m3)
SSReadily biodegradable organic substrate (g COD/m3)
toStart time
tfEnd time
BODefTotal BOD concentration
CODefTotal COD concentration
SNONitrate concentration
SPorgTotal N concentration
SPinorgTotal phosphorus concentration
TKNTotal organic N concentration
WWTP Wastewater Treatment Plant
XANitrifying organisms (g COD/m3)
XHHeterotrophic organisms (g COD/m3)
XIInert particulate organic material (g COD/m3)
XSSlowly biodegradable substrates (g COD/m3)
XPAOPoly accumulating organisms (g COD/m3)
XPHACell internal storage product of PAO’S (g COD/m3)
XPPPolyphosphate (g P/m3)
XSTOCell inner storage product of heterotopy
XTSSSuspended solids (g SS/m3)

References

  1. Matejczyk, M.; Ofman, P.; Dąbrowska, K.; Świsłocka, R.; Lewandowski, W. The study of biological activity of transformation products of diclofenac and its interaction with chlorogenic acid. J. Environ. Sci. 2020, 91, 128–141. [Google Scholar] [CrossRef] [PubMed]
  2. Solon, K. Extending Wastewater Treatment Process Models for Phosphorus Removal and Recovery: A Framework for Plant-Wide Approach of Phosphorus, Sulfur and Iron. Ph.D. Thesis, Lund University, Sweden, 2017. [Google Scholar]
  3. Asadi, A.; Verma, A.; Yang, K.; Mejabi, B. Wastewater treatment aeration process optimization: A data mining approach. J. Environ. Manag. 2017, 203, 630–639. [Google Scholar] [CrossRef] [PubMed]
  4. Olsson, G.; Newell, B. Wastewater Treatment Systems: Modelling, Diagnosis and Control; IWA Publishing: London, UK, 1999. [Google Scholar]
  5. Henze, M.; Gujer, W.; Mino, T.; van Loosdrecht, M.C.M. Activated Sludge Models ASM1, ASM2, ASM2d, and ASM3; IWA Scientific and Technical Report, No. 9; IWA Publishing: London, UK, 2000. [Google Scholar]
  6. Rieger, L.; Koch, G.; Kühni, M.; Gujer, W.; Siegrist, H. The EAWAG bio-P module for activated sludge model no. 3. Water Res. 2001, 35, 3887–3903. [Google Scholar] [CrossRef]
  7. Gernaey, K.V.; Jeppsson, U.; Vanrolleghem, P.A.; Copp, J.B. Benchmarking of Control Strategies for Wastewater Treatment Plants; IWA Scientific and Technical Report, No. 23; IWA Publishing: London, UK, 2014. [Google Scholar]
  8. Nopens, I.; Batstone, D.J.; Copp, J.B.; Jeppsson, U.; Volcke, E.; Alex, J.; Vanrolleghem, P.A. An ASM/ADM model interface for dynamic plantwide simulation. Water Res. 2009, 43, 1913–1923. [Google Scholar] [CrossRef]
  9. Jeppsson, U.; Pons, M.N.; Nopens, I.; Alex, J.; Copp, J.B.; Gernaey, K.V.; Rosén, C.; Steyer, J.P.; Vanrolleghem, P.A. Benchmark Simulation Model No. 2: General protocol and exploratory case studies. Water Sci. Technol. 2007, 56, 67–78. [Google Scholar] [CrossRef] [PubMed]
  10. Barker, P.S.; Dold, P.L. General model for biological nutrient removal activated sludge systems: Model presentation. Water Environ. Res. 1997, 69, 969–984. [Google Scholar] [CrossRef]
  11. Jeppsson, U.; Alex, J.; Batstone, D.; Benedetti, L.; Comas, J.; Copp, J.B.; Corominas, L.; Flores-Alsina, X.; Gernaey, K.V.; Nopens, I.; et al. Benchmark simulation models, quo vadis? Water Sci. Technol. 2013, 68, 1–15. [Google Scholar] [CrossRef] [Green Version]
  12. Ekama, G.A. Using bioprocess stoichiometry to build a plant-wide mass balance based steady-state WWTP model. Water Res. 2009, 43, 2101–2120. [Google Scholar] [CrossRef] [PubMed]
  13. Grau, P.; de Gracia, M.; Vanrolleghem, P.A.; Ayesa, E. A new plant-wide modelling methodology for WWTPs. Water Res. 2017, 41, 4357–4372. [Google Scholar] [CrossRef]
  14. Flores-Alsina, X.; Arnell, M.; Amerlinck, Y.; Corominas, L.; Gernaey, K.V.; Guo, L.; Lindblom, E.; Nopens, I.; Porro, J.; Shaw, A.; et al. Balancing effluent quality, economical cost and greenhouse gas emissions during the evaluation of plant-wide wastewater treatment plant control strategies. Sci. Total. Environ. 2014, 466, 616–624. [Google Scholar] [CrossRef]
  15. Flores-Alsina, X.; Saagi, R.; Lindblom, E.; Thirsing, C.; Thornberg, D.; Gernaey, K.V.; Jeppsson, U. Calibration and validation of a phenomenological influent pollutant disturbance scenario generator using full-scale data. Water Res. 2014, 51, 172–185. [Google Scholar] [CrossRef]
  16. Ruano, M.V.; Serralta, J.; Ribes, J.; Garcia-Usach, F.; Bouzas, A.; Barat, R. Application of the general model Biological Nutrient Removal Model No. 1 to upgrade two full-scale WWTPs. Environ. Technol. 2011, 33, 1005–1012. [Google Scholar] [CrossRef] [Green Version]
  17. Volcke, E.I.P.; Gernaey, K.V.; Vrecko, D.; Jeppsson, U.; van Loosdrecht, M.C.M.; Vanrolleghem, P.A. Plant-wide (BSM2) evaluation of reject water treatment with a SHARON-Anammox process. Water Sci. Technol. 2006, 54, 93–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Luca, L.; Pricopie, A.; Barbu, M.; Ifrim, G.; Caraman, S. Control strategies of phosphorus removal in wastewater treatment plants. In Proceedings of the 23rd International Conference on System Theory, Control and Computing (ICSTCC), Sinaia, Romania, 9–11 October 2019; pp. 236–241. [Google Scholar]
  19. Xu, H.; Vilanova, R. PI and fuzzy control for P-removal in wastewater treatment plant. Int. J. Comp. Commun. Control 2015, 10, 176–191. [Google Scholar] [CrossRef]
  20. Hongyang, X.; Pedret, C.; Santin, I.; Vilanova, R. Decentralized model predictive control for N and P removal in wastewater treatment plants. In Proceedings of the 22nd International Conference on System Theory, Control and Computing (ICSTCC), Sinaia, Romania, 10–12 October 2018; pp. 224–230. [Google Scholar]
  21. Shiek, A.G.; Machavolu, V.R.K.; Seepana, M.M.; Ambati, S.R. Design of control strategies for nutrient removal in a biological wastewater treatment process. Environ. Sci. Pollut. Res. 2021, 28, 12092–12106. [Google Scholar] [CrossRef] [PubMed]
  22. Maheswari, P.; Sheik, A.G.; Tejaswini, E.S.S.; Ambati, S.R. Nested control loop configuration for a three stage biological wastewater treatment process. Chem. Prod. Process Model. 2021, 16, 87–100. [Google Scholar] [CrossRef]
  23. Sheik, A.G.; Seepana, M.M.; Seshagiri Rao, A. Supervisory Control Configurations Design for Nitrogen and Phosphorus Removal in Wastewater Treatment Plants. Water Environ. Res. 2021. [Google Scholar] [CrossRef] [PubMed]
  24. Santin, I.; Pedret, C.; Meneses, M.; Vilanova, R. Artificial Neural Network for nitrogen and ammonia effluent limit violations risk detection in Wastewater Treatment Plants. In Proceedings of the International Conference on System Theory, Control and Computing (ICSTCC), Cheile Gradistei, Romania, 14–16 October 2015; pp. 589–594. [Google Scholar]
  25. Barbu, M.; Santin, I.; Vilanova, R. Applying Control Actions for the Water Line and Sludge Line to Increase the Wastewater Treatment Plants Performances. Ind. Eng. Chem. Res. 2018, 57, 5630–5638. [Google Scholar] [CrossRef]
  26. Revollar, S.; Vilanova, R.; Vega, P.; Franciscoand, M. Wastewater Treatment Plant Operation: Simple Control Schemes with a Holistic Perspective. Sustainability 2020, 12, 768. [Google Scholar] [CrossRef] [Green Version]
  27. Barbu, M.; Vilanova, R.; Meneses, M.; Santin, I. Global Evaluation of Wastewater Treatment Plants Control Strategies Including CO2 Emissions. IFAC Papers 2017, 50, 12956–12961. [Google Scholar] [CrossRef]
  28. Tejaswini, E.S.S.; Panjwani, S.; Babu, G.U.B.; Rao, A.S. Model Based Control of a Full-Scale Biological Wastewater Treatment Plant. IFAC Papers 2020, 53, 208–213. [Google Scholar] [CrossRef]
  29. Flores-Alsina, X.; Gernaey, K.V.; Jeppsson, U. Benchmarking biological nutrient removal in wastewater treatment plants: Influence of mathematical model assumptions. Water Sci. Technol. 2012, 65, 1496–1505. [Google Scholar] [CrossRef]
  30. Solon, K.; Flores-Alsina, X.; Mbamba, C.K.; Ikumi, D.; Volcke, E.I.P.; Vaneeckhaute, C.; Ekama, G.; Vanrolleghem, P.A.; Batstone, D.J.; Gernaey, K.V.; et al. Plant-wide modelling of phosphorus transformations in wastewater treatment systems: Impacts of control and operational strategies. Water Res. 2017, 113, 97–110. [Google Scholar] [CrossRef] [Green Version]
  31. Flores-Alsina, X.; Ramin, E.; Ikumi, D.; Harding, T.; Batstone, D.; Brouckaert, C.; Sotemann, S.; Gernaey, K.V. Assessment of sludge management strategies in wastewater treatment systems using a plant-wide approach. Water Res. 2021, 190, 116–714. [Google Scholar] [CrossRef]
  32. Flores-Alsina, X.; Solon, K.; Kazadi Mbamba, C.; Tait, S.; Jeppsson, U.; Gernaey, K.V.; Batstone, D.J. Modelling phosphorus, sulphur and iron interactions during the dynamic simulation of anaerobic digestion processes. Water Res. 2016, 95, 370–382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gernaey, K.V.; Flores-Alsina, X.; Rosen, C.; Benedetti, L.; Jeppsson, U. Dynamic influent pollutant disturbance scenario generation using a phenomenological modelling approach. Environ. Model. Softw. 2011, 26, 1255–1267. [Google Scholar] [CrossRef]
  34. Martin, C.; Vanrolleghem, P.A. Analysing, completing, and generating influent data for WWTP modelling: A critical review. Environ. Model. Softw. 2014, 60, 188–201. [Google Scholar] [CrossRef] [Green Version]
  35. Snip, L.J.P.; Flores-Alsina, X.; Aymerich, I.; Rodríguez-Mozaz, S.; Barcelo, D.; Plosz, B.G.; Corominas, L.; Rodriguez-Roda, I.; Jeppsson, U.; Gernaey, K.V. Generation of synthetic data to perform (micro) pollutant wastewater treatment modelling studies. Sci. Total Environ. 2016, 569, 278–290. [Google Scholar] [CrossRef] [Green Version]
  36. Otterpohl, R. Dynamische Simulation zur Unterstützung der Planung und des Betriebes von Kommunalen Kläranlagen [Dynamic Simulation to Support the Design and Operation of Municipal Wastewater Treatment Plants]. Ph.D. Thesis, Technical University of Aachen, Aachen, Germany, 1995. [Google Scholar]
  37. Takács, I.; Patry, G.G.; Nolasco, D. A dynamic model of the clarification-thickening process. Water Res. 1991, 25, 1263–1271. [Google Scholar] [CrossRef]
  38. Guerrero, J.; Flores-Alsina, X.; Guisasola, A.; Baeza, J.A.; Gernaey, K.V. Effect of nitrite, limited reactive settler and plant design configuration on the predicted performance of simultaneous C/N/P removal WWTPs. Bioresour. Technol. 2013, 136, 680–688. [Google Scholar] [CrossRef] [Green Version]
  39. Kazadi Mbamba, C.; Flores-Alsina, X.; Batstone, D.J.; Tait, S. Validation of a plant-wide phosphorus modelling approach with minerals precipitation in a full-scale WWTP. Water Res. 2016, 100, 169–183. [Google Scholar] [CrossRef] [Green Version]
  40. Flores-Alsina, X.; Ikumi, D.; Batstone, D.; Gernaey, K.V.; Brouckaert, C.; Ekama, G.A.; Jeppsson, U. Towards a plant-wide Benchmark Simulation Model with simultaneous nitrogen and phosphorus removal wastewater treatment processes. In Proceedings of the IWA 8th World Water Congress and Exhibition, Busan, South Korea, 16–21 September 2012. [Google Scholar]
  41. Henze, M.; Gujer, W.; Mino, T.; Matsuo, T.; Wentzel, M.C.; Marais, G.V.R.; van Loosdrecht, M.C.M. Activated Sludge Model No. 2d, ASM2d. Water Sci. Technol. 1999, 39, 165–182. [Google Scholar] [CrossRef]
  42. Batstone, D.J.; Keller, J.; Angelidaki, I.; Kalyuzhnyi, S.V.; Pavlostathis, S.G.; Rozzi, A.; Sanders, W.T.M.; Siegrist, H.; Vavilin, V.A. The IWA Anaerobic Digestion Model No. 1 (ADM1). Water Sci. Technol. 2002, 45, 65–73. [Google Scholar] [CrossRef]
  43. Gernaey, K.; Jeppsson, U.; Batstone, D.J.; Ingildsen, P. Impact of reactive settler models on simulated WWTP performance. Water Sci. Technol. 2006, 53, 159–167. [Google Scholar] [CrossRef] [PubMed]
  44. Grimholt, C.; Skogestad, S. Optimal PI and PID control of first-order plus delay processes and evaluation of the original and improved SIMC rules. J. Proc. Control 2018, 70, 36–46. [Google Scholar] [CrossRef]
  45. Maciejowski, J.M. Predictive Control with Constraints, 1st ed.; Pearson Education: Harlow, UK, 2002. [Google Scholar]
  46. Ljung, L. System Identification—A Theory for the User; Prentice Hall International: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  47. Ostace, G.S.; Baeza, J.A.; Guerrero, J.; Guisasola, A.; Cristea, V.M.; Agachi, P.Ş.; Lafuente, F.J. Development and economic assessment of different WWTP control strategies for optimal simultaneous removal of carbon, nitrogen and phosphorus. Comput. Chem. Eng. 2013, 53, 164–177. [Google Scholar] [CrossRef] [Green Version]
  48. Stare, A.; Vrečko, D.; Hvala, N.; Strmčnik, S. Comparison of control strategies for nitrogen removal in an activated sludge process in terms of operating costs: A simulation study. Water Res. 2007, 41, 2004–2014. [Google Scholar] [CrossRef]
  49. Guerrero, J.; Guisasola, A.; Vilanova, R.; Baeza, J.A. Improving the performance of a WWTP control system by model-based setpoint optimisation. Environ. Model. Softw. 2011, 26, 492–497. [Google Scholar] [CrossRef]
  50. Rojas, J.D.; Baeza, J.A.; Vilanova, R. Effect of the controller tuning on the performance of the BSM1 using a data driven approach. In Proceedings of the Watermatex 2011 8th IWA Symposium on Systems Analysis and Integrated Assessment, San Sebastián, Spain, 20–22 June 2011; pp. 785–792. [Google Scholar]
Figure 1. Proposed plant-wide model layout based on BSM2-P.
Figure 1. Proposed plant-wide model layout based on BSM2-P.
Energies 14 06386 g001
Figure 2. Set-point tracking for: (A) dissolved oxygen and (B) Nitrate.
Figure 2. Set-point tracking for: (A) dissolved oxygen and (B) Nitrate.
Energies 14 06386 g002
Figure 3. Control frameworks for BSM2-P: (A) PI controllers, (B) Lower level control, (C) Supervisory level control framework with a lower level.
Figure 3. Control frameworks for BSM2-P: (A) PI controllers, (B) Lower level control, (C) Supervisory level control framework with a lower level.
Energies 14 06386 g003aEnergies 14 06386 g003b
Figure 4. Dissolved oxygen tracking in the sixth reactor (PI-MPC).
Figure 4. Dissolved oxygen tracking in the sixth reactor (PI-MPC).
Energies 14 06386 g004
Figure 5. Membership functions rules: (A) Output, (B) Input, and (C) corresponding tracking performance.
Figure 5. Membership functions rules: (A) Output, (B) Input, and (C) corresponding tracking performance.
Energies 14 06386 g005aEnergies 14 06386 g005b
Figure 6. The effluent concentrations of (A) ammonia, (B) TN, (C) TP, (D) TSS, and (E) bar graph for all (SP, AE, PE, ME, HE, and consumed energy (CE)) the operational plant performance parameters are compared for four control frameworks with their corresponding average values.
Figure 6. The effluent concentrations of (A) ammonia, (B) TN, (C) TP, (D) TSS, and (E) bar graph for all (SP, AE, PE, ME, HE, and consumed energy (CE)) the operational plant performance parameters are compared for four control frameworks with their corresponding average values.
Energies 14 06386 g006aEnergies 14 06386 g006bEnergies 14 06386 g006c
Table 1. Average data of dynamic mass flow rates.
Table 1. Average data of dynamic mass flow rates.
Dynamic Mass Flow RatesAverage Data
Chemical oxygen demand (COD)8386 kg/d
Nitrogen1014 kg N/d
Phosphorus197 kg P/d
S:COD0.003 kg S Kg/COD
Table 2. Attributes of plant-wide model processes units and physical configurations.
Table 2. Attributes of plant-wide model processes units and physical configurations.
Process UnitWorking FunctionReferencesConfigurations
PSUNon-reactive[36]900 m3
SSUDouble-exponential velocity function reactive[37,38,39,40]6000 m3
ASUASM2d[32,41]4500 m3
ADUADM1[42]3400 m3
THKReactive[7,43]Underflow 30.9 m3/d
DUReactive[7,43]9.6 m3/d sludge and 168.9 m3/d reject water
SUNon-reactive[7]160 m3
Table 3. Weighting factors for EQI.
Table 3. Weighting factors for EQI.
Weighting Factors of EQI θ i
Weighting factors θ T S S θ C O D θ N K J θ N O θ B O D 5 θ P o r g θ P i n o r g
Value2130102100100
Table 4. Functioning of control strategies.
Table 4. Functioning of control strategies.
AttributesPI controller (Two Loops)PI controller (One Loop)PI (Lower Level) + MPC (Supervisory Level)PI (Lower Level) + Fuzzy (Supervisory Level)
Controlled variableSO,7 and SNO,4SO,6SNH,6SNH,6
Set-point2 gO2/m31 gN/m32 gO2/m3DO set-point is determined by higher levelDO set-point is determined by higher level
Regulating variablesKLa7 and internal recycleKLa in the last three reactorsSet-point for DO controllerSet-point for DO controller
Controller designPIPIPI and MPCPI and Fuzzy
Table 5. Linguistic functions and MF’s for control inputs and outputs.
Table 5. Linguistic functions and MF’s for control inputs and outputs.
Linguistic Variable (Output)
Linguistic ValueRangeMFCharacteristic Ranges
1LowerGaussian bell-shaped shaped0.1755.40.11
2MediumGaussian bell-shaped shaped1.065.871.36
3HigherGaussian bell-shaped shaped3.56186
Linguistic Variable (Input)
Linguistic ValueRangeMFCharacteristic Ranges
1LowerGaussian bell-shaped shaped1.899.180.034
2MediumGaussian bell-shaped shaped1.027.752.96
3HigherGaussian bell-shaped shaped8.2642.212.36
Table 6. Comparison of average concentration values of effluent discharge for four control strategies.
Table 6. Comparison of average concentration values of effluent discharge for four control strategies.
ParametersPI Controller (Two Loops)PI Controller (One Loop)PI-MPCPI-Fuzzy
SNH1.050.960.571.28
TSS15.3915.5415.3816.23
TN9.079.819.868.7
TP4.544.054.42.69
COD42.1742.0442.1741.95
BOD52.432.422.452.50
IQI97,875.7197,875.7197,875.7197,875.71
EQI14,625.9813,715.3714,39110,824.9
Average production rates
Methane1029102410351438
Hydrogen0.003920.003930.00390.0041
Carbon dioxide1504152715171640
Gas flow2630263526462722
OCI10,959.110,94911,81011,007
Average percentage of Violations
TP86.7172.570.438.4
NH2.713.120.211.28
TSS0.0250.0620.0480.58
BOD5-------0.0085
Nremoved/OCI0.081310.08000.07400.0815
Premoved/OCI0.010440.01130.00980.0139
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sheik, A.G.; Tejaswini, E.; Seepana, M.M.; Ambati, S.R.; Meneses, M.; Vilanova, R. Design of Feedback Control Strategies in a Plant-Wide Wastewater Treatment Plant for Simultaneous Evaluation of Economics, Energy Usage, and Removal of Nutrients. Energies 2021, 14, 6386. https://doi.org/10.3390/en14196386

AMA Style

Sheik AG, Tejaswini E, Seepana MM, Ambati SR, Meneses M, Vilanova R. Design of Feedback Control Strategies in a Plant-Wide Wastewater Treatment Plant for Simultaneous Evaluation of Economics, Energy Usage, and Removal of Nutrients. Energies. 2021; 14(19):6386. https://doi.org/10.3390/en14196386

Chicago/Turabian Style

Sheik, Abdul Gaffar, Eagalapati Tejaswini, Murali Mohan Seepana, Seshagiri Rao Ambati, Montse Meneses, and Ramon Vilanova. 2021. "Design of Feedback Control Strategies in a Plant-Wide Wastewater Treatment Plant for Simultaneous Evaluation of Economics, Energy Usage, and Removal of Nutrients" Energies 14, no. 19: 6386. https://doi.org/10.3390/en14196386

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop