# Feasibility of Batch Reactive Distillation with Equilibrium-Limited Consecutive Reactions in Rectifier, Stripper, or Middle-Vessel Column

Authors: T. Lukács; C. Stéger; E. Rév; M. Meyer; Z. Lelkes
Journal: International Journal of Chemical Engineering
 (2011)

---

## Abstract

A general overall feasibility methodology of batch reactive distillation of multireaction systems is developed to study all the possible configurations of batch reactive distillation. The general model equations are derived for multireaction system with any number of chemical equilibrium-limited reactions and for any number of components. The present methodology is demonstrated with the detailed study of the transesterification of dimethyl carbonate in two reversible cascade reactions in batch reactive distillation process. Pure methanol is produced as distillate, and pure diethyl carbonate is produced at the bottom simultaneously in middle-vessel column;  in each section, continuous feeding of ethanol is necessary. The results of feasibility study are successfully validated by rigorous simulations.


---

## Body

## 1. Introduction

Reactive distillation is a hybrid process integrating separation and reaction in a single unit. The main advantages of reactive distillation compared to the sequential processes include reduced investment and operating costs, reduced utility consumption, and higher conversion and selectivity [1]. Nevertheless, reactive distillation is not advantageous in every case, and systematic methods are needed to examine its feasibility for particular reaction systems.Design of complex processes usually starts with feasibility analysis including computation of limiting flows and minimum stage numbers as well. Quick, simple, and reliable method is needed to find the best feasible configurations and parameter regions applicable to produce the desired products.Several methodologies for the preliminary study of reactive distillation have been published in the past decades; short overview is in  Stéger et al. [2]. Many of them are based on graphical techniques providing a basis for developing conceptual designs (e.g., in the articles of  [3–6]). In spite of the efforts, their application is still limited because of the lack of a systematic design method applicable in general case, especially for complex column configurations [2, 7].All the possible configurations of batch reactive distillation with equilibrium-limited single reactions in rectifier with reactive boiler can be analyzed with the reliable, overall feasibility method suggested by Stéger et al. [2]. Only configurations with at most one feed are discussed there but the method can be extended for several entrainer feeds [7].In the present work, we extend our general feasibility method for any number of equilibrium-limited reactions in middle-vessel column, stripper, and rectifier. The version extended for any number of reactions and components of general model equation presented in Stéger et al. [2] for the reactive zone is derived here. Middle-vessel column, its upper part as a rectifier, and the lower part as a stripper are all included in the analysis.Multireaction system containing two or more reactions and multicomponent systems containing four or more components have already been published in the literature, but none of them deals with their overall feasibility analysis for batch processes, examining all the possible configuration.In the present work, we present the methodology of an overall feasibility study of any multicomponent and multireaction system for batch processes and any configuration of middle vessel column, stripper, and rectifier. The presented methodology can also be used to continuous processes after applying some modifications [8].The method is demonstrated in studying the feasibility of transesterification of dimethyl carbonate (DMC) in two cascade reversible reactions among five components with an intermediate methyl ethyl carbonate (MEC)  [9] in batch reactive distillation procedure. Although modeling the production of DEC in a continuous reactive distillation process has been published by Luo and Xiao [9], its overall feasibility study has not yet been published.


## 2. The Feasibility Method

### 2.1. Single-Reaction Systems

During the feasibility study, a simplified model is used with the common assumptions like instantaneous equilibrium-limited reaction in the liquid phase, neglected reaction enthalpy, and liquid hold-up on the stages, instantaneous vapor-liquid equilibrium with neglected catalyst effect [2]. Due to the simplifying assumptions, the feasible configurations can be found in a short time, and the feasible domain of the main operation parameters can also be determined.Batch extractive distillation columns can be divided into three different zones: rectifying or stripping section, extractive section, and still vessel [10–12]. The rectifying section is located below the distillate and the stripping section above the bottom product. Any extractive section is located between the middle-vessel and some feed. This terminology is applicable to batch rectifier, stripper, and middle vessel column as is shown in Figure 1. Applying heterogeneous catalyst, any section can be either reactive or nonreactive. Whereas, applying homogeneous catalyst, the sections below the catalyst feed must be reactive.Figure 1
Middle vessel column.There are several possible configurations of batch reactive distillation in a batch reactive rectifier or in a batch reactive stripper considering a reactive vessel and at most one feed, depending on the place of the reaction and feed locations. For example, a batch rectifier with at most one entrainer feed and a single contiguous reactive section including the still may consist of the still and the following alternatives listed bottom up in each version (1) a nonreactive rectifying section, (2) a reactive rectifying section, (3) a reactive rectifying section and a nonreactive rectifying section, (4) a nonreactive extractive section below an entrainer feed, and a nonreactive rectifying section, (5) a reactive extractive section below an entrainer feed, and a reactive rectifying section, (6) a reactive extractive section below an entrainer feed, and a nonreactive rectifying section, (7) a reactive extractive section below an entrainer feed, a reactive rectifying section and a nonreactive rectifying section, or (8) a reactive extractive section below a nonreactive extractive section below an entrainer feed, and a nonreactive rectifying section. Locating the feed directly to the top of the column or to the still vessel [13] is not even considered in the list above. There may be more feeds; however, and the number of contiguous reactive sections can also be higher if heterogeneous catalysis is applied. Thus, a large number of configurations can be conceived for a middle-vessel column because it is a combination of a rectifier and a stripper. One aim of the preliminary study is to select the feasible configurations.The methodology published in Stéger et al. [2] is further developed here to cope with this combinatorial diversity. Instead of considering all the configurations separately, the possible composition profiles in the potential column sections are analyzed in the proposed feasibility study. No column length constraint is taken into account, namely, infinite length is considered. The ensemble of possible profiles is constrained by product purity specifications and by the charge composition.Due to assuming instantaneous equilibrium-limited reactions, the column profiles can be represented in a transformed space. With the transformation,(1)ãi=ai-νiνrefaref,
proposed by Stéger et al. [2], where “ai” can be any liquid, vapor, feed, or product component flow rate, the reactive profiles are brought to the same form used for the nonreactive sections by Lelkes et al. [10] and Varga et al. [12] (Table 1).Table 1
Model equations of column sections for single-reaction systems.
SectionNonreactiveReactiveRectifyingdxidh=(VL)(yi-yi*) whereyi=(LV)xi+(DV)xD,i [13]dl̃idh=ṽi-ṽi* whereṽi=l̃i+d̃i  [2]Extractive  (in a rectifier)dxidh=(L′V)(yi-yi*) whereyi=(L′V)xi+(DV)xD,i-(FdV)xfd,i [13]dl̃i′dh=ṽi-ṽi* where ṽi=l̃i′+d̃i-f̃d,i  [2]Strippingdxidh=(VL)(yi-yi*) whereyi=(LV)xi-(WV)xW,i [11]dl̃idh=ṽi-ṽi* where ṽi=l̃i-w̃i [2]Extractive  (in a stripper)dxidh=(L′V)(yi-yi*) whereyi=(VL′)xi-(WV)xW,i+(FwV)xfw,i [11]dl̃i′dh=ṽi-ṽi* where  ṽi=l̃i′-w̃i+f̃w,i [2]The following transformation, proposed by Doherty and Buzad [14], is used when computing a vessel composition path:(2)Xi=xi-(νi/νref)⋅xref1-(νT/νref)⋅xref,Yi=yi-(νi/νref)⋅yref1-(νT/νref)⋅yref,
where νi is stoichiometric coefficient of component i in the reaction, νT is the sum (total) of all the stoichiometric coefficients in the reaction, and the index ref refers to a reference component. Using this transformation together with the transformation shown in   (1), proposed by Stéger et al. [2], the vessel path equation is brought to a form similar to that used for the nonreactive vessel by Lelkes et al. [10] and Varga et al. [12]. The equations have the same form if the total change of moles in the reaction is zero. This is shown in Table 2. Here U denotes the number of moles in the vessel, ξ  denotes reaction coordinate, and the sum of the reaction coordinates is made along the column. Indices D and W refer to distillate and weight product, respectively; indices fd and fw refer to feed streams entering the rectifier part and the stripper part, respectively.Table 2
Model equations of vessel path for single-reaction systems.
Nonreactivedxidt=D(xi-xD,i)+W(xi-xW,i)-Fd(xfd,i-xi)-Fw(xfw,i-xi)Reactive,νT=0dXidt=D(Xi-XD,i)+W(Xi-XW,i)-Fd(Xfd,i-Xi)-Fw(Xfw,i-Xi)Reactive,νT≠0Ũ=U(1-(νTνref)xref)  andŨ(dXidt)=D̃·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xD,ref)-XD,i)+W̃·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xW,ref)-XW,i)-F̃d·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xfd,ref)-Xfd,i)-F̃w·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xfw,ref)-Xfw,i)-Xi·(1-(νTνref)·xref)·νT·∫column(dξdt)dh+U·Xi·(νTνref)·(dxrefdt)In order to analyze all the potential column sections, the section profiles and the still path are modeled with the equations presented in Table1. A configuration is feasible if the specified products can be produced with a reasonable recovery at some operation parameter values (reflux ratio, entrainer feed ratio, etc.), that is, in a reasonable time while the vessel composition continuously changes. In each moment of this process, the column profile must connect the vessel with the specified product compositions. In order to maintain such a connection, the individual section profiles must meet at the section junctions, and the profiles in those sections neighboring the vessel must meet the vessel composition in some sense. This may not only mean simply an intersection of the profiles (and intersection of some profiles with the vessel path) but may also mean a more complicated criterion for reactive sections and reactive vessel path. The conditions of the feasibility of the junction points are presented in detail in Stéger et al. [2] and summarized briefly in Table 3.Table 3
Conditions of feasibility at column section junctions [2].
Contact of two nonreactive sectionsIntersection of the compositions profilesContact of two reactive sectionsIntersection of the compositions profilesA nonreactive section above a reactive sectionIntersection of the dew-point profile of the upper section with the composition profile of the lower sectionA reactive section above a nonreactive sectionIntersection of the compositions profilesContact between a nonreactive section and a nonreactive vesselIntersection of the composition profile and the vessel pathContact between a reactive section and a reactive vesselIntersection of the composition profile and the vessel pathContact between a nonreactive stripping section and a reactive stillIntersection of the composition profile and the vessel pathContact between a nonreactive rectification section and a reactive vesselIntersection of the dew-point profile and the reactive vessel pathIf a nonreactive section is directly situated above a reactive section, then the dew-point profile of the upper section must intersect the liquid composition profile of the lower section. If a nonreactive section is directly situated above a reactive vessel, then the dew-point profile of this section must intersect the reactive vessel path. The dew point profile is the curve of imaginary liquid compositions in equilibrium with the vapor compositions formed in the section. Note that the vapor composition profile is computed from the liquid composition profile by material balance, and; thus, the dew point profile is different from the liquid composition profile.The feasibility criterion is the intersection of the liquid composition profiles or a profile and the vessel path in all the other cases.


