Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Bayesian model updating with finite element vs surrogate models: Application to a miter gate structural system

Bayesian finite element (FE) model updating using direct model evaluations of large-scale high-fidelity FE models is extremely computationally expensive. Surrogate models can be used as fast emulators of FE models to accelerate the model calibration process. The physics/mechanics-based FE models are still the underpinning behind the surrogate models. Here, this paper evaluates the loss in accuracy and the gain in computational time while performing Bayesian model updating by using surrogate model evaluations compared to using direct FE model evaluations. This evaluation is crucial before entirely relying on surrogate models in model updating for structural health monitoring (SHM) and damage prognosis (DP) purposes. This paper also demonstrates Bayesian updating and surrogate model construction of large-scale high-fidelity FE models of infrastructure systems. In this regard, the miter gate structural system is considered as the testbed structure. Three predominant damage modes (loss of contact between gate and wall, loss of thickness due to corrosion, and loss of tension in the diagonal rods) are considered for model updating purposes. Bayesian model updating is performed using direct FE evaluations by leveraging parallel computing. Two types of surrogates, namely polynomial chaos expansion (PCE) and Gaussian process regression (GPR), are developed for the miter gate. Model updating is performed again using the trained surrogate models, and the updating results are compared with their counterparts obtained using the direct FE evaluation results. The posterior distribution of the FE model parameters obtained using the trained surrogates are sufficiently accurate with respect to the posterior obtained utilizing the direct FE evaluations. In addition, an approximate 4-fold decrease in the computational time was observed when using surrogate model evaluations instead of direct FE evaluations for model updating.

42 ENGINEERING↗
DOE OSTI2024

Repairable and Reconfigurable Structured Liquid Circuits

The advance of printed electronics is significantly bolstered by the development of liquid-state electronics that overcome the inherent limitations in flexibility and reconfigurability of solid-state electronics. By integrating the biocompatibility and conductivity of sulfonated polyaniline (S-PANI) and phytic acid (PA) with the reconfigurability of structured liquids, highly conductive all-liquid threads are developed. The dense packing and overlap of PA/S-PANI complexes at an oil/water interface promotes in-plane electron transport, and standard four-point probe measurements of PA/S-PANI interfacial assemblies demonstrate enhanced electrical properties. Notably, the rapid jetting of the ink phase into the matrix phase allows for liquid threads to be printed, enabling the fabrication of large-scale, conductive pathways between two electrodes and liquid circuits. Upon mechanical cleavage of the liquid wires, circuits can be broken, but will easily self-repair using an electric field, making this motif useful in the design of switches as well as restoring conductive pathways in series or in parallel. In conclusion, the demonstrated flexibility and reconfigurability these PA/S-PANI wires possess hold significant promise for their practical use in the design of flexible and adaptive bioelectronics that can be repaired on demand, signifying a transformative step in the evolution of liquid electronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Two-stage bulk electron heating in the diffusion region of anti-parallel symmetric reconnection

Electron bulk energization in the diffusion region during anti-parallel symmetric reconnection entails two stages. First, the inflowing electrons are adiabatically trapped and energized by an ambipolar parallel electric field. Next, the electrons gain energy from the reconnection electric field as they undergo meandering motion. These collisionless mechanisms have been described previously, and they lead to highly structured electron velocity distributions. Furthermore, a simplified control-volume analysis gives estimates for how the net effective heating scales with the upstream plasma conditions in agreement with fully kinetic simulations and spacecraft observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Here our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ~40%—through multiscale interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Enhancing quantum annealing accuracy through replication-based error mitigation *

Abstract Quantum annealers like those manufactured by D-Wave Systems are designed to find high quality solutions to optimization problems that are typically hard for classical computers. They utilize quantum effects like tunneling to evolve toward low-energy states representing solutions to optimization problems. However, their analog nature and limited control functionalities present challenges to correcting or mitigating hardware errors. As quantum computing advances towards applications, effective error suppression is an important research goal. We propose a new approach called replication based mitigation (RBM) based on parallel quantum annealing (QA). In RBM, physical qubits representing the same logical qubit are dispersed across different copies of the problem embedded in the hardware. This mitigates hardware biases, is compatible with limited qubit connectivity in current annealers, and is well-suited for currently available noisy intermediate-scale quantum annealers. Our experimental analysis shows that RBM provides solution quality on par with previous methods while being more flexible and compatible with a wider range of hardware connectivity patterns. In comparisons against standard QA without error mitigation on larger problem instances that could not be handled by previous methods, RBM consistently gets better energies and ground state probabilities across parameterized problem sets.

