Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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 469 records · Page 26

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗

Corrosion Testing of Refractory inContact with Molten Glasses Designed for Waste Vitrification - VSL Touchpoint Matrix Glasses

It is known that the predictive life of the refractory ceramic liner of nuclear waste glass melters is conservative, as demonstrated by performance of these materials such as in the Defense Waste Processing Facility (DWPF). The motivation for this task is to maximize the useful life of the melters that will be operated at the Waste Treatment and Immobilization Plant (WTP), which will in turn minimize procurement and disposal costs and melter outage times, as well as to identify maximum loadings in the waste glass of those species that corrode melter components. This task was initiated jointly with Pacific Northwest National Laboratory (PNNL) with the objective to develop a methodology and model to enable more accurate prediction of refractory service life under prototypic conditions from laboratory-scale material corrosion tests. Refractory corrosion is generally reported as physical material loss, measured in units of distance (e.g., inch) or as physical material loss rate, measured in units of distance per time (e.g., inch/day). In post-operational melters, the refractory corrosion is measured directly, sometimes reported as corrosion depth. Crucible tests are used in the laboratory to accelerate the refractory corrosion to facilitate a meaningful measurement in a commensurate amount of time. Crucible tests are particularly useful in understanding refractory corrosion across a large glass composition space, where operational testing would be prohibitive. Some of the critical parameters that are known to influence refractory corrosion by molten glass in a crucible test are temperature, system redox, molten salt phases, glass chemistry, and test duration. The majority of data collected for Monofrax® K-3 (hereafter referred to as K-3) corrosion is from crucible tests, but a small amount comes directly from scaled and production melters. Crucible test data has been collected under varying conditions, whereas data collected from operational melters is relatively fewer and represents conditions specific to the melter campaign. The result is that the published data can be grouped and analyzed in multiple ways, not all of which are readily comparable. The Standard Test Method for Isothermal Corrosion Resistance of Refractories to Molten Glass (ASTM C621) outlines the general guidelines used across industry. That method describes a sealed, static test in which the surface area of the refractory coupon and the volume of glass are fixed. A significant portion of the crucible data pertaining to nuclear waste glasses has been collected in a modified configuration; the most notable differences being the surface area of the refractory coupon to volume of the glass and use of a method for bubbling the melt. To our knowledge, the influence of those parameters on the refractory corrosion has not been quantified. In this work, it was determined that static tests and bubbled tests would be performed. Savannah River National Laboratory (SRNL) was tasked with setting up and performing static testing while PNNL was tasked with setting up and performing bubbled testing. Initial activities were performed to establish laboratory methods that reproduce data comparable to existing data sets of K-3 refractory corrosion by low activity waste (LAW) and high-level waste (HLW) glass compositions. Later activities were focused on refining the test parameters to establish a standard test practice to be used between Laboratories and collecting additional data to be used in the enhanced waste glass model development. This document serves primarily to convey the refractory loss measurement results from corrosion testing of K-3 refractory with waste glass compositions developed for use in the WTP melters. The data will be used in the enhanced property/composition models being developed for waste glass vitrification and melter operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Indoor Occupancy Sensing via Networked Nodes (2012–2022): A Review

In the past decade, different sensing mechanisms and algorithms have been developed to detect or estimate indoor occupancy. One of the most recent advancements is using networked sensor nodes to create a more comprehensive occupancy detection system where multiple sensors can identify human presence within more expansive areas while delivering enhanced accuracy compared to a system that relies on stand-alone sensor nodes. The present work reviews the studies from 2012 to 2022 that use networked sensor nodes to detect indoor occupancy, focusing on PIR-based sensors. Methods are compared based on pivotal ADPs that play a significant role in selecting an occupancy detection system for applications such as Health and Safety or occupant comfort. These parameters include accuracy, information requirement, maximum sensor failure and minimum observation rate, and feasible detection area. We briefly describe the overview of occupancy detection criteria used by each study and introduce a metric called “sensor node deployment density” through our analysis. This metric captures the strength of network-level data filtering and fusion algorithms found in the literature. It is hinged on the fact that a robust occupancy estimation algorithm requires a minimal number of nodes to estimate occupancy. This review only focuses on the occupancy estimation models for networked sensor nodes. It thus provides a standardized insight into networked nodes’ occupancy sensing pipelines, which employ data fusion strategies, network-level machine learning algorithms, and occupancy estimation algorithms. This review thus helps determine the suitability of the reviewed methods to a standard set of application areas by analyzing their gaps.

