SEARCH · Search NASA

Results for “model 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Accelerated fission rate irradiation design, pre-irradiation characterization, and adaptation of conventional PIE methods for U-10Mo and U-17Mo

Metallic U alloys have high U density and thermal conductivity and thus have been explored since the beginning of nuclear power research. Alloys of U with modest amounts of Mo, such as U-10 wt % Mo (U-10Mo), are of particular interest because the γ-U crystal structure in this alloying addition shows prolonged stability in reactor service. Historically, radiation data on U-10Mo fuels were collected in Na fast reactors or lower temperature research reactor conditions, but little is known about irradiation behavior, particularly swelling and creep, at irradiation temperatures between 250 and 500°C. This work discusses the methodology and pre-irradiation characterization results from a U-Mo irradiation campaign performed in the High Flux Isotope Reactor at Oak Ridge National Laboratory. U-10Mo and U-17Mo samples irradiations are being completed at temperatures ranging from 250 to 500°C to three targeted fission densities between 2 × 10 20 and 1.5 × 10 21 fissions per cubic centimeter. Swelling measurement of the specimen sizes studied here required development and assessment of new methods for volume determination before and after irradiation. Laser profilometry and X-ray computation tomography (XCT) were used to provide preirradiation characterization of samples to determine the error and applicability of each to determine swelling following irradiation. These outcomes are contextualized through use of BISON simulations performed to assess the predicted expansion of U-Mo fuels subjected to the irradiation conditions of this work. Use of existing BISON fuel performance models predicted a maximum of 7% swelling under the irradiation conditions of this study. Pre-irradiation characterization revealed the as-cast U-Mo fuel samples were uniformly large-grained fully cubic U crystals with small U-C/N bearing precipitates and pores distributed throughout. Samples were found to contain a bulk porosity between .4 and 3% because of the casting process. Local porosity in areas far from large, interconnected pores was found by Slice-and-View to be under .2%. Nanometer-sized precipitates rich in C and N were identified in all samples, likely because of impurities during the fabrication process. Dendritic bands were also observed throughout the samples. These bands were characterized by variable Mo content that deviated from the overall Mo content by 2–3 wt %. No other microstructural features were correlated to these bands. Mechanical properties were found to be slightly strengthened compared to literature reports of bulk U-Mo fuels due to the nano-scale precipitates throughout the sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗

DOE OSTI2022

Soil-related developments of the Biome-BGCMuSo v6.2 terrestrial ecosystem model

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from -0.047 to -0.007 m 3 m -3 for SWC and from -0.68 to -0.2 mm d -1 for soil evaporation). Simulations related to nitrogen balance and soil CO 2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N 2 O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N 2 O efflux simulation bias was between -0.13 and -0.1 mgNm -2 d -1 , and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO 2 efflux simulations bias was 0.04–0.17 gCm -2 d -1 , while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

58 GEOSCIENCES↗

DOE OSTI2020

Examination of the inversion of isobaric analogue states in mirror nuclei

A recent study by Hoff and collaborators [Nature (London) 580, 52 (2020)] presented evidence that the ground-state spin of [Math Processing Error] is different from that of its mirror, [Math Processing Error], likely due to an inversion of the ground- and first-excited states. Here, we assess the likelihood of such an inversion arising from normal isospin-symmetry breaking and, more broadly, whether such phenomena challenge our understanding of charge-symmetry-breaking forces in atomic nuclei. By placing the result within the context of previous experimental and theoretical work we demonstrate that this inversion lies entirely within the bounds of normal isospin-symmetry-breaking behavior. We further note that, in the context of isospin, neither level inversions nor the nuclear ground state hold any special significance.

59 ≤ A ≤ 89↗

DOE OSTI2023

A General Spatiotemporal Imputation Framework for Missing Sensor Data

Many applications from precision agriculture, environmental monitoring and transportation networks rely on data collected across space and time over a large geographic area. Missing data poses a significant challenge for any data-driven inference and control tasks. Data imputation or the estimation of missing data can help fill these gaps by utilizing inherent spatial relationships and temporal patterns. A variety of spatiotemporal imputation models have been developed to address missing data in spatiotemporal datasets. However, these classical methods rely on the assumption that the underlying data follows a smooth trend and fail to provide accurate estimates when there is a large number of missing points in the data. Even though there are machine learning driven tensor completion approaches such as convolutional neural network based tensor completion (CoSTCo) that capture the non-linear relationships in the dataset, the transductive nature makes the algorithm less scalable. Thus, existing approaches for estimating the missing information do not effectively capture all dimensions of the spatiotemporal data structure, resulting in erroneous predictions and poor performance. The main contributions of this paper are: (1) We propose a novel inductive framework (G-LSTM) for missing data imputation that integrates a graph neural network with LSTMs to effectively capture both spatial and temporal dependencies. (2) Experimental results on a traffic dataset demonstrate that the proposed GNN integrated with an LSTM framework achieves improved imputation and maintains steady performance even when there are extreme missing conditions in comparison with the state-of-the-art imputation framework (i.e, CoSTCo). (3) The simulation results on a traffic network show up to 69% reduction in mean absolute error and 61% reduction in root mean square error when compared to CoSTCo.

