Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2026

Spectrally accurate, reverse-mode differentiable bounce-averaging algorithm and its applications

We present a fast, spectrally (exponentially) accurate, automatically differentiable bounce-averaging algorithm that is used to simplify kinetic models. Using this algorithm, implemented in the DESC stellarator optimisation suite, we can perform efficient optimisation of many objectives to improve stellarator performance, such as the effective ripple 𝜖 eff metric for the neoclassical transport coefficient in the low collisionality regime and proxies for energetic particle confinement. For the first time, we optimise a finite-beta stellarator to directly reduce neoclassical ripple transport using reverse-mode differentiation. This ensures the computational cost of differentiation is independent of the number of controllable parameters.

fusion plasma↗
DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗
DOE OSTI2021

IBM Open Science Price - SWAP Gate Challenge

On November 30, 2020, IBM announced an open science price competition aimed at improving the fidelity of a SWAP gate on their 7 qubit Quantum system named ”Casablanca”. The details of the SWAP gate problem was outlined in a Jupyter notebook, which utilizes the IBM developed open source package Qiskit. The participants were allowed to work in teams with up to five people and restricted to only use open source software in solving the problem. Their results were due on April 16, 2021. Our approach to solving the SWAP gate challenge is based on the quantum optimal control techniques implemented in our open source packages Juqbox.jl and Quandary. The results from the optimal control approach is only as good as the accuracy in the description of the quantum system dynamics, characterized by a Hamiltonian model that was provided by IBM. Based on the calibrated control pulses IBM provide for the standard gate set, we developed a reverse engineering approach to calibrate our computational model. Techniques were developed to translate between Qiskit’s pulse representation and the B-spline formulation used in Juqbox and Quandary. The fidelity of the optimized pulse sequences were estimated using Qiskit’s randomized bench marking techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Challenges and Solutions for Leave-One-Out Biosensor Design in the Context of a Rugged Fitness Landscape

The leave-one-out (LOO) green fluorescent protein (GFP) approach to biosensor design combines computational protein design with split protein reconstitution. LOO-GFPs reversibly fold and gain fluorescence upon encountering the target peptide, which can be redefined by computational design of the LOO site. Such an approach can be used to create reusable biosensors for the early detection of emerging biological threats. Enlightening biophysical inferences for nine LOO-GFP biosensor libraries are presented, with target sequences from dengue, influenza, or HIV, replacing beta strands 7, 8, or 11. An initially low hit rate was traced to components of the energy function, manifesting in the over-rewarding of over-tight side chain packing. Also, screening by colony picking required a low library complexity, but designing a biosensor against a peptide of at least 12 residues requires a high-complexity library. This double-bind was solved using a “piecemeal” iterative design strategy. Also, designed LOO-GFPs fluoresced in the unbound state due to unwanted dimerization, but this was solved by fusing a fully functional prototype LOO-GFP to a fiber-forming protein, Drosophila ultrabithorax, creating a biosensor fiber. One influenza hemagglutinin biosensor is characterized here in detail, showing a shifted excitation/emission spectrum, a micromolar affinity for the target peptide, and an unexpected photo-switching ability.

Chemistry↗
NASA NTRS1983

Application of a viscous-inviscid interaction procedure to predict separated flows with heat transfer

A viscous-inviscid interaction procedure is described for predicting heat transfer in separated flows. The separating flow in a rearward-facing step/asymmetric channel expansion is considered. For viscous regions, the boundary-layer momentum and continuity equations are solved inversely in a coupled manner by a finite-difference numerical scheme. The streamwise convective term is altered to permit marching the solution through regions of reversed flow. The inviscid flow is computed by numerically solving the Laplace equation for stream-function in the region bounded by the displacement surfaces used in the inverse boundary-layer solution. The viscous and inviscid solutions are repeated iteratively until the edge velocities obtained from both solutions are in agreement. Predictions using this method compare favorably with experimental data and other predictions.

Hall, E. J.↗
NASA NTRS1990

University of Maryland walking robot: A design project for undergraduate students

The design and construction required that the walking robot machine be capable of completing a number of tasks including walking in a straight line, turning to change direction, and maneuvering over an obstable such as a set of stairs. The machine consists of two sets of four telescoping legs that alternately support the entire structure. A gear-box and crank-arm assembly is connected to the leg sets to provide the power required for the translational motion of the machine. By retracting all eight legs, the robot comes to rest on a central Bigfoot support. Turning is accomplished by rotating the machine about this support. The machine can be controlled by using either a user operated remote tether or the on-board computer for the execution of control commands. Absolute encoders are attached to all motors (leg, main drive, and Bigfoot) to provide the control computer with information regarding the status of the motors (up-down motion, forward or reverse rotation). Long and short range infrared sensors provide the computer with feedback information regarding the machine's relative position to a series of stripes and reflectors. These infrared sensors simulate how the robot might sense and gain information about the environment of Mars.

