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 1,297 records · Page 72

Comparative Study on the Machine Learning-Based Prediction of Adsorption Energies for Ring and Chain Species on Metal Catalyst Surfaces

Computation of adsorption and transition state energies for a large number of surface intermediates for numerous active site models pose significant computational overhead in computational screening of catalysts. Machine learning (ML) techniques can be used to predict part of these energies. To predict the energies, ML models need to be fed appropriate metal and species descriptors. For complex surface chemistries, the structures of the intermediate species can vary greatly. In this paper, working with the hydrodeoxygenation of succinic acid on six different metal surfaces, we have studied the effect of linear and non-linear ML models used along with pen-and-paper based species descriptors and two categories of metal descriptors on two different categories of intermediate species: chain and ring. More specifically, our computations include the prediction of chain species when trained on only chain species and also when trained on both chain and ring species. Similar computations were performed for predictions of ring species. In each case, results of linear ML models were compared with kernel based non-linear models. Our results indicate that ring species data does not improve the prediction of chain species. Similarly, chain species data does not improve the prediction of ring species. The use of non-linear ML models, however, did help to minimize the prediction errors compared to the linear models. Furthermore, the study also shows that electronic or adsorption energy based metal descriptors along with bond count based species fingerprints can achieve a mean absolute error (MAE) of less than 0.2 eV for complex chain molecules when used with an appropriate machine learning model.

Adsorption↗

Deep operator network surrogate for phase-field modeling of metal grain growth during solidification

A deep operator network (DeepONet) has been constructed that generates accurate representations of phase-field model simulations for evolving two dimensional metal grain morphology growing from melt. These representations serve as lower resolution, computationally efficient stand-ins for quick parameter space exploration of solutions to the the Allen-Cahn equations that dictate the phase-field model simulations. The experimental target for the phase-field model is a uranium casting system cooling a 434 g uranium charge from a maximum temperature of 1400° C at an average rate of 30° C / min , traversing the crystallographic phases of the pure metal. Experimental parameters inform the phase-field model, whose higher resolution computational model solutions are used to train the DeepONet in a given parameter space with the aim of developing a faster, more efficient method for predicting the solidifying metal's microstructure at different potential experimental values. The final DeepONet generates high accuracy, lower resolution predictions with cumulative relative approximation error over all timesteps of less than 0.5%, while ensuring solutions remain within physically feasible ranges. Further, these relative error values are comparable with other state-of-the-art DeepONet models for microstructure evolution, while significantly reducing the amount of training data required. Training a convolutional neural network simultaneously with the DeepONet, enforcing realistic values at the complex metal grain boundaries, and mathematically encoding boundary conditions into the structure of the DeepONet improved prediction accuracy and computational efficiency over a standard DeepONet model.

36 MATERIALS SCIENCE↗

Spatial representation of organic carbon and active-layer thickness of high latitude soils in CMIP5 earth system models

Soil properties such as soil organic carbon (SOC) stocks and active-layer thickness are used in earth system models (F.SMs) to predict anthropogenic and climatic impacts on soil carbon dynamics, future changes in atmospheric greenhouse gas concentrations, and associated climate changes in the permafrost regions. Accurate representation of spatial and vertical distribution of these soil properties in ESMs is a prerequisite for redudng existing uncertainty in predicting carbon-climate feedbacks. We compared the spatial representation of SOC stocks and active-layer thicknesses predicted by the coupled Modellntercomparison Project Phase 5 { CMIP5) ESMs with those predicted from geospatial predictions, based on observation data for the state of Alaska, USA. For the geospatial modeling. we used soil profile observations {585 for SOC stocks and 153 for active-layer thickness) and environmental variables (climate, topography, land cover, and surficial geology types) and generated fine-resolution (50-m spatial resolution) predictions of SOC stocks (to 1-m depth) and active-layer thickness across Alaska. We found large inter-quartile range (2.5-5.5 m) in predicted active-layer thickness of CMIP5 modeled results and small inter-quartile range (11.5-22 kg m-2) in predicted SOC stocks. The spatial coefficient of variability of active-layer thickness and SOC stocks were lower in CMIP5 predictions compared to our geospatial estimates when gridded at similar spatial resolutions (24.7 compared to 30% and 29 compared to 38%, respectively). However, prediction errors. when calculated for independent validation sites, were several times larger in ESM predictions compared to geospatial predictions. Primaly factors leading to observed differences were ( 1) lack of spatial heterogeneity in ESM predictions, (2) differences in assumptions concerning environmental controls, and (3) the absence of pedogenic processes in ESM model structures. Our results suggest that efforts to incorporate these factors in F.SMs should reduce current uncertainties associated with ESM predictions of carbon-climate feedbacks.

