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 667 records · Page 37

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗

An Optimization-Based Framework to Define the Probabilistic Design Space of Pharmaceutical Processes with Model Uncertainty

To increase manufacturing flexibility and system understanding in pharmaceutical development, the FDA launched the quality by design (QbD) initiative. Within QbD, the design space is the multidimensional region (of the input variables and process parameters) where product quality is assured. Given the high cost of extensive experimentation, there is a need for computational methods to estimate the probabilistic design space that considers interactions between critical process parameters and critical quality attributes, as well as model uncertainty. In this paper we propose two algorithms that extend the flexibility test and flexibility index formulations to replace simulation-based analysis and identify the probabilistic design space more efficiently. The effectiveness and computational efficiency of these approaches is shown on a small example and an industrial case study.

60 APPLIED LIFE SCIENCES↗

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗

Water Quality Projects Summary for the Mid-Columbia and Cumberland River Systems

Scheduling and operational control of hydropower systems is accompanied with a keen awareness of the management of water use, environmental effects, and policy, especially within the context of strict water rights policy and generation maximization. This is a multi-objective problem for many hydropower systems, including the Cumberland and Mid-Columbia river systems. Though each of these two systems have distinct operational philosophies, hydrologic characteristics, and system dynamics, they both share a responsibility to effectively manage hydropower and the environment, which requires state-of-the art improvements in the approaches and applications for water quality modeling. The Department of Energy and Oak Ridge National Laboratory have developed tools for total dissolved gas (TDG) prediction on the Mid-Columbia River and a decision-support system used for hydropower generation and environmental optimization on the Cumberland River. In conjunction with IIHR - Hydroscience & Engineering, The University of Iowa and University of Colorado s Center for Advanced Decision Support for Water and Environmental Systems (CADSWES), ORNL has managed the development of a TDG predictive methodology at seven dams along the Mid-Columbia River and has enabled the ability to utilize this methodology for optimization of operations at these projects with the commercially available software package Riverware. ORNL has also managed the collaboration with Vanderbilt University and Lipscomb University to develop a state-of-the art method for reducing high-fidelity water quality modeling results into surrogate models which can be used effectively within the context of optimization efforts to maximize generation for a reservoir system based on environmental and policy constraints. The novel contribution of these efforts is the ability to predict water quality conditions with simplified methodologies at the same level of accuracy as more complex and resource intensive computing methods. These efforts were designed to incorporate well into existing hydropower and reservoir system scheduling models, with runtimes that are comparable to existing software tools. In addition, the transferability of these tools to assess other systems is enhanced due the use of simplistic and easily attainable values for inputs, straight-forward calibration of predictive equation coefficients, and standardized comparison of traditionally familiar outputs.

54 ENVIRONMENTAL SCIENCES↗

Finding Electronic Structure Machine Learning Surrogates without Training

A myriad of phenomena in materials science and chemistry rely on quantum-level simulations of the electronic structure in matter. While moving to larger length and time scales has been a pressing issue for decades, such large-scale electronic structure calculations are still challenging despite modern software approaches and advances in high-performance computing. The silver lining in this regard is the use of machine learning to accelerate electronic structure calculations – this line of research has recently gained growing attention. The grand challenge therein is finding a suitable machine-learning model during a process called hyperparameter optimization. This, however, causes a massive computational overhead in addition to that of data generation. We accelerate the construction of machine-learning surrogate models by roughly two orders of magnitude by circumventing excessive training during the hyperparameter optimization phase. We demonstrate our workflow for Kohn-Sham density functional theory, the most popular computational method in materials science and chemistry.

36 MATERIALS SCIENCE↗

Computational Optimization of Electric Fields for Improving Catalysis of a Designed Kemp Eliminase

Here we report a computational method to improve efficiency of a de novo designed Kemp eliminase enzyme KE15, by identifying mutations that enhance electric fields and chemical positioning of the substrate that contribute to free energy stabilization of the transition state. Starting from the design that has a k cat /K M of 27 M -1 s -1 , the most improved variant introduced four computationally targeted mutations to yield a k cat /K M of 403 M -1 s -1 , with almost all of the enzyme improvement realized through a 43-fold improvement in k cat , indicative of a direct impact on the chemical step. This work raises the prospect of computationally designing enzymes that achieve better efficiency with more minimal experimental intervention using electric field optimization as guidance.

42 ENGINEERING↗

Assessment of genome annotation using gene function similarity within the gene neighborhood

