Next Article in Journal
Building Sector Issues in about 100 Years: End-Of-Life Scenarios of Carbon-Reinforced Concrete Presented in the Context of a Life Cycle Assessment, Focusing the Carbon Footprint
Next Article in Special Issue
Special Issue on “Bioethanol Production Processes”
Previous Article in Journal
Testing and Analysis of Torsional Vibration of Ship Transmission Shafting Based on Five-Point Smoothing Algorithm
Previous Article in Special Issue
Introducing a Marine Biorefinery System for the Integrated Production of Biofuels, High-Value-Chemicals, and Co-Products: A Path Forward to a Sustainable Future
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol

by
Carlos Alberto Torres Cantero
1,2,
Ricardo Pérez Zúñiga
3,
Mario Martínez García
4,*,
Silvia Ramos Cabral
4,
Manuela Calixto-Rodriguez
5,
Jorge Salvador Valdez Martínez
5,
Mayra Guadalupe Mena Enriquez
6,
Abraham Jashiel Pérez Estrada
7,
Gerardo Ortiz Torres
4,
Felipe de J. Sorcia Vázquez
4,
Azael García Rebolledo
2 and
Jesse Yoe Rumbo Morales
4
1
Tecnológico Nacional de Mexico Campus Colima, Av. Tecnológico # 1, Col. Liberación, Villa de Álvarez 28976, Colima, Mexico
2
Facultad de Ingeniería Mecánica y Eléctrica de la Universidad de Colima, Carretera Colima-Coquimatlan Km 9, Valle de las Huertas, Coquimatlán 28400, Colima, Mexico
3
Departamento en Ciencias Computacionales e Ingenierías, Sistema de Universidad Virtual de la Universidad de Guadalajara, Av. Juárez No. 976, Colonia Centro, Guadalajara 44100, Jalisco, Mexico
4
Centro Universitario de los Valles, Universidad de Guadalajara, Carretera Guadalajara-Ameca Km 45.5, Ameca 46600, Jalisco, Mexico
5
División Académica de Mecánica Industrial, Universidad Tecnológica Emiliano Zapata del Estado de Morelos, Avenida Universidad Tecnológica No. 1, Col. Palo Escrito, Emiliano Zapata 62760, Morelos, Mexico
6
Centro de Investigación Multidisciplinario en Salud, Departamento de Ciencias Biomédicas, Centro Universitario de Tonalá, Tonala 45425, Jalisco, Mexico
7
School of Engineering and Sciences, Tecnologico de Monterrey, Monterrey 64849, Nuevo Leon, Mexico
*
Author to whom correspondence should be addressed.
Processes 2022, 10(9), 1792; https://doi.org/10.3390/pr10091792
Submission received: 2 August 2022 / Revised: 24 August 2022 / Accepted: 24 August 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Bioethanol Production Processes)

Abstract

:
Extractive distillation with salts, unlike other dehydration technologies, is better due to the null toxicity that exists in the distillate, since salt cannot be evaporated. With this distillation technology, it is possible to obtain a high concentration of ethanol, however, there are still problems in the control of the distillation columns in the presence of disturbances. The present work deals with the simulation and control of an extractive distillation column using CaCl 2 as a separating agent, for which the Aspen Dynamics® simulator is used. The measurement and control of the ethanol composition are carried out by means of temperature, in addition, four control structures are evaluated and compared. These structures are L, D, LV, and DV, which are the most common in conventional distillation, and their performance is measured by means of deterministic indicators applying changes (disturbances) of composition and the flow rate in the main feed of the column. The most relevant results of this work lead to the fact that by applying a controller, it is possible to maintain the desired purity above the international purity standards (99% ethanol) that govern biofuels.

1. Introduction

In recent years, research has focused on the search for alternative sources of renewable energy. This is to reduce the dependence on fossil fuels and the greenhouse gas emissions derived from them. To reduce this dependence on fossil fuels, bioethanol has been introduced, produced through fermentation technologies from sugars and starches from raw materials, such as cane, corn, wheat, beet, barley, and sorghum, among others. This type of bioethanol production from these raw materials is known as first-generation biofuel, there are also second- and third-generation biofuels, which use materials, such as plants, trees, animal waste, algae, marine plants, or any other type of vegetable or animal waste. In addition, another type of raw material that is being used recently is the food waste that we generate throughout the food supply chain, from production to human consumption. This type of raw material is known as biomass and is divided into vegetable and lignocellulosic [1,2,3,4,5,6,7].
To produce this type of bioethanol, it is necessary to remove excess water from the ethanol-water mixture produced by fermentation. Depending on the technique used to ferment, the raw material, the bacteria used, and the environmental conditions, it can produce between 3% and 15% molar fraction [8,9]. This concentration is too low to be used as a fuel or additive, because when burned it will not have the same efficiency; that is, it will contain most of the water vapor and this will prevent the ignition and burning of the fuel, not to mention the damage that it can cause inside the engine due to the condensation of the water vapor. Therefore, to be used as a fuel or additive, bioethanol must have a purity equal to or greater than 99.5 wt%, in accordance with international standards (EN 15376; ASTMD 5798-11) [10,11].
To achieve the concentration required by the aforementioned standards, it is necessary to dehydrate the ethanol using non-ideal mixture separation technologies that allow changing or correcting the non-ideal VLE (vapor-liquid equilibrium) of the ethanol-water system to an ideal VLE. As has been mentioned and described in different articles, the ethanol-water mixture is a non-ideal mixture that at a composition of 89% molar fraction of ethanol, with an approximate temperature of 78.17 ° C and at a pressure of 1.01325 bar (1 atm), creates an azeotrope [12,13,14,15,16,17].
To break, eliminate, or move this azeotrope it is necessary to use one of the different types of distillation that can be used to dehydrate ethanol, such as: azeotropic distillation with: salts (using calcium chloride (CaCl 2 ), potassium acetate, sodium acetate, etc.), solvents, (benzene, gasoline, glycerol, etc.), vacuum distillation (which could theoretically), ionic liquids. Other separation technologies to obtain bioethanol are adsorption (using some types of zeolite), membrane processes, and reverse osmosis, or finally, the combination of processes, which are known as hybrid processes [18,19,20,21,22,23,24,25,26,27,28,29,30,31].
The present work is based on using conventional distillation column control structures with classical control techniques, PI and PID, to control the production of fuel-grade ethanol in an extractive distillation column with salts, using CaCl 2 to 16.7 wt% as a separating agent.
The purpose of these controllers is to ensure two essential aspects: product purity in accordance with international fuel ethanol quality standards in terms of water content, as well as a correct operation before input changes and oscillatory disturbances that represent more realistic behaviors in the said distillation column. Another important contribution of this work is the addition of a degree of freedom to the extractive distillation model to feed the salt through the control signals. In addition, time delays of the measured variables used in the control system are considered, which alter the dynamics of the extractive distillation column.
In addition, in this work, the extractive distillation column is simulated and controlled using the aspenONE® suite. The simulation of the technique used includes the prediction of the VLE with CaCl 2 considering experimental data provided by Nishi [32] and activity coefficients for the E-NRTL thermodynamic model for the ethanol-water-CaCl 2 system are estimated using the Aspen Properties Data Regression System® (DRS) and are reported in Cantero et al. [1]. The extractive distillation simulation uses the Aspen Plus® RadFrac block, which performs rigorous calculations for columns in ordinary distillation, absorption, desorption, extractive, and azeotropic distillation, three-stage distillation, and reactive distillation applications; it solves MESH equations, including mass balance equations (M), equilibrium relations (E), composition summation (S), and enthalpy balance (H). The distillation simulation is carried out in a steady state, where the molar flows, compositions, and temperatures are analyzed and compared with those reported in Llano-Restrepo et al. [16] and Hashemi et al. [33].
On the other hand, this process, due to the complexity of the non-linearity of the process and the non-ideal mixture, allows future work to be carried out, such as the use of different types of control (intelligent, predictive, multivariable, etc.), the use of different techniques of control structures (DB, R/F, LB, etc.), the optimization of the plant (separating agent, stages, heats, etc.), the diagnosis and isolation of faults, and most importantly, the economic evaluation, and energy consumption [1,2,3,18,34].
The document is organized in section sections, of which, Section 2 deals with the mathematical model used for VLE prediction and the distillation column. Section 3 describes the behavior of the process in a stable state. Section 4 presents the results of the analyses used to design the salt distillation column. Section 5 presents the results of the control structures and their behavior disturbances. Section 6 analyzes the error rates, according to integral error measurements and classical performance indicators, and interprets the most effective control structure. Finally, Section 7 offers a synthesis of the objective and importance of the work.

2. Mathematical Model of the Extractive Distillation Column with Salt