Active-layer thickness↗

Bedrock architecture, soil texture, and hyporheic zone characterization combining electrical resistivity and induced polarization imaging

The shallow (<15 m depth) subsurface environment of a short reach along a bedrock stream was investigated with electrical resistivity and induced polarization (IP) to map details of bedrock and soil sediments. The bedrock is mostly comprised of limestone, and is generally resistive compared to the overlying soil. The soil-bedrock interface was determined through a trial and error approach using a sharp boundary feature in the inversion model. The inferred bedrock surface determined from the inversions exhibited undulating patterns with troughs and ridges. A near continuous trough ran alongside the stream within the floodplain, and is suggestive of a paleochannel among other interpretations of this feature. The structure of the electrical resistivity above the bedrock showed small-scale elongated features. The chargeability from the IP method showed larger scale features. High values of chargeability were associated with the sediments in the floodplain, and low values were associated with bedrock, stream, and soil on the elevated banks above the stream. If the chargeability is associated with membrane polarization characteristic of clayey soils, then IP seems to highlight the mere existence of clay, while the resistivity may be more discernable of the relative proportion of clay. Vegetation differences may also explain the chargeability distribution, where parts of the survey with high chargeability had dense pine with no understory making the soils more organically rich.

54 ENVIRONMENTAL SCIENCES↗

Vector-model-supported approach in prostate plan optimization

Lengthy time consumed in traditional manual plan optimization can limit the use of step-and-shoot intensity-modulated radiotherapy/volumetric-modulated radiotherapy (S&S IMRT/VMAT). A vector model base, retrieving similar radiotherapy cases, was developed with respect to the structural and physiologic features extracted from the Digital Imaging and Communications in Medicine (DICOM) files. Planning parameters were retrieved from the selected similar reference case and applied to the test case to bypass the gradual adjustment of planning parameters. Therefore, the planning time spent on the traditional trial-and-error manual optimization approach in the beginning of optimization could be reduced. Each S&S IMRT/VMAT prostate reference database comprised 100 previously treated cases. Prostate cases were replanned with both traditional optimization and vector-model-supported optimization based on the oncologists' clinical dose prescriptions. A total of 360 plans, which consisted of 30 cases of S&S IMRT, 30 cases of 1-arc VMAT, and 30 cases of 2-arc VMAT plans including first optimization and final optimization with/without vector-model-supported optimization, were compared using the 2-sided t-test and paired Wilcoxon signed rank test, with a significance level of 0.05 and a false discovery rate of less than 0.05. For S&S IMRT, 1-arc VMAT, and 2-arc VMAT prostate plans, there was a significant reduction in the planning time and iteration with vector-model-supported optimization by almost 50%. When the first optimization plans were compared, 2-arc VMAT prostate plans had better plan quality than 1-arc VMAT plans. The volume receiving 35 Gy in the femoral head for 2-arc VMAT plans was reduced with the vector-model-supported optimization compared with the traditional manual optimization approach. Otherwise, the quality of plans from both approaches was comparable. Vector-model-supported optimization was shown to offer much shortened planning time and iteration number without compromising the plan quality.

61 RADIATION PROTECTION AND DOSIMETRY↗

THE APPLICATION OF CONTINUOUS WAVELET TRANSFORM BASED FOREGROUND SUBTRACTION METHOD IN 21 cm SKY SURVEYS

We propose a continuous wavelet transform based non-parametric foreground subtraction method for the detection of redshifted 21 cm signal from the epoch of reionization. This method works based on the assumption that the foreground spectra are smooth in frequency domain, while the 21 cm signal spectrum is full of saw-tooth-like structures, thus their characteristic scales are significantly different. We can distinguish them in the wavelet coefficient space easily and perform the foreground subtraction. Compared with the traditional spectral fitting based method, our method is more tolerant to complex foregrounds. Furthermore, we also find that when the instrument has uncorrected response error, our method can also work significantly better than the spectral fitting based method. Our method can obtain similar results with the Wp smoothing method, which is also a non-parametric method, but our method consumes much less computing time.

79 ASTRONOMY AND ASTROPHYSICS↗

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗

The Effect of the Flux Separability Approximation on Multigroup Neutron Transport

