Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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 973 records · Page 54

High Performance Computing Based Parallel HIearchical Modal Association Clustering (HPAR HMAC)

For many applications, clustering is a crucial step in order to gain insight into the makeup of a dataset. The best approach to a given problem often depends on a variety of factors, such as the size of the dataset, time restrictions, and soft clustering requirements. The HMAC algorithm seeks to combine the strengths of 2 particular clustering approaches: model-based and linkage-based clustering. One particular weakness of HMAC is its computational complexity. HMAC is not practical for mega-scale data clustering. For high-definition imagery, a user would have to wait months or years for a result; for a 16-megapixel image, the estimated runtime skyrockets to over a decade! To improve the execution time of HMAC, it is reasonable to consider an multi-core implementation that utilizes available system resources. An existing imple-mentation (Ray and Cheng 2014) divides the dataset into N partitions - one for each thread prior to executing the HMAC algorithm. This implementation benefits from 2 types of optimization: parallelization and divide-and-conquer. By running each partition in parallel, the program is able to accelerate computation by utilizing more system resources. Although the parallel implementation provides considerable improvement over the serial HMAC, it still suffers from poor computational complexity, O(N2). Once the maximum number of cores on a system is exhausted, the program exhibits slower behavior. We now consider a modification to HMAC that involves a recursive partitioning scheme. Our modification aims to exploit divide-and-conquer benefits seen by the parallel HMAC implementation. At each level in the recursion tree, partitions are divided into 2 sub-partitions until a threshold size is reached. When the partition can no longer be divided without falling below threshold size, the base HMAC algorithm is applied. This results in a significant speedup over the parallel HMAC.

Patlolla, Dilip R↗

Minimizing development costs for efficient many-core visualization using MCD3

Scientific visualization software increasingly needs to support many-core architectures. However, development time is a significant challenge due to the breadth and diversity of both visualization algorithms and architectures. With this work, we introduce a development environment for visualization algorithms on many-core devices that extends the traditional data-parallel primitive (DPP) approach with several existing constructs and an important new construct: meta-DPPs. We refer to our approach as MCD 3 — Meta-DPPs, Convenience routines, Data management, DPPs, and Devices. The twin goals of MCD 3 are to reduce developer time and to deliver efficient performance on many-core architectures, and our evaluation considers both of these goals. For development time, we study 57 algorithms implemented in the VTK-m software library and determine that MCD 3 leads to significant savings. For efficient performance, we survey ten studies looking at individual algorithms and determine that the MCD 3 hardware-agnostic approach leads to performance comparable to hardware-specific approaches: sometimes better, sometimes worse, and better in the aggregate. In total, we find that MCD 3 is an effective approach for scientific visualization libraries to support many-core architectures.

97 MATHEMATICS AND COMPUTING↗

Fast Deposition of Aligning Edge‐On Polymers for High‐Mobility Ambipolar Transistors

Abstract Fast deposition of aligning ambipolar polymers for high‐performance organic field‐effect transistors (OFETs) and inverter circuits are highly desired for both scientific studies and industry applications. Here, large‐area and ordered polymer films are prepared by a bar‐coating method at a rate of 120 mm s −1 in air. Atomic force microscopy and grazing‐incidence wide‐angle X‐ray scattering analysis indicate uniform edge‐on poly(fluoroisoindigo‐difluorobithiophene‐fluoroisoindigo‐bithiophene) (PFIBI‐BT) in 11.7 ± 1 nm film (≈5 layers). The elongated, uniformly oriented grains can reduce the adverse effects of the grain boundaries and facilitate charge transport in polymers. Furthermore, OFETs based on parallel film show high hole/electron mobilities up to 5.5/4.5 cm 2 V −1 s −1 , which are approximately nine times of the devices prepared by spin‐coating. The gain of the inverter is as high as 174, which is one of the highest values in polymer inventers currently. These results demonstrate that the excellent bipolar performance of few‐layer PFIBI‐BT can be ensured while achieving the compatibility of the experimental process with industrial preparation.

Jiang, Yingying↗

EvoGraph: On-The-Fly Efficient Mining of Evolving Graphs on GPU

With the prevalence of the World Wide Web and social networks, there has been a growing interest in high performance analytics for constantly-evolving dynamic graphs. Modern GPUs provide massive AQ1 amount of parallelism for efficient graph processing, but the challenges remain due to their lack of support for the near real-time streaming nature of dynamic graphs. Specifically, due to the current high volume and velocity of graph data combined with the complexity of user queries, traditional processing methods by first storing the updates and then repeatedly running static graph analytics on a sequence of versions or snapshots are deemed undesirable and computational infeasible on GPU. We present EvoGraph, a highly efficient and scalable GPU- based dynamic graph analytics framework.

