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 523 records · Page 29

Unsupervised Power System Event Detection and Classification Using Unlabeled PMU Data

This paper proposes a novel data-driven power system event detection and classification method based on 5TB of actual PMU measurements collected from the US western interconnect. Firstly, a set of comprehensive power quality rules are proposed to pre-filter the raw data and extract the regions of interest (ROI). Six distinct event categories are defined and corresponding patterns are chosen as references. Meanwhile, detailed characteristics of patterns are summarized to enhance our understanding of the actual events. Then, the time-independent feature vectors are generated by extracting the statistical, temporal, and spectral features from the raw time-series data. Furthermore, an ensemble model is proposed to cluster the events by combining multiple K-means clustering models using a voting strategy. Besides, both system-level and PMU-level clustering models are developed. The accuracy and robustness of the event detection method are further improved through interactive evaluation of the two-level clustering results. This paper summarizes the actual characteristics of each event category and provides a reliable basis for accurate label generation. The experiments demonstrate the effectiveness of the proposed event detection and classification method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Unsupervised Power System Event Detection and Classification Using Unlabeled PMU Data

This paper proposes a novel data-driven power system event detection and classification method based on 5TB of actual PMU measurements collected from the US western interconnect. Firstly, a set of comprehensive power quality rules are proposed to pre-filter the raw data and extract the regions of interest (ROI). Six distinct event categories are defined and corresponding patterns are chosen as references. Meanwhile, detailed characteristics of patterns are summarized to enhance our understanding of the actual events. Then, the time-independent feature vectors are generated by extracting the statistical, temporal, and spectral features from the raw time-series data. Furthermore, an ensemble model is proposed to cluster the events by combining multiple K-means clustering models using a voting strategy. Besides, both system-level and PMU-level clustering models are developed. The accuracy and robustness of the event detection method are further improved through interactive evaluation of the two-level clustering results. This paper summarizes the actual characteristics of each event category and provides a reliable basis for accurate label generation. The experiments demonstrate the effectiveness of the proposed event detection and classification method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Exploratory Visual Analysis of Anomalous Runtime Behavior in Streaming High Performance Computing Applications

Online analysis of runtime behavior is essential for performance tuning in streaming scientific workflows. Integration of anomaly detection and visualization is necessary to support human-centered analysis, such as verification of candidate anomalies utilizing domain knowledge. In this work, we propose an efficient and scalable visual analytics system for online performance analysis of scientific workflows toward the exascale scenario. Here, our approach uses a call stack tree representation to encode the structural and temporal information of the function executions. Based on the call stack tree features (e.g., execution time of the root function or vector representation of the tree structure), we employ online anomaly detection approaches to identify candidate anomalous function executions. We also present a set of visualization tools for verification and exploration in a level-of-detailed manner. General information, such as distribution of execution times, are provided in an overview visualization. The detailed structure (e.g., function invocation relations) and the temporal information (e.g., message communication) of the execution call stack of interest are also visualized. The usability and efficiency of our methods are verified in a real-world HPC application.

97 MATHEMATICS AND COMPUTING↗

Characterization and prediction of the electromechanical wear of contact tips during wire arc additive manufacturing of 316L stainless steel

Here, this study seeks to better understand the degradation of the contact tip with respect to WAAM for a 316L wire electrode as well as explore methods of monitoring the contact tip state from process data. The contact tip, a consumable component, positions the wire and serves as the electrical contact surface between the wire electrode and the welding power supply. The wear of the contact tip was characterized in terms of material loss and material contamination for a set of tips worn to discrete levels as measured by the amount of wire fed or arc time. Geometrical characterization found a 49% increase in the bore exit area at 180 meters of wire fed. Machine learning models were developed to predict the relative bore exit area of the contact tip from arc-based process data and a random forest classifier exhibited favorable performance with a cross-validated f1-score of 0.84. The regression architecture implemented a multi-layer perceptron with the ability to predict the relative exit area with an $R^2$ score of 0.75. Key features used in the prediction include the standard deviation of the voltage and the time between shorts.

Contact tip wear↗

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

Behavioral and Population Data-Driven Distribution System Load Modeling

