Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Whole-remnant and maximum-voxel SPECT/CT dosimetry in {sup 131}I-NaI treatments of differentiated thyroid cancer

Purpose: To investigate the possible differences between SPECT/CT based whole-remnant and maximum-voxel dosimetry in patients receiving radio-iodine ablation treatment of differentiated thyroid cancer (DTC). Methods: Eighteen DTC patients were administered 1.11 GBq of {sup 131}I-NaI after near-total thyroidectomy and rhTSH stimulation. Two patients had two remnants, so in total dosimetry was performed for 20 sites. Three SPECT/CT scans were performed for each patient at 1, 2, and 3–7 days after administration. The activity, the remnant mass, and the maximum-voxel activity were determined from these images and from a recovery-coefficient curve derived from experimental phantom measurements. The cumulated activity was estimated using trapezoidal-exponential integration. Finally, the absorbed dose was calculated using S-values for unit-density spheres in whole-remnant dosimetry and S-values for voxels in maximum-voxel dosimetry. Results: The mean absorbed dose obtained from whole-remnant dosimetry was 40 Gy (range 2–176 Gy) and from maximum-voxel dosimetry 34 Gy (range 2–145 Gy). For any given patient, the activity concentrations for each of the three time-points were approximately the same for the two methods. The effective half-lives varied (R = 0.865), mainly due to discrepancies in estimation of the longer effective half-lives. On average, absorbed doses obtained from whole-remnant dosimetry were 1.2 ± 0.2 (1 SD) higher than for maximum-voxel dosimetry, mainly due to differences in the S-values. The method-related differences were however small in comparison to the wide range of absorbed doses obtained in patients. Conclusions: Simple and consistent procedures for SPECT/CT based whole-volume and maximum-voxel dosimetry have been described, both based on experimentally determined recovery coefficients. Generally the results from the two approaches are consistent, although there is a small, systematic difference in the absorbed dose due to differences in the S-values, and some variability due to differences in the estimated effective half-lives, especially when the effective half-life is long. Irrespective of the method used, the patient absorbed doses obtained span over two orders of magnitude.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Bayesian Monte-Carlo Evaluation Framework of Differential and Integral Data [Slides]

This presentation discusses the Bayesian Monte Carlo evaluation framework of differential and integral data, specifically in relation to two methods that are based on Bayes theorem, which assumes that all PDFs are normal (Gaussian) and that they use a linear approximation for IBEs. The presentation states that basic components of the MC evaluation framework include the computation of random MC ensemble from ENDF File 32, the simulation of IBEs and R-matrix cross section compared to experimental data, and the computation of weighted averages. It also states that the application to U-233 indicates deviation from the conventional linear approximation. In conclusion, evaluation framework will require MCMC method e.g., Metropolis-Hastings (M.-H.).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Capacitance–voltage modeling of mid-wavelength infrared nBn detectors

Capacitance–voltage measurements are a powerful technique to determine doping profiles of semiconductor pn junctions and Schottky barrier diodes. The measurements were recently extended to III-V-based mid-wavelength nBn infrared detectors, and absorber doping densities have been extracted using the widely accepted Schottky approximation, where the potential drop across the device is assumed to be across the depleting absorber layer. However, this approach is limited to when the absorber region of the nBn is under high reverse bias and thus is only able to extract the absorber region doping profile. Here, we introduce a semi-analytical model that is capable of extracting barrier dopant polarity, doping concentration, and thickness, as well as contact and absorber layer doping concentrations, all from a capacitance–voltage measurement. Rather than solely considering the potential drop across the depleting layers, it considers the potential drop across the accumulating layer as well. This negative charge accumulation occurs for the contact and absorber layers in the case of reverse and forward biases, respectively. This allows for a single model to be applied to a capacitance–voltage curve at forward and reverse biases and it can provide regions of bias where the absorber transitions from depletion to accumulation. We compare the agreement of the semianalytical model with modeling results from commercially available finite element method software and experimental capacitance–voltage data. Finally, we show that the method is consistent with the Schottky approximation of extracting absorber doping densities at high reverse bias and discuss the model's limitations.

