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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

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↗

NASA NTRS2001

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental-iterative approach for differentiating advanced flow codes is successfully demonstrated on a 2D inviscid model problem. The method employs the reverse-mode capability of the automatic- differentiation software tool ADIFOR 3.0, and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straight-forward, black-box reverse- mode application of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-order aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoint) procedures; then, a very efficient non-iterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hessian matrices) of lift, wave-drag, and pitching-moment coefficients are calculated with respect to geometric- shape, angle-of-attack, and freestream Mach number

Taylor, Arthur C., III↗

NASA NTRS1973

Direct comparison between satellite electric field measurements and the visual aurora

Electric field data from two passes of the Injun 5 satellite, one corresponding to magnetically quiet conditions and one corresponding to substorm conditions, are compared with simultaneous all-sky-camera data from College, Alaska. In each case, a significant deviation of the electric field from the expected V x B field (where V is the satellite velocity) was evident and a distinct electric field reversal could be identified. In the region of substantial electric field equatorward of the electric field reversal a diffuse auroral arc was observed during the magnetically quiet pass and auroral patches were observed during the substorm pass. The motion of the auroral patches was consistent with the general direction and magnitude of the E x B drift computed from the satellite electric field measurements. In the substorm case the electric field reversal occurred very near a discrete auroral arc at the poleward side of the diffuse arcs and patches. Comparison of the quiet time and substorm cases suggests that the convection electric field penetrates deeper into the magnetosphere during a substorm.

Swift, D. W.↗

NASA NTRS1987

An experimental investigation of compressible three-dimensional boundary layer flow in annular diffusers

An experimental study is described in which detailed wall pressure measurements have been obtained for compressible three-dimensional unseparated boundary layer flow in annular diffusers with and without normal shock waves. Detailed mean flow-field data were also obtained for the diffuser flow without a shock wave. Two diffuser flows with shock waves were investigated. In one case, the normal shock existed over the complete annulus whereas in the second case, the shock existed over a part of the annulus. The data obtained can be used to validate computational codes for predicting such flow fields. The details of the flow field without the shock wave show flow reversal in the circumferential direction on both inner and outer surfaces. However, there is a lag in the flow reversal between the inner nad the outer surfaces. This is an interesting feature of this flow and should be a good test for the computational codes.

Om, Deepak↗

NASA NTRS1996

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗

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↗

NASA NTRS1972

Carbon monoxide oxidation rates computed for automobile thermal reactor conditions

Carbon monoxide oxidation rates in thermal reactors for exhaust manifolds are computed by integrating differential equations for system of twenty-nine reversible chemical reactions. Reactors are noncatalytic replacements for conventional exhaust manifolds and are a system for reducing carbon monoxide and hydrocarbons in automobile exhausts.

Brokaw, R. S.↗

NASA NTRS1992

Numerical simulation of STOL operations using thrust-vectoring

The flow field about a delta wing equipped with thrust reverser jets in slow speed flight near the ground has been computed. Results include the prediction of the flow about the delta wing at four fixed heights above the ground, and a simulated landing, in which the delta wing descends towards the ground. Comparison of computed and experimental lift coefficients indicates that the simulations can capture at least the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect.

Chawla, Kalpana↗

NASA NTRS1993

Numerical simulation of a powered-lift landing

The flow field about a delta wing equipped with thrust reverse jets in slow speed flight near the ground has been computed. Results include the prediction of the flow about the delta wing at four fixed heights above the ground, and simulated landing, in which the delta wing descends towards the ground. Comparison of computed and experimental lift coefficients indicates that the simulations can capture at least the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect.

Chawla, Kalpana↗

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↗

NASA NTRS1972

Effect of pictorial display configuration on the frequency of control reversals during aircraft landing approaches

A fixed-base simulator study to determine the effect of increased pictorial display realism on the frequency of control reversals made with an inside-out landing display was conducted. Display conditions included the effects of collimation and scale (head-up versus head-down presentation), horizon symbology (simple line versus white-black sky-ground surfaces), and ground plane realism (computer generated perspective versus a TV picture of a realistic model). The number of control reversals was moderately high on all displays. Control reversals to roll disturbances occurred nearly twice as frequently as reversals to either pitch or lateral rate disturbances. Though there were no significant differences among the numbers of small control reversals for the different displays, there was some evidence that this conclusion may not apply to large control reversals.

Lincoln, R.↗

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↗

DOE OSTI2020

Sequence optimizations in a high-performance computing environment

Embodiments are directed to techniques to determine dataflow graph instructions comprising one or more pick/switch instruction pairs and generate a reverse static single assignment graph based on the dataflow graph instructions, the reverse static single assignment graph comprising strongly connected components, each of the strongly connected components associated with at least one of the one or more pick/switch instruction pairs. Embodiments also include traversing the reverse static single assignment graph depth-first, and replace pick/switch instructions associated with strongly connected components having configuration values with compound instructions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Numerical analysis of base flowfield at high altitude for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a pressure based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution and turbulence model are two primary factors that influence the accuracy of the base flowfield prediction.

Wang, Ten-See↗

DOE OSTI2021

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