Olsen, Bob↗
NASA NTRS1992

Integrating aerodynamic surface modeling for computational fluid dynamics with computer aided structural analysis, design, and manufacturing

This presentation will discuss the development of a NASA Geometry Exchange Specification for transferring aerodynamic surface geometry between LeRC systems and grid generation software used for computational fluid dynamics research. The proposed specification is based on a subset of the Initial Graphics Exchange Specification (IGES). The presentation will include discussion of how the NASA-IGES standard will accommodate improved computer aided design inspection methods and reverse engineering techniques currently being developed. The presentation is in viewgraph format.

Thorp, Scott A.↗
NASA NTRS2001

Magnetized Target Fusion Driven by Plasma Liners

For practical applications of magnetized target fusion, standoff drivers to deliver the imploding momentum flux to the target plasma remotely are required. Quasi-spherically converging plasma jets have been proposed as standoff drivers for this purpose. The concept involves the dynamic formation of a quasi-spherical plasma liner by the merging of plasma jets, and the use of the liner so formed to compress a spheromak or a field reversed configuration (FRC). Theoretical analysis and computer modeling of the concept are presented. It is shown that, with the appropriate choice of the flow parameters in the liner and the target, the impact between the liner and the target plasma can be made to be shockless in the liner or to generate at most a very weak shock in the liner. Additional information is contained in the original extended abstract.

Thio, Y. C. Francis↗
NASA NTRS1974

Reversal-bounded multipushdown machines

Several representations of the recursively enumerable (r.e.) sets are presented. The first states that every r.e. set is the homomorphic image of the intersection of two linear context-free languages. The second states that every r.e. set is accepted by an on-line Turing acceptor with two pushdown stores such that in every computation, each pushdown store can make at most one reversal (that is, one change from 'pushing' to 'popping'). It is shown that this automata theoretic representation cannot be strengthened by restricting the acceptors to be deterministic multitape, nondeterministic one-tape, or nondeterministic multicounter acceptors. This provides evidence that reversal bounds are not a natural measure of computational complexity for multitape Turing acceptors.

Baker, B. S.↗
DOE OSTI2021

ARENA: Asynchronous Reconfigurable Accelerator Ring to Enable Data-Centric Parallel Computing

The next generation HPC and data centers are likely to be reconfigurable and data-centric due to the trend of hardware specialization and the emergence of data-driven applications. In this work, we propose ARENA – an asynchronous reconfigurable accelerator ring architecture as a potential scenario on how the future HPC and data centers will be like. Despite using the coarse-grained reconfigurable arrays (CGRAs) as the substrate platform, our key contribution is not only the CGRA-cluster design itself, but also the ensemble of a new architecture and programming model that enables asynchronous tasking across a cluster of reconfigurable nodes, so as to bring specialized computation to the data rather than the reverse. We presume distributed data storage without asserting any prior knowledge on the data distribution. Hardware specialization occurs at runtime when a task finds the majority of data it requires are available at the present node. In other words, we dynamically generate specialized CGRA accelerators where the data reside. The asynchronous tasking for bringing computation to data is achieved by circulating the task token, which describes the dataflow graphs to be executed for a task, among the CGRA cluster connected by a fast ring network. Evaluations on a set of HPC and data-driven applications across different domains show that ARENA can provide better parallel scalability with reduced data movement (53.9 percent). Compared with contemporary compute-centric parallel models, ARENA can bring on average 4.37× speedup. The synthesized CGRAs and their task-dispatchers only occupy 2.93mm 2 chip area under 45nm process technology and can run at 800MHz with on average 759.8mW power consumption. ARENA also supports the concurrent execution of multi-applications, offering ideal architectural support for future high-performance parallel computing and data analytics systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Computer program user's manual for advanced general aviation propeller study

A user's manual is presented for a computer program for predicting the performance (static, flight, and reverse), noise, weight and cost of propellers for advanced general aviation aircraft of the 1980 time period. Complete listings of this computer program with detailed instructions and samples of input and output are included.

Worobel, R.↗
DOE OSTI2014

Steady-state magnetohydrodynamic flow around an unmagnetized conducting sphere

