The prediction of isobaric phase equilibria for non-ideal multicomponent mixtures from heat of mixing data

A method for predicting isobaric binary and ternary vapor—liquid equilibrium data using only isothermal binary heat of mixing data and pure component vapor pressure data is presented. Three binary and two ternary hydrocarbon liquid mixtures were studied. The method consists of evaluating the parameters of the NRTL equation from isothermal heat of mixing data for the constituent binary pairs. These parameters are then used in the multicomponent NRTL equation to compute isobaric vapor—liquid equilibrium data for the ternary mixture. No ternary or higher order interaction terms are needed in the ternary calculations because of the nature of the NRTL equation. NRTL parameters derived from heat of mixing data at one temperature can be used to predict vapor—liquid equilibrium data at other temperatures up to the boiling temperature of the liquid mixture.For the systems studied this method predicted the composition of the vapor phase with a standard deviation ranging from 1–8% for the binary systems and from 4–12% for the ternary systems.     

Modeling Binary and Multicomponent Systems Containing Supercritical CO2 with Polyethylene Glycols and Compounds Relevant to the Biodiesel Production

The CPA equation of state is applied to model binary, ternary, and multicomponent mixtures that contain CO₂ with polyethylene glycols or compounds relevant to biodiesel production, such as glycerol and various triglycerides. Effort has been made to evaluate the model performance on correlating both the liquid and the vapor phase compositions, which is a demanding task, revealing the model’s and parameters’ limitations, due to the rather low concentrations of heavy compounds in the vapor phase. Initially the model’s binary parameters, which in all cases were temperature independent, were estimated using experimental data for binary systems. Those parameters were used to predict the phase behavior of supercritical CO₂ containing ternary and multicomponent mixtures. Since no parameter was adjusted to ternary or multicomponent systems’ data, the reported CPA results for such mixtures are considered as pure predictions. This is the final part of a series of studies [Tsivintzelis et al. Fluid Phase Equilibria 430 (2016) 75–92 and 504 (2020) 112337] that complete the parameterization of the CPA equation of state for systems relevant to the biodiesel production, which allows the application of the model to multicomponent mixtures of the relevant processes.     

Thermodynamics of aqueous solutions containing poly (N-isopropylacrylamide)

Hydrogels undergo reversible and discontinuous volume changes in response to variation of solution conditions such as solvent composition, temperature, salt concentration, and pH. In this contribution we focus our attention on the experimental and theoretical investigation of these swelling equilibria of aqueous cross-linked poly (N-isopropylacrylamide) solutions as well as on the connected demixing behavior of the linear polymer dissolved in water. For the experimental study of the (liquid + liquid) equilibrium an alternative method based on refractive index measurements is suggested. In order to calculate the swelling behavior a model combining an expression for the Gibbs free energy of mixing with an expression for the elastic network is applied. As a model for the Gibbs free energy of mixing the UNIQUAC-approach and the Koningsveld–Kleintjens model are used. For the elastic network contribution again two different theories, namely the phantom network theory and the affine network theory, were applied. Whereas the type of network theory has only a small influence on the calculation results, the Gibbs free energy of mixing has a large impact. Using the UNIQUAC-approach the swelling equilibria can be correlated close to the experimental data, however, this model predicts a homogeneous mixture for linear polymer chains in water. In contrast to this situation the Koningsveld–Kleintjens model does a good job in calculating the swelling equilibria as well as the demixing curve, however, the adjustable parameter must be changed slightly.     

Modeling the phase equilibria of hydrogen sulfide and carbon dioxide in mixture with hydrocarbons and water using the PCP-SAFT equation of state

The perturbed-chain polar statistical associating fluid theory (PCP-SAFT) equation of state is applied to correlate phase equilibria for mixtures of hydrogen sulfide (H₂S) and carbon dioxide (CO₂) with alkanes, with aromatics, and with water over wide temperature and pressure ranges. The binary mixtures of H₂S–methane and CO₂–methane are studied in detail including vapor–liquid, liquid–liquid and fluid–solid phase equilibria. Very satisfying results were obtained for the binary mixtures as well as for the ternary mixture of H₂S–CO₂–methane using the (constant) interaction parameters of the binary pairs.     

Vapor–liquid equilibrium of methane and methane + nitrogen and an equimolar hexane + decane mixture under isothermal conditions