The simulation of the extractive distillation technique with CaCl 2 was carried out using the aspenONE® suite, more specifically, the following suite programs. Aspen Properties® DRS v8.6 was used in the simulations of this work to obtain the VLE database. To obtain the molar flows, compositions, temperatures, heat, or enthalpy balances in a stable state, the Aspen Plus® v8.6 simulator was used. Finally, for the simulation of the dynamic state behavior of the process, together with the classic control techniques, Aspen Plus Dynamics® v8.4 was used.
The thermodynamic model, Equation (1), used in the simulation for the VLE prediction was E-NRTL [35,36]; this thermodynamic model is contained in the aspenONE® suite. m and j take the values of components 1 (water) and 2 (ethanol) in the same way as the activity coefficients for the cations, Equation (2), and anions, Equation (3). This model was chosen because it is the most recommended for this type of system [37].
G E , l c R T = m x m j x j G j m τ j m k x k G k m + c x c a x a a x a × j x j G j c , a c τ j c , a c k x k G k c , a c + a x a c x c c x c × j x j G j a , c a τ j a , c a k x k G k a , c a
1 Z c l n γ c l c = a x a a x a × k x k G k c , a c τ k c , a c k x k G k c , a c + m x m G c m k x k G k m × τ c m k x k G k m τ k m k x k G k m + a c x c c x c × x a G c a , c a k x k G k a , c a × τ c a , c a k x k G k a , c a τ k a , c a k x k G k a , c a
1 Z a l n γ a l c = c x c c x c × k x k G k a , c a τ k c , a c k x k G k a , c a + m x m G a m k x k G k m × τ a m k x k G k m τ k m k x k G k m + c a x a a x a × x a G a c , a c k x k G k c , a c × τ a c , a c k x k G k c , a c τ k c , a c k x k G k c , a c
To demonstrate that the selected model adequately represents the systems, the results obtained for molar flows, composition, and temperatures in this work were compared with the corresponding ones offered by the consulted literature. Subsequently, it was verified that the relative error obtained was less than 5%, as it is the most recommended value of this type of study.
In order to make the comparison between the simulated dehydration technologies, the molar flow of ethanol in the ethanol-water mixture was taken into account to be the same for the technologies analyzed (100 kmol/h).

Description of the Extractive Distillation Column with Salt

Extractive distillation using salt as a separating agent operates on the same principle as conventional distillation. The main difference that exists is an extractive distillation column has a second feed typically located in the second tray of the column, and this is where the salt is introduced. This salt can be fed in solid form; however, it is not recommended because the salt can clog the feeding pipe; for this reason, it is recommended to dilute the salt with the reflux. In one way or another, the salt in each tray within the distillation column associates with the ethanol-water mixture creating the salting out effect. This effect is a phenomenon that consists in modifying the thermodynamic behavior of a mixture consisting of a solute in the aqueous phase by adding salt [38,39,40].
In the case of our column, the feed of the ethanol-water mixture enters as saturated vapor; in addition, there is a reflux current, through which the salt is also fed, whose compositions for the conditions described in Table 1 are x 3 = 0.0648 on a free salt basis. This salt composition was calculated using Equation (4):
F s = 0.0648 × L × F × z
where 0.0648 is the composition of the salt; L is the reflux; F is the mixed feed; and z is the ethanol concentration.
Figure 1 shows the diagram of the column used in the simulation of extractive distillation of the ethanol-water-CaCl 2 system, ready to incorporate the control loops (or control structures) and be evaluated. Table 1 summarizes the characteristics of the said plant; it has to produce ethanol above 99% mol of purity.

3. Simulation of Extractive Distillation Column with Salt in Stable State

The study that was taken as a basis to develop this simulation was the one developed by Llano-Restrepo et al. [16]. Who simulated the process in a stable state in the “FORTRAN” programming language with the thermodynamic model E-NRTL of absolute ethanol production at atmospheric pressure extractive distillation with the addition of 16.7 wt of CaCl 2 .
Unlike the study carried out by Llano-Restrepo et al. [16], in this work, the extractive distillation column was designed, simulated, and controlled at the pilot plant level using the Aspen Plus® v8.6 simulator and the E-NRTL thermodynamic model, useful for predicting the VLE of an azeotropic mixture.
The simulation of the technique used includes the prediction of VLE-CaCl 2 , which was compared with the experimental data provided by Nishi [32]. To make this comparison, the activity coefficients had to be estimated in the E-NRTL thermodynamic model of the ethanol-water-CaCl 2 system using the Aspen Properties® DRS, since Aspen Plus® does not have all. Table 2 shows the results and the differences ( Δ ) of the compositions in apparent and salt-free molar fractions of the experimental data and those estimated by the DRS.
Table 3 shows and compares the data reported in the literature for VLE-CaCl 2 . Experimental data were reported by Nishi (1975) [32]; Llano-Restrepo et al. [16]; and Hashemi et al. [33] vs. Aspen Plus® DRS estimates. It is worth mentioning that only the vapor phase ( y E x p , E s t ) was reported, since it is the most important product to be obtained in the distillation, in addition to the absence of salt. In Figure 2, these data are represented graphically. It can be seen that the data regressed by the DRS in the vapor phase are very good approximations with respect to the experimental and temperature data; it can be seen that at the bottom of the equilibrium, the data estimated by the DRS have a maximum error of 1.74 ° C. This bottom of the equilibrium represents the reboiler of the extractive distillation column.
In addition, the extractive distillation simulation was performed using the Aspen Plus® RadFrac block, which performs rigorous calculations for columns in ordinary distillation, extractive distillation applications, and more; solving the MESH equations. Figure 3 shows the graph of the composition profile results and Figure 4 shows the graph of the simulation temperature results. These profiles were analyzed and compared with the results reported in [16,33], obtaining a relative error of less than 5%.
When the results of the VLE-CaCl 2 and the profiles of the salt distillation column were compared and validated, the next step was to design our own extractive salt distillation column using RadFrac tool of Aspen Plus®, in order to add dynamics to our plant.

4. Dynamic Simulation of Extractive Distillation Column with Salt

Dynamic simulation is used to evaluate control strategies in the face of disturbances. To perform the dynamics of the simulation carried out in Aspen plus® v8.6, it had to be exported to Aspen Plus Dynamics® v8.4 using the “pressure-driven” option. Before exporting it, it was necessary to size the equipment. For this, the procedure proposed in the book by Treybal [41] was used and compared with the application of tray-sizing of the RadFrac block of Aspen Plus®. This application calculates the characteristics of the distillation column. The levels of the reflux tanks and the one at the base of the column were calculated using the method proposed by Luyben [42] proposing 10 min of retention or hold up when they are at 50% of the liquid level. Table 4 shows the characteristics obtained by these methods. It should be noted that the total number of trays reported in Table 4 does not consider tray number 1 and the last tray, since they are the condenser and the reboiler, respectively.
Aspen Plus Dynamics® simulations use rigorous nonlinear models. The E-NRTL thermodynamic model was used due to its ability to closely match the ethanol-water-CaCl 2 VLE experimental data. The column operates at 1 atmospheric pressure with a pressure drop of 0.016 atm per tray. The feed flow rate is 0.635 kmol/h. Aspen® tray numbering notation is used (numbered stages from reflux drum to kettle). Total condensers, partial reboilers, and trays with efficiencies of 50% are used to try to make it as real as possible, these efficiencies were estimated from the O’Connell empirical correlation. A column that produces high-purity products at both ends was considered in this article.

4.1. Analysis and Location of Temperature Sensing

In the distillation columns, the compositions are measured directly by composition analyzers by means of chromatography; however, this implies a dead time of 5 to 10 min, also, its reliability is insufficient for online control. Another way to infer its composition is from temperature measurements on the trays, as this is less expensive and faster for indirect composition control, plus it usually has a lag time constant that will depend on the specifications of the sensor to use.
In distillation columns with binary isobaric systems (constant pressure), the temperature is relative to the composition. In other words, temperature measurements provide accurate information on component concentrations; however, they must be strategically located on the column tray(s). This will depend on the products or degrees of freedom that we want to control within the salt extractive distillation column.
Thanks to temperature sensors, it is possible to effectively control its composition. Furthermore, such sensors are inexpensive and reliable, and introduce only a small measure of delay to the control loop. Therefore, temperature measurements are widely used to provide inferential control of compositions [43].
The selection of the tray to control its temperature (this is the reference to infer the concentration) is not usually a matter of dynamics but changes in the stable state; for this, a balance must be maintained between the following two rules:
  • The temperature of the tray must be insensitive to changes in the composition of the feed, in other words, the composition must be constant.
  • The temperature of the tray must be sensitive to changes in the manipulated variable.
The balance of these two rules can be conducted using different selection criteria. The three criteria used in this work are described in more detail in Luyben [44].
The first criterion used was the maximum slope criterion in the temperature profile. The development was as follows:
  • From the results of the temperature profile of the simulation in the steady state with its nominal inputs, the tray or the stage of the temperature profile that varies the most in the column must be selected. This tray or stage of the column is the location where the temperature must be measured by a sensor.
  • Figure 5 deals with the temperature profile of the ethanol-water-CaCl 2 column simulation. In Figure 6, small temperature variations are observed from the top to stage 17; therefore, stages 20 and 22 were selected because they have the maximum slopes.
