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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integral Nuclear Data and Benchmarking Needs for Fusion Energy Systems

Fusion energy systems are currently being designed and optimized using radiation transport codes. To deal with the unique environment inside a fusion-based system, many of these designs incorporate novel materials able to withstand the high radiation fields, ensure adequate cooling and thermal protection, and produce tritium. Validation plays a vital role in building trust in the predictive power of these models and computational methods. Validation of a code consists of modeling documented real-world experiments and comparing the code-predicted response to the measured response. Adequate validation requires measured responses from real-world experiments, also known as integral data, that mimic the system being designed, including materials, impinging radiation, and temperature, among other variables. The most trusted integral data are experimental responses that have been through a rigorous benchmarking process that develops a recommended computational model and evaluates all experimental uncertainties. Finally, there are a few research groups around the world that have been producing integral data for fusion applications, but a substantial investment is needed to address the unique validation needs of the fusion community.

Fusion↗
DOE OSTI2019

A Logical Model of Homology for Comparative Biology

There is a growing body of research on the evolution of anatomy in a wide variety of organisms. Discoveries in this field could be greatly accelerated by computational methods and resources that enable these findings to be compared across different studies and different organisms and linked with the genes responsible for anatomical modifications. Homology is a key concept in comparative anatomy; two important types are historical homology (the similarity of organisms due to common ancestry) and serial homology (the similarity of repeated structures within an organism). We explored how to most effectively represent historical and serial homology across anatomical structures to facilitate computational reasoning. We assembled a collection of homology assertions from the literature with a set of taxon phenotypes for the skeletal elements of vertebrate fins and limbs from the Phenoscape Knowledgebase. Using seven competency questions, we evaluated the reasoning ramifications of two logical models: the Reciprocal Existential Axioms (REA) homology model and the Ancestral Value Axioms (AVA) homology model. The AVA model returned all user-expected results in addition to the search term and any of its subclasses. The AVA model also returns any superclass of the query term in which a homology relationship has been asserted. The REA model returned the user-expected results for five out of seven queries. We identify some challenges of implementing complete homology queries due to limitations of OWL reasoning. This work lays the foundation for homology reasoning to be incorporated into other ontology-based tools, such as those that enable synthetic supermatrix construction and candidate gene discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

HUNTRESS: a fast heuristic for reconstructing phylogenetic trees of tumor evolution (HUNTRESS) v0.1

We introduce HUNTRESS (Histogrammed UNion Tree REconStruction heuriStic), a computational method for tumor phylogeny reconstruction from noisy genotype matrices derived from single-cell sequencing data, whose running time is linear with the number of cells and quadratic with the number of mutations. Provided that the input genotype matrix includes no false positives, each cellular subpopulation is at least a user defined fraction of the total number of cells, and the number of cells are much bigger than the number of mutations considered, HUNTRESS computes the ground truth tumor phylogeny with high probability. On simulated data HUNTRESS is faster than available alternatives with comparable or better accuracy. Additionally, the phylogenies reconstructed by HUNTRESS on two single-cell sequencing data sets agree with the best known evolutionary scenarios for the associated tumors.

Buluc, Aydin↗
DOE OSTI2021

XCP- Computational Summer Workshop: Equation of State Overview [Slides]

An overview of Equation of State (EOS) theory, some relevant physics, experiments, computational methods, and the SESAME database will be presented. The goal will be to describe what types of data and modeling are required to make an EOS that typically spans up to ~10 9 GPa in pressure and 10 9 K in temperature.

42 ENGINEERING↗
DOE OSTI2012

Benchmark On Sensitivity Calculation (Phase III)

The sensitivities of the keff eigenvalue to neutron cross sections have become commonly used in similarity studies and as part of the validation algorithm for criticality safety assessments. To test calculations of the sensitivity coefficients, a benchmark study (Phase III) has been established by the OECD-NEA/WPNCS/EG UACSA (Expert Group on Uncertainty Analysis for Criticality Safety Assessment). This paper presents some sensitivity results generated by the benchmark participants using various computational tools based upon different computational methods: SCALE/TSUNAMI-3D and -1D, MONK, APOLLO2-MORET 5, DRAGON-SUSD3D and MMKKENO. The study demonstrates the performance of the tools. It also illustrates how model simplifications impact the sensitivity results and demonstrates the importance of 'implicit' (self-shielding) sensitivities. This work has been a useful step towards verification of the existing and developed sensitivity analysis methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Quantitative imaging and automated fuel pin identification for passive gamma emission tomography

