Speed of Sound Measurements in Helium at Pressures from 15 to 100 MPa and Temperatures from 273 to 373 K.

The speed of sound in helium was measured along five isotherms in a temperature range from 273 to 373 K at pressures from 15 to 100 MPa with a relative expanded uncertainty (k = 2) from 0.02 to 0.04%. A dual-path pulse-echo system was utilized to conduct these measurements. The data were compared with the reference equation of state developed by Ortiz Vega et al. At pressures up to 50 MPa, relative deviations were within the uncertainty of our measurements, while, at higher pressures, increasing negative deviations were observed up to -0.26%. We also compared the results with predictions based on the virial equation of state correct to the seventh virial coefficient, using the ab initio virial coefficients reported recently by Gokul et al., finding agreement to within the experimental uncertainty at all investigated states.

Electroreduction of CO2/CO to C2 Products: Process Modeling, Downstream Separation, System Integration, and Economic Analysis.

Direct electrochemical reduction of CO2 to C2 products such as ethylene is more efficient in alkaline media, but it suffers from parasitic loss of reactants due to (bi)carbonate formation. A two-step process where the CO2 is first electrochemically reduced to CO and subsequently converted to desired C2 products has the potential to overcome the limitations posed by direct CO2 electroreduction. In this study, we investigated the technical and economic feasibility of the direct and indirect CO2 conversion routes to C2 products. For the indirect route, CO2 to CO conversion in a high temperature solid oxide electrolysis cell (SOEC) or a low temperature electrolyzer has been considered. The product distribution, conversion, selectivities, current densities, and cell potentials are different for both CO2 conversion routes, which affects the downstream processing and the economics. A detailed process design and techno-economic analysis of both CO2 conversion pathways are presented, which includes CO2 capture, CO2 (and CO) conversion, CO2 (and CO) recycling, and product separation. Our economic analysis shows that both conversion routes are not profitable under the base case scenario, but the economics can be improved significantly by reducing the cell voltage, the capital cost of the electrolyzers, and the electricity price. For both routes, a cell voltage of 2.5 V, a capital cost of $10,000/m2, and an electricity price of <$20/MWh will yield a positive net present value and payback times of less than 15 years. Overall, the high temperature (SOEC-based) two-step conversion process has a greater potential for scale-up than the direct electrochemical conversion route. Strategies for integrating the electrochemical CO2/CO conversion process into the existing gas and oil infrastructure are outlined. Current barriers for industrialization of CO2 electrolyzers and possible solutions are discussed as well.

Reference Materials for Phase Equilibrium Studies. 1. Liquid-Liquid Equilibria (IUPAC Technical Report).

This article is the first of three projected IUPAC Technical Reports resulting from IUPAC Project 2011-037-2-100 (Reference Materials for Phase Equilibrium Studies). The goal of that project was to select reference systems with critically evaluated property values for the validation of instruments and techniques used in phase equilibrium studies for mixtures. This Report proposes seven systems for liquid-liquid equilibrium studies, covering the four most common categories of binary mixtures: aqueous systems of moderate solubility, non-aqueous systems, systems with low solubility, and systems with ionic liquids. For each system, the available literature sources, accepted data, smoothing equations, and estimated uncertainties are given.

Residual entropy model for predicting the viscosities of dense fluid mixtures.

In this work, we have investigated the mono-variant relationship between the reduced viscosity and residual entropy in pure fluids and in binary mixtures of hydrocarbons and hydrocarbons with dissolved carbon dioxide. The mixtures considered were octane + dodecane, decane + carbon dioxide, and 1,3-dimethylbenzene (m-xylene) + carbon dioxide. The reduced viscosity was calculated according to the definition of Bell, while the residual entropy was calculated from accurate multi-parameter Helmholtz-energy equations of state and, for mixtures, the multi-fluid Helmholtz energy approximation. The mono-variant dependence of reduced viscosity upon residual molar entropy was observed for the pure fluids investigated, and by incorporating two scaling factors (one for reduced viscosity and the other for residual molar entropy), the data were represented by a single universal curve. To apply this method to mixtures, the scaling factors were determined from a mole-fraction weighted sum of the pure-component values. This simple model was found to work well for the systems investigated. The average absolute relative deviation (AARD) was observed to be between 1% and 2% for pure components and a mixture of similar hydrocarbons. Larger deviations, with AARDs of up to 15%, were observed for the asymmetric mixtures, but this compares favorably with other methods for predicting the viscosity of such systems. We conclude that the residual-entropy concept can be used to estimate the viscosity of mixtures of similar molecules with high reliability and that it offers a useful engineering approximation even for asymmetric mixtures.

