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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

B 2 H 6 splitting on catalytic surfaces and role of BH 3 towards hydrogen spillover

A fundamental understanding of the spillover mechanism is an open and challenging problem and plays an important role in catalysis. In particular, bond-exchange spillover mechanism is considered to be effective for reversible storage and release of hydrogen at near ambient conditions. For this, three critical steps are needed: finding the right support (acceptor), the right catalyst to split H 2 , and ensuring that once H 2 is split, the H atoms can migrate on the surface with the help of secondary catalysts and eventually hydrogenate the entire material. In this paper we address these challenges using density functional theory. Here we show that BH 3 , a secondary catalyst, can be produced by symmetrically splitting its stable precursor, B 2 H 6 , on doped metal-free surfaces such as graphene and h-BN as well as on MOF5. In addition, to reduce computational cost, we develop structural descriptor and predictive model equation to effectively screen potential BH 3 binding sites. Symmetrical splitting of B 2 H 6 on different types of materials can address the hydrogen spillover challenge, making efficient storage of hydrogen possible.

08 HYDROGEN↗
DOE OSTI2023

Computation of the Strain Induced Critical Current Reduction in the 16 T Nb 3 Sn Test Facility Dipole

A test facility dipole is being developed at LBNL, targeting a 16 T field in a 144 mm wide aperture. The magnet uses a block design, with two double-pancake coils. In order to minimize motion under the large Lorentz forces, the coils are preloaded against a thick aluminum shell and iron yoke using bladder and key technology. It is then crucial to verify that the performance of the magnet is not degraded due to strain induced on the Nb 3 Sn conductor during assembly, cool-down and powering. The critical current of extracted strands was measured in a varying background magnetic field and as a function of the applied longitudinal strain. Finite element analysis was used to extract the strain state inside the superconducting strands during magnet assembly and operation. This strain was then compared to the measurements to evaluate potential reversible and irreversible effects on the magnet performances. Here the results suggest that the magnet can reach 16 T with sufficient margin, with no irreversible degradation in the high field region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Generalized Lapse Rate Formulas for Use in Entraining CAPE Calculations

Abstract Common assumptions in temperature lapse rate formulas for lifted air parcels include neglecting mixing, hydrostatic balance, the removal of all condensate once it forms (pseudoadiabatic), and/or the retention of all condensate within the parcel (adiabatic). These formulas are commonly derived from the conservation of entropy, which leads to errors when nonequilibrium mixed-phase condensate is present. To evaluate these assumptions, a new general lapse rate formula is derived from an expression for energy conservation, rather than entropy conservation. This new formula incorporates mixing of the parcel with its surroundings, relaxes the hydrostatic assumption, allows for nonequilibrium mixed-phase condensate, and can be formulated for pseudoadiabatic or adiabatic ascent. The new formula is shown to exactly conserve entropy for reversible ascent. Predictions by the new formula are compared to that of older and less general formulas. The errors in previous formulas arise from the assumption of hydrostatic balance, which results in considerable warm biases due to the neglect of the energy sink from buoyancy. Predictions of ascent with entrainment using the new formula are then compared to parcel properties along trajectories in large eddy simulations. Simulated parcel properties are better predicted by the formula using a diluted analogy to adiabatic ascent, wherein condensate is diluted at the same rate as other parcel properties, than by the diluted analogy to pseudoadiabatic ascent, wherein all condensate is removed. These results suggest that CAPE should be computed with adiabatic, rather than pseudoadiabatic, parcel ascent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2023

Quantum Dot–Organic Molecule Conjugates as Hosts for Photogenerated Spin Qubit Pairs

The inherent spin polarization present in photogenerated spin-correlated radical pairs makes them promising candidates for quantum computing and quantum sensing applications. The spin states of these systems can be probed and manipulated with microwave pulses using electron paramagnetic resonance spectrometers. However, to date, there are no reports on magnetic resonance-based spin measurements of photogenerated spin-correlated radical pairs hosted on quantum dots. In the current work, we prepare dye molecule-inorganic quantum dot conjugates and show that they can produce photogenerated spin-polarized states. The dye molecule, D131, is chosen for its ability to undergo efficient charge separation, and the nanoparticle materials, ZnO quantum dots, are chosen for their promising spin properties. Transient and steady state optical spectroscopy performed on ZnO quantum dot-D131 conjugates shows that reversible photogenerated charge separation is occurring. Transient and pulsed electron paramagnetic resonance experiments are then performed on the photogenerated radical pair, which demonstrate that (1) the radical pair is polarized at moderate temperatures and well modeled by existing theories and (2) the spin states can be accessed and manipulated with microwave pulses. In conclusion, this work opens the door to a new class of promising qubit materials that can be photogenerated in polarized states and hosted by highly tailorable inorganic nanoparticles.

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

Spectroscopy of Twisted Bilayer Graphene Correlated Insulators

