SEARCH · Search NASA

Results for “Approximation 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Techno-Economic Comparison of Molten-Salt Electrolysis and Carbothermic Reduction for the Production of Metallurgical-Grade Silicon

Metallurgical-grade silicon (MG-Si) is an important source material for many industrial applications, including the manufacture of alloys, solar photovoltaics, and electronics. The process to refine raw materials into MG-Si is energy-intensive, with the predominant method of submerged-arc furnaces requiring energy consumption of approximately 11–13 kWh/kg Si. Recent research has discussed promising methods for reducing the energy required for the silicon production process, including the use of molten-salt electrolysis (MSE), a technique that offers potential savings in energy consumption without requiring carbon inputs for the process. This paper presents a techno-economic study of a potential industrial-scale MSE plant for MG-Si production to evaluate the trade-offs between capital and operating costs of the system. Capital costs are sourced from recent MG-Si plants and an existing cost model developed for MSE processes that includes the size of the plant and the operating temperature among its inputs. The results show that MSE technology has the potential to be an economically cost-competitive option for MG-Si production if the technology successfully scales to industrial production and matures enough to allow for financing costs similar to that of a comparably sized submerged-arc furnace plant.

14 SOLAR ENERGY↗

DOE OSTI2023

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Constraining DALECv2 using multiple data streams and ecological constraints: analysis and application

We use a variational method to assimilate multiple data streams into the terrestrial ecosystem carbon cycle model DALECv2 (Data Assimilation Linked Ecosystem Carbon). Ecological and dynamical constraints have recently been introduced to constrain unresolved components of this otherwise ill-posed problem. We recast these constraints as a multivariate Gaussian distribution to incorporate them into the variational framework and we demonstrate their advantage through a linear analysis. By using an adjoint method we study a linear approximation of the inverse problem: firstly we perform a sensitivity analysis of the different outputs under consideration, and secondly we use the concept of resolution matrices to diagnose the nature of the ill-posedness and evaluate regularisation strategies. We then study the non-linear problem with an application to real data. Finally, we propose a modification to the model: introducing a spin-up period provides us with a built-in formulation of some ecological constraints which facilitates the variational approach.

58 GEOSCIENCES↗

DOE OSTI2019

Summary: Verification of the ENDF/B-VII.1 and VIII.0 AMPX 1597-group Libraries for Advanced Reactor Analysis

The SCALE/AMPX multigroup (MG) cross section processing procedure has been updated to minimize reactivity differences for various advanced thermal and fast reactor designs, as observed for the current MG libraries, resulting in excellent agreement between the calculations with the new MG libraries and the continuous-energy reference calculations. The SCALE MG calculations are widely applied to thermal spectrum light-water reactor systems, as well as fast spectrum metallic systems. Due to growing interest from industry and regulators in applying SCALE for the design of fast spectrum reactors—both sodium and molten salt—it was desirable to review and improve the SCALE/AMPX procedure for unresolved resonance self-shielded data and high-energy neutron spectra. The data were improved by generating MG unresolved resonance data based on the analytic probability table method with the narrow resonance approximation and by using very fine and intermediate group structures that are typical for fast system analysis. This study focused on verifying the improved SCALE/AMPX MG cross section processing procedure and the new AMPX 1597-group library with the ENDF/B-VII.1 and VIII.0 evaluated nuclear data. The verification was made by performing reaction rate analysis and benchmark calculations for various thermal and fast reactor systems. Results indicate that the improved SCALE/AMPX MG cross section processing and libraries provide excellent results for advanced reactor analysis.

Kim, Kang Seog↗

DOE OSTI2019

Verification of the ENDF/B-VII.1 and VIII.0 AMPX 1597-Group Libraries for Advanced Reactor Analysis

Kim, Kang Seog↗

DOE OSTI2018

Hydraulic Properties of Grout Monoliths Containing Organic Ion Exchange Resin - 18369

