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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Self-consistent continuum random-phase approximation calculations of {sup 4}He electromagnetic responses

We study the electromagnetic responses of {sup 4}He within the framework of the self-consistent continuum random-phase approximation theory. In this approach, the ground-state properties are described by a Hartree-Fock calculation. The single-particle basis constructed in this manner is used in the calculations of the continuum responses of the system. Finite-range interactions are considered in the calculations. We compare our results with photon-absorption cross sections and electron-scattering quasielastic data. From this comparison, and also from the comparison with the results of microscopic calculations, we deduce that our approach describes well the continuum excitation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Practical aspects and uncertainty analysis of biological effective dose (BED) regarding its three-dimensional calculation in multiphase radiotherapy treatment plans

Purpose: There is a growing interest in the radiation oncology community to use the biological effective dose (BED) rather than the physical dose (PD) in treatment plan evaluation and optimization due to its stronger correlation with radiobiological effects. Radiotherapy patients may receive treatments involving a single only phase or multiple phases (e.g., primary and boost). Since most treatment planning systems cannot calculate the analytical BED distribution in multiphase treatments, an approximate multiphase BED expression, which is based on the total physical dose distribution, has been used. The purpose of this paper is to reveal the mathematical properties of the approximate BED formulation, relative to the true BED. Methods: The mathematical properties of the approximate multiphase BED equation are analyzed and evaluated. In order to better understand the accuracy of the approximate multiphase BED equation, the true multiphase BED equation was derived and the mathematical differences between the true and approximate multiphase BED equations were determined. The magnitude of its inaccuracies under common clinical circumstances was also studied. All calculations were performed on a voxel-by-voxel basis using the three-dimensional dose matrices. Results: Results showed that the approximate multiphase BED equation is accurate only when the dose-per-fractions (DPFs) in both the first and second phases are equal, which occur when the dose distribution does not significantly change between the phases. In the case of heterogeneous dose distributions, which significantly vary between the phases, there are fewer occurrences of equal DPFs and hence the inaccuracy of the approximate multiphase BED is greater. These characteristics are usually seen in the dose distributions being delivered to organs at risk rather than to targets. Conclusions: The finding of this study indicates that the true multiphase BED equation should be implemented in the treatment planning systems due to the inconsistent accuracy of the approximate multiphase BED equation in most of the clinical situations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Probing Ground-State Phase Transitions through Quench Dynamics

The study of quantum phase transitions requires the preparation of a many-body system near its ground state, a challenging task for many experimental systems. The measurement of quench dynamics, on the other hand, is now a routine practice in most cold atom platforms. Here we show that quintessential ingredients of quantum phase transitions can be probed directly with quench dynamics in integrable and nearly integrable systems. As a paradigmatic example, we study global quench dynamics in a transverse-field Ising model with either short-range or long-range interactions. When the model is integrable, we discover a new dynamical critical point with a nonanalytic signature in the short-range correlators. The location of the dynamical critical point matches that of the quantum critical point and can be identified using a finite-time scaling method. We extend this scaling picture to systems near integrability and demonstrate the continued existence of a dynamical critical point detectable at prethermal timescales. We quantify the difference in the locations of the dynamical and quantum critical points away from (but near) integrability. Thus, we demonstrate that this method can be used to approximately locate the quantum critical point near integrability. Furthermore, the scaling method is also relevant to experiments with finite time and system size, and our predictions are testable in near-term experiments with trapped ions and Rydberg atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Teleoperation for Nuclear Applications Considering Communication Delay and Uncertainty - 18205

Teleoperation in remote environment is of great interest in robotics. Advancement in communication technology has enabled safe execution of complex and life threatening tasks in network environment in teleoperation. Using bilateral teleoperation, the operator receives force feedback from the slave which provides an interactive feel of forces and smoother coordination. However, time delay across the communication link in sending the command signal from the master to the slave may cause instability to the system. In this research, we use the continuous pole placement method to address the time delay problem in teleoperation. Compared to approximation based and prediction based methods, this method shows improved robustness. Simulations and experiments show a stable and well synchronized control. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2012

