Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

3D Simulations of Ion Thruster Accelerator Grid Erosion Accounting for Charge Exchange Ion Space Charge

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

Chaplin, Vernon H.↗
DOE OSTI2020

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

High Fidelity, “Faster than Real-Time” Simulator for Predicting Power System Dynamic Behavior - Final Technical Report

The “High Fidelity, Faster than Real­Time Simulator for Predicting Power System Dynamic Behavior” was designed and developed by Illinois Institute of Technology with critical contributions from Electrocon International, Argonne National Laboratory, Alstom Grid and McCoy Energy. Also essential to the project were our two utility partners: Commonwealth Edison and AltaLink. The project was a success due to several major breakthroughs in the area of large­scale power system dynamics simulation, including (1) a validated faster than real­ time simulation of both stable and unstable transient dynamics in a large­scale positive sequence transmission grid model, (2) a three­phase unbalanced simulation platform for modeling new grid devices, such as independently controlled single­phase static var compensators (SVCs), (3) the world’s first high fidelity three­phase unbalanced dynamics and protection simulator based on Electrocon’s CAPE program, and (4) a first­of­its­ kind implementation of a single­phase induction motor model with stall capability. The simulator results will aid power grid operators in their true time of need, when there is a significant risk of cascading outages. The simulator will accelerate performance and enhance accuracy of dynamics simulations, enabling operators to maintain reliability and steer clear of blackouts. In the long­term, the simulator will form the backbone of the newly conceived hybrid real­time protection and control architecture that will coordinate local controls, wide­area measurements, wide­area controls and advanced real­time prediction capabilities. The nation’s citizens will benefit in several ways, including (1) less down time from power outages due to the faster­than­real­time simulator’s predictive capability, (2) higher levels of reliability due to the detailed dynamics plus protection simulation capability, and (3) more resiliency due to the three­ phase unbalanced simulator’s ability to model three­phase and single­ phase networks and devices.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1972

Apollo oxygen tank stratification analysis, volume 2

An analysis of flight performance of the Apollo 15 cryogenic oxygen tanks was conducted with the variable grid stratification math model developed earlier in the program. Flight conditions investigated were the CMP-EVA and one passive thermal control period which exhibited heater temperature characteristics not previously observed. Heater temperatures for these periods were simulated with the math model using flight acceleration data. Simulation results (heater temperature and tank pressure) compared favorably with the Apollo 15 flight data, and it was concluded that tank performance was nominal. Math model modifications were also made to improve the simulation accuracy. The modifications included the addition of the effects of the tank wall thermal mass and an improved system flow distribution model. The modifications improved the accuracy of simulated pressure response based on comparisons with flight data.

Barton, J. E.↗
NASA NTRS1991

RMS active damping augmentation

The topics are presented in viewgraph form and include: RMS active damping augmentation; potential space station assembly benefits to CSI; LaRC/JSC bridge program; control law design process; draper RMS simulator; MIMO acceleration control laws improve damping; potential load reduction benefit; DRS modified to model distributed accelerations; accelerometer location; Space Shuttle aft cockpit simulator; simulated shuttle video displays; SES test goals and objectives; and SES modifications to support RMS active damping augmentation.

Michael G Gilbert↗
DOE OSTI2015

Building blocks for correlated superconductors and magnets

Recent efforts at Los Alamos to discover strongly correlated superconductors and hard ferromagnets are reviewed. While serendipity remains a principal engine of materials discovery, design principles and structural building blocks are beginning to emerge that hold potential for predictive discovery. Successes over the last decade with the so-called “115” strongly correlated superconductors are summarized, and more recent efforts to translate these insights and principles to novel hard magnets are discussed. While true “materials by design” remains a distant aspiration, progress is being made in coupling empirical design principles to electronic structure simulation to accelerate and guide materials design and synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Quantifying the accuracy of the tumor motion and area as a function of acceleration factor for the simulation of the dynamic keyhole magnetic resonance imaging method

