Search NASA⌕ Search

SEARCH · Search NASA

Results for “Level set method”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

High-level theoretical characterization of the vinoxy radical ( • CH 2 CHO) + O 2 reaction

Numerous processes in atmospheric and combustion chemistry produce the vinoxy radical ( • CH 2 CHO). To understand the fate of this radical and to provide reliable energies needed for kinetic modeling of such processes, we have examined its reaction with O 2 using highly reliable theoretical methods. Utilizing the focal point approach, the energetics of this reaction and subsequent reactions were obtained using coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set limit. These extrapolated energies were appended with several corrections including a treatment of full triples and connected quadruple excitations, i.e., CCSDT(Q). In addition, this study models the initial vinoxy radical + O 2 reaction for the first time with multireference methods. We predict a barrier for this reaction of approximately 0.4 kcal mol –1 . This result agrees with experimental findings but is in disagreement with previous theoretical studies. The vinoxy radical + O 2 reaction produces a 2-oxoethylperoxy radical which can undergo a number of unimolecular reactions. Abstraction of a β-hydrogen (a 1,4-hydrogen shift) and dissociation back to reactants are predicted to be competitive to each other due to their similar barriers of 21.2 and 22.3 kcal mol –1 , respectively. The minimum-energy β-hydrogen abstraction pathway produces a hydroperoxy radical (QOOH) that eventually decomposes to formaldehyde, CO, and • OH. Two other unimolecular reactions of the peroxy radical are α-hydrogen abstraction (38.7 kcal mol –1 barrier) and HO 2 • elimination (43.5 kcal mol –1 barrier). Furthermore, these pathways lead to glyoxal + • OH and ketene + HO 2 • formation, respectively, but they are expected to be uncompetitive due to their high barriers.

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

The reaction between the bromine atom and the water trimer: high level theoretical studies

Three different reaction pathways are found for the reaction of bromine atom (Br) with the lowest-energy structure of the water trimer [uud-(H 2 O) 3 ], initially using the MPW1K-DFT method. The three bromine pathways have closely related geometries and energetics, analogous to those found for the fluorine and chlorine reactions. Here, the lowest-energy pathway of the Br + uud-(H 2 O) 3 reaction was further investigated using the “gold standard” CCSD(T) method and the correlation-consistent basis sets up to cc-pVQZ(-PP). Based on the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the Br + (H 2 O) 3 reaction is endothermic by 33.3 kcal mol -1 . The classical barrier height is 29.0 kcal mol -1 between the reactants and the exit complex, and there is no barrier for the reverse reaction. The Br···(H 2 O) 3 entrance complex is found to lie 4.7 kcal mol -1 below the separated reactants, and the HBr···(H 2 O) 2 OH exit complex is bound by 6.4 kcal mol -1 relative to the separated products. This potential energy profile is further corrected by the zero point energies and spin–orbit coupling effects. Structurally, the Br + (H 2 O) 3 stationary points can be derived from those of the simpler Br + (H 2 O) 2 reaction by judiciously appending a H 2 O molecule. The Br + (H 2 O) 3 potential energy profile is compared with the Br + (H 2 O) 2 and Br + H 2 O reactions, as well as to the valence isoelectronic Cl + (H 2 O) 3 and F + (H 2 O) 3 systems. It is reasonable to expect that the reactions between the bromine atom and larger water clusters would be similar to the Br + (H 2 O) 3 reaction.

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

A scalable stochastic programming approach for the design of flexible systems

We study the problem of designing systems in order to minimize cost while meeting a given flexibility target. Flexibility is attained by enforcing a joint chance constraint, which ensures that the system will exhibit feasible operation with a given target probability level. Unfortunately, joint chance constraints are complicated mathematical objects that often need to be reformulated using mixed-integer programming (MIP) techniques. In this work, we cast the design problem as a conflict resolution problem that seeks to minimize cost while maximizing flexibility. We propose a purely continuous relaxation of this problem that provides a significantly more scalable approach relative to MIP methods and show that the formulation delivers solutions that closely approximate the Pareto set of the original joint chance-constrained problem.

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

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2019

Assessing Electronic Structure Methods for Long-Range Three-Body Dispersion Interactions: Analysis and Calculations on Well-Separated Metal Atom Trimers

