Search NASASearch

SEARCH · Search NASA

Results for “COMPUTER METHOD”

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.

777 records · Page 2

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

From zonal to nodal capacity expansion planning: Spatial aggregation impacts on a realistic test-case

Solving power system capacity expansion planning (CEP) problems at realistic spatial resolutions is computationally challenging. Thus, a common practice is to solve CEP over zonal models with low spatial resolution rather than over full-scale nodal power networks. Due to improvements in solving large-scale stochastic mixed integer programs, these computational limitations are becoming less relevant, and the assumption that zonal models are realistic and useful approximations of nodal CEP is worth revisiting. Here, this work is the first to conduct a systematic computational study on the assumption that spatial aggregation can reasonably be used for ISO-scale CEP. By considering a realistic, large-scale test network based on the state of California with over 8000 buses, we find that well-designed small spatial aggregations can yield good approximations but that coarser zonal models may result in large distortions of investment decisions, e.g., capacity under-investment of up to 41% for the lowest resolution model considered.

24 POWER TRANSMISSION AND DISTRIBUTION

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Comparison between the PISO algorithm and preconditioning methods for compressible flow

Two widely used family of algorithms, pressure-based and density-based methods, have been developed for computational fluid dynamics (CFD) problems over the years. Pressure-based methods (such as SIMPLE and PISO) use a Poisson-like equation for updating pressure instead of the continuity equation, while density-based methods use the continuity equation to update density (an equation of state is used to provide density in pressure based schemes and pressure in density based schemes). Pressure-based methods were developed originally for incompressible flows at low Reynolds numbers and were then extended to high Reynolds numbers and compressible applications. On the other hand, density based methods were originally developed for transonic flows and have been extended down to low Mach numbers through the use of preconditioning techniques. We compare these two very different approaches to solving the Navier-Stokes equations in order to gain an understanding of their similarities and differences. Specifically, we consider the PISO scheme as a representative pressure-based method and contrast it with a recently developed preconditioning scheme. We also compare the relative performance of the PISO algorithm with a Euler implicit algorithm that is employed to solve the preconditioned equations by means of a vector stability analysis.

Charles L Merkle

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING

Microscopy modality transfer of steel microstructures: Inferring scanning electron micrographs from optical microscopy using generative AI

Scanning electron microscopy (SEM) is resource intensive, which limits its throughput in some applications. As an alternative, we propose applying computer vision and machine learning to generate high-quality synthetic SEM micrographs from micrographs obtained using light optical microscopy (LOM). Working with a correlated LOM/SEM dataset of dual-phase steel images, we test generative models of various architectures, including encoder-decoder networks, generative adversarial networks (GANs), and diffusion-based models. We find that the diffusion models significantly outperform other methods on both qualitative and quantitative assessments, while preserving key metallurgical meaning. This work establishes diffusion as the state-of-the-art for microscopy modality transfer and demonstrates the potential of AI-powered microscopy to enhance LOM with micron scale structural recreation.

Computer vision

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information

CFD-Assisted Nodal Modeling of Sloshing in a Cryogenic Propellant Tank

During autogenous pressurization, tank sloshing causes a significant increase in pressurant consumption to maintain constant ullage pressure during draining of the tank. This increase in pressurant consumption is caused by a significant increase in condensation at the liquid-vapor interface. Sloshing strongly affects the liquid side heat transfer coefficient and thereby the condensation rate. Traditionally, sloshing is modeled by CFD code using the VOF (Volume of Fluid) method to track the liquid vapor interface during sloshing. CFD calculations require a very fine grid to accurately compute the heat and mass transfer at the interface. Therefore, computations are time consuming and prohibit performing many parametric studies often needed during the design of a new system. This paper describes an alternative approach by developing a CFD-assisted nodal model to predict system parameters more economically with reasonable accuracy. In this approach, a multi-node model of tank pressurization was developed using GFSSP. A multi-node model was needed to account for stratification. The model computes heat and mass transfer at the interface to calculate the condensation rate. The liquid side heat transfer is computed using the parameters of sloshing dynamics such as frequency, wave amplitude, and interface area. The parameters of sloshing dynamics are computed by the CFD code LOCI-Stream. The model predictions were compared with test data for several cases.

Cryogenic Tank Sloshing

Reduced salen in a monomeric Fe complex

We report the first monomeric formally M I transition metal salen complex without ion pairing. A combination of spectroscopy, crystallography, and computations indicates that reduction of an Fe II salen complex gives substantial reduced ligand character. We also present a computational method for accurately predicting 57 Fe Mössbauer parameters of Fe salen complexes.