Purpose: The dynamic keyhole is a new MR image reconstruction method for thoracic and abdominal MR imaging. To date, this method has not been investigated with cancer patient magnetic resonance imaging (MRI) data. The goal of this study was to assess the dynamic keyhole method for the task of lung tumor localization using cine-MR images reconstructed in the presence of respiratory motion. Methods: The dynamic keyhole method utilizes a previously acquired a library of peripheral k-space datasets at similar displacement and phase (where phase is simply used to determine whether the breathing is inhale to exhale or exhale to inhale) respiratory bins in conjunction with central k-space datasets (keyhole) acquired. External respiratory signals drive the process of sorting, matching, and combining the two k-space streams for each respiratory bin, thereby achieving faster image acquisition without substantial motion artifacts. This study was the first that investigates the impact of k-space undersampling on lung tumor motion and area assessment across clinically available techniques (zero-filling and conventional keyhole). In this study, the dynamic keyhole, conventional keyhole and zero-filling methods were compared to full k-space dataset acquisition by quantifying (1) the keyhole size required for central k-space datasets for constant image quality across sixty four cine-MRI datasets from nine lung cancer patients, (2) the intensity difference between the original and reconstructed images in a constant keyhole size, and (3) the accuracy of tumor motion and area directly measured by tumor autocontouring. Results: For constant image quality, the dynamic keyhole method, conventional keyhole, and zero-filling methods required 22%, 34%, and 49% of the keyhole size (P < 0.0001), respectively, compared to the full k-space image acquisition method. Compared to the conventional keyhole and zero-filling reconstructed images with the keyhole size utilized in the dynamic keyhole method, an average intensity difference of the dynamic keyhole reconstructed images (P < 0.0001) was minimal, and resulted in the accuracy of tumor motion within 99.6% (P < 0.0001) and the accuracy of tumor area within 98.0% (P < 0.0001) for lung tumor monitoring applications. Conclusions: This study demonstrates that the dynamic keyhole method is a promising technique for clinical applications such as image-guided radiation therapy requiring the MR monitoring of thoracic tumors. Based on the results from this study, the dynamic keyhole method could increase the imaging frequency by up to a factor of five compared with full k-space methods for real-time lung tumor MRI.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Accelerating dissipative particle dynamics simulations for soft matter systems

Dissipative particle dynamics (DPD) is a coarse-grained particle-based simulation method that offers microscopic-scale insights into soft matter systems. We present an efficient implementation of a DPD model for graphical processing units (GPUs). As implemented in the LAMMPS molecular dynamics package, it can run effectively on current-generation supercomputers which often have hybrid nodes containing multi-core CPUs and (one or more) GPUs. Using efficient communication of information between the CPUs and GPUs, DPD interactions can be computed on the GPU while other portions of a full simulation model (boundary conditions, constraints, bonded interactions, diagnostic calculations, etc.) can be performed on the CPU. Our GPU-enhanced runs show a speedup of up to 9.5× versus many-core CPU simulations, and can run scalably across thousands of compute nodes. We briefly discuss how the new GPU implementation was validated against the CPU version for thermodynamics, diffusion, and hydrodynamic behavior. We also highlight large-scale models which the faster DPD implementation has enabled, for studies of monolayer self-assembly and thin-film instabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2018

WIPP PA Sensitivity Studies in Support of CRA-2014 Completeness Determination and Compliance Decision - 18209