### 2.2. Multiple Reaction Systems

The transformation of component flow rates is extended to multireaction systems, as(3)ãn=an-A⋅B-1⋅an(ref),
where index n refers to stage n, matrix A is an C×R matrix of stoichiometric coefficients of the components in the reactions, B is an R×R matrix of stoichiometric coefficients of the reference components, C is the number of components, R is the number of reactions, an is a component flow rate vector of all the C components, whereas an(ref) is a component flow rate vector of the R reference components:(4)an≡[an,1⋮an,i⋮an,Nc]an(ref)≡[an,1(ref)⋮an,i(ref)⋮an,R(ref)].
The column profile equations given in Table 1 remain unchanged, due to the transformation.The form of the reactive vessel path equation become as follows.IfνT=0, then(5)dXidt=D(Xi-XD,i)+W(Xi-XW,i)-Fd(Xi-Xdf,i)-Fw(Xi-Xwf,i),
else (νT≠0)(6)Ũ=U(1-νTTN(ref)-1x(ref)),ŨddtXi=D̃(Xi1-bTx(ref)1-bTxD(ref)-XD,i)+W̃(Xi1-bTx(ref)1-bTxW(ref)-XW,i)-F̃d(Xdf,i1-bTx(ref)1-bTxdf(ref)-Xdf,i)-F̃w(Xwf,i1-bTx(ref)1-bTxwf(ref)Xwf,i)-νTTXi(1-bTx(ref))∫columndξdtdh+UXibTddtx(ref),
where νTT is the transpose of vector νT of total stoichiometric coefficients, of length R, N(ref) is an R×R matrix of stoichiometric coefficients of the reference components in each reaction, and bT is the transpose of vector b defined by (7)bT≡B-1νTT.
The model equations and the vessel path of multireaction systems are derived in details in Appendices A and B, respectively.


## 2.1. Single-Reaction Systems

