Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Elastic properties of perovskite ATiO{sub 3} (A = Be, Mg, Ca, Sr, and Ba) and PbBO{sub 3} (B = Ti, Zr, and Hf): First principles calculations

The mechanical properties of perovskite oxides depend on two metal oxide lattices that are intercalated. This provides an opportunity for separate tuning of hardness, Poisson's ratio (transverse expansion in response to the compression), and shear strength. The elastic constants of series of perovskite oxides were studied by first principles approach. Both A-site and B-site cations were systematically varied in order to see their effects on the elastic parameters. To study the effects of A-site cations, we studied the elastic properties of perovskite ATiO{sub 3} for A being Be, Mg, Ca, Sr, or Ba, one at a time. Similarly, for B-site cations, we studied the elastic properties of PbBO{sub 3} for B being Ti, Zr, or Hf, one at a time. The density functional first principles calculations with local density approximation (LDA) and generalized gradient approximation (GGA) were employed. It is found that the maximum C{sub 11} elastic constant is achieved when the atomic size of the cations at A-site and B-site are comparable. We also found that C{sub 12} elastic constant is sensitive to B-site cations while C{sub 44} elastic constant is more sensitive to A-site cations. Details and explanations for such dependencies are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

ROM-based Subset Selection Algorithm for Efficient Surrogate Modeling

This summary discusses a new methodology to support the construction of surrogate models for computationally expensive models with many input parameters and output responses. In realistic engineering calculations, such as thermal-hydraulics simulation, the cost of the simulation is typically expensive, especially when detailed models are employed to describe system behavior, thereby justifying the need for surrogate models. The cost of surrogate model construction is typically a function of the number of model parameters and the degree of model nonlinearity. For sufficiently complex models, the cost of constructing the surrogate model could be by itself overwhelmingly expensive, thereby diminishing its value. One key requirement of surrogate model construction is the ability to preselect the best parametric form to describe the model behavior over the expected range of parameter variations. This summary presents a new algorithm to help select the best approximation out of a template of functions using recent advances in reduced order modeling (ROM) techniques. The idea of selecting functions out of a template is typically referred to as subset selection in the statistical and applied mathematics community. Different from existing subset selection algorithms, we rely on ROM to help identify the most influential functions for the surrogate model. Previous work has shown that ROM can help reduce the effective dimensionality of the parameter space, which can help reduce the cost of surrogate model construction. In this summary, we adopt a different philosophy, where the fitting functions are treated as pseudo parameters and ROM is applied on the pseudo parameter space, with the result now being a ranked list of the pseudo parameters along with their contribution to the quality of the surrogate model fit. This approach allows the analyst to discard all terms that have negligible impact on the quantity of interest. A relap5 model is employed to demonstrate the application of the proposed methodology. (authors)

42 ENGINEERING↗
DOE OSTI2019

Nuclear properties for astrophysical and radioactive-ion-beam applications (II)

We tabulate the ground-state odd-proton and odd-neutron spins and parities, proton and neutron pairing gaps, one- and two-neutron separation energies, quantities related to β-delayed one- and two-neutron emission probabilities, average energy and average number of emitted neutrons, β-decay energy release and half-life with respect to Gamow–Teller decay with a phenomenological treatment of first-forbidden decays, one- and two-proton separation energies, and α-decay energy release and half-life for 9318 nuclei ranging from {sup 16}O to {sup 339}136 and extending from the proton drip line to the neutron drip line. This paper is a new and improved version of Atomic Data And Nuclear Data Tables [66 131 (1997)]. The starting point of our present work is the new study (FRDM(2012)) of nuclear ground-state masses and deformations based on the finite-range droplet model and folded-Yukawa single-particle potential published in a previous issue of Atomic Data And Nuclear Data Tables [109–110, 1 (2016)]. The β-delayed neutron-emission probabilities and Gamow–Teller β-decay rates are obtained from a quasi-particle random-phase approximation with single-particle levels and wave functions at the calculated nuclear ground-state shapes as input quantities. A development since 1997 is we now use a Hauser–Feshbach approach to account for (n, γ) competition and treat first-forbidden decay in a phenomenological approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