Physics↗
DOE OSTI2013

The backward phase flow method for the Eulerian finite time Lyapunov exponent computations

We propose a simple Eulerian approach to compute the moderate to long time flow map for approximating the Lyapunov exponent of a (periodic or aperiodic) dynamical system. The idea is to generalize a recently proposed backward phase flow method which is specially designed for long time level set propagation. Unlike the original phase flow method or the backward phase flow method, which is applicable only to autonomous systems, the current approach can also be applied to any time-dependent (periodic or aperiodic) flow. We will discuss the stability of the proposed method. Numerical examples will be given to demonstrate the effectiveness of the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On the evaporation of solar dark matter: spin-independent effective operators

As a part of the effort to investigate the implications of dark matter (DM)-nucleon effective interactions on the solar DM detection, in this paper we focus on the evaporation of the solar DM for a set of the DM-nucleon spin-independent (SI) effective operators. In order to put the evaluation of the evaporation rate on a more reliable ground, we calculate the non-thermal distribution of the solar DM using the Monte Carlo methods, rather than adopting the Maxwellian approximation. We then specify relevant signal parameter spaces for the solar DM detection for various SI effective operators. Based on the analysis, we determine the minimum DM masses for which the DM-nucleon coupling strengths can be probed from the solar neutrino observations. As an interesting application, our investigation also shows that evaporation effect can not be neglectd in a recent proposal aiming to solve the solar abundance problem by invoking the momentum-dependent asymmetric DM in the Sun.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A fresh look at the photoelectron spectrum of bromobenzene: A third-order non-Dyson electron propagator study

The valence-shell ionization spectrum of bromobenzene, as a representative halogen substituted aromatic, was studied using the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator. This method, also referred to as IP-ADC(3), was implemented as a part of the Q-Chem program and enables large-scale calculations of the ionization spectra, where the computational effort scales as n{sup 5} with respect to the number of molecular orbitals n. The IP-ADC(3) scheme is ideally suited for investigating low-lying ionization transitions, so fresh insight could be gained into the cationic state manifold of bromobenzene. In particular, the present IP-ADC(3) calculations with the cc-pVTZ basis reveal a whole class of low-lying low-intensity two-hole-one-particle (2h-1p) doublet and quartet states, which are relevant to various photoionization processes. The good qualitative agreement between the theoretical spectral profile for the valence-shell ionization transitions generated with the smaller cc-pVDZ basis set and the experimental photoelectron spectrum measured at a photon energy of 80 eV on the PLÉIADES beamline at the Soleil synchrotron radiation source allowed all the main features to be assigned. Some theoretical aspects of the ionization energy calculations concerning the use of various approximation schemes and basis sets are discussed.

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

Spectroscopic Properties of KPb{sub 2}Cl{sub 5} and RbPb{sub 2}Br{sub 5} Doped with Er{sup 3+} and Yb{sup 3+}

Alkali halide crystals doped with Er{sup 3+} exhibit several properties that make them promising for practical use in optical fibers and converters of the emission spectrum. In this work, we studied the luminescence spectra of crystals as a function of the excitation radiation wavelength and the temperature dependences of the luminescence intensity in the up-conversion mechanism of excitation of the most intense lines (530, 550, and 650 nm) and upon excitation at 980 nm. In KPb{sub 2}Cl{sub 5}:Er{sup 3+}, RbPb{sub 2}Br{sub 5}:Er{sup 3+}, and KPb{sub 2}Cl{sub 5}:Er{sup 3+} + Yb{sup 3+} crystals, the luminescence kinetics were measured upon excitation at 313 and 365 nm, and the theoretical estimates of lifetimes were made using the Judd–Ofelt theory. The results can be used in the development of laser-based systems and for estimating the degree of approximation of calculations by the Judd–Ofelt method for these crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmas

