No abstract

Equations of state (EoS) are essential in the modeling of a wide range of industrial and natural processes. Desired qualities of EoS are accuracy, consistency, computational speed, robustness, and predictive ability outside of the domain where they have been fitted. In this work, we review present challenges associated with established models, and give suggestions on how to overcome them in the future. The most accurate EoS available, multiparameter EoS, have a second artificial Maxwell loop in the two-phase region that gives problems in phase-equilibrium calculations and excludes them from important applications such as treatment of interfacial phenomena with mass-based density functional theory. Suggestions are provided on how this can be improved. Cubic EoS are among the most computationally efficient EoS, but they often lack sufficient accuracy. We show that extended corresponding state EoS are capable of providing significantly more accurate single-phase predictions than cubic EoS with only a doubling of the computational time. In comparison, the computational time of multiparameter EoS can be orders of magnitude larger. For mixtures in the two-phase region, however, the accuracy of extended corresponding state EoS has a large potential for improvement. The molecular-based SAFT family of EoS is preferred when predictive ability is important, for example, for systems with strongly associating fluids or polymers where few experimental data are available. We discuss some of their benefits and present challenges. A discussion is presented on why predictive thermodynamic models for reactive mixtures such as CO2–NH3 and CO2–H2O–H2S must be developed in close combination with phase- and reaction equilibrium theory, regardless of the choice of EoS. After overcoming present challenges, a next-generation thermodynamic modeling framework holds the potential to improve the accuracy and predictive ability in a wide range of applications such as process optimization, computational fluid dynamics, treatment of interfacial phenomena, and processes with reactive mixtures.

Abstract. We establish the existence of small-amplitude uni-and bimodal steady periodic gravity waves with an affine vorticity distribution, using a bifurcation argument that differs slightly from earlier theory. The solutions describe waves with critical layers and an arbitrary number of crests and troughs in each minimal period. An important part of the analysis is a fairly complete description of the local geometry of the so-called kernel equation, and of the small-amplitude solutions. Finally, we investigate the asymptotic behavior of the bifurcating solutions.

We present a perturbation theory that combines the use of a third-order Barker–Henderson expansion of the Helmholtz energy with Mie-potentials that include first- (Mie-FH1) and second-order (Mie-FH2) Feynman–Hibbs quantum corrections. The resulting equation of state, the statistical associating fluid theory for Mie potentials of variable range corrected for quantum effects (SAFT-VRQ-Mie), is compared to molecular simulations and is seen to reproduce the thermodynamic properties of generic Mie-FH1 and Mie-FH2 fluids accurately. SAFT-VRQ Mie is exploited to obtain optimal parameters for the intermolecular potentials of neon, helium, deuterium, ortho-, para-, and normal-hydrogen for the Mie-FH1 and Mie-FH2 formulations. For helium, hydrogen, and deuterium, the use of either the first- or second-order corrections yields significantly higher accuracy in the representation of supercritical densities, heat capacities, and speed of sounds when compared to classical Mie fluids, although the Mie-FH2 is slightly more accurate than Mie-FH1 for supercritical properties. The Mie-FH1 potential is recommended for most of the fluids since it yields a more accurate representation of the pure-component phase equilibria and extrapolates better to low temperatures. Notwithstanding, for helium, where the quantum effects are largest, we find that none of the potentials give an accurate representation of the entire phase envelope, and its thermodynamic properties are represented accurately only at temperatures above 20 K. Overall, supercritical heat capacities are well represented, with some deviations from experiments seen in the liquid phase region for helium and hydrogen.

In the Lennard-Jones spline (LJ/s) model, the Lennard-Jones (LJ) potential is truncated and splined so that both the pair potential and the force continuously approach zero at rc ≈ 1.74σ. It exhibits the same structural features as the LJ model, but the thermodynamic properties are different due to the shorter range of the potential. One advantage of the model is that simulation times are much shorter. In this work, we present a systematic map of the thermodynamic properties of the LJ/s model from molecular dynamics and Gibbs ensemble Monte Carlo simulations. Accurate results are presented for gas/liquid, liquid/solid and gas/solid coexistence curves, supercritical isotherms up to a reduced temperature of 2 (in LJ units), surface tensions, speed of sound, the Joule-Thomson inversion curve, and the second to fourth virial coefficients. The critical point for the LJ/s model is estimated to be T * c = 0.885 ± 0.002 and P * c = 0.075 ± 0.001, respectively. The triple point is estimated to T * tp = 0.547 ± 0.005 and P * tp = 0.0016 ± 0.0002. Despite the simplicity of the model, the acentric factor was found to be as large as 0.07±0.02. The coexistence densities, saturation pressure, and supercritical isotherms of the LJ/s model were fairly well represented by the Peng-Robison equation of state. We find that Barker-Henderson perturbation theory works much more poorly for the LJ/s model than for the LJ model. The first-order perturbation theory overestimates the critical temperature and pressure by about 10% and 90%, respectively. A second-order perturbation theory that uses the mean compressibility approximation shifts the critical point closer to simulation data, but makes the prediction of the saturation densities worse. It is hypothesized that the reason for this is that the mean compressibility approximation gives a poor representation of the second-order perturbation term for the LJ/s model and that a correction factor is needed at high densities.Our main conclusion is that we at the moment do not have a theory or model that adequately represents the thermodynamic properties of the LJ/s system.