The second criterion used was sensitivity. Its development is described below:
  • A small change in the manipulated variables must be made to analyze the changes in temperature, stage by stage.
  • For a d u a l e n d control strategy, it is necessary to measure the influence that exists on the temperature when making a small change in reflux (L) and heat in the reboiler (V).
  • Therefore, a small change in the reboiler heat and reflux (0.3% and 3%, respectively) is introduced and the variation in the temperature profile is analyzed. Figure 7 shows the results of the analysis of the gain matrix K 26 × 2 and as can be seen, stage 20 is more affected by changes in reflux than the rest of the column. On the other hand, between stages 18 and 22, it is affected more by changes in the reboiler heat, and stage 22 is selected because it is below the feed.
The third criterion was the SVD criterion.
  • The gain matrix K was factored into three matrices: K = USV, where K and U obtain information on the sensitivity and interaction of the sensors Figure 8. The absolute maximum value of vector U 1 shows the location of the most temperature-sensitive sensor, while the second position, vector U 2 , indicates less interaction (sensitivity) with respect to the first position, Luyben (2012). In the case of the ethanol-water-CaCl 2 distillation, the simulations indicate that the appropriate location of the sensors above the main feed inlet (which corresponds to the largest value of the vector U 1 ) is stage 20 and below the main feed input is stage 22 (maximum value of the vector U 2 ). The results for the d u a l e n d control are shown in Table 2.
Figure 8. Singular value decomposition (SVD). Analysis of singular vectors U 1 and U 2 vs. stage location.
Figure 8. Singular value decomposition (SVD). Analysis of singular vectors U 1 and U 2 vs. stage location.
Processes 10 01792 g008
In summary, the manipulation of the temperatures of trays 20 and 22 for the use of the structure of two control points, was the result of the application of the three criteria described above, taking into account that, in reality, the trays of the column are not affected immediately by either temperature or changes in the reflux flow. For this reason, in the simulation with Aspen Dynamics®, a delay of 6 s per tray is considered, giving a total of 2 min on the tray 20 control loop sensor, while a delay of 2.2 min was placed on the tray 22 sensor [42,44,45,46,47].

Disturbances

Disturbances are random alterations that affect our main system or plant variables. In the case of the extractive distillation column, the main feed flow (F) can be considered as a disturbance since it has a composition ( z F ) and a condition ( q f , steam). This flow is considered a disturbance because it affects the four output variables of the column ( y D , x B , M B , M D ). Furthermore, according to Luyben [48], variations in the composition of the feed are the main disturbances that exist in the control of a distillation column.
On the other hand, Figure 9 shows the extractive distillation column and its variables, the salt feed flow (S or F S ), which also has a composition ( z S ) and a condition ( q S , solid); it is not considered a disturbance because this flow does not disturb or affect the distillation column, as when it varies or there are changes in the main feed, whether in flow, composition, or reflux ratio, this feed also changes. That is why this feed flow (S) is considered an ideal control and has a relationship through Equation (4).
Figure 10 shows the continuous fluctuation disturbance. This perturbation is one of two perturbations used to evaluate the control structures.
This fluctuation was the result of the sum of three sinusoidal functions with different frequencies. The first simulations were performed with a +10% change in the composition of the ethanol in the input feed, then, a change of −10% was considered. The second series of tests was carried out in the same way but altered the nominal flow of the feed. As mentioned at the beginning of this paragraph, the third series of tests used continuous fluctuations; combining the oscillations in the feed flow and the composition of the ethanol, the different series of perturbations were introduced individually to analyze their behavior.

5. Control Strategy

In industrial practice, it is important to choose the right control structure for a distillation column. There is no “unique” structure for all columns, so some authors consider that each column should be treated independently [49]. As mentioned above, the development of the control strategy requires the conversion of the steady state model to a dynamic model to evaluate the effect of the main disturbances in the extractive distillation column.
Within the control strategies for conventional distillation columns, the control of one or two main variables is considered and this will depend on: (1) alterations in the composition of the feed and the final products; and (2) power consumption. These two types of controls were considered for the saline extractive distillation column since, as mentioned above, the extractive distillation column contains an ideal control by means of the relation of the Equation (4); therefore, the modeling of the saline extractive distillation column was treated as a binary distillation column because it considered a mixture of two pseudo-components: ethanol/salt—Water/salt. The input variables, also called manipulated variables (u) of the said column, are: reflux (L); the distillate (D); the heat in the reboiler ( Q N or V); the heat removed from the condenser ( Q 1 or V T ); and the bottom current (B). The output variables (y) or controlled variables are the distillate ( y D ) and residue ( x B ) products; as well as the levels in the reflux tanks ( M D ), column base ( M B ), and pressure (P or M V ), and can be seen in Figure 9. Consequently, the control for this distillation column is posed as a problem of 5 × 5 [50] (see Equation (5)).
u = L D B Q 1 Q N ; y = y D x B M D M B P
The output variables M D , M B , and P, are part of the secondary control loops, which ensure stability and operability in the process. This problem can be reduced to one of the following two options:
(i)
Control of a primary loop or 1 × 1 . This type of control is used when there is a mixture of two components and the interest is to maintain the purity of a product or a column stream. This structure is called the s i n g l e e n d or s i n g l e p o i n t c o n t r o l . The variables to manipulate for this work were L and D.
(ii)
Control of two primary loops or 2 × 2 . This type of control is used when there is a mixture of more than two components whose objective is to preserve the purity close to the desired specifications in the two product streams of the distillation column, both the distillate and the bottom [51]. Therefore, the structure called d u a l e n d is used. These control loops define and give the name to the control structure or strategy. The basic control structures to be used in this work are LV and DV.

5.1. Single End Structures

To evaluate the behavior of the column before the aforementioned disturbances, it was necessary to install the following control loops for structure L (Figure 11a):
  • Loop TC1: The temperature controller on tray 20. The flow rate of the reflux is used to regulate the temperature of this tray.
  • PC loop: Pressure controller. The condenser pressure is controlled by manipulating the flow of liquid from the condenser.
  • Loop LC1: Reflux tank level controller. The level of the reflux tank is regulated by manipulating the flow of the distillate.
  • Loop LC2: Waste tank level controller. The bottom tank level of the column is regulated by manipulating the flow of the residue.
  • FEED-2 loop: CaCl 2 feed. The salt feed is controlled by an ideal flow controller in order to fix the salt concentration close to 16.7 p/p in all of the trays of the column (to preserve the validity of the ELV estimate). The salt feed flow is calculated online using Equation (4).
Moreover, for structure D (Figure 11b), the control loops are the same, unlike loop TC1; temperature controller on tray 20. To regulate the temperature in this tray, the distillate flow is used.
The gains of the four controllers that were used in the simulation of the structures in Figure 11a,b are reported in Table 5, of which, the secondary loops are PIs and the main loop (the temperature one) is PID. The tunings were performed using the Ziegler-Nichols oscillation method.

5.2. Dual End Structures

For structures with two control points, the problem to be solved is the adjustment or control of the temperatures, since there is an interaction between the two loops. However, the condition number (CN) calculated from S V D of the general gain matrix G (Figure 8), suggests that the c o n t r o l d u a l of the temperature is possible. Normally, this structure is recommended by Luyben [51]:
  • For columns with medium and high concentration feeds;
  • When the control cannot be resolved with a control s i n g l e e n d ;
  • When a high purity is required in both products (distillate and bottom).
In Figure 12a, the L V control structure is shown and the control loops are clearly seen, which are defined as follows:
  • Loop TC1: Temperature controller on tray 20. The flow rate of the reflux is used to regulate the temperature of this tray.
  • Loop TC2: Temperature controller in tray 22. The regulation of the temperature in this tray is conducted by means of the heat from the reboiler.
  • PC loop: Pressure controller. Regulation of the flow of the liquid from the condenser controlled condenser pressure.
  • Loop LC1: Reflux tank level controller. The level of the reflux tank is regulated by manipulating the flow of the distillate.
  • Loop LC2: Waste tank level controller. The level at the bottom of the column is regulated by manipulating the bottom flow.
  • FEED-2 loop: CaCl 2 feed. The salt feed is controlled by an ideal flow controller, in order to fix the salt concentration close to 16.7 p/p in all the trays of the column (to preserve the validity of the ELV estimate). The salt feed flow is calculated online using Equation (4).
Likewise, the control loops for the D V structure (Figure 12b) are the same, unlike the TC1 loop; temperature controller on tray 20. To regulate the temperature in this tray, the distillate flow is used.
The gains of the five controllers that were used in the simulation of the structures in Figure 12a,b are shown in Table 6. Of which, the secondary loops are PIs and the main loops (the temperature ones) are PIDs. The tunings were performed using the Ziegler-Nichols oscillation method.

6. Checking and Assessment of Control Structures

