Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model types”

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 1,081 records · Page 60

Screening Cu-Zeolites for Methane Activation Using Curriculum-Based Training

Machine learning (ML), when used synergistically with atomistic simulations, has recently emerged as a powerful tool for accelerated catalyst discovery. However, the application of these techniques has been limited by the lack of interpretable and transferable ML models. In this work, we propose a curriculum-based training (CBT) philosophy to systematically develop reactive machine learning potentials (rMLPs) for high-throughput screening of zeolite catalysts. Our CBT approach combines several different types of calculations to gradually teach the ML model about the relevant regions of the reactive potential energy surface. The resulting rMLPs are accurate, transferable, and interpretable. We further demonstrate the effectiveness of this approach by exhaustively screening thousands of [CuOCu] 2+ sites across hundreds of Cu-zeolites for the industrially relevant methane activation reaction. Specifically, this large-scale analysis of the entire International Zeolite Association (IZA) database identifies a set of previously unexplored zeolites (i.e., MEI, ATN, EWO, and CAS) that show the highest ensemble-averaged rates for [CuOCu] 2+ -catalyzed methane activation. We believe that this CBT philosophy can be generally applied to other zeolite-catalyzed reactions and, subsequently, to other types of heterogeneous catalysts. Thus, this represents an important step toward overcoming the long-standing barriers within the computational heterogeneous catalysis community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗

Introductory guide to backgrounds in XPS spectra and their impact on determining peak intensities

Photoelectron and Auger peaks are central to most of the important uses of x-ray photoelectron spectroscopy (XPS), and thus, they receive the most attention in many types of analysis. Quantitative chemical analysis using XPS requires the assessment of the intensities of the photoemission peaks of the elements detected. Determination of peak intensities requires separation of the photoelectron peak signals from the background on which the peaks rest. For the determination of peak area intensities, the background is subtracted from overall signal intensity. The spectral background is also critical when peak fitting is used to determine intensities of overlapping peaks, and the model of background used in this process can impact the results. In addition to the impact on quantitative analysis, information about the depth distribution of elements in the near surface region can often be obtained by visual inspection of the background and quantified using appropriate modeling. Here, this introductory guide provides some basic information about backgrounds in the XPS analysis, describes the types of background models that are commonly used, suggests some of their strengths and weaknesses, and provides examples of their use and misuse. Although the fundamental nature of some components of the background signals in XPS is not understood, indicating that none of the models in use are fully correct and the area is subject to active research, appropriate good practices have been established for most routine analysis. The guide describes good practices, identifies errors that frequently appear in the literature, and uses examples to demonstrate the impacts of background selections on determinations of peak intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The dark energy survey supernova program: investigating beyond-ΛCDM

We report constraints on a variety of non-standard cosmological models using the full 5-yr photometrically classified type Ia supernova sample from the Dark Energy Survey (DES-SN5YR). Both Akaike Information Criterion (AIC) and Suspiciousness calculations find no strong evidence for or against any of the non-standard models we explore. When combined with external probes, the AIC and Suspiciousness agree that 11 of the 15 models are moderately preferred over Flat-|$\Lambda$|CDM suggesting additional flexibility in our cosmological models may be required beyond the cosmological constant. We also provide a detailed discussion of all cosmological assumptions that appear in the DES supernova cosmology analyses, evaluate their impact, and provide guidance on using the DES Hubble diagram to test non-standard models. An approximate cosmological model, used to perform bias corrections to the data holds the biggest potential for harbouring cosmological assumptions. We show that even if the approximate cosmological model is constructed with a matter density shifted by |$\Delta \Omega _{\rm m}\sim 0.2$| from the true matter density of a simulated data set the bias that arises is subdominant to statistical uncertainties. Nevertheless, we present and validate a methodology to reduce this bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Surrogate model evaluation and building energy benchmarking for commercial buildings

