Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimal transport”

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.

580 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Active Learning for Rapid Targeted Synthesis of Compositionally Complex Alloys

The next generation of advanced materials is tending toward increasingly complex compositions. Synthesizing precise composition is time-consuming and becomes exponentially demanding with increasing compositional complexity. An experienced human operator does significantly better than a novice but still struggles to consistently achieve precision when synthesis parameters are coupled. The time to optimize synthesis becomes a barrier to exploring scientifically and technologically exciting compositionally complex materials. This investigation demonstrates an active learning (AL) approach for optimizing physical vapor deposition synthesis of thin-film alloys with up to five principal elements. We compared AL-based on Gaussian process (GP) and random forest (RF) models. The best performing models were able to discover synthesis parameters for a target quinary alloy in 14 iterations. We also demonstrate the capability of these models to be used in transfer learning tasks. RF and GP models trained on lower dimensional systems (i.e., ternary, quarternary) show an immediate improvement in prediction accuracy compared to models trained only on quinary samples. Furthermore, samples that only share a few elements in common with the target composition can be used for model pre-training. We believe that such AL approaches can be widely adapted to significantly accelerate the exploration of compositionally complex materials.

Chemistry↗
NASA NTRS2026

The Contaminant Footprint of Landed Spacecraft: Toward an Inventory and Modelling Framework

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in near-vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination of solar system bodies may irrevocably degrade targets of unique value to planetary scientists, for instance, as lunar landed spacecraft introduce propellant effluents into the otherwise pristine ice of the Moon’s permanently shadowed regions. NASA’s planetary protection discipline seeks to ensure that solar system bodies are not contaminated, for scientific purposes, by terrestrial material (i.e., forward contamination). This interest aligns with planetary science interest in mitigating the transport of terrestrial contaminants onto solar system bodies and in controlling types of contamination that could compromise the scientific value of samples or measurements. NASA, and other entities that practice planetary science, have compelling and multidisciplinary interests in the preservation of special regions and sampling sites of high scientific value – including lunar permanently shadowed regions (PSRs) – from inadvertent contamination by any spacecraft, and in understanding the contamination of such regions by all spacecraft. Organic molecular contamination here represents a primary threat. Organic molecules will be introduced to solar system bodies by nominal landed spacecraft and crew processes – including by the action of descent and ascent engines; natural materials outgassing; and crew environmental and life support system sources. Molecular contaminants can also travel in the free-molecular sense at global scale across near-vacuum bodies, including into regions where they may be permanently trapped. This presentation will address a high-level study to identify sources of contaminants – in particular, organic material – generated by landed spacecraft along with the transport vectors by which these contaminants can reach sites of scientific interest on bodies like the Moon. A vision for an integrated modeling framework for the organic contamination footprint of spacecraft missions, individually and collectively, will also be described and presented along with initial conclusions related to organic molecular transport.

Gas Dynamics↗
NASA NTRSDate not supplied

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins↗
DOE OSTI2026

Weak, shallow, dry convection over Angola increases offshore stratocumulus cloud droplet number concentrations

Boundary-layer cloud interactions involving shortwave-absorbing aerosols remain one of the least understood aerosol influences on climate. Here, we find the highest stratocumulus cloud droplet number concentrations over the southeast Atlantic occur when agricultural fires coincide with synoptically-weakened surface warming over Angola, occurring June-early August. Dry convection fills a shallow continental boundary layer with smoke, and a nighttime (local solar time 2-9) land breeze transports the aerosol into the marine boundary layer. Offshore aerosol transport is strengthened by low-level easterlies from a continental pressure high southeast of Angola. Simultaneously, the South Atlantic subtropical high is weaker, allowing extensive dispersal of aerosol offshore into the boundary layer, obscuring cloud brightening from shipping. Meteorological co-variation at synoptic scales compensates for cloud brightening by the smoke. Outgoing shortwave radiation increases by 15–20% of the monthly mean in June and July when offshore droplet numbers are less but the stratocumulus deck is more developed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revealing key structures for reversible sulfur redox in amorphous polymeric sulfur