The purpose is to comply with the international standard of ethanol purity; for this, it is necessary to adjust the temperature in one or two trays during the period of operation of the distillation column. For this reason, four control structures were designed and evaluated. This evaluation was carried out using common criteria, where the temperatures of the trays were measured to infer the concentration. These criteria are defined below:
  • IAE (Integral absolute errors)—this criterion tells us how close the composition in the distillate maintains the desired purity. In other words, it shows the variation of the measured signal before a reference value.
    I A E = 0 e ( t ) d t
  • ISE (Error square integral)—this error is represented in time and indicates the errors that exist in relation to the reference; therefore, its value increases because the size of the error grows. This error does not capture small errors.
    I S E = 0 e ( t ) 2 d t
  • ITAE (Integral absolute error multiplied by time)—the result is defined over the duration of the simulation. Large values of this error are interpreted as fluctuations in the system; as shown in the simulations, disturbances in the composition of the feed cause an oscillation at the limits of a new steady state.
    I T A E = 0 t e ( t ) d t
    This ITAE is defined in two ways:
    • ITAE n e e —as the algebraic subtraction of the new steady state and the time composition signal.
    • ITAE r e f —as the algebraic subtraction of the desired signal and that of the composition in time.
In addition, it should be noted that the indirect control produces an error in the stable state ( b i a s ) and is considered an unwanted new stable state ( n e e ). So, b i a s is defined as follows:
B i a s = c o m p o s i t i o n r e f c o m p o s i t i o n n e e × 10 3
The evaluation criteria must provide information on the transient response of the system to disturbances or changes in the desired value.

6.1. Simulations Implementing the Control Structures

This section shows the graphs of the performance results of the control structures of the saline extractive distillation column, assuming changes of ±10% in the composition of the feed input. Shown in these figures is the performance of controlling the mole fraction of ethanol in the distillate and residue products, the temperature(s) controlled, and the fluxes of dissolved CaCl 2 salt (see Figure 13, Figure 14, Figure 15 and Figure 16). Then, using the common error criteria, the performance of the control is evaluated; in addition, performance measurements and the comparison of energy consumption are also used (Figure 17).
Figure 13 shows the behavior of the process using the s i n g l e e n d structure, called structure L. As can be seen, the maximum composition of the distillate is 99.84% mf and the minimum is 99.80% mf. Moreover, the temperature of tray 20 can be observed, which is the temperature that is controlled to infer and regulate the composition of the distillate; its variation is less than 1 ° C . Table 7 summarizes Figure 13 and reports the maximum and minimum values, as well as the value at which it stabilizes in a given time. On the other hand, when the disturbance is oscillatory, the maximum composition reached by the distillate is 99.85% mf with a minimum composition of 99.79% mf, giving rise to a variation of less than 3 ° C . These data are important since in the next section we will talk about performance and energy saving.
In Figure 14, in the same way, the behavior of the process is reported using the structure of a control point called structure D, and it is observed that the composition of the distillate is greater than that reported in the previous structure, with a maximum composition of 99.85% mf but with a drop that is also greater (99.78% mf), and the variation in temperature is greater than 9 ° C . Table 8 summarizes Figure 14 and reports the maximum and minimum values, as well as the value at which it stabilizes in a given time.
Figure 15 shows the behaviors of the process disturbances in the feed composition in the step and sinusoidal form but now using the two control point structures, called the LV structure. In this figure, you can see the behavior of the composition in the distillate. Table 9 reports the maximum and minimum values that this behavior has; it should be noted that the maximum composition in the distillate is 99.84 % mf and with a minimum composition of 99.80 % mf. These compositions are very similar to the one reported in structure L but with a difference in the temperature variation in tray 20, which is reported with a variation of 0.5   ° C ; in addition, the temperature of tray 22 is also regulated, with a variation of 0.4   ° C between its maximum and minimum temperature value. On the other hand, the temperature variation in trays 20 and 22 due to sinusoidal disturbances are less than 0.8   ° C and 0.1 ° C , respectively.
Figure 16 and Table 10 report the behavior of the process and its maximum and minimum values using the structures of two control points, called structure DV. In this figure, you can see the behavior of the composition in the distillate, the behavior of the temperatures in tray 20 and tray 22, as well as their flows, both in the distillate and in the residue. In Table 10, a maximum value is reported in the distillate of 99.85 % mf and with a minimum value of 99.78 % mf. In addition, with maximum temperatures of 107.52 ° C and 124.37 ° C in tray 20 and tray 22, respectively.
Figure 17 reports the behavior of the reboiler using the four different structures previously described. These figures compare the behavior of each of the structures. To make a better comparison, the maximum and minimum values of the energy consumed are reported in Table 11. Based on this table, it can be summarized that for structures L and D, there is no excess energy consumption, because the temperature of the reboiler is not controlled, it remains constant; however, for structure DV, an increase in energy consumption is reported. The energy of up to more than 200% makes this structure the worst used to control this type of process.

6.2. Analysis of Error Results

The results of the evaluation of the performance of the control structures based on the integral error criteria are summarized in Table 12, Table 13 and Table 14.
Structure D generates the highest values of ISE, IAE, and ITAE, with the exception of bias, since structure L is the one that results in the highest value, meaning that structure L has a greater deviation with respect to the reference; however, structure D shows poor performance when compared to the other errors of structure L and the other control structures. On the other hand, the step disturbances in the input variables generate higher values of ITAE r e f , which indicate that the errors in the stable state affect the plant more than the action of the controller; moreover, ITAE n e e indicates that the response of the controller is oscillating, resulting in the D and DV structures having higher overshoots. The DV and LV structures, dispute the bias, ISE, IAE, and ITAE r e f errors between intermediate to low values. However, the errors of the ITA n e e of the DV structure are close to those of the D structure, which means that the action of the controller is similar. However, the LV and DV control structures have comparable trade-offs, no matter what type of disturbance is introduced.
For all the evaluated structures, the saline extractive distillation column will behave as non-linear, in which it will produce ‘bias’ compositions in the distillate, depending on whether the disturbance corresponds to an increase or a decrease in the input variables. In addition, the positive and negative step perturbation gives much smaller ITAE n e e than the ITAE r e f , which means that the steady-state error is oscillatory.
Therefore, these structures detect such controller behaviors. The results of the simulations with positive changes in the composition of the main feed yield an ISE and an ITAE r e f , as well as that of the oscillations.
Monitoring of the heat of the closed-loop reboiler reveals that under changes in composition and feed flow, the s i n g l e e n d control structures maintain an average power equal to the nominal value of 2.23 kW, due to the reboiler heat. However, for the d u a l e n d structures, the energy consumption is proportional to the changes in the composition of the feed, giving rise to large deviations in energy consumption (but they are very close to the nominal condition).
As mentioned above, different types of simulations were carried out, but the graphs only reported disturbances in the composition of the feed. Table 13 reports the performance of the control against disturbances in the feed flow of ± 10% (from 0.5715 to 0.6985). These perturbations had the same shapes as the perturbations described in the disturbances Section 4.1. The performances of these structures were evaluated using error and energy consumption criteria.
In summary, the steady-state deviations caused by disturbances in the feed flow are small compared to those caused by disturbances in the composition of the feed. The L structures all lead to errors smaller than those of the DV structure but similar to that of the LV structure due to overshoot. On the contrary, the L and D structures have lower power consumption in the reboiler compared to the LV and DV structures, because there is no interaction of a controller in the reboiler.
Finally, the control structures were simulated in the face of simultaneous sinusoidal variations in the flow and the composition of the main feed. The results of these simulations are reported in Table 14 and Figure 17; it can be seen that the control performance fell, causing a simple effect on the composition of the distillate; however, the standards of purity are still maintained. The L and LV structures improve the transient response of the composition in the distillate since they have the lowest values of the ISE and IAE indices. In addition, the heat of the reboiler is constant.
The measurement of the energy consumption of the reboiler complements the classification of the control structures, showing that the control structure’s s i n g l e e n d saves more energy than any structure’s d u a l e n d before any type of disturbance.

7. Conclusions

In this investigation, four different classical control structures were applied and evaluated, which were tuned using the Ziegler-Nichols oscillation method. The objective of using this type of structure is to provide a reference in the control of the process using a saline distillation column, regardless of the salt used, since there are very few reports on the control of this type of column due to: its high non-linearity, a second feed, and the introduction of a solid salt; moreover, it was given a realistic effect by considering the system dynamics, size, efficiency, and common disturbances.
On the other hand, from the point of view of the international purity standard for bioethanol, it can be concluded that a good result was obtained in the application of the control structures by manipulating the temperatures, since no procedure was found in the literature for the selection of trays for the inferential control of composition for this type of extractive distillation column, so it was a good exercise to test the procedure described by Luyben for binary distillation columns; in addition, the results of the simulation of the salt extractive distillation columns are reported using the interaction parameters obtained through the Aspen Plus DRS®, indicating a difference of less than 1% in error in the composition and temperature with respect to the experimental data reported by Nishi.
In general, for this case study in which four control structures (L, D, LV, and LD) were compared, the analysis of the results of the applied error criteria shows that for this process the structures of a control point are the best options, since the flow or heat in the reboiler is not manipulated; therefore, the energy is kept constant by step disturbances in the composition and the inlet flow of the main feed, causing a lower consumption of energy during the entire operation time of the column. In particular, for the reasons mentioned above, the L structure is the one with the best performance; this conclusion is confirmed by the reported results on error rates and power consumption.
This work can be used as a case study on the dehydration of ethanol as well as a simulation of the control of structures, since it meets the necessary characteristics to put it into practice.

