Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural error”

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 541 records · Page 30

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

IDAES-PSE 2.5.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.5.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New diagnostics check for near-parallel variables and constraints. New diagnostics tools for identifying causes of infeasibility in models. New example for creating a custom model of a liquid-liquid extractor unit operation. Bug Fixes Fixed bug in Gibbs reactor that caused it to appear to have additional spurious degrees of freedom. Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Testing and Robustness Deployed the IDAES Diagnostics Toolbox to confirm that there are no structural or numerical issues in the core model libraries. Additional robustness tests for core model, and some associated improvements in the converge tester class. Fixed a number of issues that were causing unexpected warnings to be emitted during testing. Deprecations and Removals Removed examples for RIPE tool which has not been supported for a number of releases.

AS↗

Accelerating Discovery of Atomistic Defects via Machine Learning

The quantification of defects such as vacancies in crystalline structures is a cornerstone of materials science research. Traditional efforts often rely on manual detection, a process that is time-intensive, prone to human error, and challenging to scale. Here we leverage machine learning (ML) methods to identify and quantify vacancies within a crystalline lattice, aiming to expedite detection while improving accuracy. Additionally, we explore the transferability of these ML techniques, identifying characteristics of atomistic imaging data that complicate this task. We show how the integration of ML can drive innovation, providing a powerful tool that will play an increasingly crucial role in the future of materials science.

2D materials↗

Subgrid Treatment of Spacer Grids in the 2D/1D Subplane Approach

Recent efforts have focused on improving the performance of the 2D/1D subplane implementation in MPACT. One of the developments recently added, presented here, focuses on treating spacer grids using a subgrid approach within each subplane, avoiding the extra computation and potentially improving stability issues that can be present by explicitly modeling the spacer grid with a separate MOC plane. From the cases analyzed, it seems that the approach does a sufficient job of maintaining the accuracy of explicitly meshing the spacer grid (<10 pcm and <1% max pin power error) while offering relatively modest performance gains. Overall, the success of the approach with spacer grids opens up possibilities to apply similar approaches to other structural components, particularly those in the axial reflectors to reduce the number of MOC planes and overall computational burden.

Stimpson, Shane G.↗

Development and Validation of Low-Cost, High-Reflectance Composite CSP Facets: SIPS Final Report

This work investigates the various challenges associated with developing heliostat structural composite facets using 1 mm glass mirrors. Such facets are desirable for Concentrating Solar Power because 1 mm glass mirrors provide an absolute increase in reflectivity of 2-3 % over the industry standard of 4 mm glass mirrors. Prototypes of paraboloid composite facets with 1 mm glass mirrors were constructed that have a root mean square slope error on the order of 2 mrad while achieving greater than 96% reflectivity. These facets were constructed using a low-quality aluminum mold and a bill of materials that show potential to achieve cost parity with existing 4 mm glass mirrors supported by structural steel. The facets produced were able to survive up to 50 mm hail ball impacts and were robust against accelerated environmental cycling designed to expose durability concerns. Further work is needed to develop a scalable and cost-effective manufacturing process with a similar bill of materials.

14 SOLAR ENERGY↗

Orbital angular momentum constraints in the variational optimization of the two-electron reduced-density matrix

The direct variational determination of the two-electron reduced-density matrix (2-RDM) corresponding to an atomic or molecular system is usually carried out in a basis of real-valued atom-centered Gaussian basis functions, under the assumption that the 2-RDM is a real-valued quantity. Furthermore, for systems that possess orbital angular momentum symmetry, the description of states with a well-defined, nonzero z projection of the orbital angular momentum requires a computational framework generalized to include either complex basis functions or a complex-valued 2-RDM. We consider a semidefinite program suitable for the direct optimization of a complex-valued 2-RDM and explore the role of orbital angular momentum constraints in systems that possess the relevant symmetries. For atomic systems, constraints on the expectation values of the square and z projection of the orbital angular momentum operator allow one to optimize 2-RDMs for multiple orbital angular momentum states. Similarly, in linear molecules, orbital angular momentum projection constraints enable the description of multiple electronic states and, moreover, the application of such constraints is essential for a qualitatively correct description of the electronic structure. For example, in the case of molecular oxygen, we demonstrate that orbital angular momentum constraints are necessary to recover the correct energy ordering of the lowest-energy singlet and triplet states near the equilibrium geometry. However, care must still be taken in the description of the dissociation limit, because the 2-RDM-based approach is not size consistent, and the size-consistency error varies dramatically, depending on the z projections of the spin and orbital angular momenta.

74 ATOMIC AND MOLECULAR PHYSICS↗

