Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational tools”

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.

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A Wigner Monte Carlo approach to density functional theory

In order to simulate quantum N-body systems, stationary and time-dependent density functional theories rely on the capacity of calculating the single-electron wave-functions of a system from which one obtains the total electron density (Kohn–Sham systems). In this paper, we introduce the use of the Wigner Monte Carlo method in ab-initio calculations. This approach allows time-dependent simulations of chemical systems in the presence of reflective and absorbing boundary conditions. It also enables an intuitive comprehension of chemical systems in terms of the Wigner formalism based on the concept of phase-space. Finally, being based on a Monte Carlo method, it scales very well on parallel machines paving the way towards the time-dependent simulation of very complex molecules. A validation is performed by studying the electron distribution of three different systems, a Lithium atom, a Boron atom and a hydrogenic molecule. For the sake of simplicity, we start from initial conditions not too far from equilibrium and show that the systems reach a stationary regime, as expected (despite no restriction is imposed in the choice of the initial conditions). We also show a good agreement with the standard density functional theory for the hydrogenic molecule. These results demonstrate that the combination of the Wigner Monte Carlo method and Kohn–Sham systems provides a reliable computational tool which could, eventually, be applied to more sophisticated problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2019

Radial feature descriptors for cell classification and recommendation

This paper introduces computational tools for cell classification into normal and abnormal, as well as content-based-image-retrieval (CBIR) for cell recommendation. It also proposes the radial feature descriptors (RFD), which define evenly interspaced segments around the nucleus, and proportional to the convexity of the nuclear boundary. Experiments consider Herlev and CRIC image databases as input to classification via Random Forest and bootstrap; we compare 14 different feature sets by means of False Negative Rate (FNR) and Kappa (κ), obtaining FNR =0.02 and κ=0.89 for Herlev, and FNR =0.14 and κ=0.78 for CRIC. Next, we sort and rank cell images using convolutional neural networks and evaluate performance with the Mean Average Precision (MAP), achieving MAP =0.84 and MAP =0.82 for Herlev and CRIC, respectively. Finally, cell classification show encouraging results regarding RFD, including its sensitivity to intensity variation around the nuclear membrane as it bypasses cytoplasm segmentation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Utilizing machine learning to predict tensile ductility and yield strength of CoNiV-based multi-principal elements alloys

This study explores the use of machine learning (ML) as a computational tool to accelerate the design of multi-principal element alloys (MPEAs) with improved tensile elongation. An ML model was trained using available experimental data from the literature along with theoretically derived features to predict yield strength (YS) and ductility. A subset of ML-predicted compositions—CoNiVFe, CoNiVTi, CoNiVTiFe, and CoCrNiVTi—was synthesized and evaluated through tensile testing. The ML model underpredicted YS by approximately 20–30 % and overpredicted ductility by 60–70 % for Ti-containing alloys. Microstructural analysis revealed that Ti segregation at interdendritic regions contributed to early fracture, leading to discrepancies in ductility predictions. Ti segregation at these regions likely drives the increased YS due to segregation strengthening. In contrast, the CoNiVFe alloy showed good agreement with both experimental YS and elongation, with prediction errors of ∼10.2 % and ∼20.7 %, respectively. Microstructural characterization revealed minimal segregation in this alloy, suggesting that the ML model can reliably predict the properties of alloys with little to no segregation. These findings highlight the capability of ML in predicting YS with good accuracy but underscore its limitations in capturing defect-driven failure mechanisms such as segregation-induced embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An accelerated framework for predicting creep rupture lifetimes in engineering alloys

Confidently predicting high-temperature deformation, including creep and creep rupture, is paramount for the design and commercialization of candidate materials for advanced nuclear energy systems. To accelerate creep quantification, we introduce a framework that enables rapid, cost-effective, and reliable prediction of creep rupture lifetimes, minimizing reliance on time-intensive bulk creep testing. Unlike conventional creep analysis, which requires extensive time and resources, our method leverages a maximum of four short-term bulk creep tests as training data for prediction. This framework combines high-throughput nanoindentation up to 700 °C with these targeted bulk tests to inform our creep rupture model in order to predict rupture lifetimes. The strong agreement between our predictions and conventional experimental data demonstrates the effectiveness of our approach for accelerated creep analysis and lifetime prediction of structural components in high-temperature applications. Our multi-pronged approach motivates further integration of computational tools and advanced instrumentation to establish a universal framework for understanding high-temperature material responses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

GMFOLD: Subgraph matching for high-throughput DNA-aptamer secondary structure classification and machine learning interpretability

