Search NASA⌕ Search

SEARCH · Search NASA

Results for “computationally expensive function”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Probabilistic Forecast of Single–Phase Flow in Porous Media With Uncertain Properties

Uncertainty about geologic makeup and properties of the subsurface renders inadequate a unique quantitative prediction of flow and transport. Instead, multiple alternative scenarios have to be explored within the probabilistic framework, typically by means of Monte Carlo simulations (MCS). These can be computationally expensive, and often prohibitively so, especially when the goal is to compute the tails of a distribution, i.e., probabilities of rare events, which are necessary for risk assessment and decision making under uncertainty. We deploy the method of distributions to derive a deterministic equation for the cumulative distribution function (CDF) of hydraulic head in an aquifer with uncertain (random) hydraulic conductivity. The CDF equation relies on a self-consistent closure approximation, which ensures that the resulting CDF of hydraulic head has the same mean and variance as those computed with either statistical moment equation (the approach used herein) or MCS. We conduct a series of numerical experiments dealing with steady-state two-dimensional flow driven by either a natural hydraulic head gradient or a pumping well. Furthermore, these experiments reveal that the CDF method remains accurate and robust for highly heterogeneous formations with the variance of log conductivity as large as five. For the same accuracy, it is also up to four orders of magnitude faster than MCS in computing hydraulic head with a required degree of confidence (probability).

42 ENGINEERING↗
DOE OSTI2017

Modeling blur in various detector geometries for MeV radiography

Monte Carlo transport codes have been used to model the detector blur and energy deposition in various detector geometries for applications in MeV radiography. Segmented scintillating detectors, where low Z scintillators combined with a high-Z metal matrix, can be designed in which the resolution increases with increasing metal fraction. Here, the combination of various types of metal intensification screens and storage phosphor imaging plates has also been studied. A storage phosphor coated directly onto a metal intensification screen has superior performance over a commercial plate. Stacks of storage phosphor plates and tantalum intensification screens show an increase in energy deposited and detective quantum efficiency with increasing plate number, at the expense of resolution. Select detector geometries were tested by comparing simulation and experimental modulation transfer functions to validate the approach.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Creating Feature Representations of Antibody-Antigen Complexes for Fast Binding Prediction with Machine Learning

Antibody-antigen interactions underlie the function of the immune system. Numerous computational tools exist that aim to predict important quantities associated with antibodyantigen interactions, such as binding energy. Efficient and accurate prediction of such quantities would greatly accelerate vaccine (antigen) and therapeutic antibody design. Current methods and tools, however, rely on varying degrees of computationally expensive steps that limit their use on problems that require many predictions, and moreover, no single method produces consistently accurate estimates. In this project, we determined the feasibility of creating feature representations that concisely describe interacting pairs of antibody and antigen as rapidly computable numerical vectors, given the antibody-antigen co-complex. Such a featurization will allow the use of a statistical machine learning model that jointly models the results of the disparate existing computational tools available. In turn, this will greatly enhance their combined predictive performance. We developed a feature representation consisting of 86 descriptors of the antigen-antibody interface which is then reduced into a lower dimensional representation using a neural network dimensionality reduction step. We have demonstrated that this feature representation serves its intended purpose of rapidly describing antigen-antibody pairs for the task of binding strength prediction. While this feasibility study has proven that such a feature representation is feasible, additional research is likely improve its ability to support various tasks, including binding strength prediction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Quantum Mechanical Simulations of Complex Nanostructures for Photovoltaic Applications

