Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Thermodynamic properties of the system toluene — 1, 1, 1-trichloroethane

Data on molar excess enthalpy on mixing at 298.15 K and 308.15 K, vapor-liquid equilibrium, latent heats of vaporization at 91.444 kPa and vapor pressures for the system toluene – 1, 1, 1-trichloroethane are presented. A simple adiabatic calorimeter designed for molar excess enthalpy measurements is described, tested and used. On présente, dans le cas du système toluène – 1, 1, 1-trichloréthane, des résultats relatifs aux grandeurs suivantes: a) enthalpie molaire d'excès à 298.15 K et 308.15 K; b) équilibre liquid-vapeur; c) chaleurs latentes de vaporisation à une pression absolue de 91.444 kP; d) pressions de vapeur. On décrit un calorimètre adiabatique simple, conçu pour mesurer l'enthalpie molaire d'excès, dont on a fait l'essai.

Thermodynamic properties of the systems benzene-chloroethanes

This paper presents the experimental data on vapor-liquid equilibrium and heats of mixing of mixtures of benzene with 1, e-dichloroethane, 1, l, 1 -trichloroethane, and lt1,2,2-tetrachloroethane.A literature survey revealed that the heats of mixing of benzene-l,2-dichloroethane have been studied and Table I shows the extent of study on this system.

Modeling Of Hydrolytic Polymerization In A Semibatch Nylon 6 Reactor

Hydrolytic polymerization of caprolactam to Nylon 6 in a semibatch reactor is carried out by heating a mixture of water and caprolactam. Evaporation of volatiles caused by heating results in a pressure build-up. After the pressure reaches a predetermined value, vapors are vented to keep the pressure constant for some time, and thereafter, to lower the pressure to a value slightly above atmospheric in a preprogrammed manner. The characteristics of the polymer are determined by the chemical reactions and the vaporization of water and caprolactam. The semibatch operation has been simulated and the predictions have been compared with industria data. The observed temperature and pressure histories were predicted with a fair degree of accuracy. It was found that the predictions of the degree of polymerization however are sensitive to the vapor-liquid equilibrium relations. A comparison with an earlier model, which neglected mass transfer resistance, indicates that simulation using the VLE data of Giori and Hayes and accounting for mass transfer resistance is more reliable.

Redução de dimensionalidade em equilíbrio de fases fluidas

Atualmente, existem modelos matemáticos capazes de preverem acuradamente as relações entre propriedades de estado; e esta tarefa é extremamente importante no contexto da Engenharia Química, uma vez que estes modelos podem ser empregados para avaliar a performance de processos químicos. Ademais, eles são de fundamental importância para a simulação de reservatórios de petróleo e processos de separação. Estes modelos são conhecidos como equações de estado, e podem ser usados em problemas de equilíbrios de fases, principalmente em equilíbrios líquido-vapor. Recentemente, um teorema matemático foi formulado (Teorema de Redução), fornecendo as condições para a redução de dimensionalidade de problemas de equilíbrios de fases para misturas multicomponentes descritas por equações de estado cúbicas e regras de mistura e combinação clássicas. Este teorema mostra como para uma classe bem definidade de modelos termodinâmicos (equações de estado cúbicas e regras de mistura clássicas), pode-se reduzir a dimensão de vários problemas de equilíbrios de fases. Este método é muito vantajoso para misturas com muitos componentes, promovendo uma redução significativa no tempo de computação e produzindo resultados acurados. Neste trabalho, apresentamos alguns experimentos numéricos com misturas-testes usando a técnica de redução para obter pressões de ponto de orvalho sob especificação de temperaturas.

Investigation of Interfacial Tensions for Carbon Dioxide Aqueous Solutions by Perturbed-Chain Statistical Associating Fluid Theory Combined with Density-Gradient Theory

The perturbed-chain statistical associating fluid theory and density-gradient theory are used to construct an equation of state (EOS) applicable for the phase behaviors of carbon dioxide aqueous solutions. With the molecular parameters and influence parameters respectively regressed from bulk properties and surface tensions of pure fluids as input, both the bulk and interfacial properties of carbon dioxide aqueous solutions are satisfactorily correlated by adjusting the binary interaction parameter (k(ij)). Our results show that the constructed EOS is able to describe the interfacial properties of carbon dioxide aqueous solutions in a wide temperature range, and illustrate the influences of temperature, pressure, and densities in each phase on the interfacial properties.

Hydrate dissociation conditions for gas mixtures containing carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons using SAFT

A new method, a molecular thermodynamic model based on statistical mechanics, is employed to predict the hydrate dissociation conditions for binary gas mixtures with carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons in the presence of aqueous solutions. The statistical associating fluid theory (SAFT) equation of state is employed to characterize the vapor and liquid phases and the statistical model of van der Waals and Platteeuw for the hydrate phase. The predictions of the proposed model were found to be in satisfactory to excellent agreement with the experimental data.