A present key barrier for implementing large-scale hydrogen liquefaction plants is their high power consumption. The cryogenic heat exchangers are responsible for a significant part of the exergy destruction in these plants and we evaluate in this work strategies to increase their efficiency. A detailed model of a plate-fin heat exchanger is presented that incorporates the geometry of the heat exchanger, nonequilibrium ortho-para conversion and correlations to account for the pressure drop and heat transfer coefficients due to possible boiling/condensation of the refrigerant at the lowest temperatures. Based on available experimental data, a correlation for the ortho-para conversion kinetics is developed, which reproduces available experimental data with an average deviation of 2.2%. In a plate-fin heat exchanger that is used to cool the hydrogen from 47.8 K to 29.3 K with hydrogen as refrigerant, we find that the two main sources of exergy destruction are thermal gradients and ortho-para hydrogen conversion, being responsible for 69% and 29% of the exergy destruction respectively. A route to reduce the exergy destruction from the ortho-para hydrogen conversion is to use a more efficient catalyst, where we find that a doubling of the catalytic activity in comparison to ferric-oxide, as demonstrated by nickel oxide-silica catalyst, reduces the exergy destruction by 9%. A possible route to reduce the exergy destruction from thermal gradients is to employ an evaporating mixture of helium and neon at the cold-side of the heat exchanger, which reduces the exergy destruction by 7%. We find that a combination of hydrogen and helium-neon as refrigerants at high and low temperatures respectively, enables a reduction of the exergy destruction by 35%. A combination of both improved catalyst and the use of hydrogen and helium-neon as refrigerants gives the possibility to reduce the exergy destruction in the cryogenic heat exchangers by 43%. The limited efficiency of the ortho-para catalyst represents a barrier for further improvement of the efficiency.

The curvature dependence of the surface tension can be described by the Tolman length (first-order correction) and the rigidity constants (second-order corrections) through the Helfrich expansion. We present and explain the general theory for this dependence for multicomponent fluids and calculate the Tolman length and rigidity constants for a hexane-heptane mixture by use of square gradient theory. We show that the Tolman length of multicomponent fluids is independent of the choice of dividing surface and present simple formulae that capture the change in the rigidity constants for different choices of dividing surface. For multicomponent fluids, the Tolman length, the rigidity constants, and the accuracy of the Helfrich expansion depend on the choice of path in composition and pressure space along which droplets and bubbles are considered. For the hexane-heptane mixture, we find that the most accurate choice of path is the direction of constant liquid-phase composition. For this path, the Tolman length and rigidity constants are nearly linear in the mole fraction of the liquid phase, and the Helfrich expansion represents the surface tension of hexane-heptane droplets and bubbles within 0.1% down to radii of 3 nm. The presented framework is applicable to a wide range of fluid mixtures and can be used to accurately represent the surface tension of nanoscopic bubbles and droplets.

The leading order terms in a curvature expansion of the surface tension, the Tolman length (first order), and rigidities (second order) have been shown to play an important role in the description of nucleation processes. This work presents general and rigorous expressions to compute these quantities for any nonlocal density functional theory (DFT). The expressions hold for pure fluids and mixtures, and reduce to the known expressions from density gradient theory (DGT). The framework is applied to a Helmholtz energy functional based on the perturbed chain polar statistical associating fluid theory (PCP-SAFT) and is used for an extensive investigation of curvature corrections for pure fluids and mixtures. Predictions from the full DFT are compared to two simpler theories: predictive density gradient theory (pDGT), that has a density and temperature dependent influence matrix derived from DFT, and DGT, where the influence parameter reproduces the surface tension as predicted from DFT. All models are based on the same equation of state and predict similar Tolman lengths and spherical rigidities for small molecules, but the deviations between DFT and DGT increase with chain length for the alkanes. For all components except water, we find that DGT underpredicts the value of the Tolman length, but overpredicts the value of the spherical rigidity. An important basis for the calculation is an accurate prediction of the planar surface tension. Therefore, further work is required to accurately extract Tolman lengths and rigidities of alkanols, because DFT with PCP-SAFT does not accurately predict surface tensions of these fluids.