A new modified BFGS method for unconstrained optimization problems

Using chain rule, we propose a modified secant equation to get a more accurate approximation of the second curvature of the objective function. Then, based on this modified secant equation we present a new BFGS method for solving unconstrained optimization problems. The proposed method makes use of both gradient and function values, and utilizes information from two most recent steps, while the usual secant relation uses only the latest step information. Under appropriate conditions, we show that the proposed method is globally convergent without convexity assumption on the objective function. Comparative numerical results show computational efficiency of the proposed method in the sense of the Dolan–Moré performance profiles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Prediction of magnetic moment collapse in ZrFe{sub 2} under hydrostatic pressure

Electronic structure and magnetic properties of ZrFe{sub 2} in the cubic Laves phase are investigated by calculations based on density functional theory. The magnetic moment decreases with the increase of the hydrostatic pressure in an unusual way: Two-step magnetic collapse is predicted. The first one is a continuous change from 1.53 μ{sub B}/Fe to 0.63 μ{sub B}/Fe at about 3.6 GPa, and the other is from 0.25 μ{sub B}/Fe to the nonmagnetic state at about 15 GPa in a first order manner under the local spin density approximation of the exchange correlation potential. A metastable state with intermediate spin moment about 0.15 μ{sub B}/Fe may exist before that. We understand this process by the changes of density of states during it. The magnetic moment decreases under the pressure in the vicinity of the experimental lattice constant with dlnm/dp=−0.038 GPa{sup −1}. The spontaneous volume magnetostriction is 3.6%, which is huge enough to find potential applications in magnetostriction actuators and sensors. We suggest that the Invar effect of this compound may be understood when considering the magnetic moment variation according to the magnetostrictive model of Invar.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Efficient Computation of Info-Gap Robustness for Finite Element Models

A recent research effort at LANL proposed info-gap decision theory as a framework by which to measure the predictive maturity of numerical models. Info-gap theory explores the trade-offs between accuracy, that is, the extent to which predictions reproduce the physical measurements, and robustness, that is, the extent to which predictions are insensitive to modeling assumptions. Both accuracy and robustness are necessary to demonstrate predictive maturity. However, conducting an info-gap analysis can present a formidable challenge, from the standpoint of the required computational resources. This is because a robustness function requires the resolution of multiple optimization problems. This report offers an alternative, adjoint methodology to assess the info-gap robustness of Ax = b-like numerical models solved for a solution x. Two situations that can arise in structural analysis and design are briefly described and contextualized within the info-gap decision theory framework. The treatments of the info-gap problems, using the adjoint methodology are outlined in detail, and the latter problem is solved for four separate finite element models. As compared to statistical sampling, the proposed methodology offers highly accurate approximations of info-gap robustness functions for the finite element models considered in the report, at a small fraction of the computational cost. It is noted that this report considers only linear systems; a natural follow-on study would extend the methodologies described herein to include nonlinear systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2011

Generation of matter-wave solitons of the Gross-Pitaevskii equation with a time-dependent complicated potential

We examine the generation of bright matter-wave solitons in the Gross-Pitaevskii equation describing Bose-Einstein condensates with a time-dependent complex potential, which is composed of a repulsive parabolic background potential and a gravitational field. By performing a modified lens-type transformation, an explicit expression for the growth rate of a purely growing modulational instability is presented and analyzed. We point out the effects of the gravitational field, as well as of the parameter related to the feeding or loss of atoms in the condensate, on the instability growth rate. It is evident from numerical simulations that the feeding with atoms and the magnetic trap have opposite effects on the dynamics of the system. It is shown that the feeding or loss parameter can be well used to control the instability domain. Our study shows that the gravitational field changes the condensate trail of the soliton trains during the propagation. We also perform a numerical analysis to solve the Gross-Pitaevskii equation with a time-dependent complicated potential. The numerical results on the effect of both the gravitational field and the parameter of feeding or loss of atoms in the condensate agree well with predictions of the linear stability analysis. Another result of the present work is the modification of the background wave function in the Thomas-Fermi approximation during the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Time-Dependent Density Functional Theory and the Real-Time Dynamics of Fermi Superfluids