Djidjev, Hristo N. (ORCID:0000000192868824)↗
DOE OSTI2024

Probing Quantum Telecloning on Superconducting Quantum Processors

Quantum information cannot be perfectly cloned, but approximate copies of quantum information can be generated. Quantum telecloning combines approximate quantum cloning, more typically referred to as quantum cloning, and quantum teleportation. Quantum telecloning allows approximate copies of quantum information to be constructed by separate parties, using the classical results of a Bell measurement made on a prepared quantum telecloning state. Quantum telecloning can be implemented as a circuit on quantum computers using a classical coprocessor to compute classical feedforward instructions using if statements based on the results of a midcircuit Bell measurement in real time. We present universal symmetric optimal 1 → M telecloning circuits and experimentally demonstrate these quantum telecloning circuits for M = 2 up to M = 10, natively executed with real-time classical control systems on IBM Quantum superconducting processors, known as dynamic circuits. We perform the cloning procedure on many different message states across the Bloch sphere, on seven IBM Quantum processors, optionally using the error suppression technique X–X sequence digital dynamical decoupling. Two circuit optimizations are utilized: one that removes ancilla qubits for M = 2,3 , and one that reduces the total number of gates in the circuit but still uses ancilla qubits. Parallel single-qubit tomography with maximum likelihood estimation density matrix reconstruction is used in order to compute the mixed-state density matrices of the clone qubits, and clone quality is measured using quantum fidelity. These results present one of the largest and most comprehensive noisy intermediate-scale quantum computer experimental analyses on (single qubit) quantum telecloning to date. The clone fidelity sharply decreases to 0.5 for M > 5 , but for M = 2 , we are able to achieve a mean clone fidelity of up to 0.79 using dynamical decoupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Next generation of extreme-resolution electron beam lithography

In top down nanofabrication research facilities around the world, the direct-write high-resolution patterning tool of choice is overwhelmingly electron beam lithography. Remarkably small features can be written in a variety of polymeric resists [V. R. Manfrinato et al., Nano Lett. 14, 4406 (2014); V. R. Manfrinato, A. Stein, L. Zhang, Y. Nam, K. G. Yager, E. A. Stach, and C. T. Black, Nano Lett. 17, 4562 (2017)]. However, this technology, which in this article the authors will refer to as conventional electron beam lithography (CEBL), is reaching its practical resolution and precision limits [V. R. Manfrinato et al., Nano Lett. 14, 4406 (2014)]. Hydrogen depassivation lithography (HDL) [J. N. Randall, J. W. Lyding, S. Schmucker, J. R. Von Ehr, J. Ballard, R. Saini, and Y. Ding, J. Vac. Sci. Technol. B 27, 2764 (2009); J. N. Randall, J. B. Ballard, J. W. Lyding, S. Schmucker, J. R. Von Ehr, R. Saini, H. Xu, and Y. Ding, Microelectron. Eng. 87, 955 (2010)] is a different version of electron beam lithography that is not limited in resolution and precision in the way that CEBL is. It uses a cold field emitter, a scanning tunneling microscope (STM) tip, to deliver a small spot of electrons on a Si (100) 2 × 1 H-passivated surface to expose a self-developing resist that is a monolayer of H adsorbed to the Si surface. Subnanometer features [S. Chen, H. Xu, K. E. J. Goh, L. Liu, and J. N. Randall, Nanotechnology 23, 275301 (2012)], and even the removal of single H atoms can be routinely accomplished [M. A. Walsh and M. C. Hersam, Annu. Rev. Phys. Chem. 60, 193 (2009)]. It is known that the H desorption process at low biases is a multielectron process [E. Foley, A. Kam, J. Lyding, and P. Avouris, Phys. Rev. Lett. 80, 1336 (1998)], but the tunneling distribution of the electrons from the STM tip to the Si surface lattice is not known. The authors have developed a simple model that demonstrates that the combination of two highly nonlinear processes creates a much higher contrast exposure mechanism than CEBL. Currently, HDL has been used almost exclusively on the Si (100) surface and has a limited number of pattern transfer techniques including Si and Ge patterned epitaxy, selective atomic layer deposition of TiO2 followed by reactive ion etching [J. B. Ballard, T. W. Sisson, J. H. G. Owen, W. R. Owen, E. Fuchs, J. Alexander, J. N. Randall, and J. R. Von Ehr, J. Vac. Sci. Technol. B 31, 06FC01 (2013)], and selective deposition of dopant atoms for quantum devices and materials [Workshop on 2D Quantum MetaMaterials held at NIST, Gaithersburg, MD, April 25–26, 2018, edited by J. Owen and W. P. Kirk]. While the throughput of HDL is very low, going parallel in a big way appears promising [J. N. Randall, J. H. G. Owen, J. Lake, R. Saini, E. Fuchs, M. Mahdavi, S. O. R. Moheimani, and B. C. Schaefer, J. Vac. Sci. Technol. B 36, 6 (2018)]. However, the most exciting aspect of HDL is its atomic-scale resolution and precision, which is key to nanoscale research. The authors see HDL emerging as the ultimate high-resolution patterning tool in top down nanofabrication research facilities.

