Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computing”

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 253 records · Page 14

Computationally Predicted High-Throughput Free-Energy Phase Diagrams for the Discovery of Solid-State Hydrogen Storage Reactions

The design of multinary solid-state material systems that undergo reversible phase changes via changes in temperature and pressure provides a potential means of safely storing hydrogen. However, fully mapping the stabilities of known or newly targeted compounds relative to competing phases at reaction conditions has previously required many stringent experiments or computationally demanding calculations of each compound’s change in Gibbs energy with respect to temperature, G(T). Here, we have extended the approach of constructing chemical potential phase diagrams based on ΔG f (T) to enable the analysis of phase stability at non-zero temperatures. We first performed density functional theory calculations to compute the formation enthalpies of binary, ternary, and quaternary compounds within several compositional spaces of current interest for solid-state hydrogen storage. Temperature effects on solid compound stability were then accounted for using our recently introduced machine learned descriptor for the temperature-dependent contribution G δ (T) to the Gibbs energy G(T). From these Gibbs energies, we evaluated each compound’s stability relative to competing compounds over a wide range of conditions and show using chemical potential and composition phase diagrams that the predicted stable phases and H2 release reactions are consistent with experimental observations. This demonstrates that our approach rapidly computes the thermochemistry of hydrogen release reactions for compounds at sufficiently high accuracy relative to experiment to provide a powerful framework for analyzing hydrogen storage materials. This framework based on G(T) enables the accelerated discovery of active materials for a variety of technologies that rely on solid-state reactions involving these materials.

08 HYDROGEN↗

A Reversed Impedance-Based Stability Criterion for IBR Grids: Preprint

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗

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↗

Interruption as a test of the user-computer interface

In order to study the effects different logic systems might have on interrupted operation, an algebraic calculator and a reverse polish notation calculator were compared when trained users were interrupted during problem entry. The RPN calculator showed markedly superior resistance to interruption effects compared to the AN calculator although no significant differences were found when the users were not interrupted. Causes and possible remedies for interruption effects are speculated. It is proposed that because interruption is such a common occurrence, it be incorporated into comparative evaluation tests of different logic system and control/display system and that interruption resistance be adopted as a specific design criteria for such design.

Kreifeldt, J. G.↗

Reveal, A General Reverse Engineering Algorithm for Inference of Genetic Network Architectures

Given the immanent gene expression mapping covering whole genomes during development, health and disease, we seek computational methods to maximize functional inference from such large data sets. Is it possible, in principle, to completely infer a complex regulatory network architecture from input/output patterns of its variables? We investigated this possibility using binary models of genetic networks. Trajectories, or state transition tables of Boolean nets, resemble time series of gene expression. By systematically analyzing the mutual information between input states and output states, one is able to infer the sets of input elements controlling each element or gene in the network. This process is unequivocal and exact for complete state transition tables. We implemented this REVerse Engineering ALgorithm (REVEAL) in a C program, and found the problem to be tractable within the conditions tested so far. For n = 50 (elements) and k = 3 (inputs per element), the analysis of incomplete state transition tables (100 state transition pairs out of a possible 10(exp 15)) reliably produced the original rule and wiring sets. While this study is limited to synchronous Boolean networks, the algorithm is generalizable to include multi-state models, essentially allowing direct application to realistic biological data sets. The ability to adequately solve the inverse problem may enable in-depth analysis of complex dynamic systems in biology and other fields.

Liang, Shoudan↗

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↗

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↗

Pioneer 9 and OGO 5 observations of an interplanetary multiple shock ensemble on February 2, 1969

A multiple shock system was observed upstream (0.13 AU) of the earth by Pioneer 9 on February 2, 1969. The same system was observed at earth by Ogo 5 and was reported separately in the literature. This paper compares the two sets of observations in still further detail. Both magnetic-field and plasma data are used in a least-squares best-fit method to compute the characteristics of the fast forward shock wave (Pioneer 9 only) and two fast reverse shock waves. Nearly all major features (shock, piston, and tangential discontinuity) retained their characteristics during the transit of the shock ensemble from Pioneer 9 to Ogo 5. The genesis of the ensemble is believed to be due to a complex stream-stream interaction. A substantial density increase (including a large rise of alpha/proton abundance) at Ogo 5, but unobserved at Pioneer 9, is explained by a sudden meridional shift to a flow from below the ecliptic plane while the streams were en route to earth. This study demonstrates a spatial and temporal plasma inhomogeneity which is superimposed on the persistent major features.

Dryer, M.↗

Cyclic structural analysis of air-cooled gas turbine blades and vanes

The creep fatigue behavior of a fully impingement cooled blade for four cyclic cases was analyzed by using the Elas 55, finite element, nonlinear structural computer program. Expected cyclic lives were calculated by using the method of strainrange partitioning for reversed inelastic strains and time fractions for ratcheted tensile creep strains. Strainrange partitioning was also applied to previous results from a one dimensional cyclic analysis of a film impingement cooled vane. The analyses indicated that strainrange partitioning is more applicable to a constrained airfoil such as the film impingement cooled vane than to the relatively unconstrained fully impingement cooled airfoil.

Kaufman, A.↗

Cyclic structural analyses of air-cooled gas turbine blades and vanes

The creep-fatigue behavior of a fully impingement-cooled blade for four cyclic cases was analyzed by using the Elas 55, finite-element, nonlinear structural computer program. Expected cyclic lives were calculated by using the method of Strainrange Partitioning for reversed inelastic strains and time fractions for ratcheted tensile creep strains. Strainrange Partitioning was also applied to previous results from a one-dimensional cyclic analysis of a film-impingement-cooled vane. The analyses indicated that Strainrange Partitioning is more applicable to a constrained airfoil such as the film-impingement-cooled vane than to the relatively unconstrained fully impingement-cooled airfoil.

Kaufman, A.↗

Calculation of the lateral control of swept and unswept flexible wings of arbitrary stiffness

A method similar to that of NACA rep. 1000 is presented for calculating the effectiveness and the reversal speed of lateral-control devices on swept and unswept wings of arbitrary stiffness. Provision is made for using either stiffness curves and root-rotation constants or structural influence coefficients in the analysis. Computing forms and an illustrative example are included to facilitate calculations by means of the method. The effectiveness of conventional aileron configurations and the margin against aileron reversal are shown to be relatively low for swept wings at all speeds and for all wing plan forms at supersonic speeds.

Diederich, Franklin W↗

A Linear Bicharacteristic FDTD Method

The linear bicharacteristic scheme (LBS) was originally developed to improve unsteady solutions in computational acoustics and aeroacoustics [1]-[7]. It is a classical leapfrog algorithm, but is combined with upwind bias in the spatial derivatives. This approach preserves the time-reversibility of the leapfrog algorithm, which results in no dissipation, and it permits more flexibility by the ability to adopt a characteristic based method. The use of characteristic variables allows the LBS to treat the outer computational boundaries naturally using the exact compatibility equations. The LBS offers a central storage approach with lower dispersion than the Yee algorithm, plus it generalizes much easier to nonuniform grids. It has previously been applied to two and three-dimensional freespace electromagnetic propagation and scattering problems [3], [6], [7]. This paper extends the LBS to model lossy dielectric and magnetic materials. Results are presented for several one-dimensional model problems, and the FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗

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↗

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

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↗