Distillation binary, model

Fast and satisfactory mass transfer calculations are necessary since we may have to repeat such calculations many times for a rate-based distillation column model or two-phase flow with mass transfer between the phases in the design and simulation process. The generalized matrix method may be used for multicomponent mass transfer calculations. The generalized matrix method utilizes the Maxwell-Stefan model with the linearized film model for diffusion flux, assuming a constant diffusion coefficient matrix and total concentration in the diffusion region. In an isotropic medium, Fick s law may describe the multicomponent molecular mass transfer at a specified temperature and pressure, assuming independent diffusion of the species in a fluid mixture. Such independent diffusion, however, is only an approximation in the following cases (i) diffusion of a dilute component in a solvent, (ii) diffusion of various components with identical diffusion properties, and (iii) diffusion in a binary mixture. 

Earlier chapters use simplified and binary models to analyze in a very informative manner some fundamentals such as the effect of reflux ratio and feed tray location, and to delineate the differences between absorption/stripping and distillation. Following chapters concentrate on specific areas such as complex distillation, with detailed analyses of various features such as pumparounds and side-strippers, and when they should be used. Also discussed are azeotropic, extractive, and three-phase distillation operations, multi-component liquid-liquid and supercritical extraction, and reactive multistage separation. The applications are clearly explained with many practical examples. 

In Examples 5.1 and 5.2 the output variables coincide with the state variables of the two processes. Consequently, in order to develop the input-output model we need only solve the differential equations of the mass and energy balances. This is not always true. Take as an example the binary distillation column model (Example 4.13 and Figure 4.10). For this system we have ... 

Do the same as in Problem III. 1 for the equations describing the dynamic and steady-state behavior of the binary distillation column modeled in Example 4.13. 

In this study, the phase equilibrium in the binary mixtures that are expected to be found in the flash distillation was modeled with the Predictive Soave-Redlich-Kwong (PSRK) equation of state [4], using modified molecular parameters r and q. Five binary ethanol + congener mixtures were considered for new yield values for parameters r and q. The congeners considered were acetic acid, acetaldehyde, furfural, methanol, and 1-pentanol. Subsequently, the model was validated with the water + ethanol binary system, and the 1 -pentanol + ethanol + water, 1-propanol + ethanol + water, and furfural + ethanol + water ternary systems. 

Using a steady-state film model, obtain an expression for the mass transfer rate across a laminar film of thickness /. in the vapour phase for the more volatile component in a binary distillation process ... 

The conlinuous binary distillation column of Fig. 3.60 follows the same general representation as that used previously in Fig. 3.58. The modelling approach again follows closely that of Luyben (1973, 1990). 

Figure 3.60. Model representation of a continuous binary distillation column. PC is the cooling water controller, LC the reflux controller.

The variation of efficiencies is due to interaction phenomena caused by the simultaneous diffusional transport of several components. From a fundamental point of view one should therefore take these interaction phenomena explicitly into account in the description of the elementary processes (i.e. mass and heat transfer with chemical reaction). In literature this approach has been used within the non-equilibrium stage model (Sivasubramanian and Boston, 1990). Sawistowski (1983) and Sawistowski and Pilavakis (1979) have developed a model describing reactive distillation in a packed column. Their model incorporates a simple representation of the prevailing mass and heat transfer processes supplemented with a rate equation for chemical reaction, allowing chemical enhancement of mass transfer. They assumed elementary reaction kinetics, equal binary diffusion coefficients and equal molar latent heat of evaporation for each component. 

The nonlinear nature of these mixed-integer optimization problems may arise from (i) nonlinear relations in the integer domain exclusively (e.g., products of binary variables in the quadratic assignment model), (ii) nonlinear relations in the continuous domain only (e.g., complex nonlinear input-output model in a distillation column or reactor unit), (iii) nonlinear relations in the joint integer-continuous domain (e.g., products of continuous and binary variables in the schedul-ing/planning of batch processes, and retrofit of heat recovery systems). In this chapter, we will focus on nonlinearities due to relations (ii) and (iii). An excellent book that studies mixed-integer linear optimization, and nonlinear integer relationships in combinatorial optimization is the one by Nemhauser and Wolsey (1988). 

