Engineering PapersSearch

SEARCH · Engineering Papers

Results for “likelihood ratio test”

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.

916 records · Page 5

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi

Dark Energy Survey Year 6 results: Clustering redshifts and importance sampling of self-organized-maps 𝑛⁡(𝑧) realizations for 3 × 2 ⁢pt samples

This work is part of a series establishing the redshift framework for the 3 × 2 ⁢pt analysis of the Dark Energy Survey Year 6 (DES Y6). For DES Y6, photometric redshift distributions are estimated using self-organizing maps (SOMs), calibrated with spectroscopic and many-band photometric data. To overcome limitations from color-redshift degeneracies and incomplete spectroscopic coverage, we enhance this approach by incorporating clustering-based redshift constraints (clustering-z, or WZ) from angular cross-correlations with BOSS and eBOSS galaxies and eBOSS quasar samples. We define a WZ likelihood and apply importance sampling to a large ensemble of SOM-derived 𝑛⁡(𝑧) realizations, selecting those consistent with the clustering measurements to produce a posterior sample for each lens and source bin. The analysis uses angular scales corresponding to 1.5–5 Mpc to optimize signal-to-noise ratio while mitigating modeling uncertainties and marginalizes over redshift-dependent galaxy bias and other systematics informed by the N-body simulation CARDINAL . While a sparser spectroscopic reference sample limits WZ constraining power at 𝑧 >1.1, particularly for source bins, we demonstrate that combining SOM with WZ improves redshift accuracy and enhances the overall cosmological constraining power of DES Y6. As a result, we estimate an improvement in 𝑆 8 of approximately 10% for cosmic shear and 3 ×2⁢pt analysis, primarily due to the WZ calibration of the source samples.

Cosmological parameters

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.

The Origin of Palimpsests

Palimpsests are the most enigmatic impact features observed on icy worlds of the Solar System. Palimpsests are large impact features which are bright, old, mostly flat, and roughly circular. They range from >150 to ~350 km in diameter and sometimes exhibit recognizable associated radial secondary ejecta crater strings beyond their rims. Curiously, among all the Solar System's icy worlds, palimpsests are exclusively seen on the icy Galilean moons. Palimpsests are obviously different than similarly sized impact features seen on other icy satellites. For example, Evander, a 350 km diameter peak-ring impact basin on Saturn’s moon Dione, has a very low depth/diameter ratio of 0.008, indicating significant viscous relaxation. Nevertheless, it still retains short wavelength topography including a prominent central peak, peak ring, and a continuous rim, none of which are seen at palimpsests. Palimpsests display an approximately circular, and dominating, expanse of undulating plains, which is level with surrounding terrain and shows no recognizable rim. At the center of palimpsests are smooth central plains that are enclosed by a low scarp. Central plains material superposes the surrounding undulating plains and is littered with small blocks reminiscent of “small chaos” on Europa. Palimpsests have yet to be explored using modern impact simulation tools. We have begun to investigate two competing hypotheses: (A) Palimpsests form on targets which have a preexisting near-subsurface layer that is liquid or solid but very weak, and is excavated as proximal ejecta; or (B) They form in a solid ice target underlain by a steep thermal gradient that strongly influences the resulting impact feature, including the production of impact-induced melting at depth that is then a major component of proximal ejecta. Our results will have profound implications for the early thermal history of the icy Galilean satellites. We will test our two hypotheses for the formation of palimpsests on Ganymede (where they are both best observed and preserved) by means of 2D iSALE hydrodynamic simulations of impacts into ice targets that will be structured according to the scenarios for each hypothesis, with select runs in 3D. We have already conducted a pilot study using 2D iSALE that has yielded promising results, which we will report out at the meeting.

palimpsests

Influence of gradation on failure mode of saturated sand at particle scale