Feldman, Dana M. [Yale University, New Haven, CT (

Energetics of Single Substitutional Impurities in NiTi

Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.

Binary Alloys

Structured illumination for surface-resolved grazing-incidence X-ray scattering

Grazing-incidence (GI) scattering techniques are widely used to characterize thin films, offering high surface sensitivity and insight into morphology and structure. However, these approaches typically provide statistical averaged information due to elongated footprint or limited spatial resolution due to beam size. Here we introduce a method that combines structured illumination with GI X-ray scattering and leverages our computational imaging approach to resolve local structural details. We demonstrate that our method captures local features of an organic semiconductor thin film without the need for sample rotation as in tomography. The method expands GI techniques from statistical averaging to high-resolution imaging, thereby providing the capability for detailed analysis of local material properties, such as domain shape, orientation and polymorphism, which are critical for advancing material design towards more efficient and tailored materials.

Gursoy, Doga

RADAI: A Large-Scale Realistic Dataset for Radiation Detection Algorithm Development

Open, realistic datasets are essential for developing and benchmarking radiation detection algorithms, yet they remain scarce. The Radiological Anomaly Detection and Identification (RADAI) project was develop to create datasets that meet the training and testing needs for sophisticated radiation detection algorithms. The RADAI dataset is a large-scale synthetic resource that integrates high-fidelity Monte Carlo simulations with realistic urban scenarios to capture both background variability and source signatures. RADAI models construction-material NORM, people and vehicles, urban clutter, and dynamic environmental effects such as cosmic-ray and rain-induced transients, and they provide list-mode detector data with motion and response modeling suitable for algorithm training and evaluation. The RADAI project resulted in three publicly-released complementary datasets together with an online scoring portal for standardized performance assessment and an open software toolkit that supports data access, augmentation, model development, and evaluation. These resources enable reproducible comparisons across methods and promote rigorous studies at the scale required by contemporary machine learning. By grounding algorithm development in realistic, well-documented conditions, RADAI supports progress toward more robust detection, identification, and localization in complex urban environments.

Ghawaly, James M. [Division of Computer Science an

Center for Space Transportation and Applied Research Fifth Annual Technical Symposium Proceedings

This Fifth Annual Technical Symposium, sponsored by the UT-Calspan Center for Space Transportation and Applied Research (CSTAR), is organized to provide an overview of the technical accomplishments of the Center's five Research and Technology focus areas during the past year. These areas include chemical propulsion, electric propulsion, commerical space transportation, computational methods, and laser materials processing. Papers in the area of artificial intelligence/expert systems are also presented.

Source record

Smart Process Planning for Automated Fiber Placement

Many industries, including aerospace, automotive, wind energy, maritime, and sporting goods, rely on strong, lightweight materials called composites. These materials are made by layering fibers, which can come in the form of narrow strips or wider sheets, and setting them in a polymer matrix. One of the most advanced ways to make these parts is through automated fiber placement, where a machine lays down the fibers in precise patterns. This method can create very efficient and strong designs, but it is complex, expensive, and often depends heavily on the experience of skilled engineers. Today, the design, manufacturing, and inspection stages of composite production are usually handled separately. This separation means that important information, such as how a part will be built or what defects might occur, is not always shared between stages. As a result, parts may not be as lightweight, strong, or defect-free as possible, and the process can take longer and cost more. This research develops a smart process planning system that connects design, manufacturing, and inspection into one continuous process. Built as software that works with existing tools, the system can automatically plan how the fibers are placed, predicting and reducing defects while improving both manufacturability and strength. The system optimizes not only individual layers but also how defects are distributed across all layers, preventing them from stacking up in ways that weaken the final part. It also uses inspection results from completed parts to improve future designs, creating a feedback loop where each stage informs the others. The system was tested by designing a composite panel using this new approach and comparing it to a panel made with state-of-the-art manual planning methods. The results showed that the system could intentionally control where defects appeared and increase the efficiency of the planning process. By unifying design, manufacturing, and inspection, this research shows a way to make advanced composite manufacturing more efficient, consistent, and cost-effective. This approach lowers the barrier to using automated fiber placement and opens the door for its wider adoption not only in aerospace but also in industries such as automotive, wind energy, maritime, and sporting goods, where strong and lightweight structures are essential.

Computer-Aided Process Planning

Phase stability in the Hf-N and Zr-N systems

Hf and Zr nitrides are promising compounds for many technologically important areas, including high-temperature structural applications, quantum computing, and solar and optical applications. Here, this article reports on a comprehensive first-principles statistical mechanics study of phase stability in the Hf-N and Zr-N binary systems. A high solubility of nitrogen in the hcp forms of Hf and Zr is predicted. The rocksalt forms of HfN and ZrN can also tolerate a high degree of off-stoichiometry through the introduction of nitrogen and metal vacancies. The Hf-N binary favors a family of stacking faulted parent crystal structures at intermediate nitrogen concentrations that host a unique form of short-range order among nitrogen interstitials and vacancies. These phases can accommodate some degree of configurational entropy and remain ordered to temperatures as high as 1200 K.

Monte Carlo methods

From Rules to Reasoning: A Survey of Large Language Model-Based Approaches to Scientific Hypothesis and Idea Generation

Scientific hypothesis generation represents a fundamental challenge in contemporary research due to exponentially expanding literature volumes and increasing disciplinary specialization. Large language models (LLMs) have emerged as transformative tools for automated scientific discovery, moving beyond traditional rule-based and literature-mining approaches. Four paradigmatic approaches define current LLM-driven hypothesis generation: direct prompting and fine-tuning methods, knowledge-enhanced frameworks integrating retrieval-augmented generation (RAG), multi-agent collaborative systems simulating research teams, and reasoning-focused approaches implementing cognitive architectures. Domain-specific applications demonstrate statistical equivalence to human expert performance in social psychology, experimental validation in biomedical research, and near-expert quality in astronomy. Evaluation methodologies encompass human expert assessment, LLM-as-judge frameworks, and comprehensive benchmarking systems. Technical challenges include hallucination management, knowledge integration limitations, and balancing novelty with feasibility. Future directions emphasize hybrid neural-symbolic architectures and sophisticated human-AI collaboration models for responsible scientific discovery acceleration.

AI-driven discovery

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION