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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
DOE OSTI2023

Dynamic modeling and simulation of pressure swing adsorption processes using toPSAil

Pressure swing adsorption (PSA) has attracted significant recent interest for chemical process intensification due to its potential for high energy efficiency and amenability to small, modular designs. However, the lack of simulation tools that are readily available, transparent, and trusted, has been identified as a serious impediment to widespread adoption of PSA, as well as to further research on PSA modeling, numerical solution, optimization, and control. This paper presents a complete framework for dynamic modeling and simulation of PSA processes and its implementation in an open-source simulator called toPSAil. Further, the presentation is tutorial and includes many modeling and implementation details often overlooked in existing literature. Novel methods are presented for handling flow reversals and implementing various pressure–flow relationships, along with controlled boundary conditions. Finally, the code contains several innovations designed to improve efficiency and reduce the extensive trial-and-error tuning often required to produce a working PSA cycle.

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

Correlated nanoelectronics and the second quantum revolution

The growing field of correlated nanoelectronics exists at the intersection of two established fields: correlated oxide electronics and semiconductor nanoelectronics. The development of quantum technologies that exploit quantum coherence and entanglement for the purposes of computation, simulation, and sensing will require complex material properties to be controlled at nanoscale dimensions. Heterostructures and nanostructures formed at the interface between LaAlO 3 and SrTiO 3 exhibit striking behavior that arises from the ability to program the conductive behavior at extreme nanoscale dimensions. The active electronic layer, SrTiO 3 , exhibits a wide range of gate-tunable phenomena such as ferroelectricity, ferroelasticity, magnetism, superconductivity, and spin–orbit coupling, all of which can be controlled at the nanoscale using two reversible methods: conductive atomic force microscope lithography and ultra-low-voltage electron beam lithography. Mesoscopic devices such as single-electron transistors and quasi-one-dimensional electron waveguides can be “sketched” using these techniques, and the properties of these devices differ significantly from those created from traditional semiconductors, such as Si or GaAs. The strongly correlated nature of the SrTiO 3 system is evident from superconducting behavior as well as a state in which electrons are paired outside the superconducting state. A highly exotic phase was discovered in which a degenerate quantum liquid is formed from bound states of n = 2, 3, 4, … electrons. Further development of correlated nanoelectronics based on the LaAlO 3 /SrTiO 3 system can potentially lead to a general platform for quantum simulation as well as a pathway for the development of highly entangled states of multiple photons.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Numerical study of the performance of swept, curved compression surface scramjet inlets

The paper describes a computational performance-enhancement study in which systematic modifications were applied to a planar-sidewall compression scramjet inlet operating at an entrance Mach number of 4 and at a dynamic pressure of 2040 lbs/sq ft, and the performance of the modified inlet was compared to that of the original geometry by evaluating numerical flowfield solutions. The planar-sidewall compression angle was modified as a function of height using sidewall curvature and simultaneously employing both forward-swept and reverse-swept compression surfaces. The results of calculations illustrate the feasibility of obtaining enhanced performance via contour modifications of a baseline swept-sidewall compression scramjet inlet operating at an entrance Mach number of 4.

Korte, John J.↗
NASA NTRS1995

Material processing of convection-driven flow field and temperature distribution under oblique gravity

A set of mathematical formulation is adopted to study vapor deposition from source materials driven by heat transfer process under normal and oblique directions of gravitational acceleration with extremely low pressure environment of 10(exp -2) mm Hg. A series of time animation of the initiation and development of flow and temperature profiles during the course of vapor deposition has been obtained through the numerical computation. Computations show that the process of vapor deposition has been accomplished by the transfer of vapor through a fairly complicated flow pattern of recirculation under normal direction gravitational acceleration. It is obvious that there is no way to produce a homogeneous thin crystalline films with fine grains under such a complicated flow pattern of recirculation with a non-uniform temperature distribution under normal direction gravitational acceleration. There is no vapor deposition due to a stably stratified medium without convection for reverse normal direction gravitational acceleration. Vapor deposition under oblique direction gravitational acceleration introduces a reduced gravitational acceleration in vertical direction which is favorable to produce a homogeneous thin crystalline films. However, oblique direction gravitational acceleration also induces an unfavorable gravitational acceleration along horizontal direction which is responsible to initiate a complicated flow pattern of recirculation. In other words, it is necessary to carry out vapor deposition under a reduced gravity in the future space shuttle experiments with extremely low pressure environment to process vapor deposition with a homogeneous crystalline films with fine grains. Fluid mechanics simulation can be used as a tool to suggest most optimistic way of experiment with best setup to achieve the goal of processing best nonlinear optical materials.

