SEARCH · Engineering Papers
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.
Ground level runway wind statistics cover crosswind, tailwind, and headwind reversal percentage frequencies with respect to month and hour for the two major runways. Also presented are bivariate normal wind statistics for a 90 degree flight azimuth for altitudes 0 through 27 km. Wind probability distributions, synthetic vector wind profiles, and statistics for any rotation of axes are computed from five given parameters.
A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.
A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.
An intrusion of vortex edge air in D the interior of the Arctic polar vortex was observed on the January 31,2005 flight of the NASA DC-8 aircraft. This intrusion was identified as anomalously high values of ozone by the AROTAL and DIAL lidars. Our analysis shows that this intrusion formed when a blocking feature near Iceland collapsed, allowing edge air to sweep into the vortex interior. along the DC-8 flight track also shows the intrusion in both ozone and HNO3. Polar Stratospheric Clouds (PSCs) were observed by the DIAL lidar on the DC-8. The spatial variability of the PSCs can be explained using MLS HNO3 and H2O observations and meteorological analysis temperatures. We also estimate vortex denitrification using the relationship between N2O and HNO3. Reverse domain fill back trajectory calculations are used to focus on the features in the MLS data. The trajectory results improve the agreement between lidar measured ozone and MLS ozone and also improve the agreement between the HNO3 measurements PSC locations. The back trajectory calculations allow us to compute the local denitrification rate and reduction of HCl within the filament. We estimate a denitrification rate of about lO%/day after exposure to below PSC formation temperature. Analysis of Aura MLS observations made
A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.
We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.
Two-terminal memory elements, or memelements, capable of co-locating signal processing and memory via history-dependent reconfigurability at the nanoscale are vital for next-generation computing materials striving to match the brain’s efficiency and flexible cognitive capabilities. While memory resistors, or memristors, have been widely reported, other types of memelements remain underexplored or undiscovered. Here we report the first example of a volatile, voltage-controlled memcapacitor in which capacitive memory arises from reversible and hysteretic geometrical changes in a lipid bilayer that mimics the composition and structure of biomembranes. We demonstrate that the nonlinear dynamics and memory are governed by two implicitly-coupled, voltage-dependent state variables—membrane radius and thickness. Further, our system is capable of tuneable signal processing and learning via synapse-like, short-term capacitive plasticity. These findings will accelerate the development of low-energy, biomolecular neuromorphic memelements, which, in turn, could also serve as models to study capacitive memory and signal processing in neuronal membranes.
This paper provides theoretical and numerical justifications for inverse design optimization of reversed equivalent area (A e,r ) of low-boom supersonic aircraft when the cruise condition changes during the optimization iterations. A modified linear theory for steady flow around a supersonic projectile is used to establish the accuracy of A e,r -based body-of-revolution approximation of a low-boom supersonic aircraft for undertrack sonic boom analysis using computational fluid dynamics off-body pressure. As a result, designing a low-boom shape of the off-body pressure at three body lengths below the aircraft is equivalent to designing a low-boom A e,r shape. A Bezier curve with eight control points is used to define an A e,r target for inverse design optimization of A e,r of a supersonic aircraft. To make an A e,r target matchable by a low-boom supersonic aircraft, the A e,r target must have the minimal perceived level of decibels (PLdB) for sonic boom and satisfy two constraints defined by the cruise condition, effective length, and two nondimensional parameters. Numerical results are used to verify that, for different cruise conditions, all optimized A e,r targets with PLdB below 70 can be approximately generated by scaling one of them if the targets have the same nondimensional parameters and approximately the same effective length.
In this work, an accident tolerant fuel cladding system with Cr coating deposited using physical vapour deposition (PVD) method onto commercial Optimized ZIRLO™ was studied. The cladding tubes were machined into C-rings and loaded to failure under compression with real-time synchrotron micro-X-ray computed tomography (XCT) imaging at room temperature (RT), 345 °C, 650 °C and 950 °C in argon (Ar) atmosphere. The mechanical behaviour and failure processes were found strongly temperature-dependent where the Cr coating showed brittle fracture at RT and 345 °C, ductile fracture at 650 °C and a reversion to brittle fracture at 950 °C. Nanoindentation measurements and scanning electron microscopy (SEM) imaging were conducted on the materials after high temperature testing. It was found that recrystallisation of the Cr coating occurred at elevated temperatures, which significantly affected its local properties hence the failure behaviour at different temperatures. This work represents the first in situ 3D XCT observation of progressive failure processes in PVD Cr-coated ZIRLO claddings up to 950 °C providing critical insights into its brittle-to-ductile transition (BDTT) behaviour and subsequent ductile-to-brittle reversion with the increase in temperature. Moreover, results are compared with other PVD Cr-coated Zircaloy materials from open literature; the influences of the coating microstructure and local properties on the failure stress/strain and fracture processes are discussed.
There is an increasing need for sustainable energy storage solutions as fossil fuels are replaced by renewable energy sources. Multivalent batteries, specifically Mg batteries, are one energy storage technology that researchers continue to develop with hopes to surpass the performance of Li-ion batteries. However, the limited energy density and transport properties of Mg cathodes remain critical challenges preventing the realization of high-performance multivalent batteries. In this work, ABO 4 zircon materials (A = Y, Eu and B = V, Cr) are computationally and experimentally evaluated as Mg intercalation cathodes. Remarkably good Mg-ion transport properties were predicted and Mg-ion intercalation was experimentally verified in sol-gel synthesized zircon YVO 4 , EuVO 4 , and EuCrO 4 . Among them, EuVO 4 exhibited the best electrochemical performance and demonstrated repeated reversible cycling. While we believe that the one-dimensional diffusion channels and redox-active species tetragonal coordination limit the value of many zircons as high-performance cathodes, their unique structural motif of overlapping polyhedra along the diffusion pathway appears instrumental for promoting good Mg-ion mobility. The motif results in a favorable "6-5-4" change in coordination that avoids unfavorable sites with lower coordination along the diffusion pathway and a structural design metric for future Mg cathode development.
The synthesis, spectroscopic characterization, computational analysis, and photoelectrochemical behavior of a new iron-based chromophore, [(Cpy) 2 Fe- (deeb)](PF 6 ) 2 (Fe(Cpy) 2 (deeb)), where Cpy is 1-methyl-3-(2-pyridyl)imidazole and deeb is 4,4’-(CO 2 CH 2 CH 3 ) 2 -2,2’-bipyridine, is reported. Electrochemically reversible waves assigned to a metal-centered E o (Fe III/II ) = +0.48 and a ligand-centered E o (Fe 2+/+ ) = −1.47 V vs Fc +/0 reduction were evident in cyclic voltammetry measurements. The combination of a strong σ-donor and a π-acceptor lowered the energy of the metal-toligand charge-transfer (MLCT) excited state relative to the metal-centered state. Two MLCT transitions appear in the visible region at 424 and 580 nm. TDDFT calculations revealed that the lower-energy band was well formulated as Fe(II)→deeb, and the higher-energy transition was charge transfer to both the deeb and Cpy ligands. Resonance Raman spectroscopy supports these findings showing enhanced deeb vibrational modes with 532 nm excitation, both deeb and Cpy modes with 473 nm excitation, and exclusively Cpy with 405 nm excitation. Ultrafast spectroscopy reveals a short-lived (∼2 ps) MLCT excited state and a longer-lived (∼20 ps) metal-centered state. Efficient methods to deprotect the ester groups and anchor the complex to mesoporous TiO 2 (anatase) thin films in high surface coverages, Fe(Cpy) 2 (dcb)|TiO 2 σ = 3 × 10 −8 mol/cm 2 , were established. Pulsed light excitation of Fe(Cpy) 2 (dcb)|TiO 2 resulted in rapid excited state injection (k inj > 10 8 s −1 ) and formation of a charge-separated state, Fe III (Cpy) 2 (dcb)|TiO 2 (e), which persists on the millisecond time scale before returning cleanly to the ground state with secondorder kinetics. Injection yields measured 50 ns after light excitation were found to double from Φ = 0.15 with green (532 nm) light to 0.30 with blue (457 nm) light excitation. Incident photon-to-current efficiency (% IPCE) measurements as a function of excitation wavelength in a 0.5 M LiI/I 2 /CH 3 CN electrolyte provide clear evidence for band-selective “hot carrier” injection from the remote Cpy-localized excited state. Collectively, the spectroscopic and photoelectrochemical data indicate that a semiconductor can intercept hot electrons from iron chromophores even when the excited-state dipole is oriented away from the surface-anchoring ligand.
Zinc continues to garner immense interest due to its versatility as an anode material in several configurations utilizing either alkaline or mild-pH electrolytes. Current research on using mild-pH electrolytes has improved the rechargeability aspect of Zn-based batteries since Zn2+ is solely utilized for plating/stripping of Zn. Several studies have incorporated Zn metal foils, yet, dramatic improvements can be achieved by expressing Zn as a porous structure. Herein, we use a quasi-pulsed electrodeposition process to prepare a conformal Zn coating onto 3D porous copper foam. By tuning the electrodeposition parameters, we achieved an optimal Zn coating that undergoes reversible plating/stripping when tested in symmetric Zn cells, which supported a low overpotential of ~60 mV for up to 100 cycles. We further investigated changes in the surface morphology by studying the Zn surface of both foil and 3D structure using scanning electron microscopy and X-ray micro-computed tomography. Both techniques showed that the Zn foil undergoes dramatic alterations at the surface, which results in inhomogeneous deposition of Zn, whereas the 3D form exhibited minimal changes. Lastly, we paired both Zn foil and 3D Zn with vanadium oxide and demonstrated that the porous structure supports high rate capability and high specific capacity.
An MHD unsteady 1-D model is used to simulate the interaction and coalescence of two pressure waves in the outer heliosphere. Each of the two pressure waves was a compression region bounded by a shock pair. Computer simulation using Voyager data as input demonstrates the interaction and coalescence process involving one pressure wave associated with a fast stream and the other pressure wave without a fast stream. The process produced a significant change in the magnetic field and plasma signatures. The propagation of the forward and reverse shocks first widened the radial dimension of the shock compression region with increasing heliocentric distances. The shocks belonging to two neighboring compression regions eventually collided and two compression regions began to overlap with each other. This type of interaction is a dominant dynamical process in the outer heliosphere, and significantly and irreversible alters the structure of the medium.
An MHD unsteady 1-D model is used to simulate the interaction and coalescence of two pressure waves in the outer heliosphere. Each of the two pressure waves was a compression region bounded by a shock pair. Computer simulation using Voyager data as input demonstrates the interaction and coalescence process involving one pressure wave associated with a fast stream and the other pressure wave without a fast stream. The process produced a significant change in the magnetic field and plasma signatures. The propagation of the forward and reverse shocks first widened the radial dimension of the shock compression region with increasing heliocentric distances. The shocks belonging to two neighboring compression regions eventually collided and two compression regions began to overlap with each other. This type of interaction is a dominant dynamical process in the outer heliosphere, and significantly and irreversible alters the structure of the medium.
Various features of the computer codes used in the helicopter industry and by government agencies for rotorcraft aeroelastic stability analysis are compared. Mathematical rigor in modeling rotorcraft is given primarily to the rotor-system dynamic behavior; the aerodynamic modeling is still limited to strip theory and to uneven application of corrections for stall, reversed flow, yawed flow, radial flow, and unsteady aerodynamic effects. The forward-flight regime analysis is included in five of the 11 codes surveyed. However, only two of these codes are capable of a Floquet analysis for aeroelastic stability. For the hover regime, nine of the 11 codes use eigen-analysis approach. The remaining codes perform a harmonic analysis of the transient response of system.
A 6 DOF analytical aerodynamic model of a high alpha research vehicle is derived. The derivation is based on wind-tunnel model data valid in the altitude-Mach flight envelope centered at 15,000 ft altitude and 0.6 Mach number with Mach range between 0.3 and 0.9. The analytical models of the aerodynamics coefficients are nonlinear functions of alpha with all control variable and other states fixed. Interpolation is required between the parameterized nonlinear functions. The lift and pitching moment coefficients have unsteady flow parts due to the time range of change of angle-of-attack (alpha dot). The analytical models are plotted and compared with their corresponding wind-tunnel data. Piloted simulated maneuvers of the wind-tunnel model are used to evaluate the analytical model. The maneuvers considered are pitch-ups, 360 degree loaded and unloaded rolls, turn reversals, split S's, and level turns. The evaluation finds that (1) the analytical model is a good representation at Mach 0.6, (2) the longitudinal part is good for the Mach range 0.3 to 0.9, and (3) the lateral part is good for Mach numbers between 0.6 and 0.9. The computer simulations show that the storage requirement of the analytical model is about one tenth that of the wind-tunnel model and it runs twice as fast.
To date, most large-eddy simulations (LES) have been carried out with eddy viscosity subgrid scale (SGS) models, with only a few exceptions that used the mixed model. Even though the assumptions behind Smagorinsky's model are rather stringent, it has been applied successfully to a variety of turbulent flows. This success is attributed to the ability of eddy viscosity models to drain energy from large scales, thus simulating the dissipative nature of turbulence. Most SGS models are absolutely dissipative, i.e. they remove energy from the large scales at every instant. However, SGS stresses may transfer energy back to the large scales intermittently; this reverse transfer or backscatter is especially important in geophysical flows and in transition. In a fully developed channel flow, there is reverse flow of energy from small to large scales near the walls, but eddy viscosity models are unable to account for this important feature. The dynamic localization eddy viscosity model of Ghosal et al. (1995) allows backscatter by co-evolving an auxiliary equation for the SGS energy; however, the computational cost is considerably larger than for conventional SGS models (Cabot 1994). In this report, a new non-eddy viscosity model based on local approximation of total quantities in terms of filtered ones is introduced; the scale similarity model of Bardina (1983) is a special case of this model. This procedure does not require the assumption of homogeneity, permits backscatter of energy from small to large scales, and is readily implemented in finite difference codes. The results of applying the proposed model to second order finite volume simulation of plane channel flow at high Reynolds numbers (Re(sub b) = 38,000) is described in this report. Greater emphasis is placed on the high Reynolds number flow since it provides a more rigorous test of the SGS model and its potential application. The results are compared to ones produced by the conventional and dynamic Smagorinsky models and the spectral LES of Piomelli (1993).
Neodymium-iron-boron (NdFeB) magnets are the most powerful magnets per unit volume sold in the commercial market. Despite the increasing demand for clean energy applications such as electric vehicles and wind turbines, disruptive events including the COVID-19 pandemic have caused significant uncertainties in the supply and demand for NdFeB magnets. Therefore, this study aims to alleviate the risk of supply shortage for NdFeB magnets and the containing critical materials, rare-earth elements (REEs), through the development of a resilient reverse supply chain and logistics network design. We develop scenarios to model the unique impact of the COVID-19 pandemic on the proposed business, incorporating both disruption intensity and recovery rate. We formulate a chance-constrained two-stage stochastic programming model to maximize the profit while guaranteeing the network resiliency against disruption risks. To solve the problem in large-scale instances, we develop an efficient Benders decomposition algorithm that reduces the computational time by 98.5% on average compared to the default CPLEX algorithm. When applied to the United States, the model suggests the optimal facility locations, processing capacities, inventory levels, and material flows for NdFeB magnet recyclers that could meet 99.7% of the demand. To the best of our knowledge, this study is the first to incorporate the impacts of the COVID-19 pandemic to design a resilient NdFeB magnet recycling supply chain and logistics network, leveraging risk-averse stochastic programming.