Highlights: •A general approach to Monte Carlo methods for multicomponent plasmas is proposed. •We show numerical tests for the two-component (electrons and ions) case. •An optimal choice of parameters for speeding up the computations is discussed. •A rigorous estimate of the error of approximation is proved. -- Abstract: A general approach to Monte Carlo methods for Coulomb collisions is proposed. Its key idea is an approximation of Landau–Fokker–Planck equations by Boltzmann equations of quasi-Maxwellian kind. It means that the total collision frequency for the corresponding Boltzmann equation does not depend on the velocities. This allows to make the simulation process very simple since the collision pairs can be chosen arbitrarily, without restriction. It is shown that this approach includes the well-known methods of Takizuka and Abe (1977) [12] and Nanbu (1997) as particular cases, and generalizes the approach of Bobylev and Nanbu (2000). The numerical scheme of this paper is simpler than the schemes by Takizuka and Abe [12] and by Nanbu. We derive it for the general case of multicomponent plasmas and show some numerical tests for the two-component (electrons and ions) case. An optimal choice of parameters for speeding up the computations is also discussed. It is also proved that the order of approximation is not worse than O(√(ε)), where ε is a parameter of approximation being equivalent to the time step Δt in earlier methods. A similar estimate is obtained for the methods of Takizuka and Abe and Nanbu.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Hydrogen adsorption in metal-organic frameworks: The role of nuclear quantum effects

The role of nuclear quantum effects on the adsorption of molecular hydrogen in metal-organic frameworks (MOFs) has been investigated on grounds of Grand-Canonical Quantized Liquid Density-Functional Theory (GC-QLDFT) calculations. For this purpose, we have carefully validated classical H{sub 2}-host interaction potentials that are obtained by fitting Born-Oppenheimer ab initio reference data. The hydrogen adsorption has first been assessed classically using Liquid Density-Functional Theory and the Grand-Canonical Monte Carlo methods. The results have been compared against the semi-classical treatment of quantum effects by applying the Feynman-Hibbs correction to the Born-Oppenheimer-derived potentials, and by explicit treatment within the GC-QLDFT. The results are compared with experimental data and indicate pronounced quantum and possibly many-particle effects. After validation calculations have been carried out for IRMOF-1 (MOF-5), GC-QLDFT is applied to study the adsorption of H{sub 2} in a series of MOFs, including IRMOF-4, -6, -8, -9, -10, -12, -14, -16, -18, and MOF-177. Finally, we discuss the evolution of the H{sub 2} quantum fluid with increasing pressure and lowering temperature.

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

Vibrational and Thermal Properties of Oxyanionic Crystals

The vibrational and thermal properties of dolomite and alkali chlorates and perchlorates were studied in the gradient approximation of density functional theory using the method of a linear combination of atomic orbitals (LCAO). Long-wave vibration frequencies, IR and Raman spectra, and mode Gruneisen parameters were calculated. Equation-of-state parameters, thermodynamic potentials, entropy, heat capacity, and thermal expansion coefficient were also determined. The thermal expansion coefficient of dolomite was established to be much lower than for chlorates and perchlorates. The temperature dependence of the heat capacity at T > 200 K was shown to be generally governed by intramolecular vibrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Jacobi collocation scheme for variable-order fractional reaction-subdiffusion equation

We developed a numerical scheme to solve the variable-order fractional linear subdiffusion and nonlinear reaction-subdiffusion equations using the shifted Jacobi collocation method. Basically, a time-space collocation approximation for temporal and spatial discretizations is employed efficiently to tackle these equations. The convergence and stability analyses of the suggested basis functions are presented in-depth. The validity and efficiency of the proposed method are investigated and verified through numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Reducing the Cost of Energy Differences in Variational Monte Carlo with Spotlight Sampling