Building energy consumption benchmarking involves challenges associated with various energy patterns for different building types; heating, ventilating, and air-conditioning (HVAC) system types; and climates. Given significant variation in energy use patterns, accurate prediction of long-term energy use using surrogate models remains challenging. Multiple linear regression (MLR) is commonly used for building energy benchmarking because of its simple structure; however, it lacks accuracy compared to other black-box models. Although many studies have compared surrogate models and offer guidance on model selection based on metrics, they do not provide detailed analysis on improving the surrogate model accuracy. In this paper, we implement a surrogate model using polynomial ridge regression (i.e., MLR with interaction terms combined with ridge regularization) for small office and retail strip mall buildings across six HVAC system types and all climate zones, for electricity and natural gas in baseline and proposed scenarios. A simulation workflow is developed using OpenStudio TM /EnergyPlus TM to generate simulation data using measures over a wide range of efficiency inputs. Enhancements based on statistical insights are used for improving the model accuracy using filters, input transformations, and change points. Surrogate models achieved average coefficient of variation of the root mean squared error (CVRMSE) values of 2.17, 1.06, 2.05, and 3.26 for proposed electricity, proposed natural gas, baseline electricity, and baseline natural gas, respectively, with enhancements reducing CVRMSE by an average of 14.9% across all combinations. We provide model interpretation via Shapley additive explanations to determine which input variables most influence energy consumption and provide supportive arguments for enhancements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mapping causal pathways with structural modes fingerprint for perovskite oxides

Abstract Causality is innate to the determination of the fundamental mechanism controlling any physical phenomena. However, combining causality within the standard practices of computational modelling to understand structure-functionality connections is extremely rare. This work proposes a fingerprint based on key structural modes for ABO 3 -type perovskite oxides and its derivatives, combined with causal models, for predicting Kohn–Sham energies. Our study of causal models captures the inherent coupling between structural modes such as rotation, tilt and antiferroelectric displacements, responsible for phase transition, polarization, magnetization and metal–insulator transition, exhibited by these materials. Although developed for modelling specific functionality, this method is universally applicable to derive other functionalities and even different material classes while tracking hidden causal mechanisms via structural distortions.

42 ENGINEERING↗

Study on Fracture of Fiber-Reinforced Polymeric Composites Using Spiral Notch Torsion Test

Fiber-reinforced polymer-matrix composites have been widely applied in various areas due to their excellent mechanical properties. In this study, carbon fiber-reinforced polymer-matrix composites and glass fiber-reinforced polymer-matrix composites failure tests under quasi static torsion have been demonstrated by using spiral notch torsion test methodology. Three types of cylindrical specimens with different fiber orientations were prepared and tested. For both composites, after increase stage, A type exhibited limited stress variations before the sudden failure, whereas B type showed gradual failure with torque drops. Displacement contours, stress contours and energy release rates were acquired using the anisotropic elastic finite element models. Microscopic characterization exhibited that the failure of A type occurred along interfaces between fibers and matrices, whereas circular cracks in B type propagated alternatively between polymeric matrices and interfaces. Type C failure is similar to type A whereas with a curvature observed on the fracture surface.

Du, Fen↗

Mitigation of spatial nonstationarity with vision transformers

Spatial nonstationarity, the location variance of features’ statistical distributions, is ubiquitous in many natural settings. For example, in geological reservoirs rock matrix porosity varies vertically due to geomechanical compaction trends, in mineral deposits grades vary due to sedimentation and concentration processes, in hydrology rainfall varies due to the atmosphere and topography interactions, and in metallurgy crystalline structures vary due to differential cooling. Conventional geostatistical modeling workflows rely on the assumption of stationarity to be able to model spatial features for geostatistical inference. Nevertheless, this is often not a realistic assumption when dealing with nonstationary spatial data and this has motivated a variety of nonstationary spatial modeling workflows such as trend and residual decomposition, cosimulation with secondary features, and spatial segmentation and independent modeling over stationary subdomains. The advent of deep learning technologies has enabled new workflows for modeling spatial relationships. However, there is a paucity of demonstrated best practice and general guidance on mitigation of spatial nonstationarity with deep learning in the geospatial context. We demonstrate the impact of two common types of geostatistical spatial nonstationarity on deep learning model prediction performance and propose the mitigation of such impacts using self-attention (vision transformer) models. We demonstrate the utility of vision transformers for the mitigation of nonstationarity with relative errors as low as 10%, exceeding the performance of alternative deep learning methods such as convolutional neural networks. We establish best practice by demonstrating the ability of self-attention networks for modeling large-scale spatial relationships in the presence of commonly observed geospatial nonstationarity.

58 GEOSCIENCES↗

De novo atomic protein structure modeling for cryoEM density maps using 3D transformer and HMM

