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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Patient-Specific Monte Carlo-Based Dose-Kernel Approach for Inverse Planning in Afterloading Brachytherapy

Purpose: Brachytherapy planning software relies on the Task Group report 43 dosimetry formalism. This formalism, based on a water approximation, neglects various heterogeneous materials present during treatment. Various studies have suggested that these heterogeneities should be taken into account to improve the treatment quality. The present study sought to demonstrate the feasibility of incorporating Monte Carlo (MC) dosimetry within an inverse planning algorithm to improve the dose conformity and increase the treatment quality. Methods and Materials: The method was based on precalculated dose kernels in full patient geometries, representing the dose distribution of a brachytherapy source at a single dwell position using MC simulations and the Geant4 toolkit. These dose kernels are used by the inverse planning by simulated annealing tool to produce a fast MC-based plan. A test was performed for an interstitial brachytherapy breast treatment using two different high-dose-rate brachytherapy sources: the microSelectron iridium-192 source and the electronic brachytherapy source Axxent operating at 50 kVp. Results: A research version of the inverse planning by simulated annealing algorithm was combined with MC to provide a method to fully account for the heterogeneities in dose optimization, using the MC method. The effect of the water approximation was found to depend on photon energy, with greater dose attenuation for the lower energies of the Axxent source compared with iridium-192. For the latter, an underdosage of 5.1% for the dose received by 90% of the clinical target volume was found. Conclusion: A new method to optimize afterloading brachytherapy plans that uses MC dosimetric information was developed. Including computed tomography-based information in MC dosimetry in the inverse planning process was shown to take into account the full range of scatter and heterogeneity conditions. This led to significant dose differences compared with the Task Group report 43 approach for the Axxent source.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Monte Carlo Second- and Third-Order Many-Body Green’s Function Methods with Frequency-Dependent, Nondiagonal Self-Energy

We fully develop the Monte Carlo many-body Green’s function (MC-GF) method with the following enhancements: (1) The truncation order of the perturbation expansion of the Dyson self-energy is raised from the second order (MC-GF2) to the third order (MC-GF3) with the aid of a computerized procedure to enumerate and transform all 84 third-order Goldstone diagrams into Monte Carlo integrable expressions and then into central processing unit (CPU)/graphical processing unit (GPU)-parallel computer codes. (2) An efficient algorithm is proposed that computes all off-diagonal and diagonal elements of the MC-GF2 and MC-GF3 self-energy matrices by common subexpression elimination. (3) The frequency-independent approximation is lifted by introducing a method that computes frequency derivatives of the MC-GF2 and MC-GF3 self-energies up to any arbitrarily high order at nearly no additional computational cost. (4) The imaginary-time integration in the Laplace-transformed expressions of the self-energy is carried out stochastically (instead of using a quadrature in the previous implementations), resulting in a 50- to 200-fold speedup. (5) The efficiency of the redundant-walker convergence acceleration scheme is analyzed numerically, and the guidelines are established to select an optimal number of walkers for maximal efficiency. When such an optimal number is used, the cost per sample is constant of molecular size on either many CPUs or many GPUs. (6) The computational cost to obtain a binding energy within a given statistical uncertainty is observed to increase as (tentatively) O(n 4 ) and O(n 5 ) of molecular size (n is the number of electrons) with and without the diagonal approximation, respectively, as opposed to O(n 5 ) and O(n 6 ) of the corresponding deterministic algorithms. With this method applied to the electron binding energies of C 60 , we show that the third-order corrections to the self-energies are much greater in electron binding energies than in ground-state energies. As a result, they display a sign of oscillatory convergence toward experimental results, not necessarily improving the agreement with increasing perturbation order, justifying MC-GF3 and motivating even higher-order methods.

Algorithms↗
DOE OSTI2023

Dimensionality reduction using elastic measures

