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

Online Convex Optimization of Programmable Quantum Computers to Simulate Time-Varying Quantum Channels

Simulating quantum channels is a fundamental primitive in quantum computing, since quantum channels define general (trace-preserving) quantum operations. An arbitrary quantum channel cannot be exactly simulated using a finite-dimensional programmable quantum processor, making it important to develop optimal approximate simulation techniques. In this paper, we study the challenging setting in which the channel to be simulated varies adversarially with time. We propose the use of matrix exponentiated gradient descent (MEGD), an online convex optimization method, and analytically show that it achieves a sublinear regret in time. Through experiments, we validate the main results for time-varying dephasing channels using a programmable generalized teleportation processor.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Free Energy and Heat Capacity

Free energy and heat capacity of actinide elements and compounds are important properties for the evaluation of the safety and reliable performance of nuclear fuel. They are essential inputs for models that describe complex phenomena that govern the behaviour of actinide compounds during nuclear fuel fabrication and irradiation. This chapter introduces various experimental methods to measure free energy and heat capacity to serve as inputs for models and to validate computer simulations. This is followed by a discussion of computer simulation of these properties, and recent simulations of thermophysical properties of nuclear fuel are briefly reviewed.

free energy↗
DOE OSTI2012

The second-order differential phase contrast and its retrieval for imaging with x-ray Talbot interferometry

Purpose: The x-ray differential phase contrast imaging implemented with the Talbot interferometry has recently been reported to be capable of providing tomographic images corresponding to attenuation-contrast, phase-contrast, and dark-field contrast, simultaneously, from a single set of projection data. The authors believe that, along with small-angle x-ray scattering, the second-order phase derivative {Phi}{sup Double-Prime }{sub s}(x) plays a role in the generation of dark-field contrast. In this paper, the authors derive the analytic formulae to characterize the contribution made by the second-order phase derivative to the dark-field contrast (namely, second-order differential phase contrast) and validate them via computer simulation study. By proposing a practical retrieval method, the authors investigate the potential of second-order differential phase contrast imaging for extensive applications. Methods: The theoretical derivation starts at assuming that the refractive index decrement of an object can be decomposed into {delta}={delta}{sub s}+{delta}{sub f}, where {delta}{sub f} corresponds to the object's fine structures and manifests itself in the dark-field contrast via small-angle scattering. Based on the paraxial Fresnel-Kirchhoff theory, the analytic formulae to characterize the contribution made by {delta}{sub s}, which corresponds to the object's smooth structures, to the dark-field contrast are derived. Through computer simulation with specially designed numerical phantoms, an x-ray differential phase contrast imaging system implemented with the Talbot interferometry is utilized to evaluate and validate the derived formulae. The same imaging system is also utilized to evaluate and verify the capability of the proposed method to retrieve the second-order differential phase contrast for imaging, as well as its robustness over the dimension of detector cell and the number of steps in grating shifting. Results: Both analytic formulae and computer simulations show that, in addition to small-angle scattering, the contrast generated by the second-order derivative is magnified substantially by the ratio of detector cell dimension over grating period, which plays a significant role in dark-field imaging implemented with the Talbot interferometry. Conclusions: The analytic formulae derived in this work to characterize the second-order differential phase contrast in the dark-field imaging implemented with the Talbot interferometry are of significance, which may initiate more activities in the research and development of x-ray differential phase contrast imaging for extensive preclinical and eventually clinical applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Case Study on Pathogen Transport, Deposition, Evaporation and Transmission: Linking High-Fidelity Computational Fluid Dynamics Simulations to Probability of Infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose–response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations for solvent-based carbon capture. Part 2: Chemical absorption across a wetted wall column: Original Research Article: Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations

Part 1 of this paper presents a numerical model for non-reactive physical mass transfer across a wetted wall column (WWC). In Part 2, we improved the existing computational fluid dynamics (CFD) model to simulate chemical absorption occurring in a WWC as a bench-scale study of solvent-based carbon dioxide (CO 2 ) capture. In this study, to generate data for WWC model validation, CO 2 mass transfer across a monoethanolamine (MEA) solvent was first measured on a WWC experimental apparatus. The numerical model developed in this work can account for both chemical absorption and desorption of CO 2 in MEA. In addition, the overall mass transfer coefficient predicted using traditional/empirical correlations is conducted and compared with CFD prediction results for both steady and wavy falling films. A Bayesian statistical calibration algorithm is adopted to calibrate the reaction rate constants in chemical absorption/desorption of CO 2 across a falling film of MEA. The posterior distributions of the two transport properties, i.e., Henry's constant and gas diffusivity in the non-reacting nitrous oxide (N 2 O)/MEA system obtained from Part 1 of this study, serves as priors for the calibration of CO 2 reaction rate constants after using the N 2 O/CO 2 analogy method. Finally, the calibrated model can be used to predict the CO 2 mass transfer in a WWC for a wider range of operating conditions.

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

pnnl/socialsim

The goal of Computational Simulation of Online Social Behavior (SocialSim) is to develop innovative technologies for highfidelity computational simulation of online social behavior. SocialSim will focus specifically on information spread and evolution. Current computational approaches to social and behavioral simulation are limited in this regard. Top-down simulation approaches focus on the dynamics of a population as a whole, and model behavioral phenomena by assuming uniform or mostly-uniform behavior across that population. Such methods can easily scale to simulate massive populations, but can be inaccurate if there are specific, distinct variations in the characteristics of the population. In contrast, bottom-up simulation approaches treat population dynamics as an emergent property of the activities and interactions taking place within a diverse population. While such approaches can enable more accurate simulation of information spread, they do not readily scale to represent large populations. SocialSim aims to develop novel approaches to address these challenges.

Smith, William↗
DOE OSTI2013

The putative liquid-liquid transition is a liquid-solid transition in atomistic models of water. II

