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 451 records · Page 25

Strongly tunable anisotropic thermal transport in MoS 2 by strain and lithium intercalation: first-principles calculations

The possibility of tuning the vibrational properties and the thermal conductivity of layered van der Waals materials either chemically or mechanically paves the way to significant advances in nanoscale heat management. Using first-principles calculations we investigate the modulation of heat transport in MoS 2 by lithium intercalation and cross-plane strain. Here, we find that both the in-plane and cross-plane thermal conductivity (κ τ ,κ z ) of MoS 2 are extremely sensitive to both strain and electrochemical intercalation. Combining lithium intercalation and strain, the in-plane and cross-plane thermal conductivity can be tuned over one and two orders of magnitude, respectively. Furthermore, since and respond in different ways to intercalation and strain, the thermal conductivity anisotropy can be modulated by two orders of magnitude. The underlying mechanisms for such large tunability of the anisotropic thermal conductivity of MoS 2 are explored by computing and analyzing the dispersion relations, group velocities, relaxation times and mean free paths of phonons. Since both intercalation and strain can be applied reversibly, their stark effect on thermal conductivity can be exploited to design novel phononic devices, as well as for thermal management in MoS 2 -based electronic and optoelectronic systems.

36 MATERIALS SCIENCE↗

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗

Flow-Field Measurements in a Wing-Fuselage Junction Using an Embedded Particle Image Velocimetry System

A novel, embedded, 2D Particle Image Velocimetry system has been developed and implemented to obtain off-body velocity measurements in the junction of an 8% wing-body configuration in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel. Although the initial purpose for implementing the system during this test was to evaluate system performance and identify potential risks, the system worked remarkably well and a considerable amount of data were obtained in the wing-fuselage junction near the trailing edge at a Reynolds Number of 2.4 Million and angle of attack of 5 degrees. In addition to providing notable efficiencies with regard to image acquisition and test operations, the PIV system proved to be quite capable of capturing unique details of the flow separation to complement the extensive suite of measurement techniques applied during the test. Instantaneous PIV vector fields reveal that the flow separation is not stationary but rather highly dynamic. Mean flow statistics calculated from the PIV measurements highlight where reverse flow and Reynolds stresses are concentrated in the separated region and agree well with results from the embedded Laser Doppler Velocimeter system and Computational Fluid Dynamics. The comparisons and additional insight gained during this effort will help guide system improvements for follow-on tests.

Jenkins, Luther N.↗

Update on the Juncture Flow PIV Results and Future Plans

A novel, embedded, 2D Particle Image Velocimetry system has been developed and implemented to obtain off-body velocity measurements in the junction of an 8% wing-body configuration in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel. Although the initial purpose for implementing the system during this test was to evaluate system performance and identify potential risks, a considerable amount of data were obtained in the wing-fuselage junction near the trailing edge at a Reynolds Number of 2.4 Million and angle of attack of 5 degrees. In addition to providing notable efficiencies with regard to image acquisition and test operations, the PIV system captured unique details of the flow separation to complement the extensive suite of measurement techniques applied during the test. Instantaneous PIV vector fields reveal that the flow separation is not stationary but rather highly dynamic. Mean flow statistics calculated from the PIV measurements highlight where reverse flow and Reynolds stresses are concentrated in the separated region and agree well with results from the embedded Laser Doppler Velocimeter system and Computational Fluid Dynamics. The comparisons and ional insight gained during this effort will help guide system improvements for the upcoming entry in 2020.

Jenkins, Luther N.↗

An Efficient GPU-Accelerated Multi-Source Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog ∼12000 Galactic Binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA global fit↗

Reverse phase transitions in binary neutron-star systems with exotic-matter cores

Multimessenger observations of binary neutron-star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter—e.g., quark matter—the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron-star mergers, or if they appear at even higher densities that are not realized in the cosmos. Here, in this article, we perform numerical relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiraling binary component. Our simulations reveal that shortly before the merger, it is possible that such stars undergo a “reverse phase transition”—i.e., densities decrease and the quark core inside the star disappears, leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to the rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the postmerger phase transition—e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition.

79 ASTRONOMY AND ASTROPHYSICS↗

Swirling jet turbulent mixing and combustion computations

Computations are presented describing the mixing and combustion of swirling jets in a coaxial stream. It is demonstrated that the boundary layer equations represent the flow reasonably well until reversed flow is imminent. For the range of parameters investigated indications are that the edge velocity has little effect on the behavior of the flow. Furthermore, confining the flow with a constant pressure wall, or impressing a favorable pressure gradient on the coaxial flow, acts to reduce the severity of the centerline adverse pressure gradient created by the swirl decay. A simple scalar eddy viscosity model, including a potential core formulation, is shown to described the behavior of weak swirling flow in the far region but is only in fair agreement with observations in the near region. The effects of swirl on a burning hydrocarbon jet exhausting into a cold coaxial stream are shown to be intensified by the reduction of the density due to combustion. The enhanced mixing properties of high swirl flow produce rapid diffusion of the burning gases into the cold edge flow causing early cessation of the NO producing reactions. Computations show that doubling the initial jet swirl could reduce the NO production by 25 percent.