Functional annotation of bacterial genomes is an obligatory and crucially important step of information processing from the genome sequences into cellular mechanisms. However, there is a lack of computational methods to evaluate the quality of functional assignments. Results: We developed a genome-scale model that assigns Bayesian probability to each gene utilizing a known property of functional similarity between neighboring genes in bacteria. Conclusions: Our model clearly distinguished true annotation from random annotation with Bayesian annotation probability >0.95. Our model will provide a useful guide to quantitatively evaluate functional annotation methods and to detect gene sets with reliable annotations.

59 BASIC BIOLOGICAL SCIENCES↗

Simulation methods for liquid–liquid phase separation of disordered proteins

Liquid–liquid phase separation of intrinsically disordered proteins (IDPs) and other biomolecules is a highly complex but robust process used by living systems. Drawing inspiration from biology, phase separating proteins have been successfully utilized for promising applications in fields of materials design and drug delivery. These protein-based materials are advantageous due to the ability to finely tune their stimulus-responsive phase behavior and material properties, and the ability to encode biologically active motifs directly into the sequence. The number of possible protein sequences is virtually endless, which makes sequence-based design a rather daunting task, but also attractive due to the amount of control coming from exploration of this variable space. The use of computational methods in this field of research have come to the aid in several aspects, including interpreting experimental results, identifying important structural features and molecular mechanisms capable of explaining the phase behavior, and ultimately providing predictive frameworks for rational design of protein sequences. We provide an overview of computational studies focused on phase separating biomolecules and the tools that are available to researchers interested in this topic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometrically necessary dislocation fingerprints of dislocation loop absorption at grain boundaries

Here we present a numerical methodology to compute the Nye-tensor fingerprints of dislocation loop absorption at grain boundaries (GBs) for comparison with TEM observations of irradiated polycrystals. Our approach links atomistic simulations of self-interstitial atom (SIA) prismatic loops gliding toward and interacting with GBs in body-centered cubic iron with experimentally extracted geometrically necessary dislocation (GND) maps to facilitate the interpretation of damage processes. The Nye-tensor analysis is strongly mesh-size dependent—corresponding to resolution-dependent TEM observations. The method computes GND fingerprints from discretized dislocation line segments extracted from molecular dynamics simulations of dislocation loops being absorbed at a GB. Specifically, we perform MD simulation of prismatic loops of two diameters and monitor the three stages of the absorption process: loop glide, the partial, and full absorption of the loops at a [1 0 0] symmetric tilt GB. These methods provide a framework for future investigations of the nature of defect absorption by grain boundaries under irradiation conditions.

36 MATERIALS SCIENCE↗

Kernel-based global sensitivity analysis obtained from a single data set

Results from global sensitivity analysis (GSA) often guide the understanding of complicated input–output systems. Kernel-based GSA methods have recently been proposed for their capability of treating a broad scope of complex systems. In this paper, we develop a new set of kernel GSA tools when only a single set of input–output data is available. Three key advances are made: (1) A new numerical estimator is proposed that demonstrates an empirical improvement over previous procedures. (2) A computational method for generating inner statistical functions from a single data set is presented. (3) A theoretical extension is made to define conditional sensitivity indices, which reveal the degree that the inputs carry shared information about the output when inherent input–input correlations are present. Utilizing these conditional sensitivity indices, a decomposition is derived for the output uncertainty based on what is called the optimal learning sequence of the input variables, which remains consistent when correlations exist between the input variables. Further, while these advances cover a range of GSA subjects, a common single data set numerical solution is provided by a technique known as the conditional mean embedding of distributions. The new methodology is implemented on benchmark systems to demonstrate the provided insights.

42 ENGINEERING↗

Vulnerable to misinformation? Verifi!

We present Verifi2, a visual analytic system to support the investigation of misinformation on social media. On the one hand, social media platforms empower individuals and organizations by democratizing the sharing of information. On the other hand, even well-informed and experienced social media users are vulnerable to misinformation. To address the issue, various models and studies have emerged from multiple disciplines to detect and understand the effects of misinformation. However, there is still a lack of intuitive and accessible tools that help social media users distinguish misinformation from verified news. In this paper, we present Verifi2, a visual analytic system that uses state-of-the-art computational methods to highlight salient features from text, social network, and images. By exploring news on a source level through multiple coordinated views in Verifi2, users can interact with the complex dimensions that characterize misinformation and contrast how real and suspicious news outlets differ on these dimensions. To evaluate Verifi2, we conduct interviews with experts in digital media, journalism, education, psychology, and computing who study misinformation. Our interviews show promising potential for Verifi2 to serve as an educational tool on misinformation. Furthermore, our interview results highlight the complexity of the problem of combating misinformation and call for more work from the visualization community.

Karduni, Alireza↗

Knowns and unknowns in metabolomics identified by multidimensional NMR and hybrid MS/NMR methods

