Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Dark-matter harmonics beyond annual modulation

The count rate at dark-matter direct-detection experiments should modulate annually due to the motion of the Earth around the Sun. We show that higher-frequency modulations, including daily modulation, are also present and in some cases are nearly as strong as the annual modulation. These higher-order modes are particularly relevant if (i) the dark matter is light, O(10) GeV, (ii) the scattering is inelastic, or (iii) velocity substructure is present; for these cases, the higher-frequency modes are potentially observable at current and ton-scale detectors. We derive simple expressions for the harmonic modes as functions of the astrophysical and geophysical parameters describing the Earth's orbit, using an updated expression for the Earth's velocity that corrects a common error in the literature. For an isotropic halo velocity distribution, certain ratios of the modes are approximately constant as a function of nuclear recoil energy. Anisotropic distributions can also leave observable features in the harmonic spectrum. Consequently, the higher-order harmonic modes are a powerful tool for identifying a potential signal from interactions with the Galactic dark-matter halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Avoidance of a Landau pole by flat contributions in QED

We consider massless Quantum Electrodynamics in the momentum scheme and carry forward an approach based on Dyson–Schwinger equations to approximate both the β-function and the renormalized photon self-energy (Yeats, 2011). Starting from the Callan–Symanzik equation, we derive a renormalization group (RG) recursion identity which implies a non-linear ODE for the anomalous dimension and extract a sufficient but not necessary criterion for the existence of a Landau pole. This criterion implies a necessary condition for QED to have no such pole. Solving the differential equation exactly for a toy model case, we integrate the corresponding RG equation for the running coupling and find that even though the β-function entails a Landau pole it exhibits a flat contribution capable of decreasing its growth, in other cases possibly to the extent that such a pole is avoided altogether. Finally, by applying the recursion identity, we compute the photon propagator and investigate the effect of flat contributions on both spacelike and timelike photons. -- Highlights: •We present an approach to approximate both the β-function and the photon self-energy. •We find a sufficient criterion for the self-energy to entail the existence of a Landau pole. •We study non-perturbative ‘flat’ contributions that emerge within the context of our approach. •We discuss a toy model and how it is affected by flat contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗
DOE OSTI2024

Investigation of Ethane Dehydrogenation and Hydrogenolysis on Pt(111), Pt(211), and Pt(100): Bayesian Quantification and Correction of DFT-Based Enthalpic and Entropic Uncertainties

Computational investigations of heterogeneously catalyzed reactions using density functional theory (DFT) are often inaccurate, largely due to uncertainties in the choice of DFT functional (enthalpic uncertainty) and approximations for modeling adsorbate movement along the catalyst surface (entropic uncertainty). This work illustrates that both uncertainties are significant in the investigation of ethane dehydrogenation (EDH) and hydrogenolysis on Pt catalysts by considering the complete deconstruction of ethane on Pt(111), Pt(211), and Pt(100) using microkinetic modeling (MKM). Hence, this work uses both noncalibrated and Bayesian-calibrated MKMs to quantify and correct inaccuracies in macroscopic properties due to both uncertainties. A Bayesian approach to the correction of entropic errors was introduced using a “Modified Fermi Function (MFF)” to calibrate between the two bounds of entropy represented by the harmonic oscillator (HO) and free translator (FT) approximations. Regardless of enthalpic and entropic uncertainties, all three surfaces are capable of ethane activation; however, Pt(211) was found to be the most active and is largely responsible for methane production. Next, Pt(111) is largely responsible for acetylene production, and Pt(100) has the highest ethylene selectivity but is most susceptible to coking. By comparison of different calibrated models, the FT entropy approximation was found to better describe EDH under typical experimental conditions. Statistical evidence was found to support Pt(111) as the active site for EDH, assuming that one single site is responsible for the chemistry. On the three surfaces, competing second dehydrogenations to CH 2 CH 2 and CH 3 CH were observed as well as isomerization of CH 3 CH back to CH 2 CH 2 and deeper dehydrogenation of CH 3 CH. In conclusion, C–C cleavage was found to largely proceed via the CH 3 C intermediate on Pt(100) and Pt(111), while on Pt(211), it was via both CHC and CH 3 C.

Bayesian model selection↗
DOE OSTI2011

Electron-vibron coupling in halogenated acenaphthenequinone upon O K-edge soft x-ray absorption

