Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion maps”

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 613 records · Page 34

An experimental study of opposed flow diffusion flame extinction over a thin fuel in microgravity

The flame spread and flame extinction characteristics of a thin fuel burning in a low-speed forced convective environment in microgravity were examined. The flame spread rate was observed to decrease both with decreasing ambient oxygen concentration as well as decreasing free stream velocity. A new mode of flame extinction was observed, caused by either of two means: keeping the free stream velocity constant and decreasing the oxygen concentration, or keeping the oxygen concentration constant and decreasing the free stream velocity. This extinction is called quenching extinction. By combining this data together with a previous microgravity quiescent flame study and normal-gravity blowoff extinction data, a flammability map was constructed with molar percentage oxygen and characteristic relative velocity as coordinates. The Damkohler number is not sufficient to predict flame spread and extinction in the near quench limit region.

Ferkul, Paul V.↗

Exchange of Ions across the TiN/TaO x Interface during Electroformation of TaO x -Based Resistive Switching Devices

The valence change model describes the resistive switching in metal oxide-based devices as due to electroreduction of the oxide and subsequent electromigration of oxygen vacancies. Here, we present cross-sectional X-ray energy-dispersive spectroscopy elemental maps of Ta, O, N, and Ti in electroformed TiN/TaO 2.0 /TiN structures. O, N, and Ti were exchanged between the anode and the functional oxide in devices formed at high power (~1 mW), but the exchange was below the detection limit at low power (<0.5 mW). All structures exhibit a similar Ta-enriched and O-depleted filament formed by the elemental segregation in the functional oxide by the temperature gradient. The elemental interchange is interpreted as due to Fick’s diffusion caused by high temperatures in the gap of the filament and is not an essential part of electroformation.

36 MATERIALS SCIENCE↗

Premixed Gas Combustion: An Excitable System

Rotating spiral and target patterns have been observed experimentally on freely-propagating premixed gas flames in large diameter tubes at normal gravity (1-g). These modes of propagation occur in near-limit mixtures which have a Lewis number (Le, defined as the ratio of the thermal diffusivity of the cold mixture to the mass diffusivity of the scarce component into the mixture) sufficiently greater than one. However, at 1-g, buoyant flows strongly distort the flame curvature, hydrodynamics (thus stretch) and convective transport of species and heat. In turn, these alter the critical Le required for onset of instability. To isolate and better understand the mechanisms which drive the observed patterns and their dynamics, 1-g and microgravity (micro-g) experiments are being conducted to determine: (1) the structure and dynamics of the patterns, (2) a map of the critical Le and heat loss for their occurrence, (3) the relative significance of the chemical kinetics, and (4) the effect of curvature (local wave and global flame front) on wave propagation. With this in hand, we will be better prepared to discuss an additional mode, a state of 'chemical turbulence,' which seems to be the ultimate fate of many of these near-limit flames prior to extinction.

Pearlman, Howard↗

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

Evolution of Local Microstructures: Spatial Instabilities in Coarsening Clusters

Diffusion-limited capillarity-driven coarsening of precipitates is an important and intensively studied phenomenon. The classic coarsening theory developed by Lifshitz and Syozov and Wagner (LSW theory) is limited to infinitesimally small volume fractions, V(sub nu), therefore neglects all direct interparticle interactions. This work uses modeling and holographic imaging to compare coarsening rates in "high" volume fraction versus low volume fraction microstructures by observing mixed-dimensional droplets (spherical caps on a surface coarsening by two-dimensional diffusion) during ground-based investigations. The method involves filling a cell with selected homogeneous parent phase, and cooling below the consolute temperature to the isopycnic temperature in the two-phase region of a monotectic system. A microgravity holographic experiment is required for three-dimensional observations to minimize sedimentation during long-term coarsening. Determination of sizes and positions of the many droplets in the holographic images requires automation. We have developed software for automated data analysis, and demonstrated good agreement between regenerated maps and scaled photographs of the original holograms for mixed dimensional coarsening. The results of these experiments were presented in a formal microgravity Science Concept Review (SCR) on December 18, 2000.

