Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reversible computation”

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 271 records · Page 15

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

Source-to-Source Automatic Differentiation of OpenMP Parallel Loops

This article presents our work toward correct and efficient automatic differentiation of OpenMP parallel worksharing loops in forward and reverse mode. Automatic differentiation is a method to obtain gradients of numerical programs, which are crucial in optimization, uncertainty quantification, and machine learning. The computational cost to compute gradients is a common bottleneck in practice. For applications that are parallelized for multicore CPUs or GPUs using OpenMP, one also wishes to compute the gradients in parallel. Here, we propose a framework to reason about the correctness of the generated derivative code, from which we justify our OpenMP extension to the differentiation model. We implement this model in the automatic differentiation tool Tapenade and present test cases that are differentiated following our extended differentiation procedure. Performance of the generated derivative programs in forward and reverse mode is better than sequential, although our reverse mode often scales worse than the input programs.

97 MATHEMATICS AND COMPUTING↗

DAGSENS: Directed Acyclic Graph Based Direct and Adjoint Transient Sensitivity Analysis for Event-Driven Objective Functions

We present DAGSENS, a new theory for parametric transient sensitivity analysis of Differential Algebraic Equation systems (DAEs), such as SPICE-level circuits. The key ideas behind DAGSENS are, (1) to represent the entire sequence of computations, starting from DAE parameters, all the way up to the objective function whose sensitivity is needed, as a Directed Acyclic Graph (DAG) called the "sensitivity DAG", and (2) to compute the required sensitivities efficiently (with time complexity linear in the size of the sensitivity DAG) by leveraging dynamic programming techniques to traverse the DAG. DAGSENS is simple, elegant, and easy-to-understand compared to existing sensitivity analysis approaches; for example, in DAGSENS, one can switch between direct and adjoint transient sensitivities just by changing the direction of DAG traversal (i.e., topological order vs. reverse topological order). Also, DAGSENS is significantly more powerful than existing sensitivity analysis approaches because it allows one to compute the sensitivities of a much more general class of objective functions, including those defined based on "events" that occur during a transient simulation (e.g., a node voltage crossing a particular threshold, a phase-locked loop (PLL) achieving lock, a signal reaching its maximum/minimum value during a transient run, etc.). In this paper, we apply DAGSENS to compute the sensitivities of important event-driven performance metrics in several real-world electronic and biological applications, including high-speed communication (featuring sub-systems such as I/0 links and PLLs), statistical cell library characterization, and gene expression in Drosophila embryos.

42 ENGINEERING↗

Parallelizable adjoint stencil computations using transposed forward-mode algorithmic differentiation

Algorithmic differentiation (AD) is a tool for generating discrete adjoint solvers, which efficiently compute gradients of functions with many inputs, for example for use in gradient-based optimization. AD is often applied to large computations such as stencil operators, which are an important part of most structured-mesh PDE solvers. Stencil computations are often parallelized, for example by using OpenMP, and optimized by using techniques such as cache-blocking and tiling to fully utilize multicore CPUs and many-core accelerators and GPUs. Differentiating these codes with conventional reverse-mode AD results in adjoint codes that cannot be expressed as stencil operations and may not be easily parallelizable. They thus leave most of the compute power of modern architectures unused. We present a method that combines forward-mode AD and loop transformation to generate adjoint solvers that use the same memory access pattern as the original computation that they are derived from and can benefit from the same optimization techniques. The effectiveness of this method is demonstrated by generating a scalable adjoint CFD solver for multicore CPUs and Xeon Phi accelerators.

algorithmic differentiation↗

Down-selection of biomolecules to assemble “reverse micelle” with perovskites

Abstract Biological molecule-semiconductor interfacing has triggered numerous opportunities in applied physics such as bio-assisted data storage and computation, brain-computer interface, and advanced distributed bio-sensing. The introduction of electronics into biological embodiment is being quickly developed as it has great potential in providing adaptivity and improving functionality. Reciprocally, introducing biomaterials into semiconductors to manifest bio-mimetic functionality is impactful in triggering new enhanced mechanisms. In this study, we utilize the vulnerable perovskite semiconductors as a platform to understand if certain types of biomolecules can regulate the lattice and endow a unique mechanism for stabilizing the metastable perovskite lattice. Three tiers of biomolecules have been systematically tested and the results reveal a fundamental mechanism for the formation of a “reverse-micelle” structure. Systematic exploration of a large set of biomolecules led to the discovery of guiding principle for down-selection of biomolecules which extends the classic emulsion theory to this hybrid systems. Results demonstrate that by introducing biomaterials into semiconductors, natural phenomena typically observed in biological systems can also be incorporated into semiconducting crystals, providing a new perspective to engineer existing synthetic materials.