We present high resolution oxygen K near-edge x-ray absorption spectra of the acenaphthenequinone (ANQ) derivative 3,8-dibromo-5,6-dichloro-ANQ (Br{sub 2}Cl{sub 2}-ANQ). The spectral features exhibit an almost identical vibronic fine structure compared to that shown by ANQ. The unequal distances of the vibronic levels as derived from the least-squares fit analysis of the vibronic progressions allows us to determine the anharmonicity of the excited state potentials involved. We conclude that a single vibrational progression couples to the resonant excitation of O 1s core electrons preferentially. Comparison of the two ANQ derivatives gives a clear indication that the vibronic mode corresponds to a C=O stretching mode, rather than coupling to a C-H mode as suggested previously. These conclusions are supported by density functional theory calculations.

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

Link between optical spectra, crystal-field parameters, and local environments of Eu{sup 3+} ions in Eu{sub 2}O{sub 3}-doped sodium disilicate glass

Rare-earth-doped glasses are key materials for optical technology due to the luminescent properties of 4f{sup n} ions. The crystal-field model describes the effect of local environment on transitions between 4f electrons. We present a detailed modeling study of the optical spectra of sodium disilicate glass, 33Na{sub 2}O{center_dot}67SiO{sub 2}, doped with 0.2% and 1.0 mol%Eu{sub 2}O{sub 3}. This study uses very large molecular dynamics models with up to 100 Eu{sup 3+} ions, the superposition model for covalent and overlap effects on crystal-field parameters, and realistic values for homogeneous linewidth broadening. The simulated spectra are in reasonable agreement with experiment. The trends in {sup 7} F{sub J} energy levels across different Eu{sup 3+} ion sites have been examined and a very detailed analysis is presented that looks at how features of the spectra are related to features of the local environment of Eu{sup 3+} ions. Increasing the crystal-field strength S{sub total} causes the {sup 7} F{sub 0} energy level to decrease and causes the splitting of {sup 7} F{sub J} manifolds to increase, and this is due to increasing mixing of 4f wave functions. To a reasonable approximation the crystal-field strength components S{sub k} depend on angular positions of ligands independently of distances to ligands. The former are seen to be more significant in determining S{sub k}, which are closely related to the rotationally invariant bond-orientational order parameters Q{sub k}. The values of S{sub 2} are approximately linear in Q{sub 2}, and the values of Q{sub 2} are higher for fivefold than sixfold coordinated rare-earth ions. These results can be of importance for efforts to enhance the local environment of rare-earth ions in oxide glasses for optical applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Benchmark calculations on the lowest-energy singlet, triplet, and quintet states of the four-electron harmonium atom

For a wide range of confinement strengths ω, explicitly-correlated calculations afford approximate energies E(ω) of the ground and low-lying excited states of the four-electron harmonium atom that are within few μhartree of the exact values, the errors in the respective energy components being only slightly higher. This level of accuracy constitutes an improvement of several orders of magnitude over the previously published data, establishing a set of benchmarks for stringent calibration and testing of approximate electronic structure methods. Its usefulness is further enhanced by the construction of differentiable approximants that allow for accurate computation of E(ω) and its components for arbitrary values of ω. The diversity of the electronic states in question, which involve both single- and multideterminantal first-order wavefunctions, and the availability of the relevant natural spinorbitals and their occupation numbers make the present results particularly useful in research on approximate density-matrix functionals. The four-electron harmonium atom is found to possess the {sup 3}P{sub +} triplet ground state at strong confinements and the {sup 5}S{sub −} quintet ground state at the weak ones, the energy crossing occurring at ω ≈ 0.0240919.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Comparing first-principles density functionals plus corrections for the lattice dynamics of YBa 2 Cu 3 O 6

The enigmatic mechanism underlying unconventional high-temperature superconductivity, especially the role of lattice dynamics, has remained a subject of debate. Theoretical insights have long been hindered due to the lack of an accurate first-principles description of the lattice dynamics of cuprates. Recently, using the r2SCAN meta-generalized gradient approximation (meta-GGA) functional, we have been able to achieve accurate phonon spectra of an insulating cuprate YBa 2 Cu 3 O 6 and discover significant magnetoelastic coupling in experimentally interesting Cu–O bond stretching optical modes [Ning et al., Phys. Rev. B 107, 045126 (2023)]. We extend this work by comparing Perdew–Burke–Ernzerhof and r2SCAN performances with corrections from the on-site Hubbard U and the D4 van der Waals (vdW) methods, aiming at further understanding on both the materials science side and the density functional side. We demonstrate the importance of vdW and self-interaction corrections for accurate first-principles YBa 2 Cu 3 O 6 lattice dynamics. Since r2SCAN by itself partially accounts for these effects, the good performance of r2SCAN is now more fully explained. Finally, in addition, the performances of the Tao–Mo series of meta-GGAs, which are constructed in a different way from the strongly constrained and appropriately normed (SCAN) meta-GGA and its revised version r2SCAN, are also compared and discussed.

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

