Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 649 records · Page 36

Automatic Substitute Computed Tomography Generation and Contouring for Magnetic Resonance Imaging (MRI)-Alone External Beam Radiation Therapy From Standard MRI Sequences

Purpose: To validate automatic substitute computed tomography CT (sCT) scans generated from standard T2-weighted (T2w) magnetic resonance (MR) pelvic scans for MR-Sim prostate treatment planning. Patients and Methods: A Siemens Skyra 3T MR imaging (MRI) scanner with laser bridge, flat couch, and pelvic coil mounts was used to scan 39 patients scheduled for external beam radiation therapy for localized prostate cancer. For sCT generation a whole-pelvis MRI scan (1.6 mm 3-dimensional isotropic T2w SPACE [Sampling Perfection with Application optimized Contrasts using different flip angle Evolution] sequence) was acquired. Three additional small field of view scans were acquired: T2w, T2*w, and T1w flip angle 80° for gold fiducials. Patients received a routine planning CT scan. Manual contouring of the prostate, rectum, bladder, and bones was performed independently on the CT and MR scans. Three experienced observers contoured each organ on MRI, allowing interobserver quantification. To generate a training database, each patient CT scan was coregistered to their whole-pelvis T2w using symmetric rigid registration and structure-guided deformable registration. A new multi-atlas local weighted voting method was used to generate automatic contours and sCT results. Results: The mean error in Hounsfield units between the sCT and corresponding patient CT (within the body contour) was 0.6 ± 14.7 (mean ± 1 SD), with a mean absolute error of 40.5 ± 8.2 Hounsfield units. Automatic contouring results were very close to the expert interobserver level (Dice similarity coefficient): prostate 0.80 ± 0.08, bladder 0.86 ± 0.12, rectum 0.84 ± 0.06, bones 0.91 ± 0.03, and body 1.00 ± 0.003. The change in monitor units between the sCT-based plans relative to the gold standard CT plan for the same dose prescription was found to be 0.3% ± 0.8%. The 3-dimensional γ pass rate was 1.00 ± 0.00 (2 mm/2%). Conclusions: The MR-Sim setup and automatic sCT generation methods using standard MR sequences generates realistic contours and electron densities for prostate cancer radiation therapy dose planning and digitally reconstructed radiograph generation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Rapid and flexible segmentation of electron microscopy data using few-shot machine learning

Abstract Automatic segmentation of key microstructural features in atomic-scale electron microscope images is critical to improved understanding of structure–property relationships in many important materials and chemical systems. However, the present paradigm involves time-intensive manual analysis that is inherently biased, error-prone, and unable to accommodate the large volumes of data produced by modern instrumentation. While more automated approaches have been proposed, many are not robust to a high variety of data, and do not generalize well to diverse microstructural features and material systems. Here, we present a flexible, semi-supervised few-shot machine learning approach for segmentation of scanning transmission electron microscopy images of three oxide material systems: (1) epitaxial heterostructures of SrTiO 3 /Ge, (2) La 0.8 Sr 0.2 FeO 3 thin films, and (3) MoO 3 nanoparticles. We demonstrate that the few-shot learning method is more robust against noise, more reconfigurable, and requires less data than conventional image analysis methods. This approach can enable rapid image classification and microstructural feature mapping needed for emerging high-throughput characterization and autonomous microscope platforms.

36 MATERIALS SCIENCE↗

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗

Inverted initial conditions: Exploring the growth of cosmic structure and voids

We introduce and explore “paired” cosmological simulations. A pair consists of an A and B simulation with initial conditions related by the inversion δ A (x,t initial ) = –δ B (x,t initial ) (underdensities substituted for overdensities and vice versa). We argue that the technique is valuable for improving our understanding of cosmic structure formation. The A and B fields are by definition equally likely draws from ΛCDM initial conditions, and in the linear regime evolve identically up to the overall sign. As nonlinear evolution takes hold, a region that collapses to form a halo in simulation A will tend to expand to create a void in simulation B. Applications include (i) contrasting the growth of A-halos and B-voids to test excursion-set theories of structure formation, (ii) cross-correlating the density field of the A and B universes as a novel test for perturbation theory, and (iii) canceling error terms by averaging power spectra between the two boxes. Furthermore, generalizations of the method to more elaborate field transformations are suggested.

79 ASTRONOMY AND ASTROPHYSICS↗

Observational Data for Next-Generation Climate Model Evaluation: Requirements, Considerations, and Best Practices

Climate model simulations are an important source of information about our planet’s climate system and also enable informed decision-making under different future scenarios. As a new archive of results from the next generation of climate models is anticipated to become available with the Coupled Model Intercomparison Project phase 7 (CMIP7), the need to develop efficient and robust methods to evaluate models is paramount. Observations are an integral part of model evaluation, providing a means to quantify and understand the degree to which climate models can faithfully reproduce Earth system processes. Such analysis is critical for constraining climate projections, identifying areas of focus for model development, and assisting analysts in deciphering the utility of models for specific applications. Observations of Earth system come from a diversity of sources, span different space–time domains, and are produced by different communities, and each dataset features different data structures and formats, metadata standards, and its own unique uncertainties. Uncertainties in an observational dataset may stem from gaps in temporal and spatial coverage, instrumentation errors, or assumptions in retrieval and processing methods. How then does one ensure that observational data are ready for use and utilized in the most appropriate way for robust, rapid, and routine climate model evaluation? The CMIP7 Model Benchmarking Task Team with input from the broader climate modeling, model evaluation, and observational data communities present a vision and considerations for best practices toward the optimal and appropriate use of observational data to support next-generation climate model evaluation.

