Next Article in Journal
Amino Ethers of Ortho-Phosphoric Acid as Extragents for Ethanol Dehydration
Previous Article in Journal
Application of Gold Nanoparticle-Based Materials in Cancer Therapy and Diagnostics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methanol–Water Purification Control Using Multi-Loop PI Controllers Based on Linear Set Point and Disturbance Models

Department of Chemical Engineering, Faculty of Engineering, Universitas Indonesia, Depok 1642, Indonesia
*
Author to whom correspondence should be addressed.
ChemEngineering 2021, 5(4), 70; https://doi.org/10.3390/chemengineering5040070
Submission received: 9 September 2021 / Revised: 9 October 2021 / Accepted: 12 October 2021 / Published: 20 October 2021

Abstract

:
Dimethyl ether (DME), derived from methanol, has been issued as an alternative to diesel and LPG. This study will continue the research of an optimized process control design for DME purification plants, specifically the methanol purification process. This way, water can be separated, and the main product, methanol, can be recycled for DME synthesis. A setpoint-based linear model of a methanol–water purification system has been developed with four controlled variables (CV) with objectives to; maintain a separated liquid water stream in the bottom stage (Stage 30) of the distillation column for methanol–water separation at 11.22%, keep the purified liquid methanol at a condenser at 58.81 °C and 49.96% of level, and the last CV is the cooler in the distillate stream, to keep the purified methanol’s top product at 40.75 °C. To complete the model, a first order plus dead time (FOPDT) disturbance model is created against the inlet temperature and flow rate of the feed, the major cause of disturbances in the industry. Using a traditional proportional integral controller connected to each controlled variable, a multi-loop control system is formed with optimization and compared to the disturbance rejection of multivariable model predictive control (MMPC). The final improvement against the feed temperature and flow for CV1, CV2, CV3, and CV4 is shown by, respectively, Integral Absolute Error (IAE) values of 79.49%, 99.90%, 100%, and 99.99% and Integral Square Error (ISE) values of 97.1%, 100%, 99.99%, and 100%.

Graphical Abstract

1. Introduction

In 2018, the energy consumption in Indonesia reached a total of 936.33 trillion barrels of oil equivalent (BOE), with 64.46 million BOE for LPG and 450.78 million BOE for fuel [1]. Energy demands and supply projections for 2019–2050 target induction stoves and DME utilization for liquefied petroleum gas (LPG) substitution. DME features easy volatility, no toxicity, no carcinogenicity, good water and alcohol solubility, and no pollution it is miscible with many resins and solvents. At ambient conditions, DME, similar to LPG, can be liquefied at moderate pressure. With a high cetane number, DME is clean-burning and sulfur-free, with extremely low particulates [2,3]. DME is derived from methanol from coal, natural gas, organic waste, and biomass [4]. Companies such as Topsoe, Mitsubishi Co., and Total have recently focused on promoting DME as a new and renewable fuel that can replace LPG and diesel. Solichin et al., (2011) [5] designed a DME purification plant suited and adjusted for Indonesia. In addition to the close market share, raw materials and auxiliary materials can be obtained in Indonesia, resulting in a reduced dependency on imported products. The manufacturing process starts with synthetic gas (CO2 and H2O) that is further synthesized into methanol and dehydrated into DME. The conversion from methanol to DME will still produce excess methanol and water, so they would be separated in a distillation column where the purified methanol is recycled to produce more DME [6,7,8]. This is different from what has been studied by Dyrda et al. [9], who used a membrane-based separation process to separate methanol and water to improve the performance of the fuel cell. In this research, a distillation column is used for separating methanol and water.
What is left is the optimization of this plant design. Process control is the integration of a mathematical model into an operating program [10] to optimize manufacturing, and its quality results in an improvement in safety standard profits [11]. The design of a system controller for the purification process of DME and methanol has been studied several times. A performance evaluation of the traditional proportional integral (PI) controller on DME and methanol purification was conducted by Wahid and Gunawan (2015) [12]. Wahid and Brillianto (2020) then developed a multiple input–multiple output (MIMO) process model with a single controller as the central brain (centralized controller) for the purification of methanol (methanol–water separation) in a DME plant [13]. The MIMO model created is based on the interaction of a controlled variable and how the changes to its manipulated variable affect other controlled variables. Each interaction is described with a linear first order plus dead time (FOPDT) model because a linear model offers simplicity and a reliable solution for system optimization.
This study aims to enhance the centralized control of the methanol purification process into a decentralized (multi-loop) control, in which each of manipulated variables has its controller, creating a multi-loop scheme rather than a single direct control of a multivariable process. A decentralized controller produces a higher result of stability from its faster and more accurate response compared to a centralized controller [14].
Multivariable control has been studied for decades, since nearly all processes would involve a simultaneous control on the production rate (flow), inventory (level and pressure), and process environment (temperature). There are three advantages of using multiple single-loop controllers. The first is that simple algorithms are easy to use, which helps when calculations are implemented with analog computing equipment, the second is that it is easy for plant operating personnel to understand the control structure (they only need to monitor the controlled variables and adjust its related manipulated variable), and the last is that the standard control designs are developed for the common unit. Although analysis and experience are required, multiloop designs will continue to be used extensively [11]. For this reason, the performance of multi-loop PI controllers and the tuning method of Big Log-Modulus Tuning has been studied and endorsed by Vu & Lee (2007), Biyanto et al., (2018), and M. Farsi (2015) [15,16,17].
Many process control designs and implementations can benefit from linear transfer function models created using empirical approaches. FOPDT models are suitable for process control analysis and design with dead time. The linear regression identification approach is more generic and may be used in industrial processes [11]. The examples of a linear model on common units are the 2 × 2 MIMO model by Wood-Berry, Vinante-Luyben, Wardle-Wood, and Ogunnaike-Ray for a distillation column. They utilize two controlled variables and two manipulated variables. Ogunnaike-Ray also created a 3 × 3 MIMO model, and a 4 × 4 MIMO model was created by Alatiqi. All of these were studied further by Huang et al. (2003) for their simplicity and implementation [18]. The Wood-Berry distillation column was then proposed using the fractional order PID (FOPID) controller and proved to be successful in terms of performance [19]. Other research conducted by Kalpana et al. (2017), Memon and Kalhoro (2021), and Alawad and Muter (2021) has proven the success of a multivariable PID controller [20,21,22]. However, they used an existing model. Here, this study develops a disturbance model, as was done by Ogunnaike and Ray (1979) on the Wood-Berry model (1973), which was originally only the SP model [23].
Since we already have the setpoint (SP) models by Wahid and Brillianto (2020) [13], the disturbance model will be created using the same method as the FOPDT model. The used FOPDT method ensures that the actual process is best represented. To make the controller more appropriate for factory conditions in real time, we added manual disturbance and adjusted the controller to be more responsive to the disturbance.
Dietrich et al.’s (2020) study about a DME plant for power concluded that the feed composition and temperature have a significant impact on the composition of the product, and this is because they are inconsistent and prone to be unstable after derivation from syntactic gas [24]. Thus, these two variable changes will be observed with the controlled variables of the process and then described with the FOPDT model. With a MIMO process model, a disturbance model, and an efficient PI control for each manipulated variable, a comprehensive and equipped multi-loop PI controller is created for the methanol purification process in the DME production plant.