Hung, R. J.↗
DOE OSTI2022

Toward Discord : Code for Simulating Continuous Spin Systems

A new computational tool to simulate classical spin systems with frustrated crystal structures is presented. Complementary single- and cluster-spin flip algorithms are implemented to calculate the diffuse scattering patterns, spin-pair correlations, and thermodynamic quantities. Test cases of geometrically frustrated kagome, pyrochlore, and cubic systems are detailed. Two recent scientific cases are also shown here. This new method, together with recent developments of the rmc-discord package (https://github.com/zjmorgan/rmc-discord), represent integrated and strategic step in a complete forward and reverse Monte Carlo framework discord.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Catalytic N–H Bond Formation Promoted by a Ruthenium Hydride Complex Bearing a Redox-Active Pyrimidine-Imine Ligand

The synthesis of a piano-stool ruthenium hydride, [(η 5 -C 5 Me 5 )Ru(PmIm)H] (PmIm = (N-(1,3,5-trimethylphenyl)-1-(pyrimidin-2-yl)ethan-1-imine) for the dual purpose of catalytic dihydrogen activation and subsequent hydrogen atom transfer for the formation of weak chemical bonds is described. Introduction of a neutral, potentially redox-active PmIm supporting ligand was designed to eliminate the possibility of deleterious C(sp 2 )–H reductive coupling and elimination that has been identified as a deactivation pathway with related rhodium and iridium catalysts. Treatment of [(η 5 -C 5 Me 5 )RuCl 2 ] n with one equivalent PmIm ligand in the presence of zinc and sodium methoxide resulted in isolation of the diruthenium complex, [(η 5 -C 5 Me 5 )Ru(PmIm)] 2 , arising from C–C bond formation between two PmIm chelates. Addition of H 2 to the ruthenium dimer under both thermal and blue light irradiation conditions furnished the targeted hydride, [(η 5 -C 5 Me 5 )Ru(PmIm)H] which has a relatively weak DFT-calculated Ru–H bond dissociation free energy (BDFE) of 47.9 kcal/mol. Addition of TEMPO to [(η 5 -C 5 Me 5 )Ru(PmIm)H] generated the 17-electron metalloradical, [(η 5 -C 5 Me 5 )Ru(PmIm)], which was characterized by EPR spectroscopy. The C–C bond forming process was reversible as irradiation of [(η 5 -C 5 Me 5 )Ru(PmIm)] 2 generated [(η 5 -C 5 Me 5 )Ru(PmIm)H] and a piano-stool ruthenium complex containing an enamide ligand derived from H-atom ab-straction from the PmIm chelate. Equilibration studies were used to establish an experimental estimate of the effective Ru–H BDFE and a value of 50.8 kcal/mol was obtained, in agreement with the observed loss of H 2 and the DFT-computed value. Here, the ruthenium hydride was an effective catalyst for the thermal catalytic hydrogenation of TEMPO, acridine and a cobalt-imido complex and for the selective reduction of azobenzene to diphenylhydrazine using H 2 as the terminal reductant, highlighting the role of this complex is catalytic weak bond formation using H 2 as the stoichiometric reductant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

A method for determining local elastoplastic stress and strain in metallurgically bonded notched laminates subjected to a loading cycle

A semianalytical method was developed for determining elastoplastic cyclic stresses and strains at notch roots in metallurgically bonded metal laminates. The method is based on the Neuber equation, which was used with an effective stress-strain curve for the laminate. It was applied to laminates containing a circular hole which were subjected to one cycle of reversed loading. The laminates consisted of two elasto-perfectly-plastic materials with different yield strengths and with either equal or different Young's moduli. A laminate of high-strength titanium alloy with alternate layers of commercially pure titanium was also analyzed. The accuracy of the method was evaluated by comparing the stresses and strains with those calculated from a finite-element analysis. The results estimated by the simple method based on the Neuber equation agreed closely with the results computed from the more elaborate finite-element analysis.

Sova, J. A.↗
NASA NTRS1995

Kinetics of cycle length dependence of ventricular repolarization: effect of autonomic blockade

INTRODUCTION: Beat-to-beat adaptation of ventricular repolarization duration to cardiac cycle length and autonomic activity has not been previously characterized in the spontaneously beating human heart. METHODS AND RESULTS: The ECG of 14 healthy subjects was recorded from the supine and upright positions. Autonomic blockade was accomplished by atropine and propranolol. RR and RT intervals were measured by a computer algorithm, and the impulse response (h) from RR to RT computed. In the supine position the maximal adjustment of the RT interval occurred in the first beat following a change in cycle length (hpeak = 17.8 +/- 1.6 msec/sec), but continued to be detectable for 3.8 seconds (2.9-4.7 sec). Propranolol attenuated the peak impulse response to 15.8 +/- 4.0 msec/sec (P = NS). In the standing position the peak impulse response was increased to 25.2 +/- 5.0 msec/sec (P = 0.004 vs supine), and the impulse response duration (hdur) shortened to 1.4 seconds (1.3-1.6). This was reversed by beta blockade (hpeak = 10.7 +/- 3.6 [P = 0.005 vs standing]; hdur = 5.5 sec [4.8-6.1]). Parasympathetic and combined autonomic blockade resulted in too little residual heart rate variability to estimate the impulse response accurately. The slope of the regression of delta RT and delta RR in the supine position was 0.0177 +/- 0.0016, which was closely correlated with the peak impulse response (r = 0.91). CONCLUSIONS: System identification techniques can assist in characterizing the cycle dependence of ventricular repolarization and may provide new insights into conditions associated with abnormal repolarization.

Clinical Trial↗
NASA NTRS1991

Computer-based video digitizer analysis of surface extension in maize roots: kinetics of growth rate changes during gravitropism

We used a video digitizer system to measure surface extension and curvature in gravistimulated primary roots of maize (Zea mays L.). Downward curvature began about 25 +/- 7 min after gravistimulation and resulted from a combination of enhanced growth along the upper surface and reduced growth along the lower surface relative to growth in vertically oriented controls. The roots curved at a rate of 1.4 +/- 0.5 degrees min-1 but the pattern of curvature varied somewhat. In about 35% of the samples the roots curved steadily downward and the rate of curvature slowed as the root neared 90 degrees. A final angle of about 90 degrees was reached 110 +/- 35 min after the start of gravistimulation. In about 65% of the samples there was a period of backward curvature (partial reversal of curvature) during the response. In some cases (about 15% of those showing a period of reverse bending) this period of backward curvature occurred before the root reached 90 degrees. Following transient backward curvature, downward curvature resumed and the root approached a final angle of about 90 degrees. In about 65% of the roots showing a period of reverse curvature, the roots curved steadily past the vertical, reaching maximum curvature about 205 +/- 65 min after gravistimulation. The direction of curvature then reversed back toward the vertical. After one or two oscillations about the vertical the roots obtained a vertical orientation and the distribution of growth within the root tip became the same as that prior to gravistimulation. The period of transient backward curvature coincided with and was evidently caused by enhancement of growth along the concave and inhibition of growth along the convex side of the curve, a pattern opposite to that prevailing in the earlier stages of downward curvature. There were periods during the gravitropic response when the normally unimodal growth-rate distribution within the elongation zone became bimodal with two peaks of rapid elongation separated by a region of reduced elongation rate. This occurred at different times on the convex and concave sides of the graviresponding root. During the period of steady downward curvature the elongation zone along the convex side extended farther toward the tip than in the vertical control. During the period of reduced rate of curvature, the zone of elongation extended farther toward the tip along the concave side of the root. The data show that the gravitropic response pattern varies with time and involves changes in localized elongation rates as well as changes in the length and position of the elongation zone. Models of root gravitropic curvature based on simple unimodal inhibition of growth along the lower side cannot account for these complex growth patterns.

NASA Discipline Plant Biology↗
NASA NTRS1979

Gaseous NH4Cl revisited - A computational investigation of the potential surface and properties

An investigation of the potential surface and properties of the ground state of gaseous NH4Cl is presented. The calculations, which utilize a good basis set and include investigation of the valence correlation effects, result in a description of gaseous NH4Cl as a weakly bound complex with NH3 and HCl structures essentially equivalent to those of the isolated fragments. The charge distributions within the fragments are modifed in the complex, however, and there is a small amount of charge transfer. The results are in reasonable agreement with estimate of D sub 0, super 0 from high temperature mass spectroscopic investigations and vibrational frequencies from matrix experiments. The equilibrium constant for the reversible reaction NH4Cl yields NH3 + HCl indicates that the homogeneous formation of NH4Cl does not represent a significant chemical sink for HCl in the upper atmosphere

Raffenetti, R. C.↗
DOE OSTI2023

Reversible Cl/Cl - redox in a spinel Mn 3 O 4 electrode

A unique prospect of using halides as charge carriers is the possibility of the halides undergoing anodic redox behaviors when serving as charge carriers for the charge-neutrality compensation of electrodes. However, the anodic conversion of halides to neutral halogen species has often been irreversible at room temperature due to the emergence of diatomic halogen gaseous products. Here, we report that chloride ions can be reversibly converted to near-neutral atomic chlorine species in the Mn 3 O 4 electrode at room temperature in a highly concentrated chloride-based aqueous electrolyte. Notably, the Zn 2+ cations inserted in the first discharge and trapped in the Mn 3 O 4 structure create an environment to stabilize the converted chlorine atoms within the structure. Characterization results suggest that the Cl/Cl - redox is responsible for the observed large capacity, as the oxidation state of Mn barely changes upon charging. Computation results corroborate that the converted chlorine species exist as polychloride monoanions, e.g., [Cl 3 ] - and [Cl 5 ] - , inside the Zn 2+ -trapped Mn 3 O 4 , and the presence of polychloride species is confirmed experimentally. Our results point to the halogen plating inside electrode lattices as a new charge-storage mechanism.

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

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Investigation of supersonic separated flow in a compression corner by laser Doppler anemometry

A single-component, computer-operated, three dimensional traversing laser Doppler velocimetry system was designed and constructed for a supersonic wind tunnel. The model was a 10 deg compression corner, providing an example of laminar boundary layer separation. Static pressure data and color schlieren photographs were taken. The Mach number was 2.42 and the Reynolds number was 213,000. The flow was seeded with submicron sized oil droplets. The trend of decreasing upstream influence with increasing Reynolds number was confirmed. Mean velocity profiles provided experimental evidence of reversed flow. Points of separation and reattachment were consistent with those determined by schlieren and pressure scans. The flow was of the laminar type until downstream of reattachment. Individual mean velocity and turbulence profiles, as well as velocity histograms are presented.

Baroth, E. C.↗
NASA NTRS1985

The closure problem for turbulence in meteorology and oceanography

The dependent variables used for computer based meteorological predictions and in plans for oceanographic predictions are wave number and frequency filtered values that retain only scales resolvable by the model. Scales unresolvable by the grid in use become 'turbulence'. Whether or not properly processed data are used for initial values is important, especially for sparce data. Fickian diffusion with a constant eddy diffusion is used as a closure for many of the present models. A physically realistic closure based on more modern turbulence concepts, especially one with a reverse cascade at the right times and places, could help improve predictions.

Pierson, W. J., Jr.↗
NASA NTRS1982

Nacelle Design

Rohr Industries, Inc. specializes in the manufacturing of nacelles, thrust reversers, and other engine components, has made use of the COSMIC program in nacelle work for aircraft like the McDonnell Douglas DC-10, the French/German Airbus and the Boeing 727, 737 and 747 jetliners. They also manufacture complete nacelles for military and business aircraft and are supplying nacelle components for the Boeing 757. The nacelle requires research and development for each type of airplane because of complex airflows around the engine inlet and high pressures on the nacelle skin. The use of the computer program defines the airflow field around turbofan engine nacelle inlet and cowls. Pressures on the nacelle skin are estimated for various flight conditions and structural integrity of the nacelle design. Rohr estimates that use of this program saved six man-months of programmer time necessary to develop alternative software.

Source record↗
NASA NTRS2010

The Writhe of Helical Structures in the Solar Corona

Context. Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when they erupt. Aims. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. Methods. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). Results. (1) We demonstrate that the relation between writhe and projected S shape is not unique in principle, but that the ambiguity does not affect low-lying structures, thus supporting the established empirical rule which associates stable forward (reverse) S shaped structures low in the corona with positive (negative) helicity. (2) Kink-unstable erupting flux ropes are found to transform a far smaller fraction of their twist helicity into writhe helicity than often assumed. (3) Confined flux rope eruptions tend to show stronger writhe at low heights than ejective eruptions (CMEs). This argues against suggestions that the writhing facilitates the rise of the rope through the overlying field. (4) Erupting filaments which are S shaped already before the eruption and keep the sign of their axis writhe (which is expected if field of one chirality dominates the source volume of the eruption), must reverse their S shape in the course of the rise. Implications for the occurrence of the helical kink instability in such events are discussed.

WRITHE↗
NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

boundary layer transition↗