Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel in Time”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Focal Hydrothermal Ablation: Preliminary Investigation of a New Concept

PurposeTo determine whether focal tissue ablation is possible with interstitial instillation of steam.MethodsFresh swine livers were used. Through a 20 gauge needle, steam was instilled every 5 s, 3 (n = 5), 6 (n = 5), 9 (n = 5), or 12 (n = 5 + 5) times in a liver lobe. The ablated zones were sectioned parallel (n = 20) or perpendicular (n = 5) to the needle track. The longitudinal long and short axis diameters, or transverse long and short axis diameters of areas with discoloration on macroscopic examination, were measured. The experiment was repeated in vivo on a pig. Steam instillation was performed once every 5 s for 5 min in the liver (n = 3) and in muscle (n = 4), and temperature changes at three neighboring sites were monitored. Long and short axis diameters of the discolored areas were measured.ResultsA well-defined area of discoloration was invariably present at the site of steam instillation. The median longitudinal long axis diameter were 2.0, 2.5, 2.5, and 3.5 cm for 3, 6, 9, and 12 steam instillations in vitro, while median short axis diameters were 1.0, 1.5, 1.5, and 1.5 cm, respectively. Six attempts at ablation in vivo could be successfully completed. The long axis diameters of the ablated zones in the liver were 7.0 and 8.0 cm, while in muscle it ranged from 5.5 to 7.0 cm.ConclusionInstillation of steam in the liver in vitro and in vivo, and in muscle in vivo rapidly leads to circumscribed zones of coagulation necrosis.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Unusual layer-by-layer growth of epitaxial oxide islands during Cu oxidation

Elucidating metal oxide growth mechanisms is essential for precisely designing and fabricating nanostructured oxides with broad applications in energy and electronics. However, current epitaxial oxide growth methods are based on macroscopic empirical knowledge, lacking fundamental guidance at the nanoscale. Using correlated in situ environmental transmission electron microscopy, statistically-validated quantitative analysis, and density functional theory calculations, we show epitaxial Cu 2 O nano-island growth on Cu is layer-by-layer along Cu 2 O(110) planes, regardless of substrate orientation, contradicting classical models that predict multi-layer growth parallel to substrate surfaces. Growth kinetics show cubic relationships with time, indicating individual oxide monolayers follow Frank-van der Merwe growth whereas oxide islands follow Stranski-Krastanov growth. Cu sources for island growth transition from step edges to bulk substrates during oxidation, contrasting with classical corrosion theories which assume subsurface sources predominate. Our results resolve alternative epitaxial island growth mechanisms, improving the understanding of oxidation dynamics critical for advanced manufacturing at the nanoscale.

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

Glassy Dynamics in a heavy ion irradiated NbSe 2 crystal

Fascination with glassy states has persisted since Fisher introduced the vortex-glass as a new thermodynamic phase that is a true superconductor that lacks conventional long-range order. Though Fisher’s original model considered point disorder, it was later predicted that columnar defects (CDs) could also induce glassiness — specifically, a Bose-glass phase. In YBa 2 Cu 3 O 7–x (YBCO), glassy states can cause distinct behavior in the temperature (T) dependent rate of thermally activated vortex motion (S). The vortex-glass state produces a plateau in S(T) whereas a Bose-glass can transition into a state hosting vortex excitations called double-kinks that can expand, creating a large peak in S(T). Although glass phases have been well-studied in YBCO, few studies exist of other materials containing CDs that could contribute to distinguishing universal behavior. Here, we report on the effectiveness of CDs tilted ~30° from the c-axis in reducing S in a NbSe 2 crystal. The magnetization is 5 times higher and S is minimized when the field is parallel to the defects versus aligned with the c-axis. We see signatures of glassiness in both field orientations, but do not observe a peak in S(T) nor a plateau at values observed in YBCO. Lastly, we discuss the possibility that competing disorder induces a field-orientation-driven transition from a Bose-glass to an anisotropic glass involving both point and columnar disorder.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
DOE OSTI2013

Reduced Insulin/Insulin-like Growth Factor-1 Signaling and Dietary Restriction Inhibit Translation but Preserve Muscle Mass in Caenorhabditis elegans

