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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Denoising of imaginary time response functions with Hankel projections

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society 2024

Yu, Yang (ORCID:0000000186178878)↗
DOE OSTI2025

Extended Galerkin Neural Network Approximation of Singular Variational Problems with Error Control

We present extended Galerkin neural networks, a variational framework for approximating general boundary value problems (BVPs) with error control. The main contributions of this work are (1) a rigorous theory guiding the construction of new weighted least squares variational formulations suitable for use in neural network approximation of general BVPs, and (2) an “extended” feedforward network architecture which incorporates and is even capable of learning singular solution structures, thus greatly improving approximability of singular solutions. Furthermore, numerical results are presented for several problems, including steady Stokes flow around reentrant corners and in convex corners with Moffatt eddies in order to demonstrate efficacy of the method.

a posteriori error estimate↗
DOE OSTI2016

3D Monte Carlo simulation of backscattered electron signal variation across pore-solid boundaries in cement-based materials

Three-dimensional (3D) Monte Carlo simulation was used to study the variation of backscattered electron (BSE) signal across pore-solid boundaries in cement-based materials in order to enhance quantitative analysis of pore structure. The effects of pore size, depth and boundary inclination angle were investigated. It is found that pores down to 1 nm can generate sufficient contrast to be detected. Visibility improves with larger pore size, smaller beam probe size and lower acceleration voltage. However, pixels in shallow pores or near pore boundaries display higher grey values (brightness) than expected due to sampling sub-surface or neighbouring solid material. Thus, cement-based materials may appear less porous or the pores appear smaller than they actually are in BSE images. Simulated BSE images were used to test the accuracy of the Overflow pore segmentation method. Results show the method is generally valid and gives low errors for pores that are 1 μm and greater.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Random forest prediction of crystal structure from electron diffraction patterns incorporating multiple scattering

Diffraction is the most common method to solve for unknown or partially known crystal structures. However, it remains a challenge to determine the crystal structure of a new material that may have nanoscale size or heterogeneities. Here, in this study, we train an architecture of hierarchical random forest models capable of predicting the crystal system, space group, and lattice parameters from one or more unknown two-dimensional electron diffraction patterns. Our initial model correctly identifies the crystal system of a simulated electron diffraction pattern from a 20-nm-thick specimen of arbitrary orientation 67% of the time. We achieve a topline accuracy of 79% when aggregating predictions from ten patterns of the same material but different zone axes. The space group and lattice predictions range from 70% to 90% accuracy and median errors of 0.01-0.5Å, respectively, for cubic, hexagonal, trigonal, and tetragonal crystal systems while being less reliable on orthorhombic and monoclinic systems. We apply this architecture to a four-dimensional scanning transmission electron microscopy scan of gold nanoparticles, where it accurately predicts the crystal structure and lattice constants. These random forest models can be used to significantly accelerate the analysis of electron diffraction patterns, particularly in the case of unknown crystal structures. Additionally, due to the speed of inference, these models could be integrated into live transmission electron microscopy experiments, allowing real-Time labeling of a specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Heterogeneous Composition of Oxygen-Evolving Complexes in Crystal Structures of Dark-Adapted Photosystem II

Photosystem II (PSII) is a homodimeric protein complex that catalyzes water oxidation at the oxygen-evolving complex (OEC), a heterocubanoid calcium-tetra-manganese cluster. Here, we analyze the omit electron-density peaks of the OEC’s metal ions in five X-ray free-electron laser PSII structures at a resolution between 2.15 Å and 1.95 Å. The omit peaks can be described by the total number of electrons and approximated by the variance of electron density distribution when the distributions are spherically symmetric. We show that the number of electrons of metal centers are different in the two OECs of PSII dimers, implying that either the oxidation states and/or occupancies of individual metal ions are different in the two monomers. In either case, we find that the two OECs of dark-adapted PSII dimers in crystals are not fully synchronized in the same expected S 1 state. Differences in redox states of the OEC in PSII only partially account for the observation that the electron densities integrate to a smaller number of electrons than expected. Differences between the determined and expected relative electron numbers are much larger than the estimated errors, indicating heterogeneity in the OEC composition. Here, the reported findings suggest that two OECs of dark-adapted PSII crystals are typically not fully synchronized in the S 1 state and likely generate heterogeneity in the crystallographical data of PSII intermediate structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Bootstrap embedding: An internally consistent fragment-based method

Strong correlation poses a difficult problem for electronic structure theory, with computational cost scaling quickly with system size. Fragment embedding is an attractive approach to this problem. By dividing a large complicated system into smaller manageable fragments “embedded” in an approximate description of the rest of the system, we can hope to ameliorate the steep cost of correlated calculations. While appealing, these methods often converge slowly with fragment size because of small errors at the boundary between fragment and bath. We describe a new electronic embedding method, dubbed “Bootstrap Embedding,” a self-consistent wavefunction-in-wavefunction embedding theory that uses overlapping fragments to improve the description of fragment edges. We apply this method to the one dimensional Hubbard model and a translationally asymmetric variant, and find that it performs very well for energies and populations. We find Bootstrap Embedding converges rapidly with embedded fragment size, overcoming the surface-area-to-volume-ratio error typical of many embedding methods. We anticipate that this method may lead to a low-scaling, high accuracy treatment of electron correlation in large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Electronic structure of mononuclear Cu-based molecule from density-functional theory with self-interaction correction

