Next Article in Journal
Application of Bayesian Hierarchical Negative Binomial Finite Mixture Model for Cost-Benefit Analysis of Barriers Optimization, Accounting for Severe Heterogeneity
Next Article in Special Issue
Locating the Parameters of RBF Networks Using a Hybrid Particle Swarm Optimization Method
Previous Article in Journal
On the Solutions of Second-Order Differential Equations with Polynomial Coefficients: Theory, Algorithm, Application
Previous Article in Special Issue
Parallelized Swarm Intelligence Approach for Solving TSP and JSSP Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant

by
Naresh Patnana
1,†,
Swapnajit Pattnaik
1,†,
Tarun Varshney
2,† and
Vinay Pratap Singh
3,*,†
1
Electrical Engineering, National Institute of Technology, Raipur 492010, India
2
Electrical and Electronics Engineering, ABES Engineering College, Ghaziabad 201009, India
3
Electrical Engineering, Malaviya National Institute of Technology, Jaipur 302017, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Algorithms 2020, 13(11), 287; https://doi.org/10.3390/a13110287
Submission received: 16 October 2020 / Revised: 5 November 2020 / Accepted: 5 November 2020 / Published: 10 November 2020
(This article belongs to the Special Issue Swarm Intelligence Applications and Algorithms)

Abstract

:
In this investigation, self-learning salp swarm optimization (SLSSO) based proportional- integral-derivative (PID) controllers are proposed for a Doha reverse osmosis desalination plant. Since the Doha reverse osmosis plant (DROP) is interacting with a two-input-two-output (TITO) system, a decoupler is designed to nullify the interaction dynamics. Once the decoupler is designed properly, two PID controllers are tuned for two non-interacting loops by minimizing the integral-square-error (ISE). The ISEs for two loops are obtained in terms of alpha and beta parameters to simplify the simulation. Thus designed ISEs are minimized using SLSSO algorithm. In order to show the effectiveness of the proposed algorithm, the controller tuning is also accomplished using some state-of-the-art algorithms. Further, statistical analysis is presented to prove the effectiveness of SLSSO. In addition, the time domain specifications are presented for different test cases. The step responses are also shown for fixed and variable reference inputs for two loops. The quantitative and qualitative results presented show the effectiveness of SLSSO for the DROP system.

1. Introduction