Metabolomics continues to make rapid progress through the development of new and better methods and their applications to gain insight into the metabolism of a wide range of different biological systems from a systems biology perspective. Customization of NMR databases and search tools allows the faster and more accurate identification of known metabolites, whereas the identification of unknowns, without a need for extensive purification, requires new strategies to integrate NMR with mass spectrometry, cheminformatics, and computational methods. For some applications, the use of covalent and non-covalent attachments in the form of labeled tags or nanoparticles can significantly reduce the complexity of these tasks.

Bingol, Kerem↗

Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis beyond the relative scales typically employed. In principle, theoretical methods can provide such quantification but are often hindered by the innate complexity of strong electron correlation and dynamic relaxation of solvent environments. Here, we present and validate a protocol for calculating the redox potentials of ferrocene/ferrocenium redox pair in the acetonitrile. Equation-of-motion ionization potential coupled-cluster single-double (EOM-IP-CCSD) and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials (IP) as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrated the ability of modern computational methods to predict absolute redox potentials quantitatively and, more critically, quantify the correlation of dynamic effects, which underlie their origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transplatformer: translating toxicogenomic profiles between generations of platforms

Background Transcriptomic profiling technologies have advanced the analysis of biological and toxicological responses. However, substantial differences in probe design, dynamic range, gene coverage, and preprocessing pipelines across platforms introduce artifacts that limit cross-study integration and hinder the reuse of historical datasets. We aim to develop computational methods for accurate cross-platform translation to maximize the value of legacy resources. Results We present TransPlatformer a deep learning framework for translating gene expression profiles across heterogeneous toxicogenomics platforms. TransPlatformer employs a novel attention-based architecture to map high-dimensional fold-change vectors from legacy microarray technologies to current platforms. Models are trained and evaluated using DrugMatrix, spanning three technological generations. We investigate mixed-tissue, single-tissue, and cross-tissue training paradigms and benchmark performance against multilayer perceptron and matrix-completion baselines. In mixed-tissue training, TransPlatformer achieves a greater than 50% reduction in mean absolute error (0.043 vs. 0.09) and nearly doubles Pearson correlation ( ≈ 0.71 vs. 0.37) relative to baseline methods. Importantly, TransPlatformer preserves rare but biologically meaningful over- and under-expressed signals, with mean absolute error below 0.22. Single-tissue models yield further improvements for well-represented organs, such as a 10% reduction in liver mean absolute error, while underscoring the need for data augmentation strategies in low-sample tissues.ra Conclusions TransPlatformer provides an effective and scalable computational solution for cross-platform transcriptomic translation. By enabling biologically faithful harmonization of gene expression data, the proposed approach facilitates the reuse of legacy toxicogenomics datasets, enhances downstream biomarker discovery, and supports more reproducible predictive modeling in toxicology.

59 BASIC BIOLOGICAL SCIENCES↗

Integrating proteomic or transcriptomic data into metabolic models using linear bound flux balance analysis

Abstract Motivation Transcriptomics and proteomics data have been integrated into constraint-based models to influence flux predictions. However, it has been reported recently for Escherichia coli and Saccharomyces cerevisiae, that model predictions from parsimonious flux balance analysis (pFBA), which does not use expression data, are as good or better than predictions from various algorithms that integrate transcriptomics or proteomics data into constraint-based models. Results In this paper, we describe a novel constraint-based method called Linear Bound Flux Balance Analysis (LBFBA), which uses expression data (either transcriptomic or proteomic) to predict metabolic fluxes. The method uses expression data to place soft constraints on individual fluxes, which can be violated. Parameters in the soft constraints are first estimated from a training expression and flux dataset before being used to predict fluxes from expression data in other conditions. We applied LBFBA to E.coli and S.cerevisiae datasets and found that LBFBA predictions were more accurate than pFBA predictions, with average normalized errors roughly half of those from pFBA. For the first time, we demonstrate a computational method that integrates expression data into constraint-based models and improves quantitative flux predictions over pFBA. Availability and implementation Code is available in the Supplementary data available at Bioinformatics online. Supplementary information Supplementary data are available at Bioinformatics online.

Tian, Mingyuan↗

Direct Simulations of Coupled Transport and Reaction on Nano-Scale X-Ray Computed Tomography Images of Platinum Group Metal-Free Catalyst Cathodes

The nano/micro-scale geometry of polymer electrolyte fuel cell (PEFC) catalyst layers critically affects cell performance. The small length scales and complex structure of these composite layers make it challenging to analyze cell performance and physics at the particle scale by experiment. We present a computational method to simulate transport and chemical reaction phenomena at the pore/particle-scale and apply it to a PEFC cathode with platinum group metal free (PGM-free) catalyst. Here, we numerically solve the governing equations for the physics with heterogeneous oxygen diffusion coefficient and proton conductivity evaluated using the actual electrode structure and ionomer distribution obtained using nano-scale resolution X-ray computed tomography (nano-CT). Using this approach, the oxygen concentration and electrolyte potential distributions imposed by the oxygen reduction reaction are solved and the impact of the catalyst layer structure on performance is evaluated.

