Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLANE WAVE”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Calculated electron paramagnetic resonance g-tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO

Various defects in ZnO, focused on substitutional N O and N 2 in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N 2 in ZnO and substitutional N on the O-site. The g tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The g-tensor of the free $N^{+}_{2}$ and $N^{–}_{2}$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N 2 behavior in ZnO. Previous controversies on the site location of N 2 in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. Here, we find that the N 2 on the Zn site is mostly zinc-vacancy like in its spin density and g-tensor, while for the O-site, a model with the N 2 axis lying in-the basal plane and the singly occupied π g -orbital along the c axis provides good agreement with experiment. For the interstitial location, if the N 2 is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated g-tensors for N O and V Zn are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the g-tensor values.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An assessment of spall failure modes in laser powder bed fusion fabricated stainless steel 316L with low-volume intentional porosity

Here, this paper reports on spall failure and damage modes in Laser Powder Bed Fusion fabricated Stainless Steel 316L (SS316L) with intentional levels of low-volume (1–5 vol. %) porosity and pore sizes of 200, 350, and 500 μm. The fabricated specimens were subjected to uniaxial-strain plate-impact loading at ~4.5 GPa, to initiate incipient spall failure. Analysis of velocimetry profiles measured using multi-probe photon-Doppler velocimetry coupled with post-mortem analysis of soft-recovered samples reveals local suppression of spall failure (termed as spall-dominated) as a function of porosity, as the failure mechanism transitions from spall-centered tensile stress dominated to a pore-centered microstructure-dominated damage mode involving void/crack nucleation and growth at pre-existing pores. The critical porosity level where the suppression of spall failure is first observed, as well as the spall location, is dependent on both the volume fraction and the size of the initially fabricated pores. In samples of 500 μm pore size, the suppression of spall failure is observed with as little as 1 vol. % porosity, while samples with smaller pores (200 μm) still experience spall-centered tensile stress dominated failure with higher levels (5 vol. %) of porosity. In the case of pore-centered microstructure-dominated failure, spall damage can occur but the spall plane is shifted toward the rear free surface, or more generally in areas further away from the region with pores. Highly heterogeneous deformation twinning, shear banding, grain rotation, and cracking are observed in the vicinity of pre-existing pores and expected spall failure sites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

4-Electron Oxygen Reduction Reaction (ORR) with Iron Phthalocyanine (FePc) Functionalized Nanowire Templated-3D Fuzzy Graphene (NT-3DFG)

Iron phthalocyanine (FePc) is a promising alternative to platinum-based catalysts for sustainable energy devices; however, the plane-symmetry of Fe-N 4 sites, random aggregation, and poor conductivity of FePc present major barriers for their application as oxygen reduction reaction (ORR) electrocatalysts. Here, we report the synergistic effects of FePc electrocatalysts supported by a nanowire-templated three-dimensional fuzzy graphene (FePc@NT-3DFG) substrate. The in situ functionalized oxygen groups (iFOGs) at the edge of NT-3DFG localize Fe active sites in FePc under alkaline ORR conditions. With a uniform FePc distribution through many single layers of graphene, the NT-3DFG substrates improve O 2 adsorption and catalytic activity while stabilizing the electrochemical activity during reactions. The FePc@NT-3DFG catalyst exhibits fast ORR kinetics with an extremely low Tafel slope of 28.3 ± 2.7 mV dec −1 , a higher half-wave potential of 0.911 ± 0.004 V (vs RHE), and notable long-term stability at 0.5 V (vs RHE) of 96.0 ± 0.4% retention after 30 h. Surface chemistry spectra validate electronic configuration modification of Fe at the iFOGs. Density functional theory calculations indicate that the extra layers of graphene improve oxygen adsorption. Moreover, additional exploration of other transition metal phthalocyanines supports the effects of iFOGs through the transition toward 4e − ORR. This work offers an expanded strategy for active site modification through edge-based graphene substrates for 4e − ORR.

X-ray absorption spectroscopy↗
DOE OSTI2025

Anisotropic light-tailored RKKY interaction in two-dimensional 𝑑-wave altermagnets

Altermagnets are known in spintronics for their intrinsic spin-splitting and unconventional magnetic responses, particularly to magnetic impurities. However, effectively controlling the magnetic exchange interactions in altermagnets is challenging for practical applications. Here, in this work, we propose using circularly polarized light to tune the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in two-dimensional 𝑑-wave altermagnets. Using the real-space retarded Green's functions approach, our results show that while the Heisenberg and Ising exchanges dominate, a notable Dzyaloshinskii–Moriya (DM) interaction also plays a key role. Furthermore, the inherent strength of altermagnetism imprints chirp-like signatures into the magnetic responses, which can be dynamically tuned via light. We mainly demonstrate that gate-induced Rashba spin-orbit coupling is essential in response to light—light selectively and anisotropically adjusts the DM interaction without affecting the other exchanges. Our findings further indicate that rotating the altermagnet by 45° relative to the light's polarization direction generates a Dirac-like dispersion and different DM interactions. We finally extract critical thresholds where light reverses DM interactions along one axis or balances both in-plane components. The anisotropic light-driven control of RKKY interactions in altermagnets not only highlights their unique properties but also opens new avenues for engineering tailored magnetic characteristics in spintronic applications.