Here, I describe the time-dependent superfluid local density approximation, which is an adiabatic extension of the density functional theory to superfluid Fermi systems and their real-time dynamics. This new theoretical framework has been used to report on several phenomena in cold atomic gases and nuclear collective motion: excitation of the Higgs modes in strongly interacting Fermi superfluids, generation of quantized vortices, crossing and reconnection of vortices, excitation of the superflow at velocities above the critical velocity, excitation of quantum shock waves and domain walls in the collisions of superfluid atomic clouds, and excitation of collective states in nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Spectral Analysis of 3-(Adamantan-1-yl)-4-Ethyl-1-[(4-Phenylpiperazin-1-yl) Methyl]-1H-1,2,4-Triazole-5(4H)-Thione

Vibrational IR (3200–650 cm{sup –1}) and Raman spectra (3200–150 cm{sup –1}) of adamantane-containing 3-(adamantan-1-yl)-4-ethyl-1-[(4-phenylpiperazin-1-yl)methyl]-1H-1,2,4-triazole-5(4H)-thione, which is promising for drug design, were examined. The UV/Vis spectrum (450–200 nm) of the compound in EtOH was measured. Full geometry optimization using density functional theory (DFT) in the B3LYP/cc-pVDZ approximation allowed the equilibrium configuration of the molecule to be determined and IR and Raman spectra to be calculated. Based on these, the experimental vibrational IR and Raman spectra were interpreted and the biological activity indices were predicted. The UV/Vis spectrum of the title compound was simulated at the time-dependent DFT/CAM-B3LYP/cc-pVDZ level with and without solvent effects and at the ab initio multi-reference perturbation theory XMCQDPT2 level. The UV/Vis spectrum that was simulated using the multi-reference XMCQDPT2 approximation agreed very successfully with the experimental data, in contrast to the single-reference DFT method. This was probably a consequence of intramolecular charge transfer.

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

Investigation of the deuteron spin structure at short nucleon-nucleon distances in the reaction of polarized-deuteron fragmentation to cumulative pions

Experimental results on the vector (A{sub y}) and tensor (A{sub yy}) analyzing powers in the fragmentation of 5- and 9-GeV/c polarized deuterons to high-momentum pions in the kinematical region corresponding to pion production on a strongly correlated nucleon pair (cumulative meson production) are presented. The angular and momentum dependences of A{sub yy} are not described by calculations performed in the impulse approximation by using standard deuteron wave functions. An explanation for our data should be sought on the basis of models that treat the deuteron at short distances (deuteron-core region) as a multiquark state-for example, a 6q cluster, whose high orbital angular momentum (D wave) leads to the observed strong dependence of the reaction tensor analyzing power A(d-vector, {pi})X on the pion transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Calculation of Young’s Modulus of MoS{sub 2}-Based Single-Wall Nanotubes Using Force-Field and Hybrid Density Functional Theory

A force field is proposed that reproduces with a high accuracy a large number of properties of the bulk crystal MoS{sub 2} phases, monolayers, and nanotubes. The reproduced values are both the experimental results and the results of quantum chemical calculations. The elaborated interaction potential can be useful primarily for investigation of multiwall MoS{sub 2} nanotubes and their thermodynamic properties, especially, since the potential is able to reproduce the frequencies of the crystal phonon spectrum. In this study the proposed potential is applied to simulate the temperature dependence of a number of properties of the armchair and zigzag nanotubes. The calculations have been performed using molecular mechanics method within the framework of quasi harmonic approximation which is carried out through the estimation of the temperature dependence of the Helmholtz free energy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Minimal multi-element stochastic collocation for uncertainty quantification of discontinuous functions