Aptamers are oligonucleotide receptors that bind to their targets with high affinity. Here, we consider aptamers comprised of single-stranded DNA that undergo target-binding-induced conformational changes, giving rise to unique secondary and tertiary structures. Given a specific aptamer primary sequence, there are well-established computational tools (notably mfold) to predict the secondary structure via free energy minimization algorithms. While mfold generates secondary structures for individual sequences, there is a need for a high-throughput process whereby thousands of DNA structures can be predicted in real-time for use in an interactive setting, when combined with aptamer selections that generate candidate pools that are too large to be experimentally interrogated. We developed a new Python code for high-throughput aptamer secondary structure determination (GMfold). GMfold uses subgraph matching methods to group aptamer candidates by secondary structure similarities. We also improve an open-source code, SeqFold, to incorporate subgraph matching concepts. We represent each secondary structure as a lowest-energy bipartite subgraph matching of the DNA graph to itself. These new tools enable thousands of DNA sequences to be compared based on their secondary structures, using machine-learning algorithms. This process is advantageous when analyzing sequences that arise from aptamer selections via systematic evolution of ligands by exponential enrichment (SELEX). This work is a building block for future machine-learning-informed DNA-aptamer selection processes to identify aptamers with improved target affinity and selectivity and advance aptamer biosensors and therapeutics.

Aptamer↗
DOE OSTI2018

Unraveling the surface chemistry processes in lithiated and boronized plasma material interfaces under extreme conditions

The review of recent theoretical and experimental research on the complex surface chemistry processes that evolve from low-Z material conditioning on plasma-facing materials under extreme fusion plasma conditions is presented. A combination of multi-scale computational physics and chemistry modeling with real-time diagnosis of the plasma-material interface in tokamak fusion plasma edge is complemented by ex-vessel in-situ single-effect experimental facilities to unravel the evolving characteristics of low-Z components under irradiation. Effects of the lithium and boron coatings at carbon surfaces to the retention of deuterium and chemical sputtering of the plasma-facing surfaces are discussed in detail. The critical role of oxygen in the surface chemistry during hydrogen-fuel irradiation is found to drive the kinetics and dynamics of these surfaces as they interact with fusion edge plasma that ultimately could have profound effects on fusion plasma confinement behavior. Computational studies also extend in spatio-temporal scales not accessible by empirical means and therefore open the opportunity for a strategic approach at irradiation surface science studies that combined these powerful computational tools with in-vessel and ex-vessel in-situ diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Elucidating texture and grain morphology contributions to the micromechanical response of additively manufactured Inconel 625

Microstructural variation of additively manufactured (AM) metal components in comparison to wrought counterparts makes certification for critical applications a challenge. Microscale simulations leveraging modern computational tools may be used to supplement testing of AM microstructures, thus accelerating certification by reducing the number of experiments needed. However, as micromechanical response is closely tied to critical properties like fatigue-life and fracture, utilization of these simulations with macroscale experimental data alone is insufficient. One means to attain microscale experimental data is in situ diffraction data collected from synchrotron X-ray sources. In this work, such data were collected during in situ compression of AM Inconel 625 superalloy. Interpretation of experimental results was assisted by massive (8M element) complementary micromechanical simulations performed on sets of virtual microstructures generated using cellular automata. Together, micromechanical data from diffraction experiments and simulations were used to probe the effects of textured “track” microstructures generated during laser powder bed fusion and directional strength-to-stiffness on micromechanical response. Though fiber-averaged directional strength-to-stiffness ratios were expected to dominate given the high elastic anisotropy of the material, the combination of small variations in texture and specific grain configurations unique to AM microstructures lead to significant variability in micromechanical response after yield. The findings emphasize the importance of high-fidelity microstructural representation that captures key texture components and AM-specific morphology for property prediction of AM metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The carbon challenge: Design, synthesis, and chemisorption behavior of solid sorbents in direct air capture of carbon dioxide

Direct air capture (DAC) of CO 2 is a promising solution for reducing the carbon footprint through "negative emission" technology. However, the low CO 2 concentration (~400 ppm) and the dynamic nature of DAC processes present challenges in designing effective sorbent systems. Recent advancements in material design and structural engineering have led to the development of high-performance solid sorbents, offering a more stable, safe, and energy-efficient alternative to traditional liquid CO 2 capture methods. This review highlights progress in solid sorbent-based DAC, focusing on amine-modified materials, hydrogen-bonded frameworks, and ionic liquid-engineered scaffolds. The discussion covers design principles, synthesis methodologies, and their impact on CO 2 chemisorption, comparing the advantages and limitations of each approach. Characterization techniques, especially operando methods and computational tools, are reviewed to understand sorbent behavior during CO 2 integration and release. The reaction pathways and interaction mechanisms of these sorbents with CO 2 are analyzed to guide future design. Additionally, the CO 2 chemisorption behaviors, including capacity, sorption kinetics, recyclability, and durability in the presence of gaseous impurities and under humid conditions will be evaluated and compared. Further, the review offers unique insights into the physical properties, chemical structures, and surface engineering effects of these sorbents, based on comprehensive characterization and evaluation techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Modeling and simulation functional needs for molten salt reactor licensing