Author Contributions

Formal analysis, C.A.T.C., J.Y.R.M. and R.P.Z.; Investigation, C.A.T.C. and M.M.G.; Methodology, S.R.C., M.C.-R., J.S.V.M. and M.G.M.E.; Resources, A.J.P.E. and G.O.T.; Software, J.Y.R.M.; Supervision, F.d.J.S.V. and A.G.R. All authors have read and agreed to the published versionof the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Greek symbols
γ activity coefficient
τ NRTL binary interaction energy parameter
Letters
E- N R T L electrolyte non-random two-liquid
G E molar excess Gibbs energy
N R T L non-random two-liquid
Rgas constant
Tsolid phase temperature
Xeffective liquid phase mole fraction
xtrue liquid phase mole fraction base on all species, molecular and ionic
Yeffective vapor phase mole fraction
yvapor phase mole fraction
Zabsolute value of ionic change
Superscripts
l c short range contribution, represented by the local composition model
Subscripts
a , a , a anion
c , c , c cation
E s t estimated value
e x p experimental value
i , j , k any species
mmolecular species

References

  1. Cantero, C.A.T.; Lopez, G.L.; Alvarado, V.M.; Escobar Jimenez, R.F.; Morales, J.Y.; Coronado, E.M. Control structures evaluation for a salt extractive distillation pilot plant: Application to bio-ethanol dehydration. Energies 2017, 10, 1276. [Google Scholar] [CrossRef]
  2. Intan Shafinas Muhammad, N.A.; Rosentrater, K. Economic Assessment of Bioethanol Recovery Using Membrane Distillation for Food Waste Fermentation. Bioengineering 2020, 7, 15. [Google Scholar] [CrossRef] [PubMed]
  3. Li, H.; Li, S. Optimization of Continuous Solid-State Distillation Process for Cost-Effective Bioethanol Production. Energies 2020, 13, 854. [Google Scholar] [CrossRef]
  4. Champagne, P. Bioethanol from agricultural waste residues. Environ. Prog. 2008, 27, 51–57. [Google Scholar] [CrossRef]
  5. Chuck-Hernández, C.; Pérez-Carrillo, E.; Heredia-Olea, E.; Serna-Saldívar, S. Sorgo como un cultivo multifacético para la producción de bioetanol en México: Tecnologías, avances y áreas de oportunidad. Rev. Mex. Ing. Química 2011, 10, 529–549. [Google Scholar]
  6. Naik, S.; Goud, V.V.; Rout, P.K.; Dalai, A.K. Production of first and second generation biofuels: A comprehensive review. Renew. Sustain. Energy Rev. 2010, 14, 578–597. [Google Scholar] [CrossRef]
  7. Zabed, H.; Sahu, J.N.; Suely, A.; Boyce, A.N.; Faruq, G. Bioethanol production from renewable sources: Current perspectives and technological progress. Renew. Sustain. Energy Rev. 2017, 71, 475–501. [Google Scholar] [CrossRef]
  8. Inal, M.; Yiǧitoǧlu, M. Production of bioethanol by immobilized Saccharomyces cerevisiae onto modified sodium alginate gel. J. Chem. Technol. Biotechnol. 2011, 86, 1548–1554. [Google Scholar] [CrossRef]
  9. Pejin, D.J.; Mojović, L.V.; Pejin, J.D.; Grujić, O.S.; Markov, S.L.; Nikolić, S.B.; Marković, M.N. Increase in bioethanol production yield from triticale by simultaneous saccharification and fermentation with application of ultrasound. J. Chem. Technol. Biotechnol. 2012, 87, 170–176. [Google Scholar] [CrossRef]
  10. Duarte, V.M.P. Recomendaciones de especificaciones técnicas para el etanol y sus mezclas (E6), y la infraestructura para su manejo en México. Minist. Fed. Coop. Econ. Desarro. 2010, 66, 37–39. [Google Scholar]
  11. Guzmán-Martínez, C.E.; Maya-Yescas, R.; Castro-Montoya, A.J.; Nápoles Rivera, F. Dynamic simulation of control systems for bioethanol reactive dehydration: Conventional and intensified case studies. Chem. Eng. Process.-Process Intensif. 2021, 159, 108238. [Google Scholar] [CrossRef]
  12. Figueroa, J.J.; Hoss Lunelli, B.; Maciel Filho, R.; Wolf Maciel, M.R. Improvements on anhydrous ethanol production by extractive distillation using ionic liquid as solvent. Procedia Eng. 2012, 42, 1016–1026. [Google Scholar] [CrossRef] [Green Version]
  13. Gil, I.D.; Gómez, J.M.; Rodríguez, G. Control of an extractive distillation process to dehydrate ethanol using glycerol as entrainer. Comput. Chem. Eng. 2012, 39, 129–142. [Google Scholar] [CrossRef]
  14. Jeong, J.S.; Jeon, H.; Ko, K.M.; Chung, B.; Choi, G.W. Production of anhydrous ethanol using various PSA (Pressure Swing Adsorption) processes in pilot plant. Renew. Energy 2012, 42, 41–45. [Google Scholar] [CrossRef]
  15. Lei, Z.; Yi, C.; Yang, B. Design, optimization, and control of reactive distillation column for the synthesis of tert-amyl ethyl ether. Chem. Eng. Res. Des. 2013, 91, 819–830. [Google Scholar] [CrossRef]
  16. Llano-Restrepo, M.; Aguilar-Arias, J. Modeling and simulation of saline extractive distillation columns for the production of absolute ethanol. Comput. Chem. Eng. 2003, 27, 527–549. [Google Scholar] [CrossRef]
  17. Pereiro, A.B.; Araújo, J.M.M.; Esperança, J.M.S.S.; Marrucho, I.M.; Rebelo, L.P.N. Ionic liquids in separations of azeotropic systems—A review. J. Chem. Thermodyn. 2012, 46, 2–28. [Google Scholar] [CrossRef]
  18. Yin, C.; Liu, G. Optimization of Solvent and Extractive Distillation Sequence Considering Its Integration with Reactor. Processes 2021, 9, 565. [Google Scholar] [CrossRef]
  19. Singh, A.; Rangaiah, G.P. Development and optimization of a novel process of double-effect distillation with vapor recompression for bioethanol recovery and vapor permeation for bioethanol dehydration. J. Chem. Technol. Biotechnol. 2019, 94, 1041–1056. [Google Scholar] [CrossRef]
  20. Zeng, L.; Li, Z. A new process for fuel ethanol dehydration based on modeling the phase equilibria of the anhydrous MgCl2+ethanol+water system. AIChE J. 2015, 61, 664–676. [Google Scholar] [CrossRef]
  21. Cardona, C.A.; Sánchez, Ó.J. Fuel ethanol production: Process design trends and integration opportunities. Bioresour. Technol. 2007, 98, 2415–2457. [Google Scholar] [CrossRef] [PubMed]
  22. Huang, H.J.; Ramaswamy, S.; Tschirner, U.W.; Ramarao, B.V. A review of separation technologies in current and future biorefineries. Sep. Purif. Technol. 2008, 62, 1–21. [Google Scholar] [CrossRef]
  23. Khalid, A.; Aslam, M.; Qyyum, M.A.; Faisal, A.; Khan, A.L.; Ahmed, F.; Lee, M.; Kim, J.; Jang, N.; Chang, I.S.; et al. Membrane separation processes for dehydration of bioethanol from fermentation broths: Recent developments, challenges, and prospects. Renew. Sustain. Energy Rev. 2019, 105, 427–443. [Google Scholar] [CrossRef]
  24. Kumar, S.; Singh, N.; Prasad, R. Anhydrous ethanol: A renewable source of energy. Renew. Sustain. Energy Rev. 2010, 14, 1830–1844. [Google Scholar] [CrossRef]
  25. Morales, J.Y.; Lopez, G.L.; Alvarado, V.M.; Cantero, C.A.; Rivera, H.R. Optimal predictive control for a pressure oscillation adsorption process for producing bioethanol. Comput. Sist. 2019, 23, 1593–1617. [Google Scholar] [CrossRef]
  26. Soares, R.B.; Pessoa, F.L.; Mendes, M.F. Dehydration of ethanol with different salts in a packed distillation column. Process. Saf. Environ. Prot. 2015, 93, 147–153. [Google Scholar] [CrossRef]
  27. Rumbo-Morales, J.Y.; Lopez-Lopez, G.; Alvarado, V.M.; Valdez-Martinez, J.S.; Sorcia-Vázquez, F.D.; Brizuela-Mendoza, J.A. Simulation and control of a pressure swing adsorption process to dehydrate ethanol. Rev. Mex. Ing. Quim. 2018, 17, 1051–1081. [Google Scholar] [CrossRef]
  28. Rumbo Morales, J.Y.; López López, G.; Alvarado Martínez, V.M.; Sorcia Vázquez, F.d.J.; Brizuela Mendoza, J.A.; Martínez García, M. Parametric study and control of a pressure swing adsorption process to separate the water-ethanol mixture under disturbances. Sep. Purif. Technol. 2020, 236, 116214. [Google Scholar] [CrossRef]
  29. Rumbo Morales, J.Y.; Brizuela Mendoza, J.A.; Ortiz Torres, G.; Sorcia Vázquez, F.d.J.; Rojas, A.C.; Pérez Vidal, A.F. Fault-Tolerant Control implemented to Hammerstein–Wiener model: Application to Bio-ethanol dehydration. Fuel 2022, 308, 121836. [Google Scholar] [CrossRef]
  30. Núñez, A.R.L.; Morales, J.Y.R.; Villalobos, A.U.S.; Cruz-Soto, J.D.L.; Torres, G.O.; Cerda, J.C.R.; Calixto-Rodriguez, M.; Mendoza, J.A.B.; Molina, Y.A.; Durán, O.A.Z.; et al. Optimization and Recovery of a Pressure Swing Adsorption Process for the Purification and Production of Bioethanol. Fermentation 2022, 8, 293. [Google Scholar] [CrossRef]
  31. Renteria-Vargas, E.M.; Zuniga Aguilar, C.J.; Rumbo Morales, J.Y.; Vazquez, F.D.J.S.; De-La-Torre, M.; Cervantes, J.A.; Bustos, E.S.; Calixto Rodriguez, M. Neural Network-Based Identification of a PSA Process for Production and Purification of Bioethanol. IEEE Access 2022, 10, 27771–27782. [Google Scholar] [CrossRef]
  32. Nishi, Y. Vapor-liquid equilibrium relations for the system accompanied by hypothetical chemical reaction containing salt. J. Chem. Eng. Jpn. 1975, 8, 187–191. [Google Scholar] [CrossRef] [Green Version]
  33. Hashemi, N.; Pazuki, G.; Vossoughi, M.; Hemmati, S.; Saboohi, Y. Application of a New Gibbs Energy Equation to Model a Distillation Tower for Production of Pure Ethanol. Chem. Eng. Technol. 2011, 34, 1715–1722. [Google Scholar] [CrossRef]
  34. Wang, X.; Du, Z.; Zhang, Y.; Wang, J.; Wang, J.; Sun, W. Optimization of Distillation Sequences with Nonsharp Separation Columns. Processes 2019, 7, 323. [Google Scholar] [CrossRef]
  35. Chen, C.C.; Britt, H.I.; Boston, J.F.; Evans, L.B. A Local Composition Model for the Excess Gibbs Energy of Aqueous Electrolyte Systems: Part I: Single Solvent, Single Completely Dissociated Electrolyte Systems. AIChE J. 1982, 28, 588–596. [Google Scholar] [CrossRef]
  36. Chen, C.C.; Evans, L.B. A local composition model for the excess Gibbs energy of aqueous electrolyte systems. AIChE J. 1986, 32, 444–454. [Google Scholar] [CrossRef]
  37. Gil, I.D.; Uyazán, A.M.; Aguilar, J.L.; Rodríguez, G.; Caicedo, L.A. Separation of ethanol and water by extractive distillation with salt and solvent as entrainer: Process Simulation. Braz. J. Chem. Eng. 2008, 25, 207–215. [Google Scholar] [CrossRef]
  38. Furter, W.F. Salt effect in distillation: A literature review II. Can. J. Chem. Eng. 1977, 55, 229–239. [Google Scholar] [CrossRef]
  39. Johnson, A.I.; Furter, W.F. Vapor-liquid equilibrium in systems containing dissolved salts. Can. J. Chem. Eng. 1965, 43, 356–357. [Google Scholar] [CrossRef]
  40. Meranda, D.; Furter, W.F. Vapor-Liquid Equilibrium Data for System: Ethanol-Water Saturated with Potassium Acetate. Can. J. Chem. Eng. 1966, 44, 298–300. [Google Scholar] [CrossRef]
  41. Treybal, R.E. Mass Transfer Operations. McGraw-Hill: New York, NY, USA, 1980; p. 858. [Google Scholar]
  42. Luyben, W. Plantwide Dynamic Simulators in Chemical Processing and Control; CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar]
  43. Luyben, W.L. Distillation Design and Control Using Aspen Simulation; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  44. Luyben, W.L. Evaluation of criteria for selecting temperature control trays in distillation columns. J. Process Control 2006, 16, 115–134. [Google Scholar] [CrossRef]
  45. Cecil, L.S. Practical Distillation Control; John Wiley & Sons: Hoboken, NJ, USA, 2012; p. 344. [Google Scholar]
  46. Plantwide Process Control. John Wiley & Sons: Hoboken, NJ, USA, 1999; p. 395.
  47. Luyben, W.L.; Chien, I.L. Design and Control of Distillation System for Separating Azeotropes; John Wiley & Sons: Hoboken, NJ, USA, 2010; p. 107. [Google Scholar]
  48. Luyben, W.L. Steady-State Energy Conservation Aspects of Distillation Column Control System Design. Ind. Eng. Chem. Fundam. 1975, 14, 321–325. [Google Scholar] [CrossRef]
  49. Hori, E.S.; Skogestad, S. Selection of Control Structure and Temperature Location for Two-Product Distillation Columns. Chem. Eng. Res. Des. 2007, 85, 293–306. [Google Scholar] [CrossRef]
  50. Skogestad, S. The Dos and Don’ts of Distillation Column Control. Chem. Eng. Res. Des. 2007, 85, 13–23. [Google Scholar] [CrossRef]
  51. Luyben, W.L. Effect of Feed Composition on the Selection of Control Structures for High-Purity Binary Distillation. Ind. Eng. Chem. Res. 2005, 44, 7800–7813. [Google Scholar] [CrossRef]