2. Methods

2.1. Data Collection

The simulation of purification is made separately to reduce the variables that need to be controlled. The separation of methanol from water is based on the difference in boiling points: 64.7 °C and 100 °C, respectively. The data used were collected from Solichin et al.’s (2011) [5] studies on DME factory design. Table 1 and Table 2 below are the parameters for each unit.
Figure 1 and Table 3 show the simulation process and a list of controlled variables, manipulated variables, their set points, and the disturbance.

2.2. Creating the Disturbance Model through System Identification

System identification is then conducted to create a dynamic system-based mathematical model based on the input and output signals of the system called a transfer function. Empirical models are designed completely using observational empirical evidence, data input, and data output from a dynamic system. A linear empirical dynamic model of a process is calculated using experimental data, one of which is to use step-response data. Interpretations may use Marlin’s (2000) process reaction curve (PRC) method [11]. Parameters of a PRC include the dead time (ϴ), the time constant (𝞃), the magnitude of the input shift (δ), the magnitude of the steady-state change in the output (Δ), and the overall slope (S). When measurement noise is considered, the model parameters can be precisely calculated. Deriving from the tangent line, the second method of the PRC by Smith (1972) resembles the closest actual process response, as pictured in Figure 2 below [24].
Thus, the parameters resulted will be used to create a transfer function of the first order plus dead time (FOPDT) method (Equation (4)).
Kp = δ
      τ = S = 1.5 t 63 % t 28 %
θ = t 63 % τ
FOPDT G s = Kp   e θ s τ s + 1

2.3. Controller Pairing

The MIMO (multi-input–multi-output) method is a process control system that can manage process variables that communicate with a single controller at the same time [25]. A decentralized control system is a multivariable system where n inputs and n output variables are treated as n mono-variable systems. A decentralized control structure has certain advantages, such as design and hardware simplicity and ease of use [26].
The motive for using a MIMO proportional integral controller for multi-loop control is its trouble-free configuration, straightforward tuning, and its ability to attain most of the estimated control objectives [27]. In multi-loop control, each manipulated variables (n) depend on a single controller (Gc) (Figure 3).
Equation (5) below is the purification 4 × 4 MIMO process model by Wahid and Brillianto (2020). The controlled variables include the condenser temperature (𝑇cond), the cooler output temperature (𝑇cooler), the condenser level (𝐿cond), and the column level (𝐿column). The manipulated variables are the condenser load (𝑄cond), the cooler load (𝑄cooler), the distillate flow rate (𝐹D), and the lower product flow rate (𝐹B) [13].
T c o n d T c o o l e r L c o n d L c o l u m n = 0.2 e 0.68 s 0.15 s + 1 0.012 e 7.4 s 1.12 s + 1 0.584 e 0.05 s 0.03 s + 1 0.0006 e 0.57 s 1.32 s + 1     0.008 e 6.29 s 2.4 s + 1     0.502 e 0.01 s 0.33 s + 1 0.004 e 0.46 s 0.75 s + 1     0.0002 e 1.001 s 0.27 s + 1   0.001 e 1.17 s 0.75 s + 1           0.03 e 0.05 s 0.66 s + 1           0.0003 e 0.35 s 0.75 s + 1             0.000067 e 0.15 s 0.27 s + 1           0.0732 e 0.35 s 0.15 s + 1 0.049 e 1.21 s 1.32 s + 1 0.007 e 0.42 s 1.635 s + 1 0.32 e 0.02 s 0.76 s + 1 Q c o n d Q c o o l e r F D F B
Control pairings themselves are chosen based on the process interactions, dynamic responses, and sensitivity to disturbance [14]. To find which control pairing works best, a relative gain array (RGA) (λ) is one way to decide (Equation (5)). Table 4 and Table 5 explain the interaction between the variables from their RGA.
G p = K P 11 K P 1 n K P m 1 K P m n          Λ = λ 11 λ 1 n λ m 1 λ m n    Λ = G P × G p T
where
Gp = process model; m = controlled variable;
Kp = process gain; n = manipulated variable;
Λ = RGA value.

2.4. Controller Tuning

Process interactions may include undesirable interactions between two or more control loops. The special characteristics of the proportional integral controller, however, are as follows:
-
There are changes in output if the error is not zero. Because of its design, even in small error situations, this controller will eradicate errors.
-
The reset time (𝜏i) ensures that the output returns to the set stage.
PI controllers are ineffective when processes with loads shift very rapidly and processes have a substantial lag between corrective actions and the resulting actions. However, since the DME plant does not require significant changes, the PI controller (Equation (6)) is used. Major load changes and significant variations at the setpoint can be properly controlled by the PI without prolonged oscillation or a permanent offset. After the disruption has occurred, it is easily returned to its proper state.
G c s = K c 1 + 1 τ i s
Controller tuning is used for achieving closed-loop stability. A way to tune decentralized PI controllers is through the biggest long-modulus tuning (BLT) method. This method idea is to tune individual interacting loops and then to de-tune all the loops together with a common factor (F) such that the overall system is stable and gives acceptable load disturbance rejection responses. The BLT method steps are as follows:
  • Calculating initial PI parameter values;
  • BLT de-tuning:
    • Assume factor F with a typical range of 1.5 < F < 4;
    • Calculate new values of the controller’s parameter, where the detune controller gain is Kci, and the detune reset time is 𝛕ii:
      K c i = K c , W G F   and   τ i i = τ i , W G × F ;
    • Determine the log modulus function (Lc) using a multivariable Nyquist plot (W):
      W = 1 + det   1 + G P G C
      G C s = K c 11 1 + 1 τ i i 11 s K c 1 n 1 + 1 τ i i 1 n s K c m 1 1 + 1 τ i i m 1 s K c m n 1 + 1 τ i i m n s ;
    • Reassume the F factor until the bode plot of W is stable, indicators are stated in Table 6 below.

2.5. Performance Evaluation

For each process system, the sample to be taken and evaluated in this study is the CV graph monitored through setpoint tracking and disturbance rejection until the stable condition is reached again. Analysis was conducted using the performance parameters of the Integral Absolute Error (IAE) and the Integral Square Error (ISE).
The IAE is the absolute area of the difference between the area of the setpoint graph and the area of the controlled variable (CV) response graph. This makes the IAE easy to analyze. Meanwhile, the ISE value can be calculated by first squaring the existing data. IAE is better at detecting minor errors. This is because, when the value is smaller than 1, the ISE value will be smaller. IAE also does not have an error with negative values. For long-duration transients, ITAE (Integral of Time-Weighted Absolute Error) would perform best to penalize their deviations. However, the IAE and ISE perform better to correct errors precisely and with a fast response [11]. Since both hold time and deviation accountable, they are enough to underline the performance of a controller. The IAE and ISE can be used independently as a performance indicator of process control because the best IAE or ISE holds the best score [29]. Below are their equations.
I A E = S P t C V t d t
I S E = 0 S P t C V t 2 d t