The noncollisional interaction between conducting obstacles and magnetized plasma winds can be found in different scenarios, from the interaction occurring between regions inside galaxy clusters to the interaction between the solar wind and Mars, Venus, and active comets, or even the interaction between Titan and the Saturnian magnetospheric flow. These objects generate, through several current systems, perturbations in the streaming magnetic field leading to its draping around the obstacle's effective conducting surface. Recent observational results suggest that several properties associated with magnetic field draping, such as the location of the polarity reversal layer of the induced magnetotail, are affected by variations in the conditions of the streaming magnetic field. To improve our understanding of these phenomena, we perform a characterization of several magnetic field draping signatures by analytically solving an ideal problem in which a perfectly conducting magnetized plasma (with frozen-in magnetic field conditions) flows around a spherical body for various orientations of the streaming magnetic field. In particular, we compute the shift of the inverse polarity reversal layer as the orientation of the background magnetic field is changed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Differentiable McCormick relaxations

McCormick's classical relaxation technique constructs closed-form convex and concave relaxations of compositions of simple intrinsic functions. These relaxations have several properties which make them useful for lower bounding problems in global optimization: they can be evaluated automatically, accurately, and computationally inexpensively, and they converge rapidly to the relaxed function as the underlying domain is reduced in size. They may also be adapted to yield relaxations of certain implicit functions and differential equation solutions. However, McCormick's relaxations may be nonsmooth, and this nonsmoothness can create theoretical and computational obstacles when relaxations are to be deployed. This article presents a continuously differentiable variant of McCormick's original relaxations in the multivariate McCormick framework of Tsoukalas and Mitsos. Gradients of the new differentiable relaxations may be computed efficiently using the standard forward or reverse modes of automatic differentiation. Furthermore, extensions to differentiable relaxations of implicit functions and solutions of parametric ordinary differential equations are discussed. A C++ implementation based on the library MC++ is described and applied to a case study in nonsmooth nonconvex optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The role of an intramolecular hydrogen bond in the redox properties of carboxylic acid naphthoquinones

A bioinspired naphthoquinone model of the quinones in photosynthetic reaction centers but bearing an intramolecular hydrogen-bonded carboxylic acid has been synthesized and characterized electrochemically, spectroscopically, and computationally to provide mechanistic insight into the role of proton-coupled electron transfer (PCET) of quinone reduction in photosynthesis. The reduction potential of this construct is 370 mV more positive than the unsubstituted naphthoquinone. In addition to the reversible cyclic voltammetry, infrared spectroelectrochemistry confirms that the naphthoquinone/naphthoquinone radical anion couple is fully reversible. Calculated redox potentials agree with the experimental trends arising from the intramolecular hydrogen bond. Molecular electrostatic potentials illustrate the reversible proton transfer driving forces, and analysis of the computed vibrational spectra supports the possibility of a combination of electron transfer and PCET processes. The significance of PCET, reversibility, and redox potential management relevant to the design of artificial photosynthetic assemblies involving PCET processes is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Conversion of acenaphthalene to phenalene via methylation: A theoretical study

Ab initio calculations of the C 13 H 10 and C 13 H 9 potential energy surfaces related to the reaction of 1-acenaphthyl and methyl radicals and secondary isomerization of C 13 H 9 primary radical products have been performed at the chemically accurate G3(MP2,CC)//B3LYP/6-311G** level of theory to unravel the mechanism of conversion of acenaphthalene to phenalene or phenalenyl radical + H. The computed energetics and molecular parameters were utilized in Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) calculations of reaction rate constants and relative product yields. The 1-acenaphthyl + CH 3 reaction is predicted to proceed by a fast radical-radical recombination mechanism and to predominantly produce collisionally stabilized 1-methylacenaphthalene B1 or a C 13 H 9 benzylic radical A1 + H with exothermicities of 106.8 and 25.2 kcal/mol, respectively. The A1 + H channel is preferable at higher temperatures, whereas the stabilization of 1-methylacenaphthalene is favored at higher pressures. The radical A1 can nearly irreversibly interconvert to phenalenyl radical AP via a 22.7 kcal/mol exothermic isomerization process involving formal insertion of the CH 2 group into a C-C bond of the five-member ring leading to the expansion of this ring to a six-member ring. The rate constants for the H addition reaction to phenalenyl radical to form phenalene and the reverse H loss from phenalene were also computed and the strength of the weakest C-H bond in the CH 2 group of phenalene is evaluated as 62.2 kcal/mol. The analysis of the reaction kinetics allowed us to deduce a mechanism for the conversion of 1-acenaphthyl radical to phenalene or phenalenyl radical + H via methylation involving the formation of A1 via a well-skipping channel or stabilization and dissociation of (or H abstraction from) B1 followed by isomerization of A1 to phenalenyl AP, which can add an H atom producing phenalene. As a result, rate constants for the significant elementary reactions are fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the expansion of a five-member ring on a zigzag edge of PAH to a six-member ring by methylation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