Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Physics and Simulations”

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.

3,226 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

User Interface Issues in Supporting Human-Computer Integrated Scheduling

A major problem in designing user interfaces for scheduling systems is one of allowing the human to become an integral part of the system. The human role in scheduling extends beyond the simple tasks of providing the input and accepting the output. Because of the inherent intractability of most real-world scheduling problems, intelligence must be incorporated into the scheduling process in order to reach an acceptable solution in a reasonable amount of time. Artificial Intelligence research has concentrated on identifying algorithms and heuristics for this purpose. However, interfaces which allow the scheduler to take advantage of human intelligence and allow the user insight into and influence over the planning process are also needed.

Lynne P Cooper↗
DOE OSTI2026

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences↗
DOE OSTI2025

Robust Spectral Anomaly Detection in EELS Spectral Images via 3D Convolutional Variational Autoencoders

Abstract A 3D Convolutional Variational Autoencoder (3D‐CVAE) is introduced for automated anomaly detection in electron energy‐loss spectroscopy spectrum imaging (EELS‐SI) data. This approach leverages the full 3D structure of EELS‐SI data to detect subtle spectral anomalies while preserving both spatial and spectral correlations across the datacube. By employing cross‐entropy loss and training on bulk spectra, the model learns to reconstruct bulk features characteristic of the defect‐free material. In exploring methods for anomaly detection, both the 3D‐CVAE approach and principal component analysis (PCA) are evaluated, testing their performance using FeL‐edge ΔEpeak shifts designed to simulate material defects. These results show that 3D‐CVAE achieves superior anomaly detection and maintains consistent performance across various shift magnitudes. The method demonstrates clear bimodal separation between bulk and anomalous spectra, enabling reliable classification. Further analysis verifies that lower‐dimensional representations are robust to anomalies in the data. While performance advantages over PCA diminish with decreasing anomaly concentration, our method maintains high reconstruction quality even in challenging, noise‐dominated spectral regions. This approach provides a robust framework for unsupervised automated detection of spectral anomalies in EELS‐SI data, particularly valuable for analyzing complex material systems.

Chemistry↗
DOE OSTI2025

Revealing the evolution of order in materials microstructures using multi-modal computer vision

The development of high-performance materials for microelectronics, energy storage, and extreme environments depends on our ability to describe and direct property-defining microstructural order. Our present understanding is typically derived from laborious manual analysis of imaging and spectroscopy data, which is difficult to scale, challenging to reproduce, and lacks the ability to reveal latent associations needed for mechanistic models. Here, we demonstrate a multi-modal machine learning (ML) approach to describe order from electron microscopy analysis of the complex oxide La 1−x Sr x FeO 3 . We construct a hybrid pipeline based on fully and semi-supervised classification, allowing us to evaluate both the characteristics of each data modality and the value each modality adds to the ensemble. We observe distinct differences in the performance of uni- and multi-modal models, from which we draw general lessons in describing crystal order using computer vision.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using scalable computer vision to automate high-throughput semiconductor characterization

Abstract High-throughput materials synthesis methods, crucial for discovering novel functional materials, face a bottleneck in property characterization. These high-throughput synthesis tools produce 10 4 samples per hour using ink-based deposition while most characterization methods are either slow (conventional rates of 10 1 samples per hour) or rigid (e.g., designed for standard thin films), resulting in a bottleneck. To address this, we propose automated characterization (autocharacterization) tools that leverage adaptive computer vision for an 85x faster throughput compared to non-automated workflows. Our tools include a generalizable composition mapping tool and two scalable autocharacterization algorithms that: (1) autonomously compute the band gaps of 200 compositions in 6 minutes, and (2) autonomously compute the environmental stability of 200 compositions in 20 minutes, achieving 98.5% and 96.9% accuracy, respectively, when benchmarked against domain expert manual evaluation. These tools, demonstrated on the formamidinium (FA) and methylammonium (MA) mixed-cation perovskite system FA 1−x MA x PbI 3 , 0 ≤ x ≤ 1, significantly accelerate the characterization process, synchronizing it closer to the rate of high-throughput synthesis.

Science & Technology - Other Topics↗
NASA NTRSDate not supplied

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗
DOE OSTI2023

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

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↗