Figure 1. Ideal operating conditions of the saline extractive distillation column.
Figure 1. Ideal operating conditions of the saline extractive distillation column.
Processes 10 01792 g001
Figure 2. VLE expressed on a salt-free basis of ethanol-water-CaCl 2 system. (a) x y diagram; where the gray dotted line is the reference between the phases of the mixture. (b) T x y diagram.
Figure 2. VLE expressed on a salt-free basis of ethanol-water-CaCl 2 system. (a) x y diagram; where the gray dotted line is the reference between the phases of the mixture. (b) T x y diagram.
Processes 10 01792 g002
Figure 3. Comparison of the composition profiles of the saline extractive distillation column simulation. Aspen Plus® vs. Llano-Restrepo et al. and Hashemi et al. [16,33].
Figure 3. Comparison of the composition profiles of the saline extractive distillation column simulation. Aspen Plus® vs. Llano-Restrepo et al. and Hashemi et al. [16,33].
Processes 10 01792 g003
Figure 4. Comparison of the temperature profiles of the saline extractive distillation column simulation. Aspen Plus® vs. Llano-Restrepo et al. and Hashemi et al. [16,33].
Figure 4. Comparison of the temperature profiles of the saline extractive distillation column simulation. Aspen Plus® vs. Llano-Restrepo et al. and Hashemi et al. [16,33].
Processes 10 01792 g004
Figure 5. Temperature profiles of the salt extractive distillation column.
Figure 5. Temperature profiles of the salt extractive distillation column.
Processes 10 01792 g005
Figure 6. Analysis of the temperature change, stage to stage. Maximum slope criterion.
Figure 6. Analysis of the temperature change, stage to stage. Maximum slope criterion.
Processes 10 01792 g006
Figure 7. Sensor location for temperature control. Sensitivity analysis gains.
Figure 7. Sensor location for temperature control. Sensitivity analysis gains.
Processes 10 01792 g007
Figure 9. Scheme of the extractive distillation column and its variables.
Figure 9. Scheme of the extractive distillation column and its variables.
Processes 10 01792 g009
Figure 10. Perturbations of the sum of three sinusoidal functions in the feed flow and three sinusoidal functions in the feed composition.
Figure 10. Perturbations of the sum of three sinusoidal functions in the feed flow and three sinusoidal functions in the feed composition.
Processes 10 01792 g010
Figure 11. Structure of a s i n g l e e n d to regulate the purity of the distillate.
Figure 11. Structure of a s i n g l e e n d to regulate the purity of the distillate.
Processes 10 01792 g011
Figure 12. The structure’s d u a l e n d (to regulate the purity of the distillate).
Figure 12. The structure’s d u a l e n d (to regulate the purity of the distillate).
Processes 10 01792 g012aProcesses 10 01792 g012b
Figure 13. The behavior of structure L with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Figure 13. The behavior of structure L with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Processes 10 01792 g013
Figure 14. The behavior of structure D with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Figure 14. The behavior of structure D with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Processes 10 01792 g014
Figure 15. The behavior of the structure LV for the dual temperature control with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also the sinusoidal changes in the mole fraction (—red).
Figure 15. The behavior of the structure LV for the dual temperature control with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also the sinusoidal changes in the mole fraction (—red).
Processes 10 01792 g015
Figure 16. The behavior of structure DV for the dual temperature control with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Figure 16. The behavior of structure DV for the dual temperature control with changes in the main feed in a stepwise form of 10% mf (—black) and 10 % mf (—blue), also sinusoidal changes in the mole fraction (—red).
Processes 10 01792 g016
Figure 17. Power consumption of the reboiler using the control structures L, D, LV, and DV, with disturbances in the composition of the feed: (a) step of 10, (b) step of −10%, and (c) fluctuating sinusoidal input.
Figure 17. Power consumption of the reboiler using the control structures L, D, LV, and DV, with disturbances in the composition of the feed: (a) step of 10, (b) step of −10%, and (c) fluctuating sinusoidal input.
Processes 10 01792 g017
Table 1. Input data, specified parameters, and overall output data for the simulation of the saline extractive distillation column with CaCl 2 as a separating agent.
Table 1. Input data, specified parameters, and overall output data for the simulation of the saline extractive distillation column with CaCl 2 as a separating agent.
Input/Output Data and Specified ParametersSymbolValue
Input data and specified parameters
Number of equilibrium stagesN26
Main feed stage 21
Column efficiency (%) 50
Main feed thermodynamic state Saturated vapor
Main feed temperature ( ° C)T F , 2 94.5
Main feed molar flow rate (kmol/h)F 21 0.635
Main feed ethanol mole fractionx 1 , 21 0.2
Main feed water mole fractionx 2 , 21 0.8
Salt feed molar flow rate (kmol/h)F 2 0.528
Salt feed temperature ( ° C)T F , 2 78.3
Distillate molar flow rate (kmol/h) 0.127
Reflux ratio 6.42
Operating pressure (atm)P1
Output data
Ethanol molar flow rate in distillate (kmol/h)x 1 , 1 0.127
Water molar flow rate in distillate (kmol/h)x 2 , 1 0.0017
Condenser heat duty (kW)Q C −10.205
Reboiler heat duty (kW)Q R 2.2319
Table 2. VLE data estimated from experimental data of the ethanol-water-CaCl 2 system at 1 atm pressure and 16.7 wt% of CaCl 2 . Liquid phase apparent mole fractions: (a) ethanol; (b) water; (c) CaCl 2 ; vapor phase mole fraction: (d) ethanol; (e) water; (f) temperatures.
Table 2. VLE data estimated from experimental data of the ethanol-water-CaCl 2 system at 1 atm pressure and 16.7 wt% of CaCl 2 . Liquid phase apparent mole fractions: (a) ethanol; (b) water; (c) CaCl 2 ; vapor phase mole fraction: (d) ethanol; (e) water; (f) temperatures.
(a)(b)(c)
x E exp x E Est Δ x E x W exp x W Est Δ x E x CaCl 2 exp x CaCl 2 Est Δ x CaCl 2
0000.97360.973600.02640.02640
0.01950.019500.95340.953400.02710.02710
0.03890.038900.93310.933100.0280.0280
0.0970.09700.87270.872700.03030.03030
0.19310.193100.77250.772500.03440.03440
0.28850.288500.67320.673200.03830.03830
0.38310.383100.57470.574700.04220.04220
0.4770.47700.4770.47700.0460.0460
0.57010.570100.38010.380100.04980.04980
0.66250.662500.28390.283900.05360.05360
0.75410.754100.18850.188500.05740.05740
0.8450.84500.09390.093900.06110.06110
0.93520.935200000.06480.06480
(d)(e)(f)
y E exp y E Est Δ y E y W exp y W Est Δ y W T exp T Est Δ T ° C
000110102.8104.511.7116
0.2710.27106.90 × 10 6 0.7290.76340.034496.0097.7371.7370
0.3890.38902.49 × 10 6 0.6110.58720.023791.6091.3290.2709
0.5420.54202.41 × 10 6 0.4580.41500.042985.0084.2930.7061
0.6470.64703.13 × 10 7 0.3530.34390.009082.4582.3340.1158
0.6940.69401.25 × 10 7 0.3060.30030.005681.4581.8000.3500
0.7530.75294.70 × 10 7 0.2470.26860.021681.2582.3961.1466
0.7950.79494.16 × 10 7 0.2050.22760.022680.7082.3831.6831
0.8280.82791.41 × 10 7 0.1720.18190.009981.5082.0270.5274
0.8670.86694.02 × 10 8 0.1330.13600.003081.5081.7880.2887
0.9150.91507.86 × 10 9 0.0850.09040.005482.4081.7430.6561
0.9420.94208.14 × 10 8 0.0580.04380.014181.6581.0450.6044
11000082.6581.1231.5263
Table 3. Comparison between the data calculated by the D R S versus the experimental data of Nishi and those estimated by Llano and Hashemi: (a) vapor phase; (b) differences between the vapor phase data; (c) temperatures; (d) difference between the temperature data.
Table 3. Comparison between the data calculated by the D R S versus the experimental data of Nishi and those estimated by Llano and Hashemi: (a) vapor phase; (b) differences between the vapor phase data; (c) temperatures; (d) difference between the temperature data.
(a)(b)
Nishi DRS LlanoHashemi Δ y Exp
x Exp y Exp y Est y Est y Est DRS Llano Hashemi
00000000
0.020.2710.2710.2640.27100.0070
0.040.3890.3890.3940.39300.0050.004
0.100.5420.5420.5580.46600.0160.076
0.200.6470.6470.6510.59500.0040.052
0.300.6940.6940.7010.69400.0070
0.400.7530.7530.7430.76200.0100.009
0.500.7950.7950.7850.80000.0100.005
0.600.8280.8280.8290.84000.0010.012
0.700.8670.8670.8720.86800.0050.001
0.800.9150.9150.9160.91000.0010.005
0.900.9420.9420.9590.94100.0170.001
(c)(d)
Nishi DRS LlanoHashemi Δ T ° C Exp
T Exp T Est T Est T Est DRS Llano Hashemi
102.8104.51103.69108.51.710.895.70
96.0097.7495.9296.011.740.080.01
91.6091.3391.4187.670.270.193.93
85.0084.2985.1684.120.710.160.88
82.4582.3382.1082.620.120.350.17
81.4581.8081.3081.460.350.150.01
81.2582.481.1180.831.150.140.42
80.7082.3881.1480.881.680.440.18
81.5082.0381.2781.490.530.230.01
81.5081.7981.4781.840.290.030.34
82.4081.7481.7282.430.660.680.03
81.6581.0582.0281.190.600.370.46
Table 4. Characteristics of the column of the saline extractive distillation column.
Table 4. Characteristics of the column of the saline extractive distillation column.
CharacteristicsAspen Plus®TreybalLuyben
Number of stages242424
Tray spacing (m)0.100.15N.A.
Flooding factor ( C f )0.800.0424N.A.
Foaming factor0.800.80N.A.
Column diameter (m)0.1020.110.10
Total area (m 2 )0.00810.00950.0077
Active area per panel (m 2 )0.00650.0062N.A.
Downcomer area (m 2 )0.0008160.00083N.A.
Weir height (m)0.0740.077N.A.
Reflux drum total liquid volume (L)N.A.N.A.10
Reboiler total liquid volume (L)N.A.N.A.2
N.A. = Not Applicable
Table 5. Controller gain for the single temperature control.
Table 5. Controller gain for the single temperature control.
(a) Structure L(b) Structure D
ControllerPCLC1LC2TC1ControllerPCLC1LC2TC1
P K C 17.826.0132.61.44P K C 17.887.832.62.15
I τ I 1.57112.6I τ I 1.52.5126.1
D τ D 3.15D τ D 6.52
Table 6. Controller gain for the dual temperature control.
Table 6. Controller gain for the dual temperature control.
(a) Structure LV(b) Structure DV
ControllerPCLC1LC2TC1TC2ControllerPCLC1LC2TC1TC2
P K C 17.826.0132.61.4449.8P K C 17.826.0132.62.1549.8
I τ I 1.57112.64.5I τ I 1.57126.14.5
D τ D 3.151.12D τ D 6.521.12
Table 7. The behavior of the process using the structure L. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Table 7. The behavior of the process using the structure L. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Step +10%
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9981213.99 × 10 4 105.0240.1419560.5886680.0582371
Min0.9980143.33 × 10 4 104.2800.1388070.5321950.0507858
Stable0.9980833.61 × 10 4 104.8680.1397710.5471000.0518722
Step 10 %
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9984464.01 × 10 4 105.5410.1154320.5954860.0551621
Min0.9983203.22 × 10 4 104.7030.1118920.5317600.0474484
Stable0.9983723.60 × 10 4 104.8690.1142790.5745040.0537799
Sinusoidal
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9985015.33 × 10 4 106.0390.1602770.6119580.0573512
Min0.9979552.53 × 10 4 103.8730.0975080.4983450.0479607
Table 8. The behavior of the process using structure D. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Table 8. The behavior of the process using structure D. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Step +10%
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9983916.03 × 10 4 105.8010.1442840.6019450.0584786
Min0.9978662.83 × 10 4 102.3930.1378890.5285550.0495736
Stable0.9980933.61 × 10 4 104.8670.1397690.547100.0518725
Step 10 %
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9985544.74 × 10 4 107.5760.1162370.5914210.0563562
Min0.9980772.40 × 10 4 104.0140.1103990.5184550.0473115
Stable0.9983703.60 × 10 4 104.8690.1142810.5745010.0537792
Sinusoidal
Variable x D x B T s t 20 DBCaCl 2
Unitmf Ethanolmf Ethanol°Ckmol/hkmol/hkmol/h
Max0.9987641.52 × 10 4 109.4570.1655280.6229330.0604465
Min0.9974132.03 × 10 4 100.4210.09750750.4900460.0468456
Table 9. The behavior of the process using the structure LV. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Table 9. The behavior of the process using the structure LV. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Step +10%
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9981203.21 × 10 4 105.009124.1920.1433920.5865510.0582454
Min0.9980723.04 × 10 4 104.488123.7600.1393640.5025410.0516178
Stable0.9981083.16 × 10 4 104.868124.0330.1397930.5477330.0525258
Step 10 %
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9983974.02 × 10 4 105.281124.3020.1148380.6254630.0545912
Min0.9983313.54 × 10 4 104.724123.9110.1103210.5339030.0474512
Stable0.9983634.15 × 10 4 104.869124.0320.1142500.5738680.0531231
Sinusoidal
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9984896.05 × 10 4 105.873124.0800.1597200.6106310.0562294
Min0.9979612.20 × 10 4 104.288123.9760.09750560.5013910.0492466
Table 10. The behavior of the process using the structure DV. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Table 10. The behavior of the process using the structure DV. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Step +10%
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9985423.69 × 10 4 105.98124.0940.1453680.5975970.0609018
Min0.9978572.13 × 10 4 102.511123.6690.1375370.5279910.0512330
Stable0.9981073.16 × 10 4 104.868124.0330.1397920.5477330.0525263
Step 10 %
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9985546.33 × 10 4 107.564124.3740.1165040.5901630.0546879
Min0.9979883.36 × 10 4 103.930123.8990.1100410.5244050.0465062
Stable0.9983564.15 × 10 4 104.868124.0320.1142500.5738680.053119
Sinusoidal
Variable x D x B T s t 20 T s t 22 DBCaCl 2
Unitmf Ethanolmf Ethanol°C°Ckmol/hkmol/hkmol/h
Max0.9988911.44 × 10 4 110.100124.3810.1686010.6352680.0780738
Min0.9971951.06 × 10 4 100.729122.2890.09745160.4900380.0440028
Table 11. Energy consumption of the reboiler when using control structures. Reboiler energy consumption user the control structure. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Table 11. Energy consumption of the reboiler when using control structures. Reboiler energy consumption user the control structure. Maximum, minimum, and stabilization values in the event of disturbances in the feed composition.
Step +10%
StructureLDLVDV
UnitkWkWkWkW
Max2.231882.231882.468853.22535
Min2.231882.231882.224012.20217
Stable2.231882.231882.334742.33475
Step 10 %
StructureLDLVDV
UnitkWkWkWkW
Max2.231882.231882.239612.28785
Min2.231882.231882.005741.57068
Stable2.231882.231882.126972.12697
Sinusoidal
StructureLDLVDV
UnitkWkWkWkW
Max2.231882.231882.570115
Min2.231882.231881.903331.15684
Table 12. Summary of the error indices of the control structures before the disturbances in the composition of the main feed.
Table 12. Summary of the error indices of the control structures before the disturbances in the composition of the main feed.
Structure L
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)0.1526.39 × 10 5 0.4186.7820.1172.2318
Step (−10)−0.1385.11 × 10 5 0.3736.0390.2372.2318
Sinusoidal 7.88 × 10 5 0.4296.543 2.2318
Structure D
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)0.142 6.96 × 10 5 0.4246.6420.4932.2318
Step (−10)−0.135 5.57 × 10 5 0.3806.0100.3342.2318
Sinusoidal 1.55 × 10 4 0.5728.939 2.2318
Structure LV
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)0.127 4.66   ×   10 5 0.3595.7890.1582.3284
Step (−10)−0.128 4.13   ×   10 5 0.3375.5030.2792.1335
Sinusoidal 5.83   ×   10 5 0.3695.623 2.2263
Structure DV
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)0.128 6.09   ×   10 5 0.3875.8840.4322.3568
Step (−10)−0.121 4.94   ×   10 5 0.3525.4460.36082.1235
Sinusoidal 2.00   ×   10 4 0.6459.943 2.3411
Table 13. Summary of the error indices of the control structures before the disturbances in the feed flux main rate.
Table 13. Summary of the error indices of the control structures before the disturbances in the feed flux main rate.
Structure L
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)−0.005 4.61   ×   10 6 0.0520.3380.4442.2318
Step (−10)−0.011 3.06   ×   10 6 0.0530.4460.1702.2318
Sinusoidal 3.39   ×   10 6 0.0821.149 2.2318
Structure D
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)0.002 1.57   ×   10 5 0.0830.4310.3612.2318
Step (−10)−0.013 1.18   ×   10 5 0.0880.7310.3452.2318
Sinusoidal 5.67   ×   10 5 0.3194.767 2.2318
Structure LV
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)−0.019 1.80   ×   10 6 0.0590.8310.1312.3547
Step(-10)0.007 1.81   ×   10 6 0.0490.5430.2552.1170
Sinusoidal 2.50   ×   10 6 0.0691.003 2.2427
Structure DV
PerturbationBias × 10 3 ISEIAEITAE s e t ITAE n s s Average Heat Duty r e b o i l e r (kW)
Step (+10)−0.016 2.72   ×   10 5 0.1361.1110.4912.4041
Step (−10)0.021 1.72   ×   10 5 0.1090.9960.4562.0985
Sinusoidal 9.62   ×   10 5 0.4166.170 2.3017
Table 14. Summary of the error rates of the control structures in the face of simultaneous disturbances in the composition and feed flow.
Table 14. Summary of the error rates of the control structures in the face of simultaneous disturbances in the composition and feed flow.
Structure L
PerturbationISEIAEITAE s e t Average Heat Duty r e b o i l e r (kW)
Sinusoidal 8.212   ×   10 5 0.4396.6652.2318
Structure D
PerturbationISEIAEITAE s e t Average Heat Duty r e b o i l e r (kW)
Sinusoidal 2.044   ×   10 4 0.6209.0302.2318
Structure LV
PerturbationISEIAEITAE s e t Average Heat Duty r e b o i l e r (kW)
Sinusoidal 6.335   ×   10 5 0.3825.9432.2367
Structure DV
PerturbationISEIAEITAE s e t Average Heat Duty r e b o i l e r (kW)
Sinusoidal 2.965   ×   10 4 0.72810.2842.4104
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Torres Cantero, C.A.; Pérez Zúñiga, R.; Martínez García, M.; Ramos Cabral, S.; Calixto-Rodriguez, M.; Valdez Martínez, J.S.; Mena Enriquez, M.G.; Pérez Estrada, A.J.; Ortiz Torres, G.; Sorcia Vázquez, F.d.J.; et al. Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol. Processes 2022, 10, 1792. https://doi.org/10.3390/pr10091792