3. Results and Discussion

3.1. Disturbance Model

The model is obtained using the first order plus dead time (FOPDT) approach and the step response method. CV and MV changes during response changes will produce a PRC, from which we derive a transfer function that will be used to tune the controller. The disturbances that have been identified in this system are the instability of temperature and flow rate in the distillation column. Therefore, the PRC that is made is the influence of the feed-in temperature and flow rate on the liquid level and temperature in the distillation column and its condenser.
The feed-in molar flow and the temperature are increased by 5% on a separate simulation from 659.8 gmol/s to 692.78 gmol/s and from 81.74 °C to 85.827 °C. The disturbance transfer function model derived from the PRC in Figure 4 is concluded in Table 7. The excess flow causes a rise in the T-103 Condenser’s temperature and level. The higher feed-in temperature must have caused the mixture to expand more, causing a lower level in Stage 30 due to a faster condensing and an increase in condenser level. The decreasing slope results in a negative process gain value.
Thus, the final disturbance model (Gd) for methanol–water separation in DME synthesis, with FF as the feed flow rate and TF as the feed temperature, is as shown below.
G d s =   0.0001 e 1.15 s 1.00 s + 1 0.0065 e 8.32 s 3.63 s + 1 0.00003 e 1.72 s 1.38 s + 1 0.0016 e 7.95 s 3.00 s + 1       0.1066 e 2.48 s 1.63 s + 1     0.0902 e 2.48 s 1.50 s + 1 0.0148 e 1.94 s 1.32 s + 1 0.0275 e 9.45 s 3.00 s + 1   F F s T F s

3.2. Controller Pairing

The pairing results from the methanol purification 4 × 4 MIMO process model by Wahid and Brillianto (2020) [13]. RGA analysis is CV2-MV2 and CV4-MV4, as shown as bold below, from the Matlab coding program.
RGA11 =  RGA21 =  RGA31 =  RGA41 =  
0.2452  15.8404  265.6956  0.1637  
RGA12 =  RGA22 =  RGA32 =  RGA42 =  
0.0555  0.6979  −4.7820  0.1538  
RGA13 =  RGA23 =  RGA33 =  RGA43 =  
−0.7833  −50.8558  −851.0568  0.2304  
RGA14 =  RGA24 =  RGA34 =  RGA44 =  
−0.0915  −6.7712  −111.6737  1.0274  
CV1 was paired with MV1 because it has the closest value to 1. To conclude the decentralized cycle of this controller, CV3 and MV3 are treated as a pair because the manipulated variable is the desired actuator position of the T-103 Condenser, and the T-103 Condenser level is CV3, which shows correlations more than other MV available.

3.3. Multiloop PI Controller Tuning

The process of the methanol–water 4 × 4 MIMO process model with BLT was conducted. F is taken as the one closest to -8 dB because we use 4 × 4 matrices, so the stability should be 2 × 4, which is 8 dB ( L c m m a x ). With less encirclement than the others, stability is most shown by F = 1.5. This is proven by the positive magnitude and the Nyquist diagram that does not encircle (−1,0). The fine-tuning method is done by using the ‘Tune’ feature on the PID block in Simulink, MATLAB. This feature will automatically design a controller that has a balance of response and stability. Fine-tuning is done until the most optimal result is obtained and the closed loop is a table. The performance of tuning parameters published by Wahid and Gunawan (2015) [12] will also be tested and compared. Table 8 below are the final value of the controllers parameters.

3.4. Performance Test

Multi-loop controller performance was compared with multivariable model predictive control. The test was conducted through a setpoint change and the manual adding of interference. The test was done by increasing 5% of the initial value individually. For example, when the T-103 Condenser temperature (CV1) setpoint change was tested, the final value of delta (Δ) was input into the step block, whereas the rest of the step block was inactive. Table 9 below are the list of delta (Δ) inputted in the step block.
Below is the performance test results of disturbance rejection held in Simulink. Figure 5 displays the stability for all controllers. As shown on the left, the manual interference of temperature is added by increasing its original state from 81.74 °C to 85.827 °C. All controlled variables demonstrate final stability at 0. The condenser temperature (CV1) decreases when the feed-in temperature increases. It is suspected that this is caused by a slower form of top product that may result in an increase in the product temperature as well (CV2). The fine-tuned proportional-integral controller reveals a faster stabilized graphic and smaller disturbance peak than the BLT ones, but they are similar to WG. For BLT, the CV1 temperature disturbance peak reaches almost −0.32 and −0.2 for CV2, meaning the condenser temperature once reaches 58.49 °C and the MeRecy temperature decreases to 40.55 °C. They stabilize at times 150 and 50 with values back to 58.81 °C and 40.75%. In the fine-tuned controller, the CV1 peak only reaches −0.35, and CV2 and CV4 barely leave 0. All stabilize before time 100.
However, when compared with MMPC, CV1 stabilizes at a 58.167 °C condenser temperature, a 42.057 °C methanol final temperature, a 49.854% condenser level, and a 14.629% T-103 liquid volume at Stage 30. Error analysis will be further discussed in the next section.
In the middle row, the feed flow disturbance is increased from the initial value of 659.8 gmol/s to 692.72 gmol/s. Here, the fine-tuned proportional-integral controller shows the lowest oscillation. As a response towards feed molar flow increases, liquid volume in the bottom stage decreases due to a more fluid increase in heat, causing slower stage movement. The condenser temperature increases because less mixture is condensed, affecting the heat exchange process and resulting in a higher product temperature.
For MMPC results using the unisim software, the increment of CV1 stays at 59.137 °C. The BL Tune increases until 58.875 °C and fine-tunes to 58.817 °C. However, both stabilize back to 58.81 °C. The MMPC final values of CV2, CV3, and CV4 are 43.023 °C, 49.916%, and 9.9%, meaning their differences are 2.27 °C, 0.0014%, and −2.13%.
When both disturbances are added, the response is similar to when the interference is only the temperature. The condenser temperature (CV1) decreases, and the level (CV4) is barely disturbed. The MeRecy temperature increases until 0.05 °C for the BLT but barely 0.01 °C for the fine-tuning. WG has a higher oscillation than both, which shows that feed temperature interference has a higher impact, and the best controller needs to be provided. Both fine-tuned and WG PI controllers reveal a faster stabilized graphic and a smaller disturbance peak.
Their IAEs and ISEs are given in Table 10 and Table 11. The lowest error parameter values against temperature and flow disturbance are attributed to WG for CV1 and CV4, FT for CV2, and BLT for CV3. When both disturbances are added, the lowest error parameter values are attributed to fine-tuning for CV2. For CV1 and CV4, the parameters of Wahid and Gunawan (2015) have the lowest error values. Again, BLT has the lowest error values for CV3.
Other than during the setpoint change against itself, the liquid percent level of Condenser T-103 (CV3) controller is the best displayed by the big log-modulus tuned PI controller. For the total liquid volume percent of T-103 (CV4), the best controller is shown by that of Wahid and Gunawan (2015), except against the change of CV3.
The response against setpoint changes is shown in Figure 6 and Table 12, Table 13 and Table 14. In Figure 6a, CV1 or the Vessel Condenser of the T-103 temperature is increased by 5%, from 58.81 °C to 61.7505 °C. All tuned controllers of CV1 display stability at the given final value after time 50. CV2 and CV4 are stable at 0, meaning no changes. CV3 has a decline, meaning that the value of liquid percent in the T-103 condenser decreases if CV1 is increased. If the condenser temperature changes, the column level is not affected because it is a previous process. Here, the liquid percent of the condenser decreases by about 12% by BLT tuning and almost 10% by fine-tuning. Neither affect the final temperature of the purified methanol stream. CV2 remains stable at 40.75 °C.
When CV2 is tempered, there should not be any changes towards CV1 if CV2 changes because there are no directly affected connections. Thus, in Figure 6b above, when E-101 cooler temperature changes, the rest of the controlled variables remains at approximately 0, both the condenser temperature and the column level, except for CV3.
In Figure 6c, when the CV3 setpoint is increased, CV1 and CV2 spike at the beginning of time when CV3 changes. The T-103 condenser temperature must be impacted by the change of its liquid level and thus the top product temperature. The fine-tuned PI controller has a faster response than BLT.
For CV4 setpoint tracking, the condenser temperature will be affected, and so will the output product of the purified methanol’s temperature. An increase in level means more product in the condenser, causing CV3 to increase as the CV4 setpoint is changed, indicating that G34 has a strong interaction. This is proven from the result of RGA analysis.
Wahid and Gunawan (2015) displayed the smallest and shortest deviation period for the CV1 temperature controller as well. CV2’s best improvement is shown by the fine-tuned Multi-loop PI controller. Table 15 below lists the final result parameters of this research.
To corroborate this performance, the final multi-loop PI controller parameter performance was tested in Unisim. The parameter was taken back to the simulation software and input into a PI controller tool that was later attached to the process. Although the error parameters for some are not similar or even close to similar, the PI still showed improvement compared to the MMPC performance, as shown in the Table 16 below. This shows that the accountability of a process model and a disturbance model enhances the response performance of the used controller.
The control performance of the multiloop PI is better than MMPC because MMPC is for non-linear processes modeled with Unisim, whereas, in this study, the controlled process is linear.