36 MATERIALS SCIENCE↗

Subcritical Copper-Reflected-α-phase Plutonium (SCRαP) Integral Experiment Design

A Subcritical Copper-Reflected α-phase Plutonium (SCRαP) integral benchmark experiment is proposed where multiplication will be approximated using correlated neutron data from a detector system consisting of {sup 3}He tubes inside polyethylene. Measurements will be performed on various subcritical experimental configurations consisting of a weapons-grade plutonium sphere surrounded by different Cu thicknesses. In addition to the proposed base experimental configurations with Cu, additional configurations are planned with the plutonium ball nested in various thicknesses of interleaved polyethylene spherical shells mixed in with the Cu shells. The polyethylene is intended to provide fast neutron moderation and reflection, resulting in additional measurements with differing multiplication, spectra, and nuclear data sensitivity. The proposed experiments will performed at the National Criticality Experiments Research Center (NCERC). NCERC is home to a 4.5-kg α-phase stainless steel clad plutonium sphere that is referred to as the BeRP (Beryllium-Reflected Plutonium) ball due to its historical use in a beryllium-reflected critical experiment. In 2014, similar subcritical measurements were performed with the BeRP ball surrounded by nickel and tungsten. The nickel-reflected measurements were documented as a benchmark evaluation which was accepted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP). Similar to past measurements, the proposed work will help both identify deficiencies and quantify uncertainties in nuclear data and validate computational methods related to neutron multiplication inference for subcritical benchmark evaluations. Monte Carlo simulations were performed for the preliminary design of a subcritical Copper-Reflected - phase Plutonium (SCRP) integral benchmark experiment. Specifically, MCNP{sup R} version 6 was used to determine the k{sub eff} and associated nuclear data sensitivities for different preliminary configurations (reflector thickness, material type). The final experimental design process will include continuous-energy nuclear data sensitivity analysis. This will help assess the potential impact of this integral measurement as it relates to international efforts to continuously improve current libraries (ENDF, JEFF, JENDL etc) that still tend to over/under-estimate the measured results of integral experiments, sometimes significantly. Deficiencies in underlying nuclear data quantities such as Nubar have been shown to have an effect on inferred values from subcritical measurements as well. The proposed experiments and subsequent computational validation will help identify such deficiencies. It will also help validate new nuclear data evaluations, including one for Cu that was recently performed in the resolved region up to 100 MeV for {sup 65}Cu and {sup 63}Cu which is expected to improve Cu-related benchmark simulation performance. The final experimental design phase will also include simulation results that will incorporate the detectors that will be used for the subcritical measurements. For the benchmark portion of the measurements, the main detectors that are being considered are portable neutron multiplicity counters made up of {sup 3}He tubes embedded in polyethylene and lined with cadmium. These detectors employ a data acquisition system that produces list-mode data (correlated information where each event is listed sequentially) that can be analyzed using the Feynman variance-to-mean method. Other detector systems including SNAP (single {sup 3}He tube embedded in polyethylene), and a high-purity germanium (HPGe) detector are being considered for supplementary measurements only and are not intended to be used for the final benchmark measurements. Decoupling the supplementary detectors from the system eliminates any contribution to the total system multiplication, and greatly simplifies the final computational benchmark model. Therefore, an optimization of the portable neutron multiplicity counter placement in the benchmark measurements will be performed to ensure that there is minimal contribution to the total system multiplication. In addition, new temperature measurement techniques and hardware will be explored. (authors)

42 ENGINEERING↗

Analysis of buffer-IPyC separation in TRISO fuel particles

During High Temperature Gas-cooled Reactor (HTGR) operation, and due to the neutron irradiation in the reactor core, damage of the nuclear fuel coating layers occurs. The mechanism of damage formation in the TRISO fuel is explored by the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, in which the debonding process between coating layers was also investigated. The purpose of this paper is to report simulation results for two models. Firstly, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted. Secondly, a debonding enabled model, where the gap between those layers is created. The simulations were performed with the Bison code. The inputs of the simulated models are based on data from the AGR-1 experiment. Further, the research included simulations on spherical and aspherical fuel types. The simulations match results obtained by the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the behavior of the gap between buffer and IPyC layers in TRISO fuel. Based on the irradiation experiments, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. In addition, the fuel sphericity does not influence the Weibull failure probability. The results obtained for spherical and aspherical fuel are nearly identical.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