Spin currents of charged Dirac particles in rotating coordinates

The semiclassical Boltzmann transport equation of charged, massive fermions in a rotating frame of reference, in the presence of external electromagnetic fields is solved in the relaxation time approach to establish the distribution function up to linear order in the electric field in rotating coordinates, centrifugal force and the derivatives. The spin and spin current densities are calculated by means of this distribution function at zero temperature up to the first order. It is shown that the nonequilibrium part of the distribution function yields the spin Hall effect for fermions constrained to move in a plane perpendicular to the angular velocity and magnetic field. Moreover it yields an analogue of Ohm’s law for spin currents whose resistivity depends on the external magnetic field and the angular velocity of the rotating frame. Spin current densities in three-dimensional systems are also established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Damage imaging using Lamb waves for SHM applications

2-D ultrasonic arrays, due to their beam-steering capability and all azimuth angle coverage are a very promising tool for the inspection of plate-like structures using Lamb waves (LW). Contrary to the classical linear phased arrays (PAs) the 2D arrays enable unequivocal defect localization and they are even capable of mode selectivity of the received LWs . Recently, it has been shown that multistatic synthetic focusing (SF) algorithms applied for 2D arrays are much more effective than the classical phase array mode commonly used in NDT. The multistatic SF assumes multiple transmissions of elements in a transmitting aperture and off-line processing of the data acquired by a receiving aperture. In the simplest implementation of the technique, only a single multiplexed input and a number of output channels are required, which results in significant hardware simplification compared with the PA systems. On the one hand implementation of the multistatic SF to 2D arrays creates additional degrees of freedom during the design of the array topology, which complicates the array design process. On the other hand, it enables designing sparse arrays with performance similar to that of the fully populated dense arrays. In this paper we present a general systematic approach to the design and optimization of imaging systems based on the 2D array operating in the multistatic mode. We start from presenting principles of the SF schemes applied to LW imaging. Then, we outline the coarray concept and demonstrate how it can be used for reducing number of elements of the 2D arrays. Finally, efficient tools for the investigation and experimental verification of the designed 2D array prototypes are presented. The first step in the investigation is theoretical evaluation performed using frequency-dependent structure transfer function (STF), which enables approximate simulation of an array excited with a tone-burst in a dispersive medium. Finally, we show how scanning laser vibrometer, sensing waves in multiple points corresponding to the locations of the 2D receiving array elements, can be used as a tool for rapid experimental verification of the developed topologies. The presented methods are discussed in terms of the beampatterns and sparse versions of the fully populated array topologies are be presented. The effect of apodization applied to the array elements is also investigated. Both simulated and experimental results are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Cosmological angular trispectra and non-Gaussian covariance

Angular cosmological correlators are infamously difficult to compute due to the highly oscillatory nature of the projection integrals. Motivated by recent development on analytic approaches to cosmological perturbation theory, in this paper we present an efficient method for computing cosmological four-point correlations in angular space, generalizing previous works on lower-point functions. This builds on the FFTLog algorithm that approximates the matter power spectrum as a sum over power-law functions, which makes certain momentum integrals analytically solvable. The computational complexity is drastically reduced for correlators in a "separable" form—we define a suitable notion of separability for cosmological trispectra, and derive formulas for angular correlators of different separability classes. As an application of our formalism, we compute the angular galaxy trispectrum at tree level, with and without primordial non-Gaussianity. This includes effects of redshift space distortion and bias parameters up to cubic order. We also compute the non-Gaussian covariance of the angular matter power spectrum due to the connected four-point function, beyond the Limber approximation. Finally, we demonstrate that, in contrast to the standard lore, the Limber approximation can fail for the non-Gaussian covariance computation even for large multipoles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Compressed-sensing-based content-driven hierarchical reconstruction: Theory and application to C-arm cone-beam tomography