With the recent surge in big data analytics for hyperdimensional data, there is a renewed interest in dimensionality reduction techniques. In order for these methods to improve performance gains and understanding of the underlying data, a proper metric needs to be identified. This step is often overlooked, and metrics are typically chosen without consideration of the underlying geometry of the data. Here, in this paper, we present a method for incorporating elastic metrics into the t-distributed stochastic neighbour embedding (t-SNE) and Uniform Manifold Approximation and Projection (UMAP). We apply our method to functional data, which is uniquely characterized by rotations, parameterization and scale. If these properties are ignored, they can lead to incorrect analysis and poor classification performance. Through our method, we demonstrate improved performance on shape identification tasks for three benchmark data sets (MPEG-7, Car data set and Plane data set of Thankoor), where we achieve 0.77, 0.95 and 1.00 F1 score, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

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

Data-driven Distributed Learning of Multi-agent Systems: A Koopman Operator Approach

Koopman operator theory provides a model-free technique for studying nonlinear dynamical systems purely from data. Since the Koopman operator is infinite-dimensional, researchers have developed several methods that provide a finite-dimensional approximation of the Koopman operator so that it can be applied for practical use cases. One common thing with most of the methods is that their solutions are obtained by solving a centralized minimization problem. In this work, we treat the dynamical system to be a multi-agent system and propose an algorithm to compute the finite-dimensional approximation of the Koopman operator in a distributed manner using the knowledge of the topology of the underlying multi-agent system. The proposed distributed approach is shown to be equivalent to the centralized learning problem and results in a sparse Koopman whose block structure mimics the Laplacian of the multi-agent system. Extensive simulation studies illustrate the proposed framework on the network of oscillators and the IEEE 68 bus system.

Nandanoori, Sai Pushpak↗
DOE OSTI2013

BAYESIAN INFERENCE OF POLARIZED COSMIC MICROWAVE BACKGROUND POWER SPECTRA FROM INTERFEROMETRIC DATA

Detection of B-mode polarization of the cosmic microwave background (CMB) radiation is one of the frontiers of observational cosmology. Because they are an order of magnitude fainter than E-modes, it is quite a challenge to detect B-modes. Having more manageable systematics, interferometers prove to have a substantial advantage over imagers in detecting such faint signals. Here, we present a method for Bayesian inference of power spectra and signal reconstruction from interferometric data of the CMB polarization signal by using the technique of Gibbs sampling. We demonstrate the validity of the method in the flat-sky approximation for a simulation of an interferometric observation on a finite patch with incomplete uv-plane coverage, a finite beam size, and a realistic noise model. With a computational complexity of O(n {sup 3/2}), n being the data size, Gibbs sampling provides an efficient method for analyzing upcoming cosmology observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Reduced-order modeling through machine learning and graph-theoretic approaches for brittle fracture applications

Typically, thousands of computationally expensive micro-scale simulations of brittle crack propagation are needed to upscale lower length scale phenomena to the macro-continuum scale. Running such a large number of crack propagation simulations presents a significant computational challenge, making reduced-order models (ROMs) attractive for this task. Here, the ultimate goal of this research is to develop ROMs that have sufficient accuracy and low computational cost so that these upscaling simulations can be readily performed. However, constructing ROMs for these complex simulations presents its own challenge. Here, we present and compare four different approaches for reduced-order modeling of brittle crack propagation in geomaterials. These methods rely on machine learning (ML) and graph-theoretic algorithms to approximate key aspects of the brittle crack problem. These methods also incorporate different physics-based assumptions in order to reduce the training requirements while maintaining accurate physics as much as possible. Results from the ROMs are directly compared against a high-fidelity model of brittle crack propagation. Further, the strengths and weaknesses of the ROMs are discussed, and we conclude that combining smart physics-informed feature engineering with highly trainable ML models provides the best performance. The ROMs considered here have computational costs that are orders-of-magnitude less than the cost associated with high-fidelity physical models while maintaining good accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Transition metal dimer on Au(111) surface: A first principle study

The adsorption behaviour of transition metal dimers M{sub 2} (M= Cu, Ag, Au) on the Au(111) surface have been studied using the density functional theory formalism. The projector augmented wave method under the spin polarized version of generalized gradient approximation scheme was employed to calculate the total energy. The results suggest that all dimers prefer to orient in parallel to the surface plane, where two M atoms are adsorbed on two nearby threefold fcc sites. We have investigated the chemical interaction between M atoms and Au surface through electronic density of state analysis. It is found that on deposition, the electronic density of states (EDOS) of M{sub 2} dimer becomes broader in comparison to their gas phase spectrum.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Study of structural and electronic properties of ScN and ScAs in rocksalt and zincblende structure: A DFT approach

