Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational methods”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Can Electrical Resistance Tomography be used for imaging unsaturated moisture flow in cement-based materials with discrete cracks?

Previously, it has been shown that Electrical Resistance Tomography (ERT) can be used for monitoring moisture flow in undamaged cement-based materials. In this work, we investigate whether ERT could be used for imaging three-dimensional (3D) unsaturated moisture flow in cement-based materials that contain discrete cracks. Novel computational methods based on the so-called absolute imaging framework are developed and used in ERT image reconstructions, aiming at a better tolerance of the reconstructed images with respect to the complexity of the conductivity distribution in cracked material. ERT is first tested using specimens with physically simulated cracks of known geometries, and corroborated with numerical simulations of unsaturated moisture flow. Next, specimens with loading-induced cracks are imaged; here, ERT reconstructions are evaluated qualitatively based on visual observations and known properties of unsaturated moisture flow. Results indicate that ERT is a viable method of visualizing 3D unsaturated moisture flow in cement-based materials with discrete cracks. - Highlights: • 3D EIT is developed to visualize water ingress in cracked mortar. • Mortar with different size discrete cracks are used. • The EIT results are corroborated with numerical simulations. • EIT results accurately show the temporal and spatial variation of water content. • EIT is shown to be a viable method to monitor flow in cracks and matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A computationally efficient algorithm for computing convex hull prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may be required to execute side payments to participants whose costs are not covered through energy sales as determined via traditional locational marginal pricing schemes. Convex hull pricing minimizes this and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. Computing convex hull prices involves solving either a large-scale linear program or the Lagrangian dual of the corresponding non-convex scheduling problem. Further, the former approach requires explicit descriptions of market participants' convex hulls. While linear programs for computing convex hull prices are large, their structure is naturally decomposable by generators. Here, in this work, we propose and empirically analyze a Benders decomposition approach to computing convex hull prices that leverages recent advances in convex hull formulations for thermal generating units. We demonstrate across a large set of test instances that our decomposition approach only requires modest computational effort, obtaining solutions at least an order of magnitude faster than the equivalent large-scale linear programming approach. Overall, we provide a computationally feasible method for computing convex hull prices for industrial scale market clearing problems, enabling the possibility of practical adoption of this advanced pricing mechanism.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2012

Interstitial and substitutional zirconium in SrTiO3

We investigate Zr in SrTiO3 (STO) as a model for nuclear waste forms in which the fission product Sr90 eventually decays to stable Zr through beta emission. The transformation of a divalent into a tetravalent constituent is expected to affect the long-term structural and chemical stability of this solid. Computational methods of electronic structure theory, specifically the density functional theory (DFT) within the supercell model, are used to predict the thermodynamic stability and electronic states of interstitial and Sr- or Ti-substituted Zr atoms in the STO lattice. Native defects such as vacancies and antisites are also considered. When Zr replaces Sr, its most stable configuration is to simply occupy the Sr site (instead of, for example, replacing a Ti and displacing the Ti to the Sr site.) For Zr added to the lattice, its most stable configuration is to replace a Ti, making a ZrTi impurity plus a Ti interstitial (as opposed to the Zr just remaining as an interstitial atom.) ZrSr is predicted to be a double electron donor, ZrTi is electrically inactive and interstitial Zr and Ti are predicted to be quadruple donors, with all donor levels in the conduction band. ZrSr and the tetravalent interstitials are all predicted to increase the crystal volume, and the interstitials also are predicted to lead to a tetragonal distortion of the lattice. Experiments with injection of Zr atoms into STO qualitatively confirm these predictions of crystal structural changes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A quantitative lens on anaerobic life: leveraging the state-of-the-art fluxomics approach to explore clostridial metabolism

Quantitative understanding of clostridial metabolism is of longstanding interest due to the importance of Clostridia as model anaerobes, biotechnology workhorses, and contributors to evolutionary history and ecosystem. Current computational methods such as flux balance analysis-based construction of clostridial metabolism in genome scale provide a fundamental framework for metabolic analysis. However, this method alone is inadequate to characterize cellular metabolic activity. Experiment-driven approaches including isotope tracer-based fluxomics in association with genetic and biochemical methods are needed to gain a more comprehensive understanding. Here we focus on typical examples where these integrated approaches have contributed to the identification of new metabolic pathways and quantification of metabolic fluxes in Clostridia. Here, we also highlight the opportunities and challenges of cutting-edge fluxomics approaches such as machine learning modeling, deuterium tracer approach, and high throughput flux phenotyping in exploring clostridial metabolism with respect to inorganic carbon utilization, redox cofactor interconversion, and other key metabolic features.