Distribution system residential load modeling and analysis for different geographic areas within a utility or an independent system operator territory are critical for enabling small-scale, aggregated distributed energy resources to participate in grid services under Federal Energy Regulatory Commission Order No. 2222 [1]. In this study, we develop a methodology of modeling residential load profiles in different geographic areas with a focus on human behavior impact. First, we construct a behavior-based load profile model leveraging state-of-the-art appliance models. We simulate human activity and occupancy using Markov chain Monte Carlo methods calibrated with the American Time Use Survey data set. Second, we link our model with cleaned Current Population Survey data from the U.S. Census Bureau. Finally, we populate two sets of 500 households using California and Texas census data, respectively, to perform an initial analysis of the load in different geographic areas with various group features (e.g., different income levels). To distinguish the effect of population behavior differences on aggregated load, we simulate load profiles for both sets assuming fixed physical household parameters and weather data. Analysis shows that average daily load profiles vary significantly by income and income dependency varies by locality.

American Time Use Survey↗

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗

Transfer Learning Trained LSTM Models for Household Load Profile Forecasting

Grid edge renewable energy resources, such as rooftop solar photovoltaics, closely interact with consumer load profiles. Therefore, forecasting future electricity demand, ideally at the individual household level, is indispensable. In this paper, we present a transfer learning enhanced household load profile forecasting method. First, we tune a long short-term memory forecasting model to perform day-ahead prediction of household electricity load profiles. Then we improve these individualized models using transfer learning, and we use k-means clustering to create optimal source data sets. We find average improvements of 4.38% (largest improvement of 10.71%) when the entire data set was used to train the source model and 2.45% (largest improvement of 11.57%) in the mean absolute error when households were first clustered and used to train separate source models for each cluster. We find that transfer learning with clustered data can effectively boost the forecasting performance of the LSTM models. We use realistic household power measurements for 148 real residential households in Austin, Texas.

deep learning↗

Comparison Study of Conventional and Adaptive Mesh Refinement in Organic Material Decomposition Models

This study compares conventional mesh refinement techniques, specifically Uniform Mesh Refinement (UMR), with a new Adaptive Mesh Refinement (AMR) method, applied to Organic Material Decomposition (OMD) models. The proposed benefit of AMR is that only areas that require refinement, based on minimizing a specific field gradient, are refined thus decreasing model wall time compared to conventional UMR methods. This work specifically focuses on comparing UMR and AMR methods on decomposing (both No-Flow and Porous-Flow material models) Polymeric Methylene Diisocyanate (PMDI) polyurethane foam. Throughout the work, the geometry increased in complexity to assess the refinement methods performance at varying levels geometric intricacy. While AMR has been shown to work well in a variety of applications, the UMR approach proved to be computationally faster, for many of the geometries and foam decomposition models, than AMR. However, it was observed that at higher levels of refinement, greater than 3 UMR, AMR begins to be computationally better. Additionally, the settings used to perform AMR greatly impact its performance, and lessons learned, in terms of OMD models, are shared. Due to physics involved in material decomposition, specifically the evolution of state variables, these problems don’t fully benefit from the advantages of AMR.

36 MATERIALS SCIENCE↗

Method for Measuring Dry Density of Compacted Bentonite Mixture Soil Using 3D Laser Scanner data - 19119

The bentonite mixture soil is planned to be used for the low-permeability layer of low-level radioactive waste (LLW) disposal facilities. The construction of the bentonite mixture soil is similar to that of normal soil construction, and the rolling pressure and/or vibration compaction are used. The quality control of construction is the same as normal soil work, and the construction management has to confirm that the compaction degree (that is ratio of the dry density to the modified proctor curve) is above specified value. As a method of measuring the dry density after construction, the conventional methods are the core cutter method and the sand replacement method. The measurement principle of these methods is simple and the reliability of the measurement value is high, but during measurements digging holes disturb the construction place, and the measurement duration is a disadvantage in these methods. Therefore, we developed a non-destructive method that can make measurements in a short time. The developed method uses the correlation between the density change and the layer thickness change before and after compaction, and estimates the density after compaction from the height change obtained from the survey results. However, since the material is compressed while the cross-sectional area spreads, this cross-sectional area change is also taken into consideration. By using the point cloud data of the 3D laser scanner surveying, it is possible to obtain almost continuous data in a short period of time. Corrected dry density values calculated by the above procedure ρ'{sub 3D} are generally consistent with dry density after compaction ρ{sub core} by the core cutter method. This method differs depending on the boundary conditions of the bottom and the side of the construction place. Then, the reliability of the method was evaluated by applying it to two boundary conditions and comparing it with the core cutter method. The proposed estimation method showed good agreement with the core cutter method. Therefore, it was confirmed that it is possible to estimate dry density with high reliability. Although this method needs to be used in set with core cutter method, the measurement frequency of the core cutter method can be reduced, and the measurement of a vast area can be done in a short time without disturbance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Spatial DIC Errors due to Pattern-Induced Bias and Grey Level Discretization

