Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational simulation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Surveying the LLNL WSC/CP DevOps Landscape - FY23 DevOps L2: Advanced Simulation and Computing (ASC) L2 Milestone 8542, "Spack Utilization in IC Code Projects”

This L2 milestone is a multifaceted review of DevOps practices across Weapon Simulation and Computing/Computational Physics (WSC/CP) Program, which includes Weapons, ICF, and Engineering simulation codes and supporting libraries. It complements the FY22 Spack/MARBL L2: Advanced Simulation and Computing (ASC) L2 Milestone 7904, “Workflow Portability Across Cloud Services”. The overall goal is to develop a path forward for improving DevOps practices in WSC/CP which will increase developer productivity, improve software quality, and speed up our ability to deliver capabilities to end users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

SU-E-I-02: A Framework to Perform Batch Simulations of Computational Voxel Phantoms to Study Organ Doses in Computed Tomography Using a Commercial Monte Carlo Software Package

Purpose: ImpactMC (CT Imaging, Erlangen, Germany) is a Monte Carlo (MC) software package that offers a GPU enabled, user definable and validated method for 3D dose distribution calculations for radiography and Computed Tomography (CT). ImpactMC, in and of itself, offers limited capabilities to perform batch simulations. The aim of this work was to develop a framework for the batch simulation of absorbed organ dose distributions from CT scans of computational voxel phantoms. Methods: The ICRP 110 adult Reference Male and Reference Female computational voxel phantoms were formatted into compatible input volumes for MC simulations. A Matlab (The MathWorks Inc., Natick, MA) script was written to loop through a user defined set of simulation parameters and 1) generate input files required for the simulation, 2) start the MC simulation, 3) segment the absorbed dose for organs in the simulated dose volume and 4) transfer the organ doses to a database. A demonstration of the framework is made where the glandular breast dose to the adult Reference Female phantom, for a typical Chest CT examination, is investigated. Results: A batch of 48 contiguous simulations was performed with variations in the total collimation and spiral pitch. The demonstration of the framework showed that the glandular dose to the right and left breast will vary depending on the start angle of rotation, total collimation and spiral pitch. Conclusion: The developed framework provides a robust and efficient approach to performing a large number of user defined MC simulations with computational voxel phantoms in CT (minimal user interaction). The resulting organ doses from each simulation can be accessed through a database which greatly increases the ease of analyzing the resulting organ doses. The framework developed in this work provides a valuable resource when investigating different dose optimization strategies in CT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Automated platform to assess commercial off the shelf (COTS) software assurance

A computer-implemented method of verifying software is provided. The method comprises creating a number of virtual machines that simulate computing environments and running a number of software program on the virtual machines. The software programs have full access to the simulated computing environments, but the source code of the software program is unavailable. A hypervisor performs virtual machine introspection as the software programs run on the virtual machines, wherein the virtual machines and software programs are unaware the virtual machine introspection is being performed. Telemetry data is collected about the software programs, including any identified threats posed by the software programs to the simulated computing environments, and presented to a user via an interface.

Urias, Vincent↗
DOE OSTI2025

Investigation of Electrocatalytic CO 2 Reduction on MXene Materials via First-Principles Simulations

Computational studies of CO 2 reduction to yield various products were carried out on the basal plane and edges of different MXene materials. The impact of vacancies upon Mo 2 TiC 2 T x , W 2 TiC 2 T x , and Ti 3 C 2 T x (T x = O and OH) was also examined for both edge and basal sites. Initial calibrations were carried out to generate surfaces with optimal oxide/hydroxide ratios and proper termination sites upon which various vacancy sites were explored to ensure an accurate model of the surfaces’ resting states. From this work, Mo 2 TiC 2 O was determined to exhibit the lowest theoretical overpotential for methane as determined by volcano plot analyses. At a large enough vacancy concentration, the CO 2 reduction reaction (CO 2 RR) is predicted to outcompete the hydrogen evolution reaction (HER) as the predominant reaction on the surface. In conclusion, when examining the edge of Mo 2 TiC 2 O, stronger CO 2 binding was exhibited to such an extent that the reaction was predicted to terminate after the generation of formate on the edge.