Compliance of member States to the Treaty on the Non-Proliferation of Nuclear Weapons is monitored through nuclear safeguards. The Passive Gamma Emission Tomography (PGET) system is a novel instrument developed within the framework of the International Atomic Energy Agency (IAEA) project JNT 1510, which included the European Commission, Finland, Hungary and Sweden. The PGET is used for the verification of spent nuclear fuel stored in water pools. Advanced image reconstruction techniques are crucial for obtaining high-quality cross-sectional images of the spent-fuel bundle to allow inspectors of the IAEA to monitor nuclear material and promptly identify its diversion. In this work, we have developed a software suite to accurately reconstruct the spent-fuel cross sectional image, automatically identify present fuel rods, and estimate their activity. Unique image reconstruction challenges are posed by the measurement of spent fuel, due to its high activity and the self-attenuation. While the former is mitigated by detector physical collimation, we implemented a linear forward model to model the detector responses to the fuel rods inside the PGET, to account for the latter. The image reconstruction is performed by solving a regularized linear inverse problem using the fast-iterative shrinkage-thresholding algorithm. We have also implemented the traditional filtered back projection (FBP) method based on the inverse Radon transform for comparison and applied both methods to reconstruct images of simulated mockup fuel assemblies. Higher image resolution and fewer reconstruction artifacts were obtained with the inverse-problem approach, with the mean-square-error reduced by 50%, and the structural-similarity improved by 200%. We then used a convolutional neural network (CNN) to automatically identify the bundle type and extract the pin locations from the images; the estimated activity levels finally being compared with the ground truth. The proposed computational methods accurately estimated the activity levels of the present pins, with an associated uncertainty of approximately 5%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

ManiNetCluster: a novel manifold learning approach to reveal the functional links between gene networks

BACKGROUND:The coordination of genomic functions is a critical and complex process across biological systems such as phenotypes or states (e.g., time, disease, organism, environmental perturbation). Understanding how the complexity of genomic function relates to these states remains a challenge. To address this, we have developed a novel computational method, ManiNetCluster, which simultaneously aligns and clusters gene networks (e.g., co-expression) to systematically reveal the links of genomic function between different conditions. Specifically, ManiNetCluster employs manifold learning to uncover and match local and non-linear structures among networks, and identifies cross-network functional links. RESULTS:We demonstrated that ManiNetCluster better aligns the orthologous genes from their developmental expression profiles across model organisms than state-of-the-art methods (p-value <2.2×10-16). This indicates the potential non-linear interactions of evolutionarily conserved genes across species in development. Furthermore, we applied ManiNetCluster to time series transcriptome data measured in the green alga Chlamydomonas reinhardtii to discover the genomic functions linking various metabolic processes between the light and dark periods of a diurnally cycling culture. We identified a number of genes putatively regulating processes across each lighting regime. CONCLUSIONS:ManiNetCluster provides a novel computational tool to uncover the genes linking various functions from different networks, providing new insight on how gene functions coordinate across different conditions. ManiNetCluster is publicly available as an R package at https://github.com/daifengwanglab/ManiNetCluster.

Comparative network analysis↗
DOE OSTI2012

IC Workflow Project (Final Report)

The Workflow Project is a joint research task between the CASL Advanced Modeling Applications (AMA) and the Virtual Reactor Integration (VRI) focus areas, and the Industry Council. The purpose of the project is to solicit the expertise of the IC member organizations to inform CASL staff and provide feedback for the development of CASL products. The goal is for CASL staff to understand existing workflows at end-user organizations in order to achieve compatibility of CASL products with current needs and future plans. CASL is committed to developing and delivering advanced capabilities and solutions that will allow the nuclear power industry to solve complex problems and improve the operating performance and efficiency of light water reactors. Through the R&D being performed by CASL, advanced computational methods and software is being developed and embodied in the Virtual Environment for Reactor Applications (VERA). The information collected through this project will provide additional considerations that must be included in the development of VERA to ensure that it meets the needs of the nuclear industry. As will be discussed below, key additional areas that must be considered is the usability of the system in an engineering environment including the user input, quality assurance, output of results, computational runtime, and integration into existing engineering analysis processes. This document contains a summary of information gathered through a series of interviews and meetings with IC organizations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Stochastic characterization of phase detection algorithms in phase-shifting interferometry

Phase-shifting interferometry (PSI) is the preferred non-contact method for profiling sub-nanometer surfaces. Based on monochromatic light interference, the method computes the surface profile from a set of interferograms collected at separate stepping positions. Errors in the estimated profile are introduced when these positions are not located correctly. In order to cope with this problem, various algorithms that minimize the effects of certain types of stepping errors (linear, sinusoidal, etc.) have been developed. Despite the relatively large number of algorithms suggested in the literature, there is no unified way of characterizing their performance when additional unaccounted random errors are present. Here, we suggest a procedure for quantifying the expected behavior of each algorithm in the presence of independent and identically distributed (i.i.d.) random stepping errors, which can occur in addition to the systematic errors for which the algorithm has been designed. As a result, the usefulness of this method derives from the fact that it can guide the selection of the best algorithm for specific measurement situations.