During the feasibility study, a simplified model is used with the common assumptions like instantaneous equilibrium-limited reaction in the liquid phase, neglected reaction enthalpy, and liquid hold-up on the stages, instantaneous vapor-liquid equilibrium with neglected catalyst effect [2]. Due to the simplifying assumptions, the feasible configurations can be found in a short time, and the feasible domain of the main operation parameters can also be determined.Batch extractive distillation columns can be divided into three different zones: rectifying or stripping section, extractive section, and still vessel [10–12]. The rectifying section is located below the distillate and the stripping section above the bottom product. Any extractive section is located between the middle-vessel and some feed. This terminology is applicable to batch rectifier, stripper, and middle vessel column as is shown in Figure 1. Applying heterogeneous catalyst, any section can be either reactive or nonreactive. Whereas, applying homogeneous catalyst, the sections below the catalyst feed must be reactive.Figure 1
Middle vessel column.There are several possible configurations of batch reactive distillation in a batch reactive rectifier or in a batch reactive stripper considering a reactive vessel and at most one feed, depending on the place of the reaction and feed locations. For example, a batch rectifier with at most one entrainer feed and a single contiguous reactive section including the still may consist of the still and the following alternatives listed bottom up in each version (1) a nonreactive rectifying section, (2) a reactive rectifying section, (3) a reactive rectifying section and a nonreactive rectifying section, (4) a nonreactive extractive section below an entrainer feed, and a nonreactive rectifying section, (5) a reactive extractive section below an entrainer feed, and a reactive rectifying section, (6) a reactive extractive section below an entrainer feed, and a nonreactive rectifying section, (7) a reactive extractive section below an entrainer feed, a reactive rectifying section and a nonreactive rectifying section, or (8) a reactive extractive section below a nonreactive extractive section below an entrainer feed, and a nonreactive rectifying section. Locating the feed directly to the top of the column or to the still vessel [13] is not even considered in the list above. There may be more feeds; however, and the number of contiguous reactive sections can also be higher if heterogeneous catalysis is applied. Thus, a large number of configurations can be conceived for a middle-vessel column because it is a combination of a rectifier and a stripper. One aim of the preliminary study is to select the feasible configurations.The methodology published in Stéger et al. [2] is further developed here to cope with this combinatorial diversity. Instead of considering all the configurations separately, the possible composition profiles in the potential column sections are analyzed in the proposed feasibility study. No column length constraint is taken into account, namely, infinite length is considered. The ensemble of possible profiles is constrained by product purity specifications and by the charge composition.Due to assuming instantaneous equilibrium-limited reactions, the column profiles can be represented in a transformed space. With the transformation,(1)ãi=ai-νiνrefaref,
proposed by Stéger et al. [2], where “ai” can be any liquid, vapor, feed, or product component flow rate, the reactive profiles are brought to the same form used for the nonreactive sections by Lelkes et al. [10] and Varga et al. [12] (Table 1).Table 1
Model equations of column sections for single-reaction systems.
SectionNonreactiveReactiveRectifyingdxidh=(VL)(yi-yi*) whereyi=(LV)xi+(DV)xD,i [13]dl̃idh=ṽi-ṽi* whereṽi=l̃i+d̃i  [2]Extractive  (in a rectifier)dxidh=(L′V)(yi-yi*) whereyi=(L′V)xi+(DV)xD,i-(FdV)xfd,i [13]dl̃i′dh=ṽi-ṽi* where ṽi=l̃i′+d̃i-f̃d,i  [2]Strippingdxidh=(VL)(yi-yi*) whereyi=(LV)xi-(WV)xW,i [11]dl̃idh=ṽi-ṽi* where ṽi=l̃i-w̃i [2]Extractive  (in a stripper)dxidh=(L′V)(yi-yi*) whereyi=(VL′)xi-(WV)xW,i+(FwV)xfw,i [11]dl̃i′dh=ṽi-ṽi* where  ṽi=l̃i′-w̃i+f̃w,i [2]The following transformation, proposed by Doherty and Buzad [14], is used when computing a vessel composition path:(2)Xi=xi-(νi/νref)⋅xref1-(νT/νref)⋅xref,Yi=yi-(νi/νref)⋅yref1-(νT/νref)⋅yref,
where νi is stoichiometric coefficient of component i in the reaction, νT is the sum (total) of all the stoichiometric coefficients in the reaction, and the index ref refers to a reference component. Using this transformation together with the transformation shown in   (1), proposed by Stéger et al. [2], the vessel path equation is brought to a form similar to that used for the nonreactive vessel by Lelkes et al. [10] and Varga et al. [12]. The equations have the same form if the total change of moles in the reaction is zero. This is shown in Table 2. Here U denotes the number of moles in the vessel, ξ  denotes reaction coordinate, and the sum of the reaction coordinates is made along the column. Indices D and W refer to distillate and weight product, respectively; indices fd and fw refer to feed streams entering the rectifier part and the stripper part, respectively.Table 2
Model equations of vessel path for single-reaction systems.
Nonreactivedxidt=D(xi-xD,i)+W(xi-xW,i)-Fd(xfd,i-xi)-Fw(xfw,i-xi)Reactive,νT=0dXidt=D(Xi-XD,i)+W(Xi-XW,i)-Fd(Xfd,i-Xi)-Fw(Xfw,i-Xi)Reactive,νT≠0Ũ=U(1-(νTνref)xref)  andŨ(dXidt)=D̃·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xD,ref)-XD,i)+W̃·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xW,ref)-XW,i)-F̃d·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xfd,ref)-Xfd,i)-F̃w·(Xi·(1-(νT/νref)·xref)(1-(νT/νref)·xfw,ref)-Xfw,i)-Xi·(1-(νTνref)·xref)·νT·∫column(dξdt)dh+U·Xi·(νTνref)·(dxrefdt)In order to analyze all the potential column sections, the section profiles and the still path are modeled with the equations presented in Table1. A configuration is feasible if the specified products can be produced with a reasonable recovery at some operation parameter values (reflux ratio, entrainer feed ratio, etc.), that is, in a reasonable time while the vessel composition continuously changes. In each moment of this process, the column profile must connect the vessel with the specified product compositions. In order to maintain such a connection, the individual section profiles must meet at the section junctions, and the profiles in those sections neighboring the vessel must meet the vessel composition in some sense. This may not only mean simply an intersection of the profiles (and intersection of some profiles with the vessel path) but may also mean a more complicated criterion for reactive sections and reactive vessel path. The conditions of the feasibility of the junction points are presented in detail in Stéger et al. [2] and summarized briefly in Table 3.Table 3
Conditions of feasibility at column section junctions [2].
Contact of two nonreactive sectionsIntersection of the compositions profilesContact of two reactive sectionsIntersection of the compositions profilesA nonreactive section above a reactive sectionIntersection of the dew-point profile of the upper section with the composition profile of the lower sectionA reactive section above a nonreactive sectionIntersection of the compositions profilesContact between a nonreactive section and a nonreactive vesselIntersection of the composition profile and the vessel pathContact between a reactive section and a reactive vesselIntersection of the composition profile and the vessel pathContact between a nonreactive stripping section and a reactive stillIntersection of the composition profile and the vessel pathContact between a nonreactive rectification section and a reactive vesselIntersection of the dew-point profile and the reactive vessel pathIf a nonreactive section is directly situated above a reactive section, then the dew-point profile of the upper section must intersect the liquid composition profile of the lower section. If a nonreactive section is directly situated above a reactive vessel, then the dew-point profile of this section must intersect the reactive vessel path. The dew point profile is the curve of imaginary liquid compositions in equilibrium with the vapor compositions formed in the section. Note that the vapor composition profile is computed from the liquid composition profile by material balance, and; thus, the dew point profile is different from the liquid composition profile.The feasibility criterion is the intersection of the liquid composition profiles or a profile and the vessel path in all the other cases.


## 2.2. Multiple Reaction Systems

The transformation of component flow rates is extended to multireaction systems, as(3)ãn=an-A⋅B-1⋅an(ref),
where index n refers to stage n, matrix A is an C×R matrix of stoichiometric coefficients of the components in the reactions, B is an R×R matrix of stoichiometric coefficients of the reference components, C is the number of components, R is the number of reactions, an is a component flow rate vector of all the C components, whereas an(ref) is a component flow rate vector of the R reference components:(4)an≡[an,1⋮an,i⋮an,Nc]an(ref)≡[an,1(ref)⋮an,i(ref)⋮an,R(ref)].
The column profile equations given in Table 1 remain unchanged, due to the transformation.The form of the reactive vessel path equation become as follows.IfνT=0, then(5)dXidt=D(Xi-XD,i)+W(Xi-XW,i)-Fd(Xi-Xdf,i)-Fw(Xi-Xwf,i),
else (νT≠0)(6)Ũ=U(1-νTTN(ref)-1x(ref)),ŨddtXi=D̃(Xi1-bTx(ref)1-bTxD(ref)-XD,i)+W̃(Xi1-bTx(ref)1-bTxW(ref)-XW,i)-F̃d(Xdf,i1-bTx(ref)1-bTxdf(ref)-Xdf,i)-F̃w(Xwf,i1-bTx(ref)1-bTxwf(ref)Xwf,i)-νTTXi(1-bTx(ref))∫columndξdtdh+UXibTddtx(ref),
where νTT is the transpose of vector νT of total stoichiometric coefficients, of length R, N(ref) is an R×R matrix of stoichiometric coefficients of the reference components in each reaction, and bT is the transpose of vector b defined by (7)bT≡B-1νTT.
The model equations and the vessel path of multireaction systems are derived in details in Appendices A and B, respectively.


## 3. Feasibility Study

Application of the equations and criteria presented in Section1 is demonstrated here by performing the feasibility study of transesterification of dimethyl carbonate (DMC) with ethanol (EtOH) in a batch extractive distillation process.Transesterification of DMC consists of two cascade reversible reactions with an intermediate methyl ethyl carbonate (MEC) [9]. The main product is diethyl carbonate (DEC), and methanol (MeOH) is a by-product, but both of them are considered as desired products in our case as follows:(8)H3COCOOCH3DMC+C2H5OHEtOH⟺K1H3COCOOC2H5MEC+CH3OHMeOHH3COCOOC2H5MEC+C2H5OHEtOH⟺K2H5C2OCOOC2H5DEC+CH3OHMeOH.The equilibrium-limited reactions are characterized with equilibrium constantsK1=2.33 and K2=0.26 obtained by our laboratory measurements and are approximately constant in the studied temperature range. The equilibrium constants of the consecutive reactions are described by the following equations: (9)K1=[MEC]⋅[MeOH][DMC]⋅[EtOH]=2.33,K2=[DEC]⋅[MeOH][MEC]⋅[EtOH]=0.26.Total conversion cannot be reached because the reactions are equilibrium limited, but the productivity can be increased by removing the reaction products from the reaction space. The objective of this study is to examine the feasibility of simultaneously producing MeOH and DEC in middle-vessel column. Both product purities are specified to 0.95 (xd≥0.95,xw≥0.95).Vapor-liquid equilibrium is modeled with the modified Raoult-Dalton equation(10)yiP=γixipio,
where P is the actual pressure, γi is activity coefficient of component i in the liquid mixture, and pio is vapor pressure of pure component i. The activity coefficients are computed with 2-parameter Wilson model [15], and the vapor pressures with the Antoine equation in the form of(11)ln⁡pio=A-BT+C.
The model parameters are collected in Tables 4 and 5. The studied system contains 3 binary azeotropes, and there is no ternary azeotrope, as is shown in Table 6 [9].Table 4
Parameters of theWilson model (cal/mol) [7].
(λ12-λ11)(λ22–λ21)MeOH-EtOH135.8113−132.0576MeOH-DMC771.98202.44MeOH-MEC684.591203.922MeOH-DEC364.177411.055EtOH-DMC522.235450.581EtOH-MEC709.523164.462EtOH-DEC626.606331.781DMC-MEC323.814−311.777DMC-DEC−474.6591000.0MEC-DEC80.9023−83.1232Table 5
TheAntoine parameters. pio is measured in (mmHg), and T is measured in [K] [7].
ABCMEC16.822453376.596−49.461MeOH18.587523626.548−34.29EtOH18.911983800.417−41.67DEC15.561092817.835−84.304DMC16.829363253.553−44.25Table 6
VLE singularities in the system of MeOH, EtOH, DMC, MEC, and DEC.
SystemComposition (mole fractions)Temperature (°C)Type of singularity*DMC-MeOH(0.123; 0.877)63.7UNMeOH164.5SDMC-EtOH(0.31; 0.69)74.9SMEC-EtOH(0.05; 0.95)78.3SEtOH178.5SDMC190.0SMEC1109.2SDEC1126.0SN*SN: stable node; UN: unstable node; S: saddle.Feasibility study aims at determining a region of vessel compositions, called “feasible region”, from where the specified products can be reached by column profiles and the domain of column parameters such as reflux ratioR, reboil ratio S, and feed ratio(s) Fd/V and/or Fw/L. The potential combinations of column sections appropriate to produce the desired products constitute the candidate column structures.