Three-body dispersion interactions are much weaker than their two-body counterpart. However, their importance grows quickly as the number of interacting monomers rises. To explore the numerical performance of correlation methods for long-range three-body dispersion, we performed calculations on eight very simple dispersion-dominated model metal trimers: Na 3 , Mg 3 , Zn 3 , Cd 3 , Hg 3 , Cu 3 , Ag 3 , and Au 3 . One encouraging aspect is that relatively small basis sets of augmented triple-ζ size appear to be adequate for three-body dispersion in the long-range. Coupled cluster calculations were performed at high levels to assess MP3, CCSD, CCSD(T), empirical density functional theory dispersion (D3), and the many-body dispersion (MBD) approach. We found that the accuracy of CCSD(T) was generally significantly lower than for two-body interactions, with errors sometimes reaching 20% in the investigated systems, while CCSD and particularly MP3 were generally more erratic. MBD is found to perform better than D3 at large distances, whereas the opposite is true at shorter distances. When computing reference numbers for three-body dispersion, care should be taken to appropriately represent the effect of the connected triple excitations, which are significant in most cases and incompletely approximated by CCSD(T).

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

Curated Collection of More than 20,000 Experimentally Reported One-Dimensional Metal–Organic Frameworks

Herein, a collection of more than 20,000 experimentally derived crystal structures for metal–organic frameworks (MOFs) that do not have two- or three-dimensional covalently bonded networks has been developed from the materials available at the Cambridge Crystallographic Data Centre. Of these 20,000 1D MOFs, more than 12,000 structures have been verified to be solvent-free and in exact agreement with the stoichiometry of the synthesized materials. More than 10% of the complete data set comprise materials including two or more distinct metals. The band gaps of more than 12,000 1D MOFs have been computed at the density functional theory-generalized gradient approximation level, finding more than 2000 materials that have a zero band gap. Molecular simulations of CH4 adsorption in a small number of 1D MOFs indicated that adsorbate-induced deformation plays a significant role in determining adsorption isotherms in these materials. As a result, methods that have been used previously for high-throughput predictions of molecular adsorption in 3D MOFs are not suitable for 1D MOFs.

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

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗
DOE OSTI2022

Application of the ASME Boiler and Pressure Vessel Code in the Design of SRF Cavities at Fermilab

Jacketed Superconducting Radio Frequency (SRF) cavities structurally comprise of an inner niobium vessel surrounded by a liquid helium containment vessels. The pressure of the helium bath and/or its volume might be such that a jacketed SRF cavity shall be considered a system of pressure vessels. Thus, methods described in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC) should be used to analyze the structural soundness of jacketed SRF cavities. This paper will report the use of the set of rules developed at Fermilab for the design of SRF cavities, such as jacketed 1.3 GHz cavities for LCLS-II HE and jacketed Single Spoke Resonator type~2 (SSR2) for PIP-II, to ensure a similar level of safety as prescribed by the ASME BPVC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Evaluating Equivalent Monte Carlo Calculations

We address several basic issues, primarily whether two sets of MCNP tallies are stochastically equivalent (i.e., that the sets of particle history scores follow the same statistical distribution) and, if di erent, quantify the signi cance level of the di erence. The underlying topics here are of longstanding interest to the statistical community, and no \new" theoretical work needs to be carried out to address them. Related comparisons are useful to code developers in assessing code modi cations and to code users in assessing other simulations and related physical measurements. However, these statistical methods must be applied with care to avoid obtaining misleading results, as we show in several examples. The text is written for MCNP practitioners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Effect of the Flux Separability Approximation on Multigroup Neutron Transport

The angular dependence of flux-weighted multigroup cross sections is commonly neglected when generating multigroup libraries. The error of this flux separability approximation is typically not isolated from other error sources due to a lack of availability of library generation and corresponding solvers that cannot relax this approximation. These errors can now be isolated and quantified with the availability of a multigroup Monte Carlo transport and multigroup library-generation capability in the OpenMC Monte Carlo transport code. This work will discuss relevant details of the OpenMC implementation, provide an example case useful for detailing the type of errors one can expect from making the flux separability approximation, and end with more realistic problems which show the impact of the approximation and highlight how it can strongly arise from an energy-dependent resonance absorption effect. Since the angle-dependence is intrinsically linked to the energy group structure, these examples also show that relaxing the flux separability approximation with angle-dependent cross sections could be used to reduce either the fine-tuning required to set a multigroup energy structure for a specific reactor type or the number of energy groups required to obtain a desired level of accuracy for a given problem. This trade-off could increase the costs of generating multigroup cross sections, and has the potential to require more memory for storing the multigroup library during the transport calculations, but it can significantly reduce the computational time required since the runtime of a discrete ordinates or method of characteristics neutron transport solver scales roughly linearly with the number of groups.

