Vapor-liquid equilibria using the Gibbs energy and the common tangent plane criterion

Phase thermodynamics is often perceived as a difficult subject that many students never become fully comfortable with. The Gibbsian geometrical framework can help students to gain a better understanding of phase equilibria. An exercise to interpret the vapor-liquid equilibrium of a binary azeotropic mixture, using the equilibrium condition based on the common tangent plane criterion (the Gibbs stability test), is presented in this paper. From a T-composition phase diagram for the mixture, the temperature is set at different values: above, intermediate to, and below the boiling temperatures of the pure components, to intersect different regions of the system. Students prepare an Excel spreadsheet where the Gibbs energy of mixing of the vapor and liquid mixtures are calculated and represented over the whole range of compositions and then, apply the Gibbs stability test to ascertain the aggregation state of the system and to calculate the VL phase equilibrium compositions. Finally, Matlab is used to generate the 3D Gibbs energy of mixing surfaces for both phases over the whole range of temperatures which facilitates the geometrical interpretation of the vapor-liquid equilibrium.

Pitfalls on computing liquid-liquid phase equilibria using the k-value method

Poster presented in the 11th Mediterranean Congress of Chemical Engineering, Barcelona, October 21-24, 2008.

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.

Simultaneous synthesis of work exchange networks with heat integration

The optimal integration of work and its interaction with heat can represent large energy savings in industrial plants. This paper introduces a new optimization model for the simultaneous synthesis of work exchange networks (WENs), with heat integration for the optimal pressure recovery of process gaseous streams. The proposed approach for the WEN synthesis is analogous to the well-known problem of synthesis of heat exchanger networks (HENs). Thus, there is work exchange between high-pressure (HP) and low-pressure (LP) streams, achieved by pressure manipulation equipment running on common axes. The model allows the use of several units of single-shaft-turbine-compressor (SSTC), as well as stand-alone compressors, turbines and valves. Helper motors and generators are used to respond to any demand and excess of energy. Moreover, between the WEN stages the streams are sent to the HEN to promote thermal recovery, aiming to enhance the work integration. A multi-stage superstructure is proposed to represent the process. The WEN superstructure is optimized in a mixed-integer nonlinear programming (MINLP) formulation and solved with the GAMS software, with the goal of minimizing the total annualized cost. Three examples are conducted to verify the accuracy of the proposed method. In all case studies, the heat integration between WEN stages is essential to improve the pressure recovery, and to reduce the total costs involved in the process.

Integration of modular process simulators under the Generalized Disjunctive Programming framework for the structural flowsheet optimization

The optimization of chemical processes where the flowsheet topology is not kept fixed is a challenging discrete-continuous optimization problem. Usually, this task has been performed through equation based models. This approach presents several problems, as tedious and complicated component properties estimation or the handling of huge problems (with thousands of equations and variables). We propose a GDP approach as an alternative to the MINLP models coupled with a flowsheet program. The novelty of this approach relies on using a commercial modular process simulator where the superstructure is drawn directly on the graphical use interface of the simulator. This methodology takes advantage of modular process simulators (specially tailored numerical methods, reliability, and robustness) and the flexibility of the GDP formulation for the modeling and solution. The optimization tool proposed is successfully applied to the synthesis of a methanol plant where different alternatives are available for the streams, equipment and process conditions.

An alternative disjunctive optimization model for heat integration with variable temperatures

This paper presents an alternative model to deal with the problem of optimal energy consumption minimization of non-isothermal systems with variable inlet and outlet temperatures. The model is based on an implicit temperature ordering and the “transshipment model” proposed by Papoulias and Grossmann (1983). It is supplemented with a set of logical relationships related to the relative position of the inlet temperatures of process streams and the dynamic temperature intervals. In the extreme situation of fixed inlet and outlet temperatures, the model reduces to the “transshipment model”. Several examples with fixed and variable temperatures are presented to illustrate the model's performance.

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.

Biodiesel wastes: An abundant and promising source for the preparation of acidic catalysts for utilization in etherification reaction

Environmentally friendly sulfonated black carbon (BC) catalysts were prepared from biodiesel waste, glycerol. These black carbons (BCs) contain a high amount of acidic groups, mainly sulfonated and oxygenated groups. Furthermore, these catalysts show a high catalytic activity in the glycerol etherification reaction with tert-butyl alcohol, the activity being larger for the sample prepared with a higher glycerol:sulfuric acid ratio (1:3). The yield for mono-tert-butyl glycerol (MTBG), di-tert-butyl glycerol (DTBG) and tri-tert-butyl-glycerol (TTBG) were very similar to those obtained using a commercial resin, Amberlyst-15. Furthermore, experimental results show that the carbon with the lowest acidic surface group content, BC prepared in minor glycerol:sulfuric acid ratio (10:1), can be chemically treated after carbonization to achieve an improved catalytic activity. The activity of all BCs is high and very similar, about 50% and 20% for the MTBG and DTBG + TTBG, respectively.

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.

