Mixture interaction parameters

Since both ethylene and ethane have reduced temperatures nearly equal to unity at the extraction conditions of 20 C, (T =. 98) and ethylene (T = 1.04), their respective solvent capacities for butene should be about the same. This is the case as is reflected in the same values for the selectivity against butene for all pure solvent gases. One can conclude that the primary effect of the non-polar solvent is to increase the capacity of the "vapor" phase for the extracted solute near the critical. The influence of the second solvent provides only the option of modifying the physical parameters namely, pressure and temperature, under which the optimal extraction is to be conducted. The evidence for this is the effect of the ammonia on the selectivity as calculated by the EOS in Table V. The higher values for the selectivities in the ethylene mixtures are pronounced. It can be concluded that the solvent mixture interaction parameters must dominate the solubility of butene in the vapor phase. 

Lichtenthaler et al. (55) determined interaction parameters for 22 solutes in poly(dimethyl siloxane) to test several expressions of the combinatorial entropy of mixing [Eq. (7)]. The magnitude of the interaction parameter is indeed directly dependent on the evaluation of the combinatorial contribution. The combinatorial contribution was computed following both the Flory-Huggins approximation and the multiple-connected-site model recently developed by Lichtenthaler, Abrams and Prausnitz (56). This model, which retains the Flory-Huggins term, also corrects for the bulkiness of the components of the mixture. Interaction parameters were computed through both approximations, showing the sensitivity of the results to the model chosen. 

The data base contains provisions for a simple augmentation by up to eight additional compounds or substitution of other compounds for those included. Binary interaction parameters necessary for calculation of fugacities in liquid mixtures are presently available for 180 pairs. 

INTERACTION PARAMETERS FOR LIQUIO-PHASE MIXTURES COMPONENT NAMES... 

While the phase rule requires tliree components for an unsymmetrical tricritical point, theory can reduce this requirement to two components with a continuous variation of the interaction parameters. Lindli et al (1984) calculated a phase diagram from the van der Waals equation for binary mixtures and found (in accord with figure A2.5.13 that a tricritical point occurred at sufficiently large values of the parameter (a measure of the difference between the two components). 

The quantity x is called the Flory-Huggins interaction parameter It is zero for athermal mixtures, positive for endothermic mixing, and negative for exothermic mixing. These differences in sign originate from Eq. (8.39) and reaction (8.A). 

When i = J, all equations reduce to the appropriate values for a pure species. When i j, these equations define a set of interaction parameters having no physical significance. For a mixture, values of By and dBjj/dT from Eqs. (4-212) and (4-213) are substituted into Eqs. (4-183) and (4-185) to provide values of the mixture second virial coefficient B and its temperature derivative. Values of and for the mixture are then given by Eqs. (4-193) and (4-194), and values of In i for the component fugacity coefficients are given by Eq. (4-196). 

Binary interaction parameters are determined for each pq pair p q) from experimental data. Note that = k and k = k = 0. Since the quantity on the left-hand side of Eq. (4-305) represents the second virial coefficient as predicted by Eq. (4-231), the basis for Eq. (4-305) lies in Eq. (4-183), which expresses the quadratic dependence of the mixture second virial coefficient on mole fraction. 

The lattice gas has been used as a model for a variety of physical and chemical systems. Its application to simple mixtures is routinely treated in textbooks on statistical mechanics, so it is natural to use it as a starting point for the modeling of liquid-liquid interfaces. In the simplest case the system contains two kinds of solvent particles that occupy positions on a lattice, and with an appropriate choice of the interaction parameters it separates into two phases. This simple version is mainly of didactical value [1], since molecular dynamics allows the study of much more realistic models of the interface between two pure liquids [2,3]. However, even with the fastest computers available today, molecular dynamics is limited to comparatively small ensembles, too small to contain more than a few ions, so that the space-charge regions cannot be included. In contrast, Monte Carlo simulations for the lattice gas can be performed with 10 to 10 particles, so that modeling of the space charge poses no problem. In addition, analytical methods such as the quasichemical approximation allow the treatment of infinite ensembles. 

Thermodynamic models are widely used for the calculation of equilibrium and thermophysical properties of fluid mixtures. Two types of such models will be examined cubic equations of state and activity coefficient models. In this chapter cubic equations of state models are used. Volumetric equations of state (EoS) are employed for the calculation of fluid phase equilibrium and thermophysical properties required in the design of processes involving non-ideal fluid mixtures in the oil and gas and chemical industries. It is well known that the introduction of empirical parameters in equation of state mixing rules enhances the ability of a given EoS as a tool for process design although the number of interaction parameters should be as small as possible. In general, the phase equilibrium calculations with an EoS are very sensitive to the values of the binary interaction parameters. 

It is assumed that there are available NCP experimental binary critical point data. These data include values of the pressure, Pc, the temperature, Tc, and the mole fraction, xc, of one of the components at each of the critical points for the binary mixture. The vector k of interaction parameters is determined by fitting the EoS to the critical data. In explicit formulations the interaction parameters are obtained by the minimization of the following least squares objective function ... 