The modeling and simulation functional needs to license and design molten salt reactor concepts were assessed and an evaluation methodology was developed using a licensing-focused approach considering the modeling and simulation components that contribute to a molten salt reactor safety analysis report. These functional needs were prioritized by identifying the relationships between functional needs and using a quantitative weighting method to identify the most critical components to build a licensing case. This novel weighting method provided a list of functional needs prioritized relative to key advanced reactor licensing outcomes. In addition, functional need inputs were ranked to reveal those that feed the most modeling and simulation components. The findings from this prioritization effort were used to inform the modeling and simulation goals of the US Department of Energy Office of Nuclear Energy (DOE-NE) Molten Salt Reactor campaign.Modeling and simulation needs were grouped into three main categories: (1) source term definition, material accountancy, and radiological consequences, (2) safety analysis, and (3) steady-state neutronics and thermal hydraulics. Here, five near-term goals were identified among these focus areas that add value: (1) deploying a systems and safeguards analysis computational tool, (2) defining source terms for molten salt reactor analysis and interfacing with the US Nuclear Regulatory Commission on accident progression analysis, (3) developing point designs for liquid-fueled thermal-spectrum and fast-spectrum molten salt reactors, (4) building safety analysis specifications and reference solutions for Molten Salt Reactor Experiment transients, and (5) interfacing with the thermophysical and thermodynamic salt chemistry databases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

International fuel performance study of fresh fuel experiments for PCMI effects during RIA experiments

This paper presents the results of High-burnup Experiments for Reactivity-initiated Accident (HERA) Modeling & Simulation (M&S) exercise. The HERA project under the Nuclear Energy Agency (NEA) Second Framework for Irradiation Experiments (FIDES-II) program is focused on studying Light Water Reactor (LWR) fuel behavior during Reactivity-Initiated Accident (RIA) conditions. The Part I M&S cases are based on a series of tests in the Transient Reactor Test (TREAT) facility in the United States and the Nuclear Safety Research Reactor (NSRR) in Japan. The purpose of this work is to evaluate the test design to accomplish its goals in establishing clearer understanding of the effects of power pulse width during RIA conditions. Further, the blind predictions using various computational tools have been performed and compared amongst to interpret the behaviors of high burnup fuels during RIA. While many international participants evaluate the thermal–mechanical behavior of fuel rod under different conditions, a considerable scatter of outputs comes out for the cases due to the disparity between codes in predicting mechanical behaviors. In general, however, the results of thermal–mechanical analysis elaborate that nominal design conditions the shorter pulse width tests in NSRR should cause cladding failures while the TREAT tests appear to have more split prediction of failure or not. Furthermore, the sensitivity analysis varying key testing parameters reveals the considerable effect of power pulse width and total energy deposition on prediction of fuel rod failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Two-phase flow measurements and observations in a cooling panel of the reactor cavity cooling system

The analysis of the flow behavior within a cooling panel of a scaled water-cooled reactor cavity cooling system has been performed supported by high-resolution measurements of the void fraction at the risers’ outlet. The results have shown a stable, symmetric distribution of the flow through the risers during subcooled boiling conditions, characterized by small, slow bubbles with average void fraction lower than 0.3 for all risers. As the temperature increased to saturation, flow distribution appeared strongly asymmetrical and unstable, with peaks of void fraction as high as 0.9. Large, faster bubbles were observed during these peaks, immediately followed by flow of subcooled liquid water. Flow stagnation and inversion has also been observed. The high resolution void fraction results, complemented by the measurements of the coolant flow rate and temperatures, not only have further described the behavior of the system under hypothetical accident conditions, but represent a unique set of data to support the validation of system-level codes, and more advanced computational tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Influence of NOx chemistry on the prediction of natural gas end-gas autoignition in CFD engine simulations

