skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Thermodynamic analysis of anomalous region, critical point, and transition from subcritical to supercritical states: Application to van der Waals and five real fluids
All fluids exhibit large property-variations near the critical point in a region identified as the anomalous state. The anomaly starts in the liquid and extends well into the supercritical state, which can be identified thermodynamically using the Gibbs free energy (g). The specific heat, isobaric expansion, and isothermal compressibility parameters governing the transitions are: (cp/T), (vβ), and (vκ), rather cp, β, and κ. They are essentially the second-order derivatives of g and have two extrema (minimum, maximum); only maxima reported ever. When applied to the van der Waals fluid, these extrema exhibit closed loops on the phase-diagram to satisfy d3g = 0 and map the anomalous region. The predicted liquid-like to gas-like transitions are related to the ridges reported earlier, and the Widom delta falls between these loops. Evidently, in the anomalous region, both the liquid and the supercritical fluid need to be treated differently. Beyond the anomalous states, the supercritical fluids show monotonic, gradual changes in their properties. The analysis for argon, methane, nitrogen, carbon dioxide, and water validates the thermodynamic model, supports the stated observations, and identifies their delimiting pressures and temperatures for the anomalous states. It also demonstrates the applicability of the law of corresponding states. Notably, the critical point is a state where d3g = 0, the anomaly in the fluid's properties/behavior is maximal, and the governing parameters approach infinity. Also the following are presented: (a) the trajectory of the liquid–vapor line toward the melt-solid boundary and (b) a modified phase diagram (for water) exhibiting the anomalous region.  more » « less
Award ID(s):
2327571 2231393
PAR ID:
10538000
Author(s) / Creator(s):
; ;
Publisher / Repository:
AIN Publishing
Date Published:
Journal Name:
Physics of Fluids
Volume:
36
Issue:
2
ISSN:
1070-6631
Page Range / eLocation ID:
026105-1-30
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A unified one-dimensional (1D), steady-state flow and heat transfer model is presented for the pipeline transport of fluids at high pressures, including the supercritical (SC) conditions. The model includes a generalized temperature equation, presented here for the first time, and accounts for all of the important effects, including the property variation, viscous dissipation, Joule-Thomson (J-T) cooling, and heat exchange with the surrounding. With appropriate approximations, this model can yield all isothermal and nonisothermal pipe flow solutions reported thus far. A generalized multizone integral method is developed which solves the two resulting algebraic equations for pressure and temperature in conjunction with a property database, such as the National Institute of Standard and Technology (NIST) reference fluid thermodynamic and transport properties (REFPROP). With appropriately selected number and size of the zones and using property values at the mean temperature and pressure within each zone, this integral method can accurately predict the complex effects of the governing parameters, such as the pipe diameter and length, inlet and exit pressures, mass flowrate, J-T cooling, and inlet and surrounding temperatures. Its accuracy for small-to-large diameter pipes has been ascertained by a comparison with the numerical solutions of the differential form of governing equations that requires a large number of small grids along the pipe and the values of mean properties within each grid. Indeed, this integral model can be used for the pipeline transport at both subcritical and supercritical pressures as long as the fluid does not encounter its anomalous states and the phase-change. 
    more » « less
  2. As fuel injection systems advance towards higher injection pressures and the combustor environment increases in both temperature and pressure in the pursuit of improved emissions and efficiency, advanced combustion strategies are required. Injecting fuel as a supercritical fluid has the potential to improve fuel/air mixing and eliminate steps in the spray vaporization process. Experiments are carried out on a heated fuel injector in an open-air test cell using Mie scattering, Schlieren imaging, and long-distance microscopy to investigate changes in spray characteristics with varying temperature and pressure. Spray angle, spray penetration length, and vapor-liquid ratio data are collected and evaluated. Near-nozzle imaging shows distinct changes in spray morphology during the initial microseconds of spray formation. Sprays injected under conditions further into the supercritical regime exhibit increased spray angle and vapor-to-liquid ratio. Spray penetration is found to decrease with increasing temperature. A jump in vapor-to-liquid ratio is observed in the vicinity of 568 K, indicating a transition in spray behaviour trending towards more rapid fuel/air mixing across the transcritical region. Changes in the micro-scale structure of the spray during the initial microseconds of spray formation exhibit this same narrow transition region. A significantly greater fraction of the spray plume is observed to be in a vapor or vapor-like state at a given time after injection initiation as the injection conditions are advanced into the supercritical state. These findings indicate that injection fuel as a supercritical fluid has the potential to improve the mixing of a fuel/air charge, and thus, improve combustion quality. 
    more » « less
  3. Supercritical fluids, often present in modern high-performance propulsion systems, result from elevated operating pressures. When these systems utilize fluid mixtures as fuel or oxidizers, a transcritical effect often occurs. This effect can lead to misjudgments, as mixture critical points exceed those of individual components. Fluid mixing may induce phase separation, creating liquid and vapor phases due to the transcritical multi-component effect. Consequently, two-phase modeling is essential for transcritical and supercritical fluids. Traditional interface capturing methods, like Volume of Fluid (VOF) and Level Set (LS), present challenges such as computational expense and lack of conservatism. The Phase Field (PF) method, or the Diffuse Interface (DI) method which uses a phase fraction transport equation, emerges as a conservative alternative. Despite the absence of an initial interface in transcritical fluids, phase separation from mixing may form liquid droplets, necessitating multiphase modeling. To address these complexities, a Vapor-Liquid Equilibrium (VLE) model, coupled with the PR equation of state, is introduced. This model estimates phase fractions, liquid and vapor compositions, densities, and enthalpies through a flash problem solution. The conventional PF model is enhanced by replacing the phase fraction transport equation with VLE-derived values. The resulting VLE-based PF method is implemented into an OpenFOAM compressible solver, ensuring numerical stability with explicit phase field terms and a new CFL criterion. Test cases involve 1D interface convection and 2D droplet convection. In the 1D test, the VLE-based PF model adeptly captures interfaces, adjusting thickness as needed. The 2D droplet case, challenging due to a non-aligned Cartesian grid, exhibits uniform interface thickness and preserves droplet shape. The VLE-based PF model demonstrates versatility and reliability in capturing complex fluid behaviors, offering promising prospects for future research. 
    more » « less
  4. null (Ed.)
    Water + elastin-like polypeptides (ELPs) exhibit a transition temperature below which the chains transform from collapsed to expanded states, reminiscent of the cold denaturation of proteins. This conformational change coincides with liquid–liquid phase separation. A statistical-thermodynamics theory is used to model the fluid-phase behavior of ELPs in aqueous solution and to extrapolate the behavior at ambient conditions over a range of pressures. At low pressures, closed-loop liquid–liquid equilibrium phase behavior is found, which is consistent with that of other hydrogen-bonding solvent + polymer mixtures. At pressures evocative of deep-sea conditions, liquid–liquid immiscibility bounded by two lower critical solution temperatures (LCSTs) is predicted. As pressure is increased further, the system exhibits two separate regions of closed-loop of liquid–liquid equilibrium (LLE). The observation of bimodal LCSTs and two re-entrant LLE regions herald a new type of binary global phase diagram: Type XII. At high-ELP concentrations the predicted phase diagram resembles a protein pressure denaturation diagram; possible “molten-globule”-like states are observed at low concentration. 
    more » « less
  5. The authors recently reported that undercooled liquid Ag and Ag–Cu alloys both exhibit a first order phase transition from the homogeneous liquid (L-phase) to a heterogeneous solid-like G-phase under isothermal evolution. Here, we report a similar L–G transition and heterogenous G-phase in simulations of liquid Cu–Zr bulk glass. The thermodynamic description and kinetic features (viscosity) of the L-G-phase transition in Cu–Zr simulations suggest it corresponds to experimentally reported liquid–liquid phase transitions in Vitreloy 1 (Vit1) and other Cu–Zr-bearing bulk glass forming alloys. The Cu–Zr G-phase has icosahedrally ordered cores versus fcc/hcp core structures in Ag and Ag–Cu with a notably smaller heterogeneity length scale Λ . We propose the L–G transition is a phenomenon in metallic liquids associated with the emergence of elastic rigidity. The heterogeneous core–shell nano-composite structure likely results from accommodating strain mismatch of stiff core regions by more compliant intervening liquid-like medium. 
    more » « less