Emad-Ud-Din, Muhammad (ORCID:0000000279515538)↗

The structure and energetics of the HCN-HNC transition state

The optimum geometries and quadratic force constants of HCN, HNC and the transition state connecting them have been determined at the single- and double-excitation coupled-cluster (CCSD) and CCSD(T) levels of theory. Energy differences were evaluated using the CCSD and CCSD(T) methods in conjunction with large atomic natural orbital basis sets containing g-type basis functions on the heavy atoms and f-type functions on hydrogen. The most reliable structure obtained for the transition state has bond distances of 1.194, 1.188, and 1.389 A for r(CN), r(CH), and r(NH), respectively. Including a correction for zero-point vibrational energies, the transition state is predicted to be 44.6 + or - 1.0 kcal/mol above the HCN isomer, while HNC is predicted to be 14.4 + or - 1.0 kcal/mol above HCN. The latter value is in excellent agreement with the most recent experimental determination (14.8 + or - 2.0 kcal/mol).

Lee, Timothy J.↗

Integrated Design of Downwind Land-Based Wind Turbines using Analytic Gradients

Wind turbines are complex systems where component-level changes can have significant system-level effects. Effective wind turbine optimization generally requires an integrated analysis approach with a large number of design variables. Optimizing across large variable sets is orders of magnitude more efficient with gradient-based methods as compared with gradient-free method, particularly when using exact gradients. We have developed a wind turbine analysis set of over 100 components where 90% of the models provide numerically exact gradients through symbolic differentiation, automatic differentiation, and adjoint methods. This framework is applied to a specific design study focused on downwind land-based wind turbines. Downwind machines are of potential interest for large wind turbines where the blades are often constrained by the stiffness required to prevent a tower strike. The mass of these rotor blades may be reduced by utilizing a downwind configuration where the constraints on tower strike are less restrictive. The large turbines of this study range in power rating from 5-7MW and in diameter from 105m to 175m. The changes in blade mass and power production have important effects on the rest of the system, and thus the nacelle and tower systems are also optimized. For high-speed wind sites, downwind configurations do not appear advantageous. The decrease in blade mass (10%) is offset by increases in tower mass caused by the bending moment from the rotor-nacelle-assembly. For low-wind speed sites, the decrease in blade mass is more significant (25-30%) and shows potential for modest decreases in overall cost of energy (around 1-2%).

17 WIND ENERGY↗

Uncertain of uncertainties? A comparison of uncertainty quantification metrics for chemical data sets

Abstract With the increasingly more important role of machine learning (ML) models in chemical research, the need for putting a level of confidence to the model predictions naturally arises. Several methods for obtaining uncertainty estimates have been proposed in recent years but consensus on the evaluation of these have yet to be established and different studies on uncertainties generally uses different metrics to evaluate them. We compare three of the most popular validation metrics (Spearman’s rank correlation coefficient, the negative log likelihood (NLL) and the miscalibration area) to the error-based calibration introduced by Levi et al. ( Sensors 2022 , 22 , 5540). Importantly, metrics such as the negative log likelihood (NLL) and Spearman’s rank correlation coefficient bear little information in themselves. We therefore introduce reference values obtained through errors simulated directly from the uncertainty distribution. The different metrics target different properties and we show how to interpret them, but we generally find the best overall validation to be done based on the error-based calibration plot introduced by Levi et al. Finally, we illustrate the sensitivity of ranking-based methods (e.g. Spearman’s rank correlation coefficient) towards test set design by using the same toy model ferent test sets and obtaining vastly different metrics (0.05 vs. 0.65).

Rasmussen, Maria H.↗

Transport Test Problems for Radiation Detection Scenarios