High Pressure Electrochemical Reduction of CO2 to Formic Acid/Formate: A Comparison between Bipolar Membranes and Cation Exchange Membranes.

A high pressure semicontinuous batch electrolyzer is used to convert CO2 to formic acid/formate on a tin-based cathode using bipolar membranes (BPMs) and cation exchange membranes (CEMs). The effects of CO2 pressure up to 50 bar, electrolyte concentration, flow rate, cell potential, and the two types of membranes on the current density (CD) and Faraday efficiency (FE) for formic acid/formate are investigated. Increasing the CO2 pressure yields a high FE up to 90% at a cell potential of 3.5 V and a CD of ∼30 mA/cm2. The FE decreases significantly at higher cell potentials and current densities, and lower pressures. Up to 2 wt % formate was produced at a cell potential of 4 V, a CD of ∼100 mA/cm2, and a FE of 65%. The advantages and disadvantages of using BPMs and CEMs in electrochemical cells for CO2 conversion to formic acid/formate are discussed.

Thermophysical Properties and Phase Behavior of Fluids for Application in Carbon Capture and Storage Processes.

Phase behavior and thermophysical properties of mixtures of carbon dioxide with various other substances are very important for the design and operation of carbon capture and storage (CCS) processes. The available empirical data are reviewed, together with some models for the calculation of these properties. The systems considered in detail are, first, mixtures of carbon dioxide, water, and salts; second, carbon dioxide-rich nonelectrolyte mixtures; and third, mixtures of carbon dioxide with water and amines. The empirical data and the plethora of available models permit the estimation of key fluid properties required in the design and operation of CCS processes. The engineering community would benefit from the further development, and delivery in convenient form, of a small number of these models sufficient to encompass the component slate and operating conditions of CCS processes.

Atomistic Molecular Dynamics Simulations of Carbon Dioxide Diffusivity in n-Hexane, n-Decane, n-Hexadecane, Cyclohexane, and Squalane.

Atomistic molecular dynamics simulations were carried out to obtain the diffusion coefficients of CO2 in n-hexane, n-decane, n-hexadecane, cyclohexane, and squalane at temperatures up to 423.15 K and pressures up to 65 MPa. Three popular models were used for the representation of hydrocarbons: the united atom TraPPE (TraPPE-UA), the all-atom OPLS, and an optimized version of OPLS, namely, L-OPLS. All models qualitatively reproduce the pressure dependence of the diffusion coefficient of CO2 in hydrocarbons measured recently, and L-OPLS was found to be the most accurate. Specifically for n-alkanes, L-OPLS also reproduced the measured viscosities and densities much more accurately than the original OPLS and TraPPE-UA models, indicating that the optimization of the torsional potential is crucial for the accurate description of transport properties of long chain molecules. The three force fields predict different microscopic properties such as the mean square radius of gyration for the n-alkane molecules and pair correlation functions for the CO2-n-alkane interactions. CO2 diffusion coefficients in all hydrocarbons studied are shown to deviate significantly from the Stokes-Einstein behavior.

Kinetics of carbonate mineral dissolution in CO2-acidified brines at storage reservoir conditions.

