Search NASASearch

SEARCH · Search NASA

Results for “Numerical Modeling”

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.

974 records · Page 43

Space Shuttle Main Engine Propellant Path Leak Detection Using Sequential Image Processing

Initial research in this study using theoretical radiation transport models established that the occurrence of a leak is accompanies by a sudden but sustained change in intensity in a given region of an image. In this phase, temporal processing of video images on a frame-by-frame basis was used to detect leaks within a given field of view. The leak detection algorithm developed in this study consists of a digital highpass filter cascaded with a moving average filter. The absolute value of the resulting discrete sequence is then taken and compared to a threshold value to produce the binary leak/no leak decision at each point in the image. Alternatively, averaging over the full frame of the output image produces a single time-varying mean value estimate that is indicative of the intensity and extent of a leak. Laboratory experiments were conducted in which artificially created leaks on a simulated SSME background were produced and recorded from a visible wavelength video camera. This data was processed frame-by-frame over the time interval of interest using an image processor implementation of the leak detection algorithm. In addition, a 20 second video sequence of an actual SSME failure was analyzed using this technique. The resulting output image sequences and plots of the full frame mean value versus time verify the effectiveness of the system.

L Montgomery Smith

Isotopic Evidence for A Cold and Distant Origin of 3I/Atlas

Interstellar objects provide the only directly observable samples of icy planetesimals formed around other stars, and can therefore provide insight into the diversity of physical and chemical conditions occurring during exoplanet formation1−3. Here we report isotopic measurements of the interstellar comet 3I/ATLAS, which reveal an elemental composition unlike any Solar System body. The water in 3I/ATLAS is enriched in deuterium, at a level of D/H = (0.98 ± 0.06)%, which is more than an order of magnitude higher than in known comets, while its range of 12 C/ 13 C ratios (141–191 for CO 2 and 123–172 for CO) exceeds typical values found in the Solar System, as well as nearby interstellar clouds and protoplanetary disks. Such extreme isotopic signatures indicate formation at temperatures ≲ 30 K in a relatively metal-poor environment. When interpreted with respect to models for Galactic chemical evolution, the carbon isotopic composition implies that 3I/ATLAS may have accreted as long ago as 12 billion years, following a period of intense, early star formation. 3I/ATLAS thus represents a preserved fragment of an ancient planetary system.

Martin A Cordiner

Optimizing Deep Geothermal Drilling for Energy Sustainability in the Appalachian Basin

This study investigates the geological and geomechanical characteristics of the MIP 1S geothermal well in the Appalachian Basin to optimize drilling and address the wellbore stability issues encountered. Data from well logs, sidewall core analysis, and injection tests were used to derive elastic and rock strength properties, as well as stress and pore pressure profiles. A robust 1D-geomechanical model was developed and validated, correlating strongly with wellbore instability observations. This revealed significant wellbore breakout, widening the diameter from 12 ¼ inches to over 16 inches. Advanced technologies like Cerebro Force™ In-Bit Sensing were used to monitor drilling performance with high accuracy. This technology tracks critical metrics such as bit acceleration, vibration in the x, y, and z directions, Gyro RPM, stick-slip indicators, and bending on the bit. Cerebro Force™ readings identified hole drag caused by poor hole conditions, including friction between the drill string and wellbore walls and the presence of cuttings or debris. This led to higher torque and weight on bit (WOB) readings at the surface compared to downhole measurements, affecting drilling efficiency and wellbore stability. Optimal drilling parameters for future deep geothermal wells were determined based on these findings.

Environmental Sciences & Ecology

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science

Mechanics and Cooling of the CMS Phase-II Tracker Barrel

The CMS upgrade for the High-Luminosity LHC (HL-LHC) era will allow for an incredibly high rate of data collection, to be used in the next generation of precision measurements and searches for beyond the standard model physics. A key element of this upgrade is the full replacement of the silicon strip tracker to enable level-1 triggering on track level primitives, produced by pixel-strip (PS) and strip-strip (2S) modules, tiled together to form the new detector. Mechanically supporting and cooling these modules is not trivial, as material budget in the tracker volume is highly constrained, and as the modules accrue radiation damage they risk becoming entirely non-functional if the silicon sensors are not sufficiently cooled. To address this, each region of the tracker has a bespoke solution for housing and cooling their modules, mostly depending on carbon-fiber and carbon-foam structures, with dual phase CO2 cooling. For the most central part of the tracker, the Flat Barrel with PS modules (Flat TBPS) covers the pseudorapidity range out to $|\eta|=0.4$, and the mechanical structures are manufactured, assembled, and tested entirely by Fermilab and UC Davis. This presentation will describe the intense design and testing requirements for these components, and how they will ensure strong performance of the tracker through the HL-LHC era.

Cosby, Christopher [Fermilab]