For the synthesis of heterogeneous batch distillation the liquid-liquid envelope at the decanter temperature is considered in addition to the residue curve map. Therefore, the binary interaction parameters used in predicting liquid-liquid equilibrium are estimated from binary heterogeneous azeotrope or liquid-liquid equilibrium data [8,10], Table 3 shows the calculated purity of original components in each phase split at 25 °C for all heterogeneous azeotropes reported in Table 1. The thermodynamic models and binary coefficients used in the calculation of the liquid-liquid-vapour equilibrium, liquid-liquid equilibrium at 25 °C and the separatrices are reported in Table 2. 

RGA Example To illustrate use of the RGA method, consider the following steady-state version of a transfer function model for a pilot-scale, methanol-water distillation column (Wood and Berry, "Terminal Composition Control of a Binary Distillation Column, Chem. Eng. Sci., 28 1707, 1973) Kn = 12.8, Ku = -18.9, K21 = 6.6, and K22 = —19.4. It follows that X = 2 and... 

Converse and Huber (1965), Robinson (1970), Mayur and Jackson (1971), Luyben (1988) and Mujtaba (1997) used this model for simulation and optimisation of conventional batch distillation. Domenech and Enjalbert (1981) used similar model in their simulation study with the exception that they used temperature dependent phase equilibria instead of constant relative volatility. Christiansen et al. (1995) used this model (excluding column holdup) to study parametric sensitivity of ideal binary columns. 

A liquid binary mixture with B0 = 10 kmol (Hc) and xB0 = <0.6, 0.4> (xj) molefraction is subject to inverted batch distillation shown in Figure 4.12. The relative volatility of the mixture over the operating temperature range is assumed constant with a value of (a-) 2. The number of plates is, N= 10. The vapour boilup rate is, V = 10.0 kmol/hr. The total plate holdup is 0.3 kmol and the reboiler holdup is 0.1 kmol. The total batch time of operation is 4 hr with two time intervals. The first interval is of duration 1 hr and the column is operated with a reboil ratio of 0.8. The second interval is of duration 3 hrs when the column is operated with a reboil ratio of 0.9. The column operation is simulated with the type III model (section 4.3.2.1). 

Mujtaba and Hussain (1998) implemented the general optimisation framework based on the hybrid scheme for a binary batch distillation process. It was shown that the optimal control policy using a detailed process model was very close to that obtained using the hybrid model. 

The proposed modeling approach has been validated for distillation of non-reactive mixtures. For this purpose, the use is made of the total reflux distillation data for the binary mixture chlorobenzene/ethylbenzene (CB/EB) and ternary mixture methanol/acetonitrile/water (MEOH/ACN/WATER) obtained by Pelkonen (1997) as well as for the ternary mixture methanol/ethanol/water (MeOH/EtOH/WATER) measured by Mori et al. (2006). The experiments of Pelkonen (1997) were carried out in a column of 100 mm diameter, equipped with Montz-Pak A3-500 structured packing. The measured concentrations, temperature and flow rates at the condenser outlet are used as input values for simulations. 

Up to now, the proposed model has been validated using the total reflux distillation data in the column equipped with the Montz-Pak A3-500 and Montz-Pak B1-250 structured packings. A very good agreement between the calculated and experimental data for binary and ternary mixtures is found. 

The equilibrium distillation behavior of the model fuels is adequately covered in the fuel oil discussion. However, the case for the rapid droplet vaporization, which was not clearly seen for fuel oils, is more amenable to analysis for a binary system. The surface gradients are given by the following relationship,... 

Although the widely used equilibrium-stage models for distillation, described above, have proved to be quite adequate for binary and closeboiling, ideal and near-ideal multicomponent vapor-liquid mixtures. 

Strictly speaking, Eqs. (13-69) and (13-70) are valid only for describing mass transfer in binary systems under conditions where the rates of mass transfer are low. Most industrial distillation and absorption processes, however, involve more than two different chemical species. The most fundamentally sound way to model mass transfer in multi-component systems is to use the Maxwell-Stefan (MS) approach (Taylor and Krishna, op. cit.). 

---