The angular dependence of flux-weighted multigroup cross sections is commonly neglected when generating multigroup libraries. The error of this flux separability approximation is typically not isolated from other error sources due to a lack of availability of library generation and corresponding solvers that cannot relax this approximation. These errors can now be isolated and quantified with the availability of a multigroup Monte Carlo transport and multigroup library-generation capability in the OpenMC Monte Carlo transport code. This work will discuss relevant details of the OpenMC implementation, provide an example case useful for detailing the type of errors one can expect from making the flux separability approximation, and end with more realistic problems which show the impact of the approximation and highlight how it can strongly arise from an energy-dependent resonance absorption effect. Since the angle-dependence is intrinsically linked to the energy group structure, these examples also show that relaxing the flux separability approximation with angle-dependent cross sections could be used to reduce either the fine-tuning required to set a multigroup energy structure for a specific reactor type or the number of energy groups required to obtain a desired level of accuracy for a given problem. This trade-off could increase the costs of generating multigroup cross sections, and has the potential to require more memory for storing the multigroup library during the transport calculations, but it can significantly reduce the computational time required since the runtime of a discrete ordinates or method of characteristics neutron transport solver scales roughly linearly with the number of groups.

42 ENGINEERING↗

AMiBA: cluster Sunyaev-Zel'dovich effect observations with the expanded 13-element array

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3 mm to measure the Sunyaev–Zel’dovich effect (SZE) of galaxy clusters at arcminute scales. The first phase of operation—with a compact 7-element array with 0.6 m antennas (AMiBA-7)—observed six clusters at angular scales from 5{sup ′} to 23{sup ′}. Here, we describe the expansion of AMiBA to a 13-element array with 1.2 m antennas (AMiBA-13), its subsequent commissioning, and cluster SZE observing program. The most noticeable changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4 m to 4.8 m, allowing us to sample structures between 2{sup ′} and 10{sup ′}, (2) 13 new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2 m reflectors, and (3) additional correlators and six new receivers. Since the reflectors are co-mounted on and distributed over the entire six-meter CFRP platform, a refined hexapod pointing error model and phase error correction scheme have been developed for AMiBA-13. These effects—entirely negligible for the earlier central close-packed AMiBA-7 configuration—can lead to additional geometrical delays during observations. Our correction scheme recovers at least 80 ± 5% of the point-source fluxes. We, therefore, apply an upward correcting factor of 1.25 to our visibilities to correct for phase decoherence, and a ±5% systematic uncertainty is added in quadrature with our statistical errors. We demonstrate the absence of further systematics with a noise level consistent with zero in stacked uv-visibilities. From the AMiBA-13 SZE observing program, we present here maps of a subset of 12 clusters with signal-to-noise ratios above five. We demonstrate combining AMiBA-7 with AMiBA-13 observations on Abell 1689, by jointly fitting their data to a generalized Navarro–Frenk–White model. Our cylindrically integrated Compton-y values for five radii are consistent with results from the Berkeley-Illinois-Maryland Array, the Owens Valley Radio Observatory, the Sunyaev–Zel’dovich Array, and the Planck Observatory. We also report the first targeted SZE detection toward the optically selected cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗

Theory of Lamb Shift in Muonic Hydrogen

There has been for a while a large discrepancy between the values of the proton charge radius measured by the Lamb shift in muonic hydrogen and by other methods. It has already been clear that theory of muonic hydrogen is reliable at the level of this discrepancy and an error there cannot be a reason for the contradiction. Still the status of theory at the level of the uncertainty of the muonic-hydrogen experiment (which is two orders of magnitude below the discrepancy level) requires an additional clarification. Here, we revisit theory of the 2p − 2s Lamb shift in muonic hydrogen. We summarize all the theoretical contributions in order α{sup 5}m, including pure quantum electrodynamics (QED) ones as well as those which involve the proton-structure effects. Certain enhanced higher-order effects are also discussed. We basically confirm former QED calculations of other authors, present a review of recent calculations of the proton-structure effects, and treat self-consistently higher-order proton-finite-size corrections. We also overview theory of the 2p states. Eventually, we derive a value of the root-mean-square proton charge radius. It is found to be 0.840 29(55) fm, which is slightly different from that previously published in the literature (0.840 87(39) fm [Antognini et al., Science 339, 417 (2013)])

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applying machine learning and quantum chemistry to predict the glass transition temperatures of polymers