Rubel, A.↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Multistates and Polyamorphism in Phase-Change K 2 Sb 8 Se 13

The phase-change (PC) materials in the majority of optical data storage media in use today exhibit a fast, reversible crystal -> amorphous phase transition that allows them to be switched between on (1) and off (0) binary states. Solid-state inorganic materials with this property are relatively common, but those exhibiting an amorphous -> amorphous transition called polyamorphism are exceptionally rare. K 2 Sb 8 Se 13 (KSS) reported here is the first example of a material that has both amorphous -> amorphous polyamorphic transition and amorphous -> crystal transition at easily accessible temperatures (227 and 263 degrees C, respectively). The transitions are associated with the atomic coordinative preferences of the atoms, and all three states of K 2 Sb 8 Se 13 are stable in air at 25 degrees C and 1 atm. All three states of K 2 Sb 8 Se 13 exhibit distinct optical bandgaps, E-g = 1.25, 1.0, and 0.74 eV, for the amorphous-II, amorphous-I, and crystalline versions, respectively. The room-temperature electrical conductivity increases by more than 2 orders of magnitude from amorphous-I to -II and by another 2 orders of magnitude from amorphous-II to the crystalline state. This extraordinary behavior suggests that a new class of materials exist which could provide multistate level systems to enable higher-order computing logic circuits, reconfigurable logic devices, and optical switches.

amorphous-amorphous transition↗

HBMax: Optimizing Memory Efficiency for Parallel Influence Maximization on Multicore Architectures

The goal of influence maximization is to select k most-influential vertices or seeds in a network, where influence is defined by a given diffusion process. The problem has a number of important applications such as viral marketing, information spread, and epidemic control. Although computing optimal seed set is NP-Hard, due to the submodular nature of the problem efficient approximation algorithms exist. However, even state-of-the-art parallel implementations are limited by a sampling step that incurs large memory footprints. This in turn limits the problem size reach and approximation quality. In this work, we study the memory footprint of the sampling process collecting reverse reachability information in the IMM algorithm over large real-world social networks. We present an adaptive and memory-efficient optimization approach for a state-of-the-art multi-threaded parallel influence maximization algorithm. Our approach,HuffMax, uses a portion of the reverse reachable (RR) sets collected by the algorithm to learn the characteristics of the graph. Then, it compresses the intermediate reverse reachability information with Huffman coding, and queries directly on the compressed data to preserve the memory savings obtained through compression. We also propose an efficient sampling strategy based on the distribution of RR sets, which can further reduce the computation time for typical social networks with long-tail distributions. Considering a NUMA architecture, we scale up our solution on 128-core CPUs and reduce the memory footprint by up to 45.7% with negligible time overhead (or even faster) and without perceivable loss of accuracy.

Chen, Xinyu↗

Automatic Differentiation for Adjoint Stencil Loops

Stencil loops are a common motif in computations including convolutional neural networks, structured-mesh solvers for partial differential equations, and image processing. Stencil loops are easy to parallelise, and their fast execution is aided by compilers, libraries, and domain-specific languages. Reverse-mode automatic differentiation, also known as algorithmic differentiation, autodiff, adjoint differentiation, or back-propagation, is sometimes used to obtain gradients of programs that contain stencil loops. Unfortunately, conventional automatic differentiation results in a memory access pattern that is not stencil-like and not easily parallelisable. In this paper we present a novel combination of automatic differentiation and loop transformations that preserves the structure and memory access pattern of stencil loops, while computing fully consistent derivatives. The generated loops can be parallelised and optimised for performance in the same way and using the same tools as the original computation. We have implemented this new technique in the Python tool PerforAD, which we release with this paper along with test cases derived from seismic imaging and computational fluid dynamics applications.

Automatic Differentiation↗

In-situ effect in cross-ply laminates under various loading conditions analyzed with hybrid macro/micro-scale computational models

In this article, multi-scale finite element analyses based on three-dimensional (3D) hybrid macro/micro-scale computational models subjected to various loading conditions are carried out to examine the in-situ effect imposed by the neighboring plies on the failure initiation and propagation of cross-ply laminates. A detailed comparative study on crack suppression mechanisms due to the effect of embedded laminar thickness and adjacent ply orientation is presented. Furthermore, we compare the results of in-situ transverse failure strain and strength between the computational models and analytical predictions. Good agreements are generally observed, indicating the constructed computational models are highly accurate to quantify the in-situ effect. Subsequently, empirical formulas for calculating the in-situ strengths as a function of embedded ply thickness and different ply angle between embedded and adjacent plies are developed, during which several material parameters are obtained using a reverse fitting method. Finally, a new set of failure criteria for σ 22 -τ 12 , σ 22 -τ 23 , and σ 11 -τ 12 accounting for the in-situ strengths are proposed to predict laminated composites failure under multi-axial stress states. This study demonstrates an effective and efficient computational technique towards the accurate prediction of the failure behaviors and strengths of cross-ply laminates by including the in-situ effects.

