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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Augmented Adam-Gibbs model for glass melt viscosity and configuration entropy as functions of temperature and composition

As the temperature of glass melt increases, its structure approaches the state of a simple liquid while the configuration entropy approaches a maximum value. We describe this gradual change using a power law function of inverse temperature. The Adam-Gibbs model for glass viscosity as a function of temperature and glass composition augmented in this way is greatly simplified when applied to common glass families occupying moderate composition regions, such as float glass or nuclear waste glasses, on which properties can be approximated as linear functions of composition. The parsimonious model thus obtained is preferable for use in optimizing glass formulation and mathematical modeling of glass melting and forming. For multicomponent glasses with N viscosity-affecting components, the augmented Adam-Gibbs model requires 2N + 3 adjustable parameters. The model efficacy is demonstrated by fitting the model to a viscosity-temperature-composition dataset for low-activity nuclear waste glasses.

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

Bootstrap-based confidence interval estimation for thermal security region of bulk power grid

Thermal security region (TSR) is a powerful tool for monitoring and controlling the thermal security of bulk power grid with high penetration of renewable energy. One of the major challenges of TSR in practical application is to obtain the exact analytical expression of the TSR boundary (TSRB) because TSRB is a nonlinear hypersurface described by an implicit function. For this reason, the TSRB is usually approximated by hyperplane instead of getting the exact analytical expression. Traditionally, the coefficients of the hyperplane are estimated with a set of original TSRB points using least square estimation (LSE). However, LSE is a point estimation method, which is incapable of evaluating the quality of the estimation. In addition, the accuracy and computational efficiency of the hyperplane approximation are influenced by the number of TSRB points significantly. Here, a bootstrap based confidence interval estimation is proposed to estimate not only the coefficients of TSRB approximation hyperplane, but also the standard deviations and confidence intervals of the coefficients for evaluating the quality and reliability of the approximation results. First, empirical distribution functions (EDFs) of the coefficients of TSRB approximation hyperplane are approximated from a set of original TSRB points by using residual resampling bootstrap method. Then, the EDFs of the coefficients are employed to estimate the coefficients of the approximation hyperplane. Meanwhile, the standard deviations and confidence intervals of the estimated coefficients of the hyperplanes are also calculated for evaluating the quality and reliability of the approximation. The proposed approach is tested on the China Southern Power Grid (CSG). Results of simulations validate the accuracy and efficiency of the proposed approach in approximating the TSRB.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol -1 relative to the experiment. Here, the two-point QZ/5Z and TZ/QZ CBS 1/$l^3_{max}$ extrapolation schemes are inferior to their inverse power counterpart (1/($l_{max}$ + 1/2) 4 ) by 1.3 and 4.3 kJ mol -1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol -1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double- ζ/triple- ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

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

Blanco DECam Bulge Survey (BDBS) II: project performance, data analysis, and early science results

ABSTRACT The Blanco DECam Bulge Survey (BDBS) imaged more than 200 sq deg of the Southern Galactic bulge using the ugrizY filters of the Dark Energy Camera, and produced point spread function photometry of approximately 250 million unique sources. In this paper, we present details regarding the construction and collation of survey catalogues, and also discuss the adopted calibration and dereddening procedures. Early science results are presented with a particular emphasis on the bulge metallicity distribution function and globular clusters. A key result is the strong correlation (σ ∼ 0.2 dex) between (u − i)o and [Fe/H] for bulge red clump giants. We utilized this relation to find that interior bulge fields may be well described by simple closed box enrichment models, but fields exterior to b ∼ −6° seem to require a secondary metal-poor component. Applying scaled versions of the closed box model to the outer bulge fields is shown to significantly reduce the strengths of any additional metal-poor components when compared to Gaussian mixture models. Additional results include: a confirmation that the u band splits the subgiant branch in M22 as a function of metallicity, the detection of possible extratidal stars along the orbits of M 22 and FSR 1758, and additional evidence that NGC 6569 may have a small but discrete He spread, as evidenced by red clump luminosity variations in the reddest bands. We do not confirm previous claims that FSR 1758 is part of a larger extended structure.

Johnson, Christian I.↗
DOE OSTI2024

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Efficient prediction of time- and angle-resolved photoemission spectroscopy measurements on a nonequilibrium BCS superconductor