Enter the AHU (36th Chamber of ASHRAE): A Multi-site Field Study of ASHRAE G36

Despite being recognized as the best practice for advanced building controls, ASHRAE Guideline 36 (G36) has seen slow adoption in retrofit cases. Decisionmakers lack credible field evidence to justify the time and person-power investment. Most prior analyses have relied on software simulations, which overlook implementation challenges and fail to persuade owners to move from models to real-world deployment. This paper presents a multi-site field study of G36 performance, drawing on measured results from 17 projects across diverse building types and climate zones. The analysis disaggregates outcomes by the most widely adopted air handling unit (AHU) based G36 strategies, including trim-and-respond approaches to supply air temperature (SAT) and duct static pressure (DSP) reset, and economizer controls. Results for controls re programming implementations are encouraging, with HVAC savings ranging from 2% - 49%, with a median of 18%. The range aligns with simulation study findings showing 1% - 46% savings through these three strategies. These findings show that even without capital investment in control infrastructure upgrades, existing building owners are reaping significant benefits from updating HVAC sequences of operation to industry best-practice solutions. The results give practitioners and decisionmakers a reference point on what to expect, which strategies deliver, and how field performance may compare to simulation.

Deshpande, Reva

AL-CRADA-2025-01 Final Report: American-Made Wind Turbine Materials Recycling Prize Accelerate! Phase 2

This report summarizes the Phase 2 achievements of the Wind Turbine Materials Recycling Prize project led by Critical Materials Recycling (CMR) in partnership with Ames National Laboratory. The team demonstrated the technical feasibility and economic viability of recovering and reusing rare-earth Nd-Fe-B magnets from end-of-life (EOL) wind turbine generators. Through mechanical disassembly, demagnetization, and precision processing, magnets were harvested and characterized to assess structural, compositional, and magnetic properties. Results showed that recovered magnets performed on par with commercial-grade counterparts, enabling their integration into actuator and motor prototypes without redesign. System-level modeling and finite-element simulations validated operational performance, while Life Cycle Cost (LCC) analysis revealed up to 20% material cost savings and 2–3% reductions in overall motor system costs. Life Cycle Assessment (LCA) further confirmed substantial environmental benefits, including a ~94% reduction in global warming potential relative to virgin magnet production. The project’s success was enabled by a tightly integrated partnership between industry and national lab experts, illustrating a viable path for domestic magnet-to-magnet recycling. These findings support a scalable, circular solution for rare-earth material recovery, reducing dependence on imported critical minerals and advancing sustainable clean energy technologies.

Lograsso, Thomas [Ames Laboratory (AMES), Ames, IA

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)

Changes in Arctic and Antarctic Sea-Ice Properties and Processes

Although sea ice covers ~9% of the world’s oceans, its presence, and absence, vastly alter the energetic balance of Earth’s climate. Because of geographic differences, Arctic and Antarctic sea ice have different dominant properties and processes, which produce different responses to natural variability and anthropogenic change. In this Review, we survey changes in sea-ice properties and processes of the Arctic and Antarctic environments. In the Arctic, the melt season has lengthened by ~7.4 days decade −1 (1979-2024). There has been a corresponding decrease in winter sea-ice thickness (1.6 m mean total over 1980-2023), summer surface albedo (0.03 decade −1 over 1979-2020), and spring snow depth (2.5 cm decade −1 over 1954-2024). In comparison, corresponding changes in the Antarctic are small overall due to contrasting regional and temporal trends. From a dynamic perspective, sea-ice motion in both hemispheres has increased, with 0.63 cm s −1 decade −1 in the Arctic (1978-2024) and 0.69 cm s −1 decade −1 in the Antarctic (1982-2024), whereas trends in polynya occurrence differ by region. Future research should prioritize sea-ice observations that can inform how the climate system responds to natural variability and anthropogenic change, be standardized to increase uptake by the broader scientific community, and target process-oriented interactions to support model development.

Melinda A Webster

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation

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

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.

In-Situ Visualization of a Growing Brittle Crack in Aluminum Oxynitride Using Synchrotron X-Rays and the Double-Cleavage Drilled Compression Geometry

Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.

Gorske, Sara F.

Dissecting Disorder: Defect-Driven Structural Complexity in Layered Li3InCl6 Solid Electrolyte

Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.

Liu, Jue [ORNL] (ORCID:000000024453910X)

Under Pressure: The Artemis I Heatshield Char Loss Investigation and Artemis II’s Successful Entry

