Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “supercritical fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Manufacturing and Hydro Testing of a 10 MWe sCO2 Axial Turbine

Abstract Supercritical CO2 (sCO2) power cycles designed for thermodynamic efficiencies above 50% push component designs into pressure and temperature combinations that approach the limits of American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code and ASME Piping Design Codes with Turbine inlet conditions near 250 bar and 715 °C. These design conditions which limit available materials to high nickel alloys are compounded by the high fluid densities of sCO2 at these conditions which result in very compact machinery at the 10 MWe class which are currently being developed for the DOE supercritical transformational electric power (STEP) program. Design of the STEP Turbine is complete with integrated lessons learned from a previous 1 MWth subscale demonstration of a high temperature 10 MWe Turbine designed under the DOE SunShot Program. Fabrication of the STEP Turbine is nearing completion with all housings and pressure containing components in hand. This paper provides a detailed look at design intentions and limitations for the pressurized sections of the main pressure casing and the fabrication of key turbine components. Additional discussion is provided regarding postfabrication inspection and acceptance of the casing for use after issues with casing welds were identified late in the fabrication process.

Engineering↗

Simulation of Methane Oxycombustion in Supercritical Carbon Dioxide

The Allam-Fetvedt Cycle (AFC) offers the potential for electricity generation with minimal or even negative carbon emissions by using oxycombustion in supercritical carbon dioxide (sCO2). Because of the extreme high temperatures and pressures encountered in these systems, simulation of oxycombustion requires use of a complex equation of state such as the modified Soave-Redlich-Kwong (SRK) equation, coupled with equally complex correlations for transport properties. The current work simulates a model combustor using both an ideal gas equation of state and the SRK equation of state. The model combustor features a jet of CH4 mixing with and reacting with O2 in the sCO2 ambient. Two jet configurations are considered: one of pure CH4, and one of CH4 mixed with sCO2. Comparison of results shows that key output properties such as combustor temperature and flow distribution are impacted by the choice of equations of state. Inspection of simulation results shows that mixing and reaction rates, as well as the uniformity of the flow, are impacted by the choice of equation of state. These results show that even though use of a complex equation of state greatly increases the computational costs of these simulations, it is critical to accurately capturing the physics of these systems.

computational fluid dynamics↗

Detailed Simulations of a Three-Stage Supercritical Carbon Dioxide Axial Compressor with a Focus on the Shrouded Stator Cavities Flow

This paper describes the findings of detailed simulations performed to investigate the impact of seal teeth cavity leakage flow on the aerodynamic and thermal performance of a three-stage supercritical CO 2 axial compressor. The study compares a shrouded stator configuration (with cavities) to a cantilevered stator configuration (without cavities) to highlight their differences. High-fidelity computational fluid dynamics simulations were performed using non-linear harmonic (NLH) and mixing plane assumptions, considering various possible rotor/stator interface configurations for mixing plane calculations. The key performance parameters for each case were compared, and the best-performing configuration selected for further analysis. The individual stage performance parameters are also examined and compared between the cantilevered and shrouded configurations. It was observed that in the shrouded case, the leakage flow enters the cavity downstream of the stator trailing edge and gets entrained into the primary flow upstream of the stator, leading to boundary layer changes at the hub and degradation of stator and downstream rotor performance. Vortical flow structures were also observed in the stator wells, which tended to change the flow angles around the region, thereby affecting mixing and velocity distribution, resulting in a slight deterioration of compressor performance. Additionally, the study examines windage heating due to shear work from rotating walls, including the seal teeth surface. The amount of shear work done on the leakage flow and the corresponding rise in fluid temperature were quantified, tabulated, and further compared with a simple analytical model, showing good agreement between them and, hence, validating the numerical approach used.

42 ENGINEERING↗

Multiphysics Simulation of Supercritical CO2 Gasification for Hydrogen Production

A bench-scale supercritical CO2 gasifier is being considered by the U.S. Department of Energy’s National Energy Technology Laboratory (NETL). A reacting computational fluid dynamics (CFD) model was developed using ANSYS/Fluent in order to investigate the impact of various operating parameters on syngas composition. The model was validated using available data in the literature. Simulations were conducted over a range of coal-slurry loadings, reaction temperatures, system pressures, and oxygen mass flowrates with the goal of optimizing hydrogen production for the considered reactor. The results from the simulations identified a positive correlation between slurry loading and hydrogen production, while reactor temperature and pressure had a limited impact on H2 production at the conditions of interest. Additionally, the simulations indicated that oxygen mass flowrate has less of an impact on hydrogen production at higher slurry loadings.

