Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compiler techniques”

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 73 records · Page 4

Determination of β -decay feeding patterns of 88 Rb and 88 Kr using the Modular Total Absorption Spectrometer at ORNL HRIBF

We report precise determination of ground-state feeding in the β decay of fission products is an important but challenging component in modeling reactor antineutrino flux and reactor decay heat. The Modular Total Absorption Spectrometer (MTAS) is a versatile NaI(Tl) detector array that determines the true β-decay pattern free from the pandemonium effect, including precise ground-state feeding intensities. In this paper, we report MTAS results of the β feeding intensities of 88 Rb and 88 Kr, fission products with large cumulative yields in nuclear reactors. By comparing MTAS results with previous measurements, 88 Rb provides a validation of MTAS's ability to determine ground-state feedings in β decays, while the precision of 88 Kr ground-state feeding is improved when compared with the Evaluated Nuclear Structure Data File (ENSDF). The investigation of sources that contribute to β feeding branching uncertainties in MTAS experiments is discussed in detail. Lastly, the deconvolution of 88 Rb decay spectra suggests that MTAS can distinguish an allowed β spectral shape from a first forbidden unique β spectral shape.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Binary Control Pulse Optimization for Quantum Systems

Quantum control aims to manipulate quantum systems toward specific quantum states or desired operations. Designing highly accurate and effective control steps is vitally important to various quantum applications, including energy minimization and circuit compilation. In this paper we focus on discrete binary quantum control problems and apply different optimization algorithms and techniques to improve computational efficiency and solution quality. Specifically, we develop a generic model and extend it in several ways. We introduce a squared L 2 -penalty function to handle additional side constraints, to model requirements such as allowing at most one control to be active. We introduce a total variation (TV) regularizer to reduce the number of switches in the control. We modify the popular gradient ascent pulse engineering (GRAPE) algorithm, develop a new alternating direction method of multipliers (ADMM) algorithm to solve the continuous relaxation of the penalized model, and then apply rounding techniques to obtain binary control solutions. We propose a modified trust-region method to further improve the solutions. Our algorithms can obtain high-quality control results, as demonstrated by numerical studies on diverse quantum control examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Dependence of the Type Ia Supernova Host Bias on Observation or Fitting Technique

More luminous Type Ia supernovae prefer less massive hosts and regions of higher star formation. This correlation is inverted during width–color–luminosity light-curve standardization resulting in step-like biases of distance measurements with respect to host properties. Using the PMAS/PPak Integral-field Supernovahosts COmpilation (PISCO) supernova host sample and Sloan Digital Sky Survey, Galaxy Evolution Explorer, and Two Micron All Sky Survey photometry, we compare host stellar mass and specific star-formation rate (sSFR) from different observation methods, including local versus global, and fitting techniques to measure their impact on the host step biases. Mass-step measurements for all our mass samples are consistent within a 1σ significance from –0.03 ± 0.02 mag to –0.04 ± 0.02 mag. Including or excluding UV information had no effect on measured mass-step size or location. sSFR step sizes are more significant than mass-step measurements and varied from 0.05 ± 0.03 mag (Hα) and 0.06 ± 0.02 mag (UV) for a 51 host sample. The sSFR step location is influenced by the mass sample used to normalize star formation and by sSFR tracer choice. The step size is reduced to 0.04 ± 0.03 mag when using all available 73 hosts with Hα measurements. This 73 PISCO host subsample overall lacked a clear step signal, but here we are searching for whether different choices of mass or sSFR estimation can create a step signal. We find no evidence that different observation or fitting techniques choices can create a distance measurement step in either mass or sSFR.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular Symmetry in VQE: A Dual Approach for Trapped-Ion Simulations of Benzene