The existence of civilization depends upon the availability of potable water. The studies are repeatedly indicating that the potable water is acutely depleted in the last few decades. The main reason for depletion of fresh water is overuse of ground water due to growing population and industrialization [1]. In order to meet the requirements due to increased population and industrialization, it is required to explore for sources of fresh water. The desalination of sea water is one of the possible solutions since the expense of sea water is abundant [2].
The process of removing salts from saline water is called desalination [3]. Desalination is achieved by (a) thermal distillation processes such as multiple effect evaporation, multi-stage flash (MSF), vapour compression, etc., and (b) membrane processes such as nano-filtration, electro-dialysis reversal, forward osmosis, reverse osmosis (RO), etc. Among these, RO has been proved to be economical [4]. RO process also produces high quality water with minimal energy consumption. Additionally, there is low initial investment in case of RO process.
RO is a filtration process which segregates fresh water from higher concentrated saline water. In RO process, a semi-permeable membrane is used which passes the fresh water through it and rejects the concentrated solution. The various plants which use the RO process are Doha RO system [5], Chaabene RO system [6], Riverol and Pilipovik [7], etc. The various pilot-scale models [8] for RO are also presented in literature. The RO models given in literature are interacting having flux and conductivity as controlled variables generally.
The most common controller used for interacting RO models is a proportional-integral-derivative (PID) controller. The PID controller is used due to its simple structure and being able to produce satisfactory performance. The performance of PID controller depends on the quality of tuning. The classical methods [9] used for tuning are Ziegler–Nichols, gain-phase margin method, Cohen–Coon, etc. However, the main drawback with classical tuning methods is moderate performance of the controller [10]. The controller performance can be enhanced by improving the quality of tuning.
Various optimization based methods can guarantee the improved quality of tuning [11]. Recently, the intelligent optimization algorithms like fireworks algorithm [12], bacterial foraging optimization [13], elephant herding optimization [10], whale optimization algorithm [2,14], monarch butterfly optimization [15], brain storm optimization [16], etc. are proposed in the existing literature. Some techniques are proposed in the literature for tuning the PID controller for RO systems. In [11], PID controller is designed using Luus–Jaakola (LJ) optimization algorithm for the RO system. Particle swarm optimization (PSO) is utilized for tuning the PID controller parameters in [17].
Grey wolf optimization (GWO) assisted PID tuning for RO system is presented in [4]. Similarly, whale optimization algorithm and elephant herding optimization based controllers are suggested in [2,10], respectively. Other relevant controller designs can be found in [18,19]. The methods presented in [2,4,10,19] also provide better response when compared to classical methods of tuning of PID controller. However, the methods proposed in [2,4,10,19] are providing the tuning with basic algorithms only. It is well known that the basic algorithms have the problem of local minima trapping [20]. So sometimes, it becomes necessary to enhance the performance of basic algorithm by some modifications [21,22].
In this investigation, an improved algorithm, i.e., self learning salp swarm optimization (SLSSO), is utilized for tuning of PID controllers [23]. The controller tuning is presented for Doha reverse osmosis plant (DROP). DROP model is having two controlled variables namely, flux and conductivity and two manipulated variables namely, pressure and pH. These four variables are interacting in nature. This makes DROP model an interacting two-input-two-output (TITO) system. First, a decoupler is designed using feed-forward method [24] of decoupling to make the system non-interacting in nature. Thus obtained model has two non-interacting loops, one for flux and other for conductivity. Then, two PID controllers are obtained for two non-interacting loops using SLSSO algorithm. The integral-of-squared-error (ISE) is taken as tuning criterion. For the ease of simulation, the ISE is formulated in terms of alpha and beta parameters. The key features of proposed method are: (a) it is utilizing SLSSO for tuning of PID controllers, (b) a decoupler is designed for interacting Doha RO system, (c) the ISE is formulated in terms of alpha and beta parameters, and (d) the results of SLSSO are also compared with other algorithms to prove its supremacy. This study also obtained the PID controllers using artificial bee colony (ABC) algorithm, LJ algorithm, Nelder–Mead simplex (NMS) algorithm, sine cosine algorithm (SCA) and salp swarm optimization (SSO) algorithm for fair comparisons of performance of SLSSO tuned controllers. Additionally, time-domain specifications are tabulated for comparisons. Further, the statistical analysis of ISE obtained for controllers is presented. Moreover, the closed-loop responses (step, impulse and variable step responses) are also shown for different controllers.
The structure of the paper is as follows: Section 2 presents the description of reverse osmosis system in consideration. The interacting TITO system is given in the same section. In addition, the decoupler design is described in this section. The controller structure is presented in Section 3. Section 4 is presenting the SSO algorithm. The description of SLSSO algorithm is also given in the same section. Section 5 presents results and discussion where two case studies are discussed in detail. Finally, the findings of paper are concluded in Section 6.

2. Reverse Osmosis System

The schematic block diagram of Doha reverse osmosis plant (DROP) system is provided in Figure 1 [4]. The RO desalination plant comprises of four stages, namely, pre-treatment, high pressure pump, RO membrane assembly and post-treatment.
During pre-treatment stage, the saline water undergoes filtration to remove any physical or chemical impurities. This stage is important to ensure membrane longevity. The organic chemicals which can burn holes in the membrane and can cause irreparable damage to membrane are trapped using activated cartridge filters. In addition, some inhibitors are added to decrease the pH of the feed-water which helps in reducing the scaling effect. Besides this, physical impurities like mud particles and microbial growth on membrane are also removed during this stage to increase the life of membrane.
High pressure pump is utilized to provide the necessary pressure required for reverse osmosis. For brackish water, it is in the range of 15–25 bars while it is 54–80 bars for sea water. The pressurized feed-water is then fed to the membrane assembly which is made up of two or more hollow or spiral wound composite polyamide membranes. This is the main part of the RO system. It does not allow impurities to pass through it. The membrane assembly must be able to withstand the high pressure. The rejected brine from membrane assembly is sent to the discharge channel. In post-treatment stage, calcium and bi-carbonate ions are removed from the water. In addition, the pH is raised slightly above 7 and water is prepared for distribution. The two controlled variables (flux and conductivity of output water) are manipulated using two manipulated variables (feed-water pressure and feed-water pH) in case of DROP model.

