SEARCH · Search NASA
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.
The conclusions of the application of first principles model to spacecraft operations are: the first principles of Bi-phasic electrode presented model provides an explanation for many behaviors on voltage fading on LEO cycling.
We examine the accuracy of density prediction by the first principles model Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) developed by the National Center for Atmospheric Research and compare it to the accuracy of three empirical models: Jacchia 71, the Naval Research Laboratory Mass Spectrometer Incoherent Scatter Extended 2000 (NRLMSIS), Jacchia 1971, and Jacchia-Bowman 2008. Comparisons are made for three large storms: the October 2003 storm, the March 2013 storm, and the March 2015 storm. To evaluate the accuracy of these models we use tracking data for nine space objects in low Earth orbit. Additionally, we evaluate the accuracy of the TIEGCM and NRLMSIS with data from high precision accelerometers on the Challenging Minisatellite Payload (CHAMP) and Gravity field and Circulation Explorer (GOCE) satellites. The goal is to assess the use of a first principles model as a potential tool for forecasting satellite drag during large magnetic storms. For the storms considered, we found the TIEGCM, JB2008, and NRLMSIS models to be substantially more accurate than the Jacchia 71 model. The accuracies of the TIEGCM and JB2008 models were similar, but overall, the TIEGCM was more accurate. We found smaller differences for TIEGCM versus CHAMP than for NRLMIS for the Halloween Storm, and smaller differences than results published for JB2008 and the assimilative model HASDM. The empirical models are at present more practical for operational purposes, but the TIEGCM, developed as a research model, with a greater focus on operational use offers the potential for improved utility during stressing conditions.
Completeness, efficiency and autonomy are requirements for suture diagnostic reasoning systems. Methods for automating diagnostic reasoning systems include diagnosis from first principles (i.e., reasoning from a thorough description of structure and behavior) and diagnosis from experiential knowledge (i.e., reasoning from a set of examples obtained from experts). However, implementation of either as a single reasoning method fails to meet these requirements. The approach of combining reasoning from first principles and reasoning from experiential knowledge does address the requirements discussed above and can possibly ease some of the difficulties associated with knowledge acquisition by allowing developers to systematically enumerate a portion of the knowledge necessary to build the diagnosis program. The ability to enumerate knowledge systematically facilitates defining the program's scope, completeness, and competence and assists in bounding, controlling, and guiding the knowledge acquisition process.
Recent rapid progresses in physics theory and computational power have made it possible to predict the martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles [1-3]. In particular, rigorous while time-consuming thermodynamic integration has been employed to compute the anharmonic phonon free energies, which play a crucial role in determining martensitic phase transitions in SMAs. However, this approach has only been applied to simple binaries, and its accuracy is unsatisfying for certain SMAs such as the most commonly used NiTi. In this work, we report on several new developments to our method that bring first-principles theory and experiment much closer into agreement including the MTT of NiTi, and that improve the computational efficiency significantly. We have applied our refined approach to investigate the Ni0.5Ti0.5-xHfx and PdxNi0.5-xTi0.5 ternaries, and the predicted MTT for each composition is within 100K compared with experiment. We will address various techniques to overcome the difficulty encountered in studying ternaries. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties. [1] J.B. Haskins, A.E. Thompson,and J.W. Lawson, Phys. Rev B 94, 214110 (2016). [2] J.B. HaskinsandJ.W. Lawson, J. App. Phys. 121, 205103 (2017). [3] J.B. Haskins, H. Malmir, S. J. Honrao, L. A. Sandoval, and J.W. Lawson, Acta Materialia 212, 116872 (2017).
Recent rapid progresses in physics theory and computational power have made it possible to predict the martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles [1-3]. In particular, rigorous while time-consuming thermodynamic integration has been employed to compute the anharmonic phonon free energies, which play a crucial role in determining martensitic phase transitions in SMAs. However, this approach has only been applied to simple binaries, and its accuracy is unsatisfying for certain SMAs such as the most commonly used NiTi. In this work, we report on several new developments to our method that bring first-principles theory and experiment much closer into agreement including the MTT of NiTi, and that improve the computational efficiency significantly. We have applied our refined approach to investigate the Ni0.5Ti0.5-xHfx and PdxNi0.5-xTi0.5 ternaries, and the predicted MTT for each composition is within 100K compared with experiment. We will address various techniques to overcome the difficulty encountered in studying ternaries. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties.
Recent rapid progresses in physics theory and computational power have made it possible to accurately predict the phase transitions and martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles. However, previously theory and calculations [1-4] were applied only to study highly ordered stoichiometric binary alloys such as NiTi, PdTi and NiHf. Here we report on our recent first-principles investigations [5,6] on Ni 0.5 Ti 0.5-x Hf x and Pd x Ni 0.5-x Ti 0.5 ternaries and off-stoichiometric NiTi, and the predicted martensitic phase transitions in these complex SMAs are in good agreement with experimental findings. In particular, the calculated MTTs for all these compositions are within 100 K compared with the corresponding measured data, and our results also reveal the origin of the striking asymmetry in MTT of the off-stoichiometric NiTi near equiatomic compositions. We will address various techniques to overcome the difficulty encountered in studying ternaries and off-stoichiometic binaries associated with disorder and/or much lowered symmetry. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties. References: [1] J. B. Haskins, A. E. Thompson, and J. W. Lawson, Phys. Rev B 94 , 214110 (2016). [2] J. B. Haskins and J. W. Lawson, J. App. Phys. 121 , 205103 (2017). [3] J. B. Haskins, H. Malmir, S. J. Honrao, L. A. Sandoval, and J. W. Lawson, Acta Materialia 212 , 116872 (2017). [4] Z. Wu, J. W. Lawson, and O. Benafan, Phys. Rev. B 106 , L140102 (2022). [5] Z. Wu, H. Malmir, O. Benafan, and J. W. Lawson, Acta Materialia 261 , 119362 (2023). [6] Z. Wu, J. W. Lawson, and O. Benafan, Phys. Rev. B 108 , L140103 (2023).
Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.
A discrete random medium is an object in the form of a finite volume of a vacuum or a homogeneous material medium filled with quasi-randomly and quasi-uniformly distributed discrete macroscopic impurities called small particles. Such objects are ubiquitous in natural and artificial environments. They are often characterized by analyzing theoretically the results of laboratory, in situ, or remote-sensing measurements of the scattering of light and other electromagnetic radiation. Electromagnetic scattering and absorption by particles can also affect the energy budget of a discrete random medium and hence various ambient physical and chemical processes. In either case electromagnetic scattering must be modeled in terms of appropriate optical observables, i.e., quadratic or bilinear forms in the field that quantify the reading of a relevant optical instrument or the electromagnetic energy budget. It is generally believed that time-harmonic Maxwell's equations can accurately describe elastic electromagnetic scattering by macroscopic particulate media that change in time much more slowly than the incident electromagnetic field. However, direct solutions of these equations for discrete random media had been impracticable until quite recently. This has led to a widespread use of various phenomenological approaches in situations when their very applicability can be questioned. Recently, however, a new branch of physical optics has emerged wherein electromagnetic scattering by discrete and discretely heterogeneous random media is modeled directly by using analytical or numerically exact computer solutions of the Maxwell equations. Therefore, the main objective of this Report is to formulate the general theoretical framework of electromagnetic scattering by discrete random media rooted in the Maxwell- Lorentz electromagnetics and discuss its immediate analytical and numerical consequences. Starting from the microscopic Maxwell-Lorentz equations, we trace the development of the first principles formalism enabling accurate calculations of monochromatic and quasi-monochromatic scattering by static and randomly varying multiparticle groups. We illustrate how this general framework can be coupled with state-of-the-art computer solvers of the Maxwell equations and applied to direct modeling of electromagnetic scattering by representative random multi-particle groups with arbitrary packing densities. This first-principles modeling yields general physical insights unavailable with phenomenological approaches. We discuss how the first-order-scattering approximation, the radiative transfer theory, and the theory of weak localization of electromagnetic waves can be derived as immediate corollaries of the Maxwell equations for very specific and well-defined kinds of particulate medium. These recent developments confirm the mesoscopic origin of the radiative transfer, weak localization, and effective-medium regimes and help evaluate the numerical accuracy of widely used approximate modeling methodologies.
The paper presents a first-principles study of the shear modulus tensor for perfect and imperfect Coulomb solids. Allowance is made for the effects of thermal fluctuations for temperatures up to the melting conditions. The present theory treats the cases of the long-range Coulomb interaction, where volume fluctuations should be avoided in the Ewald sums.
Several battery models are under development at JPL based on first principles. The recent models are based on NiH2 and NiMH chemistries. Performance results and computation requirements are discussed.
Li-O2 batteries have traditionally used carbon based electrodes (graphite, buckypaper) as the cathode of choice due to its good electrical conductivity, stability against non-aqueous electrolytes like Dimethyl ether (DME) and ease of handling. But, the carbon cathode also leads to formation of carbonate by-products that increase overpotentials during charging leading to degradation of cathode and reduction of cyclability. In this work, we investigate some of the well-known oxides as cathodes with focus on the interface between the oxide surfaces and the discharge product: Li2O2, in the Li-O2 battery using first principles computations. Our results show that attention must be paid on choosing the appropriate surface of the oxides. We extend the analysis to suggest other possible oxide chemistries that should be investigated as cathodes in Li-O2 batteries.
A theoretical model describing orificed hollow cathode discharge is presented. The approach adopted is based on a purely analytical formulation founded on first principles. The present model predicts the emission surface temperature and plasma properties such as electron temperature, number densities and plasma potential. In general, good agreements between theory and experiment are obtained. Comparison of the results with the available related experimental data shows a maximum difference of 10 percent in emission surface temperature, 20 percent in electron temperature and 35 percent in plasma potential. In case of the variation of the electron number density with the discharge current a maximum discrepancy of 36 percent is obtained. However, in the case of the variation with the cathode internal pressure, the predicted electron number density is higher than the experimental data by a maximum factor of 2.
The imperative to quantify the Earths electromagnetic-energy budget with an extremely high accuracy has been widely recognized but has never been formulated in the framework of fundamental physics. In this paper we give a first-principles definition of the planetary electromagnetic-energy budget using the Poynting- vector formalism and discuss how it can, in principle, be measured. Our derivation is based on an absolute minimum of theoretical assumptions, is free of outdated notions of phenomenological radiometry, and naturally leads to the conceptual formulation of an instrument called the double hemispherical cavity radiometer (DHCR). The practical measurement of the planetary energy budget would require flying a constellation of several dozen planet-orbiting satellites hosting identical well-calibrated DHCRs.
Development of shape memory alloys (SMA) has been traditionally accomplished by means of extensive empirical efforts, with a very limited support from phenomenological models and related tools. We present a theoretical framework based on the Hamiltonian formalism, which connects important thermodynamic properties directly to the relative energies (the eigenvalues of the Hamiltonian), computed from the first principles using density functional theory. The estimates based on this formalism allow to establish correlations between composition and physical properties, such as relative energies, phase transition temperatures, hysteresis, thermal expansion, and other thermodynamic information. Importantly, estimates and correlations not only provide a fast answer, but also elucidate underpinning of thermodynamic features in terms of the electronic structure. Using this framework, we consider phase transformations and correlations between the selected properties and composition in the NiTi-based ternary and multicomponent shape memory alloys. Our theoretical guidance facilitates design and development of future alloys.
The SUNDIAL program offers a unique opportunity to study ionospheric behavior on a global scale. As part of this program, data pertaining to solar, interplanetary, magnetospheric, ionospheric, and thermospheric conditions are collected simultaneously from a large number of satellite and ground-based sites spread around the world. In the coming years, these data should lead to a major improvement in both empirical and first-principle ionospheric models. As a benchmark against which to compare future progress, the present state of empirical and numerical ionospheric modeling is discussed. The discussion covers the capabilities and limitations of the existing models as well as the direction of future modeling efforts.
Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.
The rate of magnetic reconnection is of the utmost importance in a variety of processes because it controls, for example, the rate energy is released in solar flares, the speed of the Dungey convection cycle in Earth’s magnetosphere, and the energy release rate in harmful geomagnetic substorms. It is known from numerical simulations and satellite observations that the rate is approximately 0.1 in normalized units, but despite years of effort, a full theoretical prediction has not been obtained. Here, we present a first-principles theory for the reconnection rate in non-relativistic electron-ion collisionless plasmas, and show that the same prediction explains why Sweet-Parker reconnection is considerably slower. The key consideration of this analysis is the pressure at the reconnection site (i.e., the x-line). We show that the Hall electromagnetic fields in antiparallel reconnection cause an energy void, equivalently a pressure depletion, at the x-line, so the reconnection exhaust opens out, enabling the fast rate of 0.1. If the energy can reach the x-line to replenish the pressure, the exhaust does not open out. In addition to heliospheric applications, these results are expected to impact reconnection studies in planetary magnetospheres, magnetically confined fusion devices, and astrophysical plasmas.
We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.