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 667 records · Page 37

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-RK Acceleration of Low-Lifecycle-Carbon Fuel (LLCF) Simulations

Designing modern combustion systems now relies on computer models that predict how changes in design will affect performance. These models have replaced older methods that relied on the designer's intuition or costly and time-consuming physical testing. By using improved models, design cycles can be shortened, and cleaner and more efficient combustion devices can be created. This project aims to improve computer simulations of low-lifecycle-carbon fuels (LLCFs) with the goal of making these simulations faster and more accurate for predicting combustion in vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Accelerating Our Understanding of Supernova Explosion Mechanism via Simulations and Visualizations with GenASiS

Core-collapse supernovae are among the most powerful explosions in the Universe, releasing about 1053 erg of energy on timescales of a few tens of seconds. These explosion events are also responsible for the production and dissemination of most of the heavy elements, making life as we know it possible. Yet exactly how they work is still unresolved. One reason for this is the sheer complexity and cost of a self-consistent, multi-physics, and multi-dimensional core-collapse supernova simulation, which is impractical, and often impossible, even on the largest supercomputers we have available today. To advance our understanding we instead must often use simplified models, teasing out the most important ingredients for successful explosions, while helping us to interpret results from higher fidelity multi-physics models. In this paper we investigate the role of instabilities in the core-collapse supernova environment. We present here simulation and visualization results produced by our code GenASiS.

Budiardja, R. D.↗

General trends of dihedral conformational transitions in a globular protein

In this paper, dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ~3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD-enhanced sampling. Side chains of the same length (particularly Nχ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. In conclusion, these general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins.

59 BASIC BIOLOGICAL SCIENCES↗

Utilizing many-core accelerators for halo and center finding within a cosmology simulation

Efficiently finding and computing statistics about “halos” (regions of high density) are essential analysis steps for N-body cosmology simulations. However, in state-of-the-art simulation codes, these analysis operators do not currently take advantage of the shared-memory data-parallelism available on multi-core and many-core architectures. The Hybrid / Hardware Accelerated Cosmology Code (HACC) is designed as an MPI+X code, but the analysis operators are parallelized only among MPI ranks, because of the difficulty in porting different X implementations (e.g., OpenMP, CUDA) across all architectures on which it is run. In this paper, we present portable data-parallel algorithms for several variations of halo finding and halo center finding algorithms. These are implemented with the PISTON component of the VTK-m framework, which uses Nvidia's Thrust library to construct data-parallel algorithms that allow a single implementation to be compiled to multiple backends to target a variety of multi-core and many-core architectures. Finally, we compare the performance of our halo and center finding algorithms against the original HACC implementations on the Moonlight, Stampede, and Titan supercomputers. The portability of Thrust allowed the same code to run efficiently on each of these architectures. On Titan, the performance improvements using our code have enabled halo analysis to be performed on a very large data set (81923 particles across 16,384 nodes of Titan) for which analysis using only the existing CPU algorithms was not feasible.

Sewell, Christopher↗