4. Conclusions

Based on the RGA, the correct installation of CV and MV are CV1-MV1, CV2-MV2, CV3-MV3, and CV4-MV4. To control the temperature is to control the heat flow, and to control the level is to control the flow. Proportional and Integral parameter values best to control CV1-MV1 pairing are 0.6667 and 0.3333, for CV2-MV2 are 4.9052 and 42.409, for CV3-MV3 are 0.0071 and 2.1603, also 13.5036 and 2.3776 for CV4-MV4.
Creating a disturbance model on the feed temperature and feed flow and focusing on each variable improved the controller performance in the methanol separation process in the DME Purification Plant by almost 100%; for CV1, CV2, CV3, and CV4, respectively, the IAE values are 79.49%, 99.90%, 100%, and 99.99, and the ISE values are 97.1%, 100%, 99.99%, and 100%. The linear model using FOPDT was proven reliable and realistic. Thus, in this study, multi-loop PI controllers have been proven effective against disturbance rejection, with big log modulus tuning, and the most suitable pairing. However, it is still unstable against setpoint changes. In addition to this drawback, this multi-loop PI controller has only been proven theoretically. A further implementation method on the physical hardware of a programmable logic controller (PLC) system needs to be explored.

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable. The study did not involve humans or animals.

Informed Consent Statement