09 BIOMASS FUELS↗
DOE OSTI2023

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used in industrial production for biofuels and bioproducts. Developing strains with better performance in biomanufacturing requires systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated in several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. In this report we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test- learn cycles, ultimately advancing the field of fungal biomanufacturing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2020

Thermally induced cracking on the massive concrete structure of the NSLS II synchrotron and its engineering remediation

Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.

42 ENGINEERING↗
DOE OSTI2019

Counter cross-flow evaporator geometries for supercritical organic Rankine cycles

With recent advancements in advanced manufacturing techniques, it is now possible to fabricate complex geometries that take advantage of well known principles of heat transfer. Therefore, unconventional configurations to enhance effectiveness beyond conventional designs can now be considered for practical application. Thermal performance and overall cost of a new design of heat exchangers in counter cross-flow configurations are studied using a simplified but accurate computational method. The new heat exchanger design was introduced and studied previously for a cross-flow configuration by Sabau et al. (2016, 2018). This new design concept uses multi-scale configurations with successive plenums for the working fluid. At the smallest scale the tubes are sized to be equal to the hydraulic entrance length of the inside fluid, in accord with constructal design. With this being said, results indicate that compared to the earlier cross-flow configuration, the counter cross-flow arrangement improves the thermal performance of the heat exchanger by as much as 17% and lowers the total cost by as much as 14%.

42 ENGINEERING↗
DOE OSTI2022

Opportunities and data requirements for data-driven prognostics and health management in liquid hydrogen storage systems

During the past decade, Prognostics and Health Management (PHM) has become an important set of tools in various areas of industry and academic reliability engineering. PHM consists of a variety of mathematical and computational methods used to support data-driven decision-making to increase the safety, availability, and reliability of complex engineering systems. In particular, PHM can provide crucial insight into reliability and safety design improvements for developing technologies where historical performance and failure data are limited. This is the case of hydrogen fueling and storage technologies. This work presents a high-level approach for designing data-driven PHM applications for bulk liquid hydrogen (LH 2 ) storage systems for hydrogen fueling stations. This paper addresses core aspects of the design, development, and implementation of data-driven PHM applications that can improve the reliability assessment of hydrogen components. The analysis focuses on the relationship between data availability and diagnostic/prognostic capabilities; potential challenges; and integration schemes for current risk mitigation measures. We identify potential condition-monitoring data sources for key components in an LH 2 storage system, including storage tanks, piping, and pumps. Further, we determine that the short-term goals for the implementation of data-driven models in PHM frameworks in hydrogen systems should focus on developing adequate data collection and analysis strategies, as well as exploring the effect on reliability, safety, and regulations for hydrogen systems.

08 HYDROGEN↗
DOE OSTI2018

A continuum dislocation-based model of wedge microindentation of single crystals

Recent Electron Backscatter Diffraction (EBSD) experiments have revealed the emergence of heterogeneous dislocation microstructures forming under a wedge indenter in fcc crystals, where micro-meter dislocation patterns challenge the predictions of traditional models of plasticity. In order to explain the formation of these features and develop a relationship between the force-displacement curve and the dislocation substructure, here we present a model of wedge indentation based on the continuum theory of dislocations. The model accounts for large deformation kinematics through the multiplicative split of the deformation gradient tensor, where the incompatible plastic component of deformation results from the flux of dislocations on different and interacting slips systems. Constitutive equations for dislocation fluxes are determined from a dissipative variational principle. As a result, each dislocation density satisfies an initial-boundary value problem with convective-diffusive character, which is coupled to the macroscopic stress and displacement fields governing the deformation process. Solution to the self-consistent continuum formulation is found using the finite element method. Computer simulations mimic the experimental conditions of wedge micro-indentation experiments into Ni single-crystals used by Kysar et al. (2010a). A comparison of overall dislocation density distribution and macroscopic mechanical response shows good overall agreement with the experimental results in terms of the detailed features of dislocation patterns and lattice rotations as well as the macroscopic force-displacement response.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Reliable thermodynamic estimators for screening caloric materials