Following the successful Artemis I skip reentry on December 11, 2022, unexpected char liberation from the Orion heat- shield prompted formation of an Anomaly Response Team to determine root cause and establish flight rationale and corrective actions for subsequent missions. Through a comprehensive investigation including detailed hardware analysis, modeling & simulation, and ground testing, the team determined that the low permeability Avcoat experienced extreme internal gas pressure buildup during the skip entry that could not adequately outgas, leading to crack formation and char liberation. Key discoveries included significant performance differences between permeable and impermeable heatshield regions and successful replication of the anomaly through ground testing. Based on extensive ground testing and analysis, Artemis II flew a modified trajectory without the skip entry to minimize crack-inducing conditions, and successfully splashed down on April 10, 2026 with significantly reduced char loss. For Artemis III and beyond, a corrective action was implemented to produce a more permeable version of Avcoat that will substantially reduce internal pressure accumulation during entry.

Artemis II

Robust Biaxial Anisotropy and Switchable Néel Vectors in LaFeO 3 Epitaxial Films

Antiferromagnets with highly stable but switchable Néel vectors are desired for antiferromagnetic spintronics with ultrafast speed and terahertz frequencies. Electrical switching of antiferromagnetic insulators has been demonstrated using binary antiferromagnets, while large families of complex antiferromagnets such as perovskites are largely unexplored. Here, we show that epitaxial LaFeO 3 thin films on SrTiO 3 (001) exhibit clear, robust biaxial anisotropy with a spin-flop field of a few tesla. Angular-dependent spin-Hall magnetoresistance (SMR) characterizations of Pt/LaFeO 3 bilayers with the current channel along SrTiO 3 [100] and [110] reveal distinct, intriguing shapes and field dependence. Simulations using a macrospin model accurately describe the main behavior and fine features of the SMR data from which key antiferromagnetic parameters are extracted. Furthermore, remanent SMR measurement confirms the high fidelity of the Néel vector along either easy axis of the biaxial anisotropy, indicating that epitaxial films of LaFeO 3 and potentially other perovskite antiferromagnets offer an attractive platform for antiferromagnetic spintronics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Design, Development, and Test of the Advanced Apollo Orbital Assembly System

As spaceflight moves toward commercial solutions for Crew Vehicles and Space Stations, opportunity exists to lower costs with novel designs. Probe and cone docking systems provide a lightweight, low cost, and high-performance docking solution. This work revisits the Apollo probe and cone design and modifies it for the requirements of today’s computer-controlled spacecraft. This new system is called the Advanced Apollo Orbital Assembly (APOA) system, and is intended to support satellite, cargo, and space station assembly missions. A simulation of the APOA was programmed, components were sized, and a design was created by scaling the probe to the size of an EELV Secondary Payload Adapter (ESPA) tunnel. Prototype test hardware was 3D printed using Fused Deposition Modeling (FDM) methods with Polylactic Acid (PLA) material. Testing of the APOA-ESPA was conducted at Marshall Space Flight Center’s (MSFC) Flat Floor, and a test-correlated simulation is used to evaluate a Monte-Carlo of Initial Contact Conditions (ICC’s) to establish baseline performance. The successful development, test, and correlation of the APOA-ESPA proves the design validity and increases Technology Readiness Level (TRL) from 2 to 4. This work opens the door to construction of an APOA-ESPA from flight like materials, and to develop a larger scale prototype APOA. When the full scale APOA is incorporated with the Common Berthing Mechanism (CBM), becomes the Hybrid Berthing System (HBS), which allows for berthing without a robotic arm.

Berthing

Lyotropic Liquid Crystalline Phase Behavior of Boron Nitride Nanotube Aqueous Dispersions

Boron nitride nanotubes (BNNTs) are gaining significant interest due to their outstanding mechanical and thermal properties, as well as their potential to serve as a model nanorod system. Processing BNNT liquid crystalline (LC) dispersions enables precise control over BNNT orientation in macroscopic assemblies, while their low absorption in the visible spectrum facilitates studying LCs at exceptionally high concentrations. Here, we investigate the behavior of BNNTs in aqueous solutions stabilized by the surfactant sodium deoxycholate (SDC), examining the effect of BNNT purity and BNNT−SDC concentrations on lyotropic LC formation. We disperse up to 15wt% BNNT in SDC solutions and use polarized light microscopy to detail the transition from an isotropic state to a biphasic regime, where isotropic and nematic domains coexist due to phase separation, to a single fully nematic phase. Cryogenic electron microscopy provides direct evidence of BNNT alignment within nematic domains. Our results show that enhanced depletion-induced attractions, driven by increased surfactant concentration, lower the threshold concentration of BNNT required to form nematic domains. In contrast, low surfactant concentrations relative to BNNT result in insufficiently coated nanotube surfaces, leading to poor dispersions and BNNT aggregation. Additionally, we fabricate well-aligned BNNT films from aqueous LC dispersions. Our findings advance the understanding of BNNT LCs, offering insight into controlling their orientation and highlighting their potential for high-performance materials.

Joe Khoury