Amorphous polymeric sulfur cathodes, such as sulfurized polyacrylonitrile (SPAN), enable high-energy lithium-sulfur batteries without cobalt or nickel, leveraging abundant sulfur. However, the limited in situ understanding of their synthesis and electrochemistry has impeded targeted optimization. Here, in this study, we integrate operando high-energy total scattering with sulfur K-edge X-ray absorption spectroscopy to monitor SPAN's formation and cycling in real time. Our results show that S-C bond formation halts further fusion of cyclized polyacrylonitrile, fostering π-π stacking and a transition from long-chain to short-chain sulfur-critical for reversible sulfur redox. These features synergistically minimize polysulfide dissolution and charge-transfer resistance, enabling optimized SPAN to achieve high capacity retention over 1,000 cycles. Operando X-ray absorption spectroscopy reveals that residual protons drive thiol-thione tautomerism, with lithium replacement during the first discharge causing ~20% irreversible capacity loss. To enhance performance, minimizing -NH groups and expanding pyridine networks are key. These findings transform SPAN optimization from empirical tuning to mechanism‑guided engineering and point the way towards sulfur loadings and energy densities competitive with state‑of‑the‑art Li‑ion cathodes.

25 ENERGY STORAGE↗
NASA NTRS2012

Returning an Entire Near-Earth Asteroid in Support of Human Exploration Beyond Low-Earth Orbit

This paper describes the results of a study into the feasibility of identifying, robotically capturing, and returning an entire Near-Earth Asteroid (NEA) to the vicinity of the Earth by the middle of the next decade. The feasibility of such an asteroid retrieval mission hinges on finding an overlap between the smallest NEAs that could be reasonably discovered and characterized and the largest NEAs that could be captured and transported in a reasonable flight time. This overlap appears to be centered on NEAs roughly 7 m in diameter corresponding to masses in the range of 250,000 kg to 1,000,000 kg. The study concluded that it would be possible to return a approx.500,000-kg NEA to high lunar orbit by around 2025. The feasibility is enabled by three key developments: the ability to discover and characterize an adequate number of sufficiently small near-Earth asteroids for capture and return; the ability to implement sufficiently powerful solar electric propulsion systems to enable transportation of the captured NEA; and the proposed human presence in cislunar space in the 2020s enabling exploration and exploitation of the returned NEA. Placing a 500-t asteroid in high lunar orbit would provide a unique, meaningful, and affordable destination for astronaut crews in the next decade. This disruptive capability would have a positive impact on a wide range of the nation's human space exploration interests. It would provide a high-value target in cislunar space that would require a human presence to take full advantage of this new resource. It would offer an affordable path to providing operational experience with astronauts working around and with a NEA that could feed forward to much longer duration human missions to larger NEAs in deep space. It represents a new synergy between robotic and human missions in which robotic spacecraft would retrieve significant quantities of valuable resources for exploitation by astronaut crews to enable human exploration farther out into the solar system. The capture, transportation, examination, and dissection of an entire NEA would provide valuable information for planetary defense activities that may someday have to deflect a much larger near-Earth object. Transportation of the NEA to lunar orbit with a total flight time of 6 to 10 years would be enabled by a ~40-kW solar electric propulsion system with a specific impulse of 3,000 s. The flight system could be launched to low-Earth orbit (LEO) on a single Atlas V-class launch vehicle, and return to lunar orbit a NEA with at least 28 times the mass launched to LEO. Longer flight times, higher power SEP systems, or a target asteroid in a particularly favorable orbit could increase the mass amplification factor from 28-to-1 to 70-to-1 or greater. The NASA GRC COMPASS team estimated the full life-cycle cost of an asteroid capture and return mission at approx.$2.6B.

asteroid mining↗
DOE OSTI2026

Advanced LiFSI-LiPF6 Electrolyte for Wide-Temperature and Thermally Stable Lithium-Ion Batteries

