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

An autoencoder compression approach for accelerating large-scale inverse problems

Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.

Mathematics↗
DOE OSTI2025

Multiclass Classification Using Bayesian Multivariate Adaptive Regression Splines

We present a new Bayesian model for the problem of multiclass classification. In this model, the probabilities of class membership of a given observation are determined by the mean of a latent Gaussian distribution. The mean functions of this latent distribution consist of combinations of highly flexible basis functions of the inputs: multivariate adaptive regression splines (MARS), first developed for multiple regression. We use reversible jump Markov chain Monte Carlo to make inference on the classification model, including the number of basis functions. We compare the probabilistic classification performance of our proposed approach to existing methods on simulated and benchmark data, and compare uncertainty estimates on simulated data. Our proposed method compares favorably with existing Bayesian and frequentist multiclass classification methods in out-of-sample probabilistic classification, and uncertainty estimation of these probabilistic classifications. We examine the fit of the proposed method to a data set of hurricane storm surge levels near Delaware Bay, US, and conclude that sea level rise is a key contributor to damage delivered by storm surge.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

Metastable metal hydrides such as AlH 3 have many attractive features as hydrogen storage media, but generally require complex reaction schemes for regeneration following H 2 release. Here in this paper, we demonstrate that the highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized Covalent Triazine Framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154 °C, with desorption complete at 250 °C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60 °C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. The resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

Coordination Chemistry↗
DOE OSTI2020

Elucidation of Active Oxygen Sites upon Delithiation of Li 3 IrO 4

Transformational increases in the storage capacity of battery cathodes could be achieved by tapping into the redox activity at oxide ligands in addition to conventional transition metal couples. However, the key signatures that govern such lattice oxygen redox (LOR) have not been ascertained. Li 3 IrO 4 has the largest reversible LOR, rendering it a unique model system. Here, X-ray spectroscopy and computational simulations reveal that LOR in Li 3 IrO 4 is selectively compensated via O sites with three lone pairs, which are activated by Li/Ir disorder. The two-electron LOR can be reversed to regenerate the initial state without unlocking competing bulk reactions observed in many other compounds. We uncover an intricate interplay between stoichiometry, O coordination, and nonbonding states in LOR and pinpoint spectroscopic signatures. This interplay is indispensable for designing materials with 3d metals that fulfill the promise of LOR to overcome the bottlenecks of current cathodes for future implementation in practical batteries.

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

An efficient level set method for tracking many materials

Here, we present an efficient level set method to track an arbitrary number of materials. The algorithm is optimal in the sense that it only needs to store a single unsigned distance-like function and a single integer indicator function, independent of the number of materials or distinct regions being tracked. Furthermore, for smooth velocity fields and smooth interface shape, arbitrarily high order solutions can be demonstrated. For interfaces that are or become kinked, the solution is limited to second-order convergence rates in the L 1 norm and first-order in the L ∞ norm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Derivations of the Total Radiation Belt Electron Content

Abstract We present multiple derivations of the Total Radiation Belt Electron Content (TRBEC), an indicator of the global number of electrons that instantaneously occupy the radiation belts. Derived from electron flux measurements, the TRBEC reduces the spatial information into a scalar quantity that concisely describes global aspects of the system. This index provides a simple, global, and long‐term assessment of the radiation belts that enables systematic analysis. In this work, we examine the TRBEC using the adiabatic invariants of which has been used in previous articles as this coordinate system removes reversible adiabatic effects. We then introduce a new expression to compute the TRBEC using the non‐adiabatic coordinates of , relevant in the contexts of energetic electron precipitation, chorus, and hiss scattering where adiabatic invariant quantities are no longer conserved. From both expressions of the TRBEC we demonstrate that an erroneous factor of that appeared in previous works using the adiabatic derivation led to an overestimate of the reported electron populations. In addition, we quantify electron loss in the outer radiation belt via a case study using the Van Allen Probes data over a 20‐day period from March 2013 specifying particle populations both in terms of the aforementioned adiabatic and non‐adiabatic variables. The total number of electrons in the outer radiation belt reached upwards of electrons at the peak of the storm, a rest mass of roughly 10 g.

Astronomy & Astrophysics↗
DOE OSTI2022

The effect of polymer stiffness on magnetization reversal of magnetorheological elastomers