HARMONIC IN-PAINTING OF COSMIC MICROWAVE BACKGROUND SKY BY CONSTRAINED GAUSSIAN REALIZATION

The presence of astrophysical emissions between the last scattering surface and our vantage point requires us to apply a foreground mask on cosmic microwave background (CMB) sky maps, leading to large cuts around the Galactic equator and numerous holes. Since many CMB analysis, in particular on the largest angular scales, may be performed on a whole-sky map in a more straightforward and reliable manner, it is of utmost importance to develop an efficient method to fill in the masked pixels in a way compliant with the expected statistical properties and the unmasked pixels. In this Letter, we consider the Monte Carlo simulation of a constrained Gaussian field and derive it CMB anisotropy in harmonic space, where a feasible implementation is possible with good approximation. We applied our method to simulated data, which shows that our method produces a plausible whole-sky map, given the unmasked pixels, and a theoretical expectation. Subsequently, we applied our method to the Wilkinson Microwave Anisotropy Probe foreground-reduced maps and investigated the anomalous alignment between quadrupole and octupole components. From our investigation, we find that the alignment in the foreground-reduced maps is even higher than the Internal Linear Combination map. We also find that the V-band map has higher alignment than other bands, despite the expectation that the V-band map has less foreground contamination than other bands. Therefore, we find it hard to attribute the alignment to residual foregrounds. Our method will be complementary to other efforts on in-painting or reconstructing the masked CMB data, and of great use to Planck surveyor and future missions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

On the Statistical Calibration of Physical Models

ABSTRACT We introduce a novel statistical calibration framework for physical models, relying on probabilistic embedding of model discrepancy error within the model. For clarity of illustration, we take the measurement errors out of consideration, calibrating a chemical model of interest with respect to a more detailed model, considered as “truth” for the present purpose. We employ Bayesian statistical methods for such model‐to‐model calibration and demonstrate their capabilities on simple synthetic models, leading to a well‐defined parameter estimation problem that employs approximate Bayesian computation. The method is then demonstrated on two case studies for calibration of kinetic rate parameters for methane air chemistry, where ignition time information from a detailed elementary‐step kinetic model is used to estimate rate coefficients of a simple chemical mechanism. We show that the calibrated model predictions fit the data and that uncertainty in these predictions is consistent in a mean‐square sense with the discrepancy from the detailed model data.

Sargsyan, K.↗
DOE OSTI2015

Study of multiband disordered systems using the typical medium dynamical cluster approximation

We generalize the typical medium dynamical cluster approximation to multiband disordered systems. Using our extended formalism, we perform a systematic study of the nonlocal correlation effects induced by disorder on the density of states and the mobility edge of the three-dimensional two-band Anderson model. We include interband and intraband hopping and an intraband disorder potential. Our results are consistent with those obtained by the transfer matrix and the kernel polynomial methods. We also apply the method to K x Fe 2-y Se 2 with Fe vacancies. Despite the strong vacancy disorder and anisotropy, we find the material is not an Anderson insulator. Moreover our results demonstrate the application of the typical medium dynamical cluster approximation method to study Anderson localization in real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Identifying Important Ions and Positions in Mass Spectrometry Imaging Data Using CUR Matrix Decompositions

Mass spectrometry imaging enables label-free, high-resolution spatial mapping of the chemical composition of complex, biological samples. Typical experiments require selecting ions and/or positions from the images: ions for fragmentation studies to identify keystone compounds and positions for follow up validation measurements using microdissection or other orthogonal techniques. Unfortunately, with modern imaging machines, these must be selected from an overwhelming amount of raw data. Existing techniques to reduce the volume of data, the most popular of which are principle component analysis and non-negative matrix factorization, have the disadvantage that they return difficult-to-interpret linear combinations of actual data elements. In this work, we show that CX and CUR matrix decompositions can be used directly to address this selection need. CX and CUR matrix decompositions use empirical statistical leverage scores of the input data to provide provably good low-rank approximations of the measured data that are expressed in terms of actual ions and actual positions, as opposed to difficult-to-interpret eigenions and eigenpositions. We show that this leads to effective prioritization of information for both ions and positions. In particular, important ions can be found either by using the leverage scores as a ranking function and using a deterministic greedy selection algorithm or by using the leverage scores as an importance sampling distribution and using a random sampling algorithm; however, selection of important positions from the original matrix performed significantly better when they were chosen with the random sampling algorithm. Also, we show that 20 ions or 40 locations can be used to reconstruct the original matrix to a tolerance of 17% error for a widely studied image of brain lipids; and we provide a scalable implementation of this method that is applicable for analysis of the raw data where there are often more than a million rows and/or columns, which is larger than SVD-based low-rank approximation methods can handle. Lastly, these results introduce the concept of CX/CUR matrix factorizations to mass spectrometry imaging, describing their utility and illustrating principled algorithmic approaches to deal with the overwhelming amount of data generated by modern mass spectrometry imaging.

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