Glass transition temperature (T g ) is important for understanding the physical and mechanical properties of a polymer material because it relates to the thermal energy required to transition between a hard glassy state and a soft rubbery one. Over the years, various models have been developed for predicting this thermal property from molecular structure to aid in designing novel polymers in selected classes. This work builds on those efforts by utilizing both machine learning (ML) and quantum chemistry (QC) techniques to develop models that can predict T g values from the molecular structure under different data availability scenarios and for a wide variety of polymer types. For the ML model, a graph convolutional network (GCN) was used to map topological polymer features; this model was trained against a dataset of more than 7500 T g values and resulted in a root mean square error (RMSE) of 38.1 °C. The QC-based regression model was trained on 83 T g values and produced an RMSE of 34.5 °C. In conclusion, this work demonstrated that while both model techniques produce accurate predictions and are suitable for different data availability scenarios, the QC-based regression model offered a more interpretable model framework with significantly less training data.

36 MATERIALS SCIENCE↗

Optimized Quantum Phase Estimation for Simulating Electronic States in Various Energy Regimes

While quantum algorithms for simulation exhibit better asymptotic scaling than their classical counterparts, they currently cannot be implemented on real-world devices. Instead, chemists and computer scientists rely on costly classical simulations of these quantum algorithms. In particular, the quantum phase estimation (QPE) algorithm is among several approaches that have attracted much attention in recent years for its genuine quantum character. However, it is memory-intensive to simulate and intractable for moderate system sizes. Herein, this paper discusses the performance and applicability of QPESIM, a new simulation of the QPE algorithm designed to take advantage of modest computational resources. In particular, we demonstrate the versatility of QPESIM in simulating various electronic states by examining the ground and core-level states of H 2 O. For these states, we also discuss the effect of the active-space size on the quality of the calculated energies. For the high-energy core-level states, we demonstrate that new QPE simulations for active spaces defined by 15 active orbitals significantly reduce the errors in core-level excitation energies compared to earlier QPE simulations using smaller active spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian calibration of irradiated graphite property models under high temperatures

Graphite under high temperatures and irradiation is central to advanced reactors. We develop a Bayesian calibration framework for graphite property models that explicitly represents model-data mismatch via a Gaussian-process discrepancy. The approach propagates uncertainty from parameters, experimental noise, and model form, with a hierarchical variance structure to capture group and cross-group noise. Using two predictive models across five grades (IG-110, NBG-18, PCEA, NBG-17, 2114) and four properties-irradiation-induced dimension change, creep, Young’s modulus change ratio, and coefficient of thermal expansion change ratio-we obtain average predictive-error reductions of 54%, 65%, 17%, and 17% when discrepancy is included. We illustrate engineering impact with a multiphysics model of a very-high-temperature reactor prismatic reflector brick, analyzing stresses under high fluence and temperature. Accounting for model discrepancy markedly improves predictive accuracy and provides a robust basis for reliable graphite component design in advanced reactors.

36 - MATERIALS SCIENCE↗

Kinetic and interaction components of the exact time-dependent correlation potential

The exact exchange-correlation (xc) potential of time-dependent density functional theory has been shown to have striking features. For example, step and peak features are generically found when the system is far from its ground-state, and these depend nonlocally on the density in space and time. We analyze the xc potential by decomposing it into kinetic and interaction components and comparing each with their exact-adiabatic counterparts, for a range of dynamical situations in model one-dimensional two-electron systems. We find that often, but not always, the kinetic contribution is largely responsible for these features that are missed by the adiabatic approximation. The adiabatic approximation often makes a smaller error for the interaction component, which we write in two parts, one being the Coulomb potential due to the time-dependent xc hole. Non-adiabatic features of the kinetic component were also larger than those of the interaction component in cases that we studied when there is negligible step structure. In ground-state situations, step and peak structures arise in cases of static correlation, when more than one determinant is essential to describe the interacting state. We investigate the time-dependent natural orbital occupation numbers and find the corresponding relation between these and the dynamical step is more complex than for the ground-state case.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THEORETICAL CLUES TO THE ULTRAVIOLET DIVERSITY OF TYPE Ia SUPERNOVAE

