Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,099 records · Page 61

A numerical investigation of hydrogen combustion in Mach 2 flow

A numerical study is conducted of a reacting flowfield generated by a single hydrogen jet injected transversely behind a rearward facing step. A nonreacting mixing study was followed by a reacting simulation. The computation simulates a current experimental condition. The validation of the numerical solution was limited to a comparison of the wall-pressure distribution along the wall. The comparison indicates agreement between the experiment and the prediction with a maximum error of 7 percent. The addition of the reaction responsible for HO2 production (H + O2 + M = HO2 + M) to an initially simplified chemical mechanism reduced the exponential growth of the free radicals and reduced the amount of heat released by combustion; as a result, the flowfield was significantly changed. A reversed flow was obtained in limited regions of the flow.

Segal, Corin↗

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method using multifidelity models for CFD-based multiobjective MDO. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

multidisciplinary optimization↗

Aeroelastic Simulation of Wind Tunnel Scale Aspect Ratio 13.5 Common Research Model with Trailing Edge Mini-Flaps

This paper presents a computation study of an aspect ratio 13.5 wind tunnel scale Common Research Model (CRM) with trailing edge mini-flaps using FUN3D. The aeroelastic simulations are conducted under free-air and inside wind tunnel conditions. The NASA Langley Research Center Transonic Dynamic Tunnel (TDT) is a closed circuit and slotted-test-section with a 16 ft ×16 ft test section. The structure mode shapes are extracted from the latest updated NASTRAN Finite-Element Model (FEM) model. A dynamic pressure sweep is performed to determine the control reversal boundary of the trailing edge mini-flaps. The wind tunnel wall interference effects on the mini-flap control derivatives are analyzed based on the simulations.

Aeroelastic↗

Closing the loop: model-predictive control for a closed-circuit reverse osmosis system

This article presents a model-predictive controller (MPC) for the maximization of the energy efficiency of a closed-circuit desalination reverse osmosis (CCRO) system. CCRO is a process for producing drinking water that is based on a cyclic operation with the following two phases: (a) filtration and (b) drain. In this article, we test model predictive control for optimal control of this process. The most important features of our approach are as follows: (a) the selection of a model structure that enables reliable forecasts of the filtration phase (up to 3 h), (b) an on-line model calibration strategy that ensures model forecast reliability, and (c) the satisfaction of equipment safety and operational constraints on the selected setpoints. We challenge this through deliberate introduction of changes in the unmeasured feed concentration and the applied constraints. Our results indicate that frequent model parameter updates are critical to maintain model reliability for MPC purposes. In addition, we illustrate that parameter identifiability is not guaranteed and that deliberate variation in flow rates is necessary even though the process never operates in steady state. Finally, MPC can compute flow rate setpoints that maximize the energy efficiency of the CCRO process while satisfying the applicable equipment and safety constraints.

closed-circuit reverse osmosis↗

DC and Structured Electric Fields Observed on the C/NOFS Satellite and Their Association with Longitude, Plasma Density, and Solar Activity