A new optimized LiFSI–LiPF6 dual-salt controlled-solvation electrolyte (E-DS) is demonstrated to enable practical graphite||LiNi0.8Mn0.1Co0.1O2 cells (˜4.0 mAh cm?²) to achieve exceptional performance and safety under extreme conditions. By optimizing anion coordination with the smaller, more dissociating FSI? anion, the E-DS forms ultrathin, dense, and inorganic-rich electrode/electrolyte interphases that dramatically suppress solvent decomposition, transition-metal dissolution, and surface reconstruction compared to the conventional LiPF6/carbonate electrolyte. Consequently, E-DS cells deliver >78% capacity retention after 300 cycles at 60 °C, retain fast discharging capacity at 30 °C, and operate effectively at -20 °C. Most strikingly, fully charged full cells with E-DS, even under overcharging to 4.8 V, show a lower heat evolution in stable formulations — transforming a traditionally unstable high-voltage/high-temperature configuration into an intrinsically safe state. This work establishes a new benchmark for carbonate-containing electrolytes, simultaneously achieving high energy density, fast-discharging capability, wide-temperature operation (-20 to 60 °C), and outstanding thermal safety in nickel-rich lithium-ion batteries.

electrode/electrolyte interphase↗
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↗
DOE OSTI2025

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗
DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗
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↗
DOE OSTI2026

Characterization of NBCMA—A DAPO Precursor

As part of our investigation into the process optimization of 2,6-diaminopyrazine-1-oxide (DAPO), small-scale sensitivity testing (impact, friction, and ESD) was conducted on the precursor, N-nitroso-bis(cyanomethyl)amine (NBCMA). NBCMA exhibited moderate sensitivity to friction, whereas no sensitivity to impact or spark was observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Spacecraft Conceptual Design for Returning Entire Near-Earth Asteroids

In situ resource utilization (ISRU) in general, and asteroid mining in particular are ideas that have been around for a long time, and for good reason. It is clear that ultimately human exploration beyond low-Earth orbit will have to utilize the material resources available in space. Historically, the lack of sufficiently capable in-space transportation has been one of the key impediments to the harvesting of near-Earth asteroid resources. With the advent of high-power (or order 40 kW) solar electric propulsion systems, that impediment is being removed. High-power solar electric propulsion (SEP) would be enabling for the exploitation of asteroid resources. The design of a 40-kW end-of-life SEP system is presented that could rendezvous with, capture, and subsequently transport a 1,000-metric-ton near-Earth asteroid back to cislunar space. The conceptual spacecraft design was developed by the Collaborative Modeling for Parametric Assessment of Space Systems (COMPASS) team at the Glenn Research Center in collaboration with the Keck Institute for Space Studies (KISS) team assembled to investigate the feasibility of an asteroid retrieval mission. Returning such an object to cislunar space would enable astronaut crews to inspect, sample, dissect, and ultimately determine how to extract the desired materials from the asteroid. This process could jump-start the entire ISRU industry.

solar electric propulsion↗
DOE OSTI2026

Scaling Real Estate Development Innovation With Emerging Local Area Developers

These case studies provide tangible examples within the development field that show how true innovation in real estate and property development hinges on shaping and optimizing diverse connections and relationships to create widely available, market-viable, and high-performance buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Probing the role of local tunnel variations in early-stage lithiation of α-MnO₂ nanowires via in situ TEM

Understanding lithium-ion transport in tunnel-structured manganese oxides is essential for designing high-performance lithium-ion battery electrode materials. Here, we elucidate the early-stage lithiation mechanism of potassium-stabilized α-MnO 2 nanowires using in situ transmission electron microscopy (TEM) coupled with electron energy-loss spectroscopy (EELS), high-resolution TEM (HRTEM), and geometric phase analysis (GPA). Real-time TEM imaging reveals clear volume expansion at the reaction front, while EELS analysis uncovers lithium-ion diffusion far beyond this region, where no visible expansion is observed, indicating fast, defect-assisted transport. GPA and HRTEM analyses show that localized tensile and compressive strain fields, originating from pre-existing local tunnel structural variations, persist after lithiation. The tensile-strained regions enable lithium-ion insertion with minimal lattice distortion, offering additional free volume that facilitates rapid lithium-ion accommodation ahead of the structural transformation. Our results demonstrate a local tunnel variation-mediated fast diffusion pathway that precedes bulk reaction, underscoring the critical role of local strain in enabling early-stage lithium transport. Given the structural versatility of MnO 2 and its ability to accommodate diverse atomic arrangements beyond the well-known tunnel phases (β-, γ-, δ-, λ-, R-phases), our findings highlight the importance of understanding and engineering local structural environments. This work provides fundamental insights into the interplay between defects, strain, and ion dynamics, and presents defect engineering as a promising approach to enhance both rate performance and structural stability in manganese-based cathodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