2.1. Control-Loops for DROP Model

The block diagram of DROP is given in Figure 2. The transfer function of interacting DROP model is given as
Y 1 ( s ) Y 2 ( s ) = D 11 ( s ) D 12 ( s ) D 21 ( s ) D 22 ( s ) R 1 ( s ) R 2 ( s )
where R 1 ( s ) and R 2 ( s ) are two manipulated variables which are feed-water pressure and feed-water pH, respectively, and Y 1 ( s ) and Y 2 ( s ) are two controlled variables which are flow and conductivity of output water, respectively. The transfer functions D 11 ( s ) , D 12 ( s ) , D 21 ( s ) and D 22 ( s ) are denoting the individual transfer functions in between two manipulated variables and two controlled variables.
The transfer functions D 11 ( s ) , D 12 ( s ) , D 21 ( s ) and D 22 ( s ) are given as:
D 11 ( s ) = k 11 ( n 11 s + 1 ) s 2 + 2 ζ 11 d 11 s + d 11
D 12 ( s ) = 0
D 21 ( s ) = k 21 ( n 21 s + 1 ) s 2 + 2 ζ 21 d 11 s + d 21
D 22 ( s ) = k 22 ( n 22 s + 1 ) s 2 + 2 ζ 22 d 22 s + d 22
It can easily be concluded from transfer functions (2)–(5) that both flow and conductivity are manipulated by feed-water pressure. However, the feed-water pH is affecting only conductivity. The values of different parameters of transfer functions (2)–(5), are given in Table 1 [4]. The interval transfer functions of (2)–(5) turn out to be (2)–(5) by considering 40% uncertainty.
D 11 ( s ) = 0.402 , 0.938 + 0.0225 , 0.0525 s 10.956 , 25.564 + 19.938 , 46.522 s + 0.60 , 1.40 s 2
D 12 ( s ) = 0
D 21 ( s ) = 3.346 , 1.434 + 1.1711 , 0.5019 s 1.302 , 3.038 + 16.65 , 38.85 s + 0.60 , 1.40 s 2
D 22 ( s ) = 113 , 57 + 42.63 , 18.27 s 0.774 , 1.806 + 1.794 , 4.186 s + 0.60 , 1.40 s 2

2.2. Design of Decoupler

A decoupler is required to be designed to nullify the effect of interaction for interacting TITO system. Figure 3 shows the decoupler along with the plant. This work employs simplified decoupling technique [17] to convert interacting TITO system into two non-interacting control loops. The transfer function of decoupler is obtained by
T = C p C d
such that
T = T 11 ( s ) 0 0 T 22 ( s )
C d = 1 x y 1
where C p and C d are, respectively, the transfer functions of plant and decoupler. From (10), it is obtained
x = D 12 ( s ) D 22 ( s )
y = D 21 ( s ) D 11 ( s )
T 11 = D 11 ( s ) D 22 ( s ) D 12 ( s ) D 21 ( s ) D 22 ( s )
T 22 = D 11 ( s ) D 22 ( s ) D 12 ( s ) D 21 ( s ) D 11 ( s )
Hence, it is clear from (13)–(16) that by choosing x and y as given in (13) and (14), respectively, the interacting system is decoupled into two non-interacting loops with transfer functions, T 11 and T 22 as given in (15) and (16), respectively. Thus, two independent PID controllers, C P I D 1 and C P I D 2 , can be obtained for two non-interacting loops having transfer functions, namely T 11 and T 22 .

3. Controller Design

3.1. Controller Structure

In literature, many structures of PID controller are existing. However, the simple and widely used parallel form given in (17) is used in this work.
u ( t ) = K p e ( t ) + K i 0 t e ( t ) d t + K d e ( t )
The parameters, K p , K i and K d are, respectively, the proportional, integral and derivative gains. However, u ( t ) and e ( t ) denote the control effort generated by controller and error signal fed to controller, respectively. In Laplace domain, (17) becomes as given in (18).
C P I D = U ( s ) E ( s ) = K p + K i s + K d s

3.2. The Determination of ISEs for DOHA Plant