Digital image correlation (DIC) is an optical metrology method widely used in experimental mechanics for full-field shape, displacement and strain measurements. The required strain resolution for engineering applications of interest mandates DIC to have a high image displacement matching accuracy, on the order of 1/100th of a pixel, which necessitates an understanding of DIC errors. In this paper, we examine two spatial bias terms that have been almost completely overlooked. They cause a persistent offset in the matching of image intensities and thus corrupt DIC results. We name them pattern-induced bias (PIB), and intensity discretization bias (IDB). We show that the PIB error occurs in the presence of an undermatched shape function and is primarily dictated by the underlying intensity pattern for a fixed displacement field and DIC settings. The IDB error is due to the quantization of the gray level intensity values in the digital camera. In this paper we demonstrate these errors and quantify their magnitudes both experimentally and with synthetic images.

47 OTHER INSTRUMENTATION↗

Luminaires for Advanced Lighting in Education

Evolving education methods and greater use of technology in the classroom are dictating the need to rethink facility designs, including classroom lighting. Advances in LED-based lighting technology have created the possibility of lighting systems that are not only cost effective and energy efficient, but also color-tunable and as durable as other facility infrastructures (a 20-30 year life expectancy). Thus, there is the opportunity that the modern lighting system can be used by educators as a tool in their teaching strategy. To meet this need, RTI International and Finelite, Inc. teamed to develop and test the Next Generation Integrated Classroom Lighting System (NICLS). The NICLS technology incorporates a high performance, color-tunable light engine into new luminaire designs (e.g., pendant, direct-indirect, downlight, troffers) that are acceptable for use in classrooms. During this project, we successfully demonstrated that the NICLS technology achieves exceptional performance and exceeds all DOE goals for the classroom of the future, including: Luminous efficacy value for NICLS luminaires in excess of 125 lpw at all CCT values; TWL range of 2,700 K to 6,500 K while maintaining a CRI of 82 or higher at all values; Capability for full-range dimming (100% to 1%) at all CCT values with flicker levels below industry guidelines; Performance of the lighting system in a classroom mock-up, incorporating daylight and occupancy sensing to provide automatic control of lighting zones to further reduce energy consumption; Rated lifetime on the system exceeding 50,000 hours with a lumen maintenance of at least 85% at 50,000 hours; and Teacher-focused UI located at the front of the classroom to operate the lighting system. A smartphone-based UI is also available to accommodate teacher movement in the classroom. A critical element of developing this technology is designing the user interface to be compatible with modern teaching methods, including increased use of icons and colors and intuitive appearance. The design of the NICLS technology and the user interface was modified with assistance from focus groups consisting of more than 80 teachers and educational professionals total. The focus groups were held in a full-sized classroom that served as a technology demonstration site for the NICLS. The NICLS technology is an advanced lighting system for educational settings that meets or exceeds all DOE photometric, electrical, and reliability goals for the COF. The NICLS technology has been demonstrated at the classroom level, and the feedback from the dozens of teachers and educational professionals who visited the demonstration site has been overwhelmingly positive. NICLS provides a state-of-the-art lighting environment that adjusts the lighting conditions—both color and illuminance levels—to the needs of students and teachers for the task at hand. Early research has suggested that such lighting conditions will improve not only teacher effectiveness but also a student’s ability to concentrate on learning activities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A search for microscopic black holes, string balls, and sphalerons in proton-proton collisions at $\sqrt{s}=13$ TeV

A search for microscopic black holes, string balls, and electroweak sphalerons using proton-proton collisions at $\sqrt{s}=13$ TeV recorded with the CMS detector at the CERN LHC during the 2016–2018 data taking, and corresponding to an integrated luminosity of 138 fb −1 , is presented. Two search strategies based on control samples in data are used. Model-independent limits on the cross section of physics phenomena with multiple energetic jets, leptons, and photons are set using a method that relies on the shape invariance of the scalar sum of the transverse momenta of all objects in the event. Model-dependent limits on black hole and sphaleron production are set using a newly introduced method that has been developed for the identification of collider events with distinct kinematic features by separating them into classes based on phase space proximity. In the context of models with large extra dimensions, semiclassical black holes and string balls with masses below 8.4–11.4 TeV and 9.0–10.7 TeV, respectively, are excluded at 95% confidence level, significantly extending the reach beyond previous searches. Results of a dedicated search for electroweak sphalerons are used to derive an upper limit of 0.0034 at 95% confidence level on the fraction of quark-quark interactions, with a center-of-mass energy above the nominal sphaleron transition energy threshold of 9 TeV, that result in a sphaleron transition.