Predictions of macro-scale fracture geometries from acoustic emission point cloud data in a hydraulic fracturing experiment

Observations of laboratory fracture testing by means of acoustic emission (AE) can provide a wealth of information regarding the fracturing process and the subsequent damage of the material or structure under load. A method for determining the structure of macro-scale fractures from a point cloud of AE events was developed and tested at the laboratory scale. An unconfined hydraulic fracturing experiment was performed on granite while monitoring acoustic emissions from six piezoelectric transducers on the surfaces of the specimen. The granite specimen dimensions were 30 × 30 × 25 cm 3 . The motivation of the AE analysis was to provide location information for a secondary wellbore placement which intersects the hydraulic fracture to complete a hydraulically connected binary well-hydraulic fracture system. Information gained from the 3D event source locations was used to optimize the location and orientation of secondary production well placement to intersect the induced hydraulic fracture. A direct method was developed to predict the location, orientation and shape of the hydraulic fracture from a dense cloud of AE events. The predicted fracture plane, which was assumed to be planar, was rotated in both the pitch and roll directions, while an average error of AE event distance between the assumed plane and the individual event locations was determined for each rotation iteration, resulting with a predicted fracture plane in the minimum error orientation. Post-test fracture observations were made and digitized in 3D space from coring and slabbing the granite specimen. The actual fracture observations were compared with the predicted fracture structure from AE and showed marked correlation. In conclusion, the simple method of determining macro-scale fracture location, orientation, and extents provided a useful tool for materials that exhibit large and disperse clouds of AE activity where coalesced fracture structure is not apparent.

42 ENGINEERING↗

Property Prediction of Organic Donor Molecules for Photovoltaic Applications Using Extremely Randomized Trees

Organic solar cells are an inexpensive, flexible alternative to traditional silicon-based solar cells but disadvantaged by low power conversion efficiency due to empirical design and complex manufacturing processes. This process can be accelerated by generating a comprehensive set of potential candidates. However, this would require a laborious trial and error method of modeling all possible polymer configurations. A machine learning model has the potential to accelerate the process of screening potential donor candidates by associating structural features of the compound using molecular fingerprints with their highest occupied molecular orbital energies. In this paper, extremely randomized tree learning models are employed for the prediction of HOMO values for donor compounds, and a web application is developed. Finally, the proposed models outperform neural networks trained on molecular fingerprints as well as SMILES, as well as other state-of-the-art architectures such as Chemception and Molecular Graph Convolution on two datasets of varying sizes.

42 ENGINEERING↗

A systematic comparison of machine learning methods for modeling of dynamic processes applied to combustion emission rate modeling

Ten established, data-driven dynamic algorithms are surveyed and a practical guide for understanding these methods generated. Existing Python programming packages for implementing each algorithm are acknowledged, and the model equations necessary for prediction are presented. A case study on a coal-fired power plant’s NO x emission rates is performed, directly comparing each modeling method’s performance on a mutual system. Each model is evaluated by its root mean squared error (RMSE) on out-of-sample future horizon predictions. Optimal hyperparameters are identified using either an exhaustive search or genetic algorithm. The top five model structures of each method are used to recursively predict future NO x emission rates over a 60-step time horizon. The RMSE at each future timestep is determined, and the recursive output prediction trends compared against measurements in time. The GRU neural network is identified as the best candidate for representing the system, demonstrating accurate and stable predictions across the future horizon by all considered models, while satisfactory performance was observed in several of the ARX/NARX formulations. Finally, these efforts have contributed 1) a concise resource of multiple proven dynamic machine learning methods, 2) a practical guide explaining the use of these methods, effectively lowering the “barrier-to-entry” of deploying such models in control systems, 3) a comparison study evaluating each method’s performance on a mutual system, 4) demonstration of accurate multi-timestep emissions modeling suitable for systems-level control, and 5) generalizable results demonstrating the suitability of each method for prediction over a multi-step future horizon to other complex dynamic systems.

42 ENGINEERING↗

Residual Setup Errors Towards the Heart After Image Guidance Linked With Poorer Survival in Lung Cancer Patients: Do We Need Stricter IGRT Protocols?

Image guided radiation therapy (IGRT) is widely used, but data directly relating set-up errors to patient outcome are scarce. This study investigates the relationship between residual IGRT shifts and overall patient survival and uses the observed relations to identify structures sensitive to radiation dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

Identifying an efficient, thermally robust inorganic phosphor host via machine learning