dynamic graphs↗

Computing the ankle-brachial index with parallel computational fluid dynamics

The ankle-brachial index (ABI), a ratio of arterial blood pressure in the ankles and upper arms, is used to diagnose and monitor circulatory conditions such as coarctation of the aorta and peripheral artery disease. Computational simulations of the ABI can potentially determine the parameters that produce an ABI indicative of ischemia or other abnormalities in blood flow. However, 0- and 1-D computational methods are limited in describing a 3-D patient-derived geometry. Thus, we present a massively parallel framework for computational fluid dynamics (CFD) simulations in the full arterial system. Additionally, using the lattice Boltzmann method to solve the Navier–Stokes equations, we employ highly parallelized and scalable methods to generate the simulation domain and efficiently distribute the computational load among processors. For the first time, we compute an ABI with 3-D CFD. In this proof-of-concept study, we investigate the dependence of ABI on the presence of stenoses, or narrowed regions of the arteries, by directly modifying the arterial geometry. As a result, our framework enables the computation a hemodynamic factor characterizing flow at the scale of the full arterial system, in a manner that is extensible to patient-specific imaging data and holds potential for treatment planning.

59 BASIC BIOLOGICAL SCIENCES↗

Backpropagation-based learning with local derivative approximation and memory replay in biologically plausible neural systems

When learning, the brain modifies individual synaptic connections to reach a desired behavior. Animal and human brains have been shown to be incredibly capable of learning complex and varied functions across a wide variety of tasks. In recent years, artificial neural networks, inspired by human and animal brains, have shown great capabilities in learning a wide variety of difficult tasks. However, artificial neural networks primarily teach themselves through the use of backpropagation, a learning method which has no clear analogue within the brain. Additionally, Artificial Neural Networks primarily use continuous activation functions, which differ significantly from the spiking neuronal behavior present in the brain. In this paper, we discuss and demonstrate a biologically plausible learning method that approximates backpropagation through two techniques on Spiking Neural Networks. First, we show that the local temporal derivatives that are necessary for backpropagation can be approximately recovered through reconstruction using spike timings. Second, we show that through learning during a sleep phase, inspired by neuroscience research into memory replay, the localized parallel feedback path can learn to approximate the derivative through the forward path weight matrix, thus solving the weight transport problem. Lastly, we demonstrate that the combination of these two methods can approach or exceed the accuracy of backpropagation-based methods for a variety of neuromorphic vision tasks while maintaining biological plausibility.

42 ENGINEERING↗

A decreasing carbon allocation to belowground autotrophic respiration in global forest ecosystems

Belowground autotrophic respiration (RAsoil) depends on carbohydrates from photosynthesis flowing to roots and rhizospheres, and is one of the most important but least understood components in forest carbon cycling. Carbon allocation plays an important role in forest carbon cycling and reflects forest adaptation to changing environmental conditions. However, carbon allocation to RAsoil has not been fully examined at the global scale. To fill this knowledge gap, first, the spatio-temporal patterns of RAsoil from 1981 to 2017 were predicted by a Random Forest (RF) algorithm using the most updated Global Soil Respiration Database (v5) with global environmental variables; second, carbon allocation from photosynthesis to RAsoil (CAB), was calculated as the ratio of RAsoil to gross primary production; and its temporal and spatial patterns were assessed in global forest ecosystems. . Globally, mean RAsoil from forests was 8.9 ± 0.08 Pg C yr-1 (mean ± standard deviation) from 1981 to 2017 with strong spatial variabilities. Temporally, RAsoil increased at a rate of 0.0059 Pg C yr-2, paralleling broader soil respiration changes and indicating increasing carbon respired by roots. Mean CAB was 0.243 ± 0.016 and decreased over time. The temporal trend of CAB varied greatly in space, reflecting uneven responses of CAB to environmental changes. This study is the first attempt to predict global CAB and analyze its temporal and spatial patterns. With the linkage of carbon use efficiency, the developed CAB offers an completely independent approach to quantify global aboveground autotropic respiration spatially and temporally, which could provide crucial insights into carbon flux partition and global carbon cycling under climate change.

Tang, Xiaolu↗

On stability criteria for kinetic magnetohydrodynamics