Modeling forms an important component of radiation detection development, allowing for testing of new detector designs, evaluation of existing equipment against a wide variety of potential threat sources, and assessing operation performance of radiation detection systems. This can, however, result in large and complex scenarios which are time consuming to model. A variety of approaches to radiation transport modeling exist with complementary strengths and weaknesses for different problems. This variety of approaches, and the development of promising new tools (such as ORNL’s ADVANTG) which combine benefits of multiple approaches, illustrates the need for a means of evaluating or comparing different techniques for radiation detection problems. This report presents a set of 9 test problems for comparing different types of radiation transport calculations, identifying appropriate tools for classes of problems, and testing and guiding the development of new methods. Cases are drawn primarily from existing or previous calculations with a preference for cases which include experimental data, or otherwise have results with a high level of confidence, are non-sensitive, and represent problem sets of interest to NA-22. From a technical perspective, cases are chosen to span a range of difficulty and to include gamma transport, neutron transport, or both and represent different important physical processes and a range of sensitivity to angular or energy fidelity. Following test case identification, existing information about geometry, measurements, and previous calculations were assembled. Monte Carlo results (MCNP decks) were reviewed or created and re-run in order to attain accurate computational times and to verify agreement with experimental data, when present. Test case information was then conveyed to ORNL in order to guide testing and development of hybrid calculations. The results of those ADVANTG calculations were then sent to PNNL for compilation. This is a report describing the details of the selected test cases and results from available transport codes. Additional results from other transport codes will be added during the next year, along with additional cases if needed for software testing. Future work will extend the test problems to include more extensive cross-code comparisons as well as more refined metrics for evaluating the tradeoffs between operator time, computational time, and solution accuracy, to assemble a formalized benchmark set.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Production Application Performance Data Streaming for System Monitoring

In this article, we present an approach to streaming collection of application performance data. Practical application performance tuning and troubleshooting in production high-performance computing (HPC) environments requires an understanding of how applications interact with the platform, including (but not limited to) parallel programming libraries such as Message Passing Interface (MPI). Several profiling and tracing tools exist that collect heavy runtime data traces either in memory (released only at application exit) or on a file system (imposing an I/O load that may interfere with the performance being measured). Although these approaches are beneficial in development stages and post-run analysis, a systemwide and low-overhead method is required to monitor deployed applications continuously. This method must be able to collect information at both the application and system levels to yield a complete performance picture. In our approach, an application profiler collects application event counters. A sampler uses an efficient inter-process communication method to periodically extract the application counters and stream them into an infrastructure for performance data collection. We implement a tool-set based on our approach and integrate it with the Lightweight Distributed Metric Service (LDMS) system, a monitoring system used on large-scale computational platforms. LDMS provides the infrastructure to create and gather streams of performance data in a low overhead manner. We demonstrate our approach using applications implemented with MPI, as it is one of the most common standards for the development of large-scale scientific applications. We utilize our tool-set to study the impact of our approach on an open source HPC application, Nalu. Our tool-set enables us to efficiently identify patterns in the behavior of the application without source-level knowledge. We leverage LDMS to collect system-level performance data and explore the correlation between the system and application events. Finally, we demonstrate how our tool-set can help detect anomalies with a low latency. We run tests on two different architectures: a system enabled with Intel Xeon Phi and another system equipped with Intel Xeon processor. Our overhead study shows our method imposes at most 0.5% CPU usage overhead on the application in realistic deployment scenarios.

59 BASIC BIOLOGICAL SCIENCES↗

Flat and Level Analysis Tool (FLAT) for real-time automated segmentation and analysis of concrete slab point clouds

In the United States, the flatness and levelness of concrete floors during construction is traditionally specified by a maximum allowable gap under a 3 meter straightedge. However, the straightedge method is inexact and rarely representative of the entire floor since the technician is free to choose any location on the floor to perform the measurement. In cases requiring a higher degree of precision and repeatability, concrete floor flatness and levelness can be measured using the standard test method ASTM E1155. With the recent introduction of advanced surveying instruments such as robotic theodolites and terrestrial laser scanners (TLS), the means now exist to modernize and expedite the measurement of floor flatness and levelness. This paper details the development and demonstration of a digital tool, named the Flat and Level Analysis Tool (FLAT), to automate and expedite the segmentation and analysis of flatness and levelness from dense point cloud data of concrete floor slabs. Segmentation algorithms were developed using unsupervised machine learning to extract the set of points belonging to the concrete floor slab from a full 360 scan of a construction site. After segmentation, automated analysis algorithms report the results according to the standard method. The developed algorithms were demonstrated on a dense point cloud captured from a concrete slab-on-grade at a construction site. Results show that the digital tool can quickly provide estimates for floor flatness and levelness with minimal human involvement with comparable accuracy to manual methods.