altermagnetism↗
DOE OSTI2020

Near-canopy horizontal concentration heterogeneity of semivolatile oxygenated organic compounds and implications for 2-methyltetrols primary emissions

Semivolatile oxygenated organic compounds (SV-OVOCs) are important atmospheric species, in particular for the production and chemistry of atmospheric particulate matter and related impacts on air quality and climate. In this study, SV-OVOCs were collected in the horizontal plane of the roughness layer over the tropical forest in the central Amazon during the wet season of 2018. A sampler mounted to a copter-type, hovering unmanned aerial vehicle was used. Underlying the collection region, a plateau forest transitioned into a slope forest across several hundred meters. The concentrations of pinonic and pinic acids, which are monoterpene oxidation products, had no statistical difference over the two forests. By comparison, across the study period, differences in the concentration of 2-methyltetrols, which are products of isoprene oxidation, ranged from -70% to +480% over the two forests. The chemical lifetime of 2-methyltetrols in the atmosphere is sufficiently long that heterogeneity in the isoprene emission rate from the two forests followed by atmospheric oxidation does not explain the concentration heterogeneity of 2-methyltetrols. Standing waves and local meteorology also cannot account for the heterogeneity. Of the possibilities considered, the most plausible explanation is the direct emission from the forest of 2-methyltetrols produced through biological processes within the plants. Under this explanation, the rate of direct SV-OVOC emissions should be modulated by forest type and related environmental stressors. Direct emissions of SV-OVOCs should be more broadly considered for constraining and improving models of atmospheric composition, transport, and chemistry over tropical forests.

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

Review of Radio Frequency Interference and Potential Impacts on the CMB-S4 Cosmic Microwave Background Survey

CMB-S4 will map the cosmic microwave background to unprecedented precision, while simultaneously surveying the millimeter-wave time-domain sky, in order to advance our understanding of cosmology and the universe. CMB-S4 will observe from two sites, the South Pole and the Atacama Desert of Chile. A combination of small- and large-aperture telescopes with hundreds of thousands of polarization-sensitive detectors will observe in several frequency bands from 20-300 GHz, surveying more than 50% of the sky to arcminute resolution with unprecedented sensitivity. CMB-54 seeks to make a dramatic leap in sensitivity while observing across a broad range of largely unprotected spectrum which is increasingly being utilized for terrestrial and satellite transmissions. Fundamental aspects of CMB instrument technology leave them vulnerable to radio frequency interference (RFI) across a wide range of frequencies, including frequencies outside of their observing bands. Ground-based CMB instruments achieve their extraordinary sensitivities by deploying large focal planes of superconducting bolometers to extremely dry, high-altitude sites, with large fractional bandwidths, wide fields of view, and years of integration time. Suitable observing sites have historically offered significant protection from RFI, both naturally through their extremely remote locations as well as through restrictions on local emissions. Since the coupling mechanisms are complex, "safe" levels or frequencies of emission that would not interfere with CMB measurements cannot always be determined through straightforward calculations. We discuss models of interference for various types of RFI relevant to CMB-54, mitigation strategies, and the potential impacts on survey sensitivity.

Cosmic Microwave Background↗
DOE OSTI2024

Energy amplification in plasma crystals due to multiple torsions

Interacting torsions are examined within a two-dimensional monolayer crystal suspended in an argon complex plasma for 1–10 W discharge powers and pressures of 135–155 mTorr. Two torsions embedded in a lattice are shown to amplify the kinetic energy and range of motion of particles located between the torsions to nearly three times that observed in single torsion systems. It is also shown that multiple torsions can interact via amplified particle energy when separated by up to 14 interparticle distances (Δ). The torsion separation distance also showed a positive linear trend with power and a slightly positive correlation with the pressure. This amplification of energy is possible due to the fact that multiple torsions in a lattice increase the interparticle distance of the lattice by 16% more than single torsion systems, leading to additional freedom of motion in the lattice plane. These combined findings show that multiple torsions heat the lattice differently depending on their separation from the other torsion. The midpoint particles between torsions absorb the majority of energy from the two torsions, and energy addition at the midpoint is nonlinear. The addition of more torsions to the lattice may lead to melting of the plasma crystal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The estimation of Poisson’s ratio by time-averaging and Cornu’s method for isotropic beams

