Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Validation of Octane Hyperboosting Phenomenon in Prenol and Structurally Related Olefinic Alcohols

Hyperboosting is a recently discovered phenomenon in which the research octane number (RON) of a blend is higher than both the neat blendstock and the neat fuel it was blended into. RON is a measure of a fuel's resistance to knock, and knock is a cause of engine inefficiency. Blends which exhibit hyperboosting are blends in which an overall improvement in engine efficiency may be expected. The first discovery of hyperboosting came from blending experiments in which prenol was blended into several different base fuels. Here, ignition delay time (IDT) measurements taken using a commercially available constant volume combustion chamber called the Advanced Fuel Ignition Delay Analyzer (AFIDA) are presented. The data show that some prenol blends have longer IDTs (lower reactivity) than either neat prenol or the base fuel, providing further evidence of hyperboosting. Additionally, more blending data is presented in which the base fuel is varied, which allows for a better understanding of hyperboosting sensitivity to chemical classes. The data indicate that aromatics may inhibit, and branched alkanes may enhance the magnitude of hyperboosting observed. Enthalpy of vaporization estimates are also given for several molecules which are blended into a 4-component surrogate. These estimates are derived from Equation of State simulations and reveal that there is no observable correlation between hyperboosting and enthalpy of vaporization. Blending data for molecules which share structural similarities with prenol are also presented. Structure property relationships are suggested, in which the double bond motif of prenol seems to play an important role in hyperboosting. This may help to understand hyperboosting and its underlying mechanism. Lastly, blending curves of surrogate blends with prenol experienced hyperboosting under lean (Homogeneous Charged Compression Ignition-HCCI) operating conditions, which validates that hyperboosting is not an artifact of the octane test methods, but inherent to the properties of prenol.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Methane mass transfer in mesoporous silica saturated with liquid hydrocarbons

Mass transfer across gas/liquid interfaces plays a central role in many industrial applications. In particular, gas dissolution and diffusion in liquid hydrocarbon mixtures, confined in nanometer-sized pores, is an essential mechanism during enhanced oil recovery (EOR) from unconventional formations. In this work, we have measured methane (C 1 ) diffusion in n-decane (C 10 ), n-hexadecane (C 16 ), and mixtures of C 10 + C 16 in a mesoporous material with an average pore size of 4 nm at 50 °C and ~ 8 MPa of gas pressure. A key conclusion of this work is that the diffusivities, measured in the bulk phase, for the relevant binary systems, are sufficient to predict the diffusion behavior of the corresponding multicomponent systems in the porous medium by using Wilke’s equation for the evaluation of effective component diffusivities, combined with an accurate equation of state (EOS) representation of the (bulk) phase behavior. This observation can facilitate accurate prediction of recovery processes in unconventional formations and, potentially, guide other applications that entail gas-liquid interface mass transfer in mesoporous materials.

42 ENGINEERING↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2021

Comparison of ablators for the polar direct drive exploding pusher platform

In this work, we examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. The higher tensile strength of these materials compared to plastic enables a thinner ablator to support higher gas pressures, which could help optimize its performance for plasma physics experiments, while ablators containing boron enable the possibility of collecting additional data to constrain models of the platform. Applying recently developed and experimentally validated equation of state models for the boron materials, we examine the performance of these materials as ablators in 2D simulations, with particular focus on changes to the ablator and gas areal density, as well as the predicted symmetry of the inherently 2D implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The STAG code: A fully relativistic super transition array calculation using Green’s functions

Calculating opacities for a wide range of plasma conditions (i.e. temperature, density, element) requires detailed knowledge of the plasma configuration space and electronic structure. For plasmas composed of heavier elements, relativistic effects are important in both the electronic structure and the details of opacity spectra. We extend our previously described superconfiguration and super transition array capabilities (Gill et al., 2023) to include a fully relativistic formalism. Further, the use of hybrid bound-continuum supershells in our superconfigurations demonstrates the importance of a consistent treatment of bound and continuum electrons in dense plasma opacities, and we expand the discussion of these consequences to include issues associated with equation of state and electron correlations between bound and continuum electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Approximation of refrigerant thermophysical properties using neural networks to speed up transient thermofluid simulations

Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.

Ma, Jiacheng↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Using X-ray spectrum imaging to quantify fission gas in high-burnup UO2 fuel

Analysis of fission gas bubbles (FGBs) in light-water reactor (LWR) fuel is needed to improve the understanding and predictive capability of fuel evolution under normal- and transient conditions. One of the most important parameters of a FGB is the pressure of the gas, primarily Xe, inside it. Bubble volume and location (inter- vs. intragranular) are important considerations as well. However, experimental analysis of such bubbles is challenging due to their small size and embedded nature, and usefulness of the data requires large numbers be analyzed. This paper proposes a method to measure the pressure of Xe bubbles using X-ray spectrum imaging (XSI) in scanning transmission electron microscopy (STEM). From X-ray generation and instrumental parameters, the number of Xe atoms yielding a given number of Xe L-series X-ray counts is estimated, and then from the estimated volume of the bubble the Xe density is estimated, which can then be converted to pressure via an equation-of-state. We apply the method to XSI datasets from high burnup (HBu) fuel from North Anna 1 reactor in as-irradiated and post-LOCA-test conditions and find Xe pressures in the FGBs clustered around 1 GPa.

Parish, Chad [ORNL] (ORCID:0000000312097439)↗
DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗
DOE OSTI2020

Effects of dissipative baryon current in heavy-ion collisions at RHIC-BES energies

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and p T spectra in heavy-ion collisions at beam energy scan energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A machine learning study on spinodal clumping in heavy ion collisions

Possible observables of baryon number clustering due to the instabilities occurring at a first order QCD phase transition are discussed. The dynamical formation of baryon clusters at a QCD phase transition can be described by numerical fluid dynamics, augmented with a gradient term and an equation of state with a mechanically unstable region. It is shown that the dynamical description of this phase transition, in nuclear collisions, will lead to the formation of dense baryon clusters at the phase boundary. State-of-the-art machine learning methods find that the coordinate space clumping leaves characteristic imprints on the spatial net density distribution in almost every event. On the other hand the momentum distributions do not show any clear event-by-event features. Lastly, it is shown that the 'third order' cumulant, the skewness, shows a peak at the beam energy where the system, created in the heavy ion collision, reaches the deconfinement phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Lastly, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Precise measurement of nuclear interaction cross sections towards neutron-skin determination with R 3 B

The R 3 B (Reactions with Relativistic Radioactive Beams) experiment as a major instrument of the NUSTAR collaboration for the research facility FAIR in Darmstadt is designed for kinematically complete studies of reactions with high-energy radioactive beams. Part of the broad physics program of R 3 B is to constrain the asymmetry term in the nuclear equation-of-state and hence improve the description of highly asymmetric nuclear matter (e.g., in neutron stars). For a precise determination of the neutron-skin thickness – an observable which is directly correlated with the symmetry energy in theoretical calculations – by measuring absolute fragmentation cross sections, it is essential to quantify the uncertainty and challenge the reaction model under stable conditions. During the successful FAIR Phase-0 campaign of R 3 B, we precisely measured the energy dependence of total interaction cross sections in 12 C + 12 C collisions, for a direct comparison with calculations based on the eikonal reaction theory.

Ponnath, L. (ORCID:0000000286742624)↗
DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