The ISE is utilized for obtaining the controller parameters. The ISE is obtained in terms of alpha and beta parameters. The ISE for two independent loops, i.e., pressure-flux and pH-conductivity loops, is given below.
(a)
ISE for pressure-flux loop:
I S E 1 = 0 e 1 2 ( t ) d t = i = 1 k β i 2 2 α i
(b)
ISE for pH-conductivity loop:
I S E 2 = 0 e 2 2 ( t ) d t = i = 1 k ¯ β ¯ i 2 2 α ¯ i
The alpha and beta parameters are obtained utilizing alpha and beta tables. The alpha and beta tables are derived from respective error in Laplace domain. Suppose the error in Laplace domain is given as
E ( s ) = U m 1 s m 1 + + U 1 s + U 0 V m s m + V m 1 s m 1 + + V 1 s + V 0
The alpha and beta tables are constructed as given in Table 2 and Table 3, respectively [4]. The first two rows of Table 2 and Table 3 are obtained from (21). The rest of the elements are obtained using elements of these two rows.

4. Salp Swarm Optimization

Salp swarm optimization (SSO) [25] is a nature inspired swarm-based meta-heuristic optimization algorithm. It is proposed recently by Mirjalili et al. [26] in 2017. This algorithm is derived from the swarming behavior of salps. Salps are found in deep ocean and are very similar to jelly fishes. The salps’ population can be categorized into (i) leaders and (ii) followers. The leaders lead the whole population while followers follow the instructions of leaders, directly or indirectly. The salps form chain to move. The salps at front of chain are known as leaders and remaining salps are known as followers. The better half of the population is treated as leader salps. However, remaining salps are considered as followers.
At first, the candidate solutions for leaders are updated. Since, followers follow the leaders, the candidate solutions for followers are updated using the solutions obtained for leaders. Suppose, X i , j , i = 1 , 2 , , M and j = 1 , 2 , , N denote initial candidate solutions for whole population where M and N represent population size and number of decision variables, respectively. The candidate solutions for leaders and followers are updated as follows.

4.1. Update of Candidate Solutions for Leaders

The candidate solutions for leaders are updated as
X i , j n e w = T j + k 1 ( ( X j + X j ) k 2 + X j ) k 3 0.5 T j k 1 ( ( X j + X j ) k 2 + X j ) k 3 < 0.5
where X i , j n e w is updated candidate solutions for X i , j , T j is the position of food source, X j and X j + represent the minimum and maximum values of decision variables, k 2 and k 3 are random numbers distributed uniformly in the range [0,1], and k 1 is modified during course of iterations as
k 1 = 2 e ( 4 t T ) 2
where T and t, respectively, are maximum number of iterations and current iteration.

4.2. Update of Candidate Solutions for Followers

The candidate solutions for followers are updated using the solutions of leaders. The mathematical expression used for updating the candidate solutions for followers is
X i , j n e w = X i , j + X i 1 , j 2
where X i , j n e w denotes updated candidate solution for follower X i , j .
After updating the whole population as suggested in (22) and (24), candidate solutions violating the minimum and maximum values of decision variables are reinitialized at respective minimum and maximum values of decision variables. The pseudo code of SSO is presented in Figure 4.

4.3. Self-Learning Salp Swarm Optimization

A self-learning rule exploits the space in close proximity of the individual position to reach its global optimum [23]. The rule provides an opportunity to each learner to enhance the individual knowledge from self-surrounding space. This phase can be mathematically expressed as:
X i , j n e w ( k ) = X i , j ( k ) ( 1 + λ ( r 0.5 ) )
where, X i , j n e w ( k ) is the new solution vector in this phase, λ is the self-learning factor which determines the self-learning capability of each individual and r [ 0 , 1 ] . The value of λ is considered 3 in this work. Update the solution vector, X i , j ( k ) , using greedy selection. Thus, the updated X i , j ( k ) after this phase takes part in the next iteration.

5. Results and Comparisons

Two case studies are performed to highlight the findings and contributions of this article. One proportional-integral-derivative (PID) controller is designed for pressure-flux loop and other PID controller is tuned for pH-conductivity loop. In both the case studies, ISE of respective loop is minimized using SLSSO algorithm [23]. For fair presentation of efficacy of SLSSO algorithm, other algorithms namely, SSO, ABC, LJ, NMS, SCA and PSO are applied to tune the PID controller.

5.1. Case Study I