Reversible, diffusionless, first-order solid-solid phase transitions accompanied by caloric effects are critical for applications in the solid-state cooling and heat-pumping devices. Accelerated discovery of caloric materials requires reliable but faster estimators for predictions and high-throughput screening of system-specific dominant caloric contributions. We assess reliability of the computational methods that provide thermodynamic properties in relevant solid phases at or near a phase transition. Here, we test the methods using the well-studied B2 FeRh alloy as a “fruit fly” in such a materials genome discovery, as it exhibits a metamagnetic transition which generates multicaloric (magneto-, elasto-, and baro-caloric) responses. For lattice entropy contributions, we find that the commonly-used linear-response and small-displacement phonon methods are invalid near instabilities that arise from the anharmonicity of atomic potentials, and we offer a more reliable and precise method for calculating lattice entropy at a fixed temperature. Then, we apply a set of reliable methods and estimators to the metamagnetic transition in FeRh (predicted 346 ± 12 K, observed 353 ± 1 K) and calculate the associated caloric properties, such as isothermal entropy and isentropic temperature changes.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Adsorption mechanism of alkyl hydroxamic acid onto bastnäsite: Fundamental steps toward rational collector design for rare earth elements

Rare earth element (REE) production is limited in part by inefficient strategies for beneficiation. Hydroxamic acid ligands are promising reagents for the selective flotation of bastnäsite [(Ce,La)FCO 3 ], a major REE ore mineral, but the mechanism and energetics of adsorption are not understood, interfering with the design of new, more efficient reagents. In this work, the adsorption of octyl hydroxamic acid onto bastnäsite was calculated using a combination of experimental and computational methods. In-situ attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy revealed changes in the hydroxamate functional group vibrational frequencies, corresponding to chelation with cerium cations at the bastnäsite surface. The results indicate a monodentate chemisorption mechanism at low surface loading that changes to bidentate chemisorption at higher concentrations. This interpretation is supported by molecular vibrational frequency shifts calculated using density functional theory (DFT), and orientation of the hydrocarbon chain measured by sum frequency generation (SFG) vibrational spectroscopy. The binding enthalpies of octyl hydroxamic acid interacting with La and Ce-bastnäsite surfaces were measured using isothermal titration calorimetry (ITC) revealing a stronger coordinating ability with bastnäsite than with a common gangue mineral, calcite (CaCO 3 ). Because octyl hydroxamate favors monodentate adsorption at low surface coverages, the relative chelating strength of metal ions could be a poor predictor for selectivity under monolayer adsorption conditions. At higher surface loadings, where the bidentate mode of adsorption is active, selectivity is likely to be limited by increased flotation of gangue ores.

58 GEOSCIENCES↗
DOE OSTI2015

Self-consistent inclusion of classical large-angle Coulomb collisions in plasma Monte Carlo simulations

Large-angle Coulomb collisions allow for the exchange of a significant proportion of the energy of a particle in a single collision, but are not included in models of plasmas based on fluids, the Vlasov–Fokker–Planck equation, or currently available plasma Monte Carlo techniques. Their unique effects include the creation of fast ‘knock-on’ ions, which may be more likely to undergo certain reactions, and distortions to ion distribution functions relative to what is predicted by small-angle collision only theories. We present a computational method which uses Monte Carlo techniques to include the effects of large-angle Coulomb collisions in plasmas and which self-consistently evolves distribution functions according to the creation of knock-on ions of any generation. The method is used to demonstrate ion distribution function distortions in an inertial confinement fusion (ICF) relevant scenario of the slowing of fusion products.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code

Numerical simulations form an indispensable tool to understand the behavior of a hot plasma that is created inside a tokamak for providing nuclear fusion energy. Various aspects of tokamak plasmas have been successfully studied through the reduced magnetohydrodynamic (MHD) model. The need for more complete modeling through the full MHD equations is addressed here. Our computational method is presented along with measures against possible problems regarding pollution, stability, and regularity. The problem of ensuring continuity of solutions in the center of a polar grid is addressed in the context of a finite element discretization of the full MHD equations. A rigorous and generally applicable solution is proposed here. Useful analytical test cases are devised to verify the correct implementation of the momentum and induction equation, the hyperdiffusive terms, and the accuracy with which highly anisotropic diffusion can be simulated. A striking observation is that highly anisotropic diffusion can be treated with the same order of accuracy as isotropic diffusion, even on non-aligned grids, as long as these grids are generated with sufficient care. This property is shown to be associated with our use of a magnetic vector potential to describe the magnetic field. Several well-known instabilities are simulated to demonstrate the capabilities of the new method. The linear growth rate of an internal kink mode and a tearing mode are benchmarked against the results of a linear MHD code. The evolution of a tearing mode and the resulting magnetic islands are simulated well into the nonlinear regime. The results are compared with predictions from the reduced MHD model. Finally, a simulation of a ballooning mode illustrates the possibility to use our method as an ideal MHD method without the need to add any physical dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Disparate Degrees of Hypervariable Loop Flexibility Control T-Cell Receptor Cross-Reactivity, Specificity, and Binding Mechanism

