Search NASASearch

SEARCH · Search NASA

Results for “Neural ordinary differential equations”

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.

Usage-based Lifing of Lithium-Ion Battery with HybridPhysics-Informed Neural Networks

Lithium-ion batteries are commonly used to power unmanned aircraft vehicles (UAVs).The ability to model and forecast the remaining useful life of these batteries enables UAV reliability assurance. Building accurate models for battery state of charge and state of health based on first principles is challenging due to the complex electrochemistry that governs battery operations and computational complexity required to solve them. Therefore, reduced order models are often used due to their ability to capture the overall battery discharge. Un-fortunately, these simplifications lead to residual discrepancy between model predictions and observed data. In this paper, we present a hybrid modeling approach merging reduced-order models and neural networks. In this approach, while most of the input-output relationship is captured by Nernst and Butler-Volmer equations, data-driven kernels reduce the gap between predictions and observations. We validate our approach using data publicly available through the NASA Prognostics Center of Excellence repository. Results showed that our hybrid battery prognosis model can be successfully calibrated, even with a limited number of observations.

Lithium-ion Battery

Li-ion Battery Aging with Hybrid Physics-Informed Neural Networks and Fleet-wide Data

In this work, we propose a hybrid model for Li-ion battery discharge and aging prediction that leverages fleet-wide data to predict future capacity drops.The model is built upon an hybrid approach merging physics-based and empirical equations, as well as neural network models in a recurrent neural network cell. The hybrid physics-informed neural network can predict voltage discharge cycles given the loading profile, and estimate the used capacity of the battery under random-loading conditions by tracking aging parameters connected to the residual capacity of the battery. By merging information on the battery aging parameters with existing fleet-wide aging data, the model can predict the future residual capacity of the battery that is being monitored, and therefore enable predictions of voltage discharge curves far ahead in the battery life cycle. We validated the approach using the NASA Prognostics Data Repository Battery data-set, which contains experimental data on Li-ion batteries discharged at random loading conditions in a controlled environment. The approach also allows the identification of discrepancies between the battery aging trend and the trend observed at the fleet level, so that batteries behaving differently from the rest of the fleet can be subject to closer monitoring and further testing to refine predictions.

PINN

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya

Extensional Flow Convecting a Reactant Undergoing a First Order Homogeneous Reaction and Diffusional Mass Transfer From a Sphere at Low to Intermediate Peclet and Damkohler Numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 ≤ Pe ≤ 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 ≤ Da ≤ 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r,Θ) = Σ c n (r)P n (cosΘ) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions P n (cosΘ), the variable coefficient PDE can be reduced to a system of N+1 second order ODEs for the radial functions C n (r), n=0,1,2, ... ,N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N+ 1) first order differential equations as a two-point boundary value problem on 1 ≤ r ≤ r •• . Convergence of the expansion for a specific value of N can thus be established and provides "spectral" behavior as well as the full concentration field c(r,Θ).

N Y Shah

Analysis of Gas Absorption to A Thin Liquid Film in the Presence of A Zero-Order Chemical Reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

S Rajagopalan

Autonomous Detection and Classification of Lunar Minerals Using a Convolutional Neural Network Based Framework for the SUCR DALI Project

NASA’s long-term goal is to deploy humans to the Moon and, from there, advance human exploration to Mars, with Artemis missions as pivotal milestones. Raman spectroscopy can uniquely identify minerals, compounds, water states, and other materials, providing distinctive fingerprints for classification. A Raman instrument has been successfully deployed and utilized on the Mars surface via the Perseverance rover, but has not yet been utilized at the lunar surface The SUCR DALI project is working towards developing a Raman spectroscopy instrument to be applied in various lunar mission concepts, including within the Artemis program. The objective of my research is to assist in the maturation of the proposed SUCR DALI lunar Raman instrument through the development of an autonomous detection and classification model capable of identifying minerals and water states on the Moon’s surface.

Convolutional Neural Networks

Differential Drag Efficacy for Close Approach Remediation

Differential drag has become a viable alternative to propulsion for satellites to avoid collisions, but there is little guidance in the literature to aid mission designers in developing a differential drag capability that verifiably meets collision avoidance efficacy standards or requirements, if such requirements were to exist. This paper proposes a differential drag efficacy determination approach based on empirical conjunctions from the NASA Conjunction Assessment Risk Analysis historical database, focusing on energy dissipation rate and change in ballistic coefficient as the key satellite parameters correlated to efficacy. The data analysis informs the discussion toward adoption of recommended differential drag requirements. A case study is presented to walk through the process to determine efficacy of a proposed mission assuming several potential requirements.

conjunction remediation

Geochemical Zoning and Early Differentiation in the Moon

