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

Efficient simulations of charge density waves in the transition metal Dichalcogenide TiSe 2

Charge density waves (CDWs) in transition metal dichalcogenides are the subject of growing scientific interest due to their rich interplay with exotic phases of matter and their potential technological applications. Here, using density functional theory with advanced meta-generalized gradient approximations (meta-GGAs) and linear response time-dependent density functional theory (TDDFT) with state-of-the-art exchange-correlation kernels, we investigate the electronic, vibrational, and optical properties in 1T-TiSe 2 with and without CDW. In both bulk and monolayer TiSe 2 , the electronic bands and phonon dispersions in either normal or CDW (semiconducting) phase are described well via meta-GGAs, which separate the valence and conduction bands just as HSE06 does but with significantly more computational feasibility. The experimentally observed humps of electron energy loss spectroscopy are successfully reproduced in TDDFT. Our work opens the door to simulating these complexities in CDW compounds from first principles by revealing meta-GGAs as an accurate low-cost alternative to HSE06.

Electronic properties and materials↗
DOE OSTI2021

Spatial correlations of entangled polymer dynamics

The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Furthermore, comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2004

Progressive Response Surfaces

Response surface functions are often used as simple and inexpensive replacements for computationally expensive computer models that simulate the behavior of a complex system over some parameter space. Progressive response surfaces are ones that are built up progressively as global information is added from new sample points in the parameter space. As the response surfaces are globally upgraded based on new information, heuristic indications of the convergence of the response surface approximation to the exact (fitted) function can be inferred. Sampling points can be incrementally added in a structured fashion, or in an unstructured fashion. Whatever the approach, at least in early stages of sampling it is usually desirable to sample the entire parameter space uniformly. At later stages of sampling, depending on the nature of the quantity being resolved, it may be desirable to continue sampling uniformly over the entire parameter space (Progressive response surfaces), or to switch to a focusing/economizing strategy of preferentially sampling certain regions of the parameter space based on information gained in early stages of sampling (Adaptive response surfaces). Here we consider Progressive response surfaces where a balanced indication of global response over the parameter space is desired.We use a variant of Moving Least Squares to fit and interpolate structured and unstructured point sets over the parameter space. On a 2-D test problem we compare response surface accuracy for three incremental sampling methods: Progressive Lattice Sampling; Simple-Random Monte Carlo; and Halton Quasi-Monte-Carlo sequences. We are ultimately after a system for constructing efficiently upgradable response surface approximations with reliable error estimates.

Romero, V. J.↗
DOE OSTI2016

Functional renormalization group study of nuclear and neutron matter

A chiral model based on nucleons interacting via boson exchange is investigated. Fluctuation effects are included consistently beyond the mean-field approximation in the framework of the functional renormalization group. The liquid-gas phase transition of symmetric nuclear matter is studied in detail. No sign of a chiral restoration transition is found up to temperatures of about 100 MeV and densities of at least three times the density of normal nuclear matter. Moreover, the model is extended to asymmetric nuclear matter and the constraints from neutron star observations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Statistics of time delay and scattering correlation functions in chaotic systems. I. Random matrix theory

We consider the statistics of time delay in a chaotic cavity having M open channels, in the absence of time-reversal invariance. In the random matrix theory approach, we compute the average value of polynomial functions of the time delay matrix Q = − iħS{sup †}dS/dE, where S is the scattering matrix. Our results do not assume M to be large. In a companion paper, we develop a semiclassical approximation to S-matrix correlation functions, from which the statistics of Q can also be derived. Together, these papers contribute to establishing the conjectured equivalence between the random matrix and the semiclassical approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

First-order derivative couplings between excited states from adiabatic TDDFT response theory

We present a complete derivation of derivative couplings between excited states in the framework of adiabatic time-dependent density functional response theory. Explicit working equations are given and the resulting derivative couplings are compared with derivative couplings from a pseudo-wavefunction ansatz. For degenerate excited states, i.e., close to a conical intersection (CI), the two approaches are identical apart from an antisymmetric overlap term. However, if the difference between two excitation energies equals another excitation energy, the couplings from response theory exhibit an unphysical divergence. This spurious behavior is a result of the adiabatic or static kernel approximation of time-dependent density functional theory leading to an incorrect analytical structure of the quadratic response function. Numerical examples for couplings close to a CI and for well-separated electronic states are given.

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

First-Principles Investigation of Magnetic Properties and Faraday Rotation of Co Doping CdTe

In this paper, we aimed to study the effect of doping of the compound CdTe with the cobalt impurity (Co), as well as the vacancy defects in Cd sites. On one hand, this leads to the investigation of the magnetic properties and the Faraday rotation effect for the studied alloy, doped with different concentrations of cobalt (0.01, 0.05, 0.1, 0.15, 0.2, 0.25). On the other hand, we have created 0.01 of vacancy defects in Cd sites. Then, we raised the vacancy defect concentration to 0.05, keeping the same concentrations of cobalt. As a result, we have investigated that there is a magnetism appearing with the cobalt doping, while the vacancy defects in Cd sites affect the stability of the magnetic states. In general, it improves the ferromagnetic state that will be well explained in the discussion. Our calculations were performed using the KKR-CPA method within the spin-polarized density functional theory (DFT) with the local density approximation (LDA). The curves of the density of states (DOS) illustrate the results of this study which has been discussed, analyzed, and explained below. In addition, the energy of each case was calculated and given in the tables below for both the DLM (half of the cobalt spins are up while the other half are down) and the ferromagnetic state in order to confirm which one of them is stable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Structural and electronic properties of TiX (X=N, As) in rock salt and zinc blende phase: A DFT study