42 ENGINEERING↗
DOE OSTI2022

Monolithic Multigrid for a Reduced-Quadrature Discretization of Poroelasticity

Advanced finite-element discretizations and preconditioners for models of poroelasticity have attracted significant attention in recent years. The equations of poroelasticity offer significant challenges in both areas, due to the potentially strong coupling between unknowns in the system, saddle-point structure, and the need to account for wide ranges of parameter values, including limiting behavior such as incompressible elasticity. This paper was motivated by an attempt to develop monolithic multigrid preconditioners for the discretization developed in [C. Rodrigo et al., Comput. Methods App. Mech. Engrg, 341 (2018), pp. 467--484]; we show here why this is a difficult task and, as a result, we modify the discretization in [Rodrigo et al.] through the use of a reduced-quadrature approximation, yielding a more “solver-friendly” discretization. Local Fourier analysis is used to optimize parameters in the resulting monolithic multigrid method, allowing a fair comparison between the performance and costs of methods based on Vanka and Braess--Sarazin relaxation. Further, numerical results are presented to validate the local Fourier analysis predictions and demonstrate efficiency of the algorithms. Finally, a comparison to existing block-factorization preconditioners is also given.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Evaluation of Plutonium Hemisphere Critical Experiments Partially Reflected by Steel and Oil

A series of 15 critical experiments performed at the Rocky Flats Critical Mass Laboratory in the late 1960s were evaluated and then determined to represent acceptable benchmark experiments for the validation of calculational methods. This series of experiments was part of a larger set of experiments performed to evaluate operational safety margins at the Rocky Flats Plant. The experiments consisted of bare plutonium metal hemishells reflected by steel hemishells of increasing thickness and motor oil. The hemishell assembly was suspended within dual aluminum tanks. Criticality was achieved by pumping oil into the tanks such that effectively infinite reflection was achieved in all directions except directly above the assembly; then the critical oil height was recorded. The results of these experiments had been initially ignored because early computational methods had been inadequate to analyze partially-reflected configurations. The dominant uncertainties include the uncertainty in the average plutonium density and the composition of materials in the gaps between the plutonium hemishells. Simple and detailed benchmark models were developed. Eigenvalue calculations using MCNP5 and ENDF/B-VII.0 were within 2s of the benchmark values. This benchmark evaluation has been added to the ICSBEP Handbook.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The Steel Equivalency Workbook: An X-Ray Transmission Calculator

The Steel Equivalency Spreadsheet was created as a robust, user-friendly method for the following: 1. Determining if it is possible to image an object. 2. Reducing time needed when selecting equipment and preparing for field imaging activities. 3. Reducing the amount of equipment taken for field imaging activities. 4. Determining a starting point for exposure settings prior to imaging an object. 5. Reducing the amount of dose deposited to complete an imaging operation (in the spirit of As Low As Reasonably Achievable, ALARA). The spreadsheet was developed using various physical models to account for different phenomena. Future implementations aim to expand beyond the Digital Radiography and Computed Tomography Single Munition Scanner (DRCT SMS) in standard configuration to include high energy XRGs (Betatrons) and sub-MeV pulsed XRGs (XRS4). Alternative computational methods that compensate for incoherent scattering when calculating relative transmission are also being pursued. Ultimately the spreadsheet exceeded the developmental goal of having less than 10% average error when comparing calculations to real-world images.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Molecular engineering of fluorescein dyes as complementary absorbers in dye co-sensitized solar cells