Tharzeen, Aabila↗

DOE OSTI2015

SU-E-T-502: Initial Results of a Comparison of Treatment Plans Produced From Automated Prioritized Planning Method and a Commercial Treatment Planning System

Purpose We developed an automated treatment planning system based on a hierarchical goal programming approach. To demonstrate the feasibility of our method, we report the comparison of prostate treatment plans produced from the automated treatment planning system with those produced by a commercial treatment planning system. Methods In our approach, we prioritized the goals of the optimization, and solved one goal at a time. The purpose of prioritization is to ensure that higher priority dose-volume planning goals are not sacrificed to improve lower priority goals. The algorithm has four steps. The first step optimizes dose to the target structures, while sparing key sensitive organs from radiation. In the second step, the algorithm finds the best beamlet weight to reduce toxicity risks to normal tissue while holding the objective function achieved in the first step as a constraint, with a small amount of allowed slip. Likewise, the third and fourth steps introduce lower priority normal tissue goals and beam smoothing. We compared with prostate treatment plans from Memorial Sloan Kettering Cancer Center developed using Eclipse, with a prescription dose of 72 Gy. A combination of liear, quadratic, and gEUD objective functions were used with a modified open source solver code (IPOPT). Results Initial plan results on 3 different cases show that the automated planning system is capable of competing or improving on expert-driven eclipse plans. Compared to the Eclipse planning system, the automated system produced up to 26% less mean dose to rectum and 24% less mean dose to bladder while having the same D95 (after matching) to the target. Conclusion We have demonstrated that Pareto optimal treatment plans can be generated automatically without a trial-and-error process. The solver finds an optimal plan for the given patient, as opposed to database-driven approaches that set parameters based on geometry and population modeling.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

DOE OSTI2016

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

An extended radar relative calibration adjustment (eRCA) technique for higher-frequency radars and range–height indicator (RHI) scans

Abstract. This study extends the relative calibration adjustment technique for calibration of weather radars to higher-frequency radars as well as range–height indicator (RHI) scans. The calibration of weather radars represents one of the most dominant sources of error for their use in a variety of fields including quantitative precipitation estimation and model comparisons. While most weather radars are routinely calibrated, the frequency of calibration is often less than required, resulting in miscalibrated time periods. While full absolute calibration techniques often require the radar to be taken offline for a period of time, there have been online calibration techniques discussed in the literature. The relative calibration adjustment (RCA) technique uses the statistics of the ground clutter surrounding the radar as a monitoring source for the stability of calibration but has only been demonstrated to work at S- and C-band for plan-position indicator (PPI) scans at a constant elevation. In this work the RCA technique is modified to work with higher-frequency radars, including Ka-band cloud radars. At higher frequencies the properties of clutter can be much more variable. This work introduces an extended clutter selection procedure that incorporates the temporal stability of clutter and helps to improve the operational stability of RCA for relatively higher-frequency radars. The technique is also extended to utilize range–height scans from radars where the elevation is varied rather than the azimuth. These types of scans are often utilized with research radars to examine the vertical structure of clouds. The newly extended technique (eRCA) is applied to four Department of Energy Atmospheric Radiation Measurement (DOE ARM) weather radars ranging in frequency from C- to Ka-band. Cross comparisons of three co-located radars with frequencies C, X, and Ka at the ARM Cloud, Aerosol, and Complex Terrain Interactions (CACTI) site show that the technique can determine changes in calibration. Using an X-band radar at the ARM Eastern North Atlantic (ENA) site, we show how the technique can be modified to be more resilient to clutter fields that show increased variability, in this case due to sea clutter. The results show that this technique is promising for a posteriori data calibration and monitoring.

47 OTHER INSTRUMENTATION↗

DOE OSTI2018

Parton Dynamics Inferred from High-Mass Drell-Yan Dimuons Induced by 120 GeV p+D Interactions