42 ENGINEERING↗
DOE OSTI2021

Asynchronous Iterative Solvers for Extreme-Scale Computing (Final Report)

This is the final report for the project: Asynchronous Iterative Solvers for Extreme-Scale Computing. This was a collaborative project. This report only covers the activities specific to Georgia Institute of Technology. The project investigated and developed iterative solvers that operate asynchronously, thereby avoiding the high cost of synchronization that is apparent when using standard, synchronous iterative solvers at extreme levels of parallelism.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Laser-Induced Ultrafast Demagnetization and Spin Switching (Final Report)

The information technology revolution demands bigger and faster magnetic storage devices. The funded research focuses on a new frontier - Laser-induced ultrafast all-optical spin switching (AOS), where a single laser pulse can switch spins from one direction to another permanently, without a magnetic field, on a time scale of 1-10 picoseconds, 2-3 orders of magnitude faster than the existing technology. The funded research has successfully resolved a long-standing issue in ultrafast demagnetization whether the lattice contributes to the demagnetization. We developed a high level first-principles theory and employed a massively parallel computer to pin down the details of the physical mechanisms that are inaccessible to other methods. This unique feature allows us to compute the force on each atom in both bulk and monolayer systems. We find that it is indeed true there is a nonzero force on each atom, but the idea of vibration contribution is not crucial to our basic understanding of physics. This paves the way to future experimental, theoretical and technological developments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Using Hardware-In-The-Loop Methodology to Develop Test Systems

Hardware in the Loop (HIL) testing methodologies have become widespread in industry. Typically, they focus on developing control algorithms for systems such as autonomous vehicles or aircraft. An oft overlooked aspect of product development is the design and fabrication of a test system for validating that the product meets requirements. Abstractly, a test system differs little from a control system—testers provide signals to the unit, monitor feedback, and base decisions on the results. While the time scales may differ, the functionalities are conceptually similar. Viewed in this light, it becomes natural to extend HIL approaches to tester development. By replacing a physical unit with a proxy model deployed to a real-time or pseudo real-time target, test systems can be developed in parallel with the design and fabrication of a first production unit. This saves considerable time in the life cycle from conceptual design to realized product. This manuscript demonstrates the process flow using a capacitive discharge unit as an exemplar.

42 ENGINEERING↗
DOE OSTI2023

The Event Horizon Telescope Image of the Quasar NRAO 530

We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5–7, when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A*. At z = 0.902, this is the most distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented angular resolution of ~20 μas, both in total intensity and in linear polarization (LP). We do not detect source variability, allowing us to represent the whole data set with static images. The images reveal a bright feature located on the southern end of the jet, which we associate with the core. The feature is linearly polarized, with a fractional polarization of ~5%–8%, and it has a substructure consisting of two components. Their observed brightness temperature suggests that the energy density of the jet is dominated by the magnetic field. The jet extends over 60 μas along a position angle ~ –28°. It includes two features with orthogonal directions of polarization (electric vector position angle), parallel and perpendicular to the jet axis, consistent with a helical structure of the magnetic field in the jet. The outermost feature has a particularly high degree of LP, suggestive of a nearly uniform magnetic field. Future EHT observations will probe the variability of the jet structure on microarcsecond scales, while simultaneous multiwavelength monitoring will provide insight into the high-energy emission origin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Relativistic Alfvén Turbulence at Kinetic Scales

Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close to k −3 or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Three practical workflow schedulers for easy maximum parallelism