This paper extends our earlier studies of free energy functions of density and crystalline order parameters for models of supercooled water, which allows us to examine the possibility of two distinct metastable liquid phases [D. T. Limmer and D. Chandler, J. Chem. Phys.135, 134503 (2011) and preprint http://arxiv.org/abs/arXiv:1107.0337 (2011)]. Low-temperature reversible free energy surfaces of several different atomistic models are computed: mW water, TIP4P/2005 water, Stillinger-Weber silicon, and ST2 water, the last of these comparing three different treatments of long-ranged forces. In each case, we show that there is one stable or metastable liquid phase, and there is an ice-like crystal phase. The time scales for crystallization in these systems far exceed those of structural relaxation in the supercooled metastable liquid. We show how this wide separation in time scales produces an illusion of a low-temperature liquid-liquid transition. The phenomenon suggesting metastability of two distinct liquid phases is actually coarsening of the ordered ice-like phase, which we elucidate using both analytical theory and computer simulation. For the latter, we describe robust methods for computing reversible free energy surfaces, and we consider effects of electrostatic boundary conditions. We show that sensible alterations of models and boundary conditions produce no qualitative changes in low-temperature phase behaviors of these systems, only marginal changes in equations of state. On the other hand, we show that altering sampling time scales can produce large and qualitative non-equilibrium effects. Recent reports of evidence of a liquid-liquid critical point in computer simulations of supercooled water are considered in this light.

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

FEMA Asteroid Impact Tabletop Exercise Simulations

We describe the computational simulations and damage assessments that we provided in support of a tabletop exercise (TTX) at the request of NASA's Near-Earth Objects Program Office. The overall purpose of the exercise was to assess leadership reactions, information requirements, and emergency management responses to a hypothetical asteroid impact with Earth. The scripted exercise consisted of discovery, tracking, and characterization of a hypothetical asteroid; inclusive of mission planning, mitigation, response, impact to population, infrastructure and GDP, and explicit quantification of uncertainty. Participants at the meeting included representatives of NASA, Department of Defense, Department of State, Department of Homeland Security/Federal Emergency Management Agency (FEMA), and the White House. The exercise took place at FEMA headquarters. Sandia's role was to assist the Jet Propulsion Laboratory (JPL) in developing the impact scenario, to predict the physical effects of the impact, and to forecast the infrastructure and economic losses. We ran simulations using Sandia's CTH hydrocode to estimate physical effects on the ground, and to produce contour maps indicating damage assessments that could be used as input for the infrastructure and economic models. We used the FASTMap tool to provide estimates of infrastructure damage over the affected area, and the REAcct tool to estimate the potential economic severity expressed as changes to GDP (by nation, region, or sector) due to damage and short-term business interruptions.

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

Visualization of Dynamic Vortex Structures in Magnetic Films with Uniaxial Anisotropy (Micromagnetic Simulation)

Three-dimensional computer simulation of dynamic processes in a moving domain boundary separating domains in a soft magnetic uniaxial film with planar anisotropy is performed by numerical solution of Landau-Lifshitz-Gilbert equations. The developed visualization methods are used to establish the connection between the motion of surface vortices and antivortices, singular (Bloch) points, and core lines of intrafilm vortex structures. A relation between the character of magnetization dynamics and the film thickness is found. The analytical models of spatial vortex structures for imitation of topological properties of the structures observed in micromagnetic simulation are constructed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Lie-algebraic classical simulations for quantum computing

The classical simulation of quantum dynamics plays an important role in our understanding of quantum complexity and in the development of quantum technologies. Efficient techniques such as those based on the Gottesman-Knill theorem for Clifford circuits, tensor networks for low entanglement-generating circuits, or Wick's theorem for fermionic Gaussian states have become central tools in quantum computing. In this work, we contribute to this body of knowledge by presenting a framework for classical simulations, dubbed “𝔤-sim”, which is based on the underlying Lie algebraic structure of the dynamical process. When the dimension of the algebra grows at most polynomially in the system size, there exist observables for which the simulation is efficient. Indeed, we show that 𝔤-sim enables new regimes for classical simulations, is able to deal with certain forms of noise in the evolution, as well as can be used to tackle several paradigmatic variational and nonvariational quantum computing tasks. For the former, we perform Lie-algebraic simulations to train and optimize parametrized quantum circuits (thus effectively showing that some variational models can be dequantized), design enhanced parameter initialization strategies, solve tasks of quantum circuit synthesis, and train a quantum-phase classifier. For the latter, we report large-scale noiseless and noisy simulations on benchmark problems. By comparing the limitations of 𝔤-sim and certain Wick's theorem-based simulations, we find that the two methods become inefficient for different types of states or observables, hinting at the existence of distinct, nonequivalent resources for classical simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

SimNet: Accurate and High-Performance Computer Architecture Simulation using Deep Learning

While cycle-accurate simulators are essential tools for architecture research, design, and development, their practicality is limited by an extremely long time-to-solution for realistic applications under investigation. This work describes a concerted effort, where machine learning (ML) is used to accelerate microarchitecture simulation. First, an ML-based instruction latency prediction framework that accounts for both static instruction properties and dynamic processor states is constructed. Then, a GPU-accelerated parallel simulator is implemented based on the proposed instruction latency predictor, and its simulation accuracy and throughput are validated and evaluated against a state-of-the-art simulator. Leveraging modern GPUs, the ML-based simulator outperforms traditional CPU-based simulators significantly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Relative plan robustness of step-and-shoot vs rotational intensity–modulated radiotherapy on repeat computed tomographic simulation for weight loss in head and neck cancer

Introduction: Interfractional anatomical alterations may have a differential effect on the dose delivered by step-and-shoot intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc therapy (VMAT). The increased degrees of freedom afforded by rotational delivery may increase plan robustness (measured by change in target volume coverage and doses to organs at risk [OARs]). However, this has not been evaluated for head and neck cancer. Materials and methods: A total of 10 patients who required repeat computed tomography (CT) simulation and replanning during head and neck IMRT were included. Step-and-shoot IMRT and VMAT plans were generated from the original planning scan. The initial and second CT simulation scans were fused and targets/OAR contours transferred, reviewed, and modified. The plans were applied to the second CT scan and doses recalculated without repeat optimization. Differences between step-and-shoot IMRT and VMAT for change in target volume coverage and doses to OARs between first and second CT scans were compared by Wilcoxon signed rank test. Results: There were clinically relevant dosimetric changes between the first and the second CT scans for both the techniques (reduction in mean D{sub 95%} for PTV2 and PTV3, D{sub min} for CTV2 and CTV3, and increased mean doses to the parotid glands). However, there were no significant differences between step-and-shoot IMRT and VMAT for change in any target coverage parameter (including D{sub 95%} for PTV2 and PTV3 and D{sub min} for CTV2 and CTV3) or dose to any OARs (including parotid glands) between the first and the second CT scans. Conclusions: For patients with head and neck cancer who required replanning mainly due to weight loss, there were no significant differences in plan robustness between step-and-shoot IMRT and VMAT. This information is useful with increased clinical adoption of VMAT.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Cold ion beam in a storage ring as a platform for large-scale quantum computers and simulators: Challenges and directions for research and development

The purpose of this paper is to evaluate the possibility of constructing a large-scale storage-ring-type ion-trap system capable of storing, cooling, and controlling a large number of ions as a platform for scalable quantum computing (QC) and quantum simulations. In such a trap, the ions form a crystalline beam moving along a circular path with a constant velocity determined by the frequency and intensity of the cooling lasers. In this paper, we consider a large leap forward in terms of the number of ions that serve as qubits in QC, from fewer than 100 available in state of the art linear ion-trap devices today to an order of 10 5 crystallized ions in the storage-ring setup. This new trap design unifies two different concepts: the storage rings of charged particles and the linear ion traps used for QC and mass spectrometry. In this paper, we use the language of particle accelerators to discuss the ion state and dynamics. We outline the differences between the above concepts, analyze challenges of the large ring with a revolving chain of ions, and propose goals for the research and development required to enable future quantum computers with 1000 times more qubits than available today. The challenge of creating such a large-scale quantum system while maintaining the necessary coherence of the qubits and the high fidelity of quantum logic operations is significant. Performing analog quantum simulations may be an achievable initial goal for such a device. Quantum calculations and simulations of complex quantum systems will move forward both the fundamental science and the applied research. Nuclear and particle physics, many-body quantum systems, lattice gauge theories, and nuclear structure calculations are just a few examples in which a large-scale quantum simulation system will become a very powerful tool to move forward our understanding of nature.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Portal dosimetry for VMAT using integrated images obtained during treatment

Purpose: Portal dosimetry provides an accurate and convenient means of verifying dose delivered to the patient. A simple method for carrying out portal dosimetry for volumetric modulated arc therapy (VMAT) is described, together with phantom measurements demonstrating the validity of the approach. Methods: Portal images were predicted by projecting dose in the isocentric plane through to the portal image plane, with exponential attenuation and convolution with a double-Gaussian scatter function. Appropriate parameters for the projection were selected by fitting the calculation model to portal images measured on an iViewGT portal imager (Elekta AB, Stockholm, Sweden) for a variety of phantom thicknesses and field sizes. This model was then used to predict the portal image resulting from each control point of a VMAT arc. Finally, all these control point images were summed to predict the overall integrated portal image for the whole arc. The calculated and measured integrated portal images were compared for three lung and three esophagus plans delivered to a thorax phantom, and three prostate plans delivered to a homogeneous phantom, using a gamma index for 3% and 3 mm. A 0.6 cm{sup 3} ionization chamber was used to verify the planned isocentric dose. The sensitivity of this method to errors in monitor units, field shaping, gantry angle, and phantom position was also evaluated by means of computer simulations. Results: The calculation model for portal dose prediction was able to accurately compute the portal images due to simple square fields delivered to solid water phantoms. The integrated images of VMAT treatments delivered to phantoms were also correctly predicted by the method. The proportion of the images with a gamma index of less than unity was 93.7% ± 3.0% (1SD) and the difference between isocenter dose calculated by the planning system and measured by the ionization chamber was 0.8% ± 1.0%. The method was highly sensitive to errors in monitor units and field shape, but less sensitive to errors in gantry angle or phantom position. Conclusions: This method of predicting integrated portal images provides a convenient means of verifying dose delivered using VMAT, with minimal image acquisition and data processing requirements.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Surrogate Model Guided Optimization of Expensive Black-Box Multi-Objective Problems: A Posteriori Methods

Many engineering applications require the simultaneous optimization of multiple conflicting objective functions. Often, these objective functions are evaluated using highly accurate computer simulations that are computationally too expensive to be evaluated hundreds or thousands of times during optimization. Thus, the goal is to find good approximations of the Pareto front using as few of these expensive simulations as possible. Here, we describe an optimization approach based on surrogate models and diverse sampling strategies to accelerate the search for the Pareto solutions. We use a separate surrogate model for approximating each objective function and then we use the surrogate models to inform where additional expensive simulations should be run. The surrogate models are updated in an active learning framework whenever new information from the expensive simulations becomes available. The sampling strategies aim at balancing local improvements of the approximate Pareto front and global exploration to identify the extrema and fill in large gaps of the approximate Pareto front. We demonstrate on a large set of benchmark problems the effectiveness of the method for finding good approximations of the Pareto front.

MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Simulation benchmarks for low-pressure plasmas: Capacitive discharges

Benchmarking is generally accepted as an important element in demonstrating the correctness of computer simulations. In the modern sense, a benchmark is a computer simulation result that has evidence of correctness, is accompanied by estimates of relevant errors, and which can thus be used as a basis for judging the accuracy and efficiency of other codes. In this paper, we present four benchmark cases related to capacitively coupled discharges. These benchmarks prescribe all relevant physical and numerical parameters. We have simulated the benchmark conditions using five independently developed particle-in-cell codes. We show that the results of these simulations are statistically indistinguishable, within bounds of uncertainty that we define. We, therefore, claim that the results of these simulations represent strong benchmarks, which can be used as a basis for evaluating the accuracy of other codes. These other codes could include other approaches than particle-in-cell simulations, where benchmarking could examine not just implementation accuracy and efficiency, but also the fidelity of different physical models, such as moment or hybrid models. We discuss an example of this kind in the Appendix. Of course, the methodology that we have developed can also be readily extended to a suite of benchmarks with coverage of a wider range of physical and chemical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Temperature trends and correlation between SQUID superparamagnetic relaxometry and dc-magnetization on model iron-oxide nanoparticles

Structural properties of iron-oxide nanoparticles deeply affect their magnetic performance in many applications such as with superparamagnetic relaxometry, when employed as cell-targeted magnetic nanoparticles for in vivo diagnostics. In this work, we present a detailed characterization of model nanoparticles for this application, with an average size of ~25 nm and a narrow size dispersion (σ < 7%). Considering the intrinsic structural properties of these model nanoparticles, the study of temperature dependence and correlation between dc-magnetization and superconducting quantum interference detector-relaxometry are discussed based on known theoretical predictions and computer simulations of the magnetic dipole moment and characteristic decay constants. Furthermore, computer simulations provide support in clarifying how important the overall collective magnetization is affected by particle size dispersion, which has a direct role on sustaining the magnetic relaxation signal in the temperature range required in preclinical and clinical settings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