The volatile elements (e.g., Rb, Pb, Tl, Bi, Cs) seem to have been depleted at the time of lunar accretion. Accordingly, it may be assumed that the moon initially accreted from refractory material. The good correlation between volatile/involatile element ratios (e.g., Cs/U, K/La, K/Zr) in both highland and maria samples means that element distribution in lunar crustal rocks is not governed by volatility differences. This and other evidence encourages the view that the moon was accreted homogeneously. A consequence of homogeneous accretion theories is that very efficient large-scale element fractionation is required to account both for the high near-surface concentrations of refractory elements (e.g., Th, U, REE, Zr, Ba, etc.) and for the Ca-Al-rich crust.

S R Taylor

High Temperature Dielectric Properties and Differential Scanning Calorimetry of Lunar Simulants

To guide development of microwave process technology that could be used during in situ construction on the Moon, we measured the high-temperature basic dielectric properties (εʹ and εʺ) of 17 lunar simulants and related materials. In order to confidently use these data one needs to understand the data’s strengths and weaknesses. Therefore, a goal of this publication is to provide insights into the comparative effects of sample composition, pre-treatments, experimental variables, high temperatures, and other factors on the measured response. The dielectric measurements were performed using the cavity perturbation method over a temperature range between room temperature to 1000 °C, or higher, and provided the real and imaginary components of permittivity at six frequencies. The utility of the original values was limited by the varying density of the pellets used in the measurement. Therefore, all of the εʹ and εʺ measurements at the frequency of 2466 MHz have been scaled to a constant density, 1.75 g/cm 3 . Here the data are presented as graphs chosen to aid analysis within and across simulant groups. To gain additional insight into the processes happening at the elevated temperatures in the dielectric measurements, heat capacity data was obtained using differential scanning calorimetry (DSC) on several of the simulant materials. Our data show that over the frequency range 397 MHz – 2985 MHz a material’s behavior does not greatly change, as compared to the scale of differences observed between lunar mare and highland simulants at high temperatures. For example at 1000 °C, the mare simulant JSC-1A absorbs 10 times more power than the highland simulant NUW-LHT-5M. We observe that as melting temperatures are reached both permittivity and dielectric loss rise non-linearly, helping to explain thermal runaway during microware heating. Our data show that even less than a few weight % of many non-lunar minerals, and the use of mixtures in simulants can affect the dielectric behavior at higher temperatures. A comparison of our results with published dielectric data for Apollo samples and with remote sensing of the Moon supports the conclusion the simulants and lunar material at room temperature have very similar dielectric values.

Differential Scanning Calorimetry

The Effect of Temperature and Pressure on the Distribution of Iron Group Elements Between Metal and Olivine Phases in the Process of Differentiation of Protoplanetary Material

The distribution patterns of Ni, Co, Mn, and Cr were studied in olivines of various origins: from meteorites (chondrites, achondrites, pallasites), which are likely analogs of the protoplanetary material, to peridotite inclusions in kimberlite pipes, which are analogs of mantle material. According to X-ray microanalysis data, each genetic group of olivines is characterized by a specific concentration of these elements. Nickel is concentrated (up to 0.34 percent) in peridotite olivines, while manganese is concentrated in meteoritic olivines. The maximum chromium content (0.2 percent) was found in ureilites, which were formed under reducing conditions. Experiments at pressures of 20 to 70 kbar and temperatures of 1100 to 2000°C have shown that in a mixture of olivine and Ni metal or NiO nickel enters the silicate phase (up to 4 percent), displacing Fe into the metallic phase. Equilibrium temperatures were estimated from the Fe, Ni distribution coefficients between the metal and olivine: 1500 K for pallasites, 1600 K for olivine-bronzite H6 chondrites, 1200 K for olivine-hypersthene L6, 900 K for LL6, and 1900 K for ureilites (at P = 1 atm). The equilibrium conditions of peridotites are close to T = 1800 K and P over 100 kbar. The distribution patterns of the transition elements are explained on the basis of physical-chemical properties. It is concluded that there is a sharp difference between the conditions of differentiation of the protoplanetary material at the time meteorites were formed and the conditions of differentiation of the planets into concentric layers.

A P Vinogradov

Costas Loop Demodulation of Suppressed Carrier BPSK Signals in the Dsn Environment: Experimental Results Obtained at TDL