In this paper, we investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C 6 H 4 S 2 ) 2 ] z in two oxidation states ( z = –2, –1) using density-functional theory (DFT) with Fermi–Löwdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3 d orbitals. The FLO-SIC method relies on optimization of Fermi–Löwdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin S = 1/2), FLO-SIC spin density originates from the Cu d and S p orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and, in particular, the 3 d orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO–LUMO gap of the dianionic state is larger than that of the monoanionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3 d -element-based systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Verification of Accuracy of CyberKnife Tumor-tracking Radiation Therapy Using Patient-specific Lung Phantoms

Purpose: To investigate the accuracy of the CyberKnife Xsight Lung Tracking System (XLTS) compared with that of a fiducial-based target tracking system (FTTS) using patient-specific lung phantoms. Methods and Materials: Three-dimensional printing technology was used to make individualized lung phantoms that closely mimicked the lung anatomy of actual patients. Based on planning computed tomographic data from 6 lung cancer patients who underwent stereotactic ablative radiation therapy using the CyberKnife, the volume above a certain Hounsfield unit (HU) was assigned as the structure to be filled uniformly with polylactic acid material by a 3-dimensional printer (3D Edison, Lokit, Korea). We evaluated the discrepancies between the measured and modeled target positions, representing the total tracking error, using 3 log files that were generated during each treatment for both the FTTS and the XLTS. We also analyzed the γ index between the film dose measured under the FTTS and XLTS. Results: The overall mean values and standard deviations of total tracking errors for the FTTS were 0.36 ± 0.39 mm, 0.15 ± 0.64 mm, and 0.15 ± 0.62 mm for the craniocaudal (CC), left–right (LR), and anteroposterior (AP) components, respectively. Those for the XLTS were 0.38 ± 0.54 mm, 0.13 ± 0.18 mm, and 0.14 ± 0.37 mm for the CC, LR, and AP components, respectively. The average of γ passing rates was 100% for the criteria of 3%, 3 mm; 99.6% for the criteria of 2%, 2 mm; and 86.8% for the criteria of 1%, 1 mm. Conclusions: The XLTS has segmentation accuracy comparable with that of the FTTS and small total tracking errors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Linear-Response Time-Dependent Embedded Mean-Field Theory

We present a time-dependent (TD) linear-response description of excited electronic states within the framework of embedded mean-field theory (EMFT). TD-EMFT allows for subsystems to be described at different mean-field levels of theory, enabling straightforward treatment of excited states and transition properties. We provide benchmark demonstrations of TD-EMFT for both local and nonlocal excitations in organic molecules, as well as applications to chlorophyll a, solvatochromic shifts of a dye in solution, and sulfur K-edge X-ray absorption spectroscopy (XAS). It is found that mixed-basis implementations of TD-EMFT lead to substantial errors in terms of transition properties; however, as previously found for ground-state EMFT, these errors are largely eliminated with the use of Fock-matrix corrections. These results indicate that TD-EMFT is a promising method for the efficient, multilevel description of excited-state electronic structure and dynamics in complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Optical properties of semiconductor nanoparticles doped with 3d and 4f block elements

Doped inorganic nanoparticles (NPs) have enabled researchers across multiple disciplines to create new materials with unique optical, electronic, magnetic, electrical, chemical, and/or compositional properties. In the quest for generating newer functional materials with controlled optical and electronic properties, the use of d- and f-block elements as dopants in the inorganic NPs remain at the forefront of research. This article summarizes experiments that identify the important control parameters, in the form of NP core and surface properties, for promoting host-sensitized dopant emission from d- and f-doped inorganic/semiconductor NPs. We conclude that energy offsets between the dopant and NP electronic levels and an understanding of the charge trapping mechanisms can be used to predict dopant luminescence, and we claim that these features can be used to choose NP host materials by design, rather than by trial and error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2017

A Pattern Language for High-Performance Computing Resilience

High-performance computing systems (HPC) provide powerful capabilities for modeling, simulation, and data analytics for a broad class of computational problems. They enable extreme performance of the order of quadrillion floating-point arithmetic calculations per second by aggregating the power of millions of compute, memory, networking and storage components. With the rapidly growing scale and complexity of HPC systems for achieving even greater performance, ensuring their reliable operation in the face of system degradations and failures is a critical challenge. System fault events often lead the scientific applications to produce incorrect results, or may even cause their untimely termination. The sheer number of components in modern extreme-scale HPC systems and the complex interactions and dependencies among the hardware and software components, the applications, and the physical environment makes the design of practical solutions that support fault resilience a complex undertaking. To manage this complexity, we developed a methodology for designing HPC resilience solutions using design patterns. We codified the well-known techniques for handling faults, errors and failures that have been devised, applied and improved upon over the past three decades in the form of design patterns. In this paper, we present a pattern language to enable a structured approach to the development of HPC resilience solutions. The pattern language reveals the relations among the resilience patterns and provides the means to explore alternative techniques for handling a specific fault model that may have different efficiency and complexity characteristics. Using the pattern language enables the design and implementation of comprehensive resilience solutions as a set of interconnected resilience patterns that can be instantiated across layers of the system stack.