Ultrasoft magnetorheological elastomers (MREs) offer convenient real-time magnetic field control of mechanical properties that provides a means to mimic mechanical cues and regulators of cells in vitro. Here, we systematically investigate the effect of polymer stiffness on magnetization reversal of MREs using a combination of magnetometry measurements and computational modeling. Poly-dimethylsiloxane-based MREs with Young’s moduli that range over two orders of magnitude were synthesized using commercial polymers Sylgard™ 527, Sylgard 184, and carbonyl iron powder. The magnetic hysteresis loops of the softer MREs exhibit a characteristic pinched loop shape with almost zero remanence and loop widening at intermediate fields that monotonically decreases with increasing polymer stiffness. A simple two-dipole model that incorporates magneto-mechanical coupling not only confirms that micrometer-scale particle motion along the applied magnetic field direction plays a defining role in the magnetic hysteresis of ultrasoft MREs but also reproduces the observed loop shapes and widening trends for MREs with varying polymer stiffnesses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Jan-Dec 2021 IC Project: w20_earthquake (Yearly Report) [Slides]

Large-scale simulations were conducted using the computational tool HOSS to investigate scientific aspects of the time reversal process. The main focus was the identification of damage/defects in wellbores for relevant subsurface applications.

58 GEOSCIENCES↗
DOE OSTI2020

Research Introduction [Slides]

In the topic of Reverse Time Imaging, we proposed a new IC to reduce computation cost but reserve image resolution for distributed sensor networks. For Induced Seismicity in Oklahoma, we analyzed fault stress state analysis at state scale, and we applied machine learning techniques to polarity picking and seismicity rate forecasting. The results provide better understanding of fault properties, stress field, and the relationships among fault, stress state, injections, and potential seismic hazards. Lastly, for Microseismic Monitoring, we detected and located 770 low-frequency events (5-50 Hz): (1) Shallow events are highly clustered, consistent with the pathway from injection well 13-10A to monitoring well; moment tensor inversion shows dominant tensile cracking; (2) Deep events are scattered and show migration pattern to the basement; moment tensors show that most events are shear cracks.

58 GEOSCIENCES↗
DOE OSTI2022

POSTER: Automatic Differentiation of Parallel Loops with Formal Methods

The accompanying poster to this short paper presents a combination of reverse mode AD and formal methods to enable efficient differentiation of (or backpropagation through) shared-memory parallel code. Compared to the state of the art, our approach can more often avoid the need for atomic updates or private data copies during the parallel derivative computation, even in the presence of unstructured or data-dependent data access patterns. This is achieved by gathering information about the memory access patterns from the input program, which is assumed to be correctly parallelized. This information is then used to build a model of assertions in a theorem prover, which can be used to check the safety of shared memory accesses during the parallel derivative computation

Automatic Differentiation↗
DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗
DOE OSTI2023

A conformational equilibrium in the nitrogenase MoFe protein with an α-V70I amino acid substitution illuminates the mechanism of H 2 formation

Study of α-V70I-substituted nitrogenase MoFe protein identified Fe6 of FeMo-cofactor (Fe 7 S 9 MoC-homocitrate) as a critical N 2 binding/reduction site. Freeze-trapping this enzyme during Ar turnover captured the key catalytic intermediate in high occupancy, denoted E 4 (4H), which has accumulated 4[e — /H + ] as two bridging hydrides, Fe2–H–Fe6 and Fe3–H–Fe7, and protons bound to two sulfurs. E 4 (4H) is poised to bind/reduce N 2 as driven by mechanistically-coupled H 2 reductive-elimination of the hydrides. This process must compete with ongoing hydride protonation (HP), which releases H 2 as the enzyme relaxes to state E 2 (2H), containing 2[e — /H + ] as a hydride and sulfur-bound proton; accumulation of E 4 (4H) in α-V70I is enhanced by HP suppression. EPR and 95 Mo ENDOR spectroscopies now show that resting-state α-V70I enzyme exists in two conformational states, both in solution and as crystallized, one with wild type (WT)-like FeMo-co and one with perturbed FeMo-co. These reflect two conformations of the Ile residue, as visualized in a reanalysis of the X-ray diffraction data of α-V70I and confirmed by computations. EPR measurements show delivery of 2[e — /H + ] to the E 0 state of the WT MoFe protein and to both α-V70I conformations generating E 2 (2H) that contains the Fe3–H–Fe7 bridging hydride; accumulation of another 2[e — /H + ] generates E 4 (4H) with Fe2–H–Fe6 as the second hydride. E 4 (4H) in WT enzyme and a minority α-V70I E 4 (4H) conformation as visualized by QM/MM computations relax to resting-state through two HP steps that reverse the formation process: HP of Fe2–H–Fe6 followed by slower HP of Fe3–H–Fe7, which leads to transient accumulation of E 2 (2H) containing Fe3–H–Fe7. In the dominant α-V70I E 4 (4H) conformation, HP of Fe2–H–Fe6 is passively suppressed by the positioning of the Ile sidechain; slow HP of Fe3–H–Fe7 occurs first and the resulting E 2 (2H) contains Fe2–H–Fe6. Importantly, it is this HP suppression in E 4 (4H) that enables α-V70I MoFe to accumulate E 4 (4H) in high occupancy. In addition, HP suppression in α-V70I E 4 (4H) kinetically unmasks hydride reductive-elimination without N 2 -binding, a process that is precluded in WT enzyme.

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

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗
DOE OSTI2021