We report experimental measurements of the dissolution rate of several carbonate minerals in CO2-saturated water or brine at temperatures between 323 K and 373 K and at pressures up to 15 MPa. The dissolution kinetics of pure calcite were studied in CO2-saturated NaCl brines with molalities of up to 5 mol kg-1. The results of these experiments were found to depend only weakly on the brine molality and to conform reasonably well with a kinetic model involving two parallel first-order reactions: one involving reactions with protons and the other involving reaction with carbonic acid. The dissolution rates of dolomite and magnesite were studied in both aqueous HCl solution and in CO2-saturated water. For these minerals, the dissolution rates could be explained by a simpler kinetic model involving only direct reaction between protons and the mineral surface. Finally, the rates of dissolution of two carbonate-reservoir analogue minerals (Ketton limestone and North-Sea chalk) in CO2-saturated water were found to follow the same kinetics as found for pure calcite. Vertical scanning interferometry was used to study the surface morphology of unreacted and reacted samples. The results of the present study may find application in reactive-flow simulations of CO2-injection into carbonate-mineral saline aquifers.

Experimental and modeling study of the phase behavior of (methane + CO2 + water) mixtures.

In this work we report phase equilibrium measurements on the system (methane + carbon dioxide + water) carried out with a high-pressure quasi-static-analytical apparatus. The measurements have been made under conditions of two-phase vapor-liquid equilibrium, three-phase vapor-liquid-liquid equilibrium (VLLE), and four-phase vapor-liquid-liquid-hydrate equilibrium. The compositions of three coexisting fluid phases have been obtained along eight isotherms at temperatures from (285.15 to 303.5) K and at pressures up to either the upper critical end point (UCEP) or up to the hydrate formation locus. Compositions of coexisting vapor and liquid phases have been obtained along three isotherms at temperatures from (323.15 to 423.15) K and pressures up to 20 MPa. The quadruple curve, along which hydrates coexist with the three fluid phases, was also measured along its entire length. The VLLE data obtained for this mixture have been compared with the predictions of the statistical associating fluid theory for potentials of variable range (SAFT-VR), implemented with the square-well potential and using parameters fitted to pure-component and binary-mixture data. Specifically, we used the SAFT-VR parameters reported by Míguez and co-workers [Míguez, J. M.; dos Ramos, M. C.; Piñeiro, M. M.; Blas, F. J. J. Phys. Chem. B 2011, 115, 9604]. The pressure along the quadruple curve was compared with the predictions of two different thermodynamic models. Furthermore, a detailed study of the ternary mixtures was carried out based on comparison with available ternary data of the type (CO2 + n-alkane + water) and available data for the constituent binary subsystems. In this way, we analyzed the observed effects on the solubility when the n-alkane component was changed or a third component was added.

Physical apparatus parameters and model for vibrating tube densimeters at pressures to 140 MPa and temperatures to 473 K.

Vibrating tube densimeters are well-established tools for measuring fluid densities precisely at elevated temperatures and pressures. However, the conventional method of calibrating them utilises a model in which the apparatus parameters are represented as polynomials of temperature and pressure that contain a variable number of terms. Here a robust, physically-based model is presented and demonstrated for six different instruments at temperatures from (273 to 473) K, pressures from (0 to 140) MPa, and densities from (0 to 1050) kg m(-3). The model's physical basis ensures that only seven apparatus parameters are required to relate the measured resonant period to fluid mass density with an average r.m.s. deviation of ±0.23 kg m(-3) across all six densimeters. Estimates for each of the apparatus parameters were made based on the geometry and material properties of the vibrating tubes, and these estimates were consistent with the parameter values determined by calibration with reference fluids. Three of the apparatus parameters describe the temperature dependence of the resonant period: for the six vibrating tubes tested, the relative standard deviations of these parameters were all within the range of values estimated from the thermoelastic properties of the Hastelloy tubes. Two distinct parameters are required to describe the pressure dependence of the vibrating tube's volume and effective spring constant, both of which are estimable from equations describing the elastic deformation of thick-walled tubes. The extensive calibrations conducted demonstrate that, for these densimeters, the variations with pressure of the tube's spring constant and its volume have a ratio that is neither 0 nor 1, as has been assumed previously. The model's physical basis allows vibrating tube densimeters to be calibrated accurately using fewer reference fluid measurements than required by the conventional method. Furthermore, use of the physically-based model reduces the uncertainty of measurements made at densities, temperatures, or pressures beyond the range of the calibration.