Observations of DC electric fields and associated E x B plasma drifts gathered by the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite are presented. We show statistical averages of the vector fields and resulting E x B plasma flows for the first three years of operations as a function of season, longitude, local time, and Fl 0.7 conditions. Magnetic field data from the VEFI science magnetometer are used to compute the plasma flows. Although typically displaying eastward and outward-directed fields during the day and westward and downward-directed fields at night, the data from DC electric field detector often reveal variations from this pattern that depend on longitude, solar activity, and plasma density. Clear "wave-4" tidal effects in both electric field components have been detected and will be presented. Zonal plasma drifts show a marked variation with solar activity and may be used as a proxy for neutral winds at night. Evidence for pre-reversal enhancements in the meridional drifts that depend on solar activity is present for some longitudes, and are corroborated by clear evidence in the plasma density data that the spacecraft journeyed below the F-peak during evenings when the rise in the ionosphere is most pronounced. In addition to DC electric fields, the data reveal considerable electric field structures at large scales (approx 100's of km) that are usually confined to the nightside. Although such electric field structures are typically associated with plasma density depletions and structures, what is surprising is the number of cases in which large amplitude, structured DC electric fields are observed without a significant plasma density counterpart structure, including their appearance at times when the ambient plasma density appears relatively quiescent. We investigate the mapping of structured electric fields along magnetic field lines from distant locations and consider tropospheric thunderstorm sources to explain some of the observations. The dependence of the structuring on season and solar activity will be reported and discussed.

Pfaff, Robert↗

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Blowing polar skyrmion bubbles in oxide superlattices

Particle-like topological structures such as skyrmions and vortices have garnered ever-increasing interests due to their rich physical insights and potential broad applications in spintronics. Here, in this work, we discover the reversible switching between polar skyrmion bubbles and ordered vortex arrays in ferroelectric superlattices under an electric field, reminiscent of the Plateau-Raleigh instability in fluid mechanics. An electric field phase diagram is constructed, showing a wide stability window for the observed polar skyrmions. A “volcano”-like pontryagin density distribution is formed, indicating the formation of a smooth circular skyrmion. The topological charge Q at different applied field is calculated, verifying the field-driven topological transition between Q = 0 and Q = ±1 states. This study is a demonstration for the computational design of field-induced topological phase transitions, giving promise for the design of next-generation nanoelectronic devices.

36 MATERIALS SCIENCE↗

Observation of plaid-like spin splitting in a noncoplanar antiferromagnet

Spatial, momentum and energy separation of electronic spins in condensed matter systems guides the development of novel devices where spin-polarized current is generated and manipulated. Recent attention on a set of previously overlooked symmetry operations in magnetic materials leads to the emergence of a new type of spin splitting, enabling giant and momentum-dependent spin polarization of energy bands on selected antiferromagnets. Despite the ever-growing theoretical predictions, the direct spectroscopic proof of such spin splitting is still lacking. Here, we provide solid spectroscopic and computational evidence for the existence of such materials. In the noncoplanar antiferromagnet MnTe 2 , the in-plane components of spin are found to be antisymmetric about the high-symmetry planes of the Brillouin zone, comprising a plaid-like spin texture in the antiferromagnetic (AFM) ground state. Such an unconventional spin pattern, further found to diminish at the high- temperature paramagnetic state, stems from the intrinsic AFM order instead of spin- orbit coupling (SOC). Furthermore, our finding demonstrates a new type of quadratic spin texture induced by time-reversal breaking, placing AFM spintronics on a firm basis and paving the way for studying exotic quantum phenomena in related materials.

Electronic devices↗

Impact of resonant magnetic perturbations on zonal flows and microturbulence

Results from a joint experimental and computational effort studying the effect of resonant magnetic perturbations (RMPs) on microturbulence levels and their connection to zonal flows in the DIII-D tokamak L-mode are presented. Beam emission spectroscopy measurements show a direct increase in density fluctuations at microturbulent scales with increasing RMP amplitude, suggesting that magnetic activity introduced by the RMP affects the regulation of microturbulence on DIII-D. This is analogous to how MHD-scale magnetic fluctuations arising from tearing modes have been observed in simulations to increase microturbulence levels in the reversed-field pinch (RFP). In the RFP, this is attributed to magnetic fluctuations eroding turbulence-limiting zonal flows; this work examines if a similar mechanism is present for DIII-D microturbulence. Gyrokinetic simulations find that the application of an RMP corresponds directly to a decrease in zonal flow levels, producing a similar increase of turbulent fluctuation levels over a range of RMP amplitudes as observed in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dopant Optimization of Donors in Semiconductor Opening Switches to Eliminate Prepulse

Semiconductor opening switches are solid-state devices capable of delivering nanosecond, hundreds of kilovolts pulses by interrupting kiloamps of current. The interruption of the current occurs in a moderately doped p-region when a high electric field region (HFR) is formed. The HFR occurs because the reverse pumping current cannot be supported by the saturation velocity and majority carrier concentration of the doping level. However, the donor profile also significantly affects the pulse performance. A secondary prepulse occurs if a secondary HFR is formed at the interface of the background n-doping and N+ doping (X n ) . By moving the location of X n deeper into the diode, the effect of the prepulse is reduced. This article investigates the effect of the donor doping profile on the performance metrics of semiconductor opening switches through technology computer-aided design (TCAD) simulations and experimental results. Through a SILVACO TCAD optimization, we designed a P + /p/n - base/n/N + where the intersection of the moderate p-region and intrinsic n-base region (X p ) is at 160 μm and X n is at 220 μm. This profile is fabricated via silicon epitaxy. Experimentally, it is shown that a deep X n (220 μm) compared with a shallow X n (300 μm) reduces the rise time by >5× . In addition, the magnitude of current density during interruption affects the prepulse foot and pulse shape. At lower current densities without the graded donor profile, high peak voltages are not achieved. Comparing the experimental results to the TCAD simulations shows that the model is predictive under high-current densities in the semiconductor opening switch (SOS) regime.

nanosecond pulse power↗

Ionospheric convection driven by NBZ currents

Computer simulations of Birkeland currents and electric fields in the polar ionosphere during periods of northward IMF were conducted. When the IMF z component is northward, an additional current system, called the NBZ current system, is present in the polar cap. These simulations show the effect of the addition of NBZ currents on ionospheric convection, particularly in the polar cap. When the total current in the NBZ system is roughly 25 to 50 percent of the net region 1 and 2 currents, convection in the central portion of the polar cap reverses direction and turns sunward. This creates a pattern of four-cell convection with two small cells located in the polar cap, rotating in an opposite direction from the larger cells. When the Birkeland currents are fixed (constant current source), the electric field is reduced in regions of relatively high conductivity, which affects the pattern of ionospheric convection. Day-night asymmetries in conductivity change convection in such a way that the two polar-cap cells are located within the large dusk cell. When ionospheric convection is fixed (constant voltage source), Birkeland currents are increased in regions of relatively high conductivity. Ionospheric currents, which flow horizontally to close the Birkeland currents, are changed appreciably by the NBZ current system. The principal effect is an increase in ionospheric current in the polar cap.

Rasmussen, C. E.↗

Alpha – omega and omega – alpha phase transformations in zirconium under hydrostatic pressure: A 3D mesoscale study

Here, a three dimensional (3D) elastoplastic phase-field model is developed for modeling the hydrostatic pressure-induced alpha – omega phase transformation and the reverse phase transformation, i.e. omega – alpha, in zirconium (Zr). Plastic deformation and strain hardening of the material are also considered in the model. The microstructure evolution during both phase transformations is studied. The transformation start pressures at different temperatures are predicted and are plotted as a phase diagram. The effect of phase transformations on the mechanical properties of the material is also studied. The input data corresponding to pure Zr are acquired from experimental studies as well as by using the CALPHAD method. Our simulations show that three different omega variants form as laths. On release of pressure, reverse phase transformation initiates at lath boundaries. We observe that both phase transformations are martensitic in nature and also occur at the same pressure, i.e. little hysteresis. The transformation start pressures and the kinetics of the transformation predicted by our model are in good agreement with experimental results.

36 MATERIALS SCIENCE↗

Note: A pure-sampling quantum Monte Carlo algorithm with independent Metropolis

Recently, Ospadov and Rothstein published a pure-sampling quantum Monte Carlo algorithm (PSQMC) that features an auxiliary Path Z that connects the midpoints of the current and proposed Paths X and Y, respectively. When sufficiently long, Path Z provides statistical independence of Paths X and Y. Under those conditions, the Metropolis decision used in PSQMC is done without any approximation, i.e., not requiring microscopic reversibility and without having to introduce any G(x → x′; τ) factors into its decision function. This is a unique feature that contrasts with all competing reptation algorithms in the literature. An example illustrates that dependence of Paths X and Y has adverse consequences for pure sampling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atmospheric chemistry of the reaction ClO + O2 reversible reaction ClO (center dot) O2: Where it stands, what needs to be done, and why?

Possible existence and chemistry of ClO (center dot) O2 was originally proposed to explain the Norrish-Neville effect that O2 suppresses chlorine photosensitized loss of ozone. It was also thought that ClO (center dot) O2 might have some atmospheric chemistry significance. Recently, doubts have been cast on this proposal, because certain laboratory data seem to imply that the equilibrium constant of the title reaction is so small that ClO (center dot) O2 may be too unstable to matter. However, those data create only a superficial illusion to that effect, because on a closer analysis they do not disprove a moderately stable and chemically significant ClO (center dot) O2. Furthermore, our state-of-the-science accurate computational chemistry calculations also suggest that ClO (center dot) O2 may be a weakly bound ClOOO radical with a reactive (2)A ground electronic state. There is therefore a need to design and perform definitive experimental tests of the existence and chemistry of the ClO (center dot) O2 species, which we discuss and which have the potential to mediate the chlorine-catalyzed stratospheric ozone depletion.

Prasad, Sheo S.↗

Negative magnetic eddy diffusivities from the test-field method and multiscale stability theory

The generation of a large-scale magnetic field in the kinematic regime in the absence of an α-effect is investigated by following two different approaches: the test-field method and the multiscale stability theory relying on the homogenization technique. Our computations of the magnetic eddy diffusivity tensor of the parity-invariant flow IV of G. O. Roberts and the modified Taylor–Green flow confirm the findings of previous studies and also explain some of their apparent contradictions. The two flows have large symmetry groups; this is used to considerably simplify the eddy diffusivity tensor. Finally, a new analytic result is presented: upon expressing the eddy diffusivity tensor in terms of solutions to auxiliary problems for the adjoint operator, we derive relations between the magnetic eddy diffusivity tensors that arise for mutually reverse small-scale flows v(x) and −v(x).

79 ASTRONOMY AND ASTROPHYSICS↗

3D numerical investigation of a free-burning argon arc with metal electrodes using a novel sheath coupling procedure

In this study, a novel coupling procedure describing the complicated plasma-electrode interaction process has been developed and applied into the 3D finite volume simulation of a direct current tungsten inert gas (TIG) welding system. This is achieved by making the space charge layer (sheath) incorporated into the computation domain and interact with both bulk plasma and cathode through the effective electrical conductivity. Both chemical and thermal nonequilibrium phenomena as well as the self-induced magnetic fields have been taken into consideration by the model to ensure a realistic numerical description of a non-thermal arc. The applicability of this coupling procedure is further improved by calculating the real electric potential, which is capable of accounting for the effects of the complicated drift and diffusion processes. Numerical results of both 100 and 200A discharge currents are presented, field reversal is obtained at near-anode regions in both cases, which is followed by the negative anode sheath potential drop. The region of the strongest electron overpopulation appears at the intersection of plasma fringes and electrode surface. Finally, numerical results of plasma temperature and voltage show good agreement with experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗

Linear system identification via backward-time observer models

Presented here is an algorithm to compute the Markov parameters of a backward-time observer for a backward-time model from experimental input and output data. The backward-time observer Markov parameters are decomposed to obtain the backward-time system Markov parameters (backward-time pulse response samples) for the backward-time system identification. The identified backward-time system Markov parameters are used in the Eigensystem Realization Algorithm to identify a backward-time state-space model, which can be easily converted to the usual forward-time representation. If one reverses time in the model to be identified, what were damped true system modes become modes with negative damping, growing as the reversed time increases. On the other hand, the noise modes in the identification still maintain the property that they are stable. The shift from positive damping to negative damping of the true system modes allows one to distinguish these modes from noise modes. Experimental results are given to illustrate when and to what extent this concept works.

Juang, Jer-Nan↗