Improving some thermodynamic consistency tests for VLE data evaluation

Thermodynamics Conference 2013 (Statistical Mechanics and Thermodynamics Group of the Royal Society of Chemistry), The University of Manchester, 3-6 September 2013.

Diseño de estrategias conjuntas UA-centros de secundaria para el fomento del conocimiento de disciplinas científico-técnicas

Continuando con la línea emprendida por el Instituto Universitario de Ingeniería de Procesos Químicos, para desarrollar actividades de fomento del conocimiento de disciplinas de ciencia y tecnología, se ha convocado el “VI Certamen de Proyectos Educativos de Ingeniería Química” durante el curso 2013-2014. Con dicha actividad se pretende promover el interés por estudios relacionados con las disciplinas científico-técnicas, entre los alumnos de ESO, con una mejor conceptualización, comprensión y caracterización de los temas. Para ello, se diseñan, planifican y desarrollan acciones entre profesorado de la UA y centros de secundaria, sobre el papel de la Química en la industria para mejorar la calidad de vida. En la presente edición del certamen se han inscrito 13 centros, que han presentado un total de 20 proyectos, relacionados con los siguientes temas: la Ingeniería Química y el medio ambiente, la Ingeniería Química y la industria alimentaria, la Ingeniería Química y el mundo de los plásticos, la Ingeniería Química y la energía y la Ingeniería Química «Verde» sostenible.

Trabajo de coordinación para la implementación del cuarto curso del Grado en Química

En el curso 2013-14 se implantará el cuarto y último curso del Grado en Química. La experiencia adquirida durante la implementación de los tres primeros cursos ha puesto de manifiesto la conveniencia de realizar un proceso de coordinación entre los profesores de las diferentes asignaturas que garantice la consecución de las competencias previstas en la memoria del título. Por ello, se ha creado en la Facultad de Ciencias de la Universidad de Alicante una red de investigación en docencia universitaria que ha estado trabajando desde el inicio del presente curso académico en este tema. Dicha red está constituida por el Vicedecano de Ordenación Académica de la Facultad de Ciencias, la Coordinadora Académica de Química y los profesores coordinadores de todas las asignaturas del 4º curso del grado (excepto Prácticas Externas y Trabajo Fin de Grado). En esta comunicación se presentarán los resultados del trabajo de investigación realizado por estos profesores que ha permitido elaborar las guías docentes de las asignaturas, planificar y coordinar las actividades a realizar para que los alumnos adquieran las competencias transversales, realizar un cronograma de actividades de evaluación y otro de prácticas de laboratorio que asegure la distribución homogénea del trabajo del alumno durante el curso académico.

A "Know-How vs. Know-What" Approach in the Teaching-Learning of Competences in Physical Chemistry

The methodological approach a teacher uses in the competence teaching-learning process determines the way students learn. Knowledge can be acquired from a series of perspectives, mainly: “know-what” (concept), where facts and descriptions of (natural or social) phenomena are pursued; “know-how” (procedure), where methods and procedures for their application are described; and “know-why” (competence), where general principles and laws that explain both the facts and their applications are sought. As all the three cases are interconnected, the boundaries between them are not fully clear and their application uses shared elements. In any case, the depth of student’s acquired competences will be directly affected by the teaching-learning perspective, traditionally aiming to a “know-why” approach for full competence acquisition. In this work, we discuss a suitable teaching-learning methodology for evaluating whether a “know-how”, “know-what” or combined approach seems better for enhancing competence learning in students. We exemplify the method using a selection of formative activities from the Physical Chemistry area in the Grades of Chemistry and Chemical Engineering.

Isobutane Alkylation Process Synthesis by means of Hybrid Simulation-Multiobjective Optimization

Multiobjective Generalized Disjunctive Programming (MO-GDP) optimization has been used for the synthesis of an important industrial process, isobutane alkylation. The two objective functions to be simultaneously optimized are the environmental impact, determined by means of LCA (Life Cycle Assessment), and the economic potential of the process. The main reason for including the minimization of the environmental impact in the optimization process is the widespread environmental concern by the general public. For the resolution of the problem we employed a hybrid simulation- optimization methodology, i.e., the superstructure of the process was developed directly in a chemical process simulator connected to a state of the art optimizer. The model was formulated as a GDP and solved using a logic algorithm that avoids the reformulation as MINLP -Mixed Integer Non Linear Programming-. Our research gave us Pareto curves compounded by three different configurations where the LCA has been assessed by two different parameters: global warming potential and ecoindicator-99.