We study how time- and angle-resolved photoemission (tr-ARPES) reveals the dynamics of BCS-type, s -wave superconducting systems with time-varying order parameters. Approximate methods are discussed, based on previous approaches to either optical conductivity or quantum dot transport, to enable computationally efficient prediction of photoemission spectra. One use of such predictions is to enable extraction of the underlying order parameter dynamics from experimental data, which is topical given the rapidly growing use of tr-ARPES in studying unconventional superconductivity. The methods considered model the two-time lesser Green's functions with an approximated lesser self-energy that describes relaxation by coupling of the system to two types of baths. The approach primarily used here also takes into consideration the relaxation of the excited states into equilibrium by explicitly including the level-broadening of the retarded and advanced Green's functions. Here, we present equilibrium and nonequilibrium calculations of tr-ARPES spectrum from our model and discuss the signatures of different types of superconducting dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

First-principles calculations of the electronic and structural properties of GaSb

In this paper, we carried out first-principles calculations in order to investigate the structural and electronic properties of the binary compound gallium antimonide (GaSb). This theoretical study was carried out using the Density Functional Theory within the plane-wave pseudopotential method. The effects of exchange and correlation (XC) were treated using the functional Local Density Approximation (LDA), generalized gradient approximation (GGA): Perdew–Burke–Ernzerhof (PBE), Perdew-Burke-Ernzerhof revised for solids (PBEsol), Perdew-Wang91 (PW91), revised Perdew–Burke–Ernzerhof (rPBE), Armiento–Mattson 2005 (AM05) and meta-generalized gradient approximation (meta-GGA): Tao–Perdew–Staroverov–Scuseria (TPSS) and revised Tao–Perdew–Staroverov–Scuseria (RTPSS) and modified Becke-Johnson (MBJ). We calculated the densities of state (DOS) and band structure with different XC potentials identified and compared them with the theoretical and experimental results reported in the literature. It was discovered that functional: LDA, PBEsol, AM05 and RTPSS provide the best results to calculate the lattice parameters (a) and bulk modulus (B{sub 0}); while for the cohesive energy (E{sub coh}), functional: AM05, RTPSS and PW91 are closer to the values obtained experimentally. The MBJ, Rtpss and AM05 values found for the band gap energy is slightly underestimated with those values reported experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A numerical method for deriving shape functions of nanoparticles for pair distribution function refinements

We report that in the structural refinement of nanoparticles, discrete atomistic modeling can be used for small nanocrystals (< 15 nm), but becomes computationally unfeasible at larger sizes, where instead unit-cell-based small-box modeling is usually employed. However, the effect of the nanocrystal's shape is often ignored or accounted for with a spherical model regardless of the actual shape due to the complexities of solving and implementing accurate shape effects. Recent advancements have provided a way to determine the shape function directly from a pair distribution function calculated from a discrete atomistic model of any given shape, including both regular polyhedra (e.g. cubes, spheres, octahedra) and anisotropic shapes (e.g. rods, discs, ellipsoids) [Olds et al. (2015). J. Appl. Cryst. 48, 1651–1659], although this approach is still limited to small size regimes due to computational demands. In order to accurately account for the effects of nanoparticle size and shape in small-box refinements, a numerical or analytical description is needed. This article presents a methodology to derive numerical approximations of nanoparticle shape functions by fitting to a training set of known shape functions; the numerical approximations can then be employed on larger sizes yielding a more accurate and physically meaningful refined nanoparticle size. Lastly, the method is demonstrated on a series of simulated and real data sets, and a table of pre-calculated shape function expressions for a selection of common shapes is provided.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Chaos control of the micro-electro-mechanical resonator by using adaptive dynamic surface technology with extended state observer

This paper addresses chaos control of the micro-electro- mechanical resonator by using adaptive dynamic surface technology with extended state observer. To reveal the mechanism of the micro- electro-mechanical resonator, the phase diagrams and corresponding time histories are given to research the nonlinear dynamics and chaotic behavior, and Homoclinic and heteroclinic chaos which relate closely with the appearance of chaos are presented based on the potential function. To eliminate the effect of chaos, an adaptive dynamic surface control scheme with extended state observer is designed to convert random motion into regular motion without precise system model parameters and measured variables. Putting tracking differentiator into chaos controller solves the ‘explosion of complexity’ of backstepping and poor precision of the first-order filters. Meanwhile, to obtain high performance, a neural network with adaptive law is employed to approximate unknown nonlinear function in the process of controller design. The boundedness of all the signals of the closed-loop system is proved in theoretical analysis. Finally, numerical simulations are executed and extensive results illustrate effectiveness and robustness of the proposed scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact of Surface Chemistry of Ultrasmall Superparamagnetic Iron Oxide Nanoparticles on Protein Corona Formation and Endothelial Cell Uptake, Toxicity, and Barrier Function