Understanding complex chemical systems—such as biomolecules, catalysts, and novel materials—is a central goal of quantum simulations. Near-term strategies hinge on the use of variational quantum eigensolver (VQE) algorithms combined with a suitable ansatz. However, straightforward application of many chemically-inspired ansatze yields prohibitively deep circuits. In this work, we employ several circuit optimization methods tailored for trapped-ion quantum devices to enhance the feasibility of intricate chemical simulations. The techniques aim to lessen the depth of the unitary coupled cluster with singles and doubles (uCCSD) ansatz's circuit compilation, a considerable challenge on current noisy quantum devices. Furthermore, we use symmetry-inspired classical post-selection methods to further refine the outcomes and minimize errors in energy measurements, without adding quantum overhead. Our strategies encompass optimal mapping from orbital to qubit, term reordering to minimize entangling gates, and the exploitation of molecular spin and point group symmetry to eliminate redundant parameters. The inclusion of error mitigation via post-selection based on known molecular symmetries improves the results to near milli-Hartree accuracy. These methods, when applied to a benzene molecule simulation, enabled the construction of an 8-qubit circuit with 69 two-qubit entangling operations, pushing the limits for variational quantum eigensolver (VQE) circuits executed on quantum hardware to date. 1 1 This manuscript has been authored in part by UT-Battelle, LLC, under contract DE-AC05-000R22725 with the US Department of Energy (DOE). The publisher acknowledges the US government license to provide public access under the DOE Public Access Plan(https://energy.gov/doe-public-access-plan).

Goings, Joshua↗

NEML2: A High Performance Library for Constitutive Modeling