Abstract Rare-earth substituted inorganic phosphors are critical for solid state lighting. New phosphors are traditionally identified through chemical intuition or trial and error synthesis, inhibiting the discovery of potential high-performance materials. Here, we merge a support vector machine regression model to predict a phosphor host crystal structure’s Debye temperature, which is a proxy for photoluminescent quantum yield, with high-throughput density functional theory calculations to evaluate the band gap. This platform allows the identification of phosphors that may have otherwise been overlooked. Among the compounds with the highest Debye temperature and largest band gap, NaBaB 9 O 15 shows outstanding potential. Following its synthesis and structural characterization, the structural rigidity is confirmed to stem from a unique corner sharing [B 3 O 7 ] 5– polyanionic backbone. Substituting this material with Eu 2+ yields UV excitation bands and a narrow violet emission at 416 nm with a full-width at half-maximum of 34.5 nm. More importantly, NaBaB 9 O 15 :Eu 2+ possesses a quantum yield of 95% and excellent thermal stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learned features from density of states for accurate adsorption energy prediction

Materials databases generated by high-throughput computational screening, typically using density functional theory (DFT), have become valuable resources for discovering new heterogeneous catalysts, though the computational cost associated with generating them presents a crucial roadblock. Hence there is a significant demand for developing descriptors or features, in lieu of DFT, to accurately predict catalytic properties, such as adsorption energies. Here, we demonstrate an approach to predict energies using a convolutional neural network-based machine learning model to automatically obtain key features from the electronic density of states (DOS). The model, DOSnet, is evaluated for a diverse set of adsorbates and surfaces, yielding a mean absolute error on the order of 0.1 eV. In addition, DOSnet can provide physically meaningful predictions and insights by predicting responses to external perturbations to the electronic structure without additional DFT calculations, paving the way for the accelerated discovery of materials and catalysts by exploration of the electronic space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometrical picture of the electron–electron correlation at the large- D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron–electron correlations and might give a bound on the correlation energy. Furthermore, the results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron–electron correlations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hierarchical modeling of molecular energies using a deep neural network

In this work, we introduce the Hierarchically Interacting Particle Neural Network (HIP-NN) to model molecular properties from datasets of quantum calculations. Inspired by a many-body expansion, HIP-NN decomposes properties, such as energy, as a sum over hierarchical terms. These terms are generated from a neural network—a composition of many nonlinear transformations—acting on a representation of the molecule. HIP-NN achieves the state-of-the-art performance on a dataset of 131k ground state organic molecules and predicts energies with 0.26 kcal/mol mean absolute error. Finally, with minimal tuning, our model is also competitive on a dataset of molecular dynamics trajectories. In addition to enabling accurate energy predictions, the hierarchical structure of HIP-NN helps to identify regions of model uncertainty.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Augmenting machine learning of energy landscapes with local structural information

We present a machine learning approach for accurately predicting formation energies of binary compounds in the context of crystal structure predictions. The success of any machine learning model depends significantly on the choice of representation used to encode the relevant physical information into machine-learnable data. We test different representation schemes based on partial radial and angular distribution functions (RDF+ADF) on Al–Ni and Cd–Te structures generated using our genetic algorithm for structure prediction. We observe a remarkable improvement in predictive accuracy upon transitioning from global to atom-centered representations, resulting in a threefold decrease in prediction errors. We show that a support vector regression model using a combination of atomic radial and angular distribution functions performs best at the formation energy prediction task, providing small root mean squared errors of 3.9 meV/atom and 10.9 meV/atom for Al–Ni and Cd–Te, respectively. We test the performance of our models against common traditional descriptors and find that RDF- and ADF-based representations significantly outperform many of those in the prediction of formation energies. The high accuracy of predictions makes our machine learning models great candidates for the exploration of energy landscapes.

Honrao, Shreyas J. (ORCID:0000000166155257)↗

Influence of fourth-order vibrational corrections on semi-experimental (reSE) structures of linear molecules

Semi-experimental structures (reSE) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level ab initio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10−5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10−4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4.

Chemistry↗

Rotation and neoclassical ripple transport in ITER

Neoclassical transport in the presence of non-axisymmetric magnetic fields causes a toroidal torque known as neoclassical toroidal viscosity (NTV). The toroidal symmetry of ITER will be broken by the finite number of toroidal field coils and by test blanket modules (TBMs). The addition of ferritic inserts (FIs) will decrease the magnitude of the toroidal field ripple. 3D magnetic equilibria in the presence of toroidal field ripple and ferromagnetic structures are calculated for an ITER steady-state scenario using the Variational Moments Equilibrium Code (VMEC). Furthermore, neoclassical transport quantities in the presence of these error fields are calculated using the Stellarator Fokker-Planck Iterative Neoclassical Conservative Solver (SFINCS).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