Hukerikar, Saurabh↗
DOE OSTI2015

Simplified Methods for Combining Natural and Mechanical Ventilation

In determining ventilation rates, it is often necessary to combine naturally-driven ventilation, such as infiltration, with mechanical systems. Modern calculation methods are sufficiently powerful that this can be done from first principles with time varying flows, but for some purposes simplified methods of combining the mechanical and natural ventilation are required—we call this “superposition”. An example of superposition would be ventilation standards that may pre-calculate some quantities within the body of the standard. When there are balanced mechanical systems, the solution is simple additivity, because a balanced system does not impact the internal pressure of the space. Unbalanced systems, however, change internal pressures and therefore can impact natural ventilation in such a way as to make it sub-additive. Several sub-additive superposition models are found in the literature. This paper presents the results of millions of hours of simulations of the physically correct solution, which span a broad range of climates, air leakage and structural conditions. This wide range of data allows for the comparison of three superposition models from the literature and eight new ones. The results show that by using the appropriate model(s) superposition errors can be reduced significantly, from the 20% over-prediction of simple linear addition to 1% or less.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Simplified Methods for Combining Natural and Mechanical Ventilation

In determining ventilation rates, it is often necessary to combine naturally-driven ventilation, such as infiltration, with mechanical systems. Modern calculation methods are sufficiently powerful that this can be done from first principles with time varying flows, but for some purposes simplified methods of combining the mechanical and natural ventilation are required—we call this “superposition”. An example of superposition would be ventilation standards that may pre-calculate some quantities within the body of the standard. When there are balanced mechanical systems, the solution is simple additivity, because a balanced system does not impact the internal pressure of the space. Unbalanced systems, however, change internal pressures and therefore can impact natural ventilation in such a way as to make it sub-additive. Several sub-additive superposition models are found in the literature. This paper presents the results of millions of hours of simulations of the physically correct solution, which span a broad range of climates, air leakage and structural conditions. This wide range of data allows for the comparison of three superposition models from the literature and eight new ones. The results show that by using the appropriate model(s) superposition errors can be reduced significantly, from the 20% over-prediction of simple linear addition to 1% or less.

42 ENGINEERING↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2017

Nature and Consequences of Biological Reductionism for the Immunological Study of Infectious Diseases

Evolution has conserved “economic” systems that perform many functions, faster or better, with less. For example, three to five leukocyte types protect from thousands of pathogens. In order to achieve so much with so little, biological systems combine their limited elements, creating complex structures. Yet, the prevalent research paradigm is reductionist. Focusing on infectious diseases, reductionist and non-reductionist views are here described. Furthermore, the literature indicates that reductionism is associated with information loss and errors, while non-reductionist operations can extract more information from the same data. When designed to capture one-to-many/many-to-one interactions—including the use of arrows that connect pairs of consecutive observations—non-reductionist (spatial–temporal) constructs eliminate data variability from all dimensions, except along one line, while arrows describe the directionality of temporal changes that occur along the line. To validate the patterns detected by non-reductionist operations, reductionist procedures are needed. Integrated (non-reductionist and reductionist) methods can (i) distinguish data subsets that differ immunologically and statistically; (ii) differentiate false-negative from -positive errors; (iii) discriminate disease stages; (iv) capture in vivo, multilevel interactions that consider the patient, the microbe, and antibiotic-mediated responses; and (v) assess dynamics. Integrated methods provide repeatable and biologically interpretable information.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

ION BEAM POLARIZATION DYNAMICS IN THE 8 GEV BOOSTER OF THE JLEIC PROJECT AT JLAB

In the Jefferson Lab’s Electron-Ion Collider (JLEIC) project, an injector of polarized ions into the collider ring is a superconducting 8 GeV booster. Both figure-8 and racetrack booster versions were considered. Our analysis showed that the figure-8 ring configuration allows one to preserve the polarization of any ion species during beam acceleration using only small longitudinal field with an integral less than 0.5 Tm. In the racetrack booster, to pre-serve the polarization of ions with the exception of deu-terons, it suffices to use a solenoidal Siberian snake with a maximum field integral of 30 Tm. To preserve deuteron polarization, we propose to use arc magnets for the race-track booster structure with a field ramp rate of the order of 1 T/s. We calculate deuteron and proton beam polari-zations in both the figure-8 and racetrack boosters includ-ing alignment errors of their magnetic elements using the Zgoubi code.

Kondratenko, A. M.↗