Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simultaneous approximation terms”

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.

29 records · Page 2

Scalable High-H 2 Flux, Robust Thin Film Solid Oxide Electrolyzer

This project was aimed at the development of proton-conducting SOEC (P-SOEC) technology that has the potential to meet key DOE H 2 production targets. A decreased proton resistance of the electrolyte and Faradaic efficiency improvements were sought to increase the fraction of consumed electrolysis power that is used to actually generate H 2 , while simultaneously decreasing cost dramatically. Moreover, the development project was intended to yield maximum durability through the use of a steam protective layer. Furthermore, sputtering was used to overcome processing challenges that have hampered P-SOEC development, while the low-temperature operation goal of 500 °C was expected to aid in mitigating thermally activated long-term degradation. The approach to high-performance, lower-temperature SOECs leveraged our existing SOFC Ni-cermet anode support and extensive thin-film sputtering layer-deposition experience. Rather than an all-in-one, reversible fuel cell approach which has many unacceptable tradeoffs, we focused on the many benefits to hydrogen generating SOECs, including the existence of synergies for reduced manufacturing costs (e.g., SOECs and SOFCs share supporting layers and overall manufacturing processing). The end result of this project was expected to increase current performance at 500 °C from approximately 0.8 A/cm 2 (at 60% Faradaic efficiency) at 1.4 V to > 1 A/cm 2 (at > 95% Faradaic efficiency) with a > 40% reduction in system cost and to enable operation of P-SOECs in steam contents >> 20% for a goal of a > 40,000 hours lifetime. To enable 500 °C operation in a very high steam atmosphere (> 20%), we proposed the use of a sputtered dense thin film (~0.1-1 µm thick) of high-stability Ba(Zr,Y)O 3 (BZY) to protect the Ba(Ce,Zr,Y,Yb)O 3 (BCZYYb) electrolyte. The BZCYYb, in turn, blocks the hole conductivity of the BZY to boost Faradaic efficiency. As FE increases, more of the consumed electricity is used in electrolysis to generate H 2 , rather than being shunted. Additionally, as cell resistance decreases, the voltage required to maintain current decreases, as well as the power required to generate the same amount of H 2 . With the proposed enhancements, these two factors result in the final 46% decrease in power needed to run the system. Likewise, a production rate of 50,000 kg H 2 /day will require 55% less active area, such that a system will need only 650 cells for an 80 cm 2 active area instead of ~1,440. Taking the 2016 DOE projected current cost and modifying the electricity cost and linearly scaling the other costs (except thermal feedstock) based on the cell area improvement, results in a 44% decrease in lifetime system cost, or a decrease from $\$$4.95/kg H 2 to $\$$2.75/kg H 2 . This is well below the 2018 DOE target of $\$$4/kg H 2 . The results from this project showed that we can create a P-SOEC with enhanced steam stability using two different electrolytes (i.e., one on top of the other) and achieve sufficiently low area specific resistance (ASR) to achieve the target performance. Unfortunately, due to extended delays at the beginning of the project and related supply chain and equipment access issues, we were not able to completely show increased Faradaic efficiency for the P-SOEC and therefore were unable to demonstrate the full proof of concept within the first budget period budget. While there are still challenges that remain to be solved, significant progress was made during this project and the concept still has merit that warrants further development.

08 HYDROGEN↗

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance Validation of Endcap Timing Layer Detector Modules for CMS Phase-2 Upgrade

The High-Luminosity Large Hadron Collider (HL-LHC) will operate up to 200 simultaneous collisions per bunch crossing, which is a significant jump from the current value of about 30 collisions per bunch crossing, producing significant pileup that challenges accurate event reconstruction. To address this, the Compact Muon Solenoid (CMS) experiment is implementing the Endcap Timing Layer (ETL), a precision timing detector designed to provide timing measurements with a resolution of approximately 50 ps per hit and 35 ps per reconstructed track. This project focuses on the characterization and validation of ETL detector modules using laboratory laser test data. ETL performance is evaluated by analyzing key quantities such as Bunch Crossing Identification (BCID), Time-of-Arrival (TOA), Time-over-Threshold (TOT), Noise Stability etc. The resulting analysis provides insight into the timing performance and operational stability of ETL modules, contributing to the quality assurance process req uired before their installation in the CMS experiment and ensure that the modules can meet the long-term reliability requirements of approximately 10 years of operation without replacement.