Molecular dynamics simulations of CO2 and brine interfacial tension at high temperatures and pressures.

Molecular dynamics simulations have been performed to study the interfacial tension of CO2 and brine for a range of temperatures between 303 and 393 K and pressures from 2 to 50 MPa. The ions involved in this study are Na(+), Ca(2+), and Cl(-). The results indicate that the interfacial tension decreases with increasing pressure under any temperature condition but increases linearly with the molality of the salt solution. The density profiles calculated from the MD simulation results also indicate a positive excess of CO2 and a negative excess of ions at the interface. The charge of the ions was found to have a larger influence than their size on the interfacial tension, a result that consistent with experimental findings.

Viscosity of liquid mixtures: the Vesovic-Wakeham method for chain molecules.

New expressions for the viscosity of liquid mixtures, consisting of chain-like molecules, are derived by means of Enskog-type analysis. The molecules of the fluid are modelled as chains of equally sized, tangentially joined, and rigid spheres. It is assumed that the collision dynamics in such a fluid can be approximated by instantaneous collisions. We determine the molecular size parameters from the viscosity of each pure species and show how the different effective parameters can be evaluated by extending the Vesovic-Wakeham (VW) method. We propose and implement a number of thermodynamically consistent mixing rules, taking advantage of SAFT-type analysis, in order to develop the VW method for chain molecules. The predictions of the VW-chain model have been compared in the first instance with experimental viscosity data for octane-dodecane and methane-decane mixtures, thus, illustrating that the resulting VW-chain model is capable of accurately representing the viscosity of real liquid mixtures.

Experimental and molecular modeling study of the three-phase behavior of (n-decane + carbon dioxide + water) at reservoir conditions.

Knowledge of the phase behavior of mixtures of oil with carbon dioxide and water is essential for reservoir engineering, especially in the processes of enhanced oil recovery and geological storage of carbon dioxide. However, for a comprehensive understanding, the study of simpler systems needs to be completed. In this work the system (n-decane + carbon dioxide + water) was studied as a model (oil + carbon dioxide + water) mixture. To accomplish our aim, a new analytical apparatus to measure phase equilibria at high pressure was designed with maximum operating temperature and pressure of 423 K and 45 MPa, respectively. The equipment relies on recirculation of two coexisting phases using a two-channel magnetically operated micropump designed during this work, with sampling and online compositional analysis by gas chromatography. The apparatus has been validated by comparison with published isothermal vapor-liquid equilibrium data for the binary system (n-decane + carbon dioxide). New experimental data have been measured for the system (n-decane + carbon dioxide + water) under conditions of three-phase equilibria. Data for the three coexisting phases have been obtained on five isotherms at temperatures from 323 to 413 K and at pressures up to the point at which two of the phases become critical. The experimental work is complemented here with a theoretical effort in which we developed models for these molecules within the framework of the statistical associating fluid theory for potentials of variable range (SAFT-VR). The phase behavior of the three binary subsystems was calculated using this theory, and where applicable, a modification of the Hudson and McCoubrey combining rules was used to treat the systems predictively. The experimental data obtained for the ternary mixture are compared to the predictions of the theory. Furthermore, a detailed analysis of the ternary mixture is carried out based on comparison with available data for the constituent binary subsystems. In this way, we analyzed the observed effects on the solubility when the third component was added.

Application of a renormalization-group treatment to the statistical associating fluid theory for potentials of variable range (SAFT-VR).