59 BASIC BIOLOGICAL SCIENCES↗

Discovery of Anion Insertion Electrochemistry in Layered Hydroxide Nanomaterials

Electrode materials which undergo anion insertion are a void in the materials innovation landscape and a missing link to energy efficient electrochemical desalination. In recent years layered hydroxides (LHs) have been studied for a range of electrochemical applications, but to date have not been considered as electrode materials for anion insertion electrochemistry. Here, we show reversible anion insertion in a LH for the first time using Co and Co-V layer hydroxides. By pairing in situ synchrotron and quartz crystal microbalance measurements with a computational unified electrochemical band-diagram description, we reveal a previously undescribed anion-insertion mechanism occurring in Co and Co-V LHs. This proof of concept study demonstrates reversible electrochemical anion insertion in LHs without significant material optimization. These results coupled with our foundational understanding of anion insertion electrochemistry establishes LHs as a materials platform for anion insertion electrochemistry with the potential for future application to electrochemical desalination.

Science & Technology - Other Topics↗

Studying rare nonadiabatic dynamics with transition path sampling quantum jump trajectories

We present a method to study rare nonadiabatic dynamics in open quantum systems using transition path sampling and quantum jump trajectories. As with applications of transition path sampling to classical dynamics, the method does not rely on prior knowledge of transition states or reactive pathways and thus can provide mechanistic insight into ultrafast relaxation processes in addition to their associated rates. In particular, we formulate a quantum path ensemble using the stochastic realizations of an unravelled quantum master equation, which results in trajectories that can be conditioned on starting and ending in particular quantum states. Because the dynamics rigorously obeys detailed balance, rate constants can be evaluated from reversible work calculations in this conditioned ensemble, allowing for branching ratios and yields to be computed in an unbiased manner. We illustrate the utility of this method with three examples: energy transfer in a donor-bridge-acceptor model, and models of photo-induced proton-coupled electron transfer and thermally activated electron transfer. These examples demonstrate the efficacy of path ensemble methods and pave the way for their use in studying complex reactive quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-temperature Majorana fermions in magnet-superconductor hybrid systems

Magnet-superconductor hybrid (MSH) structures represent one of the most promising platforms to realize, control, and manipulate Majorana modes using scanning tunneling methods. By depositing either chains or islands of magnetic atoms on the surface of a conventional, elemental superconductor such as Pb or Re, topological superconducting phases can emerge. They feature either localized Majorana bound states at the chain ends or dispersing chiral Majorana modes at the island's boundary. However, some of these experiments have not reached the spectral resolution to clearly distinguish between topological Majorana and trivial Shiba states due to very small superconducting gap sizes and experiments performed at sub-Kelvin temperatures. Here we consider superconducting substrates with unconventional spin-singlet pairing, including high-temperature d -wave and extended s -wave superconductors. We derive topological phase diagrams and compute edge states for cylinder and island geometries and discuss their properties. Several time-reversal invariant topological superconducting phases of the Zhang-Kane-Mele type are found and discussed. Addiotnally, we review one-dimensional MSH structures and show that parameters to realize topologically nontrivial magnetic chains embedded into a larger, two-dimensional substrate differ from the purely one-dimensional case. Quite generally we find that unconventional superconducting substrates work as well as the conventional s -wave substrates to realize topological phases. In particular, iron-based pnictide and chalcogenide superconductors are the most promising class of substrates for future high-temperature MSH systems.

36 MATERIALS SCIENCE↗

Hydrogen Adsorbents with High Volumetric Density: New Materials and System Projections

A high-capacity, low-cost method for storing hydrogen remains one of the primary barriers to the widespread commercialization of fuel cell vehicles. Storage via adsorption is a promising approach due to its fast kinetics, reversibility, and high gravimetric densities. An unfortunate characteristic of adsorptive storage is that high gravimetric densities typically come at the expense of volumetric density (total basis). The goal of this project is to demonstrate best-in-class adsorbents based on metal-organic frameworks (MOFs) that achieve high volumetric and gravimetric H 2 densities simultaneously, while maintaining reversibility and fast kinetics. This is accomplished using a combination ion of high-throughput computational screening and experimental synthesis and characterization. The performance of the most promising materials are projected to the system level by parameterizing system models developed by the Hydrogen Storage Engineering Center of Excellence (HSECoE).

08 HYDROGEN↗

