Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemistry - Inorganic”

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.

416 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗
DOE OSTI2026

Modern insights into the mechanisms of neptunium oxalate decomposition

Neptunium oxalate (Np(C 2 O 4 ) 2 ·6H 2 O) is a historically relevant solid phase used in nuclear processing as a precursor for neptunium dioxide (NpO 2 ). Although Np oxalate has been synthesized and used for NpO 2 production for decades, the thermal decomposition mechanism of this phase remains poorly understood and has not been evaluated in over 30 years. Conflicting reports in historical literature suggest either a direct conversion from anhydrous oxalate to NpO 2 or a decomposition that includes the formation of Np carbonate or oxidized Np intermediate phases. In this work, we reexamine the decomposition pathway of Np(C 2 O 4 ) 2 ·6H 2 O using thermal analysis coupled with evolved gas analysis and temperature-dependent Raman spectroscopy to elucidate decomposition mechanisms and intermediate phases using modern analytical techniques. Thermal analysis revealed a three-stage decomposition process, including dehydration below 200 °C, oxalate breakdown between 170 and 370 °C, and NpO 2 formation by 500 °C. However, an unidentified plateau in the thermal data was observed during measurements. Raman spectroscopy confirmed the stages of decomposition, and in the analysis of potential intermediate phases, no carbonate phases or Np 2 O 5 were identified. Raman data suggest that residual oxalate or nonstoichiometric oxide are present during decomposition before pure NpO 2 is formed. These findings clarify aspects of the Np oxalate decomposition mechanism and address longstanding discrepancies in the literature, with a specific focus on Np-specific materials chemistry.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗
DOE OSTI2026

Liquified SO 2 induced solid/cathode electrolyte interphase for lithium ion batteries

Formation of robust solid/cathode electrolyte interphases (S/CEI) is vital for long-term stability and high-performance operation of lithium-ion batteries (LIBs), particularly under high voltage regimes. However, engineering electrochemically stable S/CEIs that effectively suppress interfacial side reactions remains a key challenge. Herein, we introduce a liquefied sulfur dioxide (SO 2 )– ionic liquid complex as a fluorine-free multifunctional electrolyte additive for the first time that significantly improves the formation of sulfate/sulfite-rich S/CEI layers at both graphite and NMC811 interfaces. The unique SO 2 -N coordination with a 1,2,4-triazolide-based ionic liquid enables homogeneous SO 2 dissolution, resulting in controlled SO 2 decomposition during the initial electrochemical cycle. This decomposition yields sulfur-rich interphase species that stabilize the electrolyte-electrode interface, reduce impedance growth, and lessen electrolyte decomposition. Electrochemical tests show significantly improved cycle life, reduced polarization, and increased Coulombic efficiency for both anodes and cathodes. XPS confirms the presence of SO 2 -derived surface species that contribute to interfacial stability. In conclusion, this approach highlights a new direction for interphase engineering using liquefied gas additives and opens pathways for sulfur-based S/CEI chemistry in advanced battery systems.

Graphite↗
NASA NTRS2026

Coma Physics of an Interstellar Object: JWST Spatial-Spectral Mapping of 3I/ATLAS

We report a survey of molecular emission from cometary volatiles using the James Webb Space Telescope (JWST) toward interstellar object 3I/ATLAS carried out on UT 2025 December 22 and 23 at a heliocentric distance (\rh{}) of $2.37-2.41$ au. These measurements of CO, \ce{CO2}, \ce{H2O}, \ce{CH3OH}, and \ce{CH4} sampled molecular chemistry in 3I/ATLAS as it receded from its encounter with our Sun and entered the vicinity of the \ce{H2O} ice line --- the region between \rh{} = $2-3$ au where the temperature becomes too low for H$_2$O to vigorously sublime and CO and \ce{CO2} begin to control the overall activity. CO was the most abundant molecule, followed by \ce{H2O} and \ce{CO2}, whose molecular abundances with respect to CO were $(40.5\pm3.1)\%$ and ($41.6\pm0.3)\%$, respectively. This work presents spatial-spectral maps of column density and rotational temperature as a function of distance from the nucleus for all detected species. The spatial distributions of both quantities were highly anisotropic for the apolar species in the coma of 3I/ATLAS, yet were more nearly symmetric for the polar molecules. These results demonstrate how volatiles were segregated in the nucleus ices of 3I/ATLAS and reveal heating and cooling mechanisms in its coma. Derived maps of the ortho-to-para ratio (OPR) for \ce{H2O} were flat with increasing distance from the nucleus and consistent with a coma-averaged value $\mathrm{OPR}=2.7\pm0.2$, slightly less than the expected equilibrium value of three.

Nathan X Roth↗
DOE OSTI2026

NOx mediated formation of ε – UO 3 during thermal decomposition of uranyl nitrate hexahydrate under static air conditions