Climate models↗

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↗

Resilience Design Patterns - A Structured Approach to Resilience at Extreme Scale (version 1.1)

Reliability is a serious concern for future extreme-scale high-performance computing (HPC) systems. Projections based on the current generation of HPC systems and technology roadmaps suggest the prevalence of very high fault rates in future systems. The errors resulting from these faults will propagate and generate various kinds of failures, which may result in outcomes ranging from result corruptions to catastrophic application crashes. Therefore the resilience challenge for extreme-scale HPC systems requires management of various hardware and software technologies that are capable of handling a broad set of fault models at accelerated fault rates. Also, due to practical limits on power consumption in HPC systems future systems are likely to embrace innovative architectures, increasing the levels of hardware and software complexities. As a result the techniques that seek to improve resilience must navigate the complex trade-off space between resilience and the overheads to power consumption and performance. While the HPC community has developed various resilience solutions, application-level techniques as well as system-based solutions, the solution space of HPC resilience techniques remains fragmented. There are no formal methods and metrics to investigate and evaluate resilience holistically in HPC systems that consider impact scope, handling coverage, and performance & power efficiency across the system stack. Additionally, few of the current approaches are portable to newer architectures and software environments that will be deployed on future systems. In this document, we develop a structured approach to the management of HPC resilience using the concept of resilience-based design patterns. A design pattern is a general repeatable solution to a commonly occurring problem. We identify the commonly occurring problems and solutions used to deal with faults, errors and failures in HPC systems. Each established solution is described in the form of a pattern that addresses concrete problems in the design of resilient systems. The complete catalog of resilience design patterns provides designers with reusable design elements. We also define a framework that enhances a designer's understanding of the important constraints and opportunities for the design patterns to be implemented and deployed at various layers of the system stack. This design framework may be used to establish mechanisms and interfaces to coordinate flexible fault management across hardware and software components. The framework also supports optimization of the cost-benefit trade-offs among performance, resilience, and power consumption. The overall goal of this work is to enable a systematic methodology for the design and evaluation of resilience technologies in extreme-scale HPC systems that keep scientific applications running to a correct solution in a timely and cost-efficient manner in spite of frequent faults, errors, and failures of various types.

97 MATHEMATICS AND COMPUTING↗

Relationships among climatological vertical moisture structure, column water vapor, and precipitation over the central Amazon in observations and CMIP5 models

Bias and spread in Coupled Model Intercomparison Project Phase 5 simulated vertical specific humidity (q) structure are examined and related to both precipitation and column water vapor (cwv) near Manaus, Brazil, site of the recent Green Ocean Amazon campaign. Simulated seasonal mean q profiles are typically too dry, especially at low levels and during the local dry season, consistent with previously identified surface hydroclimate biases in the Amazon. Multimodel empirical orthogonal function analysis of the models' monthly climatological q profiles indicates two significant modes of ensemble spread in moisture vertical structure, with the leading mode peaked at low levels and the second mode in the lower free troposphere (LFT). While both modes project onto simulated cwv spread, only the first projects on precipitation, suggesting inconsistent sensitivity of simulated rainfall to LFT moisture. Finally, relative to observations, models with high cwv and low-level moisture errors tend to exhibit high precipitation error.

54 ENVIRONMENTAL SCIENCES↗

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↗

21 Tesla Fourier Transform Ion Cyclotron Resonance Mass Spectrometer Greatly Expands Mass Spectrometry Toolbox

We provide the initial performance evaluation of a 21 Tesla Fourier transform ion cyclotron resonance mass spectrometer operating at the Environmental Molecular Sciences Laboratory at Pacific Northwest National Laboratory. The spectrometer constructed for the 21T system employs a commercial dual linear ion trap mass spectrometer coupled to a FTICR spectrometer designed and built in-house. Performance gains from moving to higher magnetic field strength are exemplified by the measurement of peptide isotopic fine structure, complex natural organic matter mixtures, and large proteins. Accurate determination of isotopic fine structure was demonstrated for doubly charged substance P with minimal spectral averaging, and 8,158 molecular formulas assigned to Suwannee River Fulvic Acid standard with RMS error of 10 ppb. We also demonstrated superior performance for intact proteins; namely, broadband isotopic resolution of the entire charge state distribution of apotransferrin (78 kDa) and facile isotopic resolution of monoclonal antibody under a variety of acquisition parameters (e.g. 6 s time-domains with absorption mode processing yielded resolution of approximately 1M at m/z =2,700).

Environmental Molecular Sciences Laboratory↗

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↗