Prlina, Alexander↗

Assessing Compatibility of Natural Gas Pipeline Materials with Hydrogen, CO 2 , and Ammonia

Here, in this study, we examine the efficacy of repurposing natural gas (NG) pipelines for transporting hydrogen blends (with NG), ammonia, and CO 2 (gaseous and supercritical) from the standpoint of materials compatibility. Some information pertaining to component performance is also included, especially those components critical for pressurization and monitoring of flow. A listing of critical pipeline components and materials was developed, and their compatibilities was assessed for each fluid or gas type based on known compatibilities. Results indicate that pipeline materials should be suitable for gaseous CO 2 and anhydrous ammonia, but hydrogen blends greater than 12% may be problematic. Current compressor/regulator stations will not be suitable for use with either supercritical CO 2 or ammonia. Important knowledge gaps were identified, including (1) polymer performance with hydrogen/NG blends at low pressures, (2) compressor/regulator station polymers and epoxy coating materials with supercritical CO 2 , and (3) metal performances of hydrogen/NG blends at low pressures.

03 NATURAL GAS↗

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗

Multiphysics simulation of supercritical CO 2 gasification for hydrogen production

A bench-scale supercritical CO 2 gasifier is being considered by the U.S. Department of Energy’s National Energy Technology Laboratory (NETL). A reacting computational fluid dynamics (CFD) model was developed using ANSYS/Fluent in order to investigate the impact of various operating parameters on syngas composition. The model was validated using available data in the literature. Simulations were conducted over a range of coal-slurry loadings, reaction temperatures, system pressures, and oxygen mass flowrates with the goal of optimizing hydrogen production for the considered reactor. The results from the simulations identified a positive correlation between slurry loading and hydrogen production, while reactor temperature and pressure had a limited impact on H 2 production at the conditions of interest. Additionally, the simulations indicated that oxygen mass flowrate has less of an impact on hydrogen production at higher slurry loadings.

08 HYDROGEN↗

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗

Primordial black hole formation by vacuum bubbles. Part II

The discoveries of LIGO/Virgo black holes in recent years have revitalized the study of primordial black holes. In this work we investigate a mechanism where primordial black holes are formed by vacuum bubbles that randomly nucleate during inflation through quantum tunneling. After inflation, these bubbles typically run into the ambient radiation fluid with a large Lorentz factor. In our previous work, we assumed the bubble fields are strongly coupled to the standard model particles so that the bubble wall is impermeable. Here we complete this picture by considering bubbles interacting with the fluid only through gravity. By studying the scenario in several limits, we found that black holes could form in the either subcritical or supercritical regime. Depending on the model parameters, the resulting mass spectrum of the black holes could be wide or narrow, and may develop two peaks separated by a large mass range. Finally, with different spectra, these black holes may account for the LIGO/Virgo black holes, supermassive black holes, and may play an important role in dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of Density and Sound Speed in Mixtures Relevant to Supercritical CO2 Cycles

Abstract The objective of this research is to validate properties of mixtures relevant to supercritical carbon dioxide (sCO2) power cycles. Direct-fired sCO2 cycles are a promising technology for the future power generation systems. The working fluid of sCO2 cycles will be near and above critical point of CO2. One of the challenges is that the simulation of mixtures should consider real gas behavior. Expected operating conditions of Allam cycles reach up to 300 bar and 1000 °C. Characterizing the mixtures at the extreme conditions is an important issue in current researches and industrial applications. Thermophysical properties of mixtures may be beyond the valid range of the widely used database such as NIST REFPROP. Experimental data of mixture properties in the literature are limited which is necessary to develop high-fidelity design tools for sCO2 power cycles. We measured the density and sound speed of several multi-component mixtures. A temperature-controlled high-pressure test cell was used for the density measurements. Sound speed was measured by resonant frequency detection using an external speaker and a piezoelectric pressure sensor. Mixtures studied in this work include carbon dioxide, methane, oxygen, and water vapor. Properties of pure CO2 were measured to show the validity of our technique. Compositions were selected to be close to frozen mixtures at the inlet, mid-progress, and exhaust conditions of a model sCO2 combustor in the previous numerical simulation work. Corresponding reaction progress variables (RPV) were RPV = 0, 0.5, and 1. Temperature and pressure conditions of experiments are 310–450 K and 0–150 bar. In our study, density and sound speed from the NIST REFPROP database agree with experimental measurements within the range of our measurement uncertainties.