### 3.1. Total Reflux/Reboil Analysis

Running the system with no product removal is the limit case of any batch distillation process and is usually analyzed before finite reflux ratios and/or finite reboil ratios are studied because total reflux and/or total reboil involves the highest reachable product purity in conventional distillation systems. It is not always so with extractive and reactive distillation, but performing a total reflux/reboil analysis is worth because valuable information can be obtained this way on the behavior of the studied system.


#### 3.1.1. Nonreactive Sections

Nonreactive rectifying profiles and stripping profiles run along nonreactive residue curves in case of total reflux or total reboil, respectively. Thus, residue curves are analyzed.The studied system consists of five components. Thus, visualization of the residue curves is difficult. With the knowledge of the boiling points of the azeotropes and the pure components, the singular points and, thus, the shape of the residue curve map of the nonreactive system can be determined. In order to make the visualization easier, the system is projected to a plane (Figure2). According to this projected network, six-saddle points, one stable-node (SN), and one unstable node (UN) are found.Figure 2
The nonreactive system unfolded in 2D.The desired products are DEC as bottom product  and MeOH as distillate. Since the DEC vertex is the stable node of the system, production of pure DEC in the bottom is feasible both in batch and continuous processes. Just the contrary, pure MeOH cannot be produced at the top of a nonreactive rectifying section because the MeOH vertex is a saddle point (S) of the system. The top product would be the unstable node, that is, the DMC-MeOH azeotrope.Such a qualitative analysis of the nonreactive extractive profiles is rather difficult because their behavior cannot be visualized in a similar and clear way. Accumulation points (stable nodes) of the extractive profile bundles depend on process parameters like reflux ratioR or reboil ratio S and feed ratio F/V or F/L.


#### 3.1.2. Reactive Sections

Due to the instantaneous reaction equilibrium, dimensionality of the reactive material system is reduced to a 2D hypersurface (of the equilibrium points) in the originally 4D linear space of the mole fractions. (There are five mole fractions, their sum is unity, and; thus, only four independent mole fractions are considered.) This curved hypersurface is referred to as “the reactive surface”. Due to applying the proposed transformation, this reactive surface, represented with reactive residue curves, is projected to a 2D plane as is shown in Figure3.Figure 3
Residue curves in the transformed reactive space.The reactive space contains the same unstable node (the locus of the DMC-MeOH azeotrope) and the same stable node (DEC vertex) as in the case of the nonreactive system. Thus, without extractive section, the DMC-MeOH azeotrope and pure DEC are withdrawn as distillate and bottom products, respectively, with infinite number of theoretical stages.There is a narrow bundle of reactive residue curves running so near the MeOH vertex that they intersect the composition region of assigned purity. Thus, the assigned purity can be reached with a finite number of rectifying stages. The feasible region is the set of compositions covered by that part of this bundle which connect vertex DEC with the targeted region around vertex MeOH.Unfortunately, the targeted MeOH purity can be reached just with a rather narrow bundle of profiles running along the DEC- EtOH edge and the EtOH-MeOH edge and cannot be reached from the central part of the space where DMC-EtOH charge compositions are expected. Even with a finite number of rectifying stages, such a high EtOH concentration must be maintained in the vessel in order to keep the EtOH composition in the feasible region that makes this opportunity practically useless.Such a deficiency is usually overcome by an extractive section, according to the lessons learned from experience in extractive distillation. If a stable node of the extractive profiles is situated nearer the DEC-EtOH edge or the EtOH-MeOH edge than a candidate curve runs, then a wide bundle of extractive profiles started from the middle part of the reactive space will intersect that candidate curve. Thus, the targeted pure MeOH vertex can be approached from within a wide region of potential vessel compositions by a reactive extractive section and a reactive rectifying section above it.EtOH is a good candidate to play the role of the entrainer because, according to experience collected in earlier feasibility studies, the stable node of the extractive profiles is shifted from the lightest point of the system toward the entrainer vertex with increasing feed ratio. Thus, one expects the stable node to approach either the EtOH-MeOH edge or the DEC-EtOH edge with increasing the flow rate of the EtOH feed. Finally, EtOH would be a preferable entrainer because it is already a component of the reaction system, and; thus, no foreign component would be involved.How the stable node is shifted with the feed ratioF/V according to calculations (at total reflux) is shown in Figure 4. It is located at the DMC-MeOH azeotrope when no entrainer is fed, and is practically located in the EtOH-MeOH edge at and above F/V=0.3. The bold line in the figure represents an arbitrarily selected candidate curve. As is shown in Figure 5, all the extractive profiles started from potential charge locations run to the stable node and, thus, intersect the candidate line, that is, the candidate rectifying/stripping profile, at F/V=0.3. Thus, with a feed ratio higher than 0.3, the desired purity can be reached from the whole quadrangle, that is, from any initial composition of the charge including the stoichiometric composition.Figure 4
Shifting of the stable node of the extractive profiles at total reflux.Figure 5
Reactive extractive profiles map at total reflux and feed ratioF/V=0.3.It is also possible to withdraw DEC as bottom product in continuous process, but it is not possible to withdraw pure MeOH as distillate without entrainer feeding.


### 3.2. Finite Reflux Analysis

Due to the high dimensionality of the nonreactive system, the nonreactive column profiles are rather difficult to visualize and analyze. Therefore, only the fully reactive column configurations are analyzed here with finite reflux and finite reboil ratios.


#### 3.2.1. Effect of Reflux RatioR and Upper Feed Ratio Fd/V

Reactive rectifying profiles with finite reflux ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the MeOH vertex. At R≤4, the profile runs along the EtOH-MeOH edge but does not reach the EtOH vertex. At reflux ratios R≥5 the profiles are longer, turn toward the MEC-DEC edge, and run along it before reaching the accumulation point. R=5 can be considered as a minimum reflux ratio (Figure 6). The feasibility region without extractive section is, thus, a composition area delimited by the MEC-DEC, DEC-EtOH, and EtOH-MeOH edges and the rectifying profile computed with R=5. Producing pure MeOH as distillate without extractive section is in principle possible at R≥5 from a narrow interval of the middle-vessel composition.Figure 6
Rectifying profiles with finite reflux ratios.However, maintaining feasible region’s compositions in the vessel would require a large concentration and amount of EtOH, and this would give rise to a huge-volume vessel. Economy of the process is expected to improve by applying EtOH entrainer and an extractive section, according to the results of the total reflux analysis in Section2.1.Reactive extractive profiles with finite reflux ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the rectifying profile are shown in Figure 6 at R=5 and Fd/V=0.2.At feed ratioFd/V=0.2, all the extractive profiles run toward a point outside the quadrangle, and all the extractive profiles initiated from outside the feasible region of the rectifying profiles intersect the actual rectifying profile. When the vessel composition is situated toward the DMC-MeOH edge from the reactive rectifying profile including the stoichiometric charge composition, the products can be reached with a combination of an extractive section and a rectifying section. Thus, a wide feasible region, shown with shading in the figure, is found.