Fermilab Experiment 906/SeaQuest (E906/SeaQuest) is the latest in a well established tradition of studying leptoproduction from the annihilation of a quark and anti-quark, known as the Drell-Yan process. The broad goal of E906/SeaQuest is measuring various properties of nucleon structure in order to learn more about quarks and Quantum Chromodynamics (QCD), the mathematical description of the strong force. The present work investigated violations of the Lam-Tung relation between virtual photon polarization and quark and lepton angular momentum. The violation of Lam-Tung can be explained as the signature of quark-nucleon spin-orbit coupling through the use of the Transverse-Momentum-Dependent (TMD) framework, which assumes that the initial transverse momentum of quarks is smaller than the hard scattering scale, but also non-negligible. An analysis of the angular moments in Drell-Yan collected by E906/SeaQuest was performed with four different configurations in order to estimate the systematic errors attributed to each correction. After correction for background and error propagation, the final extraction of the azimuthal moment excluding contributions from the trigger was ν = 0.151 ± 0.88(stat.) ± 0.346(syst.) at an average transverse momentum of 0.87 ± 0.50 GeV/c and an average dimuon mass of 5.48 ± 0.70 GeV. In the future, the magnitude of the systematic errors on the extraction could potentially be reduced by improving the quality of the trigger efficiency calculation, improving the intensity dependent event reconstruction efficiency, considering the changes in acceptance due to a beam shift relative to the E906/SeaQuest spectrometer, and improving the modeling of background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Baryon acoustic oscillations from the complete SDSS-III Lyα-quasar cross-correlation function at z = 2.4

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Lyα-forest flux-transmission at a mean redshift z = 2.40. The measurement uses the complete SDSS-III data sample: 168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level allowing a more accurate measurement of this cross-correlation. We also developed the first simulations of the cross-correlation allowing us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios D H (z = 2.40)=r d = 9.01 ± 0.36 and D M (z = 2.40)=r d = 35.7 ±1.7, where the errors include marginalization over the non-linear velocity of quasars and the metal - quasar cross-correlation contribution, among other effects. These results are within 1.8σ of the prediction of the flat-ΛCDM model describing the observed CMB anisotropies.We combine this study with the Lyα-forest auto-correlation function (Bautista et al. 2017), yielding D H (z = 2.40)=r d = 8.94 ± 0.22 and D M (z = 2.40)=r d = 36.6 ± 1.2, within 2.3σ of the same flat-ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)

The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent density functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the determination of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along molecular dynamics trajectories for flexible systems or the calculation of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semi-empirical tight binding (TB) procedure similar to the well established self-consistent-charge density functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid density functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semi-empirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended atomic orbital basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the determination of atomic charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of standard time-dependent density functional theory (0.3-0.5 eV average error). An easily extendable parametrization based on coupled-cluster and density functional computed reference data for the elements H–Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semi-empirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate determination of atomic charges as well as a more accurate calculation scheme of dipole rotatory strengths within the Tamm-Dancoff approximation is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Momentum power spectrum of SDSS galaxies by massE cosmic ruler: 2.1×improvement in measure of growth rate

ABSTRACT Peculiar motion of galaxies probes the structure growth in the universe. In this study, we employ the galaxy stellar mass-binding energy (massE) relation with only two nuisance parameters to build the largest peculiar-velocity (PV) catalogue to date, consisting of 229 890 ellipticals from the main galaxy sample (MGS) of the Sloan Digital Sky Survey (SDSS). We quantify the distribution of the massE-based distances in individual narrow redshift bins (dz = 0.005), and then estimate the PV of each galaxy based on its offset from the Gaussian mean of the distribution. As demonstrated with the Uchuu-SDSS mock data, the derived PV and momentum power spectra are insensitive to accurate calibration of the massE relation itself, enabling measurements out to a redshift of 0.2, well beyond the current limit of z = 0.1 using other galaxy scaling laws. We then measure the momentum power spectrum and demonstrate that it remains almost unchanged if varying significantly the redshift bin size within which the distance is measured, as well as the intercept and slope of the massE relation, respectively. By fitting the spectra using the perturbation theory model with four free parameters, fσ8 is constrained to fσ8 = 0.459$^{+0.068}_{-0.069}$ over Δz = 0.02–0.2, 0.416$^{+0.074}_{-0.076}$ over Δz = 0.02–0.1, and 0.526$^{+0.133}_{-0.148}$ over Δz = 0.1–0.2. The error of fσ8 is 2.1 times smaller than that by the redshift space distortion (RSD) of the same sample. A Fisher matrix forecast illustrates that the constraint on fσ8 from the massE-based PV can potentially exceed that from the stage-IV RSD in late universe (z<0.5).