In this paper, we have studied the structural and electronic properties of ScN and ScAs in zincblende (ZB) and rocksalt (RS) phases. We have employed the full potential linearized augmented plane wave (FP-APW) method within the density functional theory (DFT). Generalized gradient approximation (GGA), due to Perdew-Burke-Ernzerhof (PBE) has been used to estimate the exchange-correlation functional. Our band structure results for ScN shows the metallic nature, whereas ScAs shows the semiconducting behavior. The obtained results are in excellent agreement with earlier reported data.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Angular oversampling with temporally offset layers on multilayer detectors in computed tomography

Purpose: Today’s computed tomography (CT) scanners operate at an increasingly high rotation speed in order to reduce motion artifacts and to fulfill the requirements of dynamic acquisition, e.g., perfusion and cardiac imaging, with lower angular sampling rate as a consequence. In this paper, a simple method for obtaining angular oversampling when using multilayer detectors in continuous rotation CT is presented. Methods: By introducing temporal offsets between the measurement periods of the different layers on a multilayer detector, the angular sampling rate can be increased by a factor equal to the number of layers on the detector. The increased angular sampling rate reduces the risk of producing aliasing artifacts in the image. A simulation of a detector with two layers is performed to prove the concept. Results: The simulation study shows that aliasing artifacts from insufficient angular sampling are reduced by the proposed method. Specifically, when imaging a single point blurred by a 2D Gaussian kernel, the method is shown to reduce the strength of the aliasing artifacts by approximately an order of magnitude. Conclusions: The presented oversampling method is easy to implement in today’s multilayer detectors and has the potential to reduce aliasing artifacts in the reconstructed images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Stochastic Reductions for Inertial Fluid-Structure Interactions Subject to Thermal Fluctuations

We explore the dynamics of elastic microstructures within a fluid that are subjected to thermal fluctuations. We perform analysis to obtain systematically simplified descriptions of the mechanics in the limiting regimes when (i) the coupling forces that transfer momentum between the fluid and microstructures are strong, (ii) the mass of the microstructures is small relative to the displaced mass of the fluid, and (iii) the response to stresses results in hydrodynamics that relax rapidly to a quasi-steady-state relative to the motions of the microstructure. Furthermore, we derive effective equations using a singular perturbation analysis of the backward Kolmogorov equations of the stochastic process. Our continuum mechanics description is based on the stochastic Eulerian--Lagrangian method which provides a framework for approximation of the fluid-structure interactions when subject to thermal fluctuations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics

A new time-dependent neutron transport method based on the method of characteristics (MOC) has been developed. Whereas most spatial kinetics methods treat time dependence through temporal discretization, this new method treats time dependence by defining the characteristics to span space and time. In this implementation regions are defined in space-time where the thickness of the region in time fulfills an analogous role to the time step in discretized methods. The time dependence of the local source is approximated using a truncated Taylor series expansion with high order derivatives approximated using backward differences, permitting the solution of the resulting space-time characteristic equation. To avoid a drastic increase in computational expense and memory requirements due to solving many discrete characteristics in the space-time planes, the temporal variation of the boundary source is similarly approximated. This allows the characteristics in the space-time plane to be represented analytically rather than discretely, resulting in an algorithm comparable in implementation and expense to one that arises from conventional time integration techniques. Furthermore, by defining the boundary flux time derivative in terms of the preceding local source time derivative and boundary flux time derivative, the need to store angularly-dependent data is avoided without approximating the angular dependence of the angular flux time derivative. The accuracy of this method is assessed through implementation in the neutron transport code DeCART. The method is employed with variable-order local source representation to model a TWIGL transient. The results demonstrate that this method is accurate and more efficient than the discretized method. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Attribution of PM 2.5 exposure in Beijing–Tianjin–Hebei region to emissions: implication to control strategies