A quantitative understanding of the electronic excitations in nanostructures, especially complex nanostructures, is crucial for making new-generation photovoltaic (PV) cells based on nanotechnology, which have high efficiency and low cost. Yet current quantum mechanical simulation methods are either computationally too expensive or not accurate and reliable enough, hindering the rational design of the nanoscale PV cells. The PI seeks to develop new methodologies to overcome the challenges in this very difficult and long-lasting problem, pushing the field forward so that electronic excitations can be accurately predicted for systems involving thousands of atoms. The primary objective of this project is to develop new approaches for electronic excitation calculations that are more accurate than traditional density functional theory (DFT) and are applicable to systems larger than what current beyond-DFT methods can treat. In this proposal, the PI will first address the excited-state problem within the DFT framework to obtain quasiparticle energies from both Kohn-Sham (KS) eigenvalues and orbitals; and the electron-hole binding energy will be computed based on screened Coulomb interaction of corresponding DFT orbitals. The accuracy of these approaches will be examined against many-body methods of GW/BSE and quantum Monte Carlo (QMC). The PI will also work on improving the accuracy and efficiency of the GW/BSE and QMC methods in electronic excitation computations by using better KS orbitals obtained from orbital-dependent DFT as inputs. Then an extended QMC database of ground- and excited-state properties will be developed, and this will be spot checked and supplemented with data from GW/BSE calculations. The investigation will subsequently focus on the development of an improved exchange-correlation (XC) density functional beyond the current generalized gradient approximation (GGA) level of parameterization, with parameters fitted to the QMC database. This will allow the ground-state properties of focus systems to be more precisely predicted using DFT. These new developments will then be applied to investigate a chosen set of complex nanostructures that have great potential for opening new routes in designing materials with improved transport, electronic, and optical properties for PV and other optoelectronic usages: (1) Hybrid interfaces between materials with distinct electronic and optical properties, such as organic molecules (conjugated polymers, e.g. P3HT) and inorganic semiconducting materials (Si and ZnO). Complicated interface structures, including interface bonding configurations, compositional and geometrical blending patterns, interfacial defects, and various sizes and shapes of inorganic nanomaterials, will be considered for the purpose of understanding the working mechanisms of present organic/nano PV systems and designing optimum interface structures for fast charge separation and injection. (2) Complex-structured semiconducting nanomaterials that could induce charge separation without pn- or hetero-junctions. The new methodology will allow the PI to investigate the performance of realistic semiconducting nanomaterials of internal (impurities, defects, etc.) and external (uneven surface, mechanical twisting and bending, surface chemistry, etc.) complexities on optical absorption and charge transport against charge trapping and recombination. Of particular interest is whether such structural complexity in a single material could even be beneficial for PV usage, for example, charge separation through morphology control. Successful completion of the proposed DFT methodology would have a far-reaching impact on our ability to study and exploit the nature of electronic excitations in complex materials, advancing the design of next-generation electronic and optoelectronic devices in all facets of renewable energy conversion and storage, including photovoltaics, thermoelectricity, photochemistry, etc.

14 SOLAR ENERGY↗
DOE OSTI2017

Atom-partitioned multipole expansions for electrostatic potential boundary conditions

Applications such as grid-based real-space density functional theory (DFT) use the Poisson equation to compute electrostatics. However, the expected long tail of the electrostatic potential requires either the use of a large and costly outer domain or Dirichlet boundary conditions estimated via multipole expansion. We find that the oft-used single-center spherical multipole expansion is only appropriate for isotropic mesh domains such as spheres and cubes. In this work, we introduce a method suitable for high aspect ratio meshes whereby the charge density is partitioned into atomic domains and multipoles are computed for each domain. While this approach is moderately more expensive than a single-center expansion, it is numerically stable and still a small fraction of the overall cost of a DFT calculation. The net result is that when high aspect ratio systems are being studied, form-fitted meshes can now be used in lieu of cubic meshes to gain computational speedup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2017

A near-stationary subspace for ridge approximation

Response surfaces are common surrogates for expensive computer simulations in engineering analysis. However, the cost of fitting an accurate response surface increases exponentially as the number of model inputs increases, which leaves response surface construction intractable for high-dimensional, nonlinear models. Here we describe ridge approximation for fitting response surfaces in several variables. A ridge function is constant along several directions in its domain, so fitting occurs on the coordinates of a low-dimensional subspace of the input space. We review essential theory for ridge approximation – e.g., the best mean-squared approximation and an optimal low-dimensional subspace – and we prove that the gradient-based active subspace is near-stationary for the least-squares problem that defines an optimal subspace. Motivated by the theory, we propose a computational heuristic that uses an estimated active subspace as an initial guess for a ridge approximation fitting problem. We show a simple example where the heuristic fails, which reveals a type of function for which the proposed approach is inappropriate. We then propose a simple alternating heuristic for fitting a ridge function, and we demonstrate the effectiveness of the active subspace initial guess applied to an airfoil model of drag as a function of its 18 shape parameters.