The results presented in this case study are obtained for pressure-flux loop by minimizing the performance index given in (19). Table 4 presents PID gains i.e., K p , K i , and K d obtained using SLSSO. Same table also tabulates the PID gains obtained using ABC, LJ, NMS, PSO, SCA and SSO. The mean and minimum values of ISE along with standard deviation are also presented in this table. The time domain specifications of step response obtained due to all algorithms (ABC, LJ, NMS, PSO, SCA, SSO and SLSSO) are tabulated in Table 5. The Figure 5 shows the step response obtained due to ABC, LJ, NMS, PSO, SCA, SSO and SLSSO. In addition to this, the response of pressure-flux loop for variable input due to SLSSO algorithm is plotted in Figure 6. The Figure 7 is showing the response with disturbance.
The mean value of ISE is minimum in case of SLSSO algorithm as shown in Table 4. The same is true for minimum value of ISE and standard deviation. The time domain specifications tabulated in Table 5 also supports the response obtained for pressure-flux loop. The step response of SLSSO algorithm is also reaching steady state as shown in Figure 5. For variable input and input having disturbance also, the SLSSO tuned controller is providing stable step response as shown in Figure 6 and Figure 7. Hence, it can be inferred that SLSSO algorithm is a better choice for tuning the PID gains for pressure-flux loop of DROP system.

5.2. Case Study II

This case study presents the analysis of results obtained for pH-conductivity loop. In this case, the performance index given in (20) is minimized using SLSSO algorithm. The PID gains ( K p , K i , and K d ) obtained using SLSSO are presented in Table 6. The Table 6 also provides gains obtained due to ABC, LJ, NMS, PSO, SCA and SSO. The mean and minimum values of ISE along with standard deviation are also shown in Table 6. The time domain specifications (rise time and settling time) of step response are given in Table 7 which are obtained due to ABC, LJ, NMS, PSO, SCA, SSO and SLSSO. The Figure 8 plots the step response of pH-conductivity loop due to ABC, LJ, NMS, PSO, SCA, SSO and SLSSO. Moreover, the response of pH-conductivity loop for variable step input and step input having disturbance due to SLSSO algorithm is presented in Figure 9 and Figure 10, respectively.
In this case also, the mean value of ISE is the minimum for SLSSO algorithm as shown in Table 6. The same conclusion can be derived for minimum value of ISE and standard deviation. Table 6, Figure 8, Figure 9 and Figure 10 also prove that the SLSSO algorithm is an excellent alternative for tuning of PID controller for DROP system.

6. Conclusions

This contribution presented design of proportional-integral-derivative (PID) controller for Doha reverse osmosis plant (DROP). The DROP system is interacting in nature and has two inputs and two outputs resulting in TITO system. Since, this TITO system is interacting, the interaction effect should be eliminated before designing the PID controllers for two manipulated variables. In this article, simplified decoupling technique is adopted to design decoupler for interacting TITO DROP system. After designing the decopuler, two PID controllers for pressure flux and pH conductivity loops are designed to manipulate the pressure and pH in order to generate specified flux and conductivity. For designing the PID controllers, integral-square-error (ISE) of the respective loop is minimized using lelf-learning salp swarm optimization (SLSSO) algorithm [23]. It was proved that SLSSO algorithm is an excellent choice for tuning the PID gains for two loops of the DORP system.
As far as the further extension of this work is considered, it is aimed to obtain the interval model in discrete domain. The reduced-order model for discrete domain will be designed in further research. It would also be interesting to investigate the controller design using intelligent optimization algorithms like fireworks algorithm [12], bacterial foraging optimization [13], elephant herding optimization [10], whale optimization algorithm [2,14], monarch butterfly optimization [15], brain storm optimization [16], etc.

Author Contributions

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

Funding

