Search NASA⌕ Search

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗
DOE OSTI2015

Determination of recombination radius in Si for binary collision approximation codes

Displacement damage caused by ions or neutrons in microelectronic devices can have significant effect on the performance of these devices. Therefore, it is important to predict not only the displacement damage profile, but also its magnitude precisely. Analytical methods and binary collision approximation codes working with amorphous targets use the concept of displacement energy, the energy that a lattice atom has to receive to create a permanent replacement. It was found that this “displacement energy” is direction dependent; it can range from 12 to 32 eV in silicon. Obviously, this model fails in BCA codes that work with crystalline targets, such as Marlowe. Marlowe does not use displacement energy; instead, it uses lattice binding energy only and then pairs the interstitial atoms with vacancies. Then based on the configuration of the Frenkel pairs it classifies them as close, near, or distant pairs, and considers the distant pairs the permanent replacements. Unfortunately, this separation is an ad hoc assumption, and the results do not agree with molecular dynamics calculations. After irradiation, there is a prompt recombination of interstitials and vacancies if they are nearby, within a recombination radius. In order to implement this recombination radius in Marlowe, we used the comparison of MD and Marlowe calculation in a range of ion energies in single crystal silicon target. As a result, the calculations showed that a single recombination radius of ~7.4 Å in Marlowe for a range of ion energies gives an excellent agreement with MD.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Local Hamiltonians for quantitative Green's function embedding methods

Embedding calculations that find approximate solutions to the Schrödinger equation for large molecules and realistic solids are performed commonly in a three step procedure involving (i) construction of a model system with effective interactions approximating the low energy physics of the initial realistic system, (ii) mapping the model system onto an impurity Hamiltonian, and (iii) solving the impurity problem. We have developed a novel procedure for parametrizing the impurity Hamiltonian that avoids the mathematically uncontrolled step of constructing the low energy model system. Instead, the impurity Hamiltonian is immediately parametrized to recover the self-energy of the realistic system in the limit of high frequencies or short time. The effective interactions parametrizing the fictitious impurity Hamiltonian are local to the embedded regions, and include all the non-local interactions present in the original realistic Hamiltonian in an implicit way. We show that this impurity Hamiltonian can lead to excellent total energies and self-energies that approximate the quantities of the initial realistic system very well. Moreover, we show that as long as the effective impurity Hamiltonian parametrization is designed to recover the self-energy of the initial realistic system for high frequencies, we can expect a good total energy and self-energy. Lastly, we propose two practical ways of evaluating effective integrals for parametrizing impurity models.

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

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗
DOE OSTI2018

Exact and inexact subsampled Newton methods for optimization

Abstract The paper studies the solution of stochastic optimization problems in which approximations to the gradient and Hessian are obtained through subsampling. We first consider Newton-like methods that employ these approximations and discuss how to coordinate the accuracy in the gradient and Hessian to yield a superlinear rate of convergence in expectation. The second part of the paper analyzes an inexact Newton method that solves linear systems approximately using the conjugate gradient (CG) method, and that samples the Hessian and not the gradient (the gradient is assumed to be exact). We provide a complexity analysis for this method based on the properties of the CG iteration and the quality of the Hessian approximation, and compare it with a method that employs a stochastic gradient iteration instead of the CG method. We report preliminary numerical results that illustrate the performance of inexact subsampled Newton methods on machine learning applications based on logistic regression.

Mathematics↗
DOE OSTI2015

Range-separated double-hybrid density-functional theory applied to periodic systems

Quantum chemistry methods exploiting density-functional approximations for short-range electron-electron interactions and second-order Møller-Plesset (MP2) perturbation theory for long-range electron-electron interactions have been implemented for periodic systems using Gaussian-type basis functions and the local correlation framework. The performance of these range-separated double hybrids has been benchmarked on a significant set of systems including rare-gas, molecular, ionic, and covalent crystals. The use of spin-component-scaled MP2 for the long-range part has been tested as well. The results show that the value of μ = 0.5 bohr{sup −1} for the range-separation parameter usually used for molecular systems is also a reasonable choice for solids. Overall, these range-separated double hybrids provide a good accuracy for binding energies using basis sets of moderate sizes such as cc-pVDZ and aug-cc-pVDZ.

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

The diamagnetic susceptibility of a donor in a semiconductor core shell quantum dot

The effect of Aluminium concentration, shell thickness and size of the core shell Quantum Dot on the Diamagnetic Susceptibility of a donor in the Core Shell Quantum Dot is calculated in the effective mass approximation using the variational method. The results are presented and discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2021