Accurately building 3D atomic structures from cryo-EM density maps is a crucial step in cryo-EM-based protein structure determination. Converting density maps into 3D atomic structures for proteins lacking accurate homologous or predicted structures as templates remains a significant challenge. Here, we introduce Cryo2Struct, a fully automated de novo cryo-EM structure modeling method. Cryo2Struct utilizes a 3D transformer to identify atoms and amino acid types in cryo-EM density maps, followed by an innovative Hidden Markov Model (HMM) to connect predicted atoms and build protein backbone structures. Cryo2Struct produces substantially more accurate and complete protein structural models than the widely used ab initio method Phenix. Additionally, its performance in building atomic structural models is robust against changes in the resolution of density maps and the size of protein structures.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗

Moving beyond post hoc explainable artificial intelligence: a perspective paper on lessons learned from dynamical climate modeling

AI models are criticized as being black boxes, potentially subjecting climate science to greater uncertainty. Explainable artificial intelligence (XAI) has been proposed to probe AI models and increase trust. In this review and perspective paper, we suggest that, in addition to using XAI methods, AI researchers in climate science can learn from past successes in the development of physics-based dynamical climate models. Dynamical models are complex but have gained trust because their successes and failures can sometimes be attributed to specific components or sub-models, such as when model bias is explained by pointing to a particular parameterization. We propose three types of understanding as a basis to evaluate trust in dynamical and AI models alike: (1) instrumental understanding, which is obtained when a model has passed a functional test; (2) statistical understanding, obtained when researchers can make sense of the modeling results using statistical techniques to identify input–output relationships; and (3) component-level understanding, which refers to modelers' ability to point to specific model components or parts in the model architecture as the culprit for erratic model behaviors or as the crucial reason why the model functions well. We demonstrate how component-level understanding has been sought and achieved via climate model intercomparison projects over the past several decades. Such component-level understanding routinely leads to model improvements and may also serve as a template for thinking about AI-driven climate science. Currently, XAI methods can help explain the behaviors of AI models by focusing on the mapping between input and output, thereby increasing the statistical understanding of AI models. Yet, to further increase our understanding of AI models, we will have to build AI models that have interpretable components amenable to component-level understanding. We give recent examples from the AI climate science literature to highlight some recent, albeit limited, successes in achieving component-level understanding and thereby explaining model behavior. The merit of such interpretable AI models is that they serve as a stronger basis for trust in climate modeling and, by extension, downstream uses of climate model data.

54 ENVIRONMENTAL SCIENCES↗

Application of multi-method-multi-model inference to radiation related solid cancer excess risks models for astronaut risk assessment

The impact of including model-averaged excess radiation risks (ER) into a measure of radiation attributed decrease of survival (RADS) for the outcome all solid cancer incidence and the impact on the uncertainties is demonstrated. It is shown that RADS applying weighted model averaged ER based on AIC weights result in smaller risk estimates with narrower 95% CI than RADS using ER based on BIC weights. Further a multi-method-multi-model inference approach is introduced that allows calculating one general RADS estimate providing a weighted average risk estimate for a lunar and a Mars mission. For males the general RADS estimate is found to be 0.42% (95% CI: 0.38%; 0.45%) and for females 0.67% (95% CI: 0.59%; 0.75%) for a lunar mission and 2.45% (95% CI: 2.23%; 2.67%) for males and 3.91% (95% CI: 3.44%; 4.39%) for females for a Mars mission considering an age at exposure of 40 years and an attained age of 65 years. It is recommended to include these types of uncertainties and to include model-averaged excess risks in astronaut risk assessment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Atmospheric Feedbacks Reverse the Sensitivity of Modeled Photosynthesis to Stomatal Function