These models are semiempirical and are based on the concept that intermolecular forces will cause nonrandom arrangement of molecules in the mixture. The models account for the arrangement of molecules of different sizes and the preferred orientation of molecules. In each case, the models are fitted to experimental binary vapor-liquid equilibrium data. This gives binary interaction parameters that can be used to predict multicomponent vapor-liquid equilibrium. In the case of the UNIQUAC equation, if experimentally determined vapor-liquid equilibrium data are not available, the Universal Quasi-chemical Functional Group Activity Coefficients (UNIFAC) method can be used to estimate UNIQUAC parameters from the molecular structures of the components in the mixture3. 

A distillation calculation is to be performed on a multicomponent mixture. The vapor-liquid equilibrium for this mixture is likely to exhibit significant departures from ideality, but a complete set of binary interaction parameters is not available. What factors would you consider in assessing whether the missing interaction parameters are likely to have an important effect on the calculations ... 

Example 11.1 Each component for the mixture of alkanes in Table 11.2 is to be separated into relatively pure products. Table 11.2 shows normal boiling points and relative volatilities to indicate the order of volatility and the relative difficulty of the separations. The relative volatilities have been calculated on the basis of the feed composition to the sequence, assuming a pressure of 6 barg using the Peng-Robinson Equation of State with interaction parameters set to zero (see Chapter 4). Different pressures can, in practice, be used for different columns in the sequence and if a single set of relative volatilities is to be used, the pressure at which the relative volatilities are calculated needs, as much as possible, to be chosen to represent the overall system. 

Example 11.2 Using the Underwood Equations, determine the best distillation sequence, in terms of overall vapor load, to separate the mixture of alkanes in Table 11.2 into relatively pure products. The recoveries are to be assumed to be 100%. Assume the ratio of actual to minimum reflux ratio to be 1.1 and all columns are fed with a saturated liquid. Neglect pressure drop across each column. Relative volatilities can be calculated from the Peng-Robinson Equation of State with interaction parameters assumed to be zero (see Chapter 4). Determine the rank order of the distillation sequences on the basis of total vapor load for ... 

Calculating the Characteristic Interaction Parameter of the Micellar Systems Used To perform the calculation of p12 for the systems examined, i.e. Sulfonate/Genapol/ethoxylated nonylphenol mixtures, the following assumptions were made ... 

The value of the characteristic interaction parameter of these systems (30° C), adjusted from the CMC measurements in Figure 1, was calculated by means of RST and taken equal to -2.5. This value is effectively in the range of the ones found by Graciaa for similar anionic/nonionic mixtures (8). 

The most important aspect of the simulation is that the thermodynamic data of the chemicals be modeled correctly. It is necessary to decide what equation of state to use for the vapor phase (ideal gas, Redlich-Kwong-Soave, Peng-Robinson, etc.) and what model to use for liquid activity coefficients [ideal solutions, solubility parameters, Wilson equation, nonrandom two liquid (NRTL), UNIFAC, etc.]. See Sec. 4, Thermodynamics. It is necessary to consider mixtures of chemicals, and the interaction parameters must be predictable. The best case is to determine them from data, and the next-best case is to use correlations based on the molecular weight, structure, and normal boiling point. To validate the model, the computer results of vapor-liquid equilibria could be checked against experimental data to ensure their validity before the data are used in more complicated computer calculations. 

In this equation, all molecules are divided into four groups paraffins (P), olefins (O), naphthenics (N), and aromatics (A). The v values represent the volume fractions of each component used, while the fa values are the blending values, which were calculated for each of the molecular lumps shown in Table 2. Pure component octane numbers used are designated as ON/, but one should note that in the development of the model, 57 molecular lumps were made based on GC analysis, and pure component ONs were assigned to each lump, and not necessarily each pure component. The kt values are calculated interaction parameters between paraffins, olefins, and naphthenics, and are also shown in Table 2. Based on this equation, and knowing the composition and pure octane numbers of a fuel mixture, an estimation of the blending ON may then be made. 

The optimum UNIQUAC interaction parameters u, between methylcyclohexane, methanol, and ethylbenzene were determined using the observed liquid-liquid data, where the interaction parameters describe the interaction energy between molecules i and j or between each pair of compounds. Table 4 show the calculated value of the UNIQUAC binary interaction parameters for the mixture methanol + ethylbenzene using universal values for the UNIQUAC structural parameters. The equilibrium model was optimized using an objective function, which was developed by Sorensen [15],... 

In monoethanolamine solutions the unknown interaction parameters and equilibrium constants were determined by fitting the model to data for the three component systems CC +MEA+ O and H2S+MEA+H2O. The agreement of the fitted model with Che data was found to be good. The parameters obtained in this way were then used to predict the partial pressures of mixtures of HoS and CO2 over aqueous MEA solutions. The predictions were in good agreement with experimental data, except at the higher partial pressures. 

In this case data for mixtures of H2S+CO2+DIPA+H2O were used together with the data for H2S+DIPA+H2O and CO2+DIPA+H2O to obtain the interaction parameters and equilibrium constants. The results are shown in Figures 1 and 2 to be in good agreement with the experimental data ( 3). In this case, however, in contrast to the case of MEA, the predictions use parameters evaluated from data for the four component system. 

---