Prior research suggests that direct static (e.g., uniaxial testing) and dynamic (e.g., ultrasonic wave propagation analysis) measurements differ in the estimation of Poisson’s ratio because anisotropies and heterogeneities in the sample material affect the two types of tests differently. Even assuming isotropic and homogeneous material properties, prior research further suggests that discrepancies between static/dynamic test results will exist because the error of the diagnostic techniques for the measurand are inherently different. Finally, thermodynamic effects are not present in static tests but can significantly affect dynamic test results. Given the potential for all these variables to produce discrepancies, it would be helpful to have the measurement of Poisson’s ratio obtainable from the same theory and experimental measurements by either static or dynamic testing methods. Our finite element calculations show that by combining time-averaged scanning digital holography with Cornu’s method, it is theoretically possible to estimate the effective Poisson’s ratio from the anticlastic contours at the antinode of the first out-of-plane bending mode shape. This is true regardless of frequency and therefore applicable for both static and dynamic measurements. Our results show that the estimate of Poisson’s ratio by Cornu’s method using data from simulations of mode shapes approaches the true value of Poisson’s ratio. Additionally, our research suggests that beam geometry and boundary conditions are fundamental factors limiting the convergence of the estimate of Poisson’s ratio to the true value of Poisson’s ratio regardless of performing a static or dynamic test.

42 ENGINEERING↗
DOE OSTI2021

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

Acoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically-tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near-vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be < 2%. While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

58 GEOSCIENCES↗
DOE OSTI2022

Measurement of shock roughness due to phase plate speckle imprinting relevant for x-ray diffraction experiments on 3rd and 4th generation light sources

Laser-shock compression experiments at 3rd and 4th generation light sources generally employ phase plates, which are inserted into the beamline to achieve a repeatable intensity distribution at the focal plane. Here, the laser intensity profile is characterized by a high-contrast, high-frequency laser speckle. Without sufficient smoothing, these laser non-uniformities can translate to a significant pressure distribution within the sample layer and can affect data interpretation in x-ray diffraction experiments. Here, we use a combination of one- and two-dimensional velocity interferometry to directly measure the extent to which spatial frequencies within the laser focal spot intensity pattern are smoothed out during propagation within the laser plasma and a polyimide ablator. We find that the use of thicker polyimide layers results in spatially smoother shock fronts, with the greatest degree of smoothing associated with the highest spatial frequencies. Focal spots with the smallest initial speckle separation produce the most rapid smoothing. Laser systems that employ smoothing by spectral dispersion techniques to rapidly modulate the focal plane intensity distribution are shown to be the most effective ones in producing a spatially smooth shock front. We show that a simple transport model combined with the known polyimide Hugoniot adequately describes the extent of shock smoothness as a function of polyimide thickness. Finally, our results provide a description of spatial structure smoothing across a shock front, which can be used to design targets on x-ray free electron laser facilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Deep Minima in the Triply Differential Cross Section for Ionization of Atomic Hydrogen by Electron and Positron Impact

We investigate ionization of atomic hydrogen by electron- and positron-impact. We apply the Coulomb–Born (CB1) approximation, various modified CB1 approximations, the three body distorted wave (3DW) approximation, and the time-dependent close-coupling (TDCC) method to electron-impact ionization of hydrogen. For electron-impact ionization of hydrogen for an incident energy of approximately 76.45 eV, we obtain a deep minimum in the CB1 triply differential cross section (TDCS). However, the TDCC for 74.45 eV and the 3DW for 74.46 eV gave a dip in the TDCS. For positron-hydrogen ionization (breakup) we apply the CB1 approximation and a modified CB1 approximation. We obtain a deep minimum in the TDCS and a zero in the CB1 transition matrix element for an incident energy of 100 eV with a gun angle of 56.13 °. Corresponding to a zero in the CB1 transition matrix element, there is a vortex in the velocity field associated with this element. For both electron- and positron-impact ionization of hydrogen the velocity field rotates in the same direction, which is anticlockwise. All calculations are performed for a doubly symmetric geometry; the electron-impact ionization is in-plane and the positron-impact ionization is out-of-plane.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Evolution of magnetism in the magnetic topological semimetal NdSb x Te 2–x+δ

Magnetic topological semimetals LnSbTe (Ln = Lanthanide) have attracted intensive attention because of the presence of interplay between magnetism, topological, and electron correlations depending on the choices of magnetic Ln elements. Recently, varying Sb-Te composition has been found to effectively control the electronic and magnetic states in LnSb x Te 2-x . Here, with this motivation, we report the evolution of magnetic properties with Sb-Te substitution in NdSb x Te 2-x+δ , (0 ≤ x ≤ 1). Our work reveals the interesting non-monotonic change in magnetic ordering temperature with varying composition stoichiometry. In addition, reducing the Sb content x drives the reorientation of moments from in-plane (ab-plane) to out-of-plane (c-axis) direction that results in the distinct magnetic structures for two end compounds NdTe 2 (x = 0) and NdSbTe (x = 1). Furthermore, the moment orientation in NdSb x Te 2-x+δ is also found to be strongly tunable upon application of weak magnetic field, leading to rich magnetic phases depending on the composition stoichiometry, temperature, and magnetic field. Such strong tuning of magnetism in this material establishes it as a promising platform for investigating tunable topological states and correlated topological physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nature of the magnetic stripes in fully oxygenated La 2 CuO 4+y