Not applicable. The study did not involve humans or animals.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://drive.google.com/drive/folders/18Dwsix0dJ4HAKXDzeqTJMn0WZWGGY7ye?usp=sharing, accessed on 9 October 2021.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ministry of Energy and Mineral Resources Republic of Indonesia. Handbook of Energy & Economic Statistics of Indonesia; Head of Center for Data and Information Technology on Energy and Mineral Resources: Jakarta, Indonesia, 2018.
  2. Musalaiah, M.; Engineer, T.C.; Peddintaiah, K. A Theoretical Study on Dimethyl Ether: An Alternative Fuel for Future Generations. Int. Res. J. Eng. Technol. (IRJET) 2016, 3, 238–243. [Google Scholar]
  3. Takeishi, K. Dimethyl ether (DME): A clean fuel/energy for the 21st century and the low carbon society. Int. J. Energy Environ. 2016, 10, 238–242. [Google Scholar]
  4. Falco, M.D. Dimethyl Ether (DME) Production; University UCBM: Rome, Italy, 2017. [Google Scholar]
  5. Solichin, A.; Sari, M.; Rahmiyati, W.Y. Production of DME from Syngas (for Oxygenates in Diesel Oil and Blending LPG): Plant Design Report; Department of Chemical Engineering, Faculty of Engineering, Universitas Indonesia: Depok, Indonesia, 2011. [Google Scholar]
  6. Atef, G.; Almadani, A.; Noaman, S. Manufacturing of DME from Methanol; King Fahd University of Petroleum & Minerals: Dhahran, Saudi Arabia, 2015. [Google Scholar] [CrossRef]
  7. Dieterich, V.; Buttler, A.; Hanel, A.; Spliethoff, H.; Fendt, S. Power-to-liquid via synthesis of methanol, DME or Fischer–Tropsch-fuels: A review. R. Soc. Chem. 2020, 13, 3207–3252. [Google Scholar] [CrossRef]
  8. Khaleel, A. Methanol Dehydration to Dimethyl Ether over Highly Porous Xerogel Alumina Catalyst: Flow Rate Effect. Fuel Process Technol. 2010, 91, 1505–1509. [Google Scholar] [CrossRef]
  9. Dyrda, K.M.; Wilke, V.; Haas-Santo, K.; Dittmeyer, R. Experimental Investigation of the Gas/Liquid Phase Separation Using a Membrane-Based Micro Contactor. ChemEngineering 2018, 2, 55. [Google Scholar] [CrossRef] [Green Version]
  10. Liptak, B. Instrument Engineer’s Handbook: Process Control; CRC Taylor & Francis: Boca Raton, FL, USA, 2012. [Google Scholar] [CrossRef]
  11. Marlin, T. Process Control: Designing Processes and Control System for Dynamic Performance, 2nd ed.; McGraw Hill: Singapore, 2015; ISBN 9780070393622. [Google Scholar]
  12. Wahid, A.; Gunawan, T.A. Pengendalian Proses Purifikasi DME dan Metanol pada Pabrik DME dari Gas Sintesis. Sinergi 2015, 19, 57–66. [Google Scholar] [CrossRef] [Green Version]
  13. Wahid, A.; Brillianto, Z.H. Multivariable Model Predictive Control (4 × 4) of Methanol Water Separation in Dimethyl Ether Production. AIP Conf. Proc. 2015, 2255, 030055. [Google Scholar] [CrossRef]
  14. Nandong, J. Decentralized Control Systems; Curtin University: Curtin, Malaysia, 2020. [Google Scholar] [CrossRef]
  15. Vu Truong, N.; Lee, M. Optimal Design of Multi-loop PI Controllers for Enhanced Disturbance Rejection in Multivariable Processes; School of Chemical Engineering and Technology Yeungnam University: Gyeoungbuk, Korea, 2007; Available online: http://psdc.yu.ac.kr (accessed on 22 July 2020).
  16. Biyanto, T.; Sordi, N.; Sehamat, N.; Zabiri, H. PID Multivariable Tuning System Using BLT Method for Distillation Column; IOP Publishing: Tokyo, Japan, 2018. [Google Scholar] [CrossRef]
  17. Farsi, M. Dynamic Modeling and Controllability Analysis of DME Production in an Isothermal Fixed Bed Reactor; Department of Chemical Engineering, School of Chemical and Petroleum Engineering, Shiraz University: Shiraz, Iran, 2015. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, H.; Jeng, J.; Chiang, C.; Pan, W. A Direct Method for Multi-Loop PI/PID Controller Design; Department of Chemical Engineering, National Taiwan University: Taipei, China, 2003. [Google Scholar] [CrossRef]
  19. Kumar, R.; Anand, S.; Khulbey, A.; Jha, A.N. Design of Fractional Order Controller for Wood-Berry Distillation Column; Department of Instrumentation and Control Engineering, Netaji Subhas University of Technology: New Delhi, India, 2020; ISBN 9781728169163. [Google Scholar]
  20. Kalpana, R.; Harukumar, K.; Senthukumar, J.; Balasubramanian, G.; Gour, A. Multivariable Static Output Feedback Control of a Binary Distillation Column Using Linear Matrix Inequalities and Genetic Algorithm; School of Electrical and Electronics Engineering, Sastra University: Thanjavur, India, 2017. [Google Scholar] [CrossRef]
  21. Kalhoro, A.; Memon, S. Design of Multivariable PID Controllers: A Comparative Study; Department of Electronic Engineering, Mehran University: Jamshoro, Pakistan, 2021. [Google Scholar] [CrossRef]
  22. Alawad, N.; Muter, A. Optimal Control of Model Reduction Binary Distillation Column; Department of Computer Engineering, Faculty of Engineering, Mustansiriyah University: Baghdad, Iraq, 2021. [Google Scholar] [CrossRef]
  23. Ogunnaike, B.A.; Ray, W.H. Multivarianle Design Controller Design for Linear Systems Having Multiple Time Delays; Department of Chemical Engineering, University of Winconsin: Madison, WI, USA, 1979. [Google Scholar]
  24. Smith, C. Digital Computer Process Control; Intext: Scranton, PA, USA, 1972; ISBN 0700224017/9780700224012. [Google Scholar]
  25. Camacho, E.F.; Bordons, C. Model Predictive Control, 2nd ed.; Springer: London, UK, 2007; ISBN 9780857293985. [Google Scholar]
  26. Vazquez, F.; Morilla, F. Tuning Decentralized PID Controllers for MIMO Systems with Decouplers. In Proceedings of the 15th Triennial World Congress, Barcelona, Spain, 21–26 July 2002. [Google Scholar] [CrossRef]
  27. Krishna, G.H.; Kiranmayi, R.; Rathaiah, M. Control System Design for 3 × 3 Processes Based on Effective Transfer Function and Fractional Order Filter. Int. J. Eng. Adv. Technol. (IJEAT) 2019, 8, 637–642. [Google Scholar]
  28. Chau, P.C. Chemical Process Control: A First Course with MATLAB; University of California: San Diego, CA, USA, 2001. [Google Scholar] [CrossRef]
  29. Wahid, A.; Tanuwijaya, R. Pemilihan Metode Penyetelan Pengendali PI pada Pengendalian Pabrik Regasifikasi LNG Menggunakan Metode Skor. In Proceedings of the Seminar Nasional Teknik Kimia UNPAR, Bandung, Indonesia, 19 November 2015. [Google Scholar]
