Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 253 records · Page 14

Quality Control of Silicon Sensor Modules for Particle Detectors

The High-Luminosity Large Hadron Collider (HL-LHC) will produce a higher rate of particle collisions than the current Large Hadron Collider (LHC), requiring important upgrades to the Compact Muon Solenoid (CMS) to handle an increased amount of data. An important upgrade is the Phase-2 Outer Tracker Upgrade, which consists of 13,000 silicon sensor modules made of two parallel silicon sensors and readout electronics. These modules undergo careful quality control checks both during and after module assembly to ensure precise and reliable detector performance. This project focuses on precision testing for quality control of silicon sensor modules at Fermilab. Hands-on work includes visual inspection, current-voltage testing, module testing, and ultraviolet (UV) light exposure of modules showing abnormal current-voltage behavior. The ultraviolet exposure process improves the abnormal sensor readout data by placing the selected sensor side of the module directly under the UV light inside a controlled box. In addition to laboratory testing and ultraviolet experiments, I developed a Python-based data tool that connects to a module database and allows selected testing conditions and module information to be retrieved and displayed efficiently. These different testing procedures, experimental processes, and computational tools support the broader goal of identifying module issues and improving modules that will be used in the CMS Outer Tracker Phase-2 Upgrade.

Siddiqui, Hooriya [DuPage Coll.] (ORCID:0009000151↗

Quality Control of Silicon Sensor Modules for Particle Detectors

The High-Luminosity Large Hadron Collider (HL-LHC) will produce a higher rate of particle collisions than the current Large Hadron Collider (LHC), requiring significant upgrades to the Compact Muon Solenoid (CMS) to handle the increased amount of data. An important upgrade is the Phase-2 Outer Tracker Upgrade, which consists of 13,000 silicon sensor modules made of two parallel silicon sensors and readout electronics. These modules undergo careful quality control checks both during and after module assembly to ensure precise and reliable detector performance. This project focuses on precision testing for quality control of silicon sensor modules at Fermilab. Hands-on work includes visual inspection, current-voltage testing, module testing, and ultraviolet (UV) light exposure of modules showing abnormal current-voltage behavior. The ultraviolet exposure process improves the abnormal sensor readout data by placing the selected sensor side of the module directly under the UV light inside a controlled box. In addition to laboratory testing and ultraviolet experiments, I developed a Python-based data tool that connects to a module database and allows selected testing conditions and module information to be retrieved and displayed efficiently. These different testing procedures, experimental processes, and computational tools support the broader goal of identifying module issues and improving modules that will be used in the CMS Outer Tracker Phase-2 Upgrade.

Siddiqui, Hooriya [DuPage Coll.; Fermilab] (ORCID:↗

A learning-based data-driven forecast approach for predicting future reservoir performance

Quantification of the predictive uncertainty of subsurface models has long been investigated. The traditional workflow is to calibrate prior models to match observed data, and then use the posterior models to simulate future system performance. Not only are these procedures computationally expensive, but they also have issues in maintaining geological model constraints during the calibration step. Data space inversion (DSI) was introduced recently to predict future system performance without the iterative history matching or model calibration step. In general, DSI approaches seek to establish a statistical relationship between the observed and forecast variables, as well as to quantify the predictive uncertainty of the forecast variables, by using an ensemble of uncalibrated prior models. Existing DSI approaches all require a number of complex transformation and mapping operations, which may deter their widespread use. Here in this study, we introduce a new and simpler DSI approach, the learning-based, data-driven forecast approach (LDFA), by combining dimension reduction and machine learning techniques to quickly provide accurate forecast results and reliably quantify corresponding uncertainty in the results. Our LDFA framework is demonstrated using two supervised learning algorithms, artificial neural network (ANN) and support vector regression (SVR), on two representative examples from reservoir engineering and geological carbon storage. Results suggest that our approach provides accurate forecast results (e.g., future oil production rate or cumulative injected CO 2 ) and reasonable predictive uncertainty intervals. Our framework is generic and may be applied to other surface and subsurface problems.

42 ENGINEERING↗

Computing partition functions in the one-clean-qubit model

We present a method to approximate partition functions of quantum systems using mixed-state quantum computation. For positive-semidefinite Hamiltonians, our method has an expected running-time that is almost linear in [M/( ε rel Z )] 2 , where M is the dimension of the quantum system, Z is the partition function, and ε rel is the relative precision. It is based on approximations of the exponential operator as linear combinations of certain operators related to block-encoding of Hamiltonians or Hamiltonian evolutions. The trace of each operator is estimated using a standard algorithm in the one-clean-qubit model. For large values of Z , our method may run faster than exact classical methods, whose complexities are polynomial in M . We also prove that a version of the partition function estimation problem within additive error is complete for the so-called DQC1 complexity class, suggesting that our method provides a superpolynomial speedup for certain parameter values. Overall, to attain a desired relative precision, we develop a classical procedure based on a sequence of approximations within predetermined additive errors that may be of independent interest.

97 MATHEMATICS AND COMPUTING↗

Yttrium-90 Radioembolization for Unresectable Standard-chemorefractory Intrahepatic Cholangiocarcinoma: Survival, Efficacy, and Safety Study

To assess the overall survival, efficacy, and safety of radioembolization with yttrium-90 (Y90) for unresectable standard-chemorefractory intrahepatic cholangiocarcinoma (ICC). Patients with unresectable standard-chemorefractory ICC treated with Y90 were studied. Survival was calculated from the date of first Y90 procedure. Tumor response was assessed with the Response Evaluation Criteria in Solid Tumors criteria on follow-up computed tomography or magnetic resonance imaging scans. National Cancer Institute Common Terminology Criteria (NCI CTCAE), version 3, were used for complications. Statistical analysis was performed by the Kaplan-Meier estimator by the log rank test. Nineteen patients underwent a total of 24 resin-based Y90 treatments. Median survival from the time of diagnosis and first Y90 procedure was 752 {+-} 193 [95 % confidence interval (CI) 374-1130] and 345 {+-} 128 (95 % CI 95-595) days, respectively. Median survival with Eastern Cooperative Oncology Group (ECOG) performance status 1 (n = 15) and ECOG performance status 2 (n = 4) was 450 {+-} 190 (95 % CI 78-822) and 345 {+-} 227 (95 % CI 0-790) days, respectively (p = .214). Patients with extrahepatic metastasis (n = 11) had a median survival of 404 {+-} 309 (95 % CI 0-1010) days versus 345 {+-} 117 (95 % CI 115-575) days for patients without metastasis (n = 8) (p = .491). No mortality was reported within 30 days from first Y90 radioembolization. One patient developed grade 3 thrombocytopenia as assessed by NCI CTCAE. Fatigue and transient abdominal pain were observed in 4 (21 %) and 6 (32 %) patients, respectively. Y90 radioembolization is effective for unresectable standard-chemorefractory ICC.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An empirical method for geometric calibration of a photon counting detector-based cone beam CT system

BACKGROUND: Geometric calibration is essential in developing a reliable computed tomography (CT) system. It involves estimating the geometry under which the angular projections are acquired. Geometric calibration of cone beam CTs employing small area detectors, such as currently available photon counting detectors (PCDs), is challenging when using traditional-based methods due to detectors’ limited areas. OBJECTIVE: This study presented an empirical method for the geometric calibration of small area PCD-based cone beam CT systems. METHODS: Unlike the traditional methods, we developed an iterative optimization procedure to determine geometric parameters using the reconstructed images of small metal ball bearings (BBs) embedded in a custom-built phantom. An objective function incorporating the sphericities and symmetries of the embedded BBs was defined to assess performance of the reconstruction algorithm with the given initial estimated set of geometric parameters. The optimal parameter values were those which minimized the objective function. The TIGRE toolbox was employed for fast tomographic reconstruction. To evaluate the proposed method, computer simulations were carried out using various numbers of spheres placed in various locations. Furthermore, efficacy of the method was experimentally assessed using a custom-made benchtop PCD-based cone beam CT. RESULTS: Computer simulations validated the accuracy and reproducibility of the proposed method. The precise estimation of the geometric parameters of the benchtop revealed high-quality imaging in CT reconstruction of a breast phantom. Within the phantom, the cylindrical holes, fibers, and speck groups were imaged in high fidelity. The CNR analysis further revealed the quantitative improvements of the reconstruction performed with the estimated parameters using the proposed method. CONCLUSION: Apart from the computational cost, we concluded that the method was easy to implement and robust.

Instruments & Instrumentation↗

Two-tier model reduction of viscoelastically damped finite element models

Transient simulations of linear viscoelastically damped structures necessitate critical computational resources to directly integrate the full-order finite element model with time-stepping algorithms. Traditional modal reduction techniques are not directly applicable to these systems since viscoelastic materials depend on time and frequency. A more appropriate reduction basis is obtained from the nonlinear, complex eigenvalue problem, whose eigenvectors capture the appropriate kinematics and enable frequency-based mode selection; unfortunately, the computational cost is prohibitive for computing these modes from large-scale engineering models. To address this shortcoming, this work proposes a novel two-tier reduction procedure to reduce the upfront cost of solving the complex, nonlinear eigenvalue problem. The first reduction step reduces the full-order model with real mode shapes linearized about various centering frequencies to capture the kinematics over a full range of viscoelastic material behavior (glassy, rubbery, and glass-transition zones). This tier-one reduction preserves time-temperature superposition and allows the equations to depend parametrically on operating temperature. The second-level reduction then solves the complex, nonlinear eigenmode solutions in the tier-one reduced space about a fixed temperature to further reduce the equations-of-motion. The method is demonstrated on a cantilevered sandwich plate to display its accuracy and efficiency in comparison to full-order model predictions.

97 MATHEMATICS AND COMPUTING↗

A Bloch decomposition-based stochastic Galerkin method for quantum dynamics with a random external potential

In this paper, we consider the numerical solution of the one-dimensional Schrödinger equation with a periodic lattice potential and a random external potential. This is an important model in solid state physics where the randomness results from complicated phenomena that are not exactly known. Here we generalize the Bloch decomposition-based time-splitting pseudospectral method to the stochastic setting using the generalized polynomial chaos with a Galerkin procedure so that the main effects of dispersion and periodic potential are still computed together. We prove that our method is unconditionally stable and numerical examples show that it has other nice properties and is more efficient than the traditional method. Finally, we give some numerical evidence for the well-known phenomenon of Anderson localization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗

Kernel-based global sensitivity analysis obtained from a single data set

Results from global sensitivity analysis (GSA) often guide the understanding of complicated input–output systems. Kernel-based GSA methods have recently been proposed for their capability of treating a broad scope of complex systems. In this paper, we develop a new set of kernel GSA tools when only a single set of input–output data is available. Three key advances are made: (1) A new numerical estimator is proposed that demonstrates an empirical improvement over previous procedures. (2) A computational method for generating inner statistical functions from a single data set is presented. (3) A theoretical extension is made to define conditional sensitivity indices, which reveal the degree that the inputs carry shared information about the output when inherent input–input correlations are present. Utilizing these conditional sensitivity indices, a decomposition is derived for the output uncertainty based on what is called the optimal learning sequence of the input variables, which remains consistent when correlations exist between the input variables. Further, while these advances cover a range of GSA subjects, a common single data set numerical solution is provided by a technique known as the conditional mean embedding of distributions. The new methodology is implemented on benchmark systems to demonstrate the provided insights.

42 ENGINEERING↗

Evidence Theory Representations for Loss of Assured Safety in Weak Link/Strong Link Systems.

The following topics are considered in this presentation: (i) Overview of evidence theory, (ii) Representation of loss of assured safety (LOAS) with evidence theory for a 1 SL, 1 WL system, (iii) Description of 2 SLs and 1 WL used for illustration, (iv) Plausibility and belief for LOAS and associated sampling-based verification calculations for a 2 SL, 1 WL system, (iv) Plausibility and belief for margins associated with LOAS for a 2 SL, 1 WL system, (v) Plausibility and belief for LOAS for a 2 SL, 2 WL system, (vi) Incorporation of evidence spaces for link temperature curves into LOAS calculations, (vii) Plausibility and belief for LOAS for WL/SL systems with SL subsystems, and (viii) Sampling-based procedures for the estimation of plausibility and belief. 2

42 ENGINEERING↗

Assessing the utility of phase-space-localized basis functions: Exploiting direct product structure and a new basis function selection procedure

In this paper we show that it is possible to use an iterative eigensolver in conjunction with Halverson and Poirier’s symmetrized Gaussian (SG) basis [T. Halverson and B. Poirier, J. Chem. Phys. 137, 224101 (2012)] to compute accurate vibrational energy levels of molecules with as many as five atoms. This is done, without storing and manipulating large matrices, by solving a regular eigenvalue problem that makes it possible to exploit direct-product structure. These ideas are combined with a new procedure for selecting which basis functions to use. The SG basis we work with is orders of magnitude smaller than the basis made by using a classical energy criterion. We find significant convergence errors in previous calculations with SG bases. For sum-of-product Hamiltonians, SG bases large enough to compute accurate levels are orders of magnitude larger than even simple pruned bases composed of products of harmonic oscillator functions.

97 MATHEMATICS AND COMPUTING↗

Photoconductivities from band states and a dissipative electron dynamics: Si(111) without and with adsorbed Ag clusters

A new general computational procedure is presented to obtain photoconductivities starting from atomic structures, combining ab initio electronic energy band states with populations from density matrix theory, and implemented for a specific set of materials based on Si crystalline slabs and their nanostructured surfaces without and with adsorbed Ag clusters. The procedure accounts for charge mobility in semiconductors in photoexcited states, and specifically electron and hole photomobilities at Si(111) surfaces with and without adsorbed Ag clusters using ab initio energy bands and orbitals generated from a generalized gradient functional, however with excited energy levels modified to provide correct bandgaps. Photoexcited state populations for each band and carrier type were generated using steady state solution of a reduced density matrix which includes dissipative medium effects. The present calculations provide photoexcited electronic populations and photoinduced mobilities resulting from applied electric fields and obtained from the change of driven electron energies with their electronic momentum. Extensive results for Si slabs with 8 layers, without and with adsorbed Ag clusters, show that the metal adsorbates lead to substantial increases in the photomobility and photoconductivity of electrons and holes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-G-IeP3-01: Better Kerma-Area-Product (KAP) Estimation Using the System Parameters in Radiography and Fluoroscopy

Purpose: Accurate values for Kerma-Area-Product (KAP) are needed for patient dosimetry and quality control for exams utilizing radiographic and/or fluoroscopic imaging. The KAP measured using a typical direct KAP meter built with parallel-plate transmission ionization chamber is not precise and depends on the energy spectrum of diagnostic x-rays. This study compared the accuracy and reproducibility of KAP derived from system parameters with values measured with a direct KAP meter. Methods: IEC tolerance for displayed KAP is specified up to ± 35% above 2.5 Gy-cm{sup 2} and manufacturer’s specifications are typically ± 25%. KAP values from the direct KAP meter drifts with time leading to replacement or re-calibration. More precise and consistent KAP is achievable utilizing a database of known radiation output for various system parameters. The integrated KAP meter was removed from a radiography system. A total of 48 measurements of air kerma were acquired at x-ray tube potential from 40 to 150 kVp with 10 kVp increment using ion chamber type external dosimeter at free-in-air geometry for four different types of filter combinations following the manufacturer’s service procedure. These data were used to create updated correction factors that determine air kerma computationally for given system parameters. Results of calculated KAP were evaluated against results using a calibrated ion chamber based dosimeter and a computed radiography imaging plate to measure x-ray field size. Results: The accuracy of calculated KAP from the system parameters was better within 4% deviation in all diagnostic x-ray tube potentials tested from 50 to 140 kVp. In contrast, deviations of up to 25% were measured from KAP displayed from the direct KAP meter. Conclusion: The “calculated KAP” approach provides the nominal advantage of improved accuracy and precision of displayed KAP as well as reduced cost of calibrating or replacing integrated KAP meters.

60 APPLIED LIFE SCIENCES↗

APPLICATION OF THE FAVOR-OCI FRACTURE MECHANICS COMPUTER PROGRAM TO ASME CODE SECTION XI, IWB-3610 FLAW ACCEPTANCE CRITERIA EVALUATIONS

This paper presents an overview of added features in a new version of the FAVOR (Fracture Analysis of Vessels Oak Ridge) computer code called FAVOR-OCI. The original FAVOR code was developed at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) under the sponsorship of the US Nuclear Regulatory Commission (NRC). FAVOR is applied by US and international nuclear power industries to perform deterministic and probabilistic fracture mechanics analyses of commercial nuclear reactor pressure vessels (RPVs). Appli-cations of FAVOR are focused on insuring that the structural integrity of aging, and increasingly embrittled, RPVs is main-tained throughout their licensed service life. Based on the final ORNL release of FAVOR, v16.1, FAVOR-OCI extends existing deterministic features of FAVOR while preserving all prev-iously-existing probabilistic capabilities of FAVOR.The objective of this paper is to describe new deterministic features in FAVOR-OCI that can be applied to analytical evaluations of planar flaws. These evaluations are consistent with the acceptance criteria of ASME Code, Section XI, Article IWB-3610, including Subarticles IWB-3611 (acceptance based on flaw size) and IWB-3612 (acceptance based on applied stress intensity factor). The linear elastic fracture mechanics (LEFM) capabilities of FAVOR-OCI also incorporate the analytical procedures presented in the Nonmandatory Appendix A, Analysis of Flaws, Article A-3000, Method of KI Determination, for both surface and subsurface (embedded) flaws.The paper describes a computational methodology for determining critical values of fracture-related parameters that satisfy ASME Code Section XI acceptance criteria for flaws exposed to multiple thermal-hydraulic transients. These com-pute-intensive analyses can be carried out with a single execution of FAVOR-OCI. The new methodology determines critical values by solving for either the point of tangency or point of intersection between applied KI versus time histories and a user-selected fracture toughness material property (e.g., ASME KIc, FAVOR Weibull KIc, or Master Curve Weibull KJc) for surface or subsurface flaws. Situations where warm prestress conditions apply can also be addressed. The paper highlights a need for this new capability via applications to a recent independent review of safety cases for RPVs in two Belgian nuclear power plants (NPPs). That review required ASME Section XI assessments of several thousand embedded, quasi-laminar flaws in the wall of each RPV. Analysis results provided by the new capability contributed to the technical bases compiled from several sources by the Belgian nuclear regulatory agency (FANC) and eventually used by FANC to justify the restart of these NPPs.

Dickson, Terry L.↗

An Online Generalized Multiscale Discontinuous Galerkin Method (GMsDGM) for Flows in Heterogeneous Media

Abstract Offline computation is an essential component in most multiscale model reduction techniques. However, there are multiscale problems in which offline procedure is insufficient to give accurate representations of solutions, due to the fact that offline computations are typically performed locally and global information is missing in these offline information. To tackle this difficulty, we develop an online local adaptivity technique for local multiscale model reduction problems. We design new online basis functions within Discontinuous Galerkin method based on local residuals and some optimally estimates. The resulting basis functions are able to capture the solution efficiently and accurately, and are added to the approximation iteratively. Moreover, we show that the iterative procedure is convergent with a rate independent of physical scales if the initial space is chosen carefully. Our analysis also gives a guideline on how to choose the initial space. We present some numerical examples to show the performance of the proposed method.

Physics↗

Efficient generation of sum-of-products representations of high-dimensional potential energy surfaces based on multimode expansions

The transformation of multi-dimensional potential energy surfaces (PESs) from a grid-based multimode representation to an analytical one is a standard procedure in quantum chemical programs. Within the framework of linear least squares fitting, a simple and highly efficient algorithm is presented, which relies on a direct product representation of the PES and a repeated use of Kronecker products. It shows the same scalings in computational cost and memory requirements as the POTFIT approach. In comparison to customary linear least squares fitting algorithms, this corresponds to a speed-up and memory saving by several orders of magnitude. Different fitting bases are tested, namely, polynomials, B-splines, and distributed Gaussians. Benchmark calculations are provided for the PESs of a set of small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gradient boosting machine for modeling the energy consumption of commercial buildings

Accurate savings estimations are important to promote energy efficiency projects and demonstrate their cost-effectiveness. The increasing presence of advanced metering infrastructure (AMI) in commercial buildings has resulted in a rising availability of high frequency interval data. These data can be used for a variety of energy efficiency applications such as demand response, fault detection and diagnosis, and heating, ventilation, and air conditioning (HVAC) optimization. This large amount of data has also opened the door to the use of advanced statistical learning models, which hold promise for providing accurate building baseline energy consumption predictions, and thus accurate saving estimations. The gradient boosting machine is a powerful machine learning algorithm that is gaining considerable traction in a wide range of data driven applications, such as ecology, computer vision, and biology. In the present work an energy consumption baseline modeling method based on a gradient boosting machine was proposed. To assess the performance of this method, a recently published testing procedure was used on a large dataset of 410 commercial buildings. The model training periods were varied and several prediction accuracy metrics were used to evaluate the model's performance. The results show that using the gradient boosting machine model improved the R-squared prediction accuracy and the CV(RMSE) in more than 80 percent of the cases, when compared to an industry best practice model that is based on piecewise linear regression, and to a random forest algorithm.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