Suppressed carrier binary phase-shift keyed (BPSK) signalling is currently being considered as a design alternative for future DSN telemetry in the multimegabit range. Carrier tracking of such signals is usually achieved by a Costas loop, as opposed to the ordinary phase lock loop. A Costas loop capable of demodulating BPSK signals with data rates up to 1 Msps was designed and constructed and its Doppler tracking performance with respect to a Block 3 receiver was tested at the Telecommunications Development Laboratory (TDL). The compatibility of suppressed carrier signalling with the current radiometric system, specifically Doppler tracking and ranging, was investigated. The experimental results obtained to-date with respect to Doppler tracking are presented.

R Reasoner

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG

Characterization of Ternary NiTiPt High-Temperature Shape Memory Alloys

Pt additions substituted for Ni in NiTi alloys are known to increase the transformation temperature of the alloy but only at fairly high Pt levels. However, until now only ternary compositions with a very specific stoichiometry, Ni 50-x Pt x Ti 50 , have been investigated and then only to very limited extent. In order to learn more about this potential high-temperature shape memory alloy system, a series of over twenty alloys along and on either side of a line of constant stoichiometry between NiTi and TiPt were arc melted, homogenized, and characterized in terms of their microstructure, transformation temperatures, and hardness. The resulting microstructures were examined by scanning electron microscopy and the phase compositions quantified by energy dispersive spectroscopy. "Stoichiometric" compositions along a line of constant stoichiometry between NiTi to TiPt were essentially single phase but any deviations from a stoichiometry of (Ni,Pt) 50 Ti 50 resulted in the presence of at least two different intermetallic phases, depending on the overall composition of the alloy. Essentially all alloys, whether single or two-phase, still under went a martensitic transformation. It was found that the transformation temperatures were depressed with initial Pt additions but at levels greater than 10 at.% the transformation temperature increased linearly with Pt content. Also, the transformation temperatures were relatively insensitive to alloy stoichiometry within the range of alloys examined. Finally, the dependence of hardness on Pt content for a series of Ni 50-x Pt x Ti 50 alloys showed solution softening at low Pt levels, while hardening was observed in ternary alloys containing more than about 10 at.% Pt. On either side of these "stoichiometric" compositions, hardness was also found to increase significantly.

Transformation Temperatures

Petrogenesis of Lunar Rocks: Rb-Sr Constraints and Lack of H2O

Rb and Sr isotopic data and other chemical data indicate major lunar differentiation at about 4.6 AE (AE = 10 9 years) and very limited subsequent differentiation. The constraints of limited differentiation after 4.6 AE and the apparent lack of H20 on the Moon, when applied to the derivation and petrogenesis of lunar samples, suggest the following: (1) soil samples, breccias, metaclastic rocks, and feldspathic basalts represent mixtures of repeatedly modified clastic material, which was ultimately derived from materials formed during the ~ 4.6 AE differentiation; and (2) mare basalts crystallized from melts which formed by partial melting, and which developed without equilibration between the melt and crystalline residuum.

Arden L Albee

Experimental Demonstration of a Supersonic Spike Ejector

Spike diffusers are presented as a means to improve supersonic ejector efficiency via reduction in shock losses. Design, fabrication, and testing of an articulated prototype is described. Experimental results are compared to empirical performance curves representing the loci of optimal one- and two-stage conventional ejectors driven by dry air and discharging to atmosphere. Operated at a compression ratio of ~10 with all inputs equal, the spike ejector is shown to have achieved up to 2.19X and 1.83X the secondary mass throughput of the one- and two-stage references, respectively. This suggests the potential of spike diffusers to reduce ejector operating costs by ~50%.

wind tunnel

Riddles About the Origin and Thermal History of the Moon

Analyses of lunar rocks have confirmed that the interior of the Moon has been subjected to magmatic differentiation, as has been determined by calculation of the thermal history of the Moon. However, it appears that differentiation occurred so early that it could not have been related to heating by long-lived radioactive elements, but rather indicates a high initial temperature of the Moon. The reasons for the high initial temperature remain unknown. The presence of a central, partially melted area revealed by seismic observations forces us to reject convective models of the Moon, while the dimensions of this area, together with the measured heat fluxes from the interior, serve as "boundary conditions" for any new calculations of the thermal history for a two-layered differentiated model of the Moon operative over the past 3.5 billion years. If we start with the early differentiation of the Moon and the high content of long-lived radioactive elements in it, explanation of its current thermal properties represents no difficulty. But these prerequisites themselves remain unexplained. Due to the continuing controversy concerning the true nature of the process of magmatic differentiation of the Moon, estimation of the composition of the entire Moon on the basis of analyses of surface rock remains unreliable, and use of these estimates to clarify the origin of the Moon is premature. Existing hypotheses concerning the origin of the Moon encounter dynamic difficulties, as well as difficulties with explanation of the high initial temperature, to say nothing of the difficulty involved in explaining the proposed differences in the composition of the Earth and Moon.

B Yu Levin