AMA Style

Torres Cantero CA, Pérez Zúñiga R, Martínez García M, Ramos Cabral S, Calixto-Rodriguez M, Valdez Martínez JS, Mena Enriquez MG, Pérez Estrada AJ, Ortiz Torres G, Sorcia Vázquez FdJ, et al. Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol. Processes. 2022; 10(9):1792. https://doi.org/10.3390/pr10091792

Chicago/Turabian Style

Torres Cantero, Carlos Alberto, Ricardo Pérez Zúñiga, Mario Martínez García, Silvia Ramos Cabral, Manuela Calixto-Rodriguez, Jorge Salvador Valdez Martínez, Mayra Guadalupe Mena Enriquez, Abraham Jashiel Pérez Estrada, Gerardo Ortiz Torres, Felipe de J. Sorcia Vázquez, and et al. 2022. "Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol" Processes 10, no. 9: 1792. https://doi.org/10.3390/pr10091792

APA Style

Torres Cantero, C. A., Pérez Zúñiga, R., Martínez García, M., Ramos Cabral, S., Calixto-Rodriguez, M., Valdez Martínez, J. S., Mena Enriquez, M. G., Pérez Estrada, A. J., Ortiz Torres, G., Sorcia Vázquez, F. d. J., García Rebolledo, A., & Rumbo Morales, J. Y. (2022). Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol. Processes, 10(9), 1792. https://doi.org/10.3390/pr10091792

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