Here we analytically compute the scanning tunneling microscopy (STM) signatures of integer-filled correlated ground states of the magic angle twisted bilayer graphene (TBG) narrow bands. After experimentally validating the strong-coupling approach at ± 4 electrons/moiré unit cell, we consider the spatial features of the STM signal for 14 different many-body correlated states and assess the possibility of Kekulé distortion (KD) emerging at the graphene lattice scale. Remarkably, we find that coupling the two opposite graphene valleys in the intervalley-coherent (IVC) TBG insulators does not always result in KD. As an example, we show that the Kramers IVC state and its nonchiral U(4) rotations do not exhibit any KD, while the time-reversal-symmetric IVC state does. Our results, obtained over a large range of energies and model parameters, show that the STM signal and Chern number of a state can be used to uniquely determine the nature of the TBG ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Feedback regulation of crystal growth by buffering monomer concentration

Abstract Crystallization is a ubiquitous means of self-assembly that can organize matter over length scales orders of magnitude larger than those of the monomer units. Yet crystallization is notoriously difficult to control because it is exquisitely sensitive to monomer concentration, which changes as monomers are depleted during growth. Living cells control crystallization using chemical reaction networks that offset depletion by synthesizing or activating monomers to regulate monomer concentration, stabilizing growth conditions even as depletion rates change, and thus reliably yielding desired products. Using DNA nanotubes as a model system, here we show that coupling a generic reversible bimolecular monomer buffering reaction to a crystallization process leads to reliable growth of large, uniformly sized crystals even when crystal growth rates change over time. Buffering could be applied broadly as a simple means to regulate and sustain batch crystallization and could facilitate the self-assembly of complex, hierarchical synthetic structures.

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

Mechanistic Pathways for N 2 O Elimination from trans -R 3 Sn-O-N=N-O-SnR 3 and for Reversible Binding of CO 2 to R 3 Sn-O-SnR 3 (R = Ph, Cy)

The rate and mechanism of the elimination of N 2 O from trans-R 3 Sn-O-N=N-O-SnR 3 (R = Ph ( 1 Ph ) and R = Cy ( 1 Cy )) to form R 3 Sn-O-SnR 3 (R = Ph ( 2 Ph ) and R = Cy ( 2 Cy )) have been studied using both NMR and IR techniques to monitor the reactions in the temperature range of 39–79 °C in C 6 D 6 . Activation parameters for this reaction are ΔH ‡ = 15.8 ± 2.0 kcal·mol –1 and ΔS ‡ = –28.5 ± 5 cal·mol –1 ·K –1 for 1 Ph and ΔH ‡ = 22.7 ± 2.5 kcal·mol –1 and ΔS ‡ = –12.4 ± 6 cal·mol –1 ·K –1 for 1 Cy . Addition of O 2 , CO 2 , N 2 O, or PPh 3 to sealed tube NMR experiments did not alter in a detectable way the rate or product distribution of the reactions. Computational DFT studies of elimination of hyponitrite from trans-Me 3 Sn-O-N=N-O-SnMe 3 ( 1 Me ) yield a mechanism involving initial migration of the R 3 Sn group from O to N passing through a marginally stable intermediate product and subsequent N 2 O elimination. Reactions of 1 Ph with protic acids HX are rapid and lead to formation of R 3 SnX and trans-H 2 N 2 O 2 . Reaction of 1 Ph with the metal radical •Cr(CO) 3 C 5 Me 5 at low concentrations results in rapid evolution of N 2 O. At higher •Cr(CO) 3 C 5 Me 5 concentrations, evolution of CO 2 rather than N 2 O is observed. Addition of 1 atm or less CO 2 to benzene or toluene solutions of 2 Ph and 2 Cy resulted in very rapid reaction to form the corresponding carbonates R 3 Sn-O-C(=O)-O-SnR 3 (R = Ph ( 3 Ph ) and R = Cy ( 3 Cy )) at room temperature. Evacuation results in fast loss of bound CO 2 and regeneration of 2 Ph and 2 Cy . Variable temperature data for formation of 3 Cy yield ΔH o = –8.7 ± 0.6 kcal·mol –1 , ΔS o = –17.1 ± 2.0 cal·mol –1 ·K –1 , and ΔG o 298K = –3.6 ± 1.2 kcal·mol –1 . Furthermore, DFT studies were performed and provide additional insight into the energetics and mechanisms for the reactions.

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

Dimetalloylene (M‐E‐M) Complexes of Heavier Main Group Elements Ge, Sn, Pb, Bi via Cleavage of E‐X Bonds (X=N(SiMe 3 ) 2 , O t Bu) with an Iridium Hydride