We propose a multi-element stochastic collocation method that can be applied in high-dimensional parameter space for functions with discontinuities lying along manifolds of general geometries. The key feature of the method is that the parameter space is decomposed into multiple elements defined by the discontinuities and thus only the minimal number of elements are utilized. On each of the resulting elements the function is smooth and can be approximated using high-order methods with fast convergence properties. The decomposition strategy is in direct contrast to the traditional multi-element approaches which define the sub-domains by repeated splitting of the axes in the parameter space. Such methods are more prone to the curse-of-dimensionality because of the fast growth of the number of elements caused by the axis based splitting. The present method is a two-step approach. Firstly a discontinuity detector is used to partition parameter space into disjoint elements in each of which the function is smooth. The detector uses an efficient combination of the high-order polynomial annihilation technique along with adaptive sparse grids, and this allows resolution of general discontinuities with a smaller number of points when the discontinuity manifold is low-dimensional. After partitioning, an adaptive technique based on the least orthogonal interpolant is used to construct a generalized Polynomial Chaos surrogate on each element. The adaptive technique reuses all information from the partitioning and is variance-suppressing. We present numerous numerical examples that illustrate the accuracy, efficiency, and generality of the method. When compared against standard locally-adaptive sparse grid methods, the present method uses many fewer number of collocation samples and is more accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Design and analysis of machine learning exchange-correlation functionals via rotationally invariant convolutional descriptors

Here, we explore the potential of a data-driven approach to the design of exchange-correlation (xc) functionals. The approach, inspired by convolutional filters in computer vision and surrogate functions from optimization, utilizes convolutions of the electron density to form a feature space to represent local electronic environments and neural networks to map the features to the exchange-correlation energy density. These features are orbital free, and provide a systematic route to including information at various length scales. This work shows that convolutional descriptors are theoretically capable of an exact representation of the electron density, and proposes Maxwell-Cartesian spherical harmonic kernels as a class of rotationally invariant descriptors for the construction of machine learned functionals. The approach is demonstrated using data from the B3LYP functional on a number of small molecules containing C, H, O, and N along with a neural network regression model. The machine learned functionals are compared to standard physical approximations and the accuracy is assessed for the absolute energy of each molecular system as well as formation energies. The results indicate that it is possible to reproduce the exchange-correlation portion of B3LYP formation energies to within chemical accuracy using orbital-free descriptors with a spatial extent of 0.2 Å. The findings provide empirical insight into the spatial range of electron exchange, and suggest that the combination of convolutional descriptors and machine learning regression models is a promising framework for xc functional design, although challenges remain in obtaining training data and generating models consistent with pseudopotentials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Fixed points and flow analysis on off-equilibrium dynamics in the boson Boltzmann equation

We consider fixed points of steady solutions and flow directions using the boson Boltzmann equation that is a one-dimensionally reduced kinetic equation after the angular integration. With an elastic collision integral of the two-to-two scattering process, in the dense (dilute) regime where the distribution function is large (small), the boson Boltzmann equation has approximate fixed points with a power-law spectrum in addition to the thermal distribution function. We argue that the power-law fixed point can be exact in special cases. We elaborate a graphical presentation to display evolving flow directions similarly to the renormalization group flow, which explicitly exhibits how fixed points are connected and parameter space is separated by critical lines. We discuss that such a flow diagram contains useful information on thermalization processes out of equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

First-principles study of structural stability, electronic, optical and elastic properties of binary intermetallic: PtZr

Structural, electronic, optical and elastic properties of PtZr have been studied using the full-potential linearized augmented plane wave (FP-LAPW) method within density functional theory (DFT). The energy against volume and enthalpy vs. pressure variation in three different structures i.e. B{sub 1}, B{sub 2} and B{sub 3} for PtZr has been presented. The equilibrium lattice parameter, bulk modulus and its pressure derivative have been obtained using optimization method for all the three phases. Furthermore, electronic structure was discussed to reveal the metallic character of the present compound. The linear optical properties are also studied under zero pressure for the first time. Results on elastic properties are obtained using generalized gradient approximation (GGA) for exchange correlation potentials. Ductile nature of PtZr compound is predicted in accordance with Pugh’s criteria.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