Approximation of the Frequency-Amplitude Curve Using the Homotopy Analysis Method

The Frequency-Amplitude (F-A) curve on power system oscillation under a large disturbance characterizes how a natural oscillation mode transitions to nonlinear oscillations with growing amplitudes and decaying frequencies. The existing formulation of the F-A curve is derived by solving elliptical integrals on oscillation of a single-machine-infinite-bus equivalent about the targeted oscillation mode. The formula is in a form of infinite series and needs to sum a large number of terms for satisfactory accuracy. This paper introduces an explicit, approximate expression obtained from the Homotopy Analysis Method on the F-A curve. The proposed F-A curve expression is derived from an SMIB system and verified on the IEEE 3-machine 9-bus system to show how the oscillation frequency of a dominant mode varies with oscillation amplitude under large disturbances.

frequency-amplitude (F-A) curve↗
DOE OSTI2017

Annotations on the virtual element method for second-order elliptic problems

This document contains working annotations on the Virtual Element Method (VEM) for the approximate solution of diffusion problems with variable coefficients. To read this document you are assumed to have familiarity with concepts from the numerical discretization of Partial Differential Equations (PDEs) and, in particular, the Finite Element Method (FEM). This document is not an introduction to the FEM, for which many textbooks (also free on the internet) are available. Eventually, this document is intended to evolve into a tutorial introduction to the VEM (but this is really a long-term goal).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Magnetic Properties of La-Doped ZnO(0001)–Zn Polar Surface with and without Vacancies: a First-Principle Study

The effects of La doping with and without vacancies on the magnetism of a ZnO(0001)–Zn polar surface were studied using the generalized gradient approximation + U (GGA + U) method of density functional theory to determine the magnetism source of La-doped ZnO. Results revealed that La doping with O vacancy did not cause ZnO ferromagnetism, but La-doped ZnO with Zn vacancy was ferromagnetic. The net magnetic moment was primarily related to the unpaired electron amount of the O 2p state and slightly modified by orbital and anomalous magnetic moments. The anomalous magnetic moment originated from the radiating magnetism moment via quantum electrodynamics and was confirmed to exist via theoretical calculations of magnetism. Formation energy showed that La doping easily generated Zn vacancies in ZnO, but magnetic quenching occurred when the ratio of La ions to Zn vacancies was 2:1. The unpaired O 2p electrons induced by Zn vacancies were the magnetism source of La-doped ZnO. Moreover, the Curie temperature of La-doped ZnO with Zn vacancy reached room temperature. Results showed that the ZnO(0001)–Zn polar surface can be utilized as the surface of ZnO-based diluted magnetic semiconductors because of its unipolarity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Optimal, centralized dynamic curbside parking space zoning

In this paper we formulate a dynamic mixed integer program for optimally zoning curbside parking spaces subject to transportation policy-inspired constraints and regularization terms. First, we illustrate how given some objective of curb zoning valuation as a function of zone type (paid parking, bus stop, etc.), dynamically rezoning involves unrolling this optimization program over a fixed time horizon. Second, we implement two different solution methods given an example curb zoning valuation. In the first method, we solve long horizon dynamic zoning problems via approximate dynamic programming. In the second method, we employ Dantzig-Wolfe decomposition to break-up the mixed-integer program into a master problem and several sub-problems that can be solved in parallel. This speeds up the computational solve-time of the MIP considerably. We present simulation results and comparisons of the different employed techniques on vehicle arrival-rate data obtained for a neighborhood in downtown Seattle, Washington, USA.

Nazir, Mohammad Nawaf↗
DOE OSTI2015

TU-F-CAMPUS-I-01: Investigation of the Effective Dose From Bolus Tracking Acquisitions at Different Anatomical Locations in the Chest for CT

