Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculations”

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.

298 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRS1983

Battery model for electrical power system energy balance

A model to simulate nickel-cadmium battery performance and response in a spacecraft electrical power system energy balance calculation was developed. The voltage of the battery is given as a function of temperature, operating depth-of-charge (DOD), and battery state-of-charge. Also accounted for is charge inefficiency. A battery is modeled by analysis of the results of a multiparameter battery cycling test at various temperatures and DOD's.

Hafen, D. P.↗
NASA NTRS2025

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path↗
NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

DKIN/Stiffener Debonding↗
NASA NTRS2003

Analysis of Surface and Bulk Behavior in Ni-Pd Alloys

The most salient features of the surface structure and bulk behavior of Ni-Pd alloys have been studied using the BFS method for alloys. Large-scale atomistic simulations were performed to investigate surface segregation profiles as a function of temperature, crystal face, and composition. Pd enrichment of the first layer was observed in (111) and (100) surfaces, and enrichment of the top two layers occurred for (110) surfaces. In all cases, the segregation profile shows alternate planes enriched and depleted in Pd. In addition, the phase structure of bulk Ni-Pd alloys as a function of temperature and composition was studied. A weak ordering tendency was observed at low temperatures, which helps explain the compositional oscillations in the segregation profiles. Finally, based on atom-by-atom static energy calculations, a comprehensive explanation for the observed surface and bulk features will be presented in terms of competing chemical and strain energy effects.

Ordering↗
NASA NTRSDate not supplied

Swamp Works Regolith Compaction Technologies

While the level of compaction below the lunar surface increases quickly after only a few cm of depth, in many cases during a construction mission there will be a need to excavate and transport regolith to a new location for cut-and-fill or horizontal construction of structures such as berms. In these cases to achieve high levels of bulk density, compaction must be per-formed. Additionally, in some cases surface technologies such as systems that sinter/melt the surface may desire the maximum possible compaction at the sur-face to improve melting/heating performance and the final material strength properties. Kennedy Space Center’s (KSC) Swamp Works has developed two means of compaction, lunar and mar-tian compaction. Planetary Autonomous Compaction Technology (PACT) which is part of the Multifunction End Effector for Regolith Compaction Acquisition and Transfer (MEERCAT) robotic arm end effector system’s capabilities and the Site Preparation Tooling for Operations on Mobility Platforms (STOMP) vibratory roller compactor. PACT on MEERCAT has been demonstrated to a TRL 5 and STOMP to a TRL 4 in ambient testing. The results of PACT on MEERCAT and STOMP testing will be shared with results for various simulants including BP-1, ICN-LHT-1G (aka CSM-LHT-1G), RDW-LHT-1GH (a simulant developed for the Mason Tipping Point to match characteristics of ICN-LHT-1G), and Exolith LHS-1E. This will also include discussions on methods used to verify relative density before and after compaction and means to verify density effects below depth. To calculate relative density, maximum and minimum densities for simulants were taken from literature and additional lab testing (publication in work).

redwire↗
NASA NTRS2026

Swamp Works Regolith Compaction Technologies

While the level of compaction below the lunar surface increases quickly after only a few cm of depth, in many cases during a construction mission there will be a need to excavate and transport regolith to a new location for cut-and-fill or horizontal construction of structures such as berms. In these cases to achieve high levels of bulk density, compaction must be per-formed. Additionally, in some cases surface technolo-gies such as systems that sinter/melt the surface may desire the maximum possible compaction at the sur-face to improve melting/heating performance and the final material strength properties. Kennedy Space Center’s (KSC) Swamp Works has developed two means of compaction, lunar and mar-tian compaction. Planetary Autonomous Compaction Technology (PACT) which is part of the Multifunction End Effector for Regolith Compaction Acquisition and Transfer (MEERCAT [1]) robotic arm end effector system’s capabilities and the Site Preparation Tooling for Operations on Mobility Platforms (STOMP [2]) vibratory roller compactor. PACT on MEERCAT has been demonstrated to a TRL 5 and STOMP to a TRL 4 in ambient testing. The results of PACT on MEERCAT and STOMP testing will be shared with results for various simulants including BP-1, ICN-LHT-1G (aka CSM-LHT-1G), RDW-LHT-1GH (a simulant developed for the Mason Tipping Point to match characteristics of ICN-LHT-1G), and Exolith LHS-1E. This will also include discus-sions on methods used to verify relative density before and after compaction and means to verify density effects below depth. To calculate relative density, maximum and minimum densities for simulants were taken from literature [3] [4] [5] and additional lab test-ing (publication in work).

