Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SUPERSONIC”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗
DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗
DOE OSTI2026

From dynamics to kinetics for the dissociation of •QOOH radicals derived from isopentane

Thermal dissociation rates for hydroperoxyalkyl (•QOOH) radicals are challenging to measure due to their inherent instability. Recent dynamics experiments have demonstrated that in situ synthesis of •QOOH radicals in a supersonic expansion in combination with IR action spectroscopy provides an effective route to accurately map their energy-resolved dissociation rates. Collaborative theoretical analyses provide a route for converting from those microcanonical rates to accurate temperature- and pressure-dependent thermal rate constants. Here, we present the results of this conversion for the dissociation of two isopentyl β-QOOH isomers of relevance to isopentane oxidation. The results are compared with existing literature values, as well as with results obtained in a similar fashion for related •QOOH systems. Extrapolations from benchmark calculations of barrier properties for smaller systems to accurate predictions for larger systems and multidimensional hindered rotor treatments are both shown to be important components of the theoretical analysis. For completeness, we extend the theoretical analysis to a treatment of the full RO 2 /QOOH system for isopentyl (1,1-dimethyl-propyl) radical.

Ab initio kinetics↗
DOE OSTI2025

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗
DOE OSTI2024

Effect of insulator end cap thickness on time-dependent Hartmann flow in a rotating mirror

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Identification of the Elusive Methyl-Loss Channel in the Crossed Molecular Beam Study of Gas-Phase Reaction of Dicarbon Molecules (C 2 ; X 1 Σ g + /a 3 Π u ) with 2-Methyl-1,3-butadiene (C 5 H 8 ; X 1 A')

The crossed molecular beams technique was utilized to explore the reaction of dicarbon C 2 (X 1 Σ g + /a 3 Π u ) with 2-methyl-1,3-butadiene (isoprene, CH 2 C(CH 3 )CHCH 2 ; X 1 A') at a collision energy of 28 ± 1 kJ mol⁻¹ using a supersonic dicarbon beam generated via photolysis (248 nm) of helium-seeded tetrachloroethylene (C 2 Cl 4 ). Here, experimental data combined with previous ab initio calculations provide evidence of the detection of the hitherto elusive methyl elimination channels leading to acyclic resonantly stabilized hexatetraenyl radicals: 1,2,4,5-hexatetraen-3-yl (CH 2 CC•CHCCH 2 ) and/or 1,3,4,5-hexatetraen-3-yl (CH 2 CHC•CCCH 2 ). These pathways are exclusive to the singlet potential energy surface, with the reaction initiated by the barrierless addition of dicarbon to one of the carbon-carbon double bonds in the diene. In combustion systems, both hexatetraenyl radicals can isomerize to the phenyl radical (C 6 H 5 ) through a hydrogen atom assisted isomerization – the crucial reaction intermediate and molecular mass growth species step toward the formation of polycyclic aromatic hydrocarbons (PAHs) and soot.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Photodissociation and Infrared Spectroscopy of U + (CO 2 ) n , UO + (CO 2 ) n , and UO 2+ (CO 2 ) n Cation-Molecular Complexes

Laser vaporization of uranium in a pulsed supersonic expansion of carbon dioxide is used to produce complexes of the form U + (CO 2 ) n , UO + (CO 2 ) n , and UO 2+ (CO 2 ) n . These ions are selected in a reflectron time-offlight mass spectrometer and studied with visible laser photodissociation and tunable infrared laser photodissociation spectroscopy in the region of the CO 2 antisymmetric stretch. The dissociation patterns and spectroscopy of these ions indicate that CO 2 ligands are intact molecules. Although reaction products that form oxide-carbonyl or oxalate species are predicted to be stable, there is no direct evidence in the frequency range studied for the formation of these species. There is no clear indication for the coordination numbers for singly charged uranium and its oxide complexes with CO 2 . However, there is strong support in the vibrational patterns for an eight-coordinate complex of the doubly charged UO 2+ species, i.e., UO 2+ (CO 2 ) 8 .

Cluster chemistry↗
DOE OSTI2024

Thermally Initiated Formation of Criegee Intermediate CH 2 OO in the Oxidation of Ethane

Criegee intermediates (CIs) play an important role in atmospheric chemistry as a transient source of the OH radical through their formation by the ozonolysis of unsaturated organic compounds. Here, we report thermally initiated formation of the smallest CI (CH 2 OO) in the oxidation of ethane (CH 3 CH 3 ) that may be relevant to combustion and flames. The SiO 2 /SiC oxidation microreactor is heated to 1800 K and has a short residence time of ∼100 μs. The CH 2 OO we observe is likely formed in a lower-temperature region near the microreactor’s exit. Plausible mechanisms for CH 2 OO formation and retention under these conditions mediated by methylperoxy (CH 3 OO) radicals are discussed. Pure rotational spectra of CH 2 OO and other intermediates (HO 2 , CH 3 CHO, CH 2 CHOH, c-CH 2 OCH 2 , CH 3 CH 2 CHO, CH 3 OOH, and HCOOH) are detected with a chirped-pulse Fourier transform millimeter-wave spectrometer operating in the frequency range of 60–90 GHz. In conclusion, detection occurs in a molecular beam, where the species are supersonically cooled to 5 K.

Hydrocarbons↗
DOE OSTI2025

Cation-π Bonding in Actinides: UO x + (Benzene) ( x = 0, 1, 2) Complexes Studied with Threshold Photodissociation Spectroscopy and Theory

Cation-π complexes of the form UO x + (benzene) (x = 0, 1, 2) are produced by laser vaporization and cooled in a supersonic molecular beam. These ions are mass selected and studied with UV–visible laser photodissociation spectroscopy. Each of these complexes photodissociates by elimination of the benzene ligand. Above an energetic threshold, the absorption and photodissociation are continuous, indicating a high density of strongly coupled electronic states. The thresholds for the dissociation of each of these three complexes are measured and assigned as their respective bond dissociation energies. The bond energies determined [U + –(benzene): 42.5 ± 0.3 kcal/mol; UO + –(benzene): 41.0 ± 0.3 kcal/mol; UO 2 + –(benzene): 39.7 ± 0.3 kcal/mol] are comparable to those of transition metal ion-benzene complexes. Computational studies at the DFT/B3LYP level complement the experiments, predicting dissociation energies in reasonably good agreement with the experiments. Experiments and theory agree that the U+(benzene) complex is more strongly bound than its corresponding oxide ions. This new thermochemistry on actinide cation-π bonding should stimulate higher-level computational studies on these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Bond Dissociation Energy, Ionization Energy, and Electronic Structure of Thorium Dimer

Diatomic thorium, Th 2 , has been investigated using a laser ablation, supersonic expansion source to produce the molecule and resonant two-photon ionization spectroscopy to measure its bond dissociation energy (BDE) and ionization energy (IE). The molecule has a high density of states in the vicinity of its bond dissociation energy, leading to rapid predissociation as soon as this energy is exceeded. The BDE is identified from this predissociation threshold as D 0 (Th 2 ) = 2.857(7) eV, where the assigned error limit is provided in parentheses in units of the last quoted digit. Similarly, the one-photon ionization threshold has been measured, providing the ionization energy IE(Th 2 ) = 5.042(4) eV. Together with a thermochemical cycle and the atomic ionization energy, these values provide the BDE of the cation, giving D 0 (Th 2 + ) = 4.122(8) eV. Computations show that Th 2 has three nearly degenerate low-lying electronic states (1 3 Σ u + , 1 1 Σ g + , and 1 3 Δ g ) with bonding dominated by 7s and 6d orbitals, indicating predominantly transition-metal-like behavior. The 1 3 Σ u + state exhibits a triple bond, whereas the 1 1 Σ g + and 1 3 Δ g states possess quadruple-bond character and correspondingly shorter bonds. Although 1 3 Σ u + is predicted to be the lowest state without spin–orbit coupling, the large spin–orbit stabilization of the 1 3 Δ g state makes its Ω = 1 g component the ground state. Furthermore, the calculated dissociation energy (2.840 eV) and ionization energy of Th 2 (5.098 eV) are both in excellent agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Predissociation-based measurements of bond dissociation energies: US 2 , OUS, and USe

The uranium-containing molecules US 2 , OUS, and USe have been investigated using a pulsed laser ablation supersonic beam molecular source with time-of-flight mass spectrometric detection. Spectra have been recorded using the resonant two-photon ionization method over the spectroscopic range from 277 to 238 nm. These species have a myriad of excited electronic states in this spectroscopic region, leading to spectra that are highly congested and appear quasicontinuous. Sharp predissociation thresholds are observed, allowing precise bond dissociation energies to be measured. In the case of the triatomic molecules, it was necessary to use one laser for excitation and a delayed laser for ionization in order to observe a sharp predissociation threshold that allowed a precise bond dissociation energy to be measured. The resulting thermochemical values are D 0 (SU-S) = 4.910 ± 0.003 eV, D 0 (OU-S) = 5.035 ± 0.004 eV, and D 0 (USe) = 4.609 ± 0.009 eV. These results provide the first measurement of D 0 (USe) and reduce the error limits in the previous values of D 0 (SU-S) and D 0 (OU-S) by a factor of more than 70.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Consistent solutions of the radiation diffusion equation in spherical and cylindrical geometries

We have extended the radiation diffusion model of Hammer and Rosen [Phys. Plasmas 10, 1829 (2003)] to diverging spherical and cylindrical geometries. The effect of curvilinear geometry on the supersonic, expanding wavefront increases as the internal radius of a spherical or cylindrical shell approaches zero. Small spherical geometries are important for modeling systems at the size scale of inertial confinement fusion capsules, at these scales existing quasi-analytic models for planar geometry significantly disagree with the results of simulation. With this method, the benefits of rapid iteration can be applied to common spherical systems at much smaller length scales. We present comparisons between numerical diffusion solutions and the analytic model to give ranges of applicability for the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of advanced vacuum technologies for extending plasma pulse duration on EAST

Advanced vacuum technologies, including pumping, fueling and wall conditioning, have been successfully developed or upgraded to efficiently control the fuel and impurity particles to extend the plasma pulse duration in the experimental advanced superconducting tokamak (EAST). To improve the particle exhaust rate cryopumps with a 60% increase in pumping speed and ∼2 times increase in saturation capacity have been developed, and molecular pumps with a ∼30% increase in pumping speed have been upgraded. In order to monitor the molecular pump status while avoiding bearing faults and overload accidents, a fault detection system has been built which can offer an early warning to avoid more losses within the fusion device. A series of fueling technologies have been developed including gas injection system, supersonic molecular beam injector, pellet injector (PI), massive gas injector and shattered pellet injector, installed at the midplane and divertor positions at different ports to improve fueling uniformity and efficiency. Meanwhile, routine wall conditioning such as electric and hot N 2 baking, ion cyclotron wall conditioning and glow discharge cleaning have been successfully developed to remove the impurity particles from the inner component and materials. The low Z material wall coating and real-time powder injection during plasma discharge are also designed and applied to further improve particle control capability. Finally, by using these advanced vacuum related technologies, good vacuum (<2 × 10 −6 Pa) and wall conditions are realized, and the fuel and impurity particles can be effectively and stably controlled, which promotes the achievement of the record plasma of ∼1056 s pulse duration with the line-averaged electron density of 1.8 × 10 19 m −3 on EAST. They provide a very important reference for vacuum system design and operation for future fusion devices.

EAST↗
DOE OSTI2025

Centrifugal-mirror confinement with strong azimuthal magnetic field

One practical challenge for the centrifugal-mirror confinement concept is the large radial voltage necessary to drive supersonic azimuthal rotation. In principle, the addition of a strong azimuthal field could reduce the required voltage, since the simple azimuthal E x B drift would be replaced by more rapid azimuthal trapped-particle precession. Also, if the mirror ratio is large enough, newly ionized ions are accelerated to the necessary parallel velocities in their first bounce orbit, both confining and significantly heating them. Unfortunately, MHD analysis shows that the centrifugal-force-confining plasma current is purely azimuthal. This implies that only the axial magnetic field contributes to the confining magnetic pressure, severely limiting the usefulness of the azimuthal magnetic field in a beta-limited plasma scenario.

Grad-Shafranov↗
DOE OSTI2025

Measurement of turbulent velocity and bounds for thermal diffusivity in laser shock compressed foams by x-ray photon correlation spectroscopy

Experimental benchmarking of transport coefficients under extreme conditions is required for validation of differing theoretical models. To date, measurement of transport properties of dynamically compressed samples remains a challenge with only a limited number of studies able to quantify transport in high pressure and temperature matter. x-ray photon correlation spectroscopy utilizes coherent x-ray sources to measure time correlations of density fluctuations, thus providing measurements of length and timescale-dependent transport properties. Here, we present a first-of-a-kind experiment to conduct x-ray photon correlation spectroscopy in laser shock compression experiments. We report measurement of the turbulent velocity in the wake of a laser driven supersonic shock and place an upper bound on thermal diffusivity in a solid density plasma on nanosecond timescales.

High-energy-density plasmas↗