Bridging the Gap between Ionic Liquids and Molten Salts: Group 1 Metal Salts of the Bistriflamide Anion in the Gas Phase

Fourier transform ion cyclotron resonance mass spectrometry experiments showed that liquid Group 1 metal salts of the bistriflamide anion undergoing reduced-pressure distillation exhibit a remarkable behavior that is in transition between that of the vapor-liquid equilibrium characteristics of aprotic ionic liquids and that of the Group 1 metal halides: the unperturbed vapors resemble those of aprotic ionic liquids, in the sense that they are essentially composed of discrete ion pairs. However, the formation of large aggregates through a succession of ion-molecule reactions is closer to what might be expected for Group I metal halides. Similar experiments were also carried out with bis{(trifluoromethyl)sulfonyl}amine to investigate the effect of H+, which despite being the smallest Group 1 cation, is generally regarded as a nonmetal species. In this case, instead of the complex ion-molecule reaction pattern found for the vapors of Group I metal salts, an equilibrium similar to those observed for aprotic ionic liquids was observed.

Estudo de um sistema híbrido de destilação solar para polimento de águas produzidas

The oil industry, experiencing a great economic and environmental impact, has increasingly invested in researches aiming a more satisfactory treatment of its largest effluent, i.e., produced water. These are mostly discarded at sea, without reuse and after a basic treatment. Such effluent contains a range of organic compounds with high toxicity and are difficult to remove, such as polycyclic aromatic hydrocarbons, salts, heavy metals, etc.. The main objective of this work was to study the solar distillation of produced water pre-treated to remove salts and other contaminants trough of a hybrid system with a pre-heater. This developed apparatus was called solar system, which consists of a solar heater and a conventional distillation solar still. The first device consisted of a water tank, a solar flat plate collector and a thermal reservoir. The solar distillator is of simple effect, with 1m2 of flat area and 20° of inclination. This dissertation was divided in five steps: measurements in the solar system, i.e. temperatures and distillate flow rate and weather data; modeling and simulation of the system; study of vapor-liquid equilibrium of the synthetic wastewater by the aqueous solution of p-xylene; physical and chemical analyses of samples of the feed, distillate and residue, as well as climatology pertinent variables of Natal-RN. The solar system was tested separately, with the supply water, aqueous NaCl and synthetic oil produced water. Temperature measurements were taken every minute of the thermal reservoir, water tank and distillator (liquid and vapor phases). Data of solar radiation and rainfall were obtained from INPE (National Institute for Space Research). The solar pre-heater demonstrated to be effective for the liquid systems tested. The reservoir fluid had an average temperature of 58°C, which enabled the feed to be pre-heated in the distillator. The temperature profile in the solar distillator showed a similar behavior to daily solar radiation, with temperatures near 70°C. The distillation had an average yield of 2.4 L /day, i.e., an efficiency of 27.2%. Mathematical modeling aided the identification of the most important variables and parameters in the solar system. The study of the vapor-liquid equilibrium from Total Organic Carbon (TOC) analysis indicated heteroazeotropia and the vapor phase resulted more concentrated in p-xylene. The physical-chemical analysis of pH, conductivity, Total Dissolved Solids (TDS), chlorides, cations (including heavy metals) and anions, the effluent distillate showed satisfactory results, which presents a potential for reuse. The climatological study indicates the region of Natal-RN as favorable to the operation of solar systems, but the use of auxiliary heating during periods of higher rainfall and cloud cover is also recommended

Determinação de dados de equilíbrio de fases para sistemas aquosos com eletrólitos

In the present work are established initially the fundamental relationships of thermodynamics that govern the equilibrium between phases, the models used for the description of the behavior non ideal of the liquid and vapor phases in conditions of low pressures. This work seeks the determination of vapor-liquid equilibrium (VLE) data for a series of multicomponents mixtures of saturated aliphatic hydrocarbons, prepared synthetically starting from substances with analytical degree and the development of a new dynamic cell with circulation of the vapor phase. The apparatus and experimental procedures developed are described and applied for the determination of VLE data. VLE isobarics data were obtained through a Fischer s ebulliometer of circulation of both phases, for the systems pentane + dodecane, heptane + dodecane and decane + dodecane. Using the two new dynamic cells especially projected, of easy operation and low cost, with circulation of the vapor phase, data for the systems heptane + decane + dodecane, acetone + water, tween 20 + dodecane, phenol + water and distillation curves of a gasoline without addictive were measured. Compositions of the equilibrium phases were found by densimetry, chromatography, and total organic carbon analyzer. Calibration curves of density versus composition were prepared from synthetic mixtures and the behavior excess volumes were evaluated. The VLE data obtained experimentally for the hydrocarbon and aqueous systems were submitted to the test of thermodynamic consistency, as well as the obtained from the literature data for another binary systems, mainly in the bank DDB (Dortmund Data Bank), where the Gibbs-Duhem equation is used obtaining a satisfactory data base. The results of the thermodynamic consistency tests for the binary and ternary systems were evaluated in terms of deviations for applications such as model development. Later, those groups of data (tested and approved) were used in the KijPoly program for the determination of the binary kij parameters of the cubic equations of state original Peng-Robinson and with the expanded alpha function. These obtained parameters can be applied for simulation of the reservoirs petroleum conditions and of the several distillation processes found in the petrochemistry industry, through simulators. The two designed dynamic cells used equipments of national technology for the determination of VLE data were well succeed, demonstrating efficiency and low cost. Multicomponents systems, mixtures of components of different molecular weights and also diluted solutions may be studied in these developed VLE cells