Source record↗

Localization crossover and subdiffusive transport in a classical facilitated network model of a disordered interacting quantum spin chain

Here, we consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL) and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's noninteracting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady-state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath and provides a classical analog for the MBL transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗

Effect of Fluorination on Lithium Transport and Short-Range Order in Disordered-Rocksalt-Type Lithium-Ion Battery Cathodes

Fluorine substitution is a critical enabler for improving the cycle life and energy density of disordered rocksalt (DRX) Li-ion battery cathode materials which offer prospects for high energy density cathodes, without the reliance on limited mineral resources. Due to the strong Li–F interaction, fluorine also is expected to modify the short-range cation order in these materials which is critical for Li-ion transport. In this work, density functional theory and Monte Carlo simulations are combined to investigate the impact of Li–F short-range ordering on the formation of Li percolation and diffusion in DRX materials. Additionally, the modeling reveals that F substitution is always beneficial at sufficiently high concentrations and can, surprisingly, even facilitate percolation in compounds without Li excess, giving them the ability to incorporate more transition metal redox capacity and thereby higher energy density. It is found that for F levels below 15%, its effect can be beneficial or disadvantageous depending on the intrinsic short-range order in the unfluorinated oxide, while for high fluorination levels the effects are always beneficial. Using extensive simulations, a map is also presented showing the trade-off between transition-metal capacity, Li-transport, and synthetic accessibility, and two of the more extreme predictions are experimentally confirmed.

25 ENERGY STORAGE↗

Elemental segregation to lattice defects in the CrMnFeCoNi high-entropy alloy during high temperature exposures

The influence of small plastic pre-strains on the elevated-temperature stability and microstructure of the equiatomic CrMnFeCoNi FCC solid solution is investigated. Particular attention is given to whether any of the alloy elements segregate to individual dislocations. To that end, CrMnFeCoNi samples were first deformed in tension at room temperature to plastic strains of 0.2 and 2.3%, and subsequently annealed at 973 K for 800 hours. The pre-strains activated planar slip of 1/2<110>-type dislocations on {111}-type glide planes. Interactions of this planar slip with special Σ3 grain boundaries formed a large number of dislocation segments with a <110>-type crystallographic orientation suitable for a credible end-on analysis of dislocation cores in HR-STEM. The cores of the 1/2<110> dislocations pushed up against the investigated grain boundaries were found to be close to the compact configuration. Within the sensitivity of the Super-X EDS mapping, no concentration gradient was detected near dislocations that would indicate enrichment at dislocation cores of any of the elemental constituents of the alloy after the pre-deformation and annealing. However, a Cr-rich tetragonal sigma phase nucleated and grew at grain boundary triple junctions during this anneal, processes that were not accelerated by the enhanced dislocation density present after pre-strain. A clear chromium gradient was observed in the Cr-depleted zones near grain boundaries suggesting that Cr transport occurred by relatively slow diffusion from the bulk to the grain boundaries and then by relatively fast diffusion along the grain boundaries to the precipitates. Accompanying the Cr depletion near grain boundaries is a simultaneous Ni and Mn enrichment, which promotes formation of the L1 0 NiMn phase that is observed on the grain boundaries after prolonged annealing.

36 MATERIALS SCIENCE↗

High-speed mapping of surface charge dynamics using sparse scanning Kelvin probe force microscopy