#### 3.2.2. Effect of Stripping RatioS and Lower Feed Ratio Fw/L

Reactive stripping profiles with finite reboil ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the DEC vertex. As the stripping ratio is decreased, the profiles become shorter. The stripping profile shown in Figure 7 is computed to the limit-purity specification; the profiles belonging to purer DEC run between the shown profiles and the DEC-EtOH edge and then between the final sections of the shown profile and the MeOH-EtOH edge. This is a rather narrow feasible region, which does not connect the potential stoichiometric vessel composition and the destination region.Figure 7
Reactive extractive profiles atS=10 and at Fw/L=0.3.Applying an extractive section with EtOH as entrainer may extend the feasible region even in this case. Reactive extractive profiles with finite reboil ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the stripping profile are shown in Figure 7 at S=10 and Fw/L=0.3.With feed ratioFw/L=0.3, the extractive profiles (drawn by dashed lines) run toward two different attractors somewhere outside the quadrangle and leave the reaction space through the DEC-EtOH edge near vertex DEC and through the MeOH-EtOH edge near vertex EtOH. Most of the extractive profiles intersect the stripping profile and connect the feasible region with a large interval of potential vessel compositions, including the stoichiometric one.


#### 3.2.3. Effect of Production RatioD/W to the Vessel Path

The feasible region in a middle-vessel column is the region of feasible compositions in the middle vessel. These compositions should be feasible both for the upper part and the lower part; thus, the feasible region is the intersection of the two feasible regions determined for them. Considering extractive sections both above and below the vessel, the feasible region atR=5,S=10,Fd/V=0.2, and Fw/L=0.3 is shown by shading in Figure 8.Figure 8
Feasible region of middle-vessel column with two extractive sections.The vessel path is the path of the evolving vessel composition during the process. The reactive still path is computed by integrating the equation shown in Table2 from an arbitrarily selected initial point x0(DEC,DMC,EtOH,MEC,MeOH)=(0.013,0.160,0.184,0.134,0.509). The vessel path runs in different directions according to varied product ratio D/W, as is shown in Figure 8. The vessel composition gets farthest from the DMC vertex while kept in the feasible region at about D/W=10.If the initial point contains DMC and EtOH only, then the vessel composition turns toward the MeOH-EtOH edge at middle product ratios and can be kept in the feasible region with very low product ratio only. Thus, it is easier to let it reach the border of the feasible region, parallel to the MeOH-EtOH edge, and then apply theD/W=10 value.


## 3.1. Total Reflux/Reboil Analysis

Running the system with no product removal is the limit case of any batch distillation process and is usually analyzed before finite reflux ratios and/or finite reboil ratios are studied because total reflux and/or total reboil involves the highest reachable product purity in conventional distillation systems. It is not always so with extractive and reactive distillation, but performing a total reflux/reboil analysis is worth because valuable information can be obtained this way on the behavior of the studied system.


### 3.1.1. Nonreactive Sections

Nonreactive rectifying profiles and stripping profiles run along nonreactive residue curves in case of total reflux or total reboil, respectively. Thus, residue curves are analyzed.The studied system consists of five components. Thus, visualization of the residue curves is difficult. With the knowledge of the boiling points of the azeotropes and the pure components, the singular points and, thus, the shape of the residue curve map of the nonreactive system can be determined. In order to make the visualization easier, the system is projected to a plane (Figure2). According to this projected network, six-saddle points, one stable-node (SN), and one unstable node (UN) are found.Figure 2
The nonreactive system unfolded in 2D.The desired products are DEC as bottom product  and MeOH as distillate. Since the DEC vertex is the stable node of the system, production of pure DEC in the bottom is feasible both in batch and continuous processes. Just the contrary, pure MeOH cannot be produced at the top of a nonreactive rectifying section because the MeOH vertex is a saddle point (S) of the system. The top product would be the unstable node, that is, the DMC-MeOH azeotrope.Such a qualitative analysis of the nonreactive extractive profiles is rather difficult because their behavior cannot be visualized in a similar and clear way. Accumulation points (stable nodes) of the extractive profile bundles depend on process parameters like reflux ratioR or reboil ratio S and feed ratio F/V or F/L.


### 3.1.2. Reactive Sections

Due to the instantaneous reaction equilibrium, dimensionality of the reactive material system is reduced to a 2D hypersurface (of the equilibrium points) in the originally 4D linear space of the mole fractions. (There are five mole fractions, their sum is unity, and; thus, only four independent mole fractions are considered.) This curved hypersurface is referred to as “the reactive surface”. Due to applying the proposed transformation, this reactive surface, represented with reactive residue curves, is projected to a 2D plane as is shown in Figure3.Figure 3
Residue curves in the transformed reactive space.The reactive space contains the same unstable node (the locus of the DMC-MeOH azeotrope) and the same stable node (DEC vertex) as in the case of the nonreactive system. Thus, without extractive section, the DMC-MeOH azeotrope and pure DEC are withdrawn as distillate and bottom products, respectively, with infinite number of theoretical stages.There is a narrow bundle of reactive residue curves running so near the MeOH vertex that they intersect the composition region of assigned purity. Thus, the assigned purity can be reached with a finite number of rectifying stages. The feasible region is the set of compositions covered by that part of this bundle which connect vertex DEC with the targeted region around vertex MeOH.Unfortunately, the targeted MeOH purity can be reached just with a rather narrow bundle of profiles running along the DEC- EtOH edge and the EtOH-MeOH edge and cannot be reached from the central part of the space where DMC-EtOH charge compositions are expected. Even with a finite number of rectifying stages, such a high EtOH concentration must be maintained in the vessel in order to keep the EtOH composition in the feasible region that makes this opportunity practically useless.Such a deficiency is usually overcome by an extractive section, according to the lessons learned from experience in extractive distillation. If a stable node of the extractive profiles is situated nearer the DEC-EtOH edge or the EtOH-MeOH edge than a candidate curve runs, then a wide bundle of extractive profiles started from the middle part of the reactive space will intersect that candidate curve. Thus, the targeted pure MeOH vertex can be approached from within a wide region of potential vessel compositions by a reactive extractive section and a reactive rectifying section above it.EtOH is a good candidate to play the role of the entrainer because, according to experience collected in earlier feasibility studies, the stable node of the extractive profiles is shifted from the lightest point of the system toward the entrainer vertex with increasing feed ratio. Thus, one expects the stable node to approach either the EtOH-MeOH edge or the DEC-EtOH edge with increasing the flow rate of the EtOH feed. Finally, EtOH would be a preferable entrainer because it is already a component of the reaction system, and; thus, no foreign component would be involved.How the stable node is shifted with the feed ratioF/V according to calculations (at total reflux) is shown in Figure 4. It is located at the DMC-MeOH azeotrope when no entrainer is fed, and is practically located in the EtOH-MeOH edge at and above F/V=0.3. The bold line in the figure represents an arbitrarily selected candidate curve. As is shown in Figure 5, all the extractive profiles started from potential charge locations run to the stable node and, thus, intersect the candidate line, that is, the candidate rectifying/stripping profile, at F/V=0.3. Thus, with a feed ratio higher than 0.3, the desired purity can be reached from the whole quadrangle, that is, from any initial composition of the charge including the stoichiometric composition.Figure 4
Shifting of the stable node of the extractive profiles at total reflux.Figure 5
Reactive extractive profiles map at total reflux and feed ratioF/V=0.3.It is also possible to withdraw DEC as bottom product in continuous process, but it is not possible to withdraw pure MeOH as distillate without entrainer feeding.


## 3.1.1. Nonreactive Sections

Nonreactive rectifying profiles and stripping profiles run along nonreactive residue curves in case of total reflux or total reboil, respectively. Thus, residue curves are analyzed.The studied system consists of five components. Thus, visualization of the residue curves is difficult. With the knowledge of the boiling points of the azeotropes and the pure components, the singular points and, thus, the shape of the residue curve map of the nonreactive system can be determined. In order to make the visualization easier, the system is projected to a plane (Figure2). According to this projected network, six-saddle points, one stable-node (SN), and one unstable node (UN) are found.Figure 2
The nonreactive system unfolded in 2D.The desired products are DEC as bottom product  and MeOH as distillate. Since the DEC vertex is the stable node of the system, production of pure DEC in the bottom is feasible both in batch and continuous processes. Just the contrary, pure MeOH cannot be produced at the top of a nonreactive rectifying section because the MeOH vertex is a saddle point (S) of the system. The top product would be the unstable node, that is, the DMC-MeOH azeotrope.Such a qualitative analysis of the nonreactive extractive profiles is rather difficult because their behavior cannot be visualized in a similar and clear way. Accumulation points (stable nodes) of the extractive profile bundles depend on process parameters like reflux ratioR or reboil ratio S and feed ratio F/V or F/L.