Purpose: This paper addresses the reconstruction of x-ray cone-beam computed tomography (CBCT) for interventional C-arm systems. Subsampling of CBCT is a significant issue with C-arms due to their slow rotation and to the low frame rate of their flat panel x-ray detectors. The aim of this work is to propose a novel method able to handle the subsampling artifacts generally observed with analytical reconstruction, through a content-driven hierarchical reconstruction based on compressed sensing. Methods: The central idea is to proceed with a hierarchical method where the most salient features (high intensities or gradients) are reconstructed first to reduce the artifacts these features induce. These artifacts are addressed first because their presence contaminates less salient features. Several hierarchical schemes aiming at streak artifacts reduction are introduced for C-arm CBCT: the empirical orthogonal matching pursuit approach with the ℓ{sub 0} pseudonorm for reconstructing sparse vessels; a convex variant using homotopy with the ℓ{sub 1}-norm constraint of compressed sensing, for reconstructing sparse vessels over a nonsparse background; homotopy with total variation (TV); and a novel empirical extension to nonlinear diffusion (NLD). Such principles are implemented with penalized iterative filtered backprojection algorithms. For soft-tissue imaging, the authors compare the use of TV and NLD filters as sparsity constraints, both optimized with the alternating direction method of multipliers, using a threshold for TV and a nonlinear weighting for NLD. Results: The authors show on simulated data that their approach provides fast convergence to good approximations of the solution of the TV-constrained minimization problem introduced by the compressed sensing theory. Using C-arm CBCT clinical data, the authors show that both TV and NLD can deliver improved image quality by reducing streaks. Conclusions: A flexible compressed-sensing-based algorithmic approach is proposed that is able to accommodate for a wide range of constraints. It is successfully applied to C-arm CBCT images that may not be so well approximated by piecewise constant functions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Electron capture of superheavy nuclei with realistic lepton wave functions

The superheavy nuclei push the periodic table of the elements and the chart of the nuclides to their limits, providing a unique laboratory for studies of the electron-nucleus interactions. The most important weak decay mode in known superheavy nuclei is electron capture (EC). In the standard calculations of EC, the lepton wave functions are usually considered in the lowest-order approximation. In this work, we investigate the sensitivity of EC rates on the choice of the electron wave functions by (i) assuming the single-particle approximation for the electron wave functions, and (ii) carrying out Dirac-Hartree-Fock (DHF) calculations. The nuclear response is generated based on the state-of-the-art quasiparticle random phase approximation employing relativistic nuclear energy density functional theory. Here, we show that using the improved lepton wave functions reduces the EC rates up to 40% in the superheavy nucleus oganesson (𝑍=118). Interestingly, because of screening effects, the difference between the EC rates obtained with the DHF and single-particle calculations is fairly small.

Atomic orbital↗
DOE OSTI2017

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. In this paper, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calculate solvation free energies of charged hard spheres using density functional theory interaction potentials with molecular dynamics simulation and isolate the effects of charge and cavitation, comparing to the Born (linear response) model. We show that using uncorrected Ewald summation leads to unphysical values for the single ion solvation free energy and that charging free energies for cations are approximately linear as a function of charge but that there is a small non-linearity for small anions. The charge hydration asymmetry for hard spheres, determined with quantum mechanics, is much larger than for the analogous real ions. Finally, this suggests that real ions, particularly anions, are significantly more complex than simple charged hard spheres, a commonly employed representation.

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

A new class of ensemble conserving algorithms for approximate quantum dynamics: Theoretical formulation and model problems

We develop two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in combination with the Feynman-Kleinert approximation of the density operator. These dynamics are used to improve the Feynman-Kleinert implementation of the classical Wigner approximation for the evaluation of quantum time correlation functions known as Feynman-Kleinert linearized path-integral. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics is made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that these dynamics provide a great improvement over the classical Wigner approximation, in which purely classical dynamics are used. In a special case, our first method becomes identical to centroid molecular dynamics.

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

Efficient Ensemble-Based Stochastic Gradient Methods for Optimization Under Geological Uncertainty

Ensemble-based stochastic gradient methods, such as the ensemble optimization (EnOpt) method, the simplex gradient (SG) method, and the stochastic simplex approximate gradient (StoSAG) method, approximate the gradient of an objective function using an ensemble of perturbed control vectors. These methods are increasingly used in solving reservoir optimization problems because they are not only easy to parallelize and couple with any simulator but also computationally more efficient than the conventional finite-difference method for gradient calculations. In this work, we show that EnOpt may fail to achieve sufficient improvement of the objective function when the differences between the objective function values of perturbed control variables and their ensemble mean are large. On the basis of the comparison of EnOpt and SG, we propose a hybrid gradient of EnOpt and SG to save on the computational cost of SG. We also suggest practical ways to reduce the computational cost of EnOpt and StoSAG by approximating the objective function values of unperturbed control variables using the values of perturbed ones. We first demonstrate the performance of our improved ensemble schemes using a benchmark problem. Results show that the proposed gradients saved about 30–50% of the computational cost of the same optimization by using EnOpt, SG, and StoSAG. As a real application, we consider pressure management in carbon storage reservoirs, for which brine extraction wells need to be optimally placed to reduce reservoir pressure buildup while maximizing the net present value. Results show that our improved schemes reduce the computational cost significantly.

58 GEOSCIENCES↗