An accurate prediction of phase behavior at conditions far and close to criticality cannot be accomplished by mean-field based theories that do not incorporate long-range density fluctuations. A treatment based on renormalization-group (RG) theory as developed by White and co-workers has proven to be very successful in improving the predictions of the critical region with different equations of state. The basis of the method is an iterative procedure to account for contributions to the free energy of density fluctuations of increasing wavelengths. The RG method has been combined with a number of versions of the statistical associating fluid theory (SAFT), by implementing White's earliest ideas with the improvements of Prausnitz and co-workers. Typically, this treatment involves two adjustable parameters: a cutoff wavelength L for density fluctuations and an average gradient of the wavelet function Φ. In this work, the SAFT-VR (variable range) equation of state is extended with a similar crossover treatment which, however, follows closely the most recent improvements introduced by White. The interpretation of White's latter developments allows us to establish a straightforward method which enables Φ to be evaluated; only the cutoff wavelength L then needs to be adjusted. The approach used here begins with an initial free energy incorporating only contributions from short-wavelength fluctuations, which are treated locally. The contribution from long-wavelength fluctuations is incorporated through an iterative procedure based on attractive interactions which incorporate the structure of the fluid following the ideas of perturbation theories and using a mapping that allows integration of the radial distribution function. Good agreement close and far from the critical region is obtained using a unique fitted parameter L that can be easily related to the range of the potential. In this way the thermodynamic properties of a square-well (SW) fluid are given by the same number of independent intermolecular model parameters as in the classical equation. Far from the critical region the approach provides the correct limiting behavior reducing to the classical equation (SAFT-VR). In the critical region the ß critical exponent is calculated and is found to take values close to the universal value. In SAFT-VR the free energy of an associating chain fluid is obtained following the thermodynamic perturbation theory of Wertheim from the knowledge of the free energy and radial distribution function of a reference monomer fluid. By determining L for SW fluids of varying well width a unique equation of state is obtained for chain and associating systems without further adjustment of critical parameters. We use computer simulation data of the phase behavior of chain and associating SW fluids to test the accuracy of the new equation.

A kinetic theory description of the viscosity of dense fluids consisting of chain molecules.

An expression for the viscosity of a dense fluid is presented that includes the effect of molecular shape. The molecules of the fluid are approximated by chains of equal-sized, tangentially jointed, rigid spheres. It is assumed that the collision dynamics in such a fluid can be approximated by instantaneous collisions between two rigid spheres belonging to different chains. The approach is thus analogous to that of Enskog for a fluid consisting of rigid spheres. The description is developed in terms of two molecular parameters, the diameter sigma of the spherical segment and the chain length (number of segments) m. It is demonstrated that an analysis of viscosity data of a particular pure fluid alone cannot be used to obtain independently effective values of both sigma and m. Nevertheless, the chain lengths of n-alkanes are determined by assuming that the diameter of each rigid sphere making up the chain can be represented by the diameter of a methane molecule. The effective chain lengths of n-alkanes are found to increase linearly with the number C of carbon atoms present. The dependence can be approximated by a simple relationship m=1+(C-1)3. The same relationship was reported within the context of a statistical associating fluid theory equation of state treatment of the fluid, indicating that both the equilibrium thermodynamic properties and viscosity yield the same value for the chain lengths of n-alkanes.

Vapor pressure and density of thermotropic liquid crystals: MBBA, 5CB, and novel fluorinated mesogens.

The vapor pressures and densities of six thermotropic liquid crystal-forming molecules (mesogens) have been determined experimentally as functions of temperature. The ubiquitous mesogenic compounds n-(4-methoxybenzylidene)-4-butylaniline (MBBA) and 4'-cyano-4-n-pentylbiphenyl (5CB), which both exhibit room-temperature nematic phases, are examined in this study, as are a number of trifluorinated bicyclohexyl and cyclohexylbiphenyl derivatives which find modern use in display applications. Although thermotropic mesogens are of prime importance in modern optoelectronic technologies, there is a scarcity of reliable saturation pressure data for such systems. An apparatus suitable for measurements of vapor pressures between 0.1 and 1333 Pa in the temperature range 298-523 K has been constructed. The adequacy of the apparatus has been verified by measurements on n-hexadecane at temperatures between 304 and 372 K, corresponding to vapor pressures between 0.4 and approximately 100 Pa. To our knowledge, our measurements represent the first reliable data for the saturation pressure of the fluid phase of these types of thermotropic compounds; we show that existing data for MBBA is thermodynamically inconsistent. The densities of the fluid phases of these compounds are also measured by means of a glass pycnometer at temperatures between 293 and 368 K.