Quantum-enhanced phase imaging without coincidence counting

Quantitative phase imaging provides a way to image transparent objects, such as biological cells, and measure their thickness. We report on a phase-imaging method that achieves twice the phase shift and approximately 1.7 times the spatial resolution of an equivalent spatially and temporally coherent classical quantitative phase-imaging system by using quantum interference between successive spontaneous parametric downconversion events in a nonlinear crystal. Furthermore, our method is approximately 1000 times faster than imaging the parametric downconversion photons in coincidence, which requires measurement times on the order of tens of hours. Our method may be useful for imaging sensitive transparent objects that require low illumination intensities at near-infrared and longer illumination wavelengths, such as photosensitive biological samples.

47 OTHER INSTRUMENTATION↗
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 OSTI2011

Schrieffer-Wolff transformation for quantum many-body systems

The Schrieffer-Wolff (SW) method is a version of degenerate perturbation theory in which the low-energy effective Hamiltonian H{sub eff} is obtained from the exact Hamiltonian by a unitary transformation decoupling the low-energy and high-energy subspaces. We give a self-contained summary of the SW method with a focus on rigorous results. We begin with an exact definition of the SW transformation in terms of the so-called direct rotation between linear subspaces. From this we obtain elementary proofs of several important properties of H{sub eff} such as the linked cluster theorem. We then study the perturbative version of the SW transformation obtained from a Taylor series representation of the direct rotation. Our perturbative approach provides a systematic diagram technique for computing high-order corrections to H{sub eff}. We then specialize the SW method to quantum spin lattices with short-range interactions. We establish unitary equivalence between effective low-energy Hamiltonians obtained using two different versions of the SW method studied in the literature. Finally, we derive an upper bound on the precision up to which the ground state energy of the nth-order effective Hamiltonian approximates the exact ground state energy. - Highlights: > The Schrieffer-Wolff transformation is specialized to quantum spin lattices with short-range interactions. > We provide a diagram technique for computing high-order corrections to the effective low-energy Hamiltonian. > We derive a rigorous bound on the error up to which the nth-order effective low-energy dynamics approximates the exact dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2014

Coupled cluster channels in the homogeneous electron gas

We discuss diagrammatic modifications to the coupled cluster doubles (CCD) equations, wherein different groups of terms out of rings, ladders, crossed-rings, and mosaics can be removed to form approximations to the coupled cluster method, of interest due to their similarity with various types of random phase approximations. The finite uniform electron gas (UEG) is benchmarked for 14- and 54-electron systems at the complete basis set limit over a wide density range and performance of different flavours of CCD is determined. These results confirm that rings generally overcorrelate and ladders generally undercorrelate; mosaics-only CCD yields a result surprisingly close to CCD. We use a recently developed numerical analysis [J. J. Shepherd and A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013)] to study the behaviours of these methods in the thermodynamic limit. We determine that the mosaics, on forming the Brueckner one-body Hamiltonian, open a gap in the effective one-particle eigenvalues at the Fermi energy. Numerical evidence is presented which shows that methods based on this renormalisation have convergent energies in the thermodynamic limit including mosaic-only CCD, which is just a renormalised MP2. All other methods including only a single channel, namely, ladder-only CCD, ring-only CCD, and crossed-ring-only CCD, appear to yield divergent energies; incorporation of mosaic terms prevents this from happening.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