vibration↗
NASA NTRS1983

Modeling a Li/SOCl2 battery for design purposes

A generalized code applicable to many different electrochemical systems and geometric designs is discussed. The code is to be set up so that physical property data such as thermal conductivity, viscosity, density, and configuration (e.g. physical dimensions) are the input data. Thus, by changing these parameters many different battery configurations can be handled. The outputs, as a function of time and space, are voltage, current, temperature, pressure, velocity, and species concentration.

Ernst, D. W.↗
NASA NTRS1993

Experimental and analytical ion thruster research

The results of further spectroscopic studies on the plume from a 3 cm ion source operated on an argon propellant is reported on. In particular, it is shown that it should be possible to use the spectroscopic technique to measure the plasma density of the ion plume close to the grids, where it is difficult to use electrical probe measurements. How the technique, along with electrical probe measurements in the far downstream region of the plume, can be used to characterize the operation of a three-grid, 15 cm diameter thruster from NASA JPL is outlined. Pumping speed measurements on the Vacuum Research Facility have shown that this facility should be adequate for testing the JPL thruster at pressures in the low 10(exp -5) Torr range. Finally, we describe a simple analytical model which can be used to calculate the grid impingement current which results from charge-exchange collisions in the ion plume.

Wilhelmus M. Ruyten↗
NASA NTRS2026

Recommended Practices in Laser-Induced Fluorescence (LIF) Diagnostics for Electric Propulsion

Laser-induced fluorescence (LIF) spectroscopy uses Doppler-shifted laser photons to excite a bound electron transition in an atom or ion, with subsequent fluorescence emission detected in order to measure the local velocity distribution function (VDF) and/or particle density. Due to its non-invasive nature, high spatial resolution, and acceptable difficulty of implementation, LIF has become a favored diagnostic for electric propulsion (EP) plasmas, particularly in the study of Hall thrusters. This paper presents recommended best practices for LIF measurements in electric thrusters, drawn from over three decades of implementation heritage in the EP community. Focusing on both single-photon LIF (typically used to measure VDFs) and two-photon LIF (TALIF, typically used to determine densities), the paper covers selection of atomic transitions to target, experimental setups, and interpretation of data. A number of special topics and advanced applications are discussed, including uncertainty analysis, techniques for high-speed measurements, and applications to a variety of propellants.

plasma plume↗
NASA NTRS1977

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↗
NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗
NASA NTRSDate not supplied

Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)

The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.

Adam Johnson↗
NASA NTRS1991

The Real-Time Control of Planetary Rovers Through Behavior Modification

It is not yet clear of what type, and how much, intelligence is needed for a planetary rover to function semi-autonomously on a planetary surface. Current designs assume an advanced AI system that maintains a detailed map of its journeys and the surroundings, and that carefully calculates and tests every move in advance. To achieve these abilities, and because of the limitations of space-qualified electronics, the supporting rover is quite sizable, massing a large fraction of a ton, and requiring technology advances in everything from power to ground operations. An alternative approach is to use a behavior driven control scheme. Recent research has shown that many complex tasks may be achieved by programming a robot with a set of behaviors and activation or deactivating a subset of those behaviors as required by the specific situation in which the robot finds itself. Behavior control requires much less computation than is required by tradition AI planning techniques. The reduced computation requirements allows the entire rover to be scaled down as appropriate (only down-link communications and payload do not scale under these circumstances). The missions that can be handled by the real-time control and operation of a set of small, semi-autonomous, interacting, behavior-controlled planetary rovers are discussed.

David P Miller↗
NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