Experimental vapor–liquid equilibrium data of the ternary system composed of methane and an equimolar hexane + decane mixture are reported. The experimental measurements were carried out under isothermal conditions at 258, 273, and 298 K in the pressure range 1–19 MPa. Also, experimental vapor–liquid measurements were carried out for the quaternary system methane + nitrogen and an equimolar hexane + decane mixture, at 258 K in the range 3.5–12 MPa. The results for the ternary system show that the solubility of methane in the equimolar mixture of alkanes increases when the pressure is increased at constant temperature and it increases as the temperature decreases in the whole pressure range studied. For the quaternary system with a constant amount of nitrogen, the solubility of methane in the liquid phase increases as the pressure increases at the studied temperature. The experimental results for the ternary system were satisfactorily correlated with the Peng–Robinson equation of state in the ranges of pressure and temperature studied. The equation of state was used to predict the behavior of the quaternary system using binary interaction parameters. The applicability of the principle of congruence was corroborated by comparing the vapor–liquid behavior of methane in the equimolar hexane + decane mixture with that in pure octane, at the three temperatures studied in this work.     

Molecular thermodynamics on the basis of an equation of state for argon. III. Systems of appreciably nonspherical linear molecules

A molecular model based on an equation of state for pure argon is applied to various systems involving appreciably nonspherical molecules. Since the pure components are not described adequately, their potential parameters are fitted at each temperature separately to ensure the proper vapor pressures and liquid volumes. Used in this way, the model makes reasonable predictions of the thermodynamic behavior of the binary mixtures considered. Again, the introduction of a binary parameter in the dispersion potential does not lead to any significant overall improvement.     

Deviations from ideal sublimation vapor pressure behavior in mixtures of polycyclic aromatic compounds with interacting heteroatoms

Despite the relatively small atomic fraction of a given heteroatom in a binary mixture of polycyclic aromatic compounds (PAC), the presence of heteroatoms can significantly impact mixture vapor pressure behavior over a wide range of temperatures. The vapor pressures of several binary PAC mixtures containing various heteroatoms show a range of behavior, from practically ideal behavior following Raoult’s law to significant deviations from ideality depending on the heteroatoms present. Mixtures were synthesized using the quench-cool technique with equimolar amounts of two PAC, both containing heteroatoms such as aldehyde, carboxyl, nitrogen, and sulfur substituent groups. For some mixtures, deviation from ideality is inversely related to temperature, though in other cases we see deviation from ideality increasing with temperature, whereas some appear independent of temperature. Most commonly we see lower vapor pressures than predicted by Raoult’s Law, which indicates that the interacting heteroatoms prefer the solid mixture phase as opposed to the vapor phase. Although negative deviations predominate from Raoult’s Law, the varying mixtures investigated show both higher and lower enthalpies and entropies of sublimation than predicted. In each mixture, a higher enthalpy of sublimation leads to higher entropy of sublimation than predicted, and vice versa.     

Carbon dioxide-methane mixture adsorption on activated carbon

In this work, we report new experimental data of pure and binary adsorption equilibria of carbon dioxide and methane on the activated carbon RB2 at 273 and 298 K. The pressure range studied were 0–3.5 MPa for pure gases and 0–0.1 MPa for mixtures. The combination of the generalized Dubinin model to describe the pure CO₂ and CH₄ isotherms with the IAST (Ideal Adsorbed Solution Theory) for the mixtures provide a method for the calculation of the binary adsorption equilibria. This formulation predicts with acceptable accuracy the binary adsorption data and can easily be integrated in general dynamic simulation of PSA (pressure swing adsorption process) adsorption columns. It involves only three parameters, independent of the temperature, and directly determined with only one adsorption isotherm of CO₂.     

Critical points of hydrocarbon mixtures with the Peng–Robinson,SAFT, and PC-SAFT equations of state

The phase behavior of fluids at high pressures can be rather complex, even for mixtures of relatively simple molecules, such as hydrocarbons. In this work, we use the Hicks and Young algorithm to calculate mixture critical points, comparing five modeling options: Peng–Robinson EOS: (1) original and (2) with parameters fitted from molar volume and vapor pressure data; (3) SAFT EOS; and PC-SAFT EOS: (4) original and (5) with refitted parameters to match pure component critical data. Calculations were carried out for binary hydrocarbon mixtures and 29 multicomponent mixtures. The SAFT EOS provided the worst representation of the systems tested and, interestingly, the conventional cubic EOS provided, in general, the best representation.     

Comparison of two hard-sphere reference systems for perturbation theories for mixtures

