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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Long-range current-induced spin accumulation in chiral crystals

Abstract Chiral materials, similarly to human hands, have distinguishable right-handed and left-handed enantiomers which may behave differently in response to external stimuli. Here, we use for the first time an approach based on the density functional theory (DFT)+PAOFLOW calculations to quantitatively estimate the so-called collinear Rashba–Edelstein effect (REE) that generates spin accumulation parallel to charge current and can manifest as chirality-dependent charge-to-spin conversion in chiral crystals. Importantly, we reveal that the spin accumulation induced in the bulk by an electric current is intrinsically protected by the quasi-persistent spin helix arising from the crystal symmetries present in chiral systems with the Weyl spin–orbit coupling. In contrast to conventional REE, spin transport can be preserved over large distances, in agreement with the recent observations for some chiral materials. This allows, for example, the generation of spin currents from spin accumulation, opening novel routes for the design of solid-state spintronics devices.

Roy, Arunesh↗

DOE OSTI2015

Electroluminescence from GeSn heterostructure pin diodes at the indirect to direct transition

The emission properties of GeSn heterostructure pin diodes have been investigated. The devices contain thick (400–600 nm) Ge{sub 1−y}Sn{sub y} i-layers spanning a broad compositional range below and above the crossover Sn concentration y{sub c} where the Ge{sub 1−y}Sn{sub y} alloy becomes a direct-gap material. These results are made possible by an optimized device architecture containing a single defected interface thereby mitigating the deleterious effects of mismatch-induced defects. The observed emission intensities as a function of composition show the contributions from two separate trends: an increase in direct gap emission as the Sn concentration is increased, as expected from the reduction and eventual reversal of the separation between the direct and indirect edges, and a parallel increase in non-radiative recombination when the mismatch strains between the structure components is partially relaxed by the generation of misfit dislocations. An estimation of recombination times based on the observed electroluminescence intensities is found to be strongly correlated with the reverse-bias dark current measured in the same devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Optimization of switch diagnostics on the MAIZE linear transformer driver