The thermal decomposition of uranyl nitrate hexahydrate in air at temperatures of 650 – 800 °C is expected to yield α – U 3 O 8 through intermediate UO 3 phases. Here, n this study, Raman spectroscopy complimented by powder X-ray diffraction revealed the unexpected formation of ε – UO 3 as an accompanying phase during thermal decomposition of UNH at 700 °C under static air conditions. Experiments demonstrated that ε – UO 3 does not form during the initial ramp up stage nor during the high temperature heating period; but instead forms during the cooling stage between 250 – 400 °C. Additional experiments revealed that UNH initially decomposes through an amorphous UO 3 intermediate prior to α – U 3 O 8 formation. Attempts to bypass the U 3 O 8 precursor in the formation of ε – UO 3 were unsuccessful, supporting the necessity of U 3 O 8 (α – U 3 O 8 in this study) in the formation of ε – UO 3 . The observed phase evolution is proposed to result from the NOx species generated during UNH thermal decomposition, which create localized oxidizing conditions within the furnace under static air conditions. Furthermore, a predominantly phase pure ε – UO 3 was synthesized under static air conditions through the addition of an extra plateau at 250 °C during the cooling step. These findings demonstrate the importance of gas-phase chemistry and cooling conditions in uranium oxide phase evolution and provide additional insight into ε – UO 3 formation pathways.

Epsilon uranium trioxide↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2026

Reconstruction and Dissolution of Copper Catalysts during Electrocatalytic Nitrate Reduction

Copper-based electrocatalysts are widely explored for electrochemical nitrate remediation, yet their stability under operating conditions remains poorly understood. While nanoscale and subnanoscale Cu motifs are known to restructure during the nitrate reduction reaction (NO 3 RR), how these transformations translate into irreversible material loss remains unclear. Here, we quantify Cu dissolution during NO 3 RR as a function of catalyst architecture and electrolyte chemistry using two model systems: single-atom Cu (Cu 1 ) and Cu nanoparticles (Cu NP ). Time-resolved leaching measurements, integrated with in situ X-ray absorption spectroscopy (XAS), reveal measurable Cu loss for both catalysts during NO 3 RR. Dissolution is concentrated at the initiation of electrolysis, coinciding with rapid restructuring, and depends strongly on morphology, with Cu NP consistently exhibiting greater Cu loss than Cu 1 . In situ XAS reveals that Cu 1 forms transient metallic clusters under reduction that largely redisperse upon returning to open-circuit voltage, whereas Cu NP undergoes reduction of an oxidized surface layer accompanied by sustained Cu loss during electrolysis. Notably, the presence of nitrate significantly intensifies restructuring and Cu loss, highlighting the critical role of electrolyte composition. Furthermore, these findings establish a direct link between electrochemical restructuring and Cu dissolution, unveiling electrolyte-dependent interactions as key determinants of catalyst durability in electrochemical nitrate conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mesoscale Magnetostructural Phase Separation in Fe‐deficient Fe 5 GeTe 2

Two-dimensional van der Waals ferromagnet Fe 5-x GeTe 2 (F5GT) is promising for spintronic applications due to its high Curie temperature, layered structure, and ability to host complex magnetic textures. However, the origin of its sample-dependent magnetic anisotropy remains unclear, hindering control of its magnetic behavior. Here, we use spatially resolved cryogenic scanning transmission electron microscopy (STEM) to correlatively map magnetism, lattice structure, and chemistry across atomic-to-micron scales. We reveal that only mesoscale, not nanoscale, inclusions of a Fe-deficient secondary phase significantly modify magnetic behavior, establishing a previously unrecognized critical length scale. This phase separation, induced by quenching, leads to in-plane magnetic anisotropy, while slow cooling confines separation to a few nanometers and preserves out-of-plane anisotropy. These findings reconcile prior inconsistencies and establish a predictive framework for tuning magnetism in F5GT through thermal processing, with broader implications for controlling anisotropy in other two-dimensional magnetic materials.

2D ferromagnets↗
DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗
NASA NTRSDate not supplied

Experimental Characterization of Additively Manufactured Nickel-Titanium Shape Memory Alloy Heat Pipes

Shape memory alloys (SMA) have been identified for use in spacecraft components as replacement for conventional deployment mechanisms. They may be used in thermal management components such as radiators to create self-deploying radiators. One SMA, NiTi, has also been developed for additive manufacturing processes. Heat pipes are a common way to create highly effective and lightweight spaceflight radiators, and heat pipes can also be made from NiTi and related alloys. The wick is the critical element of a functioning heat pipe, and recent progress over the past years has led to the development of additively manufactured heat pipe wicks in various materials. The combination of these efforts is the focus of this project: creating an additively manufactured, shape memory alloy self-deploying heat pipe radiator. This paper will focus on the experimental characterization of these additively manufactured NiTi heat pipes. The heat pipe coupons were additively manufactured by direct metal laser sintering (DMLS), with an integral liquid cooled condenser. Heat is input to the heat pipe via a thin film heater. Thermocouples were spot welded to the heat pipes to measure temperature at several axial locations. The heat pipes were tested with two working fluids: water and ethanol. Ethanol is not an ideal working fluid for heat pipes but is useful in characterizing them because it wets well to a wide variety of surfaces. Water is in general a superior working fluid for heat pipes, but its contact angle and therefore wicking performance strongly depends on the surface chemistry of the surface it is in contact with. A particular measurement of interest in this test is the evaporator to condenser thermal conductance, which will be compared in the full paper to recently published correlations for additively manufactured heat pipes. Experimental results for two straight geometry and one bellows geometry heat pipe will be presented. The bellows geometry is of interest for condenser of the self-deploying radiator design.

