An official website of the United States government
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
Computational Analysis of Secondary Droplet Breakup at Transcritical Conditions with Surface Tension Effects
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
- Award ID(s):
- 2237124
- PAR ID:
- 10555916
- Publisher / Repository:
- American Institute of Aeronautics and Astronautics
- Date Published:
- ISBN:
- 978-1-62410-711-5
- Format(s):
- Medium: X
- Location:
- Orlando, FL
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
This study investigates the atomization process in Respimat® Soft MistTM Inhalers (SMIs) using a validated Volume of Fluid (VOF)-to-Discrete Phase Model (DPM) to simulate the transition from colliding liquid jets to aerosolized droplets. Key parameters, including colliding jet inlet velocity, surface tension, and liquid viscosity, were systematically varied to analyze their impact on the atomization, i.e., aerosolized droplet size distributions. The VOF-to-DPM simulation results indicate that higher jet inlet velocities enhance ligament fragmentation, producing finer and more uniform droplets while reducing total atomized droplet mass. The relationship between surface tension and atomization performance in colliding jet atomization is not monotonic. Reducing surface tension plays a complex dual role in the atomization process. On the one hand, lower surface tension enhances the likelihood of liquid jet breakup into a liquid sheet, leading to the formation of smaller ligaments under the same airflow conditions and shear forces. This increases the probability of generating more secondary droplets. On the other hand, reduced surface tension also destabilizes the liquid surface shape, decreasing the formation of fine, high-sphericity droplets in regimes where surface tension is a dominant force. Viscosity also influences atomization through complex mechanisms, i.e., lower viscosity reduces resistance to ligament breakup but promotes droplet interactions and coalescence, while higher viscosity suppresses ligament fragmentation, generating larger droplets and reducing atomization efficiency. The validated VOF-to-DPM framework provides critical insights for enhancing the performance and efficiency of inhalation therapies. Future work will incorporate nozzle geometry, jet impingement angles, and surfactant effects to better understand and optimize the atomization process in SMIs, focusing on achieving preferred droplet size distributions and emitted doses for enhanced drug delivery efficiency in human respiratory systems.more » « less
-
Vapor-liquid equilibrium (VLE) is a family of first-principled thermodynamic models for transcritical multiphase flows, which can accurately capture the phase transitions at high-pressure conditions that are difficult to deal with using other models. However, VLE-based computational fluid dynamics (CFD) simulation is computationally very expensive for multi-component systems, which severely limits its applications to real-world systems. In this work, we developed a new ISAT-VLE method based on the in situ adaptive tabulation (ISAT) method to improve the computational efficiency of VLE-based CFD simulation with reduced memory usage. We developed several ISAT-VLE solvers for both fully conservative (FC) and double flux (DF) schemes. New methods are proposed to delete redundant records in the ISAT-VLE table and the ISAT-VLE method performance is further improved. To improve the convergence of the VLE solvers, a modified initial guess for equilibrium constant is also introduced. Simulations of high-pressure transcritical two-phase temporal mixing layers and shock-droplet interaction were conducted using the ISAT-VLE CFD solvers. The simulation results show that the new method obtains a speed-up factor approximately from 10 to 60 and the ISAT errors can be controlled within 1%. The shock-droplet interaction results show that the DF scheme can achieve a higher speed-up factor than the FC scheme. The two sets of simulations exhibit the phase separation at high-pressure conditions. It was found that even at supercritical pressures with respect to each component, the droplet surface could still be in a subcritical two-phase state, because the mixture critical pressure is often significantly higher than each component and hence triggers phase separation. In addition, a shock wave could partially or completely convert the droplet surface from a subcritical two-phase state to a single-phase state by raising temperature and pressure.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.2514/6.2024-2380
