skip to main content

Title: Assessing the Quality of Molecular Simulations for Vapor−Liquid Equilibria: An Analysis of the TraPPE Database
As molecular modeling and simulation techniques become increasingly important sources of thermophysical property and phase equilibrium data, the ability to assess the robustness of that data becomes more critical. Recently, the use of the compressibility factor (Z) has been suggested as a metric for testing the quality of simulation data for vapor−liquid equilibria (VLE). Here, we analyze predicted VLE data from the transferable potentials for phase equilibria (TraPPE) database and show that, apart from data entry or typographical errors, Z will always be well-behaved in Gibbs ensemble Monte Carlo (GEMC) simulations even when the simulations are not sufficiently equilibrated. However, this is not true for grand canonical Monte Carlo simulations. When the pressure is calculated from the internal forces, then pressure and density are strongly correlated for the vapor phase and, for GEMC simulations, it is recommended to treat Z as an instantaneous mechanical property. From analysis of the TraPPE VLE data, we propose a complementary metric based on the predicted vapor pressures at three neighboring temperatures and their deviation from a local Clausius−Clapeyron fit.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Journal of chemical and engineering data
Page Range / eLocation ID:
1330 - 1344
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The studies of transcritical and supercritical flow have attracted much interest in the past 30 years. However, most of them mainly focus on the single-component system, whose critical point is constant. We use the vapor-liquid equilibrium (VLE) theory to capture the thermodynamic properties of the mixture and investigate transcritical flows (i.e., supercritical CO2 oxy-combustion systems). In sCO2 oxy-combustion systems, due to the presence of water from the previous cycles, the mixture critical point increases significantly, such that the phase separation could occur in both the compressor and combustor. However, the VLE solver increases the computation cost of fluid simulation significantly, which limited the size of simulations we can conduct. Naturally, tabulation methods can be used to store the VLE solutions to avoids redundant computation. However, the size of the VLE table increases exponentially with respect to the number of components. When the number of species components is greater than three, the size of the VLE table far exceeds the RAM’s limit in today’s standard computers. In this research, an online tabulation method based on In Situ Adaptive Tabulation (ISAT) is developed to accelerate the computation of multicomponent fluids based on VLE theory. Accuracy and efficiency are analyzed and discussed. The CFD solver used in this research is based on the Pressure-Implicit with Splitting of Operators (PISO) method. Peng-Robinson equation of state (EOS) is used in the calculations of phase equilibrium. 
    more » « less
  2. The formation of subcritical methanol clusters in the vapor phase is known to complicate the analysis of nucleation measurements. Here, we investigate how this process affects the onset of binary nucleation as dilute water–methanol mixtures in nitrogen carrier gas expand in a supersonic nozzle. These are the first reported data for water–methanol nucleation in an expansion device. We start by extending an older monomer–dimer–tetramer equilibrium model to include larger clusters, relying on Helmholtz free energy differences derived from Monte Carlo simulations. The model is validated against the pressure/temperature measurements of Laksmono et al. [Phys. Chem. Chem. Phys. 13, 5855 (2011)] for dilute methanol–nitrogen mixtures expanding in a supersonic flow prior to the appearance of liquid droplets. These data are well fit when the maximum cluster size imax is 6–12. The extended equilibrium model is then used to analyze the current data. On the addition of small amounts of water, heat release prior to particle formation is essentially unchanged from that for pure methanol, but liquid formation proceeds at much higher temperatures. Once water comprises more than ∼24 mol % of the condensable vapor, droplet formation begins at temperatures too high for heat release from subcritical cluster formation to perturb the flow. Comparing the experimental results to binary nucleation theory is challenged by the need to extrapolate data to the subcooled region and by the inapplicability of explicit cluster models that require a minimum of 12 molecules in the critical cluster. 
    more » « less
  3. In recent decades, metal–organic frameworks (MOFs) have gained recognition for their potential in multicomponent gas separations. Though molecular simulations have revealed structure–property relationships of MOF–adsorbate systems, they can be computationally expensive and there is a need for surrogate models that can predict the adsorption data faster. In this work, an active learning (AL) protocol is introduced that can predict multicomponent gas adsorption in a MOF for a range of thermodynamic conditions. This methodology is applied to build a model for the adsorption of three different gas mixtures (CO2–CH4, Xe–Kr, and H2S–CO2) in the MOF Cu-BTC. A Gaussian process regression (GPR) model is used to fit the data as well to leverage its predicted uncertainty to drive the learning. The training data is generated using grand-canonical Monte Carlo (GCMC) simulations as points are iteratively added to the model to minimize the predicted uncertainty. Also, a criteria which captures the perceived performance of the GPs is introduced to terminate the AL process when the perceived accuracy threshold is met. The three systems are tested for a pressure–mole fraction (P–X), and a pressure–mole fraction–temperature (P–X–T) feature space. It is demonstrated that AL one only needs a fraction of the data from simulations to build a reliable surrogate model for predicting mixture adsorption. Further, the final GP fit from AL outperforms ideal adsorbed solution theory predictions. 
    more » « less
  4. To achieve high power density and thermodynamic cycle efficiency, the working pressures of liquid-propellant rocket engines, diesel engines, and gas turbines (based on deflagration or detonation) are continuously increasing, which could reach or go beyond the thermodynamic critical pressure of the liquid propellant. For this reason, the studies of trans- and super-critical injection are getting more and more attention. However, the simulation of transcritical phase change is still a challenging topic. The phase boundary, especially near the mixture critical point, needs to be accurately determined to investigate the multicomponent effects on transcritical injection and atomization. This work used our previously developed thermodynamic model based on the vapor-liquid equilibrium (VLE) theory, which can predict the phase separation near the mixture critical point. An \textit{in situ} adaptive tabulation (ISAT) method was developed to accelerate the computationally expensive multicomponent VLE computation such that it can be cheap enough for CFD. The new thermodynamic model was integrated into OpenFOAM to develop a VLE-based CFD solver. In this work, shock-droplet interaction and two-phase mixing simulations are conducted using our new VLE-based CFD solver. The shock-droplet interaction simulation results capture the thermodynamic condition of the surface entering the supercritical state after shock passes through. The atomization of droplets could be triggered by vorticity formed at the droplets' surface. 2D temporal mixing layer simulations show the evolution of the transcritical mixing layer and capture the phase split effect at the mixing layer. 
    more » « less
  5. Abstract

    Additive manufacturing, no longer reserved exclusively for prototyping components, can create parts with complex geometries and locally tailored properties. For example, multiple homogenous material sources can be used in different regions of a print or be mixed during printing to define properties locally. Additionally, heterogeneous composites provide an opportunity for another level of tuning properties through processing. For example, within particulate-filled polymer matrix composites before curing, the presence of an applied electric and/or magnetic fields can reorient filler particles and form hierarchical structures depending on the fields applied. Control of particle organization is important because effective material properties are highly dependent on the distribution of filler material within composites once cured. While previous work in homogenization and effective medium theories have determined properties based upon ideal analytic distributions of particle orientations and spatial location, this work expands upon these methods generating discrete distributions from quasi-Monte Carlo simulations of the electromagnetic processing event. Results of simulations provide predicted microarchitectures from which effective properties are determined via computational homogenization.

    These particle dynamics simulations account for dielectric and magnetic forces and torques in addition to hydrodynamic forces and hard particle separation. As such, the distributions generated are processing field dependent. The effective properties for a composite represented by this distribution are determined via computational homogenization using finite element analysis (FEA). This provides a path from constituents, through processing parameters to effective material properties. In this work, we use these simulations in conjunction with a multi-objective optimization scheme to resolve the relationships between processing conditions and effective properties, to inform field-assisted additive manufacturing processes.

    The constituent set providing the largest range of properties can be found using optimization techniques applied to the aforementioned simulation framework. This key information provides a recipe for tailoring properties for additive manufacturing design and production. For example, our simulation results show that stiffness for a 10% filler volume fraction can increase by 34% when aligned by an electric field as compared to a randomly distributed composite. The stiffness of this aligned sample is also 29% higher in the direction of the alignment than perpendicular to it, which only differs by 5% from the random case [1]. Understanding this behavior and accurately predicting composite properties is key to producing field processed composites and prints. Material property predictions compare favorably to effective medium theory and experimentation with trends in elastic and magnetic effective properties demonstrating the same anisotropic behavior as a result of applied field processing. This work will address the high computational expense of physics simulation based objective functions by using efficient algorithms and data structures. We will present an optimization framework using nested gradient searches for micro barium hexaferrite particles in a PDMS matrix, optimizing on composite magnetization to determine the volume fraction of filler that will provide the largest range of properties by varying the applied electric and magnetic fields.

    more » « less