Additive manufacturing↗
DOE OSTI2026

Comparing the fluorination of Pb+ and Po+: insights into p-block relativistic effects for superheavy elements

It is expected that the chemical properties of the heaviest of the superheavy elements (SHEs, = 113 – 118) do not align with what is suggested by their current positions on the periodic table. Specifically, the onset of significant relativistic effects, including increased spin–orbit splitting of the p-orbitals may lead to enhanced stability of low-oxidation states. Notably, it is predicted that flerovium ( = 114) may exhibit pseudo-noble gas behavior from a electron configuration that acts as an ‘inert pair’. Even though it is expected that the – splitting becomes pronounced toward the end of the p-block’s sixth row, there is currently limited experimental evidence to confirm its extent or impact. Here, the production of gas-phase lead ( = 82) and polonium ( = 84) fluoride cations (PbF and PoF ) were compared to elucidate differences in accessible oxidation states. The PoF and PbF species were produced and identified with the FIONA spectrometer at the Lawrence Berkeley National Laboratory 88-Inch Cylctron Facility. Polonium showed notably different fluorination chemistry than lead, producing PoF and PoF as primary products compared to lead’s PbF and PbF . The distribution of PoF products observed offer insights as to the role of spin–orbit splitting for polonium. Similar studies of superheavy elements would elucidate the accessibility of their low-oxidation states as well as to inform chemical predictions for eighth-row elements not yet discovered.

Pore, Jennifer L [Nuclear Science Division, Lawren↗
NASA NTRS2026

Scientific Highlights and Status of the Mars Science Laboratory Mission, Years 9 through 12

We present the scientific highlights and status from years 9-12 of NASA’s Mars Science Laboratory mission, with its Curiosity rover. A thermochemolysis-based sample analysis revealed a variety of organic molecules preserved in clay mineral-rich lake/lake margin deposits of the Glen Torridon region. The stratigraphically higher clay-sulfate mineral transition on Aeolis Mons was found to be associated with a lithological change from lacustrine to dry aeolian depositional environments. Water intermittently returned to the surface after the increase in aridity, including during the deposition of wave ripple strata which indicate a climate capable of sustaining an ice-free lake in the early Hesperian. The overlying Mg sulfate-bearing unit has both poly- and monohydrated Mg sulfates and is also enriched in iron carbonate of a purity consistent with chemical sedimentation and therefore atmospheric sequestration. A channel-incised canyon within Gediz Vallis and distal alluvial deposits recorded numerous fluvial and debris flow events after the deposition, lithification, and erosion of the clay- and sulfate-bearing units. Bright clasts found in Gediz Vallis canyon are composed of native sulfur. The investigation of the modern atmosphere and environment focused on the influence of the crater and mountain topography and captured high-altitude noctilucent and iridescent clouds. Radiation monitoring now spans an 11-year solar cycle. Rover systems and science instruments remain capable of addressing mission scientific objectives. All ten scientific instruments continue to return high-value data despite degradation of the Chemistry and Camera instrument’s high-voltage system, loss of multispectral imaging capability, and loss of active neutron spectroscopy. The mission has mitigated the loss of redundant brakes on two arm joints, additional wheel wear, and highly reduced capability of the redundant flight computer. Efficiency improvements in rover energy utilization and operations staffing have allowed the mission team to maintain a high level of productivity despite declining rover power generation and funding.

Mars↗
DOE OSTI2022

In Situ Tracking of Nonthermal Plasma Etching of ZIF-8 Films

Surface characterization is critical for understanding the processes used for preparing catalysts, sorbents, and membranes. Nonthermal plasma (NTP) is a process that achieves high reactivity at low temperatures and is used to tailor the surface properties of materials. In this work, we combine the capabilities of infrared reflection absorption spectroscopy (IRRAS) with NTP for the in situ interrogation of zeolitic imidazolate framework-8 (ZIF-8) thin films to probe modifications in the material induced by oxygen and nitrogen plasmas. The IRRAS measurements in oxygen plasma reveal etching of organic ligands with sequential removal of the methyl group and imidazole ring and with the formation of carbonyl moieties (C═O). In contrast, nitrogen plasma induces mild etching and grafting of nitrile groups (−C≡N). Scanning electron microscopy imaging shows that oxygen plasma, at prolonged times, significantly degrades the ZIF-8 film at the grain boundaries. Treatment of ZIF-8 membranes using mild plasma conditions yields a fivefold enhancement for H 2 /N 2 and CO 2 /CH 4 ideal selectivities and an eightfold enhancement for CO 2 /N 2 ideal selectivity. Additionally, the new tools described here can be used for spectroscopic in situ tracking of plasma-induced chemistry on thin films in general.

IRRAS↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
NASA NTRS2017

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics↗