Energy & Fuels↗

ENVIRONMENTAL DEGRADATION OF HARD AND SOFT MAGNETIC MATERIALS IN GASEOUS AND SUPERCRITICAL CO2 ENVIRONMENTS

It may be beneficial to use hermetic designs for supercritical CO2 (sCO2) cycles machinery as they would eliminate CO2 leakage through shaft end seals. This would reduce CO2 emissions and operating costs for makeup of lost process fluid. Those designs may replace traditional oil-lubricated bearings with actively controlled magnetic bearings operating at high temperatures in the process fluid environment (to reduce the need for active cooling). This paper investigates the material degradation of various types of magnets in CO2. Included are several permanent and soft magnetic materials (Alnico 9C, Alnico 5-7C, and 18-T550 grade SmCo, and Hiperco 50) with or without coatings (nickel plating or C5 coating). The materials were exposed to: (1) flowing gaseous CO2 at 1,022 °F (550 °C) and atmospheric pressure in a furnace and (2) sCO2 at 842 °F (450 °C) and 1,500 psi (103 bar) in an autoclave. The preliminary mass change measured after total exposures of 1,000 hours and 2,000 hours are included.

supercritical CO2, magnetic material↗

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dissolution and reprecipitation of amorphous silica in silica Rich shales induces Non-Monotonic evolution of porosity in acidic reactive environments

Advances in sustainable subsurface energy technologies are crucial for meeting our energy and resource needs for a climate-resilient future. Novel strategies to harness subsurface shale reservoirs for recovering valuable metals and for enabling CO 2 storage are influenced by the morphological and mineralogical heterogeneities of these materials. In this context, delineating the interactions of highly acidic solutions such as wet supercritical CO 2 on shales with varying mineralogy is crucial to inform the stability of caprock seal for CO 2 storage and enhancements in permeability for fluid transport, reactivity, and storage. The feedback chemical effects associated with the interactions of acidic solutions on the morphologies and mineralogies of shales have not been extensively investigated. These insights are crucial for assessing temporal changes in the reactivity and the fate of the fluids in subsurface environments. In this study, we investigate the effect of 1M HCl solution on the chemistry and morphology of three different shale samples with varying carbonate, clay and silica contents. An increase in the amorphous content, from 37% to 41.3%, of silica-rich and carbonate/clay lean shale is noted due to reactions with an acidic solution which is attributed to the dissolution of Si-bearing phases such as clays, accompanied by SiO 2 precipitation. In shales bearing high content of clays and carbonates, significant increase in the pore volumes and surface areas are noted. Non-monotonic changes in the micron-scale porosity of silica rich–carbonate/clay lean (e.g., Mowry shale) are noted using in-situ X-ray microtomography experiments. Due to the initial mobilization of silica and dissolution of carbonate/clay phases, the total porosity slightly increases from 6.7% to 10.7% followed by a decrease to similar to 4% caused by SiO 2 reprecipitation. In conclusion, these findings suggest that even though silica is less reactive in acidic environments, the changes in the amorphous and crystalline content due to dissolution and reprecipitation alter the porosity and fluid flow paths.

04 OIL SHALES AND TAR SANDS↗

Numerical Simulations of Carbon Dioxide Storage Efficiency in Heterogeneous Reservoir Models

The U.S. Department of Energy’s National Energy Technology Laboratory (DOE-NETL) has been developing methods and tools (the online Carbon Dioxide Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN) tool) to estimate carbon dioxide (CO2) storage potential in subsurface reservoirs. The CO2 storage efficiency terms are input in the tool to calculate storage potential in targeted reservoirs. In this effort, two CO2 storage efficiency terms were evaluated: volumetric displacement ( E V ) and microscopic displacement ( E d ). The first term deals with efficiency of CO2 propagation into an accessible reservoir volume, while the second term evaluates effectiveness of native fluid displacement with CO2. The interpreted well logs and core sample measurements were applied to create the heterogeneous reservoir models including geostatistical realizations of porosity and intrinsic permeability fields. Supercritical CO2 was injected over the course of 30 years into brine-saturated reservoir models for clastics, limestone, and dolomite lithologies and deltaic fluvial, aeolian, shallow marine, and reef depositional environments by means of varying reservoir parameters and injection scenarios. The reservoir models providing vertically heterogeneous petrophysical properties and designated as “layered reservoir models” (with homogeneous parameters along each layer of the model) were not determined to be a transition between the homogeneous and heterogeneous models in respect to storage efficiency. Another finding shows that high-efficiency factors do not necessarily mean increased CO2 storage; they rather indicate that the available volume and pore space are more fully utilized. The CO2 storage efficiency factors were evaluated dynamically at the select time points using P 10 ‐ P 50 ‐ P 90 percentiles. The results of this study show that the P 10 ‐ P 90 distribution for volumetric efficiency is wider when compared to the microscopic efficiency. It was found that where dominant buoyancy forces drive the plume to the top of a target formation, the volumetric efficiency is low. Tighter sandstone and carbonate formations show prevalence of capillary forces and better utilization of reservoir volume.