Ultrasmall superparamagnetic iron oxide nanoparticles (USPIONs) have been investigated for biomedical applications, including novel contrast agents, magnetic tracers for tumor imaging, targeted drug delivery vehicles, and magneto-mechanical actuators for hyperthermia and thrombolysis. Despite significant progress, recent clinical reports have raised concerns regarding USPION safety related to endothelial cell dysfunction; however, there is limited information on factors contributing to these clinical responses. The influence of USPION surface chemistry on nanoparticle interactions with proteins may impact endothelial cell function leading to adverse responses. Therefore, the goal of this study was to assess the effects of carboxyl-functionalized USPION (CU) or amine-functionalized USPION (AU) (approximately 30 nm diameter) on biological responses in human coronary artery endothelial cells. Increased protein adsorption was observed for AU compared with CU after exposure to serum proteins. Exposure to CU, but not AU, resulted in a concentration-dependent decrease in cell viability and perinuclear accumulation inside cytoplasmic vesicles. Internalization of CU was correlated with endothelial cell functional changes under non-cytotoxic conditions, as evidenced by a marked decreased expression of endothelial-specific adhesion proteins (eg, vascular endothelial-cadherin and platelet endothelial cell adhesion molecule-1) and increased endothelial permeability. Evaluation of downstream signaling indicated endothelial permeability is associated with actin cytoskeleton remodeling, possibly elicited by intracellular events involving reactive oxygen species, calcium ions, and the nanoparticle cellular uptake pathway. This study demonstrated that USPION surface chemistry significantly impacts protein adsorption and endothelial cell uptake, viability, and barrier function. This information will advance the current toxicological profile of USPION and improve development, safety assessment, and clinical outcomes of USPION-enabled medical products.

Toxicology↗
DOE OSTI2017

A degradation function consistent with Cocks–Ashby porosity kinetics

Here, the load carrying capacity of ductile materials degrades as a function of porosity, stress state and strain-rate. The effect of these variables on porosity kinetics is captured by the Cocks–Ashby model; however, the Cocks–Ashby model does not account for material degradation directly. This work uses a yield criteria to form a degradation function that is consistent with Cocks–Ashby porosity kinetics and is a function of porosity, stress state and strain-rate dependence. Approximations of this degradation function for pure hydrostatic stress states are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ground state wave functions for single-band Hubbard models from the Gutzwiller conjugate gradient minimisation theory

The Gutzwiller conjugate gradient minimisation (GCGM) theory is an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. Previous applications of GCGM provides satisfying accuracy of ground-state energy of Hubbard models. In the current work, we address the problem of whether the obtained wave function is a good approximation for the true ground state by comparing the correlation functions with the benchmark data. Additionally, our results confirms the accuracy of the reproduced ground state of the regular Hubbard model, but with some exception for the frustrated Hubbard model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Description of quasiparticle and satellite properties via cumulant expansions of the retarded one-particle Green's function

A cumulant-based GW approximation for the retarded one-particle Green's function is proposed, motivated by an exact relation between the improper Dyson self-energy and the cumulant generating function. We explore qualitative aspects of this method within a simple one-electron independent phonon model, where it is seen that the method preserves the energy moment of the spectral weight while also reproducing the exact Green's function in the weak-coupling limit. For the three-dimensional electron gas, this method predicts multiple satellites at the bottom of the band, albeit with inaccurate peak spacing. But, its quasiparticle properties and correlation energies are more accurate than both previous cumulant methods and standard G0W0. These results point to features that may be exploited within the framework of cumulant-based methods and suggest promising directions for future exploration and improvements of cumulant-based GW approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Diffusion quantum Monte Carlo calculations of SrFeO 3 and LaFeO 3

The equations of state, formation energy, and migration energy barrier of the oxygen vacancy in SrFeO 3 and LaFeO 3 were calculated in this paper with the diffusion quantum Monte Carlo (DMC) method. Calculations were also performed with various Density Functional Theory (DFT) approximations for comparison. DMC reproduces the measured cohesive energies of these materials with errors below 0.23(5) eV and the structural properties within 1% of the experimental values. The DMC formation energies of the oxygen vacancy in SrFeO 3 and LaFeO 3 under oxygen-rich conditions are 1.3(1) and 6.24(7) eV, respectively. Similar calculations with semi-local DFT approximations for LaFeO 3 yielded vacancy formation energies 1.5 eV lower. Comparison of charge density evaluated with DMC and DFT approximations shows that DFT tends to overdelocalize the electrons in defected SrFeO 3 and LaFeO 3 . Finally, calculations with DMC and local density approximation yield similar vacancy migration energy barriers, indicating that steric/electrostatic effects mainly determine migration barriers in these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