Shi, Yong (ORCID:0000000286146275)↗

DOE OSTI2015

A Survey of Techniques for Modeling and Improving Reliability of Computing Systems

Recent trends of aggressive technology scaling have greatly exacerbated the occurrences and impact of faults in computing systems. This has made `reliability' a first-order design constraint. To address the challenges of reliability, several techniques have been proposed. In this study, we provide a survey of architectural techniques for improving resilience of computing systems. We especially focus on techniques proposed for microarchitectural components, such as processor registers, functional units, cache and main memory etc. In addition, we discuss techniques proposed for non-volatile memory, GPUs and 3D-stacked processors. To underscore the similarities and differences of the techniques, we classify them based on their key characteristics. We also review the metrics proposed to quantify vulnerability of processor structures. Finally, we believe that this survey will help researchers, system-architects and processor designers in gaining insights into the techniques for improving reliability of computing systems.

3D-stacked processors↗

DOE OSTI2015

WE-AB-204-05: Harmonizing PET/CT Quantification in Multicenter Studies: A Case Study

Purpose: To present the implementation of a strategy to harmonize FDG PET/CT quantification (SUV), performed with different scanner models and manufacturers. Methods: The strategy was based on Boellaard (2011) and EARL FDG-PET/CT accreditation program, that propose quality control measurements for harmonizing scanner performance. A NEMA IEC Body phantom study was performed using four different devices: PHP-1 (Gemini TF Base, Philips); PHP-2 (Gemini GXL, Philips); GEH (Discovery 600, General Electric); SMS (Biograph Hi-Rez 16, Siemens). The SUV Recovery Coefficient (RC) was calculated using the clinical protocol and other clinically relevant reconstruction parameters. The most appropriate reconstruction parameters (MARP) for SUV harmonization, in each scanner, are those which achieve EARL harmonizing standards. They were identified using the lowest root mean square errors (RMSE). To evaluate the strategy’s effectiveness, the Maximum Differences (MD) between the clinical and MARP RC values were calculated. Results: The reconstructions parameters that obtained the lowest RMSE are: FBP 5mm (PHP-1); LOR-RAMLA 2i0.008l (PHP-2); VuePointHD 2i32s10mm (GEH); and FORE+OSEM 4i8s6mm (SMS). Thus, to ensure that quantitative PET image measurements are interchangeable between these sites, images must be reconstructed with the above-mentioned parameters. Although, a decoupling between the best image for PET/CT qualitative analysis and the best image for quantification studies was observed. The MD showed that the strategy was effective in reducing the variability of SUV quantification for small structures (<17mm). Conclusion: The harmonization strategy of the SUV quantification implemented with these devices was effective in reducing the variability of small structures quantification, minimizing the inter-scanner and inter-institution differences in quantification. However, it is essential that, in addition to the harmonization of quantification, the standardization of the methodology of patient preparation must be maintained, in order to minimize the SUV variability due to biological factors. Financial support by CAPES.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

A Sheaf Theoretical Approach to Uncertainty Quantification of Heterogeneous Geolocation Information

Integration of multiple, heterogeneous sensors is a challenging problem across a range of applications. Prominent among these are multi-target tracking, where one must combine observations from different sensor types in a meaningful and efficient way to track multiple targets. Because different sensors have differing error models, we seek a theoretically justified quantification of the agreement among ensembles of sensors, both overall for a sensor collection, and also at a fine-grained level specifying pairwise and multi-way interactions among sensors. We demonstrate that the theory of mathematical sheaves provides a unified answer to this need, supporting both quantitative and qualitative data. Furthermore, the theory provides algorithms to globalize data across the network of deployed sensors, and to diagnose issues when the data do not globalize cleanly. We demonstrate and illustrate the utility of sheaf-based tracking models based on experimental data of a wild population of black bears in Asheville, North Carolina. A measurement model involving four sensors deployed among the bears and the team of scientists charged with tracking their location is deployed. This provides a sheaf-based integration model which is small enough to fully interpret, but of sufficient complexity to demonstrate the sheaf’s ability to recover a holistic picture of the locations and behaviors of both individual bears and the bear-human tracking system. A statistical approach was developed in parallel for comparison, a dynamic linear model which was estimated using a Kalman filter. This approach also recovered bear and human locations and sensor accuracies. When the observations are normalized into a common coordinate system, the structure of the dynamic linear observation model recapitulates the structure of the sheaf model, demonstrating the canonicity of the sheaf-based approach. However, when the observations are not so normalized, the sheaf model still remains valid.

97 MATHEMATICS AND COMPUTING↗