TU-A-9A-03: Development and Verification of a Forward Model That Assists in Iterative Post-Processing Algorithms Used to Reduce Blur in Compton Backscatter Imaging Systems

Purpose: In Compton Backscatter Imaging (CBI), the source and detector reside on the same side of the patient. We previously demonstrated the applicability of CBI systems for medical purposes using an industrial system. To assist in post-processing images from a CBI system, a forward model based on radiation absorption and scatter principles has been developed. Methods: The forward model was developed in C++ using raytracing to track particles. The algorithm accepts phantoms of any size and resolution to calculate the fraction of incident photons scattered back to the detector, and can perform these calculations for any detector geometry and source specification. To validate the model, results were compared to MCNP-X, which is a Monte Carlo based simulation software, for various combinations of source specifications, detector geometries, and phantom compositions. Results: The model verified that the backscatter signal to the detector was based on three interaction probabilities: a) attenuation of photons going into the phantom, b) Compton scatter of photons toward the detector, and c) attenuation of photons coming out of the phantom. The results from the MCNP-X simulations and the forward model varied from 1 to 5%. This difference was less than 1% for energies higher than 30 keV, but was up to 4% for lower energies. At 50 keV, the difference was less than 1% for multiple detector widths and for both homogeneous and heterogeneous phantoms. Conclusion: As part of the optimization of a medical CBI system, an efficient and accurate forward model was constructed in C++ to estimate the output of CBI system. The model characterized individual components contributing to CBI output and increased computational efficiency over Monte Carlo simulations. It is now used in the development of novel post-processing algorithms that reduce image blur by reversing undesired contribution from outside the region of interest.

07 ISOTOPE AND RADIATION SOURCES↗

Improvement of dose distribution in breast radiotherapy using a reversible transverse magnetic field Linac-MR unit

Purpose: To investigate the improvement in dose distribution in tangential breast radiotherapy using a reversible transverse magnetic field that maintains the same direction of Lorentz force between two fields. The investigation has a potential application in future Linac-MR units. Methods: Computed tomography images of four patients and magnetic fields of 0.25–1.5 Tesla (T) were used for Monte Carlo simulation. Two patients had intact breast while the other two had mastectomy. Simulations of planning and chest wall irradiation were similar to the actual clinical process. The direction of superior-inferior magnetic field for the medial treatment beam was reversed for the lateral beam. Results: For the ipsilateral lung and heart mean doses were reduced by a mean (range) of 45.8% (27.6%–58.6%) and 26.0% (20.2%–38.9%), respectively, depending on various treatment plan setups. The mean V{sub 20} for ipsilateral lung was reduced by 55.0% (43.6%–77.3%). In addition acceptable results were shown after simulation of 0.25 T magnetic field demonstrated in dose-volume reductions of the heart, ipsilateral lung, and noninvolved skin. Conclusions: Applying a reversible magnetic field during breast radiotherapy, not only reduces the dose to the lung and heart but also produces a sharp drop dose volume histogram for planning target volume, because of bending of the path of secondary charged particles toward the chest wall by the Lorentz force. The simulations have shown that use of the magnetic field at 1.5 T is not feasible for clinical applications due to the increase of ipsilateral chest wall skin dose in comparison to the conventional planning while 0.25 T is suitable for all patients due to dose reduction to the chest wall skin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Novel Visual and Analytical Methods in Repurposing Legacy Scientific Code - A Case Study

Scientific computing is dominated by team-authored legacy code that has evolved over decades with the purpose of capturing the evolving understanding of a scientific discipline. Accumulated deprecated code, various optimization techniques, and evolving algorithms lead to convoluted source code that is impractical to reverse engineer using mainstream methods. This prevents codes from being truly repeatable or understandable, which are two of the most essential needs in scientific computing. We refactored a long-standing implementation of a common biosequence alignment algorithm in an effort to reproduce its salient behaviors in usable form. Because of the sheer size and complexity of this code base, we developed custom tools to visualize and manipulate the source code behavior under a variety of conditions. We present a case study of extracting and refactoring the algorithmic core and a novel process of discovery/prototyping/testing using a combination of openly available and custom-built tools. The result is a reduction in code size of over 2 orders of magnitude while reconstructing the key protein alignment function in BLAST

bioinformatics, scientific computing↗

Understanding and Controlling Anionic Electrochemical Activity in High-Capacity Oxides for Next Generation Li-Ion Batteries