Lund, Colton [Argonne National Laboratory (ANL), A↗
DOE OSTI2013

Estimating building simulation parameters via Bayesian structure learning

Many key building design policies are made using sophisticated computer simulations such as EnergyPlus (E+), the DOE flagship whole-building energy simulation engine. E+ and other sophisticated computer simulations have several major problems. The two main issues are 1) gaps between the simulation model and the actual structure, and 2) limitations of the modeling engine's capabilities. Currently, these problems are addressed by having an engineer manually calibrate simulation parameters to real world data or using algorithmic optimization methods to adjust the building parameters. However, some simulations engines, like E+, are computationally expensive, which makes repeatedly evaluating the simulation engine costly. This work explores addressing this issue by automatically discovering the simulation's internal input and output dependencies from ~20 Gigabytes of E+ simulation data, future extensions will use ~200 Terabytes of E+ simulation data. The model is validated by inferring building parameters for E+ simulations with ground truth building parameters. Our results indicate that the model accurately represents parameter means with some deviation from the means, but does not support inferring parameter values that exist on the distribution's tail.

Edwards, Richard E.↗
DOE OSTI2013

Estimating Building Simulation Parameters via Bayesian Structure Learning

Many key building design policies are made using sophisticated computer simulations such as EnergyPlus (E+), the DOE flagship whole-building energy simulation engine. E+ and other sophisticated computer simulations have several major problems. The two main issues are 1) gaps between the simulation model and the actual structure, and 2) limitations of the modeling engine's capabilities. Currently, these problems are addressed by having an engineer manually calibrate simulation parameters to real world data or using algorithmic optimization methods to adjust the building parameters. However, some simulations engines, like E+, are computationally expensive, which makes repeatedly evaluating the simulation engine costly. This work explores addressing this issue by automatically discovering the simulation's internal input and output dependencies from 20 Gigabytes of E+ simulation data, future extensions will use 200 Terabytes of E+ simulation data. The model is validated by inferring building parameters for E+ simulations with ground truth building parameters. Our results indicate that the model accurately represents parameter means with some deviation from the means, but does not support inferring parameter values that exist on the distribution's tail.

Energy Modeling↗
DOE OSTI2013

Development of a fourth generation predictive capability maturity model.

The Predictive Capability Maturity Model (PCMM) is an expert elicitation tool designed to characterize and communicate completeness of the approaches used for computational model definition, verification, validation, and uncertainty quantification associated for an intended application. The primary application of this tool at Sandia National Laboratories (SNL) has been for physics-based computational simulations in support of nuclear weapons applications. The two main goals of a PCMM evaluation are 1) the communication of computational simulation capability, accurately and transparently, and 2) the development of input for effective planning. As a result of the increasing importance of computational simulation to SNL s mission, the PCMM has evolved through multiple generations with the goal to provide more clarity, rigor, and completeness in its application. This report describes the approach used to develop the fourth generation of the PCMM.

Hills, Richard Guy↗
DOE OSTI2010

Dosimetric characterization of a dedicated breast computed tomography clinical prototype