NEML2, the New Engineering Material model Library, version 2, is an offshoot of NEML, an earlier material modeling code developed at Argonne National Laboratory. NEML2 extends the key philosophy of its predecessor, i.e., material models are flexible, modular, and can be built from smaller blocks. It also provides modern features that do not exist in the framework of its predecessor such as material model vectorization, automatic differentiation, device-portable just-in-time compilation, operator fusion, lazy tensor evaluation, etc. Moreover, NEML2 can seamlessly integrate with the popular machine learning package PyTorch to take advantage of modern and fast-growing machine learning techniques. In this fiscal year, the development of core library features and capabilities are complete. The purpose of this report is not to serve as a verbatim copy of the software API reference (which is available online at https://reverendbedford.github.io/neml2/). Instead, this report documents the motivation, implementation, design choices, and usage of each core capability as well as their applications in solving practical engineering problems. This report is compiled based on the NEML2 major release 2.0.0.

36 MATERIALS SCIENCE↗

Expanded analysis of machine learning models for nuclear transient identification using TPOT

Industries around the world are becoming more and more data driven. The nuclear field is no exception with several different applications being proposed. One popular area of research is the use of machine learning in transient detection. This paper seeks to build upon a previous study which made use of the AutoML package TPOT to train traditional machine learning models to classify transient events occurring with a reactor. Synthetic data was once again collected using a GPWR reactor simulator. Data on 12 different events was collected using 15 different initial conditions. Here, a dataset consisting of over 100,000 data points was compiled and used to train 7 different machine learning models using a pre-defined TPOT dictionary with 12 different preprocessing techniques. Three of the trained models were able to produce validation results in the 90s with the expanded dataset. Once the models were trained, it was possible to look into where during the simulation, misclassifications occurred. Using these three models, analysis was done to determine if TPOT could be used to train models that were effective if important features were missing. The results from this were positive with the newly trained models scoring close to the original models. Finally, to conclude this study, the three high performing models were retrained using different random states to see if there was any major variation when different states were used.

42 ENGINEERING↗

Hourly gap-filled meteorological data from PIE LTER measurements (2004-2023) used as drivers to run ELM PFLOTRAN simulations

This dataset contains continuous gap-filled precipitation, solar radiation, photosynthetically active radiation (PAR), air temperature, relative humidity, wind speed, and barometric pressure data recorded primarily at the Marshview Farm weather station within the Plum Island Long Term Ecosystems Research (PIE LTER) in Newbury Massachusetts (MA) from 2004 to 2023. We compiled the data set from published annual data packages in 15min resolution available on DataOne. Gaps were filled using different statistical techniques or available observations from the vicinity, e.g. the US-PLo and the US-PHM Ameriflux sites, also located within the PIE LTER. Flags are included in this dataset to indicate the origin of each data point. Metadata files ELMPFLOTRAN_met_dd.csv and ELMPFLOTRAN_met_flmd.csv contain more information on site locations, gap filling protocols, data variables, flags, and QA/QC methods. The data set was used in the spin up and simulations of a land surface model coupled to a biogeochemical reaction network (ELM PFLOTRAN) assessing impacts of hydrology and salinity input on methane fluxes in 2022 and 2023 (Sulman et al., 2024).

54 ENVIRONMENTAL SCIENCES↗

Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning

This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.

42 ENGINEERING↗

SPEL: Software tool for Porting E3SM Land Model with OpenACC in a Function Unit Test Framework

Most high-end computers adopt hybrid architecture, porting a large-scale scientific code onto accelerators is necessary. The paper presents a generic method for porting large-scale scientific code onto accelerators using compiler directives within a modularized function unit test platform. We have implemented the method and designed a software tool (SPEL) to port the E3SM Land Model (ELM) onto the GPUs in the Summit computer. SPEL automatically generates GPU-ready test modules for all ELM functions, such as CanopyFlux, SoilTemperature, and EcosystemDynamics. SPEL breaks the ELM into a collection of standalone unit test programs for easy code verification and further performance improvement. We further optimize several ELM test modules with advanced techniques, including memory reduction, reconstructed parallel loops, and asynchronous GPU kernel launch. We hope our study will inspire new toolkit developments that expedite large-scale scientific code porting with compiler directives.

Schwartz, Peter↗

Differential cross sections of the 16 O (𝑛,𝛼) reaction at neutron energies from 3.8 to 15 MeV

The importance of studying the 16 O (𝑛,𝛼) reaction is motivated by multiple nuclear applications which rely on an accurate nuclear reaction data library for oxygen. Discrepancies between past experimental data on the 16 O (𝑛,𝛼) 13 C reaction and its time-reverse 13 C ⁡(𝛼,𝑛) 16 O reaction have led to various different nuclear data evaluations. Here, we have measured 16 O (𝑛,𝛼) reaction cross sections using the LENZ instrument with the unmoderated white neutron source at LANSCE. Results from 2016/2017 data are discussed and used to benchmark the mcnp and geant simulations of the LENZ experimental setup. We report partial differential cross sections of 16 O (𝑛,𝛼 0 ) at 𝐸 𝑛 = 3.8–15 MeV and 16 O ⁡(𝑛,𝛼 1 +𝛼 2 +𝛼 3 ) at 𝐸 𝑛 =9–15 MeV, based on the new measurement in 2021. The resonances that we observed are in good agreement with the levels in 17 O that were previously measured. The LENZ cross sections are in better agreement with the ENDF/B-VIII.0 evaluation than the with reduced cross section found in ENDF/B-VII.1 up to 6 MeV. However, the current results for 16 O ⁡(𝑛,𝛼 0 ) 13 C g.s. appear to be in the best agreement over the entire energy range with that of JENDL/AN-2005 (ENDF/B-VI.0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

TESS Hunt for Young and Maturing Exoplanets (THYME). V. A Sub-Neptune Transiting a Young Star in a Newly Discovered 250 Myr Association

The detection and characterization of young planetary systems offer a direct path to study the processes that shape planet evolution. We report on the discovery of a sub-Neptune-sized planet orbiting the young star HD 110082 (TOI-1098). Transit events we initially detected during TESS Cycle 1 are validated with time-series photometry from Spitzer. High-contrast imaging and high-resolution, optical spectra are also obtained to characterize the stellar host and confirm the planetary nature of the transits. The host star is a late-F dwarf (M {sub ⋆} = 1.2M {sub ⊙}) with a low-mass, M dwarf binary companion (M {sub ⋆} = 0.26M {sub ⊙}) separated by nearly one arcminute (∼6200 au). Based on its rapid rotation and Lithium absorption, HD 110082 is young, but is not a member of any known group of young stars (despite proximity to the Octans association). To measure the age of the system, we search for coeval, phase-space neighbors and compile a sample of candidate siblings to compare with the empirical sequences of young clusters and to apply quantitative age-dating techniques. In doing so, we find that HD 110082 resides in a new young stellar association we designate MELANGE-1, with an age of 250{sub −70}{sup +50} Myr. Jointly modeling the TESS and Spitzer light curves, we measure a planetary orbital period of 10.1827 days and radius of R {sub p} = 3.2 ± 0.1R {sub ⊕}. HD 110082 b’s radius falls in the largest 12% of field-age systems with similar host-star mass and orbital period. This finding supports previous studies indicating that young planets have larger radii than their field-age counterparts.

47 OTHER INSTRUMENTATION↗

Laser ablation spectrometry for studies of uranium plasmas, reactor monitoring, and spent fuel safety

Nuclear security is one of the defining challenges of our time. Nuclear threats range from deliberate dispersal of radioactive material to contaminate the vital infrastructure to diversion and smuggling of special nuclear material for clandestine nuclear programs and nuclear terrorism, respectively. There is an associated need to develop and sustain the nuclear forensics capabilities, which requires the understanding of complex processes that occur in plasmas of nuclear materials. The area of nuclear safety has seen a resurgence of public interest, and there is a concomitant need to safely store used nuclear fuel and detect structural material failure in nuclear power systems, especially in innovative reactor designs envisioned for future adoption. Laser-produced plasmas are complicated extreme environments that can generate intense and rich, highly specific signatures of nuclear and radiological materials, which can then be explored in a wide range of applications. They include interdiction and rapid detection of nuclear materials, including their isotopic composition, detection over long distances, laboratory simulation of weapons effects, monitoring the condition of structural materials in dry cask storage containers, and novel instrumentation for nuclear power systems. We present a compilation of recent representative examples of the application of laser spectroscopy, and laser-induced breakdown spectroscopy in particular, to nuclear safety and security problems. A case is made that spectroscopic techniques based on laser-produced plasmas offer complementary, and sometimes unique, capabilities that motivate the continued exploration of their efficient production and understanding of the signatures they produce

(020.3260) Isotope shifts, (140.3440) Laser-induce↗

Statistical chronometry of Meteorites: II. Initial abundances and homogeneity of short-lived radionuclides

Astrophysical models of planet formation require accurate radiometric dating of meteoritic components by short-lived (Al-Mg, Mn-Cr, Hf-W) and long-lived (U-Pb) chronometers, to develop a timeline of such events in the solar nebula as formation of Ca-rich, Al-rich Inclusions (CAIs), chondrules, planetesimals, etc. CAIs formed mostly around a time (“t=0”) when the short-lived radionuclide 26 Al (t 1/2 = 0.72 Myr) was present and presumably homogeneously distributed at a known level we define as ( 26 Al/ 27 Al) SS = 5.23 × 10 -5 . The time of formation after t=0 of another object can be found by determining its initial ( 26 Al/ 27 Al) 0 ratio and comparing it to ( 26 Al/ 27 Al) SS . Dating of meteoritic objects using the Mn-Cr or Hf-W systems is hindered because the abundances ( 53 Mn/ 55 Mn) SS and ( 182 Hf/ 180 Hf) SS at t=0 are not known precisely. To constrain these quantities, we compile literature Al-Mg, Mn-Cr, Hf-W and Pb-Pb data for 13 achondrites and use novel statistical techniques to minimize the discrepancies between their times of formation across these systems. We find that for ( 53 Mn/ 55 Mn) SS = (7.80±0.36)×10 -6 , ( 182 Hf/ 180 Hf) SS = (10.41 ± 0.12) × 10-5, t SS = 4568.65 ± 0.10 Myr, and a 53 Mn half-life of 3.98 ± 0.22 Myr, these four free parameters make concordant 18 formation times recorded by the different systems in all six known volcanic achondrites (the quenched angrites D’Orbigny, Sahara 99555, Northwest Africa 1670; the psuedo-eucrites Asuka 881394 and Ibitira; and the basaltic achondrite Northwest Africa 7325). These parameters also make concordant the ages derived for chondrules from CB/CH achondrites, formed simultaneously in an impact. The other seven achondrites are not quite concordant, but are plutonic angrites or ‘carbonaceous achondrites’ for which simultaneous closure of the isotopic systems might not be expected. Our findings provide very strong support for homogeneity of 26 Al, 53 Mn, and 182 Hf in the solar nebula, and our approach offers a path for more precise chronometry

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics of qudit gates and effects of spectator modes on optimal control pulses

Qudit gates for high-dimensional quantum computing can be synthesized with high precision using numerical quantum optimal control techniques. Large circuits are broken down into modules and the tailored pulses for each module can be used as primitives for a qudit compiler. Application of the pulses of each module in the presence of extra modes may decrease their effectiveness due to crosstalk. In this paper, we address this problem by simulating qudit dynamics for circuit quantum electrodynamics systems. As a test case, we take pulses for single-qudit swap gates optimized in isolation and then apply them in the presence of spectator modes each of which are in Fock states. We provide an experimentally relevant scaling formula that can be used as a bound on the fidelity decay. Our results show that frequency shift from spectator mode populations has to be ≲ 0.1 % of the qudit's nonlinearity in order for high-fidelity single-qudit gates to be useful in the presence of occupied spectator modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MLIR loop optimizations for High-Level Synthesis: a case study

High-Level Synthesis (HLS) tools simplify the design of hardware accelerators by automatically generating Verilog/VHDL code starting from a general purpose software programming language. They include a wide range of optimization techniques in the process, most of them performed on a low-level intermediate representation (IR) of the code. Introducing optimizations on a higher level of abstraction could significantly contribute to the automated design process results; for example, polyhedral techniques for the manipulation of loops could have a significant impact on the generated accelerators when applied on a specialized IR. We use loop pipelining as a case study to explore the introduction of compiler-based transformations on top of an existing HLS process. We leverage the Multi-Level Intermediate Representation (MLIR) framework and an external scheduler to implement the required transformations, and couple them with existing HLS tools to evaluate the improvements that loop pipelining brings to the performance of generated accelerators. The proposed approach can be integrated with other high-level transformations on the MLIR representation, combining different techniques to obtain pre-optimized inputs for HLS that do not have to rely on a specific backend tool.

Curzel, Serena↗

Understanding the impacts of inorganic species in woody biomass for preprocessing and pyrolysis–A review

Woody biomass represents an abundant resource for sustainable biofuels, biochemicals, and bioproducts. Technologies for converting woody biomass have been established for decades, and research consistently highlights the critical role of inorganic species and ash plays in feedstock handling and conversion processes, including equipment plugging, corrosion, and catalyst deactivation. A thorough understanding of the variability, transport behavior, and downstream impact of inorganic species in woody biomass is essential for defining feedstock quality specifications and developing effective management strategies for conversion processes. This review compiles critical information in five main sections: 1) inorganic species concentration in woody biomass, based on anatomical fractions and their sources of variability; 2) technique features for quantifying inorganic elemental chemical analysis; 3) impacts of inorganic species on biomass preprocessing; 4) impacts of inorganic species on pyrolysis, and 5) mitigation strategies. Additionally, this review explores future challenges and opportunities in addressing the impacts of inorganic species on biomass quality. These insights aim to support the sustainable development of the biomass-to-bioenergy pipeline and ensure high-quality lignocellulosic feedstocks for efficient downstream conversions. The findings offer valuable guidance to policy makers, industry stakeholders, and researchers in developing effective strategies for managing inorganic species in woody biomass and fostering the sustainable processes for lignocellulosic biorefineries.