The MAIZE Linear Transformer Driver is made of 40 capacitor-switch-capacitor `bricks' connected in parallel. When these 40 bricks are charged to 100-kV and then discharged synchronously, the MAIZE facility generates a 1-MA current pulse with a 100-ns rise time into a matched load impedance. Discharging each of the capacitors in a brick is carried out by the breakdown of a spark-gap switch, a process which results in the emission of light. Monitoring this output light with a fiber optic coupled to a photomultiplier tube (PMT) and an oscilloscope channel provides information on switch performance and timing jitter– whether a switch red early, late, or in phase with the other switches. However, monitoring each switch with a dedicated detector- oscilloscope channel can be problematic for facilities where the number of switches to be monitored (e.g., 40 on MAIZE) greatly exceeds the number of detector-oscilloscope channels available. The technique of using fibers to monitor light emission from switches can be optimized by treating a PMT as a binary digit or bit and using a combinatorial encoding scheme, where each switch is monitored by a unique combination of fiber- PMT-oscilloscope channels simultaneously. By observing the unique combination of ber-PMT-oscilloscope channels that are turned on, the pre-firing or late-firing of a single switch on MAIZE can be identified by as few as six PMT-oscilloscope channels. The number of PMT-oscilloscope channels, N , required to monitor X switches can be calculated by 2 N = X + 1, where the number '2' is selected because the PMT-oscilloscope acts as a bit. Here, we demonstrate the use of this diagnostic technique on MAIZE. In conclusion, we also present an analysis of how this technique could be scaled to monitor the tens of thousands of switches proposed for various next generation pulsed power facilities.

47 OTHER INSTRUMENTATION↗

DOE OSTI2019

Kinetic simulation of a low-pressure helium discharge with comparison to experimental measurements

Modern computational validation efforts rely on comparison of known experimental quantities such as current, voltage, particle densities, and other plasma properties with the same values determined through simulation. A discrete photon approach for radiation transport was recently incorporated into a particle-in-cell/direct simulation Monte Carlo code. As a result, spatially and temporally resolved synthetic spectra may be generated even for non-equilibrium plasmas. The generation of this synthetic spectra lends itself to potentially new validation opportunities. In this work, initial comparisons of synthetic spectra are made with experimentally gathered optical emission spectroscopy. A custom test apparatus was constructed that contains a 0.5 cm gap distance parallel plane discharge in ultra high purity helium gas (99.9999%) at a pressure of 75 Torr. Plasma generation is initiated with the application of a fast rise-time, 100 ns full-width half maximum, 2.0 kV voltage pulse. Transient electrical diagnostics are captured along with time-resolved emission spectra. A one-dimensional simulation is run under the same conditions and compared against the experiment to determine if sufficient physics are included to model the discharge. To sync the current measurements from experiment and simulation, significant effort was undertaken to understand the kinetic scheme required to reproduce the observed features. Additionally, the role of the helium molecule excimer emission and atomic helium resonance emission on photocurrent from the cathode are studied to understand which effect dominates photo-feedback processes. Results indicate that during discharge development, atomic helium resonance emission dominates the photo-flux at the cathode even though it is strongly self-absorbed. A comparison between the experiment and simulation demonstrates that the simulation reproduces observed features in the experimental discharge current waveform. Here, the synthesized spectra from the kinetic method produces more favorable agreement with the experimental data than a simple local thermodynamic equilibrium calculation and is a first step towards using spectra generated from a kinetic method in validation procedures. The results of this study produced a detailed compilation of important helium plasma chemistry reactions for simulating transient helium plasma discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗

DOE OSTI2018

Hard-sphere-like dynamics in highly concentrated alpha-crystallin suspensions

The dynamics of concentrated suspensions of the eye-lens protein alpha crystallin have been measured using x-ray photon correlation spectroscopy. Measurements were made at wave vectors corresponding to the first peak in the hard-sphere structure factor and volume fractions close to the critical volume fraction for the glass transition. Langevin dynamics simulations were also performed in parallel to the experiments. The intermediate scattering function f(q,τ) could be fit using a stretched exponential decay for both experiments and numerical simulations. The measured relaxation times show good agreement with simulations for polydisperse hard-sphere colloids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification

Here, a massively parallel method to build large transition rate matrices from temperature-accelerated molecular dynamics trajectories is presented. Bayesian Markov model analysis is used to estimate the expected residence time in the known state space, providing crucial uncertainty quantification for higher-scale simulation schemes such as kinetic Monte Carlo or cluster dynamics. The estimators are additionally used to optimize where exploration is performed and the degree of temperature acceleration on the fly, giving an autonomous, optimal procedure to explore the state space of complex systems. The method is tested against exactly solvable models and used to explore the dynamics of C15 interstitial defects in iron. Our uncertainty quantification scheme allows for accurate modeling of the evolution of these defects over timescales of several seconds.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Towards Highly Scalable Ab Initio Molecular Dynamics (AIMD) Simulations on the Intel Knights Landing Manycore Processor

The Ab Initio Molecular Dynamics (AIMD) method allows scientists to treat the dynamics of molecular and condensed phase systems while retaining a first-principles-based description of their interactions. This extremely important method has tremendous computational requirements, because the electronic Schr¨odinger equation, approximated using Kohn-Sham Density Functional Theory (DFT), is solved at every time step. With the advent of manycore architectures, application developers have a significant amount of processing power within each compute node that can only be exploited through massive parallelism. A compute intensive application such as AIMD forms a good candidate to leverage this processing power. In this paper, we focus on adding thread level parallelism to the plane wave DFT methodology implemented in NWChem. Through a careful optimization of tall-skinny matrix products, which are at the heart of the Lagrange multiplier and nonlocal pseudopotential kernels, as well as 3D FFTs, our OpenMP implementation delivers excellent strong scaling on the latest Intel Knights Landing (KNL) processor. We assess the efficiency of our Lagrange multiplier kernels by building a Roofline model of the platform, and verify that our implementation is close to the roofline for various problem sizes. Finally, we present strong scaling results on the complete AIMD simulation for a 64 water molecules test case, that scales up to all 68 cores of the Knights Landing processor.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2017

Convolutional networks for vehicle track segmentation

Existing methods to detect vehicle tracks in coherent change detection images, a product of combining two synthetic aperture radar images taken at different times of the same scene, rely on simple, fast models to label track pixels. These models, however, are unable to capture natural track features such as continuity and parallelism. More powerful, but computationally expensive models can be used in offline settings. We present an approach that uses dilated convolutional networks consisting of a series of 3-by-3 convolutions to segment vehicle tracks. The design of our networks considers the fact that remote sensing applications tend to operate in low power and have limited training data. As a result, we aim for small, efficient networks that can be trained end-to-end to learn natural track features entirely from limited training data. We demonstrate that our 6-layer network, trained on just 90 images, is computationally efficient and improves the F-score on a standard dataset to 0.992, up from 0.959 obtained by the current state-of-the-art method.

47 OTHER INSTRUMENTATION↗

DOE OSTI2015

WE-EF-BRA-11: Precision Partial-Tumor Irradiation of Dorsal Rodent Mammary Tumors

Purpose: To introduce a pre-clinical treatment technique on a micro-irradiator to treat specific volumes of dorsal mammary tumors in BALB/c mice while sparing lungs and spine. This technique facilitates pre-clinical investigation of tumor response to sub-optimal radiation treatments in which a portion of the tumor is unirradiated, known as a “marginal miss”. In-vitro data suggests that partial tumor radiations trigger a more aggressive phenotype in non-irradiated, regional tumor cells via bystander effects. As the lung tissue is spared, the impact of marginal miss on the development of pulmonary metastasis may be assessed. Methods: End to end test was performed on three BALB/c mice as proof of concept for larger studies. 1Gy was delivered on the micro-irradiator employing previously unexplored lateral parallel-opposed diamond and/or triangle-shaped beams. The margins of the treatment beam were defined using a combination of tumor palpation, barium fiducial markers, and real-time fluoroscopic images. The dose distribution was independently verified with kilovoltage beam Monte Carlo dose calculations with 7% statistical uncertainty and double exposure images. As a final step, the technique was used in a larger pre-clinical study (15Gy, 36 BALB/c mice) and lung metastasis in response to tumor irradiation of 100%, 50% and 0% was quantified. Results: For the Monte-Carlo dose calculations, the dose volume histograms established a maximum dose within the un-irradiated and radiated portions of the mammary tumor of 0.3Gy and 1.5Gy respectively, with a sharp gradient at the boundary. 100% of the lung volume received less than 0.5Gy. This technique proved suitable for a pre-clinical marginal miss study with 50% more lung metastases in partially-radiated mouse models compared to completely. Conclusion: We have developed a novel treatment technique for partial or full irradiation of dorsal mammary tumors incorporating lung sparing.The technique will be useful for exploring the phenomena of aggressive tumor phenotype that may arise from miss-irradiation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2018

Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Here, a numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $\Omega$ on the extruded domain $\mathcal{C}=\Omega\times[0,\infty)$ following. The resulting degenerate elliptic integer order PDE is then approximated using a hybrid FEM-spectral scheme. Finite elements are used in the direction parallel to the problem domain $\Omega$, and an appropriate spectral method is used in the extruded direction. The spectral part of the scheme requires that we approximate the true eigenvalues of the integer order Laplacian over $\Omega$. We derive an a priori error estimate which takes account of the error arising from using an approximation in place of the true eigenvalues. We further present a strategy for choosing approximations of the eigenvalues based on Weyl's law and finite element discretizations of the eigenvalue problem. The system of linear algebraic equations arising from the hybrid FEM-spectral scheme is decomposed into blocks which can be solved effectively using standard iterative solvers such as multigrid and conjugate gradient. Numerical examples in two and three dimensions suggest that the approach is quasi-optimal in terms of complexity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

NANDFlashSim: High-Fidelity, Microarchitecture-Aware NAND Flash Memory Simulation [NANDFlashSim: Intrinsic Latency Variation Aware NAND Flash Memory System Modeling and Simulation at Microarchitecture Level]

As the popularity of NAND flash expands in arenas from embedded systems to high-performance computing, a high-fidelity understanding of its specific properties becomes increasingly important. Further, with the increasing trend toward multiple-die, multiple-plane architectures and high-speed interfaces, flash memory systems are expected to continue to scale and cheapen, resulting in their broader proliferation. However, when designing NAND-based devices, making decisions about the optimal system configuration is nontrivial, because flash is sensitive to a number of parameters and suffers from inherent latency variations, and no available tools suffice for studying these nuances. The parameters include the architectures, such as multidie and multiplane, diverse node technologies, bit densities, and cell reliabilities. Therefore, we introduce NANDFlashSim, a high-fidelity, latency-variation-aware, and highly configurable NAND-flash simulator, which implements a detailed timing model for 16 state-of-the-art NAND operations. Using NANDFlashSim, we notably discover the following. First, regardless of the operation, reads fail to leverage internal parallelism. Second, MLC provides lower I/O bus contention than SLC, but contention becomes a serious problem as the number of dies increases. Third, many-die architectures outperform many-plane architectures for disk-friendly workloads. Lastly, employing a high-performance I/O bus or an increased page size does not enhance energy savings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

mesoflow [SWR-22-56]

Mesoflow is a continuum scale simulation tool developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex surface morphologies (of catalysts/biomass particles among others) directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPU) compared to a single processor for representative problem sizes (2 million cell mesh).

Sitaraman, Hariswaran↗

DOE OSTI2015

Analyzing and modeling the kinetics of amyloid beta pores associated with Alzheimer’s disease pathology

Amyloid beta (Aβ) oligomers associated with Alzheimer’s disease (AD) form Ca 2+ -permeable plasma membrane pores, leading to a disruption of the otherwise well-controlled intracellular calcium (Ca 2+ ) homeostasis. The resultant up-regulation of intracellular Ca 2+ concentration has detrimental implications for memory formation and cell survival. The gating kinetics and Ca 2+ permeability of Aβ pores are not well understood. We have used computational modeling in conjunction with the ability of optical patch-clamping for massively parallel imaging of Ca 2+ flux through thousands of pores in the cell membrane of Xenopus oocytes to elucidate the kinetic properties of Aβ pores. The fluorescence time-series data from individual pores were idealized and used to develop data-driven Markov chain models for the kinetics of the Aβ pore at different stages of its evolution. Our study provides the first demonstration of developing Markov chain models for ion channel gating that are driven by optical-patch clamp data with the advantage of experiments being performed under close to physiological conditions. As a result, we demonstrate the up-regulation of gating of various Ca 2+ release channels due to Aβ pores and show that the extent and spatial range of such up-regulation increases as Aβ pores with low open probability and Ca 2+ permeability transition into those with high open probability and Ca 2+ permeability.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver

MAESTROeX is a massively parallel, finite-volume C++/F90 solver for low Mach number astrophysical flows. Our code utilizes a low Mach number equation set allowing for more efficient, long-time integration of highly subsonic flows compared to compressible approaches. The recommended range of applicability is for flows where the Mach number does not exceed ~0.1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

JAXtronomy: A JAX port of lenstronomy

Gravitational lensing is a phenomenon where light bends around massive objects, resulting in distorted images seen by an observer. Studying gravitationally lensed systems provides insights into cosmology and astrophysics, including constraints of the expansion rate of the Universe and the distribution of dark matter. Thus, we introduce JAXtronomy, a re-implementation of the gravitational lensing software package lenstronomy (Birrer, 2021; Birrer & Amara, 2018) using JAX (Bradbury et al., 2018). JAX is a Python library that uses an accelerated linear algebra (XLA) compiler to improve the performance of computing software. Our core design principle of JAXtronomy is to maintain an identical API to that of lenstronomy. The main JAX features utilized in JAXtronomy are just-in-time compilation, which can lead to significant reductions in execution time, and automatic differentiation, which allows for the implementation of gradient-based algorithms that were previously impossible. Additionally, JAX allows code to be run on GPUs or parallelized across CPU cores, further boosting the performance of JAXtronomy.

astronomy↗

DOE OSTI2017

A New Superalloy Enabling Heavy Duty Gas Turbine Wheels for Improved Combined Cycle Efficiency

The drive to increase combined cycle turbine efficiency from 62% to 65% for the next-generation advanced cycle requires a new heavy duty gas turbine wheel material capable of operating at 1200°F and above. Current wheel materials are limited by the stability of their major strengthening phase (gamma double prime), which coarsens at temperatures approaching 1200°F, resulting in a substantial reduction in strength. More advanced gamma prime superalloys, such as those used in jet engine turbine disks, are also not suitable due to size constraints; the gamma prime phase overages during the slow cooling rates inherent in processing thick-section turbine wheels. The current program addresses this need by screening two new alloy design concepts. The first concept exploits a gamma prime/gamma double prime coprecipitation reaction. Through manipulation of alloy chemistry, coprecipitation is controlled such that gamma double prime is used only to slow the growth of gamma prime during slow cooling, preventing over-aging, and allowing for subsequent heat treatment to maximize strength. In parallel, phase field modeling provides fundamental understanding of the coprecipitation reaction. The second concept uses oxide dispersion strengthening to improve on two existing alloys that exhibit excellent hold time fatigue crack growth resistance, but have insufficient strength to be considered for gas turbine wheels. Mechanical milling forces the dissolution of starting oxide powders into a metal matrix allowing for solid state precipitation of new, nanometer scale oxides that are effective at dispersion strengthening.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Machine Learning and Understanding for Intelligent Extreme Scale Scientific Computing and Discovery. DOE Workshop Report, January 7–9, 2015, Rockville, MD

Large-scale parallel simulations and data analysis drive scientific discovery across many disciplines. To drive larger and more detailed simulations, and deal with larger data volumes, exascale machines capable of 1018 operations per second are expected within the next five to ten years. However, the complexity of developing and adapting modern simulation codes for these architectures is increasing rapidly. Workloads on exascale machines will be billion-way parallel. Exploiting the full capability of the machine will require carefully assigning application tasks to cores, accelerators, deep memories, and other heterogeneous compute resources, while simultaneously optimizing for time to solution, data movement, power, and resilience. Optimizations that improve performance on one machine may slow down another. Worse, applications themselves are dauntingly complex. Production simulations comprise millions of lines of code and use sophisticated, adaptive algorithms whose performance is input-dependent. Complex workflows can couple multi-physics simulation with data preprocessing and post-processing modules.

97 MATHEMATICS AND COMPUTING↗