Purpose: To investigate the glandular dose magnitudes and characteristics resulting from image acquisition using a dedicated breast computed tomography (BCT) clinical prototype imaging system. Methods: The x-ray spectrum and output characteristics of a BCT clinical prototype (Koning Corporation, West Henrietta, NY) were determined using empirical measurements, breast phantoms, and an established spectrum model. The geometry of the BCT system was replicated in a Monte Carlo-based computer simulation using the GEANT4 toolkit and was validated by comparing the simulated results for exposure distribution in a standard 16 cm CT head phantom with those empirically determined using a 10 cm CT pencil ionization chamber and dosimeter. The computer simulation was further validated by replicating the results of a previous BCT dosimetry study. Upon validation, the computer simulation was modified to include breasts of varying sizes and homogeneous compositions spanning those encountered clinically, and the normalized mean glandular dose resulting from BCT was determined. Using the system's measured exposure output determined automatically for breasts of different size and density, the mean glandular dose for these breasts was computed and compared to the glandular dose resulting from mammography. Finally, additional Monte Carlo simulations were performed to study how the glandular dose values vary within the breast tissue during acquisition with both this BCT prototype and a typical craniocaudal (CC) mammographic acquisition. Results: This BCT prototype uses an x-ray spectrum with a first half-value layer of 1.39 mm Al and a mean x-ray energy of 30.3 keV. The normalized mean glandular dose for breasts of varying size and composition during BCT acquisition with this system ranges from 0.278 to 0.582 mGy/mGy air kerma with the reference air kerma measured in air at the center of rotation. Using the measured exposure outputs for the tube currents automatically selected by the system for the breasts of different sizes and densities, the mean glandular dose for a BCT acquisition with this prototype system varies from 5.6 to 17.5 mGy, with the value for a breast of mean size and composition being 17.06 mGy. The glandular dose throughout the breast tissue of this mean breast varies by up to {+-}50% of the mean value. During a typical CC view mammographic acquisition of an equivalent mean breast, which typically results in a mean glandular dose of 2.0-2.5 mGy, the glandular dose throughout the breast tissue varies from {approx}15% to {approx}400% of the mean value. Conclusions: Acquisition of a BCT image with the automated tube output settings for a mean breast with the Koning Corp. clinical prototype results in mean glandular dose values approximately equivalent to three to five two-view mammographic examinations for a similar breast. For all breast sizes and compositions studied, this glandular dose ratio between acquisition with this BCT prototype and two-view mammography ranges from 1.4 to 7.2. In mammography, portions of the mean-sized breast receive a considerably higher dose than the mean value for the whole breast. However, only a small portion of a breast undergoing mammography would receive a glandular dose similar to that from BCT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Emerging Computing Architectures: Simulation of Power Electronics in Power Grids

As the penetration of power electronics increases in power grids, new computing architectures needed to be evaluated for the simulation of high-fidelity models of power electronics in power grids in operations. In this paper, emerging computing architectures such as quantum processing units are evaluated for the simulation of power electronics in power grids. A hybrid algorithm based on classical computing and quantum computing is developed and tested for different use cases of electromagnetic transient (EMT) simulation of power electronics (PE)-based systems and simple circuits. The algorithms needed to simulate power electronics in emerging computing architectures are discussed and thereafter, simulation results are shown.

Debnath, Suman↗
DOE OSTI2022

Development and assessment of a reactor system prognosis model with physics-guided machine learning

Autonomous control systems provide recommendations to help operators in decision-making during plant operations ranging from normal operation to accident management. An important step of autonomous control is prognosis. In nuclear engineering domain, prognosis is the process of predicting future conditions of a system or equipment based on present signs and symptoms of a fault. The prognosis model allows predicting future reactor states for possible candidate control strategies so that the outcomes can be evaluated to determine the best control strategy. The prognosis model requires representing direct relationships between the symptoms and the predictions. In nuclear engineering, computational simulations are approximate representations of the operation of the real system. However, prognosis with computational simulations requires high computation power and time due to possible large number of scenarios. Necessary computation resources can be reduced with machine learning (ML) approach for fast predictions by building a surrogate function using the simulation data. A critical issue is, ML models are ignorant of physical knowledge, and these models approximate statistical relationships between the system variables. This ignorance can produce results that are inconsistent with physical laws, even if an optimal result is achieved from a mathematical point of view. Physics-guided machine learning (PGML) is an approach to tackle this issue. Here, this work formulates and illustrates a framework to guide development and assessment of the ML-based prognosis model for autonomous control systems. The development of the prognosis model considers the training of a ML model which consists of optimizing many aspects of the ML approach. The assessment of the prognosis model considers training data limitations and uncertainties of the ML approach. Prognosis models with standalone ML and PGML are developed and assessed on the loss-of-flow scenario of Experimental Breeder Reactor II. The results indicate that PGML based prognosis model has the best performance compared to other prognosis models.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Frequency and Clinical Significance of Previously Undetected Incidental Findings Detected on Computed Tomography Simulation Scans for Breast Cancer Patients