We present triple-axis neutron scattering studies of static and dynamic magnetic stripes in an optimally oxygen-doped cuprate superconductor, La 2 CuO 4+y , which exhibits a clean superconducting transition at T c = 42 K. Polarization analysis reveals that the magnetic stripe structure is equally represented along both of the tetragonal crystal axes and that the fluctuating stripes display significant weight for in-plane as well as out-of-plane spin components. Both static magnetic order as well as low-energy fluctuations are fully developed in zero applied magnetic field and the low-energy spin fluctuations at ℏω = 0.3 – 10 meV intensify on cooling. We interpret this as an indication that superconductivity and low-energy spin fluctuations coexist microscopically in spatial regions which are separated from domains with static magnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spatially correlated incommensurate lattice modulations in an atomically thin high-temperature Bi 2.1 Sr 1.9 CaCu 2.0 O 8+y superconductor

Strong variations in superconducting critical temperatures in different families of the cuprate perovskites, even with similar hole doping in their copper-oxygen planes, suggest the importance of lattice modulation effects. The one-dimensional incommensurate lattice modulation (ILM) of Bi 2 Sr 2 CaCu 2 O 8+y , with the average atomic positions perturbed beyond the unit cell6, offer s an ideal test ground for studying the interplay between superconductivity and the long-range incommensurate lattice fluctuations. Here we report Scanning nano X-ray Diffraction (SnXRD) imaging of incommensurate lattice modulations in Bi 2 Sr 2 CaCu 2 O 8+y Van der Waals heterostructures of thicknesses down to two-unit cells. Using SnXRD, we probe that the long-range and short-range incommensurate lattice modulations in bulk sample surface with spatial resolution below 100 nm. We find that puddle-like domains of ILM of size uniformly evolving with dimensionality. In the 2-unit cell thin sample, it is observed that the wavevectors of the long- and short-range orders become anticorrelated with emerging spatial patterns having a directional gradient. The emerging patterns, originated by tiny tuning of lattice strain, induce static mesoscopic charge density waves. Finally, our findings thus demonstrate that the strain can be used to tune and control the electromagnetic properties of two-dimensional high-temperature superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Electrical Transport Interplay with Charge Density Waves, Magnetization, and Disorder Tuned by 2D van der Waals Interface Modification via Elemental Intercalation and Substitution in ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 Crystals

Electrical transport in 2D materials exhibits unique behaviors due to reduced dimensionality, broken symmetries, and quantum confinement. It serves as both a sensitive probe for the emergence of coherent electronic phases and a tool to actively manipulate many-body correlated states. Exploring their interplay and interdependence is crucial but remains underexplored. This review integratively cross-examines the atomic and electronic structures and transport properties of van der Waals-layered crystals ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 , providing a comprehensive understanding and uncovering new discoveries and insights. A common observation from these crystals is that modifying the atomic and electronic interface structures of 2D van der Waals interfaces using heteroatoms significantly influences the emergence and stability of coherent phases, as well as phase-sensitive transport responses. In ZrTe 3 , substitution and intercalation with Se, Hf, Cu, or Ni at the 2D vdW interface alter phonon–electron coupling, valence states, and the quasi-1D interface Fermi band, affecting the onset of CDW and SC, manifested as resistance upturns and zero-resistance states. We conclude here that these phenomena originate from dopant-induced variations in the lattice spacing of the quasi-1D Te chains of the 2D vdW interface, and propose an unconventional superconducting mechanism driven by valence fluctuations at the van Hove singularity, arising from quasi-1D lattice vibrations. Short-range in-plane electronic heterostructures at the vdW interface of Cr 2 Si 2 Te 6 result in a narrowed band gap. The sharp increase in in-plane resistance is found to be linked to the emergence and development of out-of-plane ferromagnetism. The insertion of 2D magnetic layers such as Mn, Fe, and Co into the vdW gap of 2H-TaS 2 induces anisotropic magnetism and associated transport responses to magnetic transitions. Overall, 2D vdW interface modification offers control over collective electronic behavior, transport properties, and their interplays, advancing fundamental science and nanoelectronic devices.

2D van der Waals interfaces↗
DOE OSTI2022

The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations

Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.

42 ENGINEERING↗