## 3.1.2. Reactive Sections

Due to the instantaneous reaction equilibrium, dimensionality of the reactive material system is reduced to a 2D hypersurface (of the equilibrium points) in the originally 4D linear space of the mole fractions. (There are five mole fractions, their sum is unity, and; thus, only four independent mole fractions are considered.) This curved hypersurface is referred to as “the reactive surface”. Due to applying the proposed transformation, this reactive surface, represented with reactive residue curves, is projected to a 2D plane as is shown in Figure3.Figure 3
Residue curves in the transformed reactive space.The reactive space contains the same unstable node (the locus of the DMC-MeOH azeotrope) and the same stable node (DEC vertex) as in the case of the nonreactive system. Thus, without extractive section, the DMC-MeOH azeotrope and pure DEC are withdrawn as distillate and bottom products, respectively, with infinite number of theoretical stages.There is a narrow bundle of reactive residue curves running so near the MeOH vertex that they intersect the composition region of assigned purity. Thus, the assigned purity can be reached with a finite number of rectifying stages. The feasible region is the set of compositions covered by that part of this bundle which connect vertex DEC with the targeted region around vertex MeOH.Unfortunately, the targeted MeOH purity can be reached just with a rather narrow bundle of profiles running along the DEC- EtOH edge and the EtOH-MeOH edge and cannot be reached from the central part of the space where DMC-EtOH charge compositions are expected. Even with a finite number of rectifying stages, such a high EtOH concentration must be maintained in the vessel in order to keep the EtOH composition in the feasible region that makes this opportunity practically useless.Such a deficiency is usually overcome by an extractive section, according to the lessons learned from experience in extractive distillation. If a stable node of the extractive profiles is situated nearer the DEC-EtOH edge or the EtOH-MeOH edge than a candidate curve runs, then a wide bundle of extractive profiles started from the middle part of the reactive space will intersect that candidate curve. Thus, the targeted pure MeOH vertex can be approached from within a wide region of potential vessel compositions by a reactive extractive section and a reactive rectifying section above it.EtOH is a good candidate to play the role of the entrainer because, according to experience collected in earlier feasibility studies, the stable node of the extractive profiles is shifted from the lightest point of the system toward the entrainer vertex with increasing feed ratio. Thus, one expects the stable node to approach either the EtOH-MeOH edge or the DEC-EtOH edge with increasing the flow rate of the EtOH feed. Finally, EtOH would be a preferable entrainer because it is already a component of the reaction system, and; thus, no foreign component would be involved.How the stable node is shifted with the feed ratioF/V according to calculations (at total reflux) is shown in Figure 4. It is located at the DMC-MeOH azeotrope when no entrainer is fed, and is practically located in the EtOH-MeOH edge at and above F/V=0.3. The bold line in the figure represents an arbitrarily selected candidate curve. As is shown in Figure 5, all the extractive profiles started from potential charge locations run to the stable node and, thus, intersect the candidate line, that is, the candidate rectifying/stripping profile, at F/V=0.3. Thus, with a feed ratio higher than 0.3, the desired purity can be reached from the whole quadrangle, that is, from any initial composition of the charge including the stoichiometric composition.Figure 4
Shifting of the stable node of the extractive profiles at total reflux.Figure 5
Reactive extractive profiles map at total reflux and feed ratioF/V=0.3.It is also possible to withdraw DEC as bottom product in continuous process, but it is not possible to withdraw pure MeOH as distillate without entrainer feeding.


## 3.2. Finite Reflux Analysis

Due to the high dimensionality of the nonreactive system, the nonreactive column profiles are rather difficult to visualize and analyze. Therefore, only the fully reactive column configurations are analyzed here with finite reflux and finite reboil ratios.


### 3.2.1. Effect of Reflux RatioR and Upper Feed Ratio Fd/V

Reactive rectifying profiles with finite reflux ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the MeOH vertex. At R≤4, the profile runs along the EtOH-MeOH edge but does not reach the EtOH vertex. At reflux ratios R≥5 the profiles are longer, turn toward the MEC-DEC edge, and run along it before reaching the accumulation point. R=5 can be considered as a minimum reflux ratio (Figure 6). The feasibility region without extractive section is, thus, a composition area delimited by the MEC-DEC, DEC-EtOH, and EtOH-MeOH edges and the rectifying profile computed with R=5. Producing pure MeOH as distillate without extractive section is in principle possible at R≥5 from a narrow interval of the middle-vessel composition.Figure 6
Rectifying profiles with finite reflux ratios.However, maintaining feasible region’s compositions in the vessel would require a large concentration and amount of EtOH, and this would give rise to a huge-volume vessel. Economy of the process is expected to improve by applying EtOH entrainer and an extractive section, according to the results of the total reflux analysis in Section2.1.Reactive extractive profiles with finite reflux ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the rectifying profile are shown in Figure 6 at R=5 and Fd/V=0.2.At feed ratioFd/V=0.2, all the extractive profiles run toward a point outside the quadrangle, and all the extractive profiles initiated from outside the feasible region of the rectifying profiles intersect the actual rectifying profile. When the vessel composition is situated toward the DMC-MeOH edge from the reactive rectifying profile including the stoichiometric charge composition, the products can be reached with a combination of an extractive section and a rectifying section. Thus, a wide feasible region, shown with shading in the figure, is found.


### 3.2.2. Effect of Stripping RatioS and Lower Feed Ratio Fw/L

Reactive stripping profiles with finite reboil ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the DEC vertex. As the stripping ratio is decreased, the profiles become shorter. The stripping profile shown in Figure 7 is computed to the limit-purity specification; the profiles belonging to purer DEC run between the shown profiles and the DEC-EtOH edge and then between the final sections of the shown profile and the MeOH-EtOH edge. This is a rather narrow feasible region, which does not connect the potential stoichiometric vessel composition and the destination region.Figure 7
Reactive extractive profiles atS=10 and at Fw/L=0.3.Applying an extractive section with EtOH as entrainer may extend the feasible region even in this case. Reactive extractive profiles with finite reboil ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the stripping profile are shown in Figure 7 at S=10 and Fw/L=0.3.With feed ratioFw/L=0.3, the extractive profiles (drawn by dashed lines) run toward two different attractors somewhere outside the quadrangle and leave the reaction space through the DEC-EtOH edge near vertex DEC and through the MeOH-EtOH edge near vertex EtOH. Most of the extractive profiles intersect the stripping profile and connect the feasible region with a large interval of potential vessel compositions, including the stoichiometric one.


### 3.2.3. Effect of Production RatioD/W to the Vessel Path

The feasible region in a middle-vessel column is the region of feasible compositions in the middle vessel. These compositions should be feasible both for the upper part and the lower part; thus, the feasible region is the intersection of the two feasible regions determined for them. Considering extractive sections both above and below the vessel, the feasible region atR=5,S=10,Fd/V=0.2, and Fw/L=0.3 is shown by shading in Figure 8.Figure 8
Feasible region of middle-vessel column with two extractive sections.The vessel path is the path of the evolving vessel composition during the process. The reactive still path is computed by integrating the equation shown in Table2 from an arbitrarily selected initial point x0(DEC,DMC,EtOH,MEC,MeOH)=(0.013,0.160,0.184,0.134,0.509). The vessel path runs in different directions according to varied product ratio D/W, as is shown in Figure 8. The vessel composition gets farthest from the DMC vertex while kept in the feasible region at about D/W=10.If the initial point contains DMC and EtOH only, then the vessel composition turns toward the MeOH-EtOH edge at middle product ratios and can be kept in the feasible region with very low product ratio only. Thus, it is easier to let it reach the border of the feasible region, parallel to the MeOH-EtOH edge, and then apply theD/W=10 value.