Extensive research has been reported in the literature to characterize the failure mode of dry sand using various experimental techniques such as surface optical imaging, photo-elastic materials, three-dimensional (3D) computed tomography (CT), and 3D synchrotron micro-computed tomography (SMT). However, there is a limited literature about the behavior of saturated sand. This paper presents the results of axisymmetric triaxial compression (ATC) experiments that were conducted on saturated sand specimens. The behavior of specimens composed of a uniform sand with grain size between US sieves #40 and #50 is compared to specimens conducted on the same sand that has a wider gradation. 3D SMT technique was used to acquire 3D scans while shearing the specimens to probe localized events that are completely missed or misinterpreted when analyzing ATC measurements based on global standard measurements. The results show a higher effective principal stress ratio (EPSR) for the non-uniform specimen and a thicker shear band when compared to uniform specimen.

Elnur, Mohammed [University of Tennessee]

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO

Zero-field J-spectroscopy of quadrupolar nuclei

Abstract Zero- to ultralow-field nuclear magnetic resonance (ZULF NMR) allows molecular structure elucidation via measurement of electron-mediated spin-spinJ-couplings. This study examines zero-fieldJ-spectra from molecules with quadrupolar nuclei, exemplified by solutions of various isotopologues of ammonium cations. The spectra reveal differences between various isotopologues upon extracting preciseJ-coupling values from pulse-acquire measurements. A primary isotope effect,$$\triangle J=\left({\gamma }_{{}^{14}{{{{{\rm{N}}}}}}}/{\gamma }_{{}^{15}{{{{{\rm{N}}}}}}}\right){J}_{{}^{15}{{{{{\rm{N}}}}}}{{{{{\rm{H}}}}}}}-{J}_{{}^{14}{{{{{\rm{N}}}}}}{{{{{\rm{H}}}}}}}\approx -58$$ △ J = γ 14 N / γ 15 N J 15 N H − J 14 N H ≈ − 58 mHz, is deduced by analysis of the proton-nitrogenJ-coupling ratios. This study points toward further experiments with symmetric cations containing quadrupolar nuclei, promising applications in biomedicine, energy storage, and benchmarking quantum chemistry calculations.

Science & Technology - Other Topics

Density Functional Theory Analysis of Candidate Tritium Permeation Barriers for Vanadium Alloys as a Fusion Structural Material

Vanadium alloys are being considered as a candidate structural material for the breeder blanket component of future fusion reactors due to their good thermal conductivity, low activation, and chemical compatibility with FLiBe. However, vanadium has a high hydrogen solubility, resulting in hydrogen diffusion and trapping in the vanadium structural alloy and lower than expected tritium breeding ratios. To mitigate the permeation of hydrogen into the structural alloy, a tritium permeation barrier will be coated onto the vanadium alloy. Multiple materials have been proposed as candidate tritium permeation barriers, but the extent to which each of these materials impacts hydrogen diffusion and trapping is still unknown. In this project, we perform atomic-scale calculations to narrow down the list of potential tritium permeation barrier candidates and to understand the mechanisms through which hydrogen diffusion is affected in these materials.

36 MATERIALS SCIENCE

Demonstration of gold nanorod systems for enhanced total efficiency: Experimental and numerical analysis

This study investigates the photothermal performance of gold nanorods engineered to exhibit longitudinal plasmon resonances at 695 nm, 780 nm, and 970 nm. The work combines synthesis, structural characterization, extinction measurements, numerical modeling, and controlled temperature experiments to quantify how nanorod geometry, resonance tuning, concentration, and chamber shape jointly influence heat generation. Transmission electron microscopy confirms that increasing nanorod aspect ratio systematically shifts the longitudinal plasmon peak toward the near-infrared region. Extinction measurements show strong agreement with theoretical predictions, validating the numerical model across two independent datasets. Three chamber geometries were tested under laser excitation at 640 nm, 808 nm, and 980 nm: an ascending stepped base, a flat base, and a descending stepped base. Without nanorods, the ascending geometry produced the highest efficiency due to enhanced natural convection. After introducing gold nanorods, all geometries exhibited substantial thermal enhancement, with total efficiencies exceeding 20%. The strongest improvement was obtained for nanorods resonant at 780 nm with a mass concentration of 4.6 mg/mL implemented on the descending stepped-base geometry. This performance resulted from the combined effect of spectral overlapping with the 808 nm laser, the highest nanorod concentration, and localized heat accumulation that intensified buoyancy-driven flow. The findings demonstrate that total efficiency is governed by a synergistic interplay between optical resonance, nanoparticle concentration, and macroscopic chamber design, revealing the system-level coupling between nanoscale plasmonic absorption and macroscale heat-transfer phenomena. The results provide a validated framework for tuning nanoscale plasmonic absorbers and optimizing thermal systems for applications requiring efficient light-to-heat conversion.