The Carnahan—Starling equation of state for hard spheres can be extended to mixtures using either a one-fluid theory, or the generalization of scaled-particle (or Percus—Yevick theory) proposed by Boublik and by Mansoori and coworkers. The two reference systems are combined with a perturbation term of the van der Waals form; they are then used to correlate the phase behavior of binary mixtures of nonpolar molecules differing significantly in molecular size. In each case, one adjustable binary parameter (a₁₂) is used to correlate vapor—liquid equilibria over the entire composition range. Predicted Henry's constants and liquid densities for the saturated mixture are compared with experiment. The Boublik—Mansoori hard-sphere-mixture equation is superior to the Carnahan—Starling One-Fluid theory, expecially in the dilute region.     

Investigation of the surface tension of methane and n-alkane mixtures by perturbed-chain statistical associating fluid theory combined with density-gradient theory

The perturbed-chain statistical associating fluid theory (PC-SAFT) and density-gradient theory are used to construct an equation of state to describe the phase behavior of binary methane–n-alkane mixtures. With the molecular parameters and influence parameters regressed from bulk properties and surface tensions of pure fluids, respectively as input, both the bulk and interfacial properties are investigated. The surface tension of the binary systems methane–propane, methane–pentane, methane–heptane and methane–decane are predicted, and the results are satisfactory compared with the experimental data. Our results show that PC-SAFT combined with density-gradient theory is able to describe the interfacial properties of binary methane–n-alkane mixtures in wide temperature and pressure ranges, and illustrate the influence of the equilibrium bulk properties and chain length of n-alkane molecule on the interfacial properties.     

Simulation of phase separation in alcohol/water mixtures using two-body force field and standard molecular dynamics

Standard molecular dynamics simulations have been carried out on pure alcohols and alcohol/water mixtures. A simple atom-atom force field consisting of Lennard-Jones potentials plus coulombic terms over atomic point charges, but without explicit polarization terms, has been specifically fitted to reproduce several experimental properties of the pure alcohols, and has been used for mixtures by developing combination rules with the TIP3P water model. Densities, enthalpies of vaporization, radial distribution functions, self-diffusion coefficients, and rotational correlation functions of the pure alcohols are well reproduced and compare favorably with those from more sophisticated force fields. Some key aspects of the phase behaviour are correctly reproduced by the molecular dynamics simulation, showing a distinct demixing process for the n-butanol/water mixture as opposed to the stability of the t-butanol/water mixtures. The results demonstrate the ability of a molecular dynamics simulation, even in its standard form and with easily accessible time ranges, but with a carefully optimized force field, to simulate and, to a certain extent, predict the properties of binary mixtures.     

Equation of state modeling of the phase equilibria of ionic liquid mixtures at low and high pressure

Accurate design of processes based on ionic liquids (ILs) requires knowledge of the phase behavior of the systems involved. In this work, the truncated perturbed chain polar statistical associating fluid theory (tPC-PSAFT) is used to correlate the phase behavior of binary and ternary IL mixtures. Both non-polar and polar solvents are examined, while methyl imidazolium ILs are used in all cases. tPC-PSAFT accounts explicitly for weak dispersion interactions, highly directive polar interactions between permanent dipolar and quadrupolar molecules and association between hydrogen bonding molecules. For mixtures of non-polar solvents, tPC-PSAFT predicts accurately the binary mixture data. For the case of polar solvents, a binary interaction parameter is fitted to the experimental data and the agreement between experiment and correlation is very good in all cases.     

A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures

A new approach has been developed for calculating the properties of mixtures based on an equation of state explicit in reduced Helmholtz energy. This approach allows for the representation of the thermodynamic properties over a wide range of fluid states and is based on highly accurate equations of state for the pure components combined at the reduced temperature and density of the mixture. The reducing parameters used for temperature and density depend on composition. For simple mixtures (those that closely follow Raoult's law), a very accurate representation of all thermodynamic properties has been achieved with relatively simple functions. For nonideal mixtures, the reducing functions for density and temperature were modified, and a departure function was added to the equation of state. Generally, the model is able to represent liquid and vapor states with uncertainties of 0.1% in density, 1% in heat capacities and 1% in bubble point pressures if experimental data of comparable uncertainties exist. Two applications of the mixture model concepts were developed independently by the authors in the United States and Germany over the same time period. These applications include the development of individual equations for each binary system and a generalization of the model which is valid for a wide variety of mixtures. The individual approaches are presented with an explanation of the similarities and differences. Although the paper focuses mainly on binary systems, some results for ternary mixtures are also presented.     

Modeling phase equilibria for acid gas mixtures using the CPA equation of state. Part II: Binary mixtures with CO2

In Part I of this series of articles, the study of H₂S mixtures has been presented with CPA. In this study the phase behavior of CO₂ containing mixtures is modeled. Binary mixtures with water, alcohols, glycols and hydrocarbons are investigated. Both phase equilibria (vapor-liquid and liquid-liquid) and densities are considered for the mixtures involved. Different approaches for modeling pure CO₂ and mixtures are compared. CO₂ is modeled as non self-associating fluid, or as self-associating component having two, three and four association sites. Moreover, when mixtures of CO₂ with polar compounds (water, alcohols and glycols) are considered, the importance of cross-association is investigated. The cross-association is accounted for either via combining rules or using a cross-solvation energy obtained from experimental spectroscopic or calorimetric data or from ab initio calculations. In both cases two adjustable parameters are used when solvation is explicitly accounted for. The performance of CPA using the various modeling approaches for CO₂ and its interactions is presented and discussed, comparatively to various recent published investigations. It is shown that overall very good correlation is obtained for binary mixtures of CO₂ and water or alcohols when the solvation between CO₂ and the polar compound is explicitly accounted for, whereas the model is less satisfactory when CO₂ is treated as self-associating compound.     

Comparison of three models for permeation of CO2/CH4 mixtures in poly(phenylene oxide)

Two models for the permeability of pure gases have been extended to include binary gas mixtures. The first is an extension of a pure gas permeability model, proposed by Petropoulos, which is based on gradients of chemical potential. This model predicts the permeability of components in a gas mixture solely on the basis of competition for sorption sites within the polymer matrix. The second mixed gas model follows an earlier analysis by Barrer for pure gases which includes the effects of saturation of Langmuir sites on the diffusion as well as the sorption processes responsible for permeation. This generalized “competitive sorption/diffusion” model includes the effect of each gas component on the sorption and diffusion of the other component in the mixture. The flux equations from these two models have been solved numerically to predict the permeability of gas mixtures on the basis of pure gas sorption and transport parameters. Both the mixed gas Petropoulos and competitive sorption/diffusion model predictions are compared with predictions from the earlier simple competitive sorption model based on gradients of concentration. An analysis of all three models is presented for the case of CO₂/CH₄ permeability in poly(phenylene oxide) (PPO). As expected, the competitive sorption/diffusion model predicts lower permeability than either of the models which consider only competitive sorption effects. The permeability depression of both CO₂ and CH₄ predicted by the competitive sorption/diffusion model is roughly twice that predicted by the competitive sorption model, whereas the mixed gas Petropoulos model predictions for both gases lie between the other two model predictions. For the PPO/CO₂/CH₄ system, the methane permeability data lie above the predictions of all three models, whereas CO₂ data lie below the predictions of all models. Consequently, the competitive sorption/diffusion model gives the most accurate prediction for CO₂, while the simple competitive sorption model is best for methane. The effects of mixed gas sorption, fugacity, and CO₂-induced dilation were considered and do not explain the inaccuracies of any of the models. The relatively small errors in mixed gas permeability predictions using either of the three models are likely to be related to “transport plasticization” of PPO owing to high levels of CO₂ sorption and its effect on polymer segmental motions and gas diffusivity.     

Study on volumetric measurement and correlations for MTBE+1-propanol+decane at 298.15 K and atmospheric pressure

Summary Experimental densities for the ternary mixture x₁MTBE+x₂1-propanol+(1-x₁-x₂)decane and the binary mixtures xMTBE +(1-x)1-propanol and x1-propanol+(1-x)decane have been measured at 298.15 K and atmospheric pressure, using a DMA 4500 Anton Paar densimeter. Excess molar volumes were determined from the densities of the pure liquids and mixtures. Attending to the symmetry of the studied mixtures, suitable fitting equations have been used in order to correlate adequately the experimental data. For the ternary mixture, experimental data were also used to test several empirical expressions for estimating ternary properties from experimental binary results.     

On the thermodynamics of alcohol solutions. Phase equilibria of binary and ternary mixtures containing any number of alcohols

Nagata, I., 1985. On the thermodynamics of alcohol solutions. Phase equilibria of binary and ternary mixtures containing any number of alcohols. Fluid Phase Equilibria, 19: 153–174.Binary vapor—liquid and liquid—liquid equilibrium data for alcohol solutions includin one or two alcohols are correlated with the UNIQUAC associated solution theory (Nagata and Kawamura). The theory uses pure liquid association constants determined by the method of Brandani and a single value of the enthalpy of the hydrogen bond equal to ?23.2 kJ mol ^?1 for pure alcohols. For alcohol-active nonassociating component mixtures and alcohol—alcohol mixtures the theory involves additional solvation constants. The theory is extended to contain ternary mixtures with any number of alcohols. Ternary predictions of vapor—liquid and liquid—liquid equilibria are performed using only binary parameters. Good agreement is obtained between calculated and experimental results for many representative mixtures.