Beyond Standard Model↗

Fungal drops: a novel approach for macro- and microscopic analyses of fungal mycelial growth

Abstract This study presents an inexpensive approach for the macro- and microscopic observation of fungal mycelial growth. The ‘fungal drops’ method allows to investigate the development of a mycelial network in filamentous microorganisms at the colony and hyphal scales. A heterogeneous environment is created by depositing 15–20 µl drops on a hydrophobic surface at a fixed distance. This system is akin to a two-dimensional (2D) soil-like structure in which aqueous-pockets are intermixed with air-filled pores. The fungus (spores or mycelia) is inoculated into one of the drops, from which hyphal growth and exploration take place. Hyphal structures are assessed at different scales using stereoscopic and microscopic imaging. The former allows to evaluate the local response of regions within the colony (modular behaviour), while the latter can be used for fractal dimension analyses to describe the hyphal network architecture. The method was tested with several species to underpin the transferability to multiple species. In addition, two sets of experiments were carried out to demonstrate its use in fungal biology. First, mycelial reorganization of Fusarium oxysporum was assessed as a response to patches containing different nutrient concentrations. Second, the effect of interactions with the soil bacterium Pseudomonas putida on habitat colonization by the same fungus was assessed. This method appeared as fast and accessible, allowed for a high level of replication, and complements more complex experimental platforms. Coupled with image analysis, the fungal drops method provides new insights into the study of fungal modularity both macroscopically and at a single-hypha level.

Buffi, Matteo↗

Improved Guarantees for Optimal Nash Equilibrium Seeking and Bilevel Variational Inequalities

We consider a class of hierarchical variational inequality (VI) problems that subsumes VI-constrained optimization and several other problem classes, including the optimal solution selection problem and the optimal Nash equilibrium (NE) seeking problem. Our main contribution is threefold. (i) We consider bilevel VIs with monotone and Lipschitz continuous mappings and devise a single-timescale iteratively regularized extragradient method, named IR-EG 𝚖,𝚖 . We improve the existing iteration complexity results for addressing both bilevel VI and VI-constrained convex optimization problems. (ii) Under the strong monotonicity of the outer-level mapping, we develop a method named IR-EG 𝚜,𝚖 and derive faster guarantees than those in (i). We also study the iteration complexity of this method under a constant regularization parameter. These results appear to be new for both bilevel VIs and VI-constrained optimization. (iii) To our knowledge, complexity guarantees for computing the optimal NE in nonconvex settings do not exist. Motivated by this lacuna, we consider VI-constrained nonconvex optimization problems and devise an inexactly projected gradient method, named IPR-EG, where the projection onto the unknown set of equilibria is performed using IR-EG 𝚜,𝚖 with a prescribed termination criterion and an adaptive regularization parameter. We obtain new complexity guarantees in terms of a residual map and an infeasibility metric for computing a stationary point. Here, we validate the theoretical findings using preliminary numerical experiments for computing the best and the worst NEs.

bilevel optimization↗

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗

Fermilab Booster loss modelling and rebalancing using Bayesian methods

To meet PIP-II upgrade requirements, Fermilab Booster losses need to be reduced by 50% compared to present levels. So far, simulations are not good enough to predict loss patterns. Thus, an extensive Booster tune up will be necessary to achieve required performance. In this paper we present an effort to build a data-driven loss model using Bayesian techniques, and subsequently to rebalance losses for higher trip margins. We first created several sets of spatially and temporally isolated orbit and optics knobs, and trained Gaussian process models for each beam loss monitor as well as beam current. Novel techniques of uncertainty constraints and approximate GP fitting were introduced to handle safety and timing requirements. We then performed single and multi-objective tuning using scalarized objectives comprised of critical beam loss locations. We achieved significant rebalancing of losses, increasing margins by 25%, as well as an overall improvement in transmission efficiency of 0.4%. Automated data collection is being developed so that more accurate surrogate models can be trained over time.

Kuklev, N. [Fermilab]↗

Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies

Joint analyses of small-scale cosmological structure probes have been relatively unexplored yet promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model \texttt{Galacticus} to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