Chipara, Munashe [Fermilab]↗

A multi-product landscape life-cycle assessment approach for evaluating local climate mitigation potential

Increasing demand for land-based climate mitigation requires more efficient management of agricultural landscapes for competing objectives. Here we develop methods for assessing trade-offs and synergies between intensification and carbon-sequestering conservation measures in annual crop production landscapes using the DayCent ecosystem model and the Greenhouse Gases, Regulated Emissions, and Energy Use in Technologies (GREET) life-cycle assessment (LCA) model. We compiled county-scaled crop yields, fertilizer application rates, and tillage intensity for a corn–soybean farming case study landscape in the US state of Iowa. Using DayCent, we estimated a baseline soil organic carbon (SOC) accrual rate of 0.29 Mg C ha -1 y -1 driven by historical increases in crop productivity and reductions in tillage intensity. We then simulated the effects of management interventions targeted toward intensification (stover removal) and SOC sequestration (tillage intensity reduction and winter cover crop addition) individually and in combination. We propose a new multi-product landscape–LCA approach that analyzes marginal changes in corn grain, corn stover, and soybean production from the landscape in terms of their value for biofuel production (corn ethanol, soy biodiesel, and cellulosic ethanol from stover) and associated net displacement of conventional fossil-derived fuel use. This enables us to evaluate both intensification and sequestration effects in common CO 2 -equivalent mitigation units. We also used DayCent-simulated yields under the different land management scenarios to estimate farm-level costs and revenues. Our results show that intensification via collecting 30% of corn stover for biofuel production would increase the total greenhouse gas (GHG) mitigation potential of this landscape by 0.93 Mg CO 2 e ha -1 y -1 and provide $49 ha -1 y -1 of additional net revenue from biomass sales, but would reduce the baseline SOC accumulation rate by approximately 40%. In contrast, integrated approaches that include co-adoption of winter cover cropping and/or tillage intensity reduction would result in increased rates of SOC accumulation above the baseline, achieving simultaneous improvements in both farm profits and the overall GHG mitigation potential of the landscape.

54 ENVIRONMENTAL SCIENCES↗

Decoupling indentation size and strain rate effects during nanoindentation: A case study in tungsten

Materials indented at small scales may simultaneously exhibit indentation size and strain rate effects which complicate the identification of the mechanisms that control deformation and strength. Here, in this work, we explore the possibility that indentation size and rate effects in some materials can be decoupled in a simple way. Nanoindentation tests with various load-time histories were carried out to measure the hardness of a tungsten single crystal over a wide range of indentation depths (~500–3600 nm) and indentation strain rates (~5 ×10 –5 –2 ×10 –1 s –1 ). Under these conditions, this material exhibits significant indentation size and rate effects, but the size effect is, to a good approximation, independent of strain rate. It is shown that this behavior can be understood by the Nix-Gao model for the indentation size effect modified to include the effects of a strain rate dependent friction stress. As a consequence, the size and rate dependencies of the hardness can be expressed as the sum of two independent terms: H($\dot{ε}$ i , h c )=H f ($\dot{ε}$ i ) + (H 0 – H f ) $\sqrt{1+ \frac{h^*}{h_c}}$, where H($\dot{ε}$ i , h c ) is the hardness at given indentation strain rate ($\dot{ε}$ i ) and contact depth (h c ), H f ($\dot{ε}$ i ) is the hardness contributed by the rate dependent friction stress, and (H 0 – H f ) and h* are size and rate independent constants that follow from the Nix-Gao analysis. This formula, together with an expression for the rate dependence of H f , was successfully applied to decouple the indentation size and rate effects observed in tungsten. In addition, the physics underlying the rate independence of indentation size effect is discussed, which provides guidance for application of the proposed approach to other materials.

