Abstract
In this paper, a novel modeling technique has been attempted to develop the mathematical model for a bioreactor functioning at multiple operating regions. The first principle mathematical equations of the reactor are used with the POLYMATH software to generate essential data for the model development. A relative analysis is also carried out with the existing models in the literature. An optimal PID controller is then designed using a multiobjective particle swarm optimization algorithm. The controller tuning procedure is individually discussed for both the stable and unstable steady state regions. The controller tuned for each region is scheduled using a setpoint scheduler to achieve a complete control over the bioreactor. The effectiveness of the proposed scheme has been confirmed through a comparative study with the controller tuning methods proposed in the literature. The results show that, the proposed method provides enhanced performance in effective reference tracking and load disturbance rejection with minimal ISE and IAE. Finally the proposed method is validated on the nonlinear bioreactor model in the presence of a measurement noise. The results testify that the PSO tuned PID performs well in tracking the change in biomass concentration at the entire operating region.
1. Introduction
Bioreactor plays a vital role in chemical process industries to produce important chemical and biochemical compounds. In this system, living organisms also known as microbes are converted into marketable products such as beverages, antibiotics, vaccines, and industrial solvents [1–3]. The quality of the final product from a bioreactor depends mainly on the control loop employed to monitor and control the microbial growth based on the reference input. Apart from this, incidental external and internal disturbances in a reactor may result in reactor failure. Therefore, there is a strong financial inspiration to develop a finest control scheme that would facilitate rapid startup and stabilization of continuous bioreactors subject to redundant disturbances [4].
In the literature, a variety of methods have been discussed to implement a robust controller for the bioreactor operating at single or multiple steadystates. Kumar et al. have examined a bioreactor with input multiplicities. With an experimental study, the mathematical models for the different steady state operating regions are developed, and a nonlinear PI controller was implemented [5]. Sivakumaran et al. have discussed recurrent neural network (RNN) based modeling method for a nonlinear bioreactor operating at single steady state and implemented a nonlinear model predictive controller to obtain satisfactory result [6]. Nagy has proposed a neural network modelbased predictive control (NNMPC) strategy for a fermentation bioreactor model, and with simulation result they concluded that the proposed NNMPC provides enhanced performance compared to linear MPC and PID controller [7]. Giriraj Kumar et al. have discussed a geneticalgorithm (GA) based PI controller tuning for a bioreactor operating at stable steady state [8]. Bioreactor operating at unstable steady state was widely discussed by the researchers due to its complexity and instability [9–11]. A detailed analysis to realize various model based control configurations for continuous bioreactor was presented by Zhao and Skogestad [4].
In the control literature, regardless of the considerable progress in advanced process control proposals such as sliding mode control (SMC), model predictive control (MPC), and internal model control (IMC), PID controllers are still widely employed in industrial control systems because of their structural simplicity, reputation, robust behavior, and easy implementation. Along with the system’s stability, it also satisfies chief performance such as smooth reference tracking, efficient disturbance rejection, and measurement noise attenuation criteria. Most of the PID tuning approaches proposed for stable and unstable system require numerical computations to identify the optimal controller parameters [12, 13].
Recently, soft computing approachbased PID tuning has attracted the researchers due to its ability to find optimized controller parameters with a minimized computation time. The particle swarm optimization (PSO) algorithm is one of the soft computing methods introduced by Kennedy and Eberhart [14]. In recent years, it is widely considered to find optimal solutions for various engineering optimization problems [15–22]. In this paper, PSO algorithm is primarily considered for PID controller parameter optimization because (i) is an autotuning method and does not require detailed mathematical description of the process under control, and (ii) very few parameters to assign compared to other evolutionary methods such as GA and bacterial foraging optimization (BFO) [23]. PSOtuned optimal controller implementation for stable systems are reported by Korani et al. [24]. Chang and Shih proposed an improved PSO optimization algorithm to tune the PID controller for a nonlinear inverted pendulum system [15]. Zamani et al. have discussed a multiobjective PSO algorithm to tune the controller for automatic voltage regulator (AVR) problem [25]. They have proposed a multiobjective PSO algorithm to tune the H_{∞} PID controller for a nonlinear system [26]. Banu and Uma have discussed a hybrid algorithmbased PID controller implementation for a nonlinear CSTR [27]. Kanth and Latha have attempted a PSObased PID and IPD controllertuning for a class of unstable systems [28, 29]. The PSObased controller provides better servo and regulatory responses than the classical PID and modified internal model controller (IMC). They also discussed a PSObased controller implementation for a bioreactor operating at unstable operating region.
From the literature, it is observed that a nonlinear system with multiple steady states can be effectively controlled by a gainscheduled control scheme [27, 30]. This control scheme consists of a family of local controller and a scheduler (switching unit). The switching unit activates the controller based on the setpoint in order to achieve effective control over the entire operating region of the process.
The main contributions of the work are as follows. Initially, the first principle model has been developed using the POLYMATH software generated data. A multiple objective functionbased PSO algorithm has been proposed for the PID controller parameter tuning. PID tuning for the stable and unstable operating region of the bioreactor is separately addressed. The tuned PID controllers are gain scheduled to provide a complete control on the bioreactor system. The effectiveness of the proposed scheme has been validated through a simulation study by using the nonlinear bioreactor model.
2. Process Description
Bioreactor can be defined as a reactor system employed to execute a number of biological reactions in a liquid medium to form intermediate and final products. Figure 1 shows the schematic diagram of the bioreactor. The dynamic behavior of the reactor is complex, and a number of vital manufacturing processes belong to this group: The basic reaction inside the bioreactor is where “” is the reactant, “” is reaction rate constant, and “” is the product.
Biosynthesis is widely utilized to convert the living cells (biomass/microbes) into marketable chemical, pharmaceutical, food and beverage products. Equation (2) shows the operation performed during biosynthesis. In this operation, the biomass/microbes consume nutrients from the substrate (feed) to cultivate and to produce more cells and important products [1]. During this operation the bioreactor is kept under a controlled environment with constant pH, temperature, agitation rate, and dissolved oxygen tension to attain better growth of microbes. In biosynthesis, microorganisms play an essential role in the production of industrial chemicals, enzymes, and antibiotics.
The various stages of microbial growth in a bioreactor during biosynthesis are schematically shown in Figure 2. During the initialization of the process, the living cells (microbes) are placed inside the bioreactor maintained with incubated environment, and the necessary feed (substrate) is gradually applied to the microbes. In the first stage (lag phase), all the living cells are allowed to adjust their performance to the new environmental conditions after the adjustment, they are getting ready to begin reproduction. In the second stage (growth phase) the cells consume nutrients and increase in size. If the environmental conditions are favorable, each cell grows, divides into two which, in turn, grow and divide and the cycle conditions. This process is rapid and the growth rates of cells are proportional to cell concentration and the nutrient. After the growth, they reach a minimum biological space called stationary phase (third stage). At this stage, the cell growth is limited due to the lack of one or more nutrients, buildup of toxic materials during biosynthesis and organic acids generated during the growth phase. Many important fermentation products are produced in stationary space. After the stationary phase, the cells will reach death phase. Death phase is mainly due to the toxic by, products and depletion of nutrient supply. In this, a decrease in live cell concentration occurs [2, 3].
In the above three stages (Figure 2), essential work is performed in growth and stationary phase. Since it is necessary to place a controller in this region to increase productivity the cell mass concentration () and substrate concentration () are the two process state variables available for the controller design. To optimize the modelbased controller parameter, it is necessary to delineate the various stability regions in the bioreactor and to study the effect of substrate () on the biomass (). The main objective for implementing the controller for the bioreactor is to maximize the production and to minimize the waste. The stoichiometry for biomass activity is very complex since it varies with environmental conditions, microorganism, and nutrients in the feed. Due to these reasons, unstructured models are mainly considered for analysis purpose. A number of studies are available in the literature for modelbased control of bioreactor [8–10].
The following mathematical equations can describe a variety of industrial bioreactors [11]: where “” is the specific growth rate, “” the biomass concentration, “” the substrate concentration, and “” and “” are yield parameters for the product. At steady state, the variables will be , , and .
At steady state operating region, (3) will be If more than one steady state occurs, it is identified as trivial steady state and nontrivial steady state:(i)for trivial solution, when , (ii)for nontrivial solution when .
The nominal parameter and constant values considered in the mathematical equations are presented in Table 1.
3. First Principle Model Development
Figure 3 shows a generalized procedure to be followed for model development and model validation practice for nonlinear bioreactor.
Real time model development for a bioreactor is a challenging job, and sometimes it may not provide a satisfactory model due to its nonlinear nature. The operating time of the bioreactor is also large compared to the other nonlinear chemical process loops existing in the chemical industries. Generally in industries, the bioreactor may run for several days in order to convert the raw material in to products. Collecting the real time data to develop the model around the operating regions (lag, growth, and stationary phases) is time consuming. Hence, in this paper, a first principlebased model development is proposed for the bioreactor widely considered in the literature. The first principle modeling equations from (3) and (4) are considered in this work to develop the mathematical model of the bioreactor (statespace and transfer function model).
The statespace model the of the system is represented by the following equations [11]: where (where: , , , , , , and are the final steady state values of the bioreactor considered during the model development.
Table 2 shows the POLYMATH software codes developed to simulate the bioreactor data at various operating levels. Each program is separately simulated for a simulation time of 10 hours, and the corresponding plots are shown in Figures 4 to 6.
Initially the “lag phase” is simulated by considering the values as represented in Table 2. The generated data from the above program is considered to construct Figure 4, and it depicts the timerelated values of feed and product concentrations in the lag phase region. In this region, the microbes are allowed to adjust their performance to the new environmental conditions. Since the increase in the product concentration in this region is very minimal or approximately zero, in the simulation study the final values for “biomass concentration” and “substrate concentrations” are “0” and “3.850639,” respectively. From Table 1, it is observed that data attained from the simulation study is very close with the actual data (i.e., biomass concentration ; substrate concentration ).
Simulation study for growth phase is done with the values provided in Table 2. The created data from simulation is considered to construct Figure 5. It shows the timerelated values of feed and product concentrations in the growth phase region.
In this region, the microbes consume nutrients from the feed and increase its size. After reaching the maximum size, it will split into two and grows. This operation is rapid since there will be a sudden rise in the product concentration. In the simulation study the final values for “” and “” are “0.989742” and “1.528277,” respectively. From Table 1, it is observed that data attained from the simulation study is very close with the actual data (i.e., ; ).
Simulation study for stationary phase is conducted using Table 2 values, and the generated data provides Figure 6. This figure depicts the timerelated values of feed and product concentrations in the stationary phase region.
In this stage, the cell growth is limited due to the lack of the nutrients, buildup of toxic materials during biosynthesis, and organic acids generated during the growth phase. During this stage, all the available microbes are converted in to useful chemical and biochemical products. In this study the final values for “” and “” are “1.526381” and “0.159153,” respectively. From Table 1, it is observed that data attained from the simulation study is very close with the actual data (i.e., ; ).
Figure 7 represents the combinations of the product concentration achieved in lag, growth, and stationary phases. This diagram is the exact replica of the cell growth stage discussed in Section 2 (Figure 2).
From the POLYMATH simulation, the following values are obtained for the final feed and product concentrations:(i)lag phase: g/lit; g/lit,(ii)growth phase: g/lit; g/lit,(iii)stationary phase: g/lit; g/lit.
These values are implemented in (7), to develop the first principle model of the bioreactor for the growth and stationary phases, and the developed models are tabulated in Table 3. From the Table 3 values, it is confirmed that the developed model by this procedure is approximately similar to the model existing in the literature [11].
The developed mathematical models shown in Table 3 are utilized along with a delay time of “0.1” in the proposed study (the delay time = measurement of the product concentration, converting the measured concentration into the value acceptable by the control loop) for the PSObased PID controller tuning.
The POLYMATH program considered in this simulation study also provides the following information about the bioreactor considered in this study.
Figure 8 shows the specific growth rate of microbes with respect to the substrate (feed) concentration. Each microbe (cell) has its own size, and after reaching the maximum growth, the microbe starts the reproduction. Later it consumes the feed and continues the reproduction until it reache the elimination phase.
Figure 9 shows the relationship between the specific growth rate and the biomass concentration. Initially the growth rate rapidly increases due to the transition of the microbe from the lag phase to the growth phase. After reaching the maximum growth rate, due to the reproduction operation, the biomass concentration increases.
Figure 10 depicts the relationship between the feed and the biomass concentration. Initially the biomass concentration linearly increases with respect to the feed concentration. After reaching the maximum growth rate in microbes, even though the feed is less, due to the reproduction operations, the concentration in the biomass increases.
Figure 11 shows the specific growth rate of microbes with respect to time. In this a simulation period of 10 hrs is shown. In growth phase, the rate of growth linearly increases with respect to time. When it attains the maximum growth rate, the microbes researches their saturation level. Following the saturation, some of the microbe (cell) may expire due to the reasons such as lack of food ageing. Hence, after saturation, the growth rate of microbe rapidly decreases with respect to time.
Figure 12 shows the relationship between the specific growth rate and the “pH” value inside the bioreactor. Form the diagram, it is observed that the microbes can effectively grow when the surrounding “pH” is 6 < pH < 7.5.
4. PSO Algorithm
Particle swarm optimization (PSO) algorithm is a populationbased evolutionary computation technique developed by the inspiration of the social behavior in bird flocking or fish schooling [14]. It has become one of the most powerful soft computing methods for solving optimization problems. It attempts to mimic the natural process of group communication of individual knowledge, to achieve some optimum property. In PSO algorithm, a population of swarm is initialized to move in a “” dimensional search space. Each particle in swarm has a position represented by a vector “ = () ” and velocity represented by a vector “ = () ”. At the beginning, each particle in the swarm population is scattered randomly throughout the entire search space “” and with the guidance of the performance criterion, the flying particles dynamically regulate their velocity and position according to their own flying experience and their companions flying experience. Each particle remembers its best position obtained so far, which is denoted (). It also receives the globally best position achieved by any particle in the population, which is denoted as () [17, 20, 21, 24]. The updated velocity of each particle can be calculated using the present velocity and the distances from and . The updated velocity and the position are given in (6) and (8), respectively. Equation (7) shows the inertia weight: where and are positive constants known as acceleration constants. “” is the cognitive learning rate and “” is the global learning rate. and are random numbers in the range 01. The parameter “” is inertia weight that increases the overall performance of PSO. The larger value of “” (i.e., ) can favor the global widerange search, and lower value of “” (i.e., ) implies a higher ability for local nearby search (see Pseudocode 1).

5. PSObased PID Tuning
The PID tuning process is to find that the optimal values for , , and form the search space that minimizes the objective function (11). During this search, the performance criterion “” guides the optimization algorithm to get appropriate value for the controller parameters.
In the literature, there is no clear guide line to assign the algorithm parameters for the PSO algorithm. In this study, before proceeding with the PSObased controller tuning, we assigned the following parameters.
Dimension of search space “” is three (i.e., , , ); the number of swarm and bird steps is considered as twelve; the cognitive () and global () search parameter is assigned the value of 2, the inertia weight “” and “” is assigned the value 0.2 and 1, respectively.
In this study, a noninteracting form of parallel PID controller is considered to control the nonlinear system.
Parallel PID structure where: , .
5.1. Controller Tuning
The controller tuning process is employed to find the best possible values for , , and . In order to achieve the superior accuracy during the optimization search, it is necessary to assign appropriate performance index which guides the PSO algorithm.
In this work we considered the following performance criterion (11) with six parameters, ISE, IAE, , , , and : where where isrise time (time required for to reach 100% of its setpoint at the first instant), issettling time – time required for to reach an stay at [i.e.,], issteady stateerror, is simulation time, where are weighting functions used to set the priority of the multiobjective performance index parameters, and the value of “” varies from 0 to 10. The performance criterion guides the PFO algorithm to get appropriate values for the controller parameters.
Figure 13 shows the block diagram of the MOPSObased PID controller tuning proposed for the bioreactor model.
5.2. Optimization Search
Prior to the optimization search, it is necessary to assign the parameters for the PSO algorithm and the multiobjective performance index criterion.
In this study, the following values are assigned:(i)dimension of the search space (i.e., , , ),(ii)number of swarm and bird steps is considered as twelve, the cognitive () and global () search parameter is assigned the value of 2, the inertia weight “” and “” is assigned the value 0.2 and 1, respectively, (iii)the maximum iteration for generation (Iter) is set to 250,(iv)boundaries for the three dimensional search space is assigned as (v)the weighting values are assigned as and ,(vi) maximum simulation time is selected as 100 sec,(vii)the “” is preferred as <25% of the maximum simulation time. The simulation time should be selected based on the process time delay,(viii)the overshoot () range is selected as <100% of the reference signal,(ix) the “” is preferred as <50% of the maximum simulation time,(x)the is assigned as zero,(xi) the reference signal is considered as unity (i.e., ).
5.3. Comparative Study
In order to evaluate the performance of the proposed algorithm, a comparative analysis is done with the most successful soft computing methods like PSO and BFO.
PSO
The simulation is carried out by using the PSO algorithm attempted by RajniKanth and Latha [28]. The following algorithm parameters are considered: dimension of search space is three (i.e., , , ); number of swarm and bird steps is considered as 25; the cognitive () and global () search parameter is assigned the value of 2 and 1.5, respectively. The inertia weight “” is fixed as 0.7.
BFO
For the basic BFO algorithm, the following values are considered: dimension of search space is three; number of bacteria is chosen as ten; the number of chemotaxis steps is set to five; number of reproduction steps and length of a swim is considered as four; number of eliminationdispersal events is two; number of bacteria reproduction is assigned as five; probability for eliminationdispersal has a value of 0.2 [29].
6. Results and Discussions
6.1. SetPoint Scheduling
A complete control scheme of the nonlinear bioreactor operating at multiple steady states is schematically presented in Figure 14. This system has a family of classical PID controller (PID1 for unstable region and PID2 for stable region) and a gain scheduler switch. The switch considered in this study is an electronic switch, and the switching time by this device is negligible.
In the bioreactor, during the lag phase, a manual control is employed to monitor and control the biomass activity (0 to 10 hr in Figure 7). During the growth phase, the cell growth is very rapid, and its concentration is the function of the time and the substrate feed rate (from 10 to 20 hr in Figure 7). At this region, in order to increase the production rate and to reduce the wastage, it is essential to implement an optimally tuned controller. All the necessary products are produced in the stationary phase (from 20 to 30 hr in Figure 7).
In order to maintain the product quality, it is necessary to implement another controller in this region since the PID tuned for growth phase (unstable state) will not provide the robust performance in the stationary phase (stable steady state). This can be avoided by implementing a gain scheduler along with a separate PID controller based on the operating regions.
The setpointbased gain scheduler considered in this work is a comparator assisted electronic switch, which activates the corresponding PID controller based on the available setpoint value. The scheduler unit continuously monitors the reference signal (biomass concentration) applied to the bioreactor and selects the suitable local controller unit to execute a complete control over the entire operating region as depicted in Figure 14.
6.2. Growth Phase Model
The secondorder unstable transfer function model of the bioreactor with a measurement delay of “0.1” is shown in (12):
This process model has two stable poles and an unstable pole. The PSObased PID controller tuning is proposed for the model with the multiobjective performance criterion represented in (11). To examine the performance of the MOPSO algorithm, five trials are performed. The best solution among the trials is chosen as the optimal solution. The final convergence of the optimized controller parameters is shown in Figure 15, and the , , values are summarized in Table 3.
It shows that the proposed MOPSObased tuning has less number of iteration compared to BFO algorithm [29]. Figure 16 depicts the servo response of the growth phase model. From this response, it is observed that, the proposed method provides a reduced value of “” and “” compared to PSO and BFO tuned PID controller. Through a simulation time of 100 min, the MOPSO algorithm provides an enhanced result in reference tracking with reduced ISE and IAE than a single objective PSO algorithm (Table 4) [28].
The regulatory response is studied with a load disturbance of 0.3 (30% of setpoint) introduced at 75th min as in Figure 17. The controller successfully rejects disturbance at 87th min and allows the system to track the setpoint from 88th min onwards. The controller also provides a smooth output. From the result, the observation is that the proposed PID controller provides a very stable and smooth response for both the setpoint tracking and load disturbance elimination.
6.3. Stationary Phase Model
The secondorder stable transfer function model of the bioreactor with a measurement delay of 0.1 hr is depicted in (13):
The MOPSObased PID controller tuning is proposed for the above model, and the ultimate convergence of the controller parameters (best solution among the five trials) is shown in Figure 18, and the optimized values are presented in Table 5.
A step change of 1 g/lit is introduced to the bioreactor model, and the corresponding reference tracking performance of the controller (PID2) is presented in Figure 19. A relative analysis with the previously published work is also carried out [8]. The proposed controller provides improved performance with reduced values for , , ISE, and IAE compared to other methods such as ZN, GA, and Skogestad.
A load disturbance of 0.3 (30% of setpoint) is introduced at 75th min to test the regulatory performance of the controller. From Figure 20, it can be inferred that the PSOtuned PID controller is able to reject the load disturbance quickly (77th min) and maintains the biomass concentration based on the given setpoint. A summary of the performance comparisons between the proposed and the existing methods are presented in Table 4. The MOPSOtuned PID provides superior performance with reduced ISE and IAE values in entire operating region. (In Table 5 the suffix “” is for setpoint and “” is the load disturbance.)
6.4. Controller Implementation on Bioreactor Model
The performance of the proposed method is examined on the nonlinear bioreactor model, developed using the first principle equations from (3) to (5). Simulation studies have been carried out on the nonlinear model to demonstrate the reference tracking and disturbance rejection performance of the proposed MOPSOtuned PID controller in the growth phase and stationary phase of the biomass.
Initially, a step change in the biomass concentration of 0.9951 g/lit has been introduced in the system. To study the regulatory performance, a 20% change in the dilution rate is applied at 300th sampling instant, and the corresponding change in system parameters such as substrate concentration, dilution rate, and controller output (PID1) is observed (Figures 21 and 22). The observation is that change in substrate concentration is more compared to other variables. Later a step change of 1.5302 g/lit has been applied at 700th sampling instant to study the performance in the stationary phase. This change in setpoint, activates the local PID controller (PID 2) with the help of the setpoint scheduler. From Figure 21, the observation is that though there is a small overshoot due to the process nonlinearity, the proposed controller provides a smooth setpoint tracking in the entire operating region. The regulatory performance is then studied by introducing 20% change in the dilution rate at 1000th sampling instant. The corresponding change in the substrate concentration and the dilution rate is described in Figure 22. Due to the presence disturbance, the dilution rate and the substrate concentration are more oscillatory from 1000th to 1150th sampling instant. Later the response reaches a very smooth value up to 1400th sampling instant.
The transition between the controller (from PID1 to PID2) is initiated by the setpoint scheduler during the 700th sampling instant (Figure 21). Due to this changeover, a spiky change in the controller output is observed at 700th sampling instant. This effect disappears quickly and the controller output reaches a steadystate value at 783rd sampling instant. The robustness of the proposed control scheme is then tested by introducing a measurement noise (noise power of 0.001; sampling time of 0.1 sec).
Figure 23 shows the variations of biomass concentration and Figures 24 and 25 depict the variation of substrate concentration, dilution rate, and controller signal in the presence of measurement noise. The reference tracking response of the setpoint scheduled bioreactor shows that the proposed scheme works well in the noisy environment. In the presence of measurement noise, the dilution rate and the substrate concentration are more oscillatory from 0th to 500th sampling instant. After 500th sampling instant, the system enters into stationaty phase since the effect of measurement noise is minimal. From these observations, it is concluded that the proposed controller scheme is robust, and it helps to provide a better output for setpoint tracking in the entire operating region of the bioreactor.
7. Conclusion
Modeling and the modelbased robust controller implementation for the nonlinear chemical system is a challenging work. This paper presents a novel method for modeling and modelbased PID controller implementation for a nonlinear bioreactor. First principle equationbased model for the bioreactor is developed using POLYMATH software. The developed model is then considered for PSO algorithmbased PID controller tuning. The proposed method implements a time domain associated multiobjective performance criterion to guide the PSO algorithm in order to obtain optimized PID parameters. Two different local PID controllers are tuned separately for both the stable and unstable operating regions and implemented using a setpointassisted gain scheduler to track the biomass concentration based on the reference signal. From the extensive simulation studies, it can be observed that the proposed scheme provides enhanced result in setpoint tracking and disturbance rejection. This controller also shows a robust performance in the presence of measurement noise. From the study, it can be concluded that the proposed method helps to provide an approximated mathematical model using the modelling equations governing the systems. The proposed controller scheme can be considered as an alternative to other gain scheduled controllers widely employed in nonlinear system control.
Nomenclature:
Positive constants (0–2)  
:  Feed concentration 
:  Reactor concentration 
:  Error 
:  Global best position 
IAE:  Integral absolute error 
ISE:  Integral squared error 
:  Proportional gain 
:  Integral gain 
:  Derivative gain 
:  Filter constant (10) 
:  Local best position 
:  Measured output 
:  Inlet flow rate 
:  Random number (01) 
:  Reference input 
:  Position of particle 
:  Setpoint 
:  Integral time constant 
:  Derivative time constant 
:  PID controller output 
:  Velocity of particle 
:  Inertia weight of particle 
:  Measured variable 
:  Process output. 
:  Iteration number 
:  Updated iteration number. 