Quantum mechanical first principle calculations have been performed to study the electronic and structural properties of TiN and TiAs in zinc blende (ZB) and rock salt (RS) structures. The full-potential linearized augmented plane wave (FP-LAPW) method has been used within the framework of density functional theory (DFT). The exchange correlation functional has been solved employing generalized gradient approximation (GGA). Our predicted results for lattice constants are in good agreement with the earlier findings. The electronic band structures of TiX are metallic in both the phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Ab initio theory of the equation of state for compressed rare gas crystals

Nonempirical equations of state of compressed rare gas crystals Ne, Ar, Kr, and Xe are studied on the basis of the earlier-obtained ab initio adiabatic potential. The paired and three-body short-range repulsive potentials are calculated by the Hartree–Fock method in the basis of localized functions with their exact mutual orthogonalization and do not contain experimentally determined parameters. The theory is compared with the experiment and results of calculations by other authors. Analysis of the proposed equations of state for large compressions has shown the importance of taking into account the three-body interaction and the terms of the higher order in the overlap integral in compressed neon and the sufficiency of the quadratic approximation in the orthogonalization of functions in heavy rare gas crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2013

On Using Intensity Interferometry for Feature Identification and Imaging of Remote Objects

We derive an approximation to the intensity covariance function of two scanning pinhole detectors, facing a distant source (e.g., a star) being occluded partially by an absorptive object (e.g., a planet). We focus on using this technique to identify or image an object that is in the line-of-sight between a well-characterized source and the detectors. We derive the observed perturbation to the intensity covariance map due to the object, showing that under some reasonable approximations it is proportional to the real part of the Fourier transform of the source's photon-flux density times the Fourier transform of the object's intensity absorption. We highlight the key parameters impacting its visibility and discuss the requirements for estimating object-related parameters, e.g., its size, velocity or shape. We consider an application of this result to determining the orbit inclination of an exoplanet orbiting a distant star. Finally, motivated by the intrinsically weak nature of the signature, we study its signal-to-noise ratio and determine the impact of system parameters.

quantum optics↗
NASA NTRS1998

A Study of Transport in the Near-Earth Plasma Sheet During A Substorm Using Time-Dependent Large Scale Kinetics

In this study we investigate the transport of H+ ions that made up the complex ion distribution function observed by the Geotail spacecraft at 0740 UT on November 24, 1996. This ion distribution function, observed by Geotail at approximately 20 R(sub E) downtail, was used to initialize a time-dependent large-scale kinetic (LSK) calculation of the trajectories of 75,000 ions forward in time. Time-dependent magnetic and electric fields were obtained from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interaction with the solar wind and the interplanetary magnetic field (IMF) as observed during the interval of the observation of the distribution function. Our calculations indicate that the particles observed by Geotail were scattered across the equatorial plane by the multiple interactions with the current sheet and then convected sunward. They were energized by the dawn-dusk electric field during their transport from Geotail location and ultimately were lost at the ionospheric boundary or into the magnetopause.

El-Alaoui, M.↗
DOE OSTI2014

Electronic and optical properties of NbO{sub 2}

In the present study, we combine theoretical and experimental approaches in order to gain insight into the electronic properties of both the high-temperature, rutile (metallic) and low-temperature, body-centered tetragonal (insulating) phase of niobium dioxide (NbO{sub 2}) as well as the optical properties of the low-temperature phase. Theoretical calculations performed at the level of the local density approximation, Hubbard U correction, and hybrid functional are complemented with the spectroscopic ellipsometry (SE) of epitaxial films grown by molecular beam epitaxy. For the rutile phase, the local density approximation (LDA) gives the best description and predicts Fermi surface nesting consistent with wave vectors that lead to niobium-niobium dimerization during the phase transition. For the insulating phase, LDA provides a good quantitative description of the lattice, but only a qualitative description for the band gap. Including a Hubbard U correction opens the band gap at the expense of correctly describing the valence band and lattice of both phases. The hybrid functional slightly overestimates the band gap. Ellipsometric measurement is consistent with insulating behavior with a 1.0 eV band gap. Comparison with the theoretical dielectric functions, obtained utilizing a scissors operator to adjust the LDA band gap to reproduce the ellipsometry data, allows for identification of the SE peak features.

36 MATERIALS SCIENCE↗
NASA NTRS2003

On the Applicability of High-frequency Approximations to Lilley's Equation

Three forms of the high-frequency asymptotic Green's function for Lilley's equation are reviewed and compared to the exact solution over wide range of Strouhal numbers. The asymmetric approximation, which applies to sources away form the jet axis, and the quasi-symmetric approximation, which is arrived at by making a near-axis source assumption, are both obtained for parallel round jets from a formal Fourier-transform solution. The ray-theory solution, which is the only high-frequency approximation that can be applied to more general mean flows, follows from a WKB ansatz and is shown to be closely related to the asymmetric approximation. The comparisons show that the best overall prediction of the exact Green's function is given by the asymmetric approximation which remains accurate down to a Strouhal number of 1/2. The close relationship between the asymmetric and ray-theory approximations suggests that the high-frequency asymptotic Green's function for more general mean flows would be similarly successful.

Wundrow, David W.↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1943

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