09 BIOMASS FUELS↗

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗

Coulomb excitation of 102 Ru with 12 C and 16 O

Here, the Coulomb excitation of 102 Ru was performed with beams of 12 C and 16 O ions. The beam particles scattered at forward angles were momentum analyzed with a magnetic spectrograph. The resolution achieved enabled the populations of the $2^{+}_{1}$ state, the unresolved $2^{+}_{2}$/ $4^{+}_{1}$, and $2^{+}_{4}$/$3^{-}_{1}$, doublets of states, and the $3^{-}_{2}$ state to be determined as a function of the scattering angle. These populations are compared with gosia calculations, yielding B(E2; $2^{+}_{1}$→$0^{+}_{1}$) = 41.5 ± 2.3 W.u., B(E2; $2^{+}_{2}$→$0^{+}_{1}$) = 1.75 ± 0.11 W.u., B(E3; $3^{-}_{1}$→$0^{+}_{1}$) = 31.5 ± 3.5 W.u., and B(E3; $3^{-}_{2}$→$0^{+}_{1}$) = 6.8 ± 0.5 W.u. The B(E3; $3^{-}_{1}$→$0^{+}_{1}$) value is significantly larger than previously measured. The weakly populated $2^{+}_{3}$ state, presumed to be a member of the band built on the $0^{+}_{2}$ state, was observed clearly for a single angle only, and a fit to its population results in B(E2; $2^{+}_{3}$→$0^{+}_{1}$) = 0.053 ± 0.011 W.u. Using the known γ -ray branching ratios for the $2^{+}_{3}$ level, the B(E2; $2^{+}_{3}$→$0^{+}_{2}$) value is calculated to be 18 ± 4 W.u., substantially less than the B(E2; $2^{+}_{1}$→$0^{+}_{1}$) . This suggests that the deformation of the $0^{+}_{2}$ state is lower than that of the $0^{+}_{1}$ state. The results are compared with beyond-mean-field calculations with the Gogny-D1S interaction using the symmetry-conserving configuration-mixing method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