Figure 1. Methanol purification process simulation.
Figure 1. Methanol purification process simulation.
Chemengineering 05 00070 g001
Figure 2. Process reaction curve (PRC) of (a) Method I, (b) Method II, and (c) Method III (Source: Smith, 1972; Wade, 2004).
Figure 2. Process reaction curve (PRC) of (a) Method I, (b) Method II, and (c) Method III (Source: Smith, 1972; Wade, 2004).
Chemengineering 05 00070 g002
Figure 3. Block diagram of the 2 × 2 MIMO process system multi-loop controller. (Source: Nandong, 2020) [14].
Figure 3. Block diagram of the 2 × 2 MIMO process system multi-loop controller. (Source: Nandong, 2020) [14].
Chemengineering 05 00070 g003
Figure 4. PRC of methanol purification controlled variables after the addition of (a) feed-in flow disturbance and (b) feed-in temperature disturbance.
Figure 4. PRC of methanol purification controlled variables after the addition of (a) feed-in flow disturbance and (b) feed-in temperature disturbance.
Chemengineering 05 00070 g004
Figure 5. Responses of the methanol purification process using BLT, FT, or WG (top to bottom) with disturbance(s) of (a) feed temperature, (b) feed flow, and (c) feed temperature and feed flow (from left to right). Vessel temperature of Column T-103 (CV1) is colored blue, temperature of the MeRecy stream (CV2) is colored yellow, the liquid percent level of Condenser T-103 (CV3) is colored green, and the total volume percent of T-103 (CV4) is colored red. Step of CV1 is colored orange, Step of CV2 is colored purple, Step of CV3 is colored light blue, and Step of CV4 is colored blue. Note that the CVs are analyzed, so none are manually increasing, and thus there are no changes in the step lines.
Figure 5. Responses of the methanol purification process using BLT, FT, or WG (top to bottom) with disturbance(s) of (a) feed temperature, (b) feed flow, and (c) feed temperature and feed flow (from left to right). Vessel temperature of Column T-103 (CV1) is colored blue, temperature of the MeRecy stream (CV2) is colored yellow, the liquid percent level of Condenser T-103 (CV3) is colored green, and the total volume percent of T-103 (CV4) is colored red. Step of CV1 is colored orange, Step of CV2 is colored purple, Step of CV3 is colored light blue, and Step of CV4 is colored blue. Note that the CVs are analyzed, so none are manually increasing, and thus there are no changes in the step lines.
Chemengineering 05 00070 g005
Figure 6. Responses against setpoint tracking using BLT, FT, and WG (top to bottom) in terms of (a) the vessel temperature of Column T-103 (CV1), (b) the temperature of the MeRecy stream (CV2), (c) the liquid percent level of Condenser T-103 (CV3), and (d) the total volume percent of T-103 (CV4). CV1 is colored blue, CV2 is colored yellow, CV3 is colored green, and CV4 is colored red. Step of CV1 is colored orange, Step of CV2 is colored purple, Step of CV3 is colored light blue, and Step of CV4 is colored blue.
Figure 6. Responses against setpoint tracking using BLT, FT, and WG (top to bottom) in terms of (a) the vessel temperature of Column T-103 (CV1), (b) the temperature of the MeRecy stream (CV2), (c) the liquid percent level of Condenser T-103 (CV3), and (d) the total volume percent of T-103 (CV4). CV1 is colored blue, CV2 is colored yellow, CV3 is colored green, and CV4 is colored red. Step of CV1 is colored orange, Step of CV2 is colored purple, Step of CV3 is colored light blue, and Step of CV4 is colored blue.
Chemengineering 05 00070 g006
Table 1. Cooler and pump process parameters.
Table 1. Cooler and pump process parameters.
UnitCodeObjectivesInitialFinal
CoolerE-101Cools methanol distillates240 ℉110 ℉
PumpP-101To increase the recycle flow pressure200 kPa3500 kPa
(Source: Solichin et al. 2011) [5].
Table 2. Distillation column process parameters.
Table 2. Distillation column process parameters.
Process ParametersT-103
Number of Stages30
Condenser Pressure14.69 psia
Reboiler Pressure14.69 psia
Stream Inlet Stage19
(Source: Solichin et al. 2011) [5].
Table 3. CV and MV of the methanol purification process.
Table 3. CV and MV of the methanol purification process.
NoUnit ProcessController TypeControlled VariableManipulated VariableSet PointActuatorDisturbance
1Condenser T-103Temperature controlTemperatureHeat flow58.81 °CControl Valve Qcond1Changes in the inlet feed (methanol + water) temperature and molar flow
2Cooler E-101Heat flow40.75 °CControl Valve QcoolMeRecy
3Condenser T-103Level
control
Column LevelMethanol to recycle Flow49.96%Actuator Desired Position VLV- 105
4Column T-103Wastewater flow11.22%Actuator Desired Position VLV- 106
Table 4. Relative gain array effect.
Table 4. Relative gain array effect.
λPossible Pairing
λmn < 0Unstable interactions
λmn = 1No interaction exists
λmn = 0.5Strong interactions
λmn > 1Interactions in the opposite direction of the variables, and stronger as value increases
λmn < 1Interactions in the same direction of the variables, and stronger as the value decreases
Table 5. Quick summaries of several possibilities.
Table 5. Quick summaries of several possibilities.
ΛPossible Pairing
Λ = 1 0 0 1 No interaction. The controller design is a SISO system.
Λ = 1 δ δ 1 Strong interaction if 𝛿 is close to 1; weak interaction if 𝛿 << 1.
(Source: Chau, 2001) [28].
Table 6. Bode plot stability parameters.
Table 6. Bode plot stability parameters.
StableMarginally StableUnstable
Phase over frequency > gain cross over frequencyPhase over frequency = gain cross over frequencyPhase over frequency < gain cross over frequency
Gain margin > 1Gain margin = 1Gain margin < 1
Phase margin = positivePhase margin = 0Phase margin = negative
(Source: Marlin, 2015) [11].
Table 7. FOPDT of the methanol purification disturbance model.
Table 7. FOPDT of the methanol purification disturbance model.
DisturbanceVariable Process
CV1CV2CV3CV4
Inlet Flow   0.0001 e 1.15 s 1.00 s + 1 0.0065 e 8.32 s 3.63 s + 1 0.00003 e 1.72 s 1.38 s + 1 0.0016 e 7.95 s 3.00 s + 1
Inlet Temperature   0.1066 e 2.48 s 1.63 s + 1 0.0902 e 2.48 s 1.50 s + 1 0.0148 e 1.94 s 1.32 s + 1 0.0275 e 9.45 s 3.00 s + 1
Table 8. Proportional integral (PI) controller parameters of the methanol purification process.
Table 8. Proportional integral (PI) controller parameters of the methanol purification process.
Controlled Variable (CV)Controller ParameterBig Log-Modulus TuningFine
Tuning
Wahid and Gunawan (2015)
1Vessel Temperature of Condenser T-103P0.670.241
I0.331.640.5
2MeRecy TemperatureP0.674.901
I0.3342.410.5
3Liquid Percent Level of Condenser T-103P0.01−67.010.01
I2.166.75 × 10−191.80
4Total Liquid Volume Percent of T-103P9.003.4913.50
I1.593.562.38
Table 9. Final value of Simulink’s step block.
Table 9. Final value of Simulink’s step block.
CVSet PointIncrease 5%Delta (Δ)
158.81 ℃61.75 ℃2.94
240.75 ℃42.79 ℃2.04
349.96%52.46%2.5
411.22%11.78%0.56
Disturbance
Feed Temperature81.74 °C85.83 °C4.09
Feed Flow659.8 gmol/s692.72 gmol/s32.98
Table 10. IAE comparison of MMPC and multi-loop PI against disturbance rejection.
Table 10. IAE comparison of MMPC and multi-loop PI against disturbance rejection.
CVIAEImprovement
MMPC
(Wahid and Brillianto 2020)
Multi-Loop PIMulti-Loop PI
Big Log-Modulus TuningWahid & Gunawan (2015)Fine-TuningBig Log-Modulus TuningWahid & Gunawan (2015)Fine-Tuning
Feed Temperature Disturbance
13.70 × 1011.05 × 1014.716.6571.58%87.27%82.01%
27.96 × 1011.242.98 × 101.62 × 10−198.45%96.26%99.80%
36.186.26 × 10−96.84 × 1027.93 × 10−5100.00%−1.10 × 104%100.00%
42.04 × 1022.47 × 10−21.10 × 10−22.82 × 10−299.99%99.99%99.97%
Feed Flow Disturbance
13.75 × 1013.11 × 10−11.53 × 10−12.09 × 10−199.17%99.59%99.44%
21.55 × 1024.97 × 10−19.20 × 10−11.18 × 10−299.68%99.41%99.99%
33.321.46 × 10−61.62 × 1011.87 × 10−6100.00%−385.80%100.00%
47.41 × 1011.16 × 10−25.15 × 10−31.33 × 10−299.98%99.99%99.98%
Feed Flow and Temperature Disturbance
12.24 × 1011.02 × 1014.596.4854.19%79.49%71.02%
21.64 × 1021.472.661.70 × 10−199.11%98.38%99.90%
31.24 × 1016.12 × 10−56.69 × 1027.74 × 10−5100.00%−5.31 × 103%100.00%
42.04 × 1022.47 × 10−21.10 × 10−22.82 × 10−299.98%99.99%99.98%
Table 11. ISE comparison of MMPC and multi-loop PI against disturbance rejection.
Table 11. ISE comparison of MMPC and multi-loop PI against disturbance rejection.
CVISEImprovement
MMPC
(Wahid and Brillianto 2020)
Multi-Loop PIMulti-Loop PI
Big Log-Modulus TuningWahid & Gunawan (2015)Fine-TuningBig Log-Modulus TuningWahid & Gunawan (2015)Fine-Tuning
Feed Temperature Disturbance
12.30 × 1012.038.01 × 10−11.4191.20%96.52%93.89%
21.06 × 1029.90 × 10−21.58 × 10−17.03 × 10−499.91%99.85%100.00%
36.41 × 10−18.31 × 10−59.53 × 1021.29 × 10−499.99%−1.49 × 105%99.98%
46.94 × 1021.01 × 10−42.39 × 10−51.63 × 10−4100.00%100.00%100.00%
Feed Flow Disturbance
16.62 × 1011.33 × 10−35.47 × 10−41.02 × 10−3100.00%100.00%100.00%
24.43 × 1022.17 × 10−24.20 × 10−11.62 × 10−5100.00%99.91%100.00%
37.68 × 10−14.66 × 10−85.41 × 10−17.31 × 10−8100.00%29.57%100.00%
49.64 × 1012.23 × 10−55.35 × 10−63.70 × 10−5100.00%100.00%100.00%
Feed Flow and Temperature Disturbance
12.68 × 1011.957.79 × 10−11.3692.71%97.09%94.91%
23.85 × 1021.08 × 10−11.29 × 10−18.38 × 10−499.97%99.97%100.00%
31.837.93 × 10−59.10 × 1021.22 × 10−4100.00%−4.96 × 104%99.99%
46.94 × 1021.01 × 10−42.39 × 10−51.63 × 10−4100.00%100.00%100.00%
Table 12. IAE comparison of MMPC and multi-loop PI.
Table 12. IAE comparison of MMPC and multi-loop PI.
CVIAEImprovement
MMPC
(Wahid and Brillianto 2020)
Multi-Loop PIMulti-Loop PI
Big Log-Modulus TuningWahid & Gunawan (2015)Fine-TuningBig Log-Modulus TuningWahid & Gunawan (2015)Fine-Tuning
CV1 Setpoint Change
11.81 × 1026.61 × 1012.96 × 1014.17 × 10163.51%83.66%77.02%
24.288.979.261.04−109.63%−116.24%75.70%
34.364.08 × 1034.16 × 1034.82 × 103−93,386.24%−9.53 × 104%−1.10 × 105%
42.00 × 10−22.29 × 10−31.37 × 10−31.65 × 10−288.54%93.16%17.70%
CV2 Setpoint Change
14.10 × 10−17.34 × 10−13.28 × 10−14.60 × 10−1−79.02%19.95%−12.15%
25.64 × 1011.83 × 1018.202.71 × 10−167.48%85.46%99.52%
33.20 × 10−13.59 × 1013.73 × 1014.35 × 101−11,112.50%−11,546.88%−13,493.75%
42.00 × 10−21.68 × 10−48.12 × 10−56.45 × 10−499.16%99.59%96.78%
Table 13. IAE comparison of MMPC and multi-loop PI (Cont’d).
Table 13. IAE comparison of MMPC and multi-loop PI (Cont’d).
CVIAEImprovement
MMPC
(Wahid and Brillianto 2020)
Multi-Loop PIMulti-Loop PI
Big Log-Modulus TuningWahid & Gunawan (2015)Fine-TuningBig Log-Modulus TuningWahid & Gunawan (2015)Fine-Tuning
CV3 Setpoint Change
17.327.515.806.97−2.60%20.71%4.77%
24.679.21 × 10−12.203.80 × 10−180.29%52.89%91.86%
35.59 × 1014.16 × 1023.92 × 1021.34 × 103−643.43%−601.59%−2301.22%
45.00 × 10−21.43 × 10−14.61 × 10−11.27−185.20%−821.40%−2430.00%
CV4 Setpoint Change
15.20 × 10−17.191.311.84−1282.31%151.15%−254.42%
25.30 × 10−12.98 × 10−17.42 × 10−14.31 × 10−243.70%−39.91%91.87%
35.791.73 × 1021.89 × 1022.19 × 102−2886.18%−3.17 × 103%−3.68 × 103%
45.74 × 1011.42 × 10−11.03 × 10−12.84 × 10−199.75%99.82%99.51%
Table 14. ISE comparison of MMPC and multi-loop PI.
Table 14. ISE comparison of MMPC and multi-loop PI.
CVISEImprovement
MMPC
(Wahid and Brillianto 2020)
Multi-Loop PIMulti-Loop PI
Big Log-Modulus TuningWahid & Gunawan (2015)Fine-TuningBig Log-Modulus TuningWahid & Gunawan (2015)Fine-Tuning
CV1 Setpoint Change
14.35 × 1028.90 × 1013.92 × 1016.42 × 10179.54%90.99%85.23%
24.00 × 10−21.72 × 10−11.77 × 10−13.37 × 10−2−328.75%−343.50%15.80%
36.00 × 10−23.39 × 1043.49 × 1044.71 × 104−5.65 × 107%−5.82 × 107%−7.84 × 107%
45.36 × 10−79.18 × 10−84.40 × 10−89.48 × 10−682.88%91.79%−1669.03%
CV2 Setpoint Change
15.00 × 10−47.45 × 10−33.07 × 10−36.63 × 10−3−1389.60%514.60%−1226.20%
27.11 × 1011.44 × 1015.942.49 × 10−179.68%91.6425%99.65%
32.00 × 10−42.732.883.911.37 × 106%−1.44 × 106%−1.95 × 106%
48.30 × 10−74.43 × 10−93.29 × 10−92.17 × 10−799.47%99.60%73.86%
CV3 Setpoint Change
11.00 × 10−24.86 × 10−11.851.90−4759%−18,370.00%−18,930.00%
26.00 × 10−27.18 × 10−31.70 × 10−17.38 × 10−288.04%−182.83%−23.00%
33.10 × 1013.71 × 1023.19 × 1021.22 × 103−1099.23%−930.68%−3845.75%
45.80 × 10−64.27 × 10−25.69 × 10−11.59−7.36 × 105%−9.82 × 106%−2.73 × 107%
CV4 Setpoint Change
14.00 × 10−44.79 × 10−19.19 × 10−21.35 × 10−1−1.20 × 105%−2.29 × 104%−3.37 × 104%
25.00 × 10−46.56 × 10−39.32 × 10−36.75 × 10−5−1.21 × 103%−1764.80%86.51%
31.30 × 10−16.29 × 1017.23 × 1019.73 × 101−48,30%−5.55 × 104%−7.47 × 104%
42.50 × 1014.20 × 10−22.80 × 10−27.99 × 10−299.83%99.89%99.68%
Table 15. Final controller parameter values.
Table 15. Final controller parameter values.
Controlled Variable (CV)Tuning MethodPI
1Vessel Temperature of Condenser T-103WG0.670.33
2MeRecy TemperatureFT4.9042.41
3Liquid Percent Level of Condenser T-103BLT0.012.16
4Total Liquid Volume Percent of T-103WG13.502.38
Table 16. IAE comparison of multi-loop PI in Simulink and Unisim.
Table 16. IAE comparison of multi-loop PI in Simulink and Unisim.
CVIAEImprovementISEImprovement
MMPC (Wahid & Brillianto 2020)Multi-LoopMMPC (Wahid & Brillianto 2020)Multi-Loop
Tuning
Method
SimulinkUnisimSimulinkUnisimSimulinkUnisimSimulinkUnisim
CV1 Setpoint Change
11.81 × 102WG2.96 × 1012.08 × 10−383.66%100.00%4.35 × 1023.92 × 1012.65 × 10−490.99%100.00%
24.28FT1.041.15 × 10−275.70%99.73%4.00 × 10−23.37 × 10−21.92 × 10−415.80%99.52%
34.36BLT4.08 × 1031.34 × 10−2−93,386.24%99.69%6.00 × 10−23.39 × 1041.83 × 10−4−5.65 × 107%99.69%
42.00 × 10−2WG1.37 × 10−31.37 × 10−393.16%93.16%5.36 × 10−74.40 × 10−82.89 × 10−791.79%46.18%
CV2 Setpoint Change
14.10 × 10−1WG3.28 × 10−11.27 × 10−219.95%96.91%5.00 × 10−43.07 × 10−32.45 × 10−4−514.60%51.05%
25.64 × 101FT2.71 × 10−12.76 × 10−199.52%99.51%7.11 × 1012.49 × 10−16.43 × 10−199.65%99.10%
33.20 × 10−1BLT3.59 × 1013.25 × 10−3−11,112.50%98.99%2.00 × 10−42.731.48 × 10−5−1.37 × 106%92.59%
42.00 × 10−2WG8.12 × 10−51.36 × 10−399.59%93.22%8.30 × 10−73.29 × 10−92.78 × 10−799.60%66.51%
CV3 Setpoint Change
17.32WG5.801.25 × 10−220.71%99.83%1.00 × 10−21.852.40 × 10−4−18370.00%97.60%
24.67FT3.80 × 10−11.11 × 10−291.86%99.76%6.00 × 10−27.38 × 10−21.83 × 10−4−23.00%99.70%
35.59 × 101BLT4.16 × 1022.66 × 10−3−643.43%100.00%3.10 × 1013.71 × 1021.47 × 10−5−1099.23%100.00%
45.00 × 10−2WG4.61 × 10−18.26 × 10−4−821.40%98.35%5.80 × 10−65.69 × 10−11.10 × 10−6−9.82 × 106%81.06%
CV4 Setpoint Change
15.20 × 10−1WG1.311.28 × 10−2−151.15%97.53%4.00 × 10−49.19 × 10−22.51 × 10−4−2.29 × 104%37.25%
25.30 × 10−1FT4.31 × 10−21.14 × 10−291.87%97.84%5.00 × 10−46.75 × 10−51.91 × 10−486.51%61.80%
35.79BLT1.73 × 1023.26 × 10−3−2886.18%99.94%1.30 × 10−16.29 × 1011.50 × 10−5−4.83 × 104%99.99%
45.74 × 101WG1.03 × 10−14.04 × 10−499.82%100.00%2.50 × 1012.80 × 10−22.83 × 10−799.89%100.00%
Feed Temperature Disturbance
13.70 × 101WG4.711.61 × 10−187.27%99.56%2.30 × 1018.01 × 10−12.78 × 10−296.52%99.88%
27.96 × 101FT1.62 × 10−11.07 × 10−299.80%99.99%1.06 × 1027.03 × 10−41.14 × 10−4100.00%100.00%
36.18BLT6.26 × 10−99.44 × 10−3100.00%99.85%6.41 × 10−18.31 × 10−59.06 × 10−599.99%99.99%
42.04 × 102WG1.10 × 10−22.52 × 10−399.99%100.00%6.94 × 1022.39 × 10−51.17 × 10−5100.00%100.00%
Feed Flow Disturbance
13.75 × 101WG1.53 × 10−11.67 × 10−199.59%99.56%6.62 × 1015.47 × 10−44.22 × 10−5100.00%100.00%
21.55 × 102FT1.18 × 10−21.38 × 10−299.99%99.99%4.43 × 1021.62 × 10−51.15 × 10−4100.00%100.00%
33.32BLT1.46 × 10−67.36 × 10−4100.00%99.98%7.68 × 10−14.66 × 10−86.73 × 10−7100.00%100.00%
47.41 × 101WG5.15 × 10−39.42 × 10−399.99%99.99%9.64 × 1015.35 × 10−69.45 × 10−5100.00%100.00%
Feed Flow and Temperature Disturbance
12.24 × 101WG4.591.70 × 10−179.49%99.24%2.68 × 1017.79 × 10−13.01 × 10−297.09%99.89%
21.64 × 102FT1.70 × 10−11.14 × 10−299.90%99.99%3.85 × 1028.38 × 10−41.14 × 10−4100.00%100.00%
31.24 × 101BLT6.12 × 10−53.42 × 10−3100.00%99.97%1.837.93 × 10−52.96 × 10−5100.00%100.00%
42.04 × 102WG1.10 × 10−23.27 × 10−299.99%99.98%6.94 × 1022.39 × 10−51.43 × 10−3100.00%100.00%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wahid, A.; Rikimata, M. Methanol–Water Purification Control Using Multi-Loop PI Controllers Based on Linear Set Point and Disturbance Models. ChemEngineering 2021, 5, 70. https://doi.org/10.3390/chemengineering5040070

AMA Style

Wahid A, Rikimata M. Methanol–Water Purification Control Using Multi-Loop PI Controllers Based on Linear Set Point and Disturbance Models. ChemEngineering. 2021; 5(4):70. https://doi.org/10.3390/chemengineering5040070

Chicago/Turabian Style

Wahid, Abdul, and Monica Rikimata. 2021. "Methanol–Water Purification Control Using Multi-Loop PI Controllers Based on Linear Set Point and Disturbance Models" ChemEngineering 5, no. 4: 70. https://doi.org/10.3390/chemengineering5040070

Article Metrics

Back to TopTop