Reduced signaling through the C. elegans insulin/IGF1 like tyrosine kinase receptor daf2 and dietary restriction via bacterial dilution are two well-characterized lifespan-extending interventions that operate in parallel or through (partially) independent mechanisms. Using accurate mass and time tag LCMS/MS quantitative proteomics we detected that the abundance of a large number of ribosomal subunits is decreased in response to dietary restriction as well as in the daf2(e1370) insulin/IGF1 receptor mutant. In addition, general protein synthesis levels in these long-lived worms are repressed. Surprisingly, ribosomal transcript levels were not correlated to actual protein abundance, suggesting that posttranscriptional regulation determines ribosome content. Proteomics also revealed increased presence of many structural muscle cell components in long-lived worms, which appears to result from prioritized preservation of muscle cell volume in nutrient-poor conditions or low insulin-like signaling. Activation of DAF16, but not diet-restriction, stimulates mRNA expression of muscle-related genes to prevent muscle atrophy. Important daf2 specific proteome changes include overexpression of aerobic metabolism enzymes and a general activation of stress responsive and immune defense systems, while increased abundance of many protein subunits of the proteasome core complex is a DR-specific characteristic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

High-order-harmonic generation of vibrating H 2 + and D 2 +

In this work, we study the isotope effect on high-order-harmonic generation (HHG) of vibrating H 2 + and D 2 + molecular ions aligned parallel to the polarization of the laser field. The time-dependent Schrödinger equation is solved accurately and efficiently by the time-dependent generalized pseudospectral and Fourier grid methods in three spatial coordinates, one of them being the internuclear separation and the other two describing the electronic motion. The laser pulses have a carrier wavelength of 800 nm and duration of 10 or 16 optical cycles. The peak intensities used in the calculations are 2×10 14 and 3×10 14 W/cm 2 . The effect of nuclear vibration is visible in both H 2 + and D 2 + but more pronounced in the lighter H 2 + molecule. Striking differences from the fixed nuclei case are a total disappearance of the traditional plateau in the HHG spectrum at the higher intensity and significant redshift of the harmonic peaks in the central part of the spectrum. These phenomena are explained based on the analysis of the dynamics of the nuclear vibrational wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Parallel general purpose multiobjective optimization framework with application to electron beam dynamics

Particle accelerators are invaluable tools for research in the basic and applied sciences, such as materials science, chemistry, the biosciences, particle physics, nuclear physics and medicine. The design, commissioning, and operation of accelerator facilities is a nontrivial task, due to the large number of control parameters and the complex interplay of several conflicting design goals. The Argonne Wakefield Accelerator facility has some unique challenges resulting from its purpose to carry out advanced accelerator R&D. Individual experiments often have challenging beam requirements, and the physical configuration of the beam lines is often changed to accommodate the variety of supported experiments. The need for rapid deployment of different operational settings further complicates the optimization work that must be done for multiple constraints and challenging operational regimes. One example of this is an independently staged two-beam acceleration experiment which requires the construction of an additional beam line (this is now in progress). The high charge drive beam, well into the space charge regime, must be threaded through small aperture (17.6 mm) decelerating structures. In addition, the bunch length must be sufficiently short to maximize power generation in the decelerator. We propose to tackle this problem by means of multiobjective optimization algorithms which also facilitate a parallel deployment. In order to compute solutions in a meaningful time frame, a fast and scalable software framework is required. In this paper, we present a general-purpose framework for simulation-based multiobjective optimization methods that allows the automatic investigation of optimal sets of machine parameters. Using evolutionary algorithms as the optimizer and opal as the forward solver, validation experiments and results of multiobjective optimization problems in the domain of beam dynamics are presented. Optimized solutions for the new high charge drive beam line found by the framework were used to finish the design of a two beam acceleration experiment. The selected solution along with the associated beam parameters is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

A Unit Testing Framework for Scientific Legacy Code