{alpha}{beta} T-cell receptors (TCRs) recognize multiple antigenic peptides bound and presented by major histocompatibility complex molecules. TCR cross-reactivity has been attributed in part to the flexibility of TCR complementarity-determining region (CDR) loops, yet there have been limited direct studies of loop dynamics to determine the extent of its role. Here we studied the flexibility of the binding loops of the {alpha}{beta} TCR A6 using crystallographic, spectroscopic, and computational methods. A significant role for flexibility in binding and cross-reactivity was indicated only for the CDR3{alpha} and CDR3{beta} hypervariable loops. Examination of the energy landscapes of these two loops indicated that CDR3{beta} possesses a broad, smooth energy landscape, leading to rapid sampling in the free TCR of a range of conformations compatible with different ligands. The landscape for CDR3{alpha} is more rugged, resulting in more limited conformational sampling that leads to specificity for a reduced set of peptides as well as the major histocompatibility complex protein. In addition to informing on the mechanisms of cross-reactivity and specificity, the energy landscapes of the two loops indicate a complex mechanism for TCR binding, incorporating elements of both conformational selection and induced fit in a manner that blends features of popular models for TCR recognition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Predicting metal-binding proteins and structures through integration of evolutionary-scale and physics-based modeling

Metals are essential elements in all living organisms, binding to approximately 50% of proteins. They serve to stabilize proteins, catalyze reactions, regulate activities, and fulfill various physiological and pathological functions. While there have been many advancements in determining the structures of protein-metal complexes, numerous metal-binding proteins still need to be identified through computational methods and validated through experiments. Here, to address this need, we have developed the ESMBind workflow, which combines evolutionary scale modeling (ESM) for metal-binding prediction and physics-based protein-metal modeling. Our approach utilizes the ESM-2 and ESM-IF models to predict metal-binding probability at the residue level. In addition, we have designed a metal-placement method and energy minimization technique to generate detailed 3D structures of protein-metal complexes. Our workflow outperforms other models in terms of residue and 3D-level predictions. To demonstrate its effectiveness, we applied the workflow to 142 uncharacterized fungal pathogen proteins and predicted metal-binding proteins involved in fungal infection and virulence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Type IV failure in weldment of creep resistant ferritic alloys: I. Micromechanical origin of creep strain localization in the heat affected zone

Creep strength enhanced ferritic (CSEF) steels containing 9-12wt% chromium have been extensively used in fossil-fuel-fired power plants. Despite their excellent creep resistance at high temperatures, premature failures (especially Type IV cracking) are often found in the fine-grained heat affected zone (HAZ) or intercritical HAZ of the welded components. This failure mode is preceded by the strain localization in the HAZ, as measured by the Digital Image Correlation (DIC) technique. The present work aims to develop a finite-element based computational method to determine the micromechanical and microstructural origin of the strain localization phenomenon. We construct a two-dimensional digital microstructure based on the actual microstructure of ferritic steel weldments by using the Voronoi-tessellation method, to account for the effects of its large grain-size gradients. A mechanism-based finite element method is developed for modeling the high temperature deformation resulting from a synergy of thermally activated dislocation movements, diffusional flow and grain boundary sliding. The numerical results agree well with the strain measurements by our DIC technique, particularly revealing the effect of pre-welding tempering on the evolution of strain localization in HAZ of creep resistant steel weldments. It is found that the diffusional creep with dependence on grain sizes, dislocation creep with dependence on material strength, and more importantly, grain boundary sliding, contribute synergistically to the creep strain accumulation in the HAZ, and their relative degree of significance is quantified. The creep rupture life will be investigated in the companion paper.

36 MATERIALS SCIENCE↗