Interaction parameters fitted

Figure 4. Three-phase equilibrium LtL2V in the system carbon dioxide-water-1-propanol at 333 K and 13.1 MPa exp., this work — Calculated with Peng-Robinson EOS using Panagiotopoulos and Reid mixing rule, left side prediction from pure component and binary data alone, right side interaction parameters fitted to ternary three-phase equilibria at temperatures between 303 and 333 K...

In Great Britain, the National Engineering Laboratory (NEL, formerly a government agency but now privatized) has prodnced a database for thermodynamic and transport properties. PPDS contains correlations for properties of a large number of pure components these are based on evaluated experimental data where possible but also include some estimated properties. For mixtures, the database contains binary interaction parameters fitted to data for use with common equation-of-state and liqnid-activity methods for calculating phase eqnilibria. Information is available at their Web site [14]. 

Predictive methods make possible to treat the non-ideality of a liquid mixture without the knowledge of binary interaction parameters fitted from experimental data. Obviously, the predictive should be used only for exploratory purposes. Here we present two approaches. The first one, called the regular solution theory, requires information only about pure components. The second one, UNIFAC, is based on group contributions, and makes use indirectly of experimental data. 

A more difficult and industrially important test for an equation of state is the ability to predict the VLE of ternary and multicomponent mixtures using interactions parameters fitted to corresponding binary data. Tsonopoulos and co-workers analyzed experimental VLE data for 6 ternary mixtures and showed that Redlich-Kwong-Joffe-Zudkevitch is again slightly more precise than the Soave-Redlich-Kwong and Peng-Robinson equations of state. 

Fig. 4. UNIFAC group interaction parameter matrix, the ISi represents parameters fit and parameters not available (168). A represents an aromatic...

To illustrate the use of the gas sorption mode , we show in Figure 7 results of the supercritical ethylene sorption in low-density polyethylene (12,16). As seen in Figure 7, the theory is capable of fitting the ethylene sorption data. In this instance, the data at three temperatures can be fit within experimental precision using interaction parameters (p o) of 3235 atm, 3178 atm, or 3101 atm at 126°C, 140 0, and 155 C, respectively. 

It is well known that cubic equations of state have inherent limitations in describing accurately the fluid phase behavior. Thus our objective is often restricted to the determination of a set of interaction parameters that will yield an "acceptable fit" of the binary VLE data. The following implicit least squares objective function is suitable for this purpose... 

When the fit is judged to be excellent the statistically best interaction parameters can be efficiently obtained by performing implicit ML estimation. This was found to be the case with the methane-methanol and the nitrogen-ethane systems presented later in this chapter. 

Data at two temperatures were obtained from Zeck and Knapp (1986) for the nitrogen-ethane system. The implicit LS estimates of the binary interaction parameters are ka=0, kb=0, kc=0 and kd=0.0460. The standard deviation of kd was found to be equai to 0.0040. The vapor liquid phase equilibrium was computed and the fit was found to be excellent (Englezos et al. 1993). Subsequently, implicit ML calculations were performed and a parameter value of kd=0.0493 with a standard deviation equal to 0.0070 was computed. Figure 14.2 shows the experimental phase diagram as well as the calculated one using the implicit ML parameter estimate. 

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. 

This expression can be modified to apply directly to any of various techniques used to measure the interaction parameter, including membrane and vapor osmometry, freezing point depression, light scattering, viscometry, and inverse gas chromatography [89], A polynomial curve fit is typically used for the concentration dependence of %, while the temperature dependence can usually be fit over a limited temperature range to the form [47]... 

The fit was again good, providing that the interaction parameter was positive for both types of adsorbed hydrogen, indicating an attractive interaction between adsorbed particles and thus a decreasing AHads with increasing 0. 

In contrast to the NRTL-SAC model, the UNIFAC model developed by Fredenslund et. al. [29] divides each molecule into a set of functional groups that interact with each other on a binaiy basis and whose interactions are combined together to describe the global liquid phase interaction between molecules. Because the segments in UNIFAC are based on functional groups it is possible to model a system provided that all of the molecular structures are known. The problem with pharmaceutical sized molecules is that existing UNIFAC parameter tables do not contain many of the group interaction parameters that are necessary, and even when they do, the interactions are fitted to a database of chemicals that are much smaller and simpler than pharmaceuticals, and typically fail to represent them adequately. 

First a database of solute-solvent properties are created in SoluCalc. The database needs the melting point, the enthalpy of fusion and the Hildebrand solubility parameter of the solute (Cimetidine) and the solvents for which solubility data is available. Using the available data, SoluCalc first prepares a list of the most sensitive group interactions and fits sequentially, the solubility data for the minimum set of group interaction parameters that best represent the total data set. For a small set of solvents, the fitted values from SoluCalc are shown in Table 9. It can be noted that while the correlation is very good, the local model is more like a UNIQUAC model than a group contribution model... 

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. 

Figure 6. Dimensionless interaction parameters as determined from isothermal fits at various temperatures ((- -) KBr-water ( ) NaCl—water)...

The interaction parameters are weak, linear functions of temperature, as shown in Table 5, Table 6 and Figure 6. These tables and figure show the results of isothermal fits for activity coefficient data of aqueous NaCl and KBr at various temperatures. The Pitzer equation parameters are, however, strongly dependent on temperature (Silvester and Pitzer, (23)). 

A second type of ternary electrolyte systems is solvent -supercritical molecular solute - salt systems. The concentration of supercritical molecular solutes in these systems is generally very low. Therefore, the salting out effects are essentially effects of the presence of salts on the unsymmetric activity coefficient of molecular solutes at infinite dilution. The interaction parameters for NaCl-C02 binary pair and KCI-CO2 binary pair are shown in Table 8. Water-electrolyte binary parameters were obtained from Table 1. Water-carbon dioxide binary parameters were correlated assuming dissociation of carbon dioxide in water is negligible. It is interesting to note that the Setschenow equation fits only approximately these two systems (Yasunishi and Yoshida, (24)). 

We fitted the interaction parameters with equations as represented in the graphs(8) and given in Table II. Values do not differ significantly from those of the equations of EMNP(3). 

J. As with the alkane - water systems, the interaction parameters for the aqueous liquid phase were found to be temperature - dependent. However, the compositions for the benzene - rich phases could not be accurately represented using any single value for the constant interaction parameter. The calculated water mole fractions in the hydrocarbon - rich phases were always greater than the experimental values as reported by Rebert and Kay (35). The final value for the constant interaction parameter was chosen to fit the three phase locus of this system. Nevertheless, the calculated three-phase critical point was about 9°C lower than the experimental value. 

The interdependence of the Gibbs energy of adsorption and the molecular interaction parameter was recently discussed in detail by Karol-czak, who used a six-parameter model. Contrary to the rather general Damaskin model, no relation between the molecular interaction parameter A and AGads was assumed. It was suggested that this is an arbitrary relation dependent on the theoretical model used in fitting experimental data within acceptable experimental errors. 

D = D° exp(-ac ), where D is the diffusion, D represents the zero-concentration limit, c is the concentration, a and v are parameters, fits the data from a wide variety of probes and matrix polymers ( ). Several theoretical justifications for this behavior have been presented (97-1011. but it is not possible to tell yet which, if any, is uniquely correct. The treatments range from simple physical considerations (98) to treatments of hydrodynsumical interaction of probe and matrix (97,991. Other more complex and general treatments (1001 do not explicitly arrive at the stretched exponential form, but do closely fit the available data. Much more work needs to be done on probe diffusion in such transient networks. Beyond enhancing the arsenal of gel characterization, the problem is quite fundamental to a number of other important processes. 

---