Computational and experimental search for potential polyanionic K-ion cathode materials

Discovering high-energy cathode materials is critical to construct K-ion batteries for practical applications. Owing to the great success of layered oxides in Li- and Na-ion systems, K layered cathodes have also been investigated in recent years. However, the much larger size of K + compared to Li or Na introduces strong K + –K + interaction within the layer, which results in a sloped voltage profile thereby limiting the specific capacity and operating voltage. In contrast, polyanionic materials with a three-dimensional K + arrangement can effectively mitigate K + –K + interaction. In this work, ten K polyanionic compounds with theoretical capacity >100 mA h g –1 are screened from the Inorganic Crystal Structure Database as potential cathode materials for K-ion batteries. Among the ten proposed compounds, K 2 MnP 2 O 7 , K 2 Mn 2 P 2 O 7 F 2 , K 2 Fe 2 P 2 O 7 F 2 , and K 6 V 2 (PO 4 ) 4 with average voltage <4.5 V are synthesized and evaluated electrochemically. While the re-insertion of K into these compounds is not fully reversible, it may be related to the very high migration barrier that we compute for K ions. In addition, we show the successful synthesis of a series of K 3 V 3–x Cr x (PO 4 ) 4 (x = 0, 1, 2, 3) compounds. Among these, K 3 V 2 Cr(PO 4 ) 4 exhibits the largest reversible capacity, as revealed by the in situ investigation. Lastly, we find that the redox couples in many of these compounds sit at remarkably high potential, even higher than in equivalent Li compounds, which brings both opportunities and challenges in the future research of K polyanion cathodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Design, synthesis, and biological testing of biphenylmethyloxazole inhibitors targeting HIV-1 reverse transcriptase

In this article, we report non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) using a biphenylmethyloxazole pharmacophore. A crystal structure of benzyloxazole 1 was obtained and suggested the potential viability of biphenyl analogues. In particular, 6a, 6b, and 7 turned out to be potent NNRTIs with low-nanomolar activity in enzyme inhibition and infected T-cell assays, and with low cytotoxicity. Though modeling further suggested that analogues with fluorosulfate and epoxide warheads might provide covalent modification of Tyr188, synthesis and testing did not find evidence for this outcome.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Thermodynamics and kinetics of precipitation and austenite reversion during aging of Ti-free grade 300 maraging steel manufactured by laser powder bed fusion (LPBF)

Austenite reversion and intermetallic precipitation are two key phase transformations in maraging steels which determine the strength and ductility. Usually, an optimum aging temperature/time is determined using trial and error approach with inputs from thermodynamic databases and a-posteriori microscopy/hardness evaluations. In this study, a thermo-kinetic model is developed to predict the precipitation and austenite reversion kinetics during various post-fabrication heat treatments in a Ti-free variant of grade 300 maraging steel manufactured using laser powder bed fusion. The thermo-kinetic model was based on the principles of Scheil solidification and the classical nucleation theory. The predictions from the thermo-kinetic model were verified using suitable microscopy techniques, and the results from the thermo-kinetic model agree well with the experimental measurements. Overall, results highlight the fact that an Integrated Computational Materials Engineering approach can be used to predict the precipitation and austenite reversion kinetics, and be used as an input for an efficient design of heat treatments for achieving optimum mechanical properties.

36 MATERIALS SCIENCE↗