Hayes, Nolan↗

An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: Tetrahedral molecules

A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB{sub 4} molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q){sup −2} type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH{sub 4} molecule is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving Comfort in Hot-Humid Climates with a Whole-House Dehumidifier, Windermere, Florida [Fact Sheet]

Maintaining comfort in a home can be challenging in hot-humid climates. At the common summer temperature set point of 75°F, the perceived air temperature can vary by 11°F because higher indoor humidity reduces comfort. Often the air conditioner (AC) thermostat set point is lower than the desirable cooling level to try to increase moisture removal so that the interior air is not humid or "muggy." However, this method is not always effective in maintaining indoor relative humidity (RH) or comfort. In order to quantify the performance of a combined whole-house dehumidifier (WHD) AC system, researchers from the U.S. Department of Energy's Building America team Consortium of Advanced Residential Buildings (CARB) monitored the operation of two Lennox AC systems coupled with a Honeywell DH150 TrueDRY whole-house dehumidifier for a six-month period. By using a WHD to control moisture levels (latent cooling) and optimizing a central AC to control temperature (sensible cooling), improvements in comfort can be achieved while reducing utility costs. Indoor comfort for this study was defined as maintaining indoor conditions at below 60% RH and a humidity ratio of 0.012 lbm/lbm while at common dry bulb set point temperatures of 74°-80°F. In addition to enhanced comfort, controlling moisture to these levels can reduce the risk of other potential issues such as mold growth, pests, and building component degradation. Because a standard AC must also reduce dry bulb air temperature in order to remove moisture, a WHD is typically needed to support these latent loads when sensible heat removal is not desired.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Improved Data Interpretation through Identification of Time Series Periodicity Changes

Analysis and interpretation of time series data is easiest when the data values occur at uniform intervals in time, but actual data may have differing data sampling frequencies, such as monthly and daily readings. Applying data analysis techniques, such as smoothing, to such a data set may not give a representative result between time segments. The ability to automatically distinguish time segments of differing data frequency would provide a means for applying data analysis independently to each segment, though a suitable blending at segment boundaries would be required. A method for detecting frequency changes was developed and applied to Gaussian and median smoothing of hydraulic head data from groundwater wells at the U.S. Department of Energy Hanford Site in southeastern Washington state. The process identifies time segments of high-frequency (daily) or low-frequency (greater than daily) data using adjusted-bandwidth Gaussian kernel density estimation and a threshold value, which are further refined to address small blocks of low-frequency data within larger blocks of high-frequency data. User-selectable levels of smoothing are then applied independently to the time segments prior to combining the segment results for a single smoothed data set. This time segment identification approach provides effective low- and high-frequency data separation, which provides a method to apply data analysis independently to each time segment.

97 MATHEMATICS AND COMPUTING↗

Radiation Detection Computational Benchmark Scenarios

