GE models and algorithms for condensed phase equilibrium data regression in ternary systems: limitations and proposals

Phase equilibrium data regression is an unavoidable task necessary to obtain the appropriate values for any model to be used in separation equipment design for chemical process simulation and optimization. The accuracy of this process depends on different factors such as the experimental data quality, the selected model and the calculation algorithm. The present paper summarizes the results and conclusions achieved in our research on the capabilities and limitations of the existing GE models and about strategies that can be included in the correlation algorithms to improve the convergence and avoid inconsistencies. The NRTL model has been selected as a representative local composition model. New capabilities of this model, but also several relevant limitations, have been identified and some examples of the application of a modified NRTL equation have been discussed. Furthermore, a regression algorithm has been developed that allows for the advisable simultaneous regression of all the condensed phase equilibrium regions that are present in ternary systems at constant T and P. It includes specific strategies designed to avoid some of the pitfalls frequently found in commercial regression tools for phase equilibrium calculations. Most of the proposed strategies are based on the geometrical interpretation of the lowest common tangent plane equilibrium criterion, which allows an unambiguous comprehension of the behavior of the mixtures. The paper aims to show all the work as a whole in order to reveal the necessary efforts that must be devoted to overcome the difficulties that still exist in the phase equilibrium data regression problem.

Correlation of the liquid–liquid equilibrium data for specific ternary systems with one or two partially miscible binary subsystems

This paper presents the results of a liquid–liquid equilibrium data correlation for 11 ternary systems which have not been previously fitted using the NRTL model or, when they have, the results presented in the literature are inconsistent with the experimental behavior of the system. These ternary systems include mixtures with one or two partially miscible pairs. During the correlation process, new restrictions were imposed on the values for the NRTL binary parameters to ensure correct prediction of the total or partial miscibility for the binary pairs involved. In addition, topological concepts related to the Gibbs stability test have been applied in order to validate the results in the whole range of compositions.

Modelling liquid–liquid equilibria for island type ternary systems

Closed miscibility gaps in ternary liquid mixtures, at constant temperature and pressure, are obtained if phase separations occur only in the ternary region, whilst all binary mixtures involved in the system are completely miscible. This type of behaviour, although not very frequent, has been observed for a certain number of systems. Nevertheless, we have found no information about the applicability of the common activity coefficient models, as NRTL and UNIQUAC, for these types of ternary systems. Moreover, any of the island type systems published in the most common liquid–liquid equilibrium data collections, are correlated with any model. In this paper, the applicability of the NRTL equation to model the LLE of island type systems is assessed using topological concepts related to the Gibbs stability test. A first attempt to correlate experimental LLE data for two island type ternary systems is also presented.

Some examples of modeling in chemical engineering: thermal treatment of polymers, phase equilibrium calculations and process design

Presentation submitted to PSE Seminar, Chemical Engineering Department, Center for Advanced Process Design-making (CAPD), Carnegie Mellon University, Pittsburgh (USA), October 2012.

Geometrical study of the three-dimensional temperature–composition diagrams. An important aid to understand the behavior of the liquid–vapour equilibrium in ternary systems

In this work the usefulness of qualitatively studying and drawing three-dimensional temperature–composition diagrams for ternary systems is pointed out to understand and interpret the particular behavior of the liquid–vapour equilibrium of non-ideal ternary systems. Several examples have been used in order to highlight the interest and the possibilities of this tool, which should be an interesting support not only for lecturers, but also for researchers interested in experimental equilibrium data determination.

Liquid–Liquid, Vapor–Liquid, and Vapor–Liquid–Liquid Equilibrium Data for the Water–n-Butanol–Cyclohexane System at Atmospheric Pressure: Experimental Determination and Correlation

The temperature and the composition of the vapor–liquid–liquid equilibrium (VLLE) and the vapor–liquid equilibrium (VLE) of a ternary mixture of water–n-butanol–cyclohexane were measured at atmospheric pressure (101.32 kPa) in a modified dynamic recirculating still. As found in the literature, the experimental data obtained reveal a ternary azeotrope at 341.86 K with a mole fraction composition of 0.281, 0.034, and 0.685 water, n-butanol, and cyclohexane, respectively. The liquid–liquid equilibrium (LLE) compositions were measured at a constant temperature of 313.15 K and compared with data in the literature collected at other temperatures. Thermodynamic consistency of all the experimental data was demonstrated. The universal quasichemical (UNIQUAC) and the nonrandom two-liquid (NRTL) thermodynamic models were used to correlate the VLE and LLE data, while the original universal functional (UNIFAC) model was used to compare the predicted data.

Solubility, Density, Refractive Index, and Viscosity for the Polyhydric Alcohol + CsBr + H2O Ternary Systems at Different Temperatures

The solubility, density, refractive index, and viscosity data for the ethylene glycol + CsBr + H2O, 1,2-propanediol + CsBr + H2O, and glycerin + CsBr + H2O ternary systems have been determined at (288.15, 298.15, and 308.15) K. In all cases, the solubility of CsBr in aqueous solutions was decreased significantly due to the presence of polyhydric alcohol. The liquid–solid equilibrium experimental data were correlated using the NRTL (nonrandom two-liquid) activity coefficient model, considering nondissociation of the dissolved salt in the liquid phase, and new interaction parameters were estimated. The mean deviations between calculated and experimental compositions were low, showing the good descriptive quality and applicability of the NRTL model. The refractive indices, densities, and viscosities for the unsaturated solutions of the three ternary systems have also been measured at three temperatures. Values for all of the properties were correlated with the salt concentrations and proportions of polyhydric alcohol in the solutions.