42 ENGINEERING↗
DOE OSTI2015

3D RISM theory with fast reciprocal-space electrostatics

The calculation of electrostatic solute-solvent interactions in 3D RISM (“three-dimensional reference interaction site model”) integral equation theory is recast in a form that allows for a computational treatment analogous to the “particle-mesh Ewald” formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Design Paradigm for Modular Multilevel Converter-Based Generator Rectifier Systems

Modular Multilevel Converters (MMC) are being widely considered for medium to high voltage DC applications. Designing such converters through multi-objective optimization is of interest because such an approach allows the trade-off between competing objectives (for example mass and loss) to be explicitly and quantitatively identified. In this work, an optimization based design paradigm for MMC based generator rectifier systems is proposed. Such development typically requires detailed component design and simulation models for the electric machine and converter which are computationally expensive. As an alternative, the proposed work utilizes an electric machine metamodel, inductor metamodel, and high-speed steady-state simulation model for the MMC to facilitate multi-objective optimization minimizing system metrics of interest while satisfying system constraints. A case study is undertaken to demonstrate the functionality of the proposed design paradigm.

42 ENGINEERING↗
DOE OSTI2014

Assessing the standard Molybdenum projector augmented wave VASP potentials

Density Functional Theory (DFT) based Equation of State (EOS) construction is a prominent part of Sandia’s capabilities to support engineering sciences. This capability is based on augmenting experimental data with information gained from computational investigations, especially in those parts of the phase space where experimental data is hard, dangerous, or expensive to obtain. A key part of the success of the Sandia approach is the fundamental science work supporting the computational capability. Not only does this work enhance the capability to perform highly accurate calculations but it also provides crucial insight into the limitations of the computational tools, providing high confidence in the results even where results cannot be, or have not yet been, validated by experimental data. This report concerns the key ingredient of projector augmented-wave (PAW) potentials for use in pseudo-potential computational codes. Using the tools discussed in SAND2012-7389 we assess the standard Vienna Ab-initio Simulation Package (VASP) PAWs for Molybdenum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Determination of Optimal Operating Schemes for a Multi-Reservoir System Under Environmental Constraints

In strati ed reservoirs, both dam tailwater discharge and thermal plant intake water quality and temperature can be highly dependent on structure depth. A two- dimensional laterally-averaged model allows for better prediction of water quality over time at specfi c depths. Because high- fidelity models are typically too computationally expensive for direct inclusion within optimization algorithms, water quality is incorporated using one dimensional models are simple flow requirements. Water quality predictions can be incorporated within the optimization process through using surrogate modeling methods, in this application artificial neural network (ANN) models. ANNs are flexible machine learning tools for function approximation composed of a structure of neurons assembled within a multi-layer architecture. They are capable of handling large amounts of training data and modeling nonlinear dynamic systems, making ANNs a well-suited method for this application. This report illustrates the development of ANN models to emulate the hydrodynamic and water quality modeling capabilities of the high- fidelity, two-dimensional CE-QUAL-W2 (W2) model, as well as a linked riverine reservoir system optimization process which accounts for energy generation, water balance and hydraulics, and compliance point water quality. A process for hourly hydropower generation planning is demonstrated on a pair of reservoirs linked in series. The two reservoirs are U.S. Army Corps of Engineers projects with hydropower capabilities on the Cumberland River near Nashville, Tennessee, USA.

13 HYDRO ENERGY↗
DOE OSTI2017

Why wait? Let's start computing while data is still on the wire