The existence of a potential energy functional in the zero-Larmor-radius collisionless plasma theory of Kruskal & Oberman (Phys. Fluids, vol. 1, 1958 p. 275), Rosenbluth & Rostoker (Phys. Fluids, vol. 2, 1959, p. 23) allows us to derive easily sufficient conditions for linear stability. However, this kinetic magnetohydrodynamics (KMHD) theory does not have a self-adjointness property, making it difficult to derive necessary conditions. In particular, the standard methods to prove that an instability follows if some trial perturbation makes the incremental potential energy negative, which rely on the self-adjointness of the force operator or on the existence of a complete basis of normal modes, are not applicable to KMHD. This paper investigates KMHD linear stability criteria based on the time evolution of initial-value solutions, without recourse to the classic bounds or comparison theorems of Kruskal–Oberman and Rosenbluth–Rostoker for the KMHD potential energy. The adopted approach does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Most importantly, the investigation of a necessary condition for stability does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes. Thus, it is thereby shown that stability in isothermal ideal-MHD is a sufficient condition for stability in KMHD and that, with a proviso on the long-time behaviour of oscillations about stable equilibria, stability in the double-adiabatic fluid theory, including the variation of the parallel fluid displacement, would be a necessary condition for stability in KMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Single Small-Scale Plutonium Redox Reaction System Yields Three Crystallographically-Characterizable Organoplutonium Complexes