The effect of metallicity on the observed light of Type Ia supernovae (SNe Ia) could lead to systematic errors as the absolute magnitudes of local and distant SNe Ia are compared to measure luminosity distances and determine cosmological parameters. The UV light may be especially sensitive to metallicity, though different modeling methods disagree as to the magnitude, wavelength dependence, and even the sign of the effect. The outer density structure, {sup 56}Ni, and to a lesser degree asphericity, also impact the UV. We compute synthetic photometry of various metallicity-dependent models and compare to UV/optical photometry from the Swift Ultra-Violet/Optical Telescope. We find that the scatter in the mid-UV to near-UV colors is larger than predicted by changes in metallicity alone and is not consistent with reddening. We demonstrate that a recently employed method to determine relative abundances using UV spectra can be done using UVOT photometry, but we warn that accurate results require an accurate model of the cause of the variations. The abundance of UV photometry now available should provide constraints on models that typically rely on UV spectroscopy for constraining metallicity, density, and other parameters. Nevertheless, UV spectroscopy for a variety of supernova explosions is still needed to guide the creation of accurate models. A better understanding of the influences affecting the UV is important for using SNe Ia as cosmological probes, as the UV light may test whether SNe Ia are significantly affected by evolutionary effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of small variations of speed of sound in optoacoustic tomographic imaging

Purpose: Speed of sound difference in the imaged object and surrounding coupling medium may reduce the resolution and overall quality of optoacoustic tomographic reconstructions obtained by assuming a uniform acoustic medium. In this work, the authors investigate the effects of acoustic heterogeneities and discuss potential benefits of accounting for those during the reconstruction procedure. Methods: The time shift of optoacoustic signals in an acoustically heterogeneous medium is studied theoretically by comparing different continuous and discrete wave propagation models. A modification of filtered back-projection reconstruction is subsequently implemented by considering a straight acoustic rays model for ultrasound propagation. The results obtained with this reconstruction procedure are compared numerically and experimentally to those obtained assuming a heuristically fitted uniform speed of sound in both full-view and limited-view optoacoustic tomography scenarios. Results: The theoretical analysis showcases that the errors in the time-of-flight of the signals predicted by considering the straight acoustic rays model tend to be generally small. When using this model for reconstructing simulated data, the resulting images accurately represent the theoretical ones. On the other hand, significant deviations in the location of the absorbing structures are found when using a uniform speed of sound assumption. The experimental results obtained with tissue-mimicking phantoms and a mouse postmortem are found to be consistent with the numerical simulations. Conclusions: Accurate analysis of effects of small speed of sound variations demonstrates that accounting for differences in the speed of sound allows improving optoacoustic reconstruction results in realistic imaging scenarios involving acoustic heterogeneities in tissues and surrounding media.

60 APPLIED LIFE SCIENCES↗

Damage Pseudo Cross Section Generator

This software is a python script for generating group damage pseudo-cross-sections for specific elements, and tabulating them in forms useful to python or to office-open workbooks (xlsx). The damage pseudo-cross-section is needed for calculating how much neutron radiation damage a structural material receives in a reactor measured in displacements per atom (DPA). The process of generating these pseudo-cross-sections is laborious as it requires running the nuclear data processing code, NJOY, for hundreds of isotopes, and then post-processing the data. This also allows many opportunities for human error to be involved. This script eliminates these problems by fully automating the process.

Gale, MicahD.↗

RELIABILITY BASED DESIGN OF FIXED FOUNDATION WIND TURBINES

Recent analysis of offshore wind turbine foundations using both applicable API and IEC standards show that the total load demand from wind and waves is greatest in wave driven storms. Further, analysis of overturning moment loads (OTM) reveal that impact forces exerted by breaking waves are the largest contributor to OTM in big storms at wind speeds above the operating range of 25 m/s. Currently, no codes or standards for offshore wind power generators have been adopted by the Bureau of Ocean Energy Management Regulation and Enforcement (BOEMRE) for use on the Outer Continental Shelf (OCS). Current design methods based on allowable stress design (ASD) incorporate the uncertainty in the variation of loads transferred to the foundation and geotechnical capacity of the soil and rock to support the loads is incorporated into a factor of safety. Sources of uncertainty include spatial and temporal variation of engineering properties, reliability of property measurements applicability and sufficiency of sampling and testing methods, modeling errors, and variability of estimated load predictions. In ASD these sources of variability are generally given qualitative rather than quantitative consideration. The IEC 61400‐3 design standard for offshore wind turbines is based on ASD methods. Load and resistance factor design (LRFD) methods are being increasingly used in the design of structures. Uncertainties such as those listed above can be included quantitatively into the LRFD process. In LRFD load factors and resistance factors are statistically based. This type of analysis recognizes that there is always some probability of failure and enables the probability of failure to be quantified. This paper presents an integrated approach consisting of field observations and numerical simulation to establish the distribution of loads from breaking waves to support the LRFD of fixed offshore foundations.

17 WIND ENERGY↗