Purpose: To determine the frequency and clinical significance of previously undetected incidental findings found on computed tomography (CT) simulation images for breast cancer patients. Methods and Materials: All CT simulation images were first interpreted prospectively by radiation oncologists and then double-checked by diagnostic radiologists. The official reports of CT simulation images for 881 consecutive postoperative breast cancer patients from 2009 to 2010 were retrospectively reviewed. Potentially important incidental findings (PIIFs) were defined as any previously undetected benign or malignancy-related findings requiring further medical follow-up or investigation. For all patients in whom a PIIF was detected, we reviewed the clinical records to determine the clinical significance of the PIIF. If the findings from the additional studies prompted by a PIIF required a change in management, the PIIF was also recorded as a clinically important incidental finding (CIIF). Results: There were a total of 57 (6%) PIIFs. The 57 patients in whom a PIIF was detected were followed for a median of 17 months (range, 3-26). Six cases of CIIFs (0.7% of total) were detected. Of the six CIIFs, three (50%) cases had not been noted by the radiation oncologist until the diagnostic radiologist detected the finding. On multivariate analysis, previous CT examination was an independent predictor for PIIF (p = 0.04). Patients who had not previously received chest CT examinations within 1 year had a statistically significantly higher risk of PIIF than those who had received CT examinations within 6 months (odds ratio, 3.54; 95% confidence interval, 1.32-9.50; p = 0.01). Conclusions: The rate of incidental findings prompting a change in management was low. However, radiation oncologists appear to have some difficulty in detecting incidental findings that require a change in management. Considering cost, it may be reasonable that routine interpretations are given to those who have not received previous chest CT examinations within 1 year.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A line-source method for aligning on-board and other pinhole SPECT systems

Purpose: In order to achieve functional and molecular imaging as patients are in position for radiation therapy, a robotic multipinhole SPECT system is being developed. Alignment of the SPECT system—to the linear accelerator (LINAC) coordinate frame and to the coordinate frames of other on-board imaging systems such as cone-beam CT (CBCT)—is essential for target localization and image reconstruction. An alignment method that utilizes line sources and one pinhole projection is proposed and investigated to achieve this goal. Potentially, this method could also be applied to the calibration of the other pinhole SPECT systems.Methods: An alignment model consisting of multiple alignment parameters was developed which maps line sources in three-dimensional (3D) space to their two-dimensional (2D) projections on the SPECT detector. In a computer-simulation study, 3D coordinates of line-sources were defined in a reference room coordinate frame, such as the LINAC coordinate frame. Corresponding 2D line-source projections were generated by computer simulation that included SPECT blurring and noise effects. The Radon transform was utilized to detect angles (α) and offsets (ρ) of the line-source projections. Alignment parameters were then estimated by a nonlinear least squares method, based on the α and ρ values and the alignment model. Alignment performance was evaluated as a function of number of line sources, Radon transform accuracy, finite line-source width, intrinsic camera resolution, Poisson noise, and acquisition geometry. Experimental evaluations were performed using a physical line-source phantom and a pinhole-collimated gamma camera attached to a robot.Results: In computer-simulation studies, when there was no error in determining angles (α) and offsets (ρ) of the measured projections, six alignment parameters (three translational and three rotational) were estimated perfectly using three line sources. When angles (α) and offsets (ρ) were provided by the Radon transform, estimation accuracy was reduced. The estimation error was associated with rounding errors of Radon transform, finite line-source width, Poisson noise, number of line sources, intrinsic camera resolution, and detector acquisition geometry. Statistically, the estimation accuracy was significantly improved by using four line sources rather than three and by thinner line-source projections (obtained by better intrinsic detector resolution). With five line sources, median errors were 0.2 mm for the detector translations, 0.7 mm for the detector radius of rotation, and less than 0.5° for detector rotation, tilt, and twist. In experimental evaluations, average errors relative to a different, independent registration technique were about 1.8 mm for detector translations, 1.1 mm for the detector radius of rotation (ROR), 0.5° and 0.4° for detector rotation and tilt, respectively, and 1.2° for detector twist.Conclusions: Alignment parameters can be estimated using one pinhole projection of line sources. Alignment errors are largely associated with limited accuracy of the Radon transform in determining angles (α) and offsets (ρ) of the line-source projections. This alignment method may be important for multipinhole SPECT, where relative pinhole alignment may vary during rotation. For pinhole and multipinhole SPECT imaging on-board radiation therapy machines, the method could provide alignment of SPECT coordinates with those of CBCT and the LINAC.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Thermal Neutron Scattering Evaluation Methods Development at ORNL [Slides]