Abstract Reactions of the Ir V hydride [ Me BDI Dipp ]IrH 4 {BDI=(Dipp)NC(Me)CH(Me)CN(Dipp); Dipp=2,6‐ i Pr 2 C 6 H 3 } with E[N(SiMe 3 ) 2 ] 2 (E=Sn, Pb) afforded the unusual dimeric dimetallotetrylenes ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 in good yields. Moreover, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Ge) 2 was formed in situ from thermal decomposition of [ Me BDI Dipp ]Ir(H) 2 Ge[N(SiMe 3 ) 2 ] 2 . These reactions are accompanied by liberation of HN(SiMe 3 ) 2 and H 2 through the apparent cleavage of an E−N(SiMe 3 ) 2 bond by Ir−H. In a reversal of this process, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 reacted with excess H 2 to regenerate [ Me BDI Dipp ]IrH 4 . Varying the concentrations of reactants led to formation of the trimeric ([ Me BDI Dipp ]IrH 2 ) 3 ( μ 2 ‐E) 3 . The further scope of this synthetic route was investigated with group 15 amides, and ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Bi) 2 was prepared by the reaction of [ Me BDI Dipp ]IrH 4 with Bi(NMe 2 ) 3 or Bi(O t Bu) 3 to afford the first example of a “naked” two‐coordinate Bi atom bound exclusively to transition metals. A viable mechanism that accounts for the formation of these products is proposed. Computational investigations of the Ir 2 E 2 (E=Sn, Pb) compounds characterized them as open‐shell singlets with confined nonbonding lone pairs at the E centers. In contrast, Ir 2 Bi 2 is characterized as having a closed‐shell singlet ground state.

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

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Signal amplification in a solid-state sensor through asymmetric many-body echo

Electronic spins of nitrogen–vacancy centres in diamond constitute a promising system for micro- and nanoscale magnetic sensing, because of their operation under ambient conditions, ease of placement in close proximity to sensing targets and biological compatibility. At high densities, the electronic spins interact through dipolar coupling, which typically limits but can also potentially enhance sensing performance. Here we report the experimental demonstration of many-body signal amplification in a solid-state, room-temperature quantum sensor. Our approach uses time-reversed two-axis-twisting interactions, engineered through dynamical control of the quantization axis and Floquet engineering in a two-dimensional ensemble of nitrogen–vacancy centres. We observe that optimal amplification occurs when the backward evolution time equals twice the forward evolution time, in sharp contrast to the conventional Loschmidt echo. These observations can be understood as resulting from an underlying time-reversed mirror symmetry of the microscopic dynamics, providing key insights into signal amplification and opportunities for practical entanglement-enhanced quantum sensing.

quantum information↗
DOE OSTI2023

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

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

Arctic Tipping Points Triggering Global Change (LDRD Final Report)

The Arctic is warming and feedbacks in the coupled Earth system may be driving the Arctic to tipping events that could have critical downstream impacts for the rest of the globe. In this project we have focused on analyzing sea ice variability and loss in the coupled Earth system Summer sea ice loss is happening rapidly and although the loss may be smooth and reversible, it has significant consequences for other Arctic systems as well as geopolitical and economic implications. Accurate seasonal predictions of sea ice minimum extent and long-term estimates of timing for a seasonally ice-free Arctic depend on a better understanding of the factors influencing sea ice dynamics and variation in this strongly coupled system. Under this project we have investigated the most influential factors in accurate predictions of September Arctic sea ice extent using machine learning models trained separately on observational data and on simulation data from five E3SM historical ensembles. Monthly averaged data from June, July, and August for a selection of ice, ocean, and atmosphere variables were used to train a random forest regression model. Gini importance measures were computed for each input feature with the testing data. We found that sea ice volume is most important earlier in the season (June) and sea ice extent became a more important predictor closer to September. Results from this study provide insight into how feature importance changes with forecast length and illustrates differences between observational data and simulated Earth system data. We have additionally performed a global sensitivity analysis (GSA) using a fully coupled ultra- low resolution configuration E3SM. To our knowledge, this is the first global sensitivity analysis involving the fully-coupled E3SM Earth system model. We have found that parameter variations show significant impact on the Arctic climate state and atmospheric parameters related to cloud parameterizations are the most significant. We also find significant interactions between parameters from different components of E3SM. The results of this study provide invaluable insight into the relative importance of various parameters from the sea ice, atmosphere and ocean components of the E3SM (including cross-component parameter interactions) on various Arctic-focused quantities of interest (QOIs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Non-topotactic reactions enable high rate capability in Li-rich cathode materials

High-rate cathode materials for Li-ion batteries require fast Li transport kinetics, which typically rely on topotactic Li intercalation/de-intercalation because it minimally disrupts Li transport pathways. In contrast to this conventional view, here we demonstrate that the rate capability in a Li-rich cation-disordered rocksalt cathode can be significantly improved when the topotactic reaction is replaced by a non-topotactic reaction. The fast non-topotactic lithiation reaction is enabled by facile and reversible transition metal octahedral-to-tetrahedral migration, which improves rather than impedes Li transport. Using this concept, we show that high-rate performance can be achieved in Mn- and Ni-based cation-disordered rocksalt materials when some of the transition metal content can reversibly switch between octahedral and tetrahedral sites. Overall, this study provides a new perspective on the design of high-performance cathode materials by demonstrating how the interplay between Li and transition metal migration in materials can be conducive to fast non-topotactic Li intercalation/de-intercalations.

25 ENERGY STORAGE↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