The Waste Isolation Pilot Plant (WIPP) has been developed by the U.S. Department of Energy (DOE) for the geologic (deep underground) disposal of transuranic (TRU) waste. Containment of TRU waste at the WIPP is regulated by the U.S. Environmental Protection Agency (EPA). The DOE demonstrates compliance with the containment requirements by means of performance assessment (PA) calculations. WIPP PA calculations estimate the probability and consequence of potential radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. Since TRU waste was first received at the WIPP in 1999, compliance recertification applications (CRAs) have been required by the EPA on five-year cycle. The latest CRA (CRA-2014) was submitted by the DOE to the EPA in 2014. As part of the application, a formal set of PA calculations (CRA-2014 PA) was performed and results were documented in detail. The CRA-2014 PA incorporated several parameter updates based on new data, as well as the implementation of a different model for panel closures, barriers planned for the opening of each waste panel. The results presented as part of the CRA-2014 showed via the required complementary cumulative distribution function (CCDF) release curves that the estimated releases were within compliance limits. As part of their process for reviewing the completeness of the CRA-2014 [1], the EPA formally requested that the DOE perform sensitivity studies based on EPA-requested changes to model parameters. The DOE performed and documented the results of three sensitivity studies described below: CRA14-SEN2 [2], CRA14-SEN3 [3], and CRA14-SEN4 [4]. In each case, the new analysis results were compared to CRA- 2014 PA (CRA14) results to determine the impact of the EPA-requested changes. CRA14-SEN2: the application of EPA-requested modified parameters in the operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model resulted in increased pressures and decreased brine saturations in waste areas and increased pressures and brine saturations in the operations and experimental areas. However, the overall effects on total high-probability mean releases from the repository were entirely insignificant. It was concluded that the modeling assumptions associated with the operations and experimental areas of the repository have an insignificant effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN3: the application of EPA-requested modified parameters in the panel closure areas and operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model predominantly result in increased pressures and decreased brine saturations in waste areas, increased pressures and brine saturations in the operations and experimental areas, and decreased brine and gas flows across all panel closures. The pressure increases in repository waste regions yielded slightly decreased brine saturations (on average) in those areas. However, the overall effects on total high-probability mean releases from the repository are entirely insignificant. It was concluded that the modeling assumptions made for this sensitivity study which were associated with the panel closures and north end areas of the repository have only a small effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN4: the application of several changes to the current WIPP PA model (CRA14), including: parameter changes, use of an updated code, and a revised computational grid for the BRAGFLO code. The modifications to the repository model resulted in increased releases in all primary release mechanisms. The impacts of each EPA-requested change to CRA14 were analyzed with respect to each release mechanism in the CRA14-SEN4 study. Overall, total high-probability predicted mean releases from the repository were increased by about 15%. Total low-probability predicted mean releases were increased by about 107%. It was concluded that the EPA-requested changes to the CRA14 result in increases to the predicted total releases from the repository. However, releases calculated in the CRA14-SEN4 analysis remained below regulatory limits, demonstrating continued compliance of the WIPP. In January 2017, the EPA determined that CRA-2014 was complete and then began a technical evaluation of CRA-2014 to determine whether the WIPP continued to comply with the radiation protection standards for disposal. During their evaluation, the EPA utilized the sensitivity study results to inform their ultimate determination in July 2017 that the WIPP continues to comply with regulatory standards. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

A plasma–vacuum interface tracking algorithm for magnetohydrodynamic simulations of coaxial plasma accelerators

The resistive Magneto-hydrodynamic (MHD) model describes the behavior of a strongly ionized plasma in the presence of external electric and magnetic fields. For problems involving the plasma with an open boundary to very low pressure/vacuum regions, the continuum assumption is no longer valid in the entire domain of interest. For example, this is seen in a plasma accelerator where a dense plasma expands into a vacuum background. A common practice to deal with this issue is to assign small values of density and pressure to the vacuum regions and proceed to solve the continuum based MHD equations throughout the domain. We show that this approach fails to produce solutions consistent with the physics of the problem and can give rise to unacceptable artifacts such as spurious shocks. Here, we develop a plasma–vacuum interface tracking approach to mitigate this problem. The plasma–vacuum interface is tracked in a physically consistent manner and the MHD equations are solved only in the regions that contain the plasma. The interface tracking is achieved using a face-flux formulation derived from the theoretical solution to a 1D free expansion problem. Coupled with a threshold based approach, the interface tracking is implemented for both explicit and implicit time stepping frameworks on generalized unstructured grids. Simulations of magnetized thermal plasma jets expanding into a vacuum background indicate plume profiles devoid of the unphysical shock obtained using the small background density approach. In the context of resistive MHD simulations, the interface tracking approach overcomes the numerical stiffness induced by specifying a large background resistivity in the vacuum regions, resulting in significant wall-clock time gains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Accelerated Snowmelt Protocol to Simulate Climate Change Induced Impacts on Snowpack Dependent Ecosystems