Fluorescein dye derivatives exhibit extended optical absorption up to 500 nm, rendering these compounds suitable as co-absorbers in dye-sensitized solar cells (DSCs). A molecular engineering approach is presented, which embraces this intrinsic optical attribute of fluoresceins, while modifying the dye chemistry to enhance their light harvesting efficiency, in order to effectively tailor them for DSC applications. This approach first realizes relationships between the molecular structure and the optoelectronic properties for a series of five a priori known (parent) fluorescein dyes: 5-carboxyfluorescein (1), a mixture of m -carboxyfluorescein where m = 5 or 6 (2), 5-carboxyfluorescein diacetate (3), 6-carboxyfluorescein diacetate (4), a mixture of n -carboxy-2',7'-dichlorofluorescein diacetate where n = 5 or 6 (5). The first step in this approach combines, where available, experimental and computational methods so that electronic structure calculations can also be validated for representative fluorescein dyes. Such calculations can then be used reliably to predict the structure and properties of fluorescein dyes for cases where experimental data are lacking. Structure-function relationships established from this initial step inform the selection of parent dye 1 that is taken forward to the second step in molecular engineering: in silico chemical derivation to re-functionalize 1 for DSC applications. For this purpose, computational calculations are used to extend the charge conjugation in 1 between its donor and acceptor moieties. These structural modifications result in a bathochromic shift of the lowest excitation by ~1.3-1.9 eV (100-170 nm), making the dye optically absorb in the visible region. Further calculations on dye molecules adsorbed onto the surface of a TiO 2 cluster are used to investigate the dye sensitization behavior via dye adsorption energies and anchoring modes. The results of this theoretical investigation lead to two molecularly engineered fluoresceins being proposed to act as co-sensitizers together with a rhodamine dye. This combination of three dyes ensures chemical compatibility, panchromatic absorption, and restores optical absorption dipping otherwise observed in a DSC device at ~350-400 nm owing to the I-/I- 3 electrolyte. Altogether, the results of this study demonstrate that molecular engineering can be used to identify suitable chemical modifications for organic dyes with improved light harvesting properties for photovoltaic applications.

14 SOLAR ENERGY↗
DOE OSTI2014

Data free inference with processed data products

Here, we consider the context of probabilistic inference of model parameters given error bars or confidence intervals on model output values, when the data is unavailable. We introduce a class of algorithms in a Bayesian framework, relying on maximum entropy arguments and approximate Bayesian computation methods, to generate consistent data with the given summary statistics. Once we obtain consistent data sets, we pool the respective posteriors, to arrive at a single, averaged density on the parameters. This approach allows us to perform accurate forward uncertainty propagation consistent with the reported statistics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Accelerated basis adaptation in homogeneous chaos spaces

Polynomial chaos expansions (PCE) provide an efficient approach to uncertainty quantification (UQ) and have been adapted to diverse applications across the spectrum of science and engineering. For situations involving large stochastic parameterizations, the curse of dimensionality renders PCE-based methods computationally prohibitive. A basis adaptation approach for PCE was proposed by Tipireddy and Ghanem (2014) which transforms the input random variables through an isometry such that sufficient probabilistic characterization of specific quantities of interest (QoI) is concentrated in an algebraic manifold embedded in the linear span of the dominant transformed variables. While quite versatile, that original version of the PCE basis adaptation exhibited slow convergence for a number of problems of practical significance. In the present paper, we propose two novel methods to accelerate the convergence of the original basis adaptation approach, thus expanding its reach while also providing insight into its performance. In the first method, information gained from a pilot PCE representation is used to correct the mean and Gaussian coefficients in the adapted space. By taking advantages of probabilistic information in higher dimensional adaptation gleaned from an initial adaptation, the second method updates the rotation matrix used to identify the dominant transformed variables. In this manner, the new rotation matrix concentrates even more probabilistic information in its first few dimensions. These two method can be combined to achieve even better performance, the combined method is referred to as sequentially optimized adaptation method. Here, the methods are demonstrated on an analytical test function and a model of a space structure with several sub-components and a non-smooth quantity of interest representing the maximum acceleration over time. Both methods achieve accelerated convergence of the basis adaptation approach with negligible additional costs.

Accelerated convergence↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2011

Computational Design of Proteins Targeting the Conserved Stem Region of Influenza Hemagglutinin

We describe a general computational method for designing proteins that bind a surface patch of interest on a target macromolecule. Favorable interactions between disembodied amino acid residues and the target surface are identified and used to anchor de novo designed interfaces. The method was used to design proteins that bind a conserved surface patch on the stem of the influenza hemagglutinin (HA) from the 1918 H1N1 pandemic virus. After affinity maturation, two of the designed proteins, HB36 and HB80, bind H1 and H5 HAs with low nanomolar affinity. Further, HB80 inhibits the HA fusogenic conformational changes induced at low pH. The crystal structure of HB36 in complex with 1918/H1 HA revealed that the actual binding interface is nearly identical to that in the computational design model. Such designed binding proteins may be useful for both diagnostics and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Optimization Algorithms for Hierarchical Problems, with Application to Nanoporous Materials

This project investigated computational methods for optimizing the design of nanoporous materials. These materials hold promise in energy storage applications such as electrical supercapacitors and hydrogen storage systems. A significant issue for nanoporous storage devices is the design of the channel network within these materials to allow rapid diffusion while retaining storage capacity.

97 MATHEMATICS AND COMPUTING↗