Runtime scheduling and workflow systems are an increasingly popular algorithmic component in HPC because they allow full system utilization with relaxed synchronization requirements. There are so many special-purpose tools for task scheduling, one might wonder why more are needed. Use cases seen on the Summit supercomputer needed better integration with MPI and greater flexibility in job launch configurations. Preparation, execution, and analysis of computational chemistry simulations at the scale of tens of thousands of processors revealed three distinct workflow patterns. A separate job scheduler was implemented for each one using extremely simple and robust designs: file-based, task-list based, and bulk-synchronous. Comparing to existing methods shows unique benefits of this work, including simplicity of design, suitability for HPC centers, short startup time, and well-understood per-task overhead. All three new tools have been shown to scale to full utilization of Summit, and have been made publicly available with tests and documentation. This work presents a complete characterization of the minimum effective task granularity for efficient scheduler usage scenarios. Here, these schedulers have the same bottlenecks, and hence similar task granularities as those reported for existing tools following comparable paradigms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]

Here, we present a numerical method for modeling the mechanical effects of nonlinearly-compliant joints in elasto-plastic media. The method uses a series of strain-rate and stress update algorithms to determine joint closure, slip, and solid stress within computational cells containing multiple “embedded” joints. This work facilitates efficient modeling of nonlinear wave propagation in large spatial domains containing a large number of joints that affect bulk mechanical properties. We implement the method within the massively parallel Lagrangian code GEODYN-L and provide verification and examples. We highlight the ability of our algorithms to capture joint interactions and multiple weakness planes within individual computational cells, as well as its computational efficiency. We also discuss the motivation for developing the proposed technique: to simulate large-scale wave propagation during the Source Physics Experiments (SPE), a series of underground explosions conducted at the Nevada National Security Site (NNSS).

42 ENGINEERING↗
DOE OSTI2016

High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations

We study Chebyshev filter diagonalization as a tool for the computation of many interior eigenvalues of very large sparse symmetric matrices. In this technique the subspace projection onto the target space of wanted eigenvectors is approximated with filter polynomials obtained from Chebyshev expansions of window functions. After the discussion of the conceptual foundations of Chebyshev filter diagonalization we analyze the impact of the choice of the damping kernel, search space size, and filter polynomial degree on the computational accuracy and effort, before we describe the necessary steps towards a parallel high-performance implementation. Because Chebyshev filter diagonalization avoids the need for matrix inversion it can deal with matrices and problem sizes that are presently not accessible with rational function methods based on direct or iterative linear solvers. To demonstrate the potential of Chebyshev filter diagonalization for large-scale problems of this kind we include as an example the computation of the 10{sup 2} innermost eigenpairs of a topological insulator matrix with dimension 10{sup 9} derived from quantum physics applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Orientation of Ferroelectric Domains and Disappearance upon Heating Methylammonium Lead Triiodide Perovskite from Tetragonal to Cubic Phase

We study the spontaneous polarization of the archetypal semiconducting halide perovskite methylammonium lead triiodide (CH 3 NH 3 PbI 3 ) that is currently being investigated for use in thin film solar cells and light-emitting diodes. Using both lateral and vertical piezoresponse force microscopy (PFM) to image polycrystalline thin films, we observed domains in the piezoresponse that reversibly appear and disappear below and above the tetragonal-to-cubic phase transition temperature. Importantly, we observe these domains to exhibit a piezoresponse that is predominantly in-plane for films with the (110) plane oriented parallel to the substrate, providing a measure of the polarization associated with specific crystal planes. We characterize the polarization and its temporal response using both local switching spectroscopy and time-dependent PFM spectra. These data show hysteresis loops with the polarization switching with bias but relaxing back on time scales of several minutes. Furthermore, our results suggest the existence of ferroelectric behavior due to off-center displacement of the Pb 2+ cation, although the local polarization response is complicated by the presence of local ionic and electronic conductivity. Understanding the nature of these domains paves the way for further optimization of optoelectronic devices using CH 3 NH 3 PbI 3 perovskite material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