Unraveling local dynamic charge processes is vital for progress in diverse fields, from microelectronics to energy storage. This relies on the ability to map charge carrier motion across multiple length- and timescales and understanding how these processes interact with the inherent material heterogeneities. Towards addressing this challenge, we introduce high-speed sparse scanning Kelvin probe force microscopy, which combines sparse scanning and image reconstruction. This approach is shown to enable sub-second imaging (>3 frames per second) of nanoscale charge dynamics, representing several orders of magnitude improvement over traditional Kelvin probe force microscopy imaging rates. Bridging this improved spatiotemporal resolution with macroscale device measurements, we successfully visualize electrochemically mediated diffusion of mobile surface ions on a LaAlO 3 /SrTiO 3 planar device. Such processes are known to impact band-alignment and charge-transfer dynamics at these heterointerfaces. Furthermore, we monitor the diffusion of oxygen vacancies at the single grain level in polycrystalline TiO 2 . Through temperature-dependent measurements, we identify a charge diffusion activation energy of 0.18 eV, in good agreement with previously reported values and confirmed by DFT calculations. Together, these findings highlight the effectiveness and versatility of our method in understanding ionic charge carrier motion in microelectronics or nanoscale material systems.

25 ENERGY STORAGE↗

Lateral charge transport from heavy-ion tracks in integrated circuit chips

A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.

Zoutendyk, J. A.↗

STAR-X: Survey and Time-domain Astrophysical Research eXplorer

STAR-X is a MIDEX mission proposed to NASA in December 2021. Comprising an X-ray telescope (XRT) provided by GSFC and MIT, a UV telescope (UVT) provided by the University of Colorado, and a spacecraft (SC) provided by Ball Aerospace, STAR-X is designed to conduct time-domain survey and to respond rapidly to transient events discovered by other observatories such as LIGO, Rubin LSST, Roman WFIRST, and SKA. The key features of the XRT are its excellent PSF (2.5 arc-seconds half-power diameter), large effective area (1,800 cm2 at 1 keV), and large field of view (1 deg2), making it more than an order of magnitude more capable and more sensitive than Chandra and Swift/XRT to conduct survey and to find and study transient sources. The UVT has a 30-cm aperture with 5 filters covering the 160nm to 340nm band, providing simultaneous spectral coverage with the XRT and enabling reverberation mapping of accretion disk geometries around black holes. The SC is highly autonomous and is capable of fast slewing, enabling efficient raster scans and time-domain surveys. In particular, in combination with a state-of-the-art mission operations center at the University of Colorado, the SC can respond to targets of opportunity within 2 hours 90% of the time. With its nearly equatorial low-earth orbit, STAR-X’s telescopes will have low particle background, enabling them to have unprecedented sensitivity for measuring faint diffuse emissions from clusters of galaxies. STAR-X is a timely response to Astro2020’s recommendation for a space-based, sustaining time-domain and multi-messenger program.

William W. Zhang↗

The Second Radio Synchrotron Background Workshop: Conference Summary and Report

We summarize the second radio synchrotron background workshop, which took place on 2022 June 15–17 in Barolo, Italy. This meeting was convened because available measurements of the diffuse radio zero level continue to suggest that it is several times higher than can be attributed to known Galactic and extragalactic sources and processes, rendering it the least well-understood electromagnetic background at present and a major outstanding question in astrophysics. The workshop agreed on the next priorities for investigations of this phenomenon, which include searching for evidence of the radio Sunyaev–Zel’dovich effect, carrying out cross-correlation analyses of radio emission with other tracers, and supporting the completion of the 310 MHz absolutely calibrated sky map project.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗

Satellite-derived surface energy balance estimates in the Alaskan sub-Arctic

The capabilities of the heat capacity mapping mission (HCMM) satellite are described, together with the numerical model used for data analysis. The HCMM carries the HCM radiometer for remote sensing in the visible and 0.55-1.1 micron wavelength regions, as well as in the 10.5-12.5 micron interval, with the swath being 720 km and each pixel being a square 500 m on a side. The HCMM is intended to aid in hydrological studies of soil moisture, runoff, and evapotranspiration estimates. Data are taken of the albedos and temperatures of vegetation, assuming that all reflection is diffuse. Corrections are made in the algorithm to account for sun angle and the spacecraft distance from the earth. Sample calculations are provided from scans of a coastal plain, mountains, and plateaus of Alaska.

Gurney, R. J.↗

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