The Waste Treatment and Immobilization Plant (WTP) at Hanford is being constructed to treat 56 million gallons of radioactive waste currently stored in underground tanks at the Hanford site. Operation of the WTP will generate several solid secondary waste (SSW) streams including spent spherical resorcinol formaldehyde cation exchange resin (IXr). The current plan for safe disposal of IXr is solidification in cementitious material (CM) typically comprised of cement, fly ash and blast furnace slag. Three CM formulations were tested to stabilize IXr for long-term land disposal. Two IXr loading rates, 0.1 and 0.3 v/v IXr:grout were used with the three CM formulations to prepare simulated waste form monoliths for hydraulic property testing. Following curing, the monoliths were removed from molds and tested to determine saturated hydraulic conductivity (K{sub sat}) and water retention properties. K{sub sat} was measured using a standard flexible permeameter method ASTM D-5084-10. A constant vapor pressure (CVP) method was used to determine van Genuchten parameters for desorption moisture retention curves for each of the SSW grout formulations. For the CVP method, vacuum saturated, intact wafers approximately 1.25 cm thick were placed above a saturated salt solution inside a sealed container for equilibration. Each saturated salt solution has an associated vapor pressure related to relative humidity in the container. Visual examination and three-dimensional x-ray micro computed tomography (μCT) indicated that the IXr was uniformly distributed in the cured waste form. Saturated hydraulic conductivity was within the range previously reported by Flach et al. in a study of the hydraulic properties of 98 pastes and grouts. The solidification of IXr in CM resulted in a waste form with water retention properties similar to cementitious pastes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2016

An improved multilevel Monte Carlo method for estimating probability distribution functions in stochastic oil reservoir simulations

In this paper, we develop an improved multilevel Monte Carlo (MLMC) method for estimating cumulative distribution functions (CDFs) of a quantity of interest, coming from numerical approximation of large-scale stochastic subsurface simulations. Compared with Monte Carlo (MC) methods, that require a significantly large number of high-fidelity model executions to achieve a prescribed accuracy when computing statistical expectations, MLMC methods were originally proposed to significantly reduce the computational cost with the use of multifidelity approximations. The improved performance of the MLMC methods depends strongly on the decay of the variance of the integrand as the level increases. However, the main challenge in estimating CDFs is that the integrand is a discontinuous indicator function whose variance decays slowly. To address this difficult task, we approximate the integrand using a smoothing function that accelerates the decay of the variance. In addition, we design a novel a posteriori optimization strategy to calibrate the smoothing function, so as to balance the computational gain and the approximation error. The combined proposed techniques are integrated into a very general and practical algorithm that can be applied to a wide range of subsurface problems for high-dimensional uncertainty quantification, such as a fine-grid oil reservoir model considered in this effort. The numerical results reveal that with the use of the calibrated smoothing function, the improved MLMC technique significantly reduces the computational complexity compared to the standard MC approach. Finally, we discuss several factors that affect the performance of the MLMC method and provide guidance for effective and efficient usage in practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Green’s function coupled cluster formulations utilizing extended inner excitations

Here, we analyze new approximations of the Green’s function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i( n, m ) method, where the excitation level of the inner auxiliary operators ( m ) used to describe the ionization potential and electron affinity effects in the N-1 and N+1 particle spaces is higher than the excitation level ( n ) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called “n + 1” rule in this category [or the GFCC-i( n, n + 1) method], which states that in order to maintain size-extensivity of the Green’s function matrix elements, the excitation level of inner auxiliary operators X p (ω) and Y q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N 2 and CO molecules in the inner and outer valence regimes. Lastly, in comparison with the Green’s function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Explicitly correlated ring-coupled-cluster-doubles theory

The connection between the random-phase approximation and the ring-coupled-cluster-doubles method bridges the gap between density-functional and wave-function theories and the importance of the random-phase approximation lies in both its broad applicability and this linking role in electronic-structure theory. In this contribution, we present an explicitly correlated approach to the random-phase approximation, based on the direct ring-coupled-cluster-doubles ansatz, which overcomes the problem of slow basis-set convergence, inherent to the random-phase approximation. Benchmark results for a test set of 106 molecules and a selection of 10 organic complexes from the S22 test set demonstrate that convergence to within 99% of the basis-set limit is reached for triple-zeta basis sets for atomisation energies, while quadruple-zeta basis sets are required for interaction energies. Corrections due to single excitations into the complementary auxiliary space reduce the basis-set incompleteness error by one order of magnitude, while contributions due to the coupling of conventional and geminal amplitudes are in general negligible. We find that a non-iterative explicitly correlated correction to first order in perturbation theory exhibits the best ratio of accuracy to computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Information-theoretic discrepancy based iterative reconstructions (IDIR) for polychromatic x-ray tomography