Here, an approach to obtaining substantial amounts of data from a hazardous starting material that can only be obtained and handled in small quantities is demonstrated by the investigation of a single small-scale reaction of cyclooctatetraene, C 8 H 8 , with a solution obtained from the reduction of Cp' 3 Pu (Cp' = C 5 H 4 SiMe 3 ) with potassium graphite. This one reaction coupled with oxidation of a product has provided single-crystal X-ray structural data on three organoplutonium compounds as well as information on redox chemistry thereby demonstrating an efficient route to new reactivity and structural information on this highly radioactive element. The crystal structures were obtained from the reduction of C 8 H 8 by a putative Pu(II) complex, (Cp' 3 Pu II ) 1– , generated in situ, to form the Pu(III) cyclooctatetraenide complex, [K(crypt)][(C 8 H 8 ) 2 Pu III ], 1-Pu , and the tetra(cyclopentadienyl) Pu(III) complex, [K(crypt)][Cp' 4 Pu III ], 2-Pu . Oxidation of the sample of 1-Pu with Ag(I) afforded a third organoplutonium complex that has been structurally characterized for the first time, (C 8 H 8 ) 2 Pu IV , 3-Pu . Complexes 1-Pu and 3-Pu contain Pu sandwiched between parallel (C 8 H 8 ) 2– rings. The (Cp' 4 Pu III ) – anion in 2-Pu features three η 5 -Cp' rings and one η 1 -Cp' ring, which is a rare example of a formal Pu–C η 1 -bond. In addition, this study addresses the challenge of small-scale synthesis imparted by radiological and material availability of transuranium isotopes, in particular that of pure metal samples. A route to an anhydrous Pu(III) starting material from the more readily available Pu IV O 2 was developed to facilitate reproducible syntheses and allow complete spectroscopic analysis of 1-Pu and 2-Pu . Pu IV O 2 was converted to Pu III Br 3 (DME) 2 (DME = CH 3 OCH 2 CH 2 OCH 3 ) and subsequently Pu III Br 3 (THF) x , which was used to independently synthesize 1-Pu , 2-Pu , and 3-Pu .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Impacts of chemical gradients on microbial community structure

Succession of redox processes is sometimes assumed to define a basic microbial community structure for ecosystems with oxygen gradients. In this paradigm, aerobic respiration, denitrification, fermentation and sulfate reduction proceed in a thermodynamically determined order, known as the ‘redox tower’. Here, we investigated whether redox sorting of microbial processes explains microbial community structure at low-oxygen concentrations. We subjected a diverse microbial community sampled from a coastal marine sediment to 100 days of tidal cycling in a laboratory chemostat. Oxygen gradients (both in space and time) led to the assembly of a microbial community dominated by populations that each performed aerobic and anaerobic metabolism in parallel. This was shown by metagenomics, transcriptomics, proteomics and stable isotope incubations. Effective oxygen consumption combined with the formation of microaggregates sustained the activity of oxygen-sensitive anaerobic enzymes, leading to braiding of unsorted redox processes, within and between populations. Analyses of available metagenomic data sets indicated that the same ecological strategies might also be successful in some natural ecosystems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toroidal effect on runaway vortex and avalanche growth rate

The momentum-space dynamics of runaway electrons in a slab geometry, in terms of both the geometry and topological transition of the runaway vortex when synchrotron radiative damping is taking into account, have recently been shown to play a crucial role in runaway mitigation and avoidance. In a tokamak geometry, magnetic trapping arises from parallel motion along the magnetic field that scales as 1/R in strength with R the major radius. Since the transit time for a runaway electron moving along the field is of order 10 -8 sec while the collisional time is of ~ 0.01 sec in ITER-like plasmas, a bounce-averaged formulation can drastically reduce computational cost. Here the LAPS-RFP code’s implementation of a bounce-averaged relativistic Fokker- Planck model, along with the essential physics of synchrotron radiation damping and knock-on collisions, is described. It is found that the magnetic trapping can reduce the volume of the runaway vortex as the momentum-space fluxes are strongly modified inside the trapped-region. As a result, the avalanche growth rate is reduced at off-axis locations. In addition to benchmarking with previous calculations that did not take into account radiation damping, we also clarify how synchrotron radiation damping modifies the avalanche growth rate in a tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coherent electronic-vibrational dynamics in deuterium bromide probed via attosecond transient-absorption spectroscopy

Ultrafast laser excitation can trigger multiplex coherent dynamics in molecules. Here, we report attosecond transient-absorption experiments addressing simultaneous probing of electronic and vibrational dynamics in a prototype molecule, deuterium bromide (DBr), following its strong-field ionization. Electronic and vibrational coherences in the ionic X$^2\Pi_{3/2}$ and X$^2\Pi_{1/2}$ states are characterized in the Br-$3d$ core-level absorption spectra via quantum beats with 12.6-fs and 19.9-fs periodicities, respectively. Polarization scans reveal that the phase of the electronic quantum beats depends on the probe direction, experimentally showing that the coherent electronic motion corresponds to the oscillation of the hole density along the ionization-field direction. The vibrational quantum beats are found to maintain a relatively constant amplitude, whereas the electronic quantum beats exhibit a partial decrease in time. Quantum wave-packet simulations show that the decoherence effect from the vibrational motion is insignificant because of the parallel relation between the X$^2\Pi_{3/2}$ and X$^2\Pi_{1/2}$ potentials. Finally, a comparison between the DBr and HBr results suggests that rotation motion is responsible for the decoherence since it leads to initial alignment prepared by the strong-field ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗

Broadband light-trapping enhancement of graphene absorptivity

Low optical absorptivity of graphene limits the quantum efficiency of graphene-based optoelectronic devices. We demonstrate a broadband enhancement of absorptivity in an architecture of graphene intercalation inside a simple-cubic woodpile photonic crystal. Here, significant light-trapping effect is achieved due to spectrally broad and dense resonances in the architecture, especially modes propagating nearly parallel to the interface between the photonic crystal and the background material. We demonstrate an overall absorptivity ~10 to 17 times larger than the free-standing bare graphene monolayer with bandwidth to central frequency ratio 0.3, and ~5 to 11 times larger with bandwidth to central frequency ratio 1.33. A distributed intercalation of three graphene monolayers proves a ~70% overall absorptivity with bandwidth to central frequency ratio 0.3.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An extraction scheme for future CEBAF FFA based energy upgrade

Jefferson lab is considering an energy increase from current 12 GeV to 22 GeV for its CEBAF accelerator. This will be accomplished by recirculating 5-6 additional turns through two parallel CEBAF LINACs using an FFA arc at each end of the racetrack. The total recirculation turns would be 10 times, the first four turns use present conventional arcs to make the 180-degree bends from one LINAC to the other. However, the last 5-6 turns will all share a single beam line inside two FFA arcs. This reduces the footprint and the cost of the project significantly. On the other hand, having the trajectories of last 5-6 recirculating beams close to each other makes it challenging to extract beams from different passes with different energies. In this paper we will explain our present extraction system for 12 GeV, our challenges and limitations, and a possible extraction solution for the 22 GeV upgrade with the goal of extracting beam at different turns/energies to different experimental halls.

Accelerator Physics↗

Extraction scheme for future CEBAF FFA based energy upgrade

Jefferson lab is considering an energy increase from current 12 GeV to 22 GeV for its CEBAF accelerator. This will be accomplished by recirculating 5-6 additional turns through two parallel CEBAF LINACs using an FFA arc at each end of the racetrack. The total recirculation turns would be 10 times, the first four turns use present conventional arcs to make the 180-degree bends from one LINAC to the other. However, the last 5-6 turns will all share a single beam line inside two FFA arcs. This reduces the footprint and the cost of the project significantly. On the other hand, having the trajectories of last 5-6 recirculating beams close to each other makes it challenging to extract beams from different passes with different energies. In this paper we will explain our present extraction system for 12 GeV, our challenges and limitations, and a possible extraction solution for the 22 GeV upgrade with the goal of extracting beam at different turns/energies to different experimental halls.

Kazimi, Reza↗

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