Implementing an Equation of State without Derivatives: teqp

Efficient Approximation with Space Filling Quadtrees: Application to Phase Equilibria in Binary Mixtures

A relatively succinct set of property calculations can be used to construct efficient (in memory and speed) numerical approximations of that function in a rectangular domain. Through the use of adaptive subdivision with quadtrees and bi-Chebyshev expansions in each leaf, the function can be practically represented to the order of the noise in the function to be approximated. A further benefit is that evaluation of the approximation is noniterative and thus cannot fail to converge (a relatively common problem in thermophysical property libraries, especially for mixtures). Evaluation of the approximation function requires only a few bisection steps to identify the leaf of interest such that evaluation of the approximation data structure takes less than a microsecond. The technique is demonstrated by application to the vapor-liquid-equilibria evaluated with two different models (COSMO-SAC activity coefficient model and multifluid model). For the more expensive COSMO-SAC case, the approximation function is more than 2000 times faster to evaluate, and deviations in pressure are less than a part in 108 which is practically equal to the iteration convergence criterion.

Influence of repulsion on entropy scaling and density scaling of monatomic fluids

Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.

An Analysis of the Critical Region of Multiparameter Equations of State

In this work, two classes of defects with multiparameter equations of state are investigated. In the first, it is shown that the critical point provided by equation of state developers often does not exactly meet the criticality conditions based on the first two density derivatives of the pressure being zero at the critical point. Based on the more accurate locations of the critical points given in the first part, the scaling of the densities along the binodal and spinodal in the critical region are investigated, and we find that the vast majority of equations have reasonable behavior but a few do not.

The NIST REFPROP Database for Highly Accurate Properties of Industrially Important Fluids

The NIST REFPROP software program is a powerful tool for calculating thermophysical properties of industrially important fluids, and this manuscript describes the models implemented in, and features of, this software. REFPROP implements the most accurate models available for selected pure fluids and their mixtures that are valid over the entire fluid range including gas, liquid, and supercritical states, with the goal of uncertainties approaching the level of the underlying experimental data. The equations of state for thermodynamic properties are primarily of the Helmholtz energy form; a variety of models are implemented for the transport properties. We document the models for the 147 fluids included in the current version. A graphical user interface generates tables and provides extensive plotting capabilities. Properties can also be accessed through third-party apps or user-written code via the core property subroutines compiled into a shared library. REFPROP disseminates international standards in both the natural gas and refrigeration industries, as well as standards for water/steam.