The goal of the project was to develop a procedure to generate covariance matrix for S( α,β ) data that incorporated both computational simulations and experimental data. Experimental data and computer simulation fit is achieved using the Unified Monte Carlo (UMC) method. The framework is material & simulation method independent, and this is demonstrated using light water. A methodology was developed for generalized generation of covariances for thermal scattering files. The framework was subsequently tested on light water by generating covariance of pDOS with encouraging results. Preliminary results suggest temperature effects of polyethylene may warrant further investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Stochastic modeling of catalytic processes in nanoporous materials: Beyond mean-field approach

Transport and reaction in zeolites and other porous materials, such as mesoporous silica particles, has been a focus of interest in recent years. This is in part due to the possibility of anomalous transport effects (e.g. single-file diffusion) and its impact in the reaction yield in catalytic processes. Computational simulations are often used to study these complex nonequilibrium systems. Computer simulations using Molecular Dynamics (MD) techniques are prohibitive, so instead coarse grained one-dimensional models with the aid of Kinetic Monte Carlo (KMC) simulations are used. Both techniques can be computationally expensive, both time and resource wise. These coarse-grained systems can be exactly described by a set of coupled stochastic master equations, that describe the reaction-diffusion kinetics of the system. The equations can be written exactly, however, coupling between the equations and terms within the equations make it impossible to solve them exactly; approximations must be made. One of the most common methods to obtain approximate solutions is to use Mean Field (MF) theory. MF treatments yield reasonable results at high ratios of reaction rate k to hop rate h of the particles, but fail completely at low k=h due to the over-estimation of fluxes of particles within the pore. We develop a method to estimate fluxes and intrapore diffusivity in simple one- dimensional reaction-diffusion models at high and low k=h, where the pores are coupled to an equilibrated three-dimensional fluid. We thus successfully describe analytically these simple reaction-diffusion one-dimensional systems. Extensions to models considering behavior with long range steric interactions and wider pores require determination of multiple boundary conditions. We give a prescription to estimate the required parameters for these simulations. For one dimensional systems, if single-file diffusion is relaxed, additional parameters to describe particle exchange have to be introduced. We use Langevin Molecular Dynamics (MD) simulations to assess these parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Pathway Controlled Penetration (PcP) [Slides]

The technical approach employs advanced computational simulation tools to demonstrate how current assets can destroy RWK-RFI-12-0001's HDBT, a tunnel complex with two portals built into the base of a granite mountain. The granite over layer is assumed to be 60 meters thick over both portals and 80 meters over the facility's mission space. Key S&T is the completed development of a highly innovative viscoplastic fracture material model, 3D parallel gas-fracture capabilities into FDEM, and a stochastic handling of the material properties. Phase I - Develop and validate code simulation tools: (1) develop, incorporate and validate AZ-Frac material model for granite; and (2) Develop and incorporate gas-driven-fracture modeling into LANL's FDEM MUNROU code; (3) Develop and incorporate stochastic features into FDEM modeling. Phase II - Conduct PcP analysis on above HDBT: (1) Acquire HDBT design data, develop simulation model; and (2) Evaluate and select most promising defeat alternative. Phase III - Deliver code, train Service target analysts, and conduct simulations against real world HDBTs. PcP uses advanced computer simulations to enhance HDBT functional defeat efforts. Newly developed material models that account for fractural energy coupled with the finite discrete element methodology (FDEM) will provide targeting packages that will create penetration avenues for current or future lethality options. This novel computational approach requires full 3D geologic and structure characterization as well as significant high performance computing capabilities. The goal is to distinctively alter the targeting paradigm by leveraging critical DoD assets along with insitu geologic strata. In other words, assets will utilize underground rock structure to their benefit by creating rubbilization zones that will allow pathway controlled penetration.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Applying Quantum Computing to Simulate Power System Dynamics

Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. This paper demonstrates the potential use of quantum computing algorithms to model the power system dynamics. Leveraging a symbolic programming framework, we equivalently convert the power system dynamics’ differential algebraic equations (DAEs) into ordinary differential equations (ODEs), where the data of the state vector can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion, the quantum state tensor, and Hamiltonian simulation, whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can simulate the dynamics of the power system with high accuracy, whereas its complexity is polynomial in the logarithm of the system dimension. Our work also illustrates the use of scientific machine learning tools for implementing scientific computing concepts, e.g., Taylor expansion, DAEs/ODEs transform, and quantum computing solver, in the field of power engineering.

Tran, Huynh↗