The Beijing–Tianjin–Hebei (BTH) region is one of the most heavily polluted regions in China, with both high PM 2.5 concentrations and a high population density. A quantitative source-receptor relationship can provide valuable insights that can inform effective emission control strategies. Both source apportionment (SA) and source sensitivity (SS) can provide such information from different perspectives. In this study, both methods are applied in northern China to identify the most significant emission categories and source regions for PM 2.5 exposure in BTH in 2013. Despite their differences, both models show similar distribution patterns for population and simulated PM 2.5 concentrations, resulting in overall high PM 2.5 exposure values (approximately 110 μg/m 3 ) and particularly high exposure values during the winter (approximately 200 μg/m 3 ). Both methods show that local emissions play a dominant role (70%), with some contribution from surrounding provinces (e.g., Shandong) via regional transport. The two methods also agree on the priority of local emission controls: both identify industrial, residential, and agricultural emissions as the top three categories that should be controlled locally. In addition, the effect of controlling agricultural ammonia emissions is approximately doubled when the co-benefits of reducing nitrate are considered. In conclusion, the synthesis of SA and SS for addressing specific categories of emissions provides a quantitative basis for the development of emission control strategies and policies for controlling PM 2.5 in China.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

Sensitivity analysis for dose deposition in radiotherapy via a Fokker–Planck model

In this paper, we study the sensitivities of electron dose calculations with respect to stopping power and transport coefficients. We focus on the application to radiotherapy simulations. We use a Fokker–Planck approximation to the Boltzmann transport equation. Equations for the sensitivities are derived by the adjoint method. The Fokker–Planck equation and its adjoint are solved numerically in slab geometry using the spherical harmonics expansion (P N ) and an Harten-Lax-van Leer finite volume method. Our method is verified by comparison to finite difference approximations of the sensitivities. Finally, we present numerical results of the sensitivities for the normalized average dose deposition depth with respect to the stopping power and the transport coefficients, demonstrating the increase in relative sensitivities as beam energy decreases. In conclusion, this in turn gives estimates on the uncertainty in the normalized average deposition depth, which we present.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Multiphase Optimal and Non-Singular Power Flow by Successive Linear Approximations: Preprint

We propose an optimal power flow method for multiphase power systems. The method works for both radial and meshed networks, and is able to accommodate different types of load/source connections such as wye, delta, etc. Moreover, the method generates a sequence of non-singular points in the electrical state space and a corresponding sequence of power injections, which can be used in real-time control settings. Under certain conditions, the sequences attain the local minimum. The method has been tested using several typical networks, including the IEEE 37-bus and 123-bus test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing

Angular two-point statistics of large-scale structure observables are important cosmological probes. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation; the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. The new method significantly simplifies the calculation and improves the numerical speed and stability. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. We find that using the Limber approximation for galaxy clustering in future analyses like LSST Year 1 and DES Year 6 may cause significant biases in cosmological parameters, indicating that going beyond the Limber approximation is necessary for these analyses.

weak gravitational lensing↗
DOE OSTI2014

Cholesky-decomposed density MP2 with density fitting: Accurate MP2 and double-hybrid DFT energies for large systems

Our recently developed QQR-type integral screening is introduced in our Cholesky-decomposed pseudo-densities Møller-Plesset perturbation theory of second order (CDD-MP2) method. We use the resolution-of-the-identity (RI) approximation in combination with efficient integral transformations employing sparse matrix multiplications. The RI-CDD-MP2 method shows an asymptotic cubic scaling behavior with system size and a small prefactor that results in an early crossover to conventional methods for both small and large basis sets. We also explore the use of local fitting approximations which allow to further reduce the scaling behavior for very large systems. The reliability of our method is demonstrated on test sets for interaction and reaction energies of medium sized systems and on a diverse selection from our own benchmark set for total energies of larger systems. Timings on DNA systems show that fast calculations for systems with more than 500 atoms are feasible using a single processor core. Parallelization extends the range of accessible system sizes on one computing node with multiple cores to more than 1000 atoms in a double-zeta basis and more than 500 atoms in a triple-zeta basis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