77 NANOSCIENCE AND NANOTECHNOLOGY

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics

Wild Blue Yonder Propulsion Schemes

This paper will include a discussion of the ORION concept, several gaseous core nuclear rockets, thermonuclear propulsion utilizing superconducting magnets, and finally a lightweight radioisotope power generation system for electric propulsion. With the exception of the latter concept, all of these schemes have much in common. The initial vehicle weights would be very large - on the order of several million pounds. The payload fractions are high - on the order of 25 to 50 percent of the takeoff weight - for near Earth missions. The development problems would be severe, and, correspondingly, the development costs would be extreme - on the order of many billions of dollars. In addition, the launching problems from Earth would be fantastic - with nuclear radiation hazards and political overtones added for good measure. However, the reward for success would be great. One can contemplate large payload fractions propelled on space missions - with thrust-to-weight ratios, at least in some cases, greater than unity and with specific impulses of several thousand seconds. The mission transportation costs would run in terms of dollars per pound of payload with clear opportunities for reasonable manned expeditions across the solar system. This is the carrot that leads the endorsement of such gigantic projects.

John C Evvard

Tunable magnetic excitations in the honeycomb antiferromagnet (Co,Ni)⁢TiO3

The solid solution of the honeycomb antiferromagnet (AFM) Co0.5⁢Ni0.5⁢TiO3 (CNTO) was synthesized by mixing a 1:1 molar ratio of CoTiO3 (CTO) and NiTiO3 (NTO). Powder neutron scattering was used to determine the structure and magnetic spectrum, where the nuclear and magnetic structures resemble those of the end members. The Néel temperature (T𝑁 = 27 ± 1 K), the ordered magnetic moment (M = 2.29 ± 0.03 𝜇𝐵 per magnetic ion), and lattice parameters are intermediate between those of the two. From inelastic neutron scattering, it is observed that the magnon density of states (MDOS) extends up to 13 meV, and the interaction can be described by an XXZ-type magnetic Hamiltonian. Above 15 meV, spin-orbit excitons (SOEs) are observed similarly to those in CTO, but with modified spin-orbit crystal-field splitting because of the substitution of Ni at Co sites, and the disorder induced by the Ni substitution modifies the spin-orbit crystal field. Thus, the magnetic dynamics of the CNTO cannot be described as a simple average of the two parent compounds.