42 ENGINEERING↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

A Unification Framework for Euclidean and Hyperbolic Graph Neural Networks

Hyperbolic neural networks have recently gained significant attention due to their promising results on several graph problems including node classification and link prediction. The primary reason for this success is the effectiveness of hyperbolic space in capturing the inherent hierarchy of graph datasets. However, they are limited in terms of generalization, scalability, and have inferior performance when applied to non-hierarchical datasets. In this paper, we take a completely different perspective for modeling hyperbolic networks and answer the following question: is an Euclidean model able to approximate a function or behavior in the hyperbolic space? Extending the universal approximation theory developed for Euclidean models, We draw an analogy from the hyperbolic components to the Euclidean counterparts and conclude that, in order to capture hierarchical features, it is possible to generalize hyperbolic models to be a special case of Euclidean models with the proposed Pseudo-Poincaré technique. We applied our non-linear hyperbolic normalization to the current state-of-the-art homogeneous and multi-relational graph networks and demonstrate significant improvements in performance compared to both Euclidean and hyperbolic counterparts. The primary impact of this work lies in its ability to capture hierarchical features in the Euclidean space, and thus, can replace hyperbolic networks without any loss in performance metrics while simultaneously leveraging the power of Euclidean networks such as interpretability and efficient execution of various model components.

Khatir, Mehrdad↗

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Land of Opportunity: Potential for Renewable Energy on Federal Lands

Renewable energy (RE) in the United States has historically been deployed primarily on private lands, but the growing interest in RE development, generally across the country and specifically on federal lands, raises questions about the potential for RE on public lands. This study seeks to estimate RE technical potential on federal lands and project the amount of RE to be developed on federal lands under decarbonization scenarios for the contiguous United States. The study applied a combination of high-resolution geospatial analysis and power sector modeling and relied on multiple partner federal agencies and land administrators from the Bureau of Land Management, the U.S. Fish and Wildlife Service, the U.S. Forest Service, the U.S. Department of Defense, and the U.S. Department of Energy. In our reference siting access case, we estimate 44 million acres of federal land across the contiguous United States is potentially suitable for UPV development, which corresponds to a capacity technical potential of 5,750 gigawatts (GW) (Figure ES-1). Federal land area available for wind development is similar to UPV but wind’s generating capacity technical potential is lower (875 GW). The technical potential is estimated for two geothermal technologies, hydrothermal and enhanced geothermal systems (EGS), both of which have a smaller amount of federal land available for development (12 and 27 million acres, respectively). However, in terms of capacity, the technical potential for EGS (975 GW) is approximately equal to wind’s technical potential and there is an estimated 130 GW of hydrothermal potential. We also developed cases with more-limited land available for UPV and wind (for federal and non-federal lands) resulting in 96% reduction of wind capacity potential and 70% reduction for UPV. In a case with additional constraints applied to non-federal lands only (Limited Private), the overall (federal and non-federal) technical potential declines but the share of that technical potential on federal lands is higher than in the other siting cases. The technical potential is the maximum amount that could be developed, but only a small fraction would be developed or needed in the future. Across seven scenarios that achieve 100% carbon-free electricity by 2035, we estimate 26 GW to 270 GW of RE capacity could be deployed on federal lands by 2035. The three central scenarios have 51–84 GW of RE deployed by 2035 and requiring about 500,000 to 1,000,000 acres of total land area. Direct land consumption is less than total land use required, thus enabling co-use opportunities. The large technical potential estimates and the increasing deployment projections from the collection of scenarios show the opportunities for RE development on federal lands. Capturing these opportunities-while minimizing conflicts with other land uses, federal department or agency missions, and public interest, and simultaneously maximizing the economic, grid, and social value of the projects-would require collaborative planning among federal land administrators, grid planners, project developers, the public, and other stakeholders.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