Gibbs Energy of Mixing Function: Topological Analysis in Azeotropic Systems

ESAT 2014. 27th European Symposium on Applied Thermodynamics, Eindhoven University of Technology, July 6-9, 2014.

Isothermal (liquid + liquid) equilibrium data at T = 313.15 K and isobaric (vapor + liquid + liquid) equilibrium data at 101.3 kPa for the ternary system (water + 1-butanol + p-xylene)

The (vapor + liquid), (liquid + liquid) and (vapor + liquid + liquid) equilibria of the ternary system (water + 1-butanol + p-xylene) have been determined. (Water + 1-butanol + p-xylene) is a type 2 heterogeneous ternary system with partially miscible (water + 1-butanol) and (water + p-xylene) pairs. By contrast, (1-butanol + p-xylene) is totally miscible under atmospheric conditions. This paper examines the (vapor + liquid) equilibrium in both heterogeneous and homogeneous regions at 101.3 kPa of pressure. (Liquid + liquid) equilibrium data at T = 313.15 K have also been determined, and for comparison, the obtained experimental data have been calculated by means of several thermodynamic models: UNIQUAC, UNIFAC and NRTL. Some discrepancies were found between the (vapor + liquid + liquid) correlations; however, the models reproduced the (liquid + liquid) equilibrium data well. The obtained data reveal a ternary heterogeneous azeotrope with mole fraction composition: 0.686 water, 0.146 1-butanol and 0.168 p-xylene.

Procedure for calculating the chloride diffusion coefficient and surface concentration from a profile having a maximum beyond the concrete surface

Chlorides induce local corrosion in the steel reinforcements when reaching the bar surface. The measurement of the rate of ingress of these ions, is made by mathematically fitting the so called “error function equation” into the chloride concentration profile, obtaining so the diffusion coefficient and the chloride concentration at the concrete surface. However, the chloride profiles do not always follow Fick’s law by having the maximum concentration at the concrete surface, but often the profile shows a maximum concentration more in the interior, which indicates a different composition and performance of the most external concrete layer with respect to the internal zones. The paper presents a procedure prepared during the time of the RILEM TC 178-TMC: “Testing and modeling chloride penetration in concrete”, which suggests neglecting the external layer where the chloride concentration increases and using the maximum as an “apparent” surface concentration, called C max and to fit the error function equation into the decreasing concentration profile towards the interior. The prediction of evolution should be made also from the maximum.

An R script to plot the conditioned Gini (or skewness) coefficient in the bivariate case

In the cs.index.zip file we provide an R script which let us plot the conditioned Gini (or skewness) coefficient used in the working paper entitled "On conditional skewness with applications in environmental data" submitted to Environmental and Ecological Statistics. On the other hand, the ReadMe.pdf explains how to use the cs.index.R script.

New methodology for describing the equilibrium beach profile applied to the Valencia's beaches

Mathematical models used for the understanding of coastal seabed morphology play a key role in beach nourishment projects. These projects have become the fundamental strategy for coastal maintenance during the last few years. Accordingly, the accuracy of these models is vital to optimize the costs of coastal regeneration projects. Planning of such interventions requires methodologies that do not generate uncertainties in their interpretation. A study and comparison of mathematical simulation models of the coastline is carried out in this paper, as well as elements that are part of the model that are a source of uncertainty. The equilibrium profile (EP) and the offshore limit corresponding to the depth of closure (DoC) have been analyzed taking into account different timescale ranges. The results have thus been compared using data sets from three different periods which are identified as present, past and future. Accuracy in data collection for the beach profiles and the definition of the median grain size calculation using collected samples are the two main factors that have been taken into account in this paper. These data can generate high uncertainties and can produce a lack of accuracy in nourishment projects. Together they can generate excessive costs due to possible excess or shortage of sand used for the nourishment. The main goal of this paper is the development of a new methodology to increase the accuracy of the existing equilibrium beach profile models, providing an improvement to the inputs used in such models and in the fitting of the formulae used to obtain seabed shape. This new methodology has been applied and tested on Valencia's beaches.

Isobaric vapor–liquid–liquid–solid equilibrium of the water + NaCl + 1-butanol system at 101.3 kPa

A mixture of water + NaCl + 1-butanol at 101.3 kPa is studied in order to determine the influence of salt on its experimental vapor–liquid–liquid–solid equilibrium. A detailed analysis of the evolution with temperature of the different equilibrium regions is carried out. The study is conducted at a constant pressure of 101.3 kPa in a recirculating still that has been modified by our research group. The changes in the 1-butanol/water composition ratio in the vapor phase that are provoked by the salt are studied as a function of equilibrium region. In addition, the mutual solubility of 1-butanol and water is assessed in the liquid–liquid and solid–liquid regions.