This research was funded by SERB, grant number ECR/2017/000212.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rathore, N.S.; Singh, V.; Phuc, B.D.H. A modified controller design based on symbiotic organisms search optimization for desalination system. J. Water Supply Res. Technol. 2019, 68, 337–345. [Google Scholar] [CrossRef]
  2. Rathore, N.S.; Singh, V. Whale optimisation algorithm-based controller design for reverse osmosis desalination plants. Int. J. Intell. Eng. Inform. 2019, 7, 77–88. [Google Scholar] [CrossRef]
  3. Rathore, N.S.; Singh, V. Design of optimal PID controller for the reverse osmosis using teacher-learner-based-optimization. Membr. Water Treat. 2018, 9, 129–136. [Google Scholar]
  4. Rathore, N.S.; Singh, V.; Kumar, B. Controller design for doha water treatment plant using grey wolf optimization. J. Intell. Fuzzy Syst. 2018, 35, 5329–5336. [Google Scholar] [CrossRef]
  5. Alatiqi, I.; Ghabris, A.; Ebrahim, S. System identification and control of reverse osmosis desalination. Desalination 1989, 75, 119–140. [Google Scholar] [CrossRef]
  6. Chaaben, A.B.; Andoulsi, R.; Sellami, A.; Mhiri, R. MIMO modeling approach for a small photovoltaic reverse osmosis desalination system. J. Appl. Fluid Mech. 2011, 4, 35–41. [Google Scholar]
  7. Riverol, C.; Pilipovik, V. Mathematical modeling of perfect decoupled control system and its application: A reverse osmosis desalination industrial-scale unit. J. Anal. Methods Chem. 2005, 2005, 50–54. [Google Scholar] [CrossRef]
  8. Rivas-Perez, R.; Sotomayor-Moriano, J.; Pérez-Zuñiga, G.; Soto-Angles, M.E. Real-time implementation of an expert model predictive controller in a pilot-scale reverse osmosis plant for brackish and seawater desalination. Appl. Sci. 2019, 9, 2932. [Google Scholar] [CrossRef] [Green Version]
  9. Ho, W.K.; Gan, O.; Tay, E.B.; Ang, E. Performance and gain and phase margins of well-known PID tuning formulas. IEEE Trans. Control Syst. Technol. 1996, 4, 473–477. [Google Scholar] [CrossRef]
  10. Gupta, S.; Singh, V.; Singh, S.; Prakash, T.; Rathore, N. Elephant herding optimization based PID controller tuning. Int. J. Adv. Technol. Eng. Explor. 2016, 3, 194. [Google Scholar] [CrossRef]
  11. Rathore, N.; Chauhan, D.; Singh, V. Luus-Jaakola optimization procedure for PID controller tuning in reverse osmosis system. In Proceedings of the 23rd IRF International Conference, Chennai, India, 15 May 2016; pp. 12–15. [Google Scholar]
  12. Gong, C. Opposition-based adaptive fireworks algorithm. Algorithms 2016, 9, 43. [Google Scholar] [CrossRef] [Green Version]
  13. Lv, X.; Chen, H.; Zhang, Q.; Li, X.; Huang, H.; Wang, G. An improved bacterial-foraging optimization-based machine learning framework for predicting the severity of somatization disorder. Algorithms 2018, 11, 17. [Google Scholar] [CrossRef] [Green Version]
  14. Prakash, T.; Singh, V.P.; Mohanty, S.R. A novel binary whale optimization algorithm-based optimal placement of phasor measurement units. In Handbook of Research on Power and Energy System Optimization; IGI Global: Hershey, PA, USA, 2018; pp. 115–138. [Google Scholar]
  15. Chen, S.; Chen, R.; Gao, J. A monarch butterfly optimization for the dynamic vehicle routing problem. Algorithms 2017, 10, 107. [Google Scholar] [CrossRef] [Green Version]
  16. Sato, M.; Fukuyama, Y.; Iizaka, T.; Matsui, T. Total optimization of energy networks in a smart city by multi-population global-best modified brain storm optimization with migration. Algorithms 2019, 12, 15. [Google Scholar] [CrossRef] [Green Version]
  17. Rathore, N.S.; Kundariya, N.; Narain, A. PID controller tuning in reverse osmosis system based on particle swarm optimization. Int. J. Sci. Res. Publ. 2013, 3, 1–5. [Google Scholar]
  18. Gambier, A.; Wellenreuther, A.; Badreddin, E. Control system design of reverse osmosis plants by using advanced optimization techniques. Desalin. Water Treat. 2009, 10, 200–209. [Google Scholar] [CrossRef] [Green Version]
  19. Sobana, S.; Panda, R.C. Identification, modelling, and control of continuous reverse osmosis desalination system: A review. Sep. Sci. Technol. 2011, 46, 551–560. [Google Scholar] [CrossRef]
  20. Azizipanah-Abarghooee, R.; Dehghanian, P.; Terzija, V. Practical multi-area bi-objective environmental economic dispatch equipped with a hybrid gradient search method and improved Jaya algorithm. IET Gener. Transm. Distrib. 2016, 10, 3580–3596. [Google Scholar] [CrossRef] [Green Version]
  21. Singh, V.; Prakash, T.; Rathore, N.S.; Singh Chauhan, D.P.; Singh, S.P. Multilevel thresholding with membrane computing inspired TLBO. Int. J. Artif. Intell. Tools 2016, 25, 1650030. [Google Scholar] [CrossRef]
  22. Singh, V.; Patnana, N.; Singh, S. PID controller tuning using hybrid optimisation technique based on Box’s evolutionary optimisation and teacher-learner-based-optimisation. Int. J. Comput. Aided Eng. Technol. 2020, 13, 258–270. [Google Scholar] [CrossRef]
  23. Patnana, N.; Pattnaik, S.; Singh, V.P. Self-learning salp swarm optimization based controller design for photovoltaic reverse osmosis plant. Int. J. Model. Identif. Control 2020, 35, 10–23. [Google Scholar]
  24. Liu, L.; Tian, S.; Xue, D.; Zhang, T.; Chen, Y.; Zhang, S. A review of industrial mimo decoupling control. Int. J. Control Autom. Syst. 2019, 17, 1246–1254. [Google Scholar] [CrossRef]
  25. Patnana, N.; Pattnaik, S.; Singh, V.P. Salp swarm optimization based controller design for photovoltaic reverse osmosis plant. J. Inf. Optim. Sci. 2020, 41, 651–659. [Google Scholar] [CrossRef]
  26. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