Rathnayaka, Srimal [University of Virginia, Charlo

Valorization of consolidated bioprocessing residues for bioplastics

This study demonstrates an organic solvent-free processing strategy to valorize consolidated bioprocessing (CBP) residues, from switchgrass and poplar biomass, into functional poly(butylene succinate) (PBS)-based biocomposites using high-shear homogenization (HSH). HSH transformed the switchgrass and poplar CBP residues (CBP-R) into fine, uniformly distributed particles and microfibers. The composites of PBS with homogenized switchgrass residues (H-CBP-R-SG) or homogenized poplar residues (H-CBP-R-P) at a 70/30 weight ratio exhibited improved processability and mechanical integrity, with the Young's modulus for the PBS/H-CBP-R-SG and PBS/H-CBP-R-P nearly doubling to 0.66 ± 0.07 GPa and 0.65 ± 0.04 GPa, respectively, compared to neat PBS (0.36 ± 0.02 GPa). Dynamic Mechanical Analysis (DMA) reveals a significant suppression of the tan δ peak magnitude, indicating that HSH-mediated physical activation facilitates stress transfer in composites typical of covalent chemical grafting systems. While the transition to a stiffness-dominated profile reduces ductility, the resulting composites exhibit the dimensional stability and resistance to thermal warping required for high-fidelity FDM 3D printing and injection molding. Beyond material performance, comprehensive techno-economic analysis (TEA) and life cycle assessment (LCA) confirmed that diverting CBP residues into composite products can improve the economic viability of the biorefinery without substantially increasing biorefinery global warming potential (GWP). At a 30 wt% blend ratio, incorporating residuals into PBS yielded a minimum selling price for the composite of $\$4.07$ per kg compared to the conventional bioplastic price of $\$5.00$ per kg. This approach aligns with circular bioeconomy principles by converting waste streams into value-added products. Furthermore, this innovative strategy addresses key challenges in bioplastic development, including cost, compatibility, and performance, while simultaneously advancing waste minimization strategies for sustainable manufacturing systems.

09 BIOMASS FUELS

Evaluation of an Ejector Ramjet Based Propulsion System for Air-Breathing Hypersonic Flight

A Rocket Based Combined Cycle (RBCC) engine system is designed to combine the high thrust to weight ratio of a rocket along with the high specific impulse of a ramjet in a single, integrated propulsion system. This integrated, combined cycle propulsion system is designed to provide higher vehicle performance than that achievable with a separate rocket and ramjet. The RBCC engine system studied in the current program is the Aerojet strutjet engine concept, which is being developed jointly by a government-industry team as part of the Air Force HyTech program pre-PRDA activity. The strutjet is an ejector-ramjet engine in which small rocket chambers are embedded into the trailing edges of the inlet compression struts. The engine operates as an ejector-ramjet from takeoff to slightly above Mach 3. Above Mach 3 the engine operates as a ramjet and transitions to a scramjet at high Mach numbers. For space launch applications the rockets would be re-ignited at a Mach number or altitude beyond which air-breathing propulsion alone becomes impractical. The focus of the present study is to develop and demonstrate a strutjet flowpath using hydrocarbon fuel at up to Mach 7 conditions.

Scott R Thomas

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up

Single‐photon emitters in PECVD‐grown silicon nitride films: from material growth to photophysical properties

Abstract Silicon nitride (SiN) is a key material for quantum photonics due to its wide transparency window, high refractive index, low optical losses, and semiconductor foundry compatibility. We study the formation of single‐photon emitters in SiN films grown by plasma‐enhanced chemical vapor deposition (PECVD), exploring their photophysical properties and dependence on growth conditions. Emitters were observed across the entire range of nitrogen‐to‐silicon precursor ratios, from silicon‐rich to nitrogen‐rich conditions, enabled by the low background fluorescence. We demonstrate single‐photon emitters in SiN films with a higher refractive index (1.8–1.9) compared to our previous reports (∼1.7). Notably, nitrogen‐rich, thinner films yield particularly bright emitters with shorter emission lifetimes, likely due to more efficient annealing. Silicon‐rich SiN films exhibit red‐shifted emission, suggesting that composition may provide a mechanism for wavelength tuning. These findings establish the feasibility of emitters formation in foundry standard PECVD tools, advancing the scalability and lab‐to‐fab transition of SiN‐based quantum photonic technologies.

Materials Science

Air-Broadened Line Shapes in the 2n3 R Branch of 12ch4 Between 6014 and 6100 cm-1

Complete and accurate information on line shape parameters of 2?3 methane transitions for air broadening as a function of temperature is critical not only for the correct interpretation of the observed atmospheric spectra but also for the development of a reliable theoretical model. For this reason, we obtained a series of high-resolution, high S/N spectra of high-purity 12CH4 and 12CH4 broadened with dry air at temperatures in the 130 to 295 K range using the Bruker IFS 125HR Fourier transform spectrometer at JPL. Two absorption cells were used in the experiment, a White cell with path length of 13 m for room temperature spectra and a 21 m Herriott cella for cold sample spectra. The 15 spectra used in the analysis consisted of 3 low pressure (0.26 to 2.57 Torr) spectra with pure 12CH4 and 12 air-broadened spectra with total sample pressures of 79-805 Torr and volume mixing ratios of methane between 0.23 and 1%. A multispectrum least-squares fitting techniqueb was employed to fit all 15 spectra simultaneously. Preliminary results for select R(J) manifolds will be presented.

V Malathy Devi