42 ENGINEERING↗
DOE OSTI2025

An MLCommons Scientific Benchmarks Ontology

Scientific machine learning research spans diverse domains and data modalities, yet existing benchmark efforts remain siloed and lack standardization. This makes novel and transformative applications of machine learning to critical scientific use-cases more fragmented and less clear in pathways to impact. This paper introduces an ontology for scientific benchmarking developed through a unified, community-driven effort that extends the MLCommons ecosystem to cover physics, chemistry, materials science, biology, climate science, and more. Building on prior initiatives such as XAI-BENCH, FastML Science Benchmarks, PDEBench, and the SciMLBench framework, our effort consolidates a large set of disparate benchmarks and frameworks into a single taxonomy of scientific, application, and system-level benchmarks. New benchmarks can be added through an open submission workflow coordinated by the MLCommons Science Working Group and evaluated against a six-category rating rubric that promotes and identifies high-quality benchmarks, enabling stakeholders to select benchmarks that meet their specific needs. The architecture is extensible, supporting future scientific and AI/ML motifs, and we discuss methods for identifying emerging computing patterns for unique scientific workloads. The MLCommons Science Benchmarks Ontology provides a standardized, scalable foundation for reproducible, cross-domain benchmarking in scientific machine learning. A companion webpage for this work has also been developed as the effort evolves: https://mlcommons-science.github.io/benchmark/

Hawks, Ben [Fermilab] (ORCID:0000000157000288)↗
DOE OSTI2023

Spatiotemporal Super-Resolution with Generative Machine Learning for Creating Renewable Energy Resource Data Under Climate Change Scenarios

As we plan for a future with higher penetrations of renewables and increasing electrification, it becomes more important to understand how the electricity grid will operate under a variety of weather events. We must also consider that the weather our future grid will experience will be different and possibly more extreme than the historical weather that we have extensive data for. We can use data from global climate models (GCMs) to help understand how our climate may change over the next several decades, but there is often a significant gap between the low-resolution GCM data and the high-resolution weather data required to study power systems under specific weather events. Therefore, our objective in this work is to develop tools that can bridge this gap by using low-resolution GCM data to create realistic high-resolution weather datasets that can be used to study renewable energy generation and electricity demand. To accomplish this objective, we have developed a set of generative machine learning models that can rapidly downscale GCM daily average output data at an approximate grid resolution of 100km to hourly data at an approximate 4 km grid resolution. The models can be used to create high resolution data from nearly any GCM included in the Coupled Model Intercomparison Project (CMIP) Phase 5 or 6. Our methods include all datasets regularly used to study the integration of wind and solar power plants as well as changes in electricity demand due to heating and cooling loads. These models and datasets enable power systems modelers to study climate change-influenced weather events and their impact on the grid. We have downscaled and validated wind, solar, temperature, and humidity data with very promising results. The generative machine learning methods are computationally efficient and produce data that has similar statistical characteristics to current state-of-the-art historical datasets. We have trained initial generative models and produced an initial dataset collectively referred to as Sup3rCC: Super-Resolved Renewable Energy Resource Data with Climate Change Impacts. The data covers a (mostly) historical period from 2015-2025 and a future period from 2050-2059. We have also taken hypothetical high-electrification load data and scaled the heating and cooling loads with respect to the 2050-2059 high-resolution Sup3rCC meteorology. The results show how future levels of renewable energy generation and electrified load may be impacted by climate change, setting the stage for capacity expansion models to consider a dynamic climate through model years.

climate change↗
DOE OSTI2025

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗
DOE OSTI2017

Automating cell detection and classification in human brain fluorescent microscopy images using dictionary learning and sparse coding