Stomata mediate fluxes of carbon and water between terrestrial plants and the atmosphere. These fluxes are governed by stomatal function and can be modulated in many Earth system models by an empirical parameter within the calculation of stomatal conductance, the stomatal slope (𝑔 1⁢𝑀 ). Intuitively, 𝑔 1⁢𝑀 represents the marginal water cost of carbon, relating it to the emergent plant property of water use efficiency. Observations show that 𝑔 1⁢𝑀 can range widely across and within plant types in varying environments, and this distribution of 𝑔 1⁢𝑀 is not captured within Earth system models which represent each plant type with a single 𝑔 1⁢𝑀 value. Here we examine how 𝑔 1⁢𝑀 influences photosynthesis using coupled Earth system model simulations by perturbing 𝑔 1⁢𝑀 to observed 5⁢t⁢h and 95⁢t⁢h percentiles for each plant type. We find that high 𝑔 1⁢𝑀 reduces photosynthesis nearly everywhere, while low 𝑔 1⁢𝑀 has regionally dependent responses. Under fixed atmospheric conditions, low 𝑔 1⁢𝑀 increases photosynthesis in the Amazon and central North America but decreases photosynthesis in boreal Canada. These responses reverse when the atmosphere responds interactively due to spatially differing sensitivity to increases in temperature and vapor pressure deficit. Choice of 𝑔 1⁢𝑀 also influences photosynthetic response to changes in atmospheric carbon dioxide (CO 2 ), with lower and higher 𝑔1⁢𝑀 modifying total global response to elevated 2x preindustrial CO 2 by +6.4% and −9.6%, respectively. Our work demonstrates that atmospheric feedbacks are critical for determining the photosynthetic response to 𝑔 1⁢𝑀 assumptions and some regions are particularly sensitive to choice of 𝑔 1⁢𝑀 .

Liu, Amy X. [University of Washington, Seattle, WA↗

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An extended polynomial chaos expansion for PDF characterization and variation with aleatory and epistemic uncertainties

This paper presents an extended polynomial chaos formalism for epistemic uncertainties and a new framework for evaluating sensitivities and variations of output probability density functions (PDF) to uncertainty in probabilistic models of input variables. An ”extended” polynomial chaos expansion (PCE) approach is developed that accounts for both aleatory and epistemic uncertainties, modeled as random variables, thus allowing a unified treatment of both types of uncertainty. We explore in particular epistemic uncertainty associated with the choice of prior probabilistic models for input parameters. A PCE-based Kernel Density (KDE) construction provides a composite map from the PCE coefficients and germ to the PDF of quantities of interest (QoI). Here, the sensitivities of these PDF with respect to the input parameters are then evaluated. Input parameters of the probabilistic models are considered. By sampling over the epistemic random variable, a family of PDFs is generated and the failure probability is itself estimated as a random variable with its own PCE. Integrating epistemic uncertainties within the PCE framework results in a computationally efficient paradigm for propagation and sensitivity evaluation. Two typical illustrative examples are used to demonstrate the proposed approach.

Aleatory uncertainty↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids: Preprint

Increased deployment of Distributed Energy Resources (DERs) and microgrids are significantly changing the operation of distribution grid. Increasing penetration of power electronics based DERs also raise several operational problems such as protection coordination issues, voltage raise, recloser synchronization, stability, and frequency control. Legacy systems are designed for unidirectional power flow. Bidirectional power flows, and varying short circuit levels are some of the major challenges in implementing these legacy protection schemes for future power grid with high penetration of power electronics based DERs. Inverter based DER’s can contribute to a maximum of 2 p.u. fault current. Interconnection standards might require even lower fault current contribution. If these protection related challenges are not addressed, it will result in limiting the penetration of DERs in distribution feeders. In this paper, the feasibility of high frequency signature generated by transients is evaluated to identify and locate faults in a distribution system. The paper will discuss the need for an accurate line modeling approach to accurately capture the high frequency signatures. These high frequency signatures are naturally occurring electromagnetic transients for events such as switching events in the network, faults events, etc. This is different than that implemented in transmission systems as the proposed approach utilizes advanced signal processing techniques to retrieve valuable information from the high frequency signatures generated during faults and does not just rely on time delay measurements between reflections of the high frequency signals. In the full paper, information of a modified IEEE distribution feeder test system to accurately model the high frequency signatures will be provided. The distribution lines will be modeled using universal line modeling (ULM) technique in an electromagnetic transients program (EMTP) modeling platform. Different test cases will be simulated and these test cases will include different fault types, multiple fault locations, capacitor bank switching, etc. The feeder system modeled for these simulations will also include power electronics based DERs. Additionally, experimental results demonstrating the high frequency signatures generated due to a transient event in a medium voltage (13.8 kV) distribution line will also be presented in the full paper. Furthermore, the analysis combined with simulation and experimental results in the full paper will help demonstrate the feasibility and efficacy of the proposed protection scheme. The developed protection scheme will not just help solve problems generated by bidirectional current flow, varying short-circuit current capacity, but also help in proliferation of DERs in distribution systems.

41 EE - Solar Energy Technologies Office (EE-4S)↗