Here, we investigate an approximate sampling scheme that can significantly reduce the cost scaling of variational Monte Carlo when it is employed to predict the energy differences associated with local chemical changes. Inspired by side-chaining and embedding methods, this spotlight sampling approach adopts an approximate fragmented Hamiltonian and correlated sampling to reduce cost scaling to the point that it is essentially linear with system size, with the potential to go sublinear if certain conditions are met. In tests on bond stretching energies in alcohols, hydrogen dimer chains, and molecules with various degrees of π-system delocalization, we observe the anticipated linear scaling and an explicit cost crossover with standard variational Monte Carlo.

Bumann, Sonja [University of California, Berkeley,↗
DOE OSTI2015

Ab initio study of the structural, electronic and optical properties of ZnTe compound

Structural, electronic and optical properties of ZnTe compound were calculated using Density Functional Theory (DFT) based on the pseudopotentials and planewaves (PP-PW) method as implemented in the ABINIT computer code, where the exchange–correlation functional is approximated using the local density approximation (LDA) and the generalized gradient approximation (GGA). The obtained results from either LDA or GGa calculation for lattice parameter, energy band gap and optical parameters, such as the fundamental absorption edge, the peaks observed in the imaginary part of the dielectric function, the macroscopic dielectric constants and the optical dielectric constant, are compared with the available theoretical results and experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Optimal Energy Management for Microgrid with Stationary and Mobile Storage

This paper studies energy management in a Microgrid (MG) with solar generation, Battery Energy Management System (BESS) and gridable (V2G) EVs. A two-stage stochastic optimization method is proposed to capture the intermittent solar generation and random EV user behaviors. It is subsequently formulated as a Mixed Integer Linear Programming (MILP) problem. To evalutate the proposed method, real solar generation, loads, BESS and EV data is used in Sample Average Approximation (SAA). Computational results show the correctness of the proposed method as well as steady and tightly bounded optimality gap. Comparisons demonstrate that the proposed stochastic method outperforms its deterministic counterpart at the expense of higher computational cost. It is also observed that moderate number of EVs helps to reduce the overall operational cost of the MG, which sheds light on future EV integration to the smart grid.

Battery Energy Management System (BESS)↗
DOE OSTI2024

nuclear-score-maximization v1.0

This software library presents efficient and multithreaded implementations of matrix low rank approximation via column selection in C++17 code. The algorithms are described in Fornace, Mark, and Michael Lindsey. "Column and row subset selection using nuclear scores: algorithms and theory for Nystro m approximation, CUR decomposition, and graph Laplacian reduction." arXiv preprint arXiv:2407.01698 (2024). The presented methods are by-and-large ver novel, have provable approximation guarantees, multiple use-cases, and exhibit higher quality approximations on a variety of studied examples.

Fornace, Mark↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2016

Electronic and optical properties of mixed perovskites CsSn{sub x}Pb{sub (1−x)}I{sub 3}

The electronic structure and optical properties of the new solar cells absorber material: mixed perovskites CsSn{sub x}Pb{sub (1−x)}I{sub 3} are studied by the first-principle calculations with mBJ + SOC (modified Beak Johnson approximation plus spin-orbit coupling) method. The band gap of the serial of compounds almost quasi-linearly reduces with increasing Sn content from 0.96 eV (x = 0) to 0.16 eV (x = 1). Optical absorption coefficient revealed a progressive red shift with the increment of the Sn content, accompanying with the absorption edge broadening. The absorption coefficient and Ideal Power Absorption Coefficient (IPAC) increase greatly with the Pb atoms being partially substituted by Sn atoms. The pure CsSnI{sub 3} has the highest IPAC, but it is unstable in the air because the Sn{sup 2+} will be oxidized to Sn{sup 4+}. So our results indicate that partially substituted CsSn{sub x}Pb{sub (1−x)}I{sub 3} might be the good solar cell absorption material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