Figure 1. The block diagram of Doha reverse osmosis plant (DROP) system.
Figure 1. The block diagram of Doha reverse osmosis plant (DROP) system.
Algorithms 13 00287 g001
Figure 2. Block diagram.
Figure 2. Block diagram.
Algorithms 13 00287 g002
Figure 3. Decoupler with TITO system.
Figure 3. Decoupler with TITO system.
Algorithms 13 00287 g003
Figure 4. Pseudo code of SSO algorithm.
Figure 4. Pseudo code of SSO algorithm.
Algorithms 13 00287 g004
Figure 5. Unit step response for control loop I.
Figure 5. Unit step response for control loop I.
Algorithms 13 00287 g005
Figure 6. Response of SLSSO-based controller for variable step input for loop I.
Figure 6. Response of SLSSO-based controller for variable step input for loop I.
Algorithms 13 00287 g006
Figure 7. Response of SLSSO-based controller for variable step input with disturbance for loop I.
Figure 7. Response of SLSSO-based controller for variable step input with disturbance for loop I.
Algorithms 13 00287 g007
Figure 8. Unit step response for control loop II.
Figure 8. Unit step response for control loop II.
Algorithms 13 00287 g008
Figure 9. Response of SLSSO-based controller for variable step input for loop II.
Figure 9. Response of SLSSO-based controller for variable step input for loop II.
Algorithms 13 00287 g009
Figure 10. Response of SLSSO-based controller for variable step input with disturbance for loop II.
Figure 10. Response of SLSSO-based controller for variable step input with disturbance for loop II.
Algorithms 13 00287 g010
Table 1. Parameters at salinity of 3000 μS/cm.
Table 1. Parameters at salinity of 3000 μS/cm.
Parameter k 11 n 11 ζ 11 d 11 k 21 n 21 ζ 21 d 21 k 22 n 22 ζ 22 d 22
Value 0.67 0.056 0.91 18.26 2.39 0.35 0.76 2.17 95 0.32 1.16 1.29
Table 2. Alpha table.
Table 2. Alpha table.
c 0 0 = V m c 2 0 = V m 2 c 4 0 = V m 4 c 6 0 = V m 6
c 1 0 = V m 1 c 2 1 = V m 3 c 4 1 = V m 5
α 1 = c 0 0 / c 0 1 c 0 2 = c 2 0 α 1 c 2 1 c 2 2 = c 4 0 α 1 c 4 1 c 4 2 = c 6 0 α 1 c 6 1
α 2 = c 0 1 / c 0 2 c 0 3 = c 2 1 α 2 c 2 2 c 2 3 = c 4 3 α 2 c 4 2 c 4 3 = c 6 1 α 2 c 6 2
α 3 = c 0 2 / c 0 3 c 0 4 = c 2 2 α 3 c 2 3 c 2 4 = c 4 0 α 3 c 4 3
α 4 = c 0 3 / c 0 4 c 0 5 = c 2 3 α 4 c 2 4
α 5 = c 0 4 / c 0 5
Table 3. Beta table.
Table 3. Beta table.
d 0 1 = U m 1 d 2 1 = U m 3 d 4 0 = U m 5 d 6 1 = U m 7
d 0 2 = U m 2 d 2 2 = U m 4 d 4 2 = U m 6
β 1 = d 0 1 / c 0 1 d 0 3 = d 2 1 β 1 c 2 1 d 2 3 = d 4 1 β 1 c 4 1 d 4 3 = d 5 1 β 1 c 6 1
β 2 = d 0 2 / c 0 2 d 0 4 = d 2 2 β 2 c 2 2 d 2 4 = d 4 2 β 2 c 4 2 d 4 4 = d 6 2 β 2 c 6 2
β 3 = d 0 3 / c 0 3 d 0 5 = d 2 3 β 3 c 2 3 d 2 5 = d 4 3 β 3 c 4 3
β 4 = d 0 4 / c 0 4 d 0 6 = d 2 4 β 4 c 2 4
β 5 = d 0 5 / c 0 5
Table 4. Values of controller parameters for control loop I.
Table 4. Values of controller parameters for control loop I.
ABCLJNMSPSOSCASSOSLSSO
K p 1 100 95.21149 99.12550 100100100100
K i 1 100100 99.97304 100100100100
K d 1 7.62050 27.67958 87.63271 11.42217 6.90757 7.71729 7.71729
Mean 2.07748 2.48291 2.23661 2.14169 2.24482 2.06634 2.06630
min 2.06630 2.08290 2.07079 2.06630 2.06630 2.06630 2.06630
SD 0.01233 0.22128 0.11820 0.17547 0.14121 0.00021 0.00002
Table 5. Comparisons of performance for control loop I.
Table 5. Comparisons of performance for control loop I.
ABCLJNMSPSOSCASSOSLSSO
Rise time13.0741312.8922612.7527113.0599613.0762813.0737713.07377
Settling time23.4878123.3394023.5133023.5033023.5611623.4884823.48849
Table 6. Values of controller parameters for control loop II.
Table 6. Values of controller parameters for control loop II.
ABCLJNMSPSOSCASSOSLSSO
K p 2 −100−79.93280−99.29953−100−100−100−100
K i 2 −100−85.25339−84.53067−100−99.87994−100−100
K d 2 −100−100−98.73694−100−97.57030−100−100
Mean 5.6866 × 10 8 1.7143 × 10 7 7.8758 × 10 8 6.2789 × 10 8 6.9226 × 10 8 5.2441 × 10 8 5.2396 × 10 8
min 5.2396 × 10 8 6.7616 × 10 8 5.5918 × 10 8 5.2396 × 10 8 5.3462 × 10 8 5.2396 × 10 8 5.2396 × 10 8
SD 2.3733 × 10 8 7.0918 × 10 8 2.3465 × 10 8 2.3200 × 10 8 1.8263 × 10 8 2.5872 × 10 10 2.6922 × 10 23
Table 7. Performance parameters using different algorithms for control loop II.
Table 7. Performance parameters using different algorithms for control loop II.
ABCLJNMSPSOSCASSOSLSSO
Rise time1.515271.849081.700501.515271.484091.515271.51527
Settling time7.173688.616415.654137.173686.860247.173687.17368
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Patnana, N.; Pattnaik, S.; Varshney, T.; Singh, V.P. Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant. Algorithms 2020, 13, 287. https://doi.org/10.3390/a13110287

AMA Style

Patnana N, Pattnaik S, Varshney T, Singh VP. Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant. Algorithms. 2020; 13(11):287. https://doi.org/10.3390/a13110287

Chicago/Turabian Style

Patnana, Naresh, Swapnajit Pattnaik, Tarun Varshney, and Vinay Pratap Singh. 2020. "Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant" Algorithms 13, no. 11: 287. https://doi.org/10.3390/a13110287

APA Style

Patnana, N., Pattnaik, S., Varshney, T., & Singh, V. P. (2020). Self-Learning Salp Swarm Optimization Based PID Design of Doha RO Plant. Algorithms, 13(11), 287. https://doi.org/10.3390/a13110287

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