Purpose: X-ray photons generated from a typical x-ray source for clinical applications exhibit a broad range of wavelengths, and the interactions between individual particles and biological substances depend on particles' energy levels. Most existing reconstruction methods for transmission tomography, however, neglect this polychromatic nature of measurements and rely on the monochromatic approximation. In this study, we developed a new family of iterative methods that incorporates the exact polychromatic model into tomographic image recovery, which improves the accuracy and quality of reconstruction.Methods: The generalized information-theoretic discrepancy (GID) was employed as a new metric for quantifying the distance between the measured and synthetic data. By using special features of the GID, the objective function for polychromatic reconstruction which contains a double integral over the wavelength and the trajectory of incident x-rays was simplified to a paraboloidal form without using the monochromatic approximation. More specifically, the original GID was replaced with a surrogate function with two auxiliary, energy-dependent variables. Subsequently, the alternating minimization technique was applied to solve the double minimization problem. Based on the optimization transfer principle, the objective function was further simplified to the paraboloidal equation, which leads to a closed-form update formula. Numerical experiments on the beam-hardening correction and material-selective reconstruction were conducted to compare and assess the performance of conventional methods and the proposed algorithms.Results: The authors found that the GID determines the distance between its two arguments in a flexible manner. In this study, three groups of GIDs with distinct data representations were considered. The authors demonstrated that one type of GIDs that comprises “raw” data can be viewed as an extension of existing statistical reconstructions; under a particular condition, the GID is equivalent to the Poisson log-likelihood function. The newly proposed GIDs of the other two categories consist of log-transformed measurements, which have the advantage of imposing linearized penalties over multiple discrepancies. For all proposed variants of the GID, the aforementioned strategy was used to obtain a closed-form update equation. Even though it is based on the exact polychromatic model, the derived algorithm bears a structural resemblance to conventional methods based on the monochromatic approximation. The authors named the proposed approach as information-theoretic discrepancy based iterative reconstructions (IDIR). In numerical experiments, IDIR with raw data converged faster than previously known statistical reconstruction methods. IDIR with log-transformed data exhibited superior reconstruction quality and faster convergence speed compared with conventional methods and their variants.Conclusions: The authors' new framework for tomographic reconstruction allows iterative inversion of the polychromatic data model. The primary departure from the traditional iterative reconstruction was the employment of the GID as a new metric for quantifying the inconsistency between the measured and synthetic data. The proposed methods outperformed not only conventional methods based on the monochromatic approximation but also those based on the polychromatic model. The authors have observed that the GID is a very flexible means to design an objective function for iterative reconstructions. Hence, the authors expect that the proposed IDIR framework will also be applicable to other challenging tasks.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2021

Efficient calculation of carrier scattering rates from first principles

Abstract The electronic transport behaviour of materials determines their suitability for technological applications. We develop a computationally efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 23 semiconductors, including the large 48 atom CH 3 NH 3 PbI 3 hybrid perovskite, and comparing the results against experimental measurements and more detailed scattering simulations. The Spearman rank coefficient of mobility against experiment ( r s = 0.93) improves significantly on results obtained using a constant relaxation time approximation ( r s = 0.52). We find our approach offers similar accuracy to state-of-the art methods at approximately 1/500th the computational cost, thus enabling its use in high-throughput computational workflows for the accurate screening of carrier mobilities, lifetimes, and thermoelectric power.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches

Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and their induced asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. In this work, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic timescale separation; the Franck-Condon method, which includes dynamical effects but only at the harmonic level; as well as the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Multi-element least square HDMR methods and their applications for stochastic multiscale model reduction