In this era of Big Data, computing useful information from data is becoming increasingly complicated, particularly due to the ever increasing volumes of data that need to travel over the network to data centers to be stored and processed, all highly expensive operations in the long haul. In this paper we suggest that we can do computing and analysis of data "on the wire," i.e., while data is still in transit. The nature of these computations include analysis, visualization, pattern recognition, and prediction, or forecasting, on the streaming data. We follow a service function chaining architecture to implement this, assuming that the data packets arrive within a single network administrative domain. As a demonstration of this new computing paradigm, we present three examples. Firstly, we demonstrate pattern recognition and data visualization on streaming forex data, which can be used for lucrative trading in the forex market. In our second example, we analyze and learn user buying patterns from clickstream data streaming from multiple websites. Finally, we monitor solar sensors for a zero reading while the packets are still on their way to the data center, to schedule any maintenance and requisite repairs with no time delay.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ML-based Dimension Reduction Strategies

Deep learning (DL)--based surrogate models have achieved success in various applications in carbon capture and storage (CCS). However, the model training on high-dimensional spaces is computationally expensive and impractical for large-scale and complex geological models, because the models usually contain hundreds of thousands to millions of grid cells, each with a set of parameters. Furthermore, the high cost of generating training data with sufficient variation is another limitation of model training on high-dimensional spaces, which may result in overfitting and reduce the model efficiency and prediction performance. We proposed the workflow incorporating dimension reduction methods and deep learning models, which aim to extract the latent variables of input parameters and output state variables, and then build the mapping function at the latent spaces. The proposed workflow can significantly reduce the computational complexity in solving both forward and inverse problems compared to models trained on high-dimensional spaces. Dimensionality reduction models showed great potential in workflows for fast reservoir simulation, history matching, prior model generation, visualization, and more, ultimately enhancing DL model performance in related SMART Work Packages.

Hosseini, Seyyed↗
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Furthermore, our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

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

Improved reference system for the corrected rigid spheres equation of state model

The Corrected Rigid Spheres (CRIS) EOS model, developed from fluid perturbation theory using a hard sphere reference system, has been successfully used to calculate the EOS of many materials, including gases and metals. The radial distribution function (RDF) plays a pivotal role in choosing the sphere diameter, through a variational principle, as well as the thermodynamic response. Despite its success, the CRIS model has some shortcomings in that it predicts too large a temperature for liquid-vapor critical points, can break down at large compression, and is computationally expensive. In this paper, we first demonstrate that an improved analytic representation of the hard sphere RDF does not alleviate these issues. Relaxing the strict adherence of the RDF to hard spheres allows an accurate fit to the isotherms and vapor dome of the Lennard-Jones fluid using an arbitrary reference system. Second order correction are eliminated, limiting the breakdown at large compression and significantly reducing the computation cost. The transferrability of the new model to real systems is demonstrated on argon, with an improved vapor dome compared to the original CRIS model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

TGCM: (T)rait, (G)ene, and (C)rop Growth (M)odel Directed Targeted Gene Characterization in Sorghum (Final Technical Report)

Understanding which genes control important crop traits could help scientists develop better bioenergy and food crops more efficiently. However, plant genomes contain tens of thousands of genes, and testing each one individually is expensive and time-consuming. This project developed computational tools to predict which genes are most likely to matter, allowing researchers to focus their efforts where they will have the greatest impact. This project developed and validated integrated approaches combining machine learning, quantitative genetics, and crop growth modeling to improve the efficiency of functional gene characterization in sorghum (Sorghum bicolor), a critical bioenergy and food security crop. The research addressed a fundamental challenge in plant biology: the majority of genes in plant genomes lack experimentally validated functions, making it difficult to prioritize which genes to study using resource-intensive reverse genetics approaches.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Algorithm 1028: VTMOP: Solver for Blackbox Multiobjective Optimization Problems

VTMOP is a Fortran 2008 software package containing two Fortran modules for solving computationally expensive bound-constrained blackbox multiobjective optimization problems. VTMOP implements the algorithm of [32], which handles two or more objectives, does not require any derivatives, and produces well-distributed points over the Pareto front. The first module contains a general framework for solving multiobjective optimization problems by combining response surface methodology, trust region methodology, and an adaptive weighting scheme. The second module features a driver subroutine that implements this framework when the objective functions can be wrapped as a Fortran subroutine. Lastly, support is provided for both serial and parallel execution paradigms, and VTMOP is demonstrated on several test problems as well as one real-world problem in the area of particle accelerator optimization.

97 MATHEMATICS AND COMPUTING↗