## 3.2.1. Effect of Reflux RatioR and Upper Feed Ratio Fd/V

Reactive rectifying profiles with finite reflux ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the MeOH vertex. At R≤4, the profile runs along the EtOH-MeOH edge but does not reach the EtOH vertex. At reflux ratios R≥5 the profiles are longer, turn toward the MEC-DEC edge, and run along it before reaching the accumulation point. R=5 can be considered as a minimum reflux ratio (Figure 6). The feasibility region without extractive section is, thus, a composition area delimited by the MEC-DEC, DEC-EtOH, and EtOH-MeOH edges and the rectifying profile computed with R=5. Producing pure MeOH as distillate without extractive section is in principle possible at R≥5 from a narrow interval of the middle-vessel composition.Figure 6
Rectifying profiles with finite reflux ratios.However, maintaining feasible region’s compositions in the vessel would require a large concentration and amount of EtOH, and this would give rise to a huge-volume vessel. Economy of the process is expected to improve by applying EtOH entrainer and an extractive section, according to the results of the total reflux analysis in Section2.1.Reactive extractive profiles with finite reflux ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the rectifying profile are shown in Figure 6 at R=5 and Fd/V=0.2.At feed ratioFd/V=0.2, all the extractive profiles run toward a point outside the quadrangle, and all the extractive profiles initiated from outside the feasible region of the rectifying profiles intersect the actual rectifying profile. When the vessel composition is situated toward the DMC-MeOH edge from the reactive rectifying profile including the stoichiometric charge composition, the products can be reached with a combination of an extractive section and a rectifying section. Thus, a wide feasible region, shown with shading in the figure, is found.


## 3.2.2. Effect of Stripping RatioS and Lower Feed Ratio Fw/L

Reactive stripping profiles with finite reboil ratio are computed with the equations shown in Table1. Numerical integration of the differential equation is started from an initial point inside the targeted purity region near the DEC vertex. As the stripping ratio is decreased, the profiles become shorter. The stripping profile shown in Figure 7 is computed to the limit-purity specification; the profiles belonging to purer DEC run between the shown profiles and the DEC-EtOH edge and then between the final sections of the shown profile and the MeOH-EtOH edge. This is a rather narrow feasible region, which does not connect the potential stoichiometric vessel composition and the destination region.Figure 7
Reactive extractive profiles atS=10 and at Fw/L=0.3.Applying an extractive section with EtOH as entrainer may extend the feasible region even in this case. Reactive extractive profiles with finite reboil ratio and feed ratio are computed with the equations shown in Table1. The initial compositions are selected near an edge of the reaction quadrangle. A few computed extractive profiles together with the stripping profile are shown in Figure 7 at S=10 and Fw/L=0.3.With feed ratioFw/L=0.3, the extractive profiles (drawn by dashed lines) run toward two different attractors somewhere outside the quadrangle and leave the reaction space through the DEC-EtOH edge near vertex DEC and through the MeOH-EtOH edge near vertex EtOH. Most of the extractive profiles intersect the stripping profile and connect the feasible region with a large interval of potential vessel compositions, including the stoichiometric one.


## 3.2.3. Effect of Production RatioD/W to the Vessel Path

The feasible region in a middle-vessel column is the region of feasible compositions in the middle vessel. These compositions should be feasible both for the upper part and the lower part; thus, the feasible region is the intersection of the two feasible regions determined for them. Considering extractive sections both above and below the vessel, the feasible region atR=5,S=10,Fd/V=0.2, and Fw/L=0.3 is shown by shading in Figure 8.Figure 8
Feasible region of middle-vessel column with two extractive sections.The vessel path is the path of the evolving vessel composition during the process. The reactive still path is computed by integrating the equation shown in Table2 from an arbitrarily selected initial point x0(DEC,DMC,EtOH,MEC,MeOH)=(0.013,0.160,0.184,0.134,0.509). The vessel path runs in different directions according to varied product ratio D/W, as is shown in Figure 8. The vessel composition gets farthest from the DMC vertex while kept in the feasible region at about D/W=10.If the initial point contains DMC and EtOH only, then the vessel composition turns toward the MeOH-EtOH edge at middle product ratios and can be kept in the feasible region with very low product ratio only. Thus, it is easier to let it reach the border of the feasible region, parallel to the MeOH-EtOH edge, and then apply theD/W=10 value.


## 4. Validation by Rigorous Simulation

The results of the feasibility study are validated by rigorous simulations performed using ProSimBatch, a dynamic simulator for batch distillation and reactive processes. This software product works by numerically integrating differential equations, and; thus, it can also compute the start-up phase of the process. Unfortunately, the method of the start-up procedure is fixed in the software, and is not a selectable option.This simulator allows performing simulation of batch-reactive distillation processes in a rectifier, in a stripper, or in a middle-vessel column. In the case of a middle-vessel column, an additional withdrawal has to be defined at the bottom of the column. The middle vessel is defined as a tray with large holdup; its leaving liquid flow rate and its heat duty can be specified as well.Instantaneously equilibrium-limited reaction cannot be specified during the simulations. Instead, four kinetically controlled reactions are specified. The reaction rate coefficients are set large (to achieve rather fast reactions) in such a way that their ratios are equal to the equilibrium constants of the equilibrium limited reactions:(12)K1=k+1k-1=2.33,K2=k+2k-2=0.26.
For simulation, stage numbers should also have to be specified. In order to obtain initial estimates for the needed number of stages, preliminary calculations are performed. For this aim, not the simulator but a simple stage-by-stage calculating device is used. The calculations are started from the vessel composition, performed along the supposed extractive profile, and,  when a composition similar to that obtained by the feasibility study is achieved, then the calculation is continued with the supposed rectifier or stripping profile until the specified purity of the product is reached. Considering the initial and final vessel compositions, the calculated minimum and maximum stage numbers are 30 and 56 in the rectifier section, 20 and 20 in the stripper section, 6 and 16 in the upper extractive section, and 5 and 10 in the lower extractive section.


### 4.1. Validation of the Reactive Profiles

