SEARCH · Engineering Papers
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.
We are making available a multi-petabyte data set for open science via a data hackathon. The participants bring their science to the computing and data facility, alleviating the need to transfer and manage the large volume of data. This accelerates the process of discovery science while increasing scientific throughput.
Abstract Experimental measurements in deep-inelastic scattering and lepton-pair production on deuterium targets play an important role in the flavor separation of u and d (anti)quarks in global QCD analyses of the parton distribution functions (PDFs) of the nucleon. We investigate the impact of theoretical corrections accounting for the light-nuclear structure of the deuteron upon the fitted u , d -quark, gluon, and other PDFs in the CJ15 and CT18 families of next-to-leading order CTEQ global analyses. The investigation is done using the $$L_2$$ L 2 sensitivity statistical method, which provides a common metric to quantify the strength of experimental constraints on various PDFs and ratios of PDFs in the two distinct fitting frameworks. Using the $$L_2$$ L 2 sensitivity and other approaches, we examine the compatibility of deuteron data sets with other fitted experiments under varied implementations of the deuteron corrections. We find that freely-fitted deuteron corrections modify the PDF uncertainty at large momentum fractions and will be relevant for future PDFs affecting electroweak precision measurements.
The large data rates at the LHC require an online trigger system to select relevant collisions. Rather than compressing individual events, we propose to compress an entire data set at once. We use a normalizing flow as a deep generative model to learn the probability density of the data online. The events are then represented by the generative neural network and can be inspected offline for anomalies or used for other analysis purposes. We demonstrate our new approach for a toy model and a correlation-enhanced bump hunt.
Deep learning (DL) models are being deployed at medical centers to aid radiologists for diagnosis of lung conditions from chest radiographs. Such models are often trained on a large volume of publicly available labeled radiographs. These pre-trained DL models’ ability to generalize in clinical settings is poor because of the changes in data distributions between publicly available and privately held radiographs. In chest radiographs, the heterogeneity in distributions arises from the diverse conditions in X-ray equipment and their configurations used for generating the images. In the machine learning community, the challenges posed by the heterogeneity in the data generation source is known as domain shift, which is a mode shift in the generative model. In this work, we introduce a domain-shift detection and removal method to overcome this problem. Our experimental results show the proposed method’s effectiveness in deploying a pre-trained DL model for abnormality detection in chest radiographs in a clinical setting.
SAND2023-06610O The Nuclear Incident Response Program (NIRP) Core Software Suite is a set of code that supports multiple applications. It includes miscellaneous base code for data objects, mathematic equations, and user interface components that are used for developing large applications, including those used for data management. Among the types of applications: • Analyst Manager is used for managing contact information, which is also often included in reports and identifies the owners of calculations. • Radionuclide Viewer provides access to the Dose Coefficient File Package (DCFPAK) radiological data, and it complements the Mixture Manager tool. • Mixture Manager creates and manages radionuclides mixtures that are commonly used in other applications. • High Explosive Manager manages explosives and properties. • Chart Viewer is used for data charts such as meteorology charts and other charts specific to data needs. • Unit Converter converts numeric values to/from different units. This package includes object-oriented base code such as data objects describing the date, time, and location of an event; weather conditions; data types that support predefined ranges of values; a unit library, and others. The package also includes mathematic equations such as determining the distance between two points; erf functions; factorial functions; incidence fractions of probit curves; and others. Other features are customized user interface components for text fields; unit selectors (combo boxes); menus; panels; bug-tracking; and automated checks for software updates. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
ABSTRACT We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales ($\lt 30\, h^{-1}$ Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the ‘lensing is low’ tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets.
Here we show how the galaxy four-point correlation function can test for cosmological parity violation. The detection of cosmological parity violation would reflect previously unknown forces present at the earliest moments of the Universe. Recent developments both in rapidly evaluating galaxy N-point correlation functions and in determining the corresponding covariance matrices make the search for parity violation in the four-point correlation function possible in current and upcoming surveys such as those undertaken by Dark Energy Spectroscopic Instrument, the Euclid satellite, and the Vera C. Rubin Observatory. We estimate the limits on cosmic parity violation that could be set with these data.
Abstract We report on the design and performance of the B icep3 instrument and its first three-year data set collected from 2016 to 2018. B icep3 is a 52 cm aperture refracting telescope designed to observe the polarization of the cosmic microwave background (CMB) on degree angular scales at 95 GHz. It started science observation at the South Pole in 2016 with 2400 antenna-coupled transition-edge sensor bolometers. The receiver first demonstrated new technologies such as large-diameter alumina optics, Zotefoam infrared filters, and flux-activated SQUIDs, allowing ∼10× higher optical throughput compared to the Keck design. B icep3 achieved instrument noise equivalent temperatures of 9.2, 6.8, and 7.1 μ K CMB s and reached Stokes Q and U map depths of 5.9, 4.4, and 4.4 μ K arcmin in 2016, 2017, and 2018, respectively. The combined three-year data set achieved a polarization map depth of 2.8 μ K arcmin over an effective area of 585 square degrees, which is the deepest CMB polarization map made to date at 95 GHz.
A search is performed for a heavy scalar resonance X decaying to a Higgs boson (H) and a Higgs-like scalar boson (Y) in the two bottom quark (H → $b\bar{b}$) and four quark (Y → VV → 4q) final state, where V denotes a W or Z boson. Masses of the X between 900 and 4000 GeV and the Y between 60 and 2800 GeV are considered. The search is performed in data collected by the CMS experiment at the CERN LHC from proton-proton collisions at 13 TeV center-of-mass energy, with a data set corresponding to a total integrated luminosity of 138 fb −1 . It targets the Lorentz-boosted regime, in which the products of the H → $b\bar{b}$ decay can be reconstructed as a single large-area jet, and those from the Y → VV → 4q decay as either one Y → 4q or two jets V → H → $q\bar{q}$. Jet identification and mass reconstruction exploit machine-learning tools, including a novel attention-based “particle transformer” for Y → 4q identification. No significant excess is observed in the data above the standard model background expectation. Upper limits on the product of production cross section and branching fraction as low as 0.2 fb are derived at 95% confidence level for various mass hypotheses. This is the first search at the LHC for scalar resonances in the all-hadronic H → $b\bar{b}$VV decay channel.
This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.
Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.
Cryo-electron tomography provides detailed views of macromolecules in situ. However, imaging a large field of view to provide more cellular context requires reducing magnification during data collection, which in turn restricts the resolution. To circumvent this trade-off between field of view and resolution, we have developed a montage data collection scheme that uniformly distributes the dose throughout the specimen. In this approach, sets of slightly overlapping circular tiles are collected at high magnification and stitched to form a composite projection image at each tilt angle. These montage tilt-series are then reconstructed into massive tomograms with a small pixel size but a large field of view. For proof-of-principle, we applied this method to the thin edge of HeLa cells. Thon rings to better than 10 Å were detected in the montaged tilt-series, and diverse cellular features were observed in the resulting tomograms. These results indicate that the additional dose required by this technique is not prohibitive to performing structural analysis to intermediate resolution across a large field of view. We anticipate that montage tomography will prove particularly useful for lamellae, increase the likelihood of imaging rare cellular events, and facilitate visual proteomics.
Turbine generator power from simulations using Actuator Line Models and Actuator Disk Models with a Filtered Lifting Line Correction are compared to field data of a V27 turbine. Preliminary results of the wake characteristics are also presented. Turbine quantities of interest from traditional ALM and ADM with the Gaussian kernel (..epsilon..) set at the optimum value for matching power production and that resolve the kernel at all mesh sizes are also presented. The atmospheric boundary layer is simulated using Nalu-Wind, a Large Eddy Simulation code which is part of the ExaWind code suite. The effect of mesh resolution on quantities of interest is also examined.
We develop a method for the determination of the top quark mass using the distribution of the decay length of the B -hadrons originating from its decay. This technique is based on our earlier observation regarding the location of the peak of the b quark energy distribution. Such “energy-peak” methods enjoy a greater degree of model-independence with respect to the kinematics of top quark production compared to earlier proposals. The CMS experiment has implemented the energy-peak method using associated b -jet energy as an approximation for b quark energy. The new method uses B -hadron decay lengths, which are related to b quark energies by convolution. The advantage of the new decay length method is that it can be applied in a way that evades jet-energy scale (JES) uncertainties. Indeed, CMS has measured the top quark mass using B -hadron decay lengths, but they did not incorporate the energy-peak method. Therefore, mismodeling of top quark transverse momentum remains a large uncertainty in their result. We demonstrate that, using energy-peak methods, this systematic uncertainty can become negligible. We show that with the current LHC data set, a sub-GeV statistical uncertainty on the top quark mass can be attained with this method. To achieve a comparable systematic uncertainty as is true for many methods based on exclusive or semi-inclusive observables using hadrons, we find that the quark-hadron transition needs to be described significantly better than is the case with current fragmentation functions and hadronization models.
Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.
Adiabatic quantum computers are a promising platform for efficiently solving challenging optimization problems. Therefore, many are interested in using these computers to train computationally expensive machine learning models. We present a quantum approach to solving the balanced k-means clustering training problem on the D-Wave 2000Q adiabatic quantum computer. In order to do this, we formulate the training problem as a quadratic unconstrained binary optimization (QUBO) problem. Unlike existing classical algorithms, our QUBO formulation targets the global solution to the balanced k-means model. We test our approach on a number of small problems and observe that despite the theoretical benefits of the QUBO formulation, the clustering solution obtained by a modern quantum computer is usually inferior to the solution obtained by the best classical clustering algorithms. Nevertheless, the solutions provided by the quantum computer do exhibit some promising characteristics. We also perform a scalability study to estimate the run time of our approach on large problems using future quantum hardware. Finally, as a final proof of concept, we used the quantum approach to cluster random subsets of the Iris benchmark data set.
Soft-chemistry techniques provide a versatile approach to synthesizing inorganic materials under mild conditions, enabling access to compositions and structures that are challenging to achieve through traditional thermodynamically driven solid-state methods. However, these solution-based routes often result in phase competition, requiring precise control over reaction conditions to achieve selective product formation. While one-variable-at-a-time (OVAT) approaches have traditionally been used for phase selection, data-driven strategies are emerging as more efficient methods for navigating complex synthetic spaces. Ternary metal halides, such as cesium cadmium bromides (Cs–Cd–Br), are of growing interest due to their potential in wide and ultrawide band gap applications. Unlike the well-studied cesium lead halide phases, the compositional diversity and solution-based synthesis of ternary Cs–Cd–Br phases remain largely unexplored. This study systematically investigates the synthetic phase space of the Cs–Cd–Br system by constructing a data-driven phase map. Using a common set of precursors and a standardized experimental procedure, we successfully synthesize all four known Cs–Cd–Br phases—CsCdBr 3 , Cs 2 CdBr 4 , Cs 3 CdBr 5 , and Cs 7 Cd 3 Br 13 —each exhibiting distinct structures, morphologies, and optical properties. Our findings highlight the potential of soft-chemistry methods for expanding the library of ternary metal halides and provide key insights into the thermodynamic and kinetic factors governing phase formation.
The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.