Stochastic multiscale modeling has become a necessary approach to quantify uncertainty and characterize multiscale phenomena for many practical problems such as flows in stochastic porous media. The numerical treatment of the stochastic multiscale models can be very challengeable as the existence of complex uncertainty and multiple physical scales in the models. To efficiently take care of the difficulty, we construct a computational reduced model. To this end, we propose a multi-element least square high-dimensional model representation (HDMR) method, through which the random domain is adaptively decomposed into a few subdomains, and a local least square HDMR is constructed in each subdomain. These local HDMRs are represented by a finite number of orthogonal basis functions defined in low-dimensional random spaces. The coefficients in the local HDMRs are determined using least square methods. We paste all the local HDMR approximations together to form a global HDMR approximation. To further reduce computational cost, we present a multi-element reduced least-square HDMR, which improves both efficiency and approximation accuracy in certain conditions. To effectively treat heterogeneity properties and multiscale features in the models, we integrate multiscale finite element methods with multi-element least-square HDMR for stochastic multiscale model reduction. This approach significantly reduces the original model's complexity in both the resolution of the physical space and the high-dimensional stochastic space. We analyze the proposed approach, and provide a set of numerical experiments to demonstrate the performance of the presented model reduction techniques. - Highlights: • Multi-element least square HDMR is proposed to treat stochastic models. • Random domain is adaptively decomposed into some subdomains to obtain adaptive multi-element HDMR. • Least-square reduced HDMR is proposed to enhance computation efficiency and approximation accuracy in certain conditions. • Integrating MsFEM and multi-element least square HDMR can significantly reduce computation complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Computational Homogenization of Polycrystals

This paper reviews the current state of the art in the simulation of the mechanical behavior of polycrystalline materials by means of computational homogenization. The key ingredients of this modeling strategy are presented in detail starting with the parameters needed to describe polycrystalline microstructures and the digital representation of such microstructures in a suitable format to perform computational homogenization. The different crystal plasticity frameworks that can describe the physical mechanisms of deformation in single crystals (dislocation slip and twinning) at the microscopic level are presented next. This is followed by the description of computational homogenization methods based on mean-field approximations by means of the viscoplastic self-consistent approach, or on the full-field simulation of the mechanical response of a representative polycrystalline volume element by means of the finite element method or the fast Fourier transform-based method. Multiscale frameworks based on the combination of mean-field homogenization and the finite element method are presented next to model the plastic deformation of polycrystalline specimens of arbitrary geometry under complex mechanical loading. Examples of application to predict the strength, fatigue life, damage, and texture evolution under different conditions are presented to illustrate the capabilities of the different models. Lastly, current challenges and future research directions in this field are summarized.

36 MATERIALS SCIENCE↗

DOE OSTI2015

The Hubble flow of plateau inflation

In the absence of CMB precision measurements, a Taylor expansion has often been invoked to parametrize the Hubble flow function during inflation. The standard ''horizon flow'' procedure implicitly relies on this assumption. However, the recent Planck results indicate a strong preference for plateau inflation, which suggests the use of Padé approximants instead. We propose a novel method that provides analytic solutions of the flow equations for a given parametrization of the Hubble function. This method is illustrated in the Taylor and Padé cases, for low order expansions. We then present the results of a full numerical treatment scanning larger order expansions, and compare these parametrizations in terms of convergence, prior dependence, predictivity and compatibility with the data. Finally, we highlight the implications for potential reconstruction methods.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Analysis of the static magnetic field-dependent optical transmission of Ni nanorod colloidal suspensions

The magnetic field-dependent optical transmission of dilute Ni nanorod aqueous suspensions was investigated. A series of four samples of nanorods were synthesized using the AAO template method and processed to stable colloids. The distributions of their length and diameter were characterized by analysis of TEM images and revealed average diameters of ∼25 nm and different lengths in the range of 60 nm–1100 nm. The collinear magnetic and optical anisotropy was studied by static field-dependent transmission measurements of linearly polarized light parallel and perpendicular to the magnetic field direction. The experimental results were modelled assuming the field-dependent orientation distribution function of a superparamagnetic ensemble for the uniaxial ferromagnetic nanorods in liquid dispersion and extinction cross sections for longitudinal and transversal optical polarization derived from different approaches, including the electrostatic approximation and the separation of variables method, both applied to spheroidal particles, as well as finite element method simulations of spheroids and capped cylindrical particles. The extinction cross sections were compared to reveal the differences associated with the approximations of homogeneous polarization and/or particle shape. The consequences of these approximations for the quantitative analysis of magnetic field-dependent optical transmission measurements were investigated and a reliable protocol derived. Furthermore, the changes in optical cross sections induced by electromagnetic interaction between two nanorods in parallel end-to-end and side-by-side configuration as a function of their separation were studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