Five column configurations are considered with different reflux and reboil ratios and with different feed ratios in order to compare the results of the simulated and the calculated profiles. The parameters applied in the simulations are collected in Table7.Table 7
Parameters of the simulations with ProSimBatch.
Simulation #12345R∞∞5∞5S∞∞∞1010Fd/V00.30.200.2Fw/L0000.30.3Nrect3030303030Nupperext1616161616Nstrip1010101010Nlowerext1010101010Qcond4880 W4880 W4880 W1035 W1035 WQreboil1020 W1020 W1020 W5180 W5180 WHold-up (tray)0.35 mol0.35 mol0.35 mol30 cm330 cm3Hold-up (MV)350 mol350 mol350 mol3000 cm33000 cm3Charge4000 mol4000 mol4000 mol5125 cm35125 cm3The simulated composition profiles along the column height are represented by thick lines in Figures9 to 13. Some profile curves calculated in the feasibility study are represented by thin lines in these figures in order to compare the simulation results with these curve bundles. Middle-vessel compositions are labeled by letters MV.Figure 9
Simulated and calculated residue curves.Two cases are simulated with total reflux and reboil, no feed, and two arbitrarily selected vessel compositions, as shown in Figure9. The dotted and the continuous thick lines represent the sections below and above the middle vessel, respectively. The residue curves calculated during the feasibility study (drawn by thin continuous lines) are also shown in the figure. The simulated reactive profiles are well approximated by the calculated curves.Two composition profiles simulated with total reflux and reboil, and with an average feed ratio (Fd/V=0.3) in the upper part of the middle vessel column, and with two arbitrarily selected reboiler compositions, are shown in Figure 10. The dotted, the dashed, and the continuous lines correspond, respectively, to the section below the middle vessel and the sections below and above the feed in the upper part of the middle vessel column. The (dashed) extractive profiles approach the stable node of the extractive profile curves’ bundle, as is also shown by two thin dotted lines calculated during the feasibility study. The thick continuous lines, representing the simulated reactive rectifying profiles, run almost together along the MeOH-EtOH edge and the EtOH-DMC edge (see also Section 3.1.2). According to this Figure, the results of the feasibility study provide a good estimation for the composition profile along the column.Figure 10
Simulated and calculated profiles at total reflux and reboil ratio, and atFd/V=0.3.A simulated column profile with total reboil, finite reflux ratio (R=5), and with an average feed ratio (Fd/V=0.2) is shown in Figure 11. The dotted, dashed, and the continuous thick lines correspond, respectively, to the stripping section below the middle vessel, and the sections below and above the feed in the upper part of the middle vessel column. The continuous and the dotted thin lines represent the reactive rectifying profile curve and a reactive upper extractive curve, respectively, as calculated in the feasibility study. The thin lines again provide good estimation to the simulated ones. The simulated rectifying and extractive profiles meet in a point nearer the EtOH vertex than anticipated because the simulated distillate composition is a little bit different from that specified in the feasibility calculations.Figure 11
Simulated and calculated profiles at total reboil ratio, atR=5 and at Fd/V=0.2.Two simulated column profiles with total reflux, finite reboil ratio (S=10) and with an average feed ratio (Fw/L=0.3) are shown in Figure 12. The dotted, the dashed, and the continuous thick lines, respectively, represent the rectifying section and the sections above and below the feed in the lower part of the middle vessel column. Two reactive stripping profile curves calculated in the feasibility study are shown by thin continuous lines, and two lower reactive extractive profile curves are also shown by thin dotted lines in this figure. The tendency of the calculated profile curves and the simulated profiles agree in this case, too.Figure 12
Simulated and calculated profiles at total reflux, atS=10, and at Fw/L=0.3.Figure 13
Simulated and   calculated profiles during the feasibility study atR=5, at Fd/V=0.2, at S=10 and at Fw/L=0.3.Two simulated column profiles with finite reflux ratio (R=5), finite reboil ratio (S=10), upper feed ratio Fd/V=0.2, and lower feed ratio (Fw/L=0.3) are shown in Figure 13. The simulated rectifying and stripping profiles and the simulated extractive profiles in the upper part and in the lower part are represented by thick continuous, dotted, and dashed lines, respectively. Profile curves calculated in the feasibility study are shown with thin lines. The calculated profile curves well approximate the simulated profiles in the lower part of the column but show deviation in the upper part because the simulated distillate composition is a little bit different from that specified in the feasibility calculations. However, the calculated profile curves give valid information on the qualitative shape of the simulated profiles even in this case.


## 4.1. Validation of the Reactive Profiles

Five column configurations are considered with different reflux and reboil ratios and with different feed ratios in order to compare the results of the simulated and the calculated profiles. The parameters applied in the simulations are collected in Table7.Table 7
Parameters of the simulations with ProSimBatch.
Simulation #12345R∞∞5∞5S∞∞∞1010Fd/V00.30.200.2Fw/L0000.30.3Nrect3030303030Nupperext1616161616Nstrip1010101010Nlowerext1010101010Qcond4880 W4880 W4880 W1035 W1035 WQreboil1020 W1020 W1020 W5180 W5180 WHold-up (tray)0.35 mol0.35 mol0.35 mol30 cm330 cm3Hold-up (MV)350 mol350 mol350 mol3000 cm33000 cm3Charge4000 mol4000 mol4000 mol5125 cm35125 cm3The simulated composition profiles along the column height are represented by thick lines in Figures9 to 13. Some profile curves calculated in the feasibility study are represented by thin lines in these figures in order to compare the simulation results with these curve bundles. Middle-vessel compositions are labeled by letters MV.Figure 9
Simulated and calculated residue curves.Two cases are simulated with total reflux and reboil, no feed, and two arbitrarily selected vessel compositions, as shown in Figure9. The dotted and the continuous thick lines represent the sections below and above the middle vessel, respectively. The residue curves calculated during the feasibility study (drawn by thin continuous lines) are also shown in the figure. The simulated reactive profiles are well approximated by the calculated curves.Two composition profiles simulated with total reflux and reboil, and with an average feed ratio (Fd/V=0.3) in the upper part of the middle vessel column, and with two arbitrarily selected reboiler compositions, are shown in Figure 10. The dotted, the dashed, and the continuous lines correspond, respectively, to the section below the middle vessel and the sections below and above the feed in the upper part of the middle vessel column. The (dashed) extractive profiles approach the stable node of the extractive profile curves’ bundle, as is also shown by two thin dotted lines calculated during the feasibility study. The thick continuous lines, representing the simulated reactive rectifying profiles, run almost together along the MeOH-EtOH edge and the EtOH-DMC edge (see also Section 3.1.2). According to this Figure, the results of the feasibility study provide a good estimation for the composition profile along the column.Figure 10
Simulated and calculated profiles at total reflux and reboil ratio, and atFd/V=0.3.A simulated column profile with total reboil, finite reflux ratio (R=5), and with an average feed ratio (Fd/V=0.2) is shown in Figure 11. The dotted, dashed, and the continuous thick lines correspond, respectively, to the stripping section below the middle vessel, and the sections below and above the feed in the upper part of the middle vessel column. The continuous and the dotted thin lines represent the reactive rectifying profile curve and a reactive upper extractive curve, respectively, as calculated in the feasibility study. The thin lines again provide good estimation to the simulated ones. The simulated rectifying and extractive profiles meet in a point nearer the EtOH vertex than anticipated because the simulated distillate composition is a little bit different from that specified in the feasibility calculations.Figure 11
Simulated and calculated profiles at total reboil ratio, atR=5 and at Fd/V=0.2.Two simulated column profiles with total reflux, finite reboil ratio (S=10) and with an average feed ratio (Fw/L=0.3) are shown in Figure 12. The dotted, the dashed, and the continuous thick lines, respectively, represent the rectifying section and the sections above and below the feed in the lower part of the middle vessel column. Two reactive stripping profile curves calculated in the feasibility study are shown by thin continuous lines, and two lower reactive extractive profile curves are also shown by thin dotted lines in this figure. The tendency of the calculated profile curves and the simulated profiles agree in this case, too.Figure 12
Simulated and calculated profiles at total reflux, atS=10, and at Fw/L=0.3.Figure 13
Simulated and   calculated profiles during the feasibility study atR=5, at Fd/V=0.2, at S=10 and at Fw/L=0.3.Two simulated column profiles with finite reflux ratio (R=5), finite reboil ratio (S=10), upper feed ratio Fd/V=0.2, and lower feed ratio (Fw/L=0.3) are shown in Figure 13. The simulated rectifying and stripping profiles and the simulated extractive profiles in the upper part and in the lower part are represented by thick continuous, dotted, and dashed lines, respectively. Profile curves calculated in the feasibility study are shown with thin lines. The calculated profile curves well approximate the simulated profiles in the lower part of the column but show deviation in the upper part because the simulated distillate composition is a little bit different from that specified in the feasibility calculations. However, the calculated profile curves give valid information on the qualitative shape of the simulated profiles even in this case.


## 5. Conclusion

The feasibility method of batch reactive distillation published in Stéger et al. [2] is extended for dealing with consecutive reactions. General model equations are derived for multireaction system with any number of chemical equilibrium-limited reactions. Since the presented general methodology is of graphical nature, performing the feasibility study of systems with more than 3 degrees of freedom might be difficult in practice. Thus, further investigation is needed for such, high complexity, systems.The present methodology, applying some modifications, might be used to continuous processes, as well.The presented methodology is demonstrated with the detailed study of the transesterification of dimethyl carbonate. Since the system contains 5 components, and the degree of freedom of any nonreactive section is 4; only the total reflux analysis is performed for the nonreactive sections. The feasibility analysis of the reactive sections is successfully performed.  The proposed new transformation is used to present the graphical methodology applicable to multireaction and multicomponent system.According to the presented feasibility study, at least one of the fully reactive middle-vessel column configurations is feasible. Extractive stages are needed both in the upper and lower parts of the column, with the following limiting parameters:R≥5,Fd/V≥0.2,S≥10,Fw/L≥0.3. However, only the configurations proposed by Stéger et al. [2] have been investigated in the present work, and there may be other feasible configurations depending on the place and number of the applied feeds and reactive sections. Any of them can be studied with the presented method.Effect of production ratio (D/W) to the vessel path and feasible region of reactive vessel is also examined.The results of the feasibility study are successfully validated by rigorous simulations.



---

*Source: 231828-2011-10-26.xml*