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: Surface tension and evaporation behavior of liquid fuel droplets at transcritical conditions: Towards bridging the gap between molecular dynamics and continuum simulations
The phase transition from subcritical to supercritical conditions, referred to as transcritical behavior, significantly impacts the evaporation and fuel–air mixing in high-pressure liquid-fuel propulsion systems. Transcritical behavior is characterized as a transition from classical two-phase evaporation to a single-phase gas-like diffusion regime as surface tension and latent heat of vaporization reduce. However, the interfacial behavior represented by the surface tension coefficient and evaporation rate during this transition which are crucial inputs for Computational Fluid Dynamics (CFD) simulations of practical transcritical fuel spray is still missing. This study aims at developing new evaporation rate and surface tension models for transcritical n-dodecane droplets using molecular dynamics (MD) simulations irrespective of the droplet size. As MD simulations are primarily limited to the nanoscale, the new models can bridge the gap between MD and continuum simulations and enable the direct application of these findings to microscopic droplets. A new characteristic timescale, i.e., “undroplet time,” is defined which marks the transition from classical two-phase evaporation to single-phase gas-like diffusion behavior. The undroplet time indicates the onset of droplet core disintegration and penetration of nitrogen molecules into the droplet, which occurs after the vanishment of the surface tension. By normalizing the time with respect to the undroplet time, the rate of surface tension decay, evaporation rate, and the rate of droplet mass depletion become independent of the droplet size. Calculation of pairwise correlation coefficients for the entire MD results shows that both surface tension coefficient and evaporation rate are strongly correlated with the background temperature, while pressure and droplet size play a less significant role past the critical point. Therefore, new models for surface tension coefficient and evaporation rate spanning from sub- to supercritical conditions are developed as a function of background pressure and temperature, which can be used in continuum simulations. The identified phase change behavior based on the undroplet time shows a good agreement with the phase change regime maps obtained using microscale experiments and nanoscale MD predictions.  more » « less
Award ID(s):
2237124
PAR ID:
10474222
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Fuel
ISSN:
0016-2361
Page Range / eLocation ID:
130187
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In the pursuit of enhanced engine performance and reduced emissions, the design of liquid-fueled propulsion systems is shifting towards much higher combustor pressures, surpassing the nominal critical pressure of the fuel and air. This trend leads to the adoption of supercritical conditions, wherein the liquid fuel is injected into the ambient air at supercritical pressure and temperature, causing the fuel temperature to exceed its nominal critical point. This transition from a liquid-like to a gas-like behavior, known as "transcritical behavior," is a crucial aspect governing the operation of modern high-pressure propulsion and energy conversion systems. In these systems, the primary liquid jet breakup and the subsequent break-up of the resulting droplets into smaller droplets, namely secondary breakup, significantly impact mixing and combustion processes. Despite its importance, there has been a limited focus on droplet breakup at supercritical conditions, particularly at higher flow speeds relevant to high-speed liquid-fuel propulsion systems. Surface tension effects are often neglected in the simulation of transcritical flow, assuming surface tension vanishes beyond the critical point, while recent experiments and molecular dynamics simulations suggest that surface tension effects persist at transcritical conditions. To gain insight into the effects of surface tension on transcritical flows, we have developed a fully compressible multiphaseDirect Numerical Simulation (DNS) approach that accounts for decaying surface effects. The diffuse interface method is employed to represent transcritical interfaces, accounting for surface tension effects calculated using molecular dynamics simulations. This approach is employed to investigate the behavior of subcritical n-dodecane droplets in a supercritical nitrogen environment interacting with a shockwave, aiming to identify the governing breakup regimes at transcritical conditions. The development of quantitative measures enables the generalization of droplet breakup modes for transcritical droplets. The insights gained from this study contribute to advancing the understanding of transcritical liquid breakup, providing valuable knowledge for designing and optimizing high-speed propulsion systems 
    more » « less
  2. Fuels when sprayed under superheated and elevated fuel pressure show different behavior than traditional fuel injection sprays. In this work optical diagnostics were used to study the behavior of Jet A-1 under subcritical, transcritical, and supercritical sprays into open air ambience. Five different temperatures were tested, and the resultant spray images were processed to obtain quantitative measurements such as spray penetrations, and spray cone angle for each case. The spray structure transition with changing parameters from subcritical, transcritical, and supercritical states were also studied. The transition between the three different states are shown in this study and the resulting spray cone angles and penetrations are compared for the fuel. The results show that a transcritical spray has a measurable variation in the spray cone formation and penetration process for a fixed injection pressure. At this state the spray cone shows a bimodal spray angle relationship with increasing penetration. Flash boiling of the fuel is observed near the nozzle of the injector. Increasing the temperature further into the supercritical regime, the spray plume shows a thinning of the jet near the nozzle with a reduced overall penetration compared to lower temperatures. 
    more » « less
  3. 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
  4. This paper numerically investigates the evaporation characteristics of a single n-dodecane fuel droplet in high-pressure nitrogen environment relevant to rotating detonation engines. A validated computational fluid dynamics solver coupled with real-fluid thermophysical models is utilized. The effects of pressure, droplet temperature, and ambient gas temperature on the evaporation rate are analyzed by tracking the droplet diameter evolution. Two interface tracking techniques, namely a mean density-based method and a novel vapor-liquid equilibrium-based method, are implemented and compared. The results show appreciable deviations from the classical d2-law for droplet evaporation. Increasing the ambient temperature and droplet temperature (toward critical point) substantially accelerate the evaporation process. Meanwhile, higher pressures decrease the evaporation rate owing to slower species/thermal diffusions. At certain conditions, discernible differences are observed between the two interface tracking methods indicating deficiencies in the simple mean density approach. The paper demonstrates an effective computational framework for transcritical droplet evaporation simulations. And the generated high-pressure droplet evaporation datasets can inform sub-model development for spray combustion modeling. 
    more » « less
  5. 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