Background Immunofluorescence (IF) plays a major role in quantifying protein expression in situ and understanding cell function. It is widely applied in assessing disease mechanisms and in drug discovery research. Automation of IF analysis can transform studies using experimental cell models. However, IF analysis of postmortem human tissue relies mostly on manual interaction, often subjected to low-throughput and prone to error, leading to low inter and intra-observer reproducibility. Human postmortem brain samples challenges neuroscientists because of the high level of autofluorescence caused by accumulation of lipofuscin pigment during aging, hindering systematic analyses. We propose a method for automating cell counting and classification in IF microscopy of human postmortem brains. Our algorithm speeds up the quantification task while improving reproducibility. New method Dictionary learning and sparse coding allow for constructing improved cell representations using IF images. These models are input for detection and segmentation methods. Classification occurs by means of color distances between cells and a learned set. Results Our method successfully detected and classified cells in 49 human brain images. We evaluated our results regarding true positive, false positive, false negative, precision, recall, false positive rate and F1 score metrics. We also measured user-experience and time saved compared to manual countings. Comparison with existing methods We compared our results to four open-access IF-based cell-counting tools available in the literature. Lastly, our method showed improved accuracy for all data samples. Conclusion The proposed method satisfactorily detects and classifies cells from human postmortem brain IF images, with potential to be generalized for applications in other counting tasks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The entropy of galaxy spectra: how much information is encoded?

Abstract The inverse problem of extracting the stellar population content of galaxy spectra is analysed here from a basic standpoint based on information theory. By interpreting spectra as probability distribution functions, we find that galaxy spectra have high entropy, thus leading to a rather low effective information content. The highest variation in entropy is unsurprisingly found in regions that have been well studied for decades with the conventional approach. We target a set of six spectral regions that show the highest variation in entropy – the 4000 Å break being the most informative one. As a test case with real data, we measure the entropy of a set of high-quality spectra from the Sloan Digital Sky Survey, and contrast entropy-based results with the traditional method based on line strengths. The data are classified into star-forming (SF), quiescent (Q), and active galactic nucleus (AGN) galaxies, and show – independently of any physical model – that AGN spectra can be interpreted as a transition between SF and Q galaxies, with SF galaxies featuring a more diverse variation in entropy. The high level of entanglement complicates the determination of population parameters in a robust, unbiased way, and affects traditional methods that compare models with observations, as well as machine learning (especially deep learning) algorithms that rely on the statistical properties of the data to assess the variations among spectra. Entropy provides a new avenue to improve population synthesis models so that they give a more faithful representation of real galaxy spectra.

Ferreras, Ignacio (ORCID:0000000345843127)↗
DOE OSTI2022

Joint Determination of Reactor Antineutrino Spectra from U 235 and Pu 239 Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of 235U and 239Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of 235U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant 235U and 239Pu isotopes and improves the uncertainty of the 235U spectral shape to about 3%. The 235U and 239Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the 235U and 239Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Heat Exchanger Design Optimization for Mitigation of Diesel Engine Exhaust Fouling

Abstract The design optimization of a diesel exhaust-coupled heat and mass exchanger that drives a 2.71 kW cooling capacity absorption heat pump is presented in this study. Fouling layer thermal resistance and pressure drops from single-tube experiments are used to develop a thermodynamic, heat transfer, and pressure drop model for the exhaust-coupled desorber. A parametric study is performed to select a desorber design that meets system performance while minimizing footprint. Experimental heat duties and pressure drops are within 10% and 3%, respectively, of the model predictions. Thus, large data sets from single-tube experiments with representative geometries are successful in accounting for fouling effects at the component level. Desorber design optimization based on this approach ensures continued heat pump performance after fouling. This study, along with the single-tube experiments, presents a systematic approach to design exhaust-coupled heat exchangers while considering the effects of fouling. These results are applicable for a wide range of waste-heat recovery applications, and this method can be extended to different geometries and operating conditions.

Energy & Fuels↗
DOE OSTI2018

VisIRR: A Visual Analytics System for Information Retrieval and Recommendation for Large-Scale Document Data

In this paper, we present an interactive visual information retrieval and recommendation system, called VisIRR, for large-scale document discovery. VisIRR effectively combines the paradigms of (1) a passive pull through query processes for retrieval and (2) an active push that recommends items of potential interest to users based on their preferences. Equipped with an efficient dynamic query interface against a large-scale corpus, VisIRR organizes the retrieved documents into high-level topics and visualizes them in a 2D space, representing the relationships among the topics along with their keyword summary. In addition, based on interactive personalized preference feedback with regard to documents, VisIRR provides document recommendations from the entire corpus, which are beyond the retrieved sets. Such recommended documents are visualized in the same space as the retrieved documents, so that users can seamlessly analyze both existing and newly recommended ones. This article presents novel computational methods, which make these integrated representations and fast interactions possible for a large-scale document corpus. We illustrate how the system works by providing detailed usage scenarios. Finally, we present preliminary user study results for evaluating the effectiveness of the system.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