Rechargeable Li-ion batteries with higher energy density are in urgent demand to address the global challenge of energy storage. In comparison with anode materials, the relatively low capacity of cathode oxides, which exhibit classical cationic redox activity, has become one of the major bottlenecks to reach higher energy density. Recently, anionic activity, such as oxygen redox reaction, has been discovered in the electrochemical processes, providing extra reversible capacity for certain transition-metal oxides. Consequently, a more complete understanding and precise controlling on anionic electrochemical activity in these high-capacity oxides have become a flourishing, yet challenging subject. Here, this perspective highlights (1) key features of the anionic electrochemical activities; (2) computational and experimental tools to characterize and quantify the anionic activity; and (3) design principles that correlate the chemical and structural compositions with high reversible capacity to accelerate the discovery of novel cathode oxides for next generation Li-ion batteries.

25 ENERGY STORAGE↗

Diffusioosmotic flow reversals due to ion–ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion–ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion–ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion–ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. Furthermore, the mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin Dynamics of SCRPs in ZnO Quantum Dot–Organic Molecule Conjugates Studied with Pulsed EPR

Photogenerated spin-correlated radical pairs (SCRPs) are emerging as promising new candidates for Quantum Information Science applications, where they act as electron spin qubit pairs (SQPs). Historically, these pairs have been mostly investigated in natural photosynthesis and in molecular organic donor–(linker)–acceptor systems. Recently, we have shown that these spin pairs can also be observed by time-resolved electron-paramagnetic resonance (EPR) spectroscopy in hybrid inorganic–organic conjugates. The current study builds on recent work (ACS Nano, 2025, 19, 12194–12207), in which we systematically prepared hybrid inorganic–organic molecule systems with tunable geometries that can host SCRPs/SQPs. Here, we demonstrate using pulsed EPR spectroscopy that we can generate, study, and manipulate the spins in these photogenerated SCRPs hosted on highly tailorable inorganic–organic hybrid systems. Microwave pulse-based spin manipulation is the first step toward developing quantum computing and quantum sensing applications with these promising spin-based qubit materials, which can be reversibly photogenerated in highly spin-polarized states at moderate temperatures. An attractive feature of this new class of hybrid materials is that the g-factor of the unpaired electron spin in the ZnO QD (as part of the SCRP/SQP) can be adjusted using the quantum size effect in these QDs, which allows for the selective addressability of each spin in the SCRP using microwave pulses. Additionally, we also observed a difference between the g-values of the unpaired electron in the ZnO QD in the transient SCRP state and stable photoreduced state. We hypothesize that the unpaired electron in the SCRP is delocalized as a band-like conduction electron, while in the stable photoreduced state, the electron is delocalized in a shallow electron trap state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D Nanomagnetism in Low Density Interconnected Nanowire Networks

Free-standing, interconnected metallic nanowire networks with densities as low as 40 mg/cm 3 have been achieved over centimeter-scale areas, using electrodeposition into polycarbonate membranes that have been ion-tracked at multiple angles. Networks of interconnected magnetic nanowires further provide an exciting platform to explore 3-dimensional nanomagnetism, where their structure, topology, and frustration may be used as additional degrees of freedom to tailor the materials properties. New magnetization reversal mechanisms in cobalt networks are captured by the first-order reversal curve method, which demonstrate the evolution from strong demagnetizing dipolar interactions to intersection-mediated domain wall pinning and propagation, and eventually to shape-anisotropy dominated magnetization reversal. These findings open up new possibilities for 3-dimensional integrated magnetic devices for memory, complex computation, and neuromorphics.

36 MATERIALS SCIENCE↗

Study of wrap mode impact on Pseudomonas aeruginosa motion in the chemotactic field of fungi

An agent-based discrete computational model biologically calibrated to Pseudomonas aeruginosa migration is used to explore the impacts of bacterial reversals and wrap mode on the efficiency of motion in different environments, both with and without chemotaxis. It is first shown that wrap mode increases the exploration of continuous multimodal chemotactic profiles such as those produced by biologically relevant fungal networks. For cells undergoing a run-reverse pattern, it is shown that the bacteria are likely to remain at the first local chemoattractant maximal production site on a hypha they find. However, with wrap mode, the bacteria can more easily escape these local sites to further explore their neighboring environment along the fungi, suggesting that wrap mode may be beneficial for migration along the fungi in liquid. In a different set of simulations of bacterial motion close to an isolated chemotactic source, wrap mode is shown to increase the ability of a bacterium to reorient toward the source while reducing the overall motion required for similar chemotactic efficiency as a run-reverse strategy, suggesting a potential metabolic benefit. In contrast, model simulations show that wrap mode can increase the rate of dispersal of P. aeruginosa in a nonchemotactic environment.

Hansen, Austin [University of California, Riversid↗

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