Field studies that simulate the effects of climate change are important for a predictive understanding of ecosystem responses to a changing environment. Among many concerns, regional warming can result in advanced timing of spring snowmelt in snowpack dependent ecosystems, which could lead to longer snow-free periods and drier summer soils. Past studies investigating these impacts of climate change have manipulated snowmelt with a variety of techniques that include manual snowpack alteration with a shovel, infrared radiation, black sand and fabric covers. Within these studies however, sufficient documentation of methods is limited, which can make experimental reproduction difficult. Here, we outline a detailed plot-scale protocol that utilizes a permeable black geotextile fabric deployed on top of an isothermal spring snowpack to induce advanced snowmelt. The method offers a reliable and cost-effective approach to induce snowmelt by passively increasing solar radiation absorption at the snow surface. In addition, control configurations with no snowpack manipulation are paired adjacent to the induced snowmelt plot for experimental comparison. Past and ongoing deployments in Colorado subalpine ecosystems indicate that this approach can accelerate snowmelt by 14-23 days, effectively mimicking snowmelt timing at lower elevations. This protocol can be applied to a variety of studies to understand the hydrological, ecological, and geochemical impacts of regional warming in snowpack dependent ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Smart-PGSim: Using Neural Network to Accelerate AC-OPF Power Grid Simulation

In this work we address the problem of accelerating complex power-grid simulation through machine learning ( ML). Specifically, we develop a framework, Smart-PGSim,which generates multitask-learning (MTL) neural network (NN)models to predict the initial values of variables critical to the problem convergence. MTL models allow information sharing when predicting multiple dependent variables while including customized layers to predict individual variables. We show that,to achieve the required accuracy, it is paramount to embed domain-specific constraints derived from the specific power-grid components in the MTL model. Smart-PGSim then employs the predicted initial values as a high-quality initial condition for the power-grid numerical solver (warm start), resulting in both higher performance compared to state-of-the-art solutions while maintaining the required accuracy. Smart-PGSim brings 2.60×speedup on average (up to 3.28×) computed over 10,000 problems, without losing solution optimality.

machine learning, neural networks↗
DOE OSTI2017

Influence of Alloying Elements and Effect of Stress on Anisotropic Hydrogen Diffusion in Zr-Based Alloys Predicted by Accelerated Kinetic Monte Carlo Simulations

The formation of hydrides degrades the toughness of Zr-based claddings, and thus is one of the primary concerns regarding cladding integrity during fuel operation and used fuel storage. However, the formation mechanism of hydrides, particularly the nucleation mechanism, is yet to be fully understood. To explore the nucleation process, atomic scale studies including density functional theory, molecular dynamics and kinetic Monte Carlo have been used to investigate H diffusion, segregation, and the possible nucleation paths of hydrides. Here, density functional theory calculations are used to parameterize a kinetic Monte Carlo model, which can accurately predict the effective H diffusivity in various Zr-based claddings including Zircaloy, Zirlo and Zr-2.5Nb. Assisted by density functional theory calculations, a recently developed charge-optimized-many-body potential is further improved and used to investigate the possible nucleation paths of hydrides. It’s shown that the formation of face-centered hydrides starts with coherent zeta-phase, shearing of which in the basal plane changes the hcp stacking into fcc and leads to the formation of gamma-hydride nuclei. The shearing involves a negligible barrier and negligible strain accumulation in the matrix by adopting three equivalent shear partials. The results obtained at the atomic scale are consistent with previous theory and experiments, and improves the current understanding on how and how fast hydride forms in Zr based cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