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.

Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

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.

Partial collapse of a ventilated stone façade: Diagnosis and analysis of the anchorage system

The mechanical behaviour of transventilated façades performed by natural stone is necessarily based on the correct execution of both anchoring elements on the stone cladding as in the ones corresponding to the enclosure support, either with brick masonry walls or reinforced concrete walls. In the case studied in the present work, the origin of the damages suffered on the façade of a building located in Alcoy has been analyzed, where the detachment of part of the outer enclosure occurred. This enclosure is a transventilated façade formed by Bateig Blue stone tiles. To this end, “in situ” tests of the anchoring systems employed have been performed, as well as laboratory tests of mechanical characterization of the material and of different types of anchor, comparing these results with those obtained in both the simplified analytical models of continuum mechanics as developed by the Finite Element Method (FEM).

Difference schemes for numerical solutions of lagging models of heat conduction

Non-Fourier models of heat conduction are increasingly being considered in the modeling of microscale heat transfer in engineering and biomedical heat transfer problems. The dual-phase-lagging model, incorporating time lags in the heat flux and the temperature gradient, and some of its particular cases and approximations, result in heat conduction modeling equations in the form of delayed or hyperbolic partial differential equations. In this work, the application of difference schemes for the numerical solution of lagging models of heat conduction is considered. Numerical schemes for some DPL approximations are developed, characterizing their properties of convergence and stability. Examples of numerical computations are included.

An improved procedure for obtaining and maintaining well characterized partial water saturation states on concrete samples to be used for mass transport tests

A conditioning procedure is proposed allowing to install into the concrete specimens any selected value of water saturation degree with homogeneous moisture distribution. This is achieved within the least time and the minimum alteration of the concrete specimens. The protocol has the following steps: obtaining basic drying data at 50 °C (water absorption capacity and drying curves); unidirectional drying of the specimens at 50 °C until reaching the target saturation degree values; redistribution phase in closed containers at 50 °C (with measurement of the quasi-equilibrium relative humidities); storage into controlled environment chambers until and during mass transport tests, if necessary. A water transport model is used to derive transport parameters of the tested materials from the drying data, i.e., relative permeabilities and apparent water diffusion coefficients. The model also allows calculating moisture profiles during isothermal drying and redistribution phases, thus allowing optimization of the redistribution times for obtaining homogeneous moisture distributions.

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.

Disipation Functions of Flow Equations in Models of Complex Systems

In an open system, each disequilibrium causes a force. Each force causes a flow process, these being represented by a flow variable formally written as an equation called flow equation, and if each flow tends to equilibrate the system, these equations mathematically represent the tendency to that equilibrium. In this paper, the authors, based on the concepts of forces and conjugated fluxes and dissipation function developed by Onsager and Prigogine, they expose the following hypothesis: Is replaced in Prigogine’s Theorem the flow by its equation or by a flow orbital considering conjugate force as a gradient. This allows to obtain a dissipation function for each flow equation and a function of orbital dissipation.

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.

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.

A Quantile-Based Probabilistic Mean Value Theorem

For non-negative random variables with finite means we introduce an analogous of the equilibrium residual-lifetime distribution based on the quantile function. This allows us to construct new distributions with support (0, 1), and to obtain a new quantile-based version of the probabilistic generalization of Taylor's theorem. Similarly, for pairs of stochastically ordered random variables we come to a new quantile-based form of the probabilistic mean value theorem. The latter involves a distribution that generalizes the Lorenz curve. We investigate the special case of proportional quantile functions and apply the given results to various models based on classes of distributions and measures of risk theory. Motivated by some stochastic comparisons, we also introduce the “expected reversed proportional shortfall order”, and a new characterization of random lifetimes involving the reversed hazard rate function.