Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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 OSTI2017

On the use of reverse Brownian motion to accelerate hybrid simulations

Multiscale and multiphysics simulations are two rapidly developing fields of scientific computing. Efficient coupling of continuum (deterministic or stochastic) constitutive solvers with their discrete (stochastic, particle-based) counterparts is a common challenge in both kinds of simulations. We focus on interfacial, tightly coupled simulations of diffusion that combine continuum and particle-based solvers. The latter employs the reverse Brownian motion (rBm), a Monte Carlo approach that allows one to enforce inhomogeneous Dirichlet, Neumann, or Robin boundary conditions and is trivially parallelizable. We discuss numerical approaches for improving the accuracy of rBm in the presence of inhomogeneous Neumann boundary conditions and alternative strategies for coupling the rBm solver with its continuum counterpart. Numerical experiments are used to investigate the convergence, stability, and computational efficiency of the proposed hybrid algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental validation benchmark data for CFD of transient convection from forced to natural with flow reversal on a vertical flat plate

Transient convection has been investigated experimentally for the purpose of providing Computational Fluid Dynamics (CFD) validation benchmark data. A specialized facility for validation benchmark experiments called the Rotatable Buoyancy Tunnel was used to acquire thermal and velocity measurements of flow over a smooth, vertical heated plate. The initial condition was forced convection downward with subsequent transition to mixed convection, ending with natural convection upward after a flow reversal. Data acquisition through the transient was repeated for ensemble-averaged results. With simple flow geometry, validation data were acquired at the benchmark level. All boundary conditions (BCs) were measured and their uncertainties quantified. Temperature profiles on all four walls and the inlet were measured, as well as as-built test section geometry. Inlet velocity profiles and turbulence levels were quantified using Particle Image Velocimetry. System Response Quantities (SRQs) were measured for comparison with CFD outputs and include velocity profiles, wall heat flux, and wall shear stress. Extra effort was invested in documenting and preserving the validation data. Details about the experimental facility, instrumentation, experimental procedure, materials, BCs, and SRQs are made available through this paper. As a result, the latter two are available for download and the other details are included in this work.

42 ENGINEERING↗
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 OSTI2015

Equation-of-motion coupled-cluster method for doubly ionized states with spin-orbit coupling

In this work, we report implementation of the equation-of-motion coupled-cluster method for doubly ionized states (EOM-DIP-CC) with spin-orbit coupling (SOC) using a closed-shell reference. Double ionization potentials (DIPs) are calculated in the space spanned by 2h and 3h1p determinants with the EOM-DIP-CC approach at the CC singles and doubles level (CCSD). Time-reversal symmetry together with spatial symmetry is exploited to reduce computational effort. To circumvent the problem of unstable dianion references when diffuse basis functions are included, nuclear charges are scaled. Effect of this stabilization potential on DIPs is estimated based on results from calculations using a small basis set without diffuse basis functions. DIPs and excitation energies of some low-lying states for a series of open-shell atoms and molecules containing heavy elements with two unpaired electrons have been calculated with the EOM-DIP-CCSD approach. Results show that this approach is able to afford a reliable description on SOC splitting. Furthermore, the EOM-DIP-CCSD approach is shown to provide reasonable excitation energies for systems with a dianion reference when diffuse basis functions are not employed.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2018

Statistical relationship between metabolic decomposition and chemical uptake predicts bioconcentration factor data for diverse chemical exposures

Background: A challenge of in vitro to in vivo extrapolation (IVIVE) is to predict the physical state of organisms exposed to chemicals in the environment from in vitro exposure assay data. Although toxicokinetic modeling approaches promise to bridge in vitro screening data with i Danio rerio ) based upon particle swarm optimizations (PSO) of the chemical uptake and degradation rates that predict bioconcentration factors (BCF) for a broad range of chemicals. PSO reveals a relationship between chemical uptake and decomposition parameter values that predicts chemical-specific BCF values with moderate statistical agreement to a limited yet diverse chemical dataset, and all without a need to retrain the model to new data. Conclusions: The presented model requires only the octanol-water partitioning ratio to predict BCFs to a fidelity consistent with existing QSAR models. This success begs re-evaluation of the modeling assumptions; specifically, it suggests that chemical uptake into arterial blood may be limited by transport across gill membranes (diffusion) rather than by counter-current flow between gill lamellae (convection). Therefore, more detailed molecular modeling of aquatic respiration may further improve predictive accuracy of the rTK approach.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2013

Evaluation of the robustness of the preprocessing technique improving reversible compressibility of CT images: Tested on various CT examinations

Purpose: To modify the preprocessing technique, which was previously proposed, improving compressibility of computed tomography (CT) images to cover the diversity of three dimensional configurations of different body parts and to evaluate the robustness of the technique in terms of segmentation correctness and increase in reversible compression ratio (CR) for various CT examinations.Methods: This study had institutional review board approval with waiver of informed patient consent. A preprocessing technique was previously proposed to improve the compressibility of CT images by replacing pixel values outside the body region with a constant value resulting in maximizing data redundancy. Since the technique was developed aiming at only chest CT images, the authors modified the segmentation method to cover the diversity of three dimensional configurations of different body parts. The modified version was evaluated as follows. In randomly selected 368 CT examinations (352 787 images), each image was preprocessed by using the modified preprocessing technique. Radiologists visually confirmed whether the segmented region covers the body region or not. The images with and without the preprocessing were reversibly compressed using Joint Photographic Experts Group (JPEG), JPEG2000 two-dimensional (2D), and JPEG2000 three-dimensional (3D) compressions. The percentage increase in CR per examination (CR{sub I}) was measured.Results: The rate of correct segmentation was 100.0% (95% CI: 99.9%, 100.0%) for all the examinations. The median of CR{sub I} were 26.1% (95% CI: 24.9%, 27.1%), 40.2% (38.5%, 41.1%), and 34.5% (32.7%, 36.2%) in JPEG, JPEG2000 2D, and JPEG2000 3D, respectively.Conclusions: In various CT examinations, the modified preprocessing technique can increase in the CR by 25% or more without concerning about degradation of diagnostic information.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

First-order finite-Larmor-radius fluid modeling of tearing and relaxation in a plasma pinch

Drift and Hall effects on magnetic tearing, island evolution, and relaxation in pinch configurations are investigated using a non-reduced first-order finite-Larmor-radius (FLR) fluid model with the nonideal magnetohydrodynamics (MHD) with rotation, open discussion (NIMROD) code [C.R. Sovinec and J. R. King, J. Comput. Phys. 229, 5803 (2010)]. An unexpected result with a uniform pressure profile is a drift effect that reduces the growth rate when the ion sound gyroradius ({rho}{sub s}) is smaller than the tearing-layer width. This drift is present only with warm-ion FLR modeling, and analytics show that it arises from {nabla}B and poloidal curvature represented in the Braginskii gyroviscous stress. Nonlinear single-helicity computations with experimentally relevant {rho}{sub s} values show that the warm-ion gyroviscous effects reduce saturated-island widths. Computations with multiple nonlinearly interacting tearing fluctuations find that m = 1 core-resonant-fluctuation amplitudes are reduced by a factor of two relative to single-fluid modeling by the warm-ion effects. These reduced core-resonant-fluctuation amplitudes compare favorably to edge coil measurements in the Madison Symmetric Torus (MST) reversed-field pinch [R. N. Dexter et al., Fusion Technol. 19, 131 (1991)]. The computations demonstrate that fluctuations induce both MHD- and Hall-dynamo emfs during relaxation events. The presence of a Hall-dynamo emf implies a fluctuation-induced Maxwell stress, and the simulation results show net transport of parallel momentum. The computed magnitude of force densities from the Maxwell and competing Reynolds stresses, and changes in the parallel flow profile, are qualitatively and semi-quantitatively similar to measurements during relaxation in MST.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