Assessing the activity coefficients of water in cholinium-based ionic liquids: experimental measurements and COSMO-RS modeling

The vapor liquid-equilibrium of water + ionic liquids is relevant for a wide range of applications of these compounds. It is usually measured by ebulliometric techniques, but these are time consuming and expensive. In this work it is shown that the activity coefficients of water in a series of cholinium-based ionic liquids can be reliably and quickly estimated at 298.15K using a humidity meter instrument. The cholinium based ionic liquids were chosen to test this experimental methodology since data for water activities of quaternary ammonium salts are available in the literature allowing the validation of the proposed technique. The COSMO-RS method provides a reliable description of the data and was also used to understand the molecular interactions occurring on these binary systems. The estimated excess enthalpies indicate that hydrogen bonding between water and ionic liquid anion is the dominant interaction that governs the behavior of water and cholinium-based ionic liquids systems, while the electrostatic-misfit and van der Walls forces have a minor contribution to the total excess enthalpies.

Mapping Binary Liquid-Vapor or Liquid-Liquid-Vapor Equilibria Regions, including the Different Azeotropic Behaviours, as a Function of the NRTL Binary Parameters

13th Mediterranean Congress of Chemical Engineering (Sociedad Española de Química Industrial e Ingeniería Química, Fira Barcelona, Expoquimia), Barcelona, September 30-October 3, 2014

Calorimetry In Vapor Absorption And Binary Mixture Vapor Compression Refrigeration And Heatpump Systems

The use of binary fluid systems in thermally driven vapour absorption and mechanically driven vapour compression refrigeration and heatpump cycles has provided an impetus for obtaining experimental date on caloric properties of such fluid mixtures. However, direct measurements of these properties are somewhat scarce in spite of the calorimetric techniques described in the literature being quite adequate. Most of the design data are derived through calculations using theoretical models and vapour-liquid equilibrium data. This article addresses the choice of working fluids and the current status on the data availability vis-a-vis engineering applications. Particular emphasis is on organic working fluid pairs.

Deterpenation of Bergamot Essential Oil Using Liquid-Liquid Extraction: Equilibrium Data of Model Systems at 298.2 K

The deterpenation of bergamot essential oil can be performed by liquid liquid extraction using hydrous ethanol as the solvent. A ternary mixture composed of 1-methyl-4-prop-1-en-2-yl-cydohexene (limonene), 3,7-dimethylocta-1,6-dien-3-yl-acetate (linalyl acetate), and 3,7-dimethylocta-1,6-dien-3-ol (linalool), three major compounds commonly found in bergamot oil, was used to simulate this essential oil. Liquid liquid equilibrium data were experimentally determined for systems containing essential oil compounds, ethanol, and water at 298.2 K and are reported in this paper. The experimental data were correlated using the NRTL and UNIQUAC models, and the mean deviations between calculated and experimental data were lower than 0.0062 in all systems, indicating the good descriptive quality of the molecular models. To verify the effect of the water mass fraction in the solvent and the linalool mass fraction in the terpene phase on the distribution coefficients of the essential oil compounds, nonlinear regression analyses were performed, obtaining mathematical models with correlation coefficient values higher than 0.99. The results show that as the water content in the solvent phase increased, the kappa value decreased, regardless of the type of compound studied. Conversely, as the linalool content increased, the distribution coefficients of hydrocarbon terpene and ester also increased. However, the linalool distribution coefficient values were negatively affected when the terpene alcohol content increased in the terpene phase.

Vapor–liquid equilibrium of binary mixtures containing isopropyl acetate and alkanols at 101.32 kPa

Isobaric vapor-liquid equilibria of binary mixtures of isopropyl acetate plus an alkanol (1-propanol, 2-propanol, 1-butanol, or 2-butanol) were measured at 101.32 kPa, using a dynamic recirculating still. An azeotropic behavior was observed only in the mixtures of isopropyl acetate + 2-propanol and isopropyl acetate + 1-propanol. The application of four thermodynamic consistency tests (the Herington test, the Van Ness test, the infinite dilution test, and the pure component test) showed the high quality of the experimental data. Finally, both NRTL and UNIQUAC activity coefficient models were successfully applied in the correlation of the measured data, with the average absolute deviations in vapor phase composition and temperature of 0.01 and 0.16 K, respectively.