Myshakin, Evgeniy M.↗

Permeability Evolution of Fractures in Shale in the Presence of Supercritical CO 2

We experimentally examined the impact of supercritical carbon dioxide (scCO 2 ) on fracture permeability and fracture surface characteristics in shales of various compositions. We measured permeability and fracture normal displacement at different effective stresses using both argon and scCO 2 as pore fluids. Both natural fractures and saw cuts in intact samples were used in our study. The permeability and fracture normal displacement decrease after multiple loading/unloading cycles, apparently as the result of inelastic compaction. In samples with high carbonate content, we observe an increase in permeability and fracture surface degradation after more than 3.5 days of exposure to scCO 2 , apparently due to carbonate dissolution. We demonstrate that the sensitivity of permeability to effective normal stress correlates well with the fracture normal displacement. The dependence of permeability on effective normal stress increases after exposure to scCO 2 as the fracture surface became more compliant.

58 GEOSCIENCES↗

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗

Supercritical Reforming of Wet Ethanol for High Efficiency Direct-injection Heavy-duty Compression-ignition Engines

The purpose of the research was to investigate the potential of using a mature bio-fuel in heavy-duty compression ignition engines. By co-optimizing both the fuel characteristics and engine system the potential for a superior outcome was demonstrated. The use of wet bio-ethanol eliminates the majority of the energy intensive distilling and dehydrating fuel production processes, which moves the fuel towards carbon neutral and also lowers the fuel costs. The relatively high water content of the resulting fuel is leveraged in the proposed novel combustion system by incorporating an integrated high efficiency exhaust waste heat recovery system. This results in significantly higher thermal efficiency. In addition the combustion system features low criteria pollutant emissions and the potential to reduce the initial cost of the engine system hardware. Substantial societal benefits are demonstrated through the co-optimization of the fuel and engine system. A computational proof-of-concept study has been performed to demonstrate the potential benefits of a novel wet ethanol heavy-duty compression ignition combustion system featuring integrated exhaust waste heat recovery. A combined in-cylinder closed cycle 3D computational fluid dynamics (CFD) - 1D engine system simulation approach was used. The models were validated to baseline engine data using diesel fuel and then applied to the wet ethanol study. The original concept was to maximize exhaust waste heat recovery through supercritical reforming of the wet ethanol fuel. Phase I simulation results indicated that the optimal solution for maximum engine efficiency gains were realized through maximizing thermo-mechanical recuperation with negligible fuel reformation. The results show the potential to achieve impressive gains in brake thermal efficiency (BTE) over the base diesel engine. The potential to increase BTE up to 20.9% over the base diesel engine was demonstrated, with even larger gains possible through reduced in-cylinder heat transfer losses. The majority of the efficiency gains were realized through integrated high efficiency exhaust waste heat recovery. The concept also has the potential to achieve future ultra-low NO x emissions standards and negligible engine-out soot emissions. The mixing controlled combustion of high temperature wet ethanol features relatively low engine-out NO x emissions without the need for exhaust gas recirculation (EGR). The soot free combustion enabled by the relatively high oxygen content of ethanol also allows for the use of stoichiometric mixing controlled combustion, which is not practical with diesel fuel. When stoichiometric combustion is used a simple passive 3-way catalyst can be used for exhaust emissions after-treatment and near zero tailpipe emissions. The Phase I simulations results have defined the system layout and requirements in preparation for the Phase II experimental proof-of-concept study. The potential applications of the research include most current applications of diesel engines. The Phase I study focused on heavy-duty on-highway class 8 trucks. However, virtually any application that requires highly efficient clean power generation would benefit from the novel engine system proposed. The results indicate substantial fuel cost savings and reduced greenhouse gas emissions with similar or reduced initial system hardware costs compared to modern diesel engine systems.

09 BIOMASS FUELS↗

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