Large-scale scientific applications play important roles in supporting research. However, it is often very expensive and time-consuming to make changes to, maintain and evolve the scientific code due to its complexity and poor programming skills of researchers. Therefore, in order to visualize scientific code architecture to optimize software design, understand undocumented source code, and analyze software flow and functionality, we first introduce a unit testing framework (UTF). Then, because such infrastructure's performance is very crucial in practical use since the scientific legacy applications simulate instances in a long period of time, we improve the UTF by applying Message Passing based Parallelization and parallel I/O operations. Furthermore, due to the scientific code has enormous state data and the I/O capacity on the server is limited, we apply in situ data analysis method to encounter fewer resource limitations, and adopt signal processing to greatly reduce data transfer. Last, we demonstrated the correctness and high-efficiency of our framework for legacy Earth model on Titan supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Accelerating Bilevel Optimization With Hierarchical Many-Threaded Parallel Differential Evolution

Bilevel optimization is encountered in many relevant real-world applications. The main feature of this type of problem is that an upper-level optimization problem is constrained by a nested lower-level optimization problem. Because of this nested structure, bilevel problems (BLPs) are usually computationally expensive to solve. Differential evolution (DE) has demonstrated promising results in solving BLPs of relatively small scales. As the problem scale increases, the decision space becomes intrinsically larger, requiring a growing number of function evaluations for the method to work properly. In this context, heavy parallelization and high-performance computing techniques are indispensable to enable the resolution of more complex and challenging optimization problems. Hence, we propose a hierarchical many-threaded parallel DE approach for BLPs, where both levels are parallelized. The computational experiments demonstrate that the parallel implementation achieved runtime speeds ranging from 44 to 2559 times faster than the sequential version on a well-known scalable SMD benchmark test problem when executed on an NVIDIA A100 GPU. The findings indicate that the algorithm’s convergence is strongly influenced by the number of both upper- and lower-level generations. Moreover, the success of experiments with large-scale problems is closely linked to the choice of small population sizes.

Dufek, Amanda S↗
DOE OSTI2026

Fast and Accurate Intersections on a Sphere

We introduce a fast, high-precision algorithm for calculating intersections between great circle arcs and lines of constant latitude on the unit sphere. We first propose a simplified intersection point formula with improved speed and numerical robustness over the ones traditionally implemented in geoscience software. We then show how algorithms based on the concept of error-free transformations (EFT) can be applied to evaluate this formula within a relative error bound that is on the order of machine precision. Here, we demonstrate that, with a vectorized and parallelized implementation, this enhanced accuracy is achieved with no compute time overhead compared to a direct calculation in hardware floating point, making our algorithm suitable for performance-sensitive applications like regridding of high-resolution climate data. In contrast, evaluating our formula using high-precision data types like quadruple precision and arbitrary precision, or using the robust intersection computation routines from the Computational Geometry Algorithms Library, leads to significant computational overhead, especially since these alternatives inhibit vectorization. More generally, our work demonstrates how EFT techniques can be combined and extended to implement nontrivial geometric calculations with high accuracy and speed.

Environmental sciences↗
DOE OSTI2018

Integrated Model, Batch and Domain Parallelism in Training Neural Networks

We propose a new integrated method of exploiting model, batch and domain parallelism for the training of deep neural networks (DNNs) on large distributed-memory computers using minibatch stochastic gradient descent (SGD). Our goal is to find an efficient parallelization strategy for a fixed batch size using $P$ processes. Our method is inspired by the communication-avoiding algorithms in numerical linear algebra. We see $P$ processes as logically divided into a $P_r \times P_c$ grid where the $P_r$ dimension is implicitly responsible for model/domain parallelism and the $P_c$ dimension is implicitly responsible for batch parallelism. In practice, the integrated matrix-based parallel algorithm encapsulates these types of parallelism automatically. We analyze the communication complexity and analytically demonstrate that the lowest communication costs are often achieved neither with pure model nor with pure data parallelism. We also show how the domain parallel approach can help in extending the theoretical scaling limit of the typical batch parallel method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

KinCat v.1.0

SAND2024-02099O The software is designed to allow researchers to perform kinetic Monte Carlo (KMC) simulations of catalytic reactions on a 2D lattice. The code is written efficiently to run on a variety of shared memory computing architectures (e.g. GPU, multi-core) and to natively express the full complexity of lateral interactions on reaction rates. The software allows researchers to perform KMC simulations of catalytic reactions on a 2D lattice. It uses parallel shared-memory computing architectures to reduce run-times and allows for simultaneous simulation of multiple independent runs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Najm, Habib↗