Natural gas (NG) represents a promising low-cost/low-emission alternative to diesel fuel when used in high-efficiency internal combustion engines. Advanced combustion strategies utilizing high EGR rates and controlled end-gas autoignition can be implemented with NG to achieve diesel-like efficiencies; however, to support the design of these next-generation NG ICEs, computational tools, including single- and multi-dimensional simulation packages will need to account for the complex chemistry that can occur between the reactive species found in EGR (including NOx) and the fuel. Research has shown that NOx plays an important role in the promotion/inhibition of large hydrocarbon autoignition and when accounted for in CFD engine simulations, can significantly improve the prediction of end-gas autoignition for these fuels. However, reduced NOx-enabled NG mechanisms for use in CFD engine simulations are lacking, and as a result, the influence of NOx chemistry on NG engine operation remains unknown. Here, we analyze the effects of NOx chemistry on the prediction of NG/oxidizer/EGR autoignition and generate a reduced mechanism of a suitable size to be used in engine simulations. Results indicate that NG ignition is sensitive to NOx chemistry, where it was observed that the addition of EGR, which included NOx, promoted NG autoignition. The modified mechanism captured well all trends and closely matched experimentally measured ignition delay times for a wide range of EGR rates and NG compositions. Here, the importance of C2-C3 chemistry is noted, especially for wet NG compositions containing high fractions of ethane and propane. Finally, when utilized in CFD simulations of a Cooperative Fuels Research (CFR) engine, the new reduced mechanism was able to predict the knock onset crank angle (KOCA) to within one crank angle degree of experimental data, a significant improvement compared to previous simulations without NOx chemistry.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Asynchronous Two-Level Checkpointing Scheme for Large-Scale Adjoints in the Spectral-Element Solver Nek5000

Adjoints are an important computational tool for large-scale sensitivity evaluation, uncertainty quantification, and derivative-based optimization. An essential component of their performance is the storage/recomputation balance in which efficient checkpointing methods play a key role. We introduce a novel asynchronous two-level adjoint checkpointing scheme for multistep numerical time discretizations targeted at large-scale numerical simulations. The checkpointing scheme combines bandwidth-limited disk checkpointing and binomial memory checkpointing. Based on assumptions about the target petascale systems, which we later demonstrate to be realistic on the IBM Blue Gene/Q system Mira, we create a model of the expected performance of our checkpointing approach and validate it using the highly scalable Navier-Stokes spectralelement solver Nek5000 on small to moderate subsystems of the Mira supercomputer. In turn, this allows us to predict optimal algorithmic choices when using all of Mira. We also demonstrate that two-level checkpointing is significantly superior to single-level checkpointing when adjoining a large number of time integration steps. To our knowledge, this is the first time two-level checkpointing had been designed, implemented, tuned, and demonstrated on fluid dynamics codes at large scale of 50k+ cores.

Adjoints↗
DOE OSTI2021

Construct design for CRISPR/Cas-based genome editing in plants

CRISPR construct design is a key step in the practice of genome editing, which includes identification of appropriate Cas proteins, design and selection of guide RNAs (gRNAs), and selection of regulatory elements to express gRNAs and Cas proteins. Here, we review the choices of CRISPR-based genome editors suited for different needs in plant genome editing applications. We consider the technical aspects of gRNA design and the associated computational tools. We also discuss strategies for the design of multiplex CRISPR constructs for high-throughput manipulation of complex biological processes or polygenic traits. We provide recommendations for different elements of CRISPR constructs and discuss the remaining challenges of CRISPR construct optimization in plant genome editing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

HEPOM: Using Graph Neural Networks for the Accelerated Predictions of Hydrolysis Free Energies in Different pH Conditions

Hydrolysis is a fundamental family of chemical reactions where water facilitates the cleavage of bonds. The process is ubiquitous in biological and chemical systems, owing to water’s remarkable versatility as a solvent. However, accurately predicting the feasibility of hydrolysis through computational techniques is a difficult task, as subtle changes in reactant structure like heteroatom substitutions or neighboring functional groups can influence the reaction outcome. Furthermore, hydrolysis is sensitive to the pH of the aqueous medium, and the same reaction can have different reaction properties at different pH conditions. In this work, we have combined reaction templates and high-throughput ab initio calculations to construct a diverse data set of hydrolysis free energies. The developed framework automatically identifies reaction centers, generates hydrolysis products, and utilizes a trained graph neural network (GNN) model to predict ΔG values for all potential hydrolysis reactions in a given molecule. The long-term goal of the work is to develop a data-driven, computational tool for high-throughput screening of pH-specific hydrolytic stability and the rapid prediction of reaction products, which can then be applied in a wide array of applications including chemical recycling of polymers and ion-conducting membranes for clean energy generation and storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