36 MATERIALS SCIENCE↗

Observable measurement-induced transitions

One of the main postulates of quantum mechanics is that measurements destroy quantum coherence (wave function collapse). Recently it was discovered that in a many-body system dilute local measurements still preserve some coherence across the entire system. As the measurement density is increased, a phase transition occurs that is characterized by the disentanglement of different parts of the system. Unfortunately, this transition is impossible to observe experimentally for macroscopic systems because it requires an exponentially costly full tomography of the many-body wave function or a comparison with the simulation on an oracle classical computer. In this work we report the discovery of another measurement-induced phase transition that can be observed experimentally if quantum dynamics can be reversed. On one side of this phase transition the quantum information encoded in some part of the Hilbert space is fully recovered after the time inversion. On the other side, all quantum information is corrupted. This transition also manifests itself as the change in the behavior of the probability to observe the same measurement outcome in the process that consists of identical blocks repeated many times. In each block the unitary evolution is followed by the measurement. On one side of the transition the probability decreases exponentially with the number of repetitions, on the other it tends to a constant as the number of repetitions is increased. We confirm the existence of the proposed phase transition through numerical simulations of realistic quantum circuits and analytical calculations using an effective random-matrix theory model.

Measurement-induced phase transitions↗

Excitations in the field-induced quantum spin liquid state of α-RuCl 3

The celebrated Kitaev quantum spin liquid (QSL) is the paradigmatic example of a topological magnet with emergent excitations in the form of Majorana Fermions and gauge fluxes. Upon breaking of time-reversal symmetry, for example in an external magnetic field, these fractionalized quasiparticles acquire non-Abelian exchange statistics, an important ingredient for topologically protected quantum computing. Consequently, there has been enormous interest in exploring possible material realizations of Kitaev physics and several candidate materials have been put forward, recently including α-RuCl 3 . In the absence of a magnetic field this material orders at a finite temperature and exhibits low-energy spin wave excitations. However, at moderate energies, the spectrum is unconventional and the response shows evidence for fractional excitations. Here in this paper, we use time-of-flight inelastic neutron scattering to show that the application of a sufficiently large magnetic field in the honeycomb plane suppresses the magnetic order and the spin waves, leaving a gapped continuum spectrum of magnetic excitations. Our comparisons of the scattering to the available calculations for a Kitaev QSL show that they are consistent with the magnetic field induced QSL phase.

36 MATERIALS SCIENCE↗

Ultrasonic Doppler measurement of renal artery blood flow

Studies were made of (1) blood flow redistribution during lower body negative pressure (LBNP), (2) the profile of blood flow across the mitral annulus of the heart (both perpendicular and parallel to the commissures), (3) testing and evaluation of a number of pulsed Doppler systems, (4) acute calibration of perivascular Doppler transducers, (5) redesign of the mitral flow transducers to improve reliability and ease of construction, and (6) a frequency offset generator designed for use in distinguishing forward and reverse components of blood flow by producing frequencies above and below the offset frequency. Finally methodology was developed and initial results were obtained from a computer analysis of time-varying Doppler spectra.

Freund, W. R.↗

Overlapping emission peaks in the solar C I multiplets at lambda 1560 and lambda 1657

Observations of the C I multiplets at 1560 and 1657 A made with the University of Colorado spectrometer on the OSO 8 satellite are presented and compared with computed profiles for the Vernazza-Avrett-Loeser solar atmosphere. These are optically thick emission lines formed in the solar chromosphere that show the central reversals typical of such lines. In each multiplet there is an interesting case of overlapping emission peaks which shows that such peaks do not constructively combine but instead weaken. This behavior is easily understood and reproduced with an optically thick, non-LTE mode of formation for these lines and is not consistent with an optically thin mechanism. We also find that the shapes of these blends are very sensitive to the magnitude of the nonthermal microvelocities.

Shine, R. A.↗

Observation of gravity waves during the extreme tornado outbreak of 3 April 1974

A continuous wave-spectrum high-frequency radiowave Doppler sounder array was used to observe upper-atmospheric disturbances during an extreme tornado outbreak. The observations indicated that gravity waves with two harmonic wave periods were detected at the F-region ionospheric height. Using a group ray path computational technique, the observed gravity waves were traced in order to locate potential sources. The signals were apparently excited 1-3 hours before tornado touchdown. Reverse ray tracing indicated that the wave source was located at the aurora zone with a Kp index of 6 at the time of wave excitation. The summation of the 24-hour Kp index for the day was 36. The results agree with existing theories (Testud, 1970; Titheridge, 1971; Kato, 1976) for the excitation of large-scale traveling ionospheric disturbances associated with geomagnetic activity in the aurora zone.

Hung, R. J.↗

Biomimetic robots using EAP as artificial muscles - progress and challenges

Biology offers a great model for emulation in areas ranging from tools, computational algorithms, materials science, mechanisms and information technology. In recent years, the field of biomimetics, namely mimicking biology, has blossomed with significant advances enabling the reverse engineering of many animals' functions and implementation of some of these capabilities.

robotics↗