Modeling forms an important component of radiation detection development, allowing for testing of new detector designs, evaluation of existing equipment against a wide variety of potential threat sources, and assessing operation performance of radiation detection systems. This can, however, result in large and complex scenarios which are time consuming to model. A variety of approaches to radiation transport modeling exist with complementary strengths and weaknesses for different problems. This variety of approaches, and the development of promising new tools (such as ORNL’s ADVANTG) which combine benefits of multiple approaches, illustrates the need for a means of evaluating or comparing different techniques for radiation detection problems. This report presents a set of 9 benchmark problems for comparing different types of radiation transport calculations, identifying appropriate tools for classes of problems, and testing and guiding the development of new methods. The benchmarks were drawn primarily from existing or previous calculations with a preference for scenarios which include experimental data, or otherwise have results with a high level of confidence, are non-sensitive, and represent problem sets of interest to NA-22. From a technical perspective, the benchmarks were chosen to span a range of difficulty and to include gamma transport, neutron transport, or both and represent different important physical processes and a range of sensitivity to angular or energy fidelity. Following benchmark identification, existing information about geometry, measurements, and previous calculations were assembled. Monte Carlo results (MCNP decks) were reviewed or created and re-run in order to attain accurate computational times and to verify agreement with experimental data, when present. Benchmark information was then conveyed to ORNL in order to guide testing and development of hybrid calculations. The results of those ADVANTG calculations were then sent to PNNL for compilation. This is a report describing the details of the selected Benchmarks and results from various transport codes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On predicting contamination levels of HALOE optics aboard UARS using direct simulation Monte Carlo

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flowfield and surface conditions and geometric orientations in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. Problems resolving species outgassing and vent flux rates that varied over many orders of magnitude were handled using species weighting factors. Results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface are presented, along with data related to code performance. Using procedures developed in standard contamination analyses, the cumulative level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated to be about 2700A.

Woronowicz, Michael S.↗

A Data-Driven Framework for Power System Event Type Identification via Safe Semi-Supervised Techniques

Herein this paper investigates the use of phasor measurement unit (PMU) data with deep learning techniques to construct real-time event identification models for transmission networks. Increasing penetration of distributed energy resources represents a great opportunity to achieve decarbonization, as well as challenges in systematic situational awareness. When high-resolution PMU data and sufficient manually recorded event labels are available, the power event identification problem is defined as a statistical classification problem that can be solved by numerous cutting-edge classifiers. However, in real grids, collecting tremendous high-quality event labels is quite expensive. Utilities frequently have a large number of event records without in-depth details (i.e., unlabeled events). To bridge this gap, we propose a novel semi-supervised learning-based method to improve the performance of event classifiers trained with a limited number of labeled events by exploiting the information from massive unlabeled events. In other words, compared to existing data-driven methods, our method requires only a small portion of labeled data to achieve a similar level of accuracy. Meanwhile, this work discusses and addresses the performance degradation caused by class distribution mismatch between the training set and the real applications. Based on the proposed safe learning mechanism, our model does not directly use all unlabeled events during model training, but selectively uses them through a comprehensive evaluation procedure. Numerical studies on a sizable PMU dataset have been used to validate the performance of the proposed method.

42 ENGINEERING↗

A theoretical investigation of the (0001) covellite surfaces

We report on the properties of the (0001) covellites surfaces, which we investigate by periodic slab density functional theory calculations. The absolute surface energies have been computed for all bulk terminations, showing that surfaces terminated by the flat CuS layer are associated with the lowest surface energy. Cleavage is predicted to occur across the [0001] interlayer Cu–S bond. The surfaces obtained by lowest energy cleavage are analyzed in terms of the atomic vertical relaxation, workfunction, and surface band structure. Our study predicts the presence of a shallow p{sub z}-derived surface state located 0.26 eV below the Fermi level, which is set to play an important role in the surface reactivity of covellite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Failure of Correlation to Describe Out-of-Plane Carbon=Carbon Bending

Carbon-carbon multiply bonded systems are improperly described with standard correlation methods and basis sets. For computations of vibrational modes, the out-of- plane bends can be reported as imaginary at worst or simply too low at best. Utilizing the simplest of aromatic structures (cyclopropenylidene) and various levels of theory, this work diagnoses this known behavior for the first time. A combined 1-particle and n-particle basis set effect conspire to produce these non-physical results. When moving from sp2 to sp3 hybridization in the carbon atoms, the larger number of basis functions overcorrects the energy. This is exacerbated by correlation methods. These allow for occupation of the π and π∗ orbitals in the expanded wave function that combine with the hydrogen s orbitals. As a result, the improperly described space can be further and non-physically stabilized by post-Hartree-Fock correlation. This represents a fundamental problem with at least Hartree-Fock based methods of all flavors in describing carbon. Beyond being a flaw in quantum chemical theory, other repercussions will be present in computations regarding spectroscopy as well as energy and environmental studies where highly-accurate hydrocabon vibrational transitions or thermochemical data are needed.

Fortenberry, Ryan C.↗