Purpose: Stationary table acquisitions (Bolus tracking) in X-ray Computed Tomography (CT) can Result in dose length products (DLP) comparable to spiral scans. It is today unclear whether or not the effective dose (E) for Bolus Tracking can be approximated using target region specific conversion factors (E/DLP). The purpose of this study was to investigate how E depends on the anatomical location of the Bolus Tracking in relation to Chest CT scans with the same DLP. Methods: Effective doses were approximated for the ICRP 110 adult Reference Male (AM) and adult Reference Female (FM) computational voxel phantoms using software for CT dose approximations (pre-simulated MC data). The effective dose was first approximated for a Chest CT scan using spiral technique and a CTDIvol (32 cm) of 6 mGy. The effective dose from the spiral scan was then compared to E approximated for contiguous Bolus Tracking acquisitions (1 cm separation), with a total collimation of 1 cm, over different locations of the chest of the voxel phantoms. The number of rotations used for the Bolus Tracking acquisitions was adjusted to yield the same DLP (32 cm) as the spiral scan. Results: Depending on the anatomical location of the Bolus Tracking, E ranged by factors of 1.3 to 6.8 for the AM phantom and 1.4 to 3.3 for the AF phantom, compared to the effective dose of the spiral scans. The greatest E for the Bolus Tracking acquisitions was observed for anatomical locations coinciding with breast tissue. This can be expected as breast tissue has a high tissue weighting factor in the calculation of E. Conclusion: For Chest CT scans, the effective dose from Bolus Tracking is highly dependent on the anatomical location where the scan is administered and will not always accurately be represented using target region specific conversion factors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-IeP3-04: Effective Dose Measurements in Fast Kvp Switch Dual Energy Computed Tomography

Purpose: The objective of this study was two-fold: (a) to test a new approach to approximating organ dose by using the effective energy of the combined 80kV/140kV beam in dual-energy (DE) computed tomography (CT), and (b) to derive the effective dose (ED) in the abdomen-pelvis protocol in DECT. Methods: A commercial dual energy CT scanner was employed using a fast-kV switch abdomen/pelvis protocol alternating between 80 kV and 140 kV. MOSFET detectors were used for organ dose measurements. First, an experimental validation of the dose equivalency between MOSFET and ion chamber (as a gold standard) was performed using a CTDI phantom. Second, the ED of DECT scans was measured using MOSFET detectors and an anthropomorphic phantom. For ED calculations, an abdomen/pelvis scan was used using ICRP 103 tissue weighting factors; ED was also computed using the AAPM Dose Length Product (DLP) method and compared to the MOSFET value. Results: The effective energy was determined as 42.9 kV under the combined beam from half-value layer (HVL) measurement. ED for the dual-energy scan was calculated as 16.49 ± 0.04 mSv by the MOSFET method and 14.62 mSv by the DLP method. Conclusion: Tissue dose in the center of the CTDI body phantom was 1.71 ± 0.01 cGy (ion chamber) and 1.71 ± 0.06 (MOSFET) respectively; this validated the use of effective energy method for organ dose estimation. ED from the abdomen-pelvis scan was calculated as 16.49 ± 0.04 mSv by MOSFET and 14.62 mSv by the DLP method; this suggests that the DLP method provides a reasonable approximation to the ED.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Structural transformation features in titanium crystallite under mechanical loading

The features of defect structure generation and development in titanium crystallites were studied on the basis of molecular dynamics method. Interatomic interaction was described using many-body potentials calculated in the approximation of the embedded atom method. It is shown that local structural changes begin to occur in the crystallite at the achievement of a threshold strain value, which is accompanied by a dramatic decrease of potential energy. The features of the formation of local structural changes, which are precursors of classical stacking faults, were studied.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗
DOE OSTI2018

Generalized Jacobi–Galerkin method for nonlinear fractional differential algebraic equations

In this paper, we provide an approximate approach based on the Galerkin method to solve a class of nonlinear fractional differential algebraic equations. The fractional derivative operator in the Caputo sense is utilized and the generalized Jacobi functions are employed as trial functions. The existence and uniqueness theorem as well as the asymptotic behavior of the exact solution are provided. It is shown that some derivatives of the solutions typically have singularity at origin dependence on the order of the fractional derivative. The influence of the perturbed data on the exact solutions along with the convergence analysis of the proposed scheme is also established. Some illustrative examples provided to demonstrate that this novel scheme is computationally efficient and accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Structural, electronic and elastic properties of REIr{sub 2} (RE=La and Ce) Laves phase compounds

REIr{sub 2} (RE = La and Ce) Laves phase intermetallic compounds were investigated with respect to their structural, electronic and elastic properties using full potential linearized augmented plane wave (FP-LAPW) method within generalized gradient approximation (GGA) as implemented in WIEN2k code. The ground state properties such as lattice constants (a{sub 0}), bulk modulus (B), pressure derivative of bulk modulus (B′) and density of state at Fermi level N(E{sub F}) have been obtained by optimization method. The electronic structure (BS, TDOS and PDOS) reveals that these Laves phase compounds are metallic in nature. The calculated elastic constants indicate that these compounds are mechanically stable at ambient pressure and found to be ductile in nature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