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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Development of TDLAS for Water Vapor Concentration Measurement during VHTR Steam Ingress

Very-high-temperature reactors (VHTR) utilize a graphite-moderated, helium-cooled design that provides higher efficiencies and a low-carbon source of process heat (>700 deg. C) for many industrial applications. Safety concerns have led to extensive accident condition investigations. Steam generator heat-exchanger tube rupture accidents (SGTR) can lead to increased reactor core steam density, which causes neutron slowing-down and chemical corrosion of the graphite fuel elements and structural columns. It is important to be able to computationally model the propagation of steam throughout the lower plenum and nuclear reactor core, but experimental data is needed to validate such models. The primary objective of this project is to instrument a scaled down experiment of a VHTR designed by Ohio State University (OSU), with tunable diode laser absorption spectroscopy (TDLAS) diagnostics to measure steam concentration changes and map how water vapor propagates through the test facility's lower plenum and core. The basic principle of absorption spectroscopy is based on a molecular species ability to absorb photons of specific energies. A collimated laser beam of known energy is passed through a test section containing a sample gas. When a transition of the absorbing gas species resonates with the laser frequency, a certain amount of energy will be absorbed based on the rovibrational characteristics of the molecule. The intensity of the transmitted laser light can be measured by a detector after exiting the test section. The Beer-Lambert Law of absorbance is used to recover this information relating beam attenuation to the system dimensions, gas composition, pressure, and temperature. The line-strength is a function of temperature alone and relates the fraction of molecules in a given rovibrational state through the partition function. The partition function describes the statistical properties of a system in thermodynamic equilibrium where the states obey a Boltzmann distribution. Thus, the number density of molecules in a particular state can be determined from the total internal partition function at a particular temperature. The lineshape function of a given absorption transition is derived from the physical interaction between the absorbing species and is dependent on temperature and pressure. Random thermal motion of absorbing molecules result in broadening of the linewidth, called Doppler broadening, which follows a Gaussian distribution. Molecular collisions result in broadening and shifting of the linewidth, called pressure broadening, which follows a Lorentzian distribution. The overall lineshape results in a Voigt profile, which is a convolution of the Gaussian and Lorentzian distributions. The Voigt profile does not have a simple analytical form and is most often calculated numerically or approximated by complex functions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

First-Principles-Based Machine Learning Models for Phase Behavior and Transport Properties of CO 2

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN and SCAN-rvv10 approximations of density functional theory. Here we employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

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

Spectral-density estimation with the Gaussian integral transform

The spectral-density operator $\hat{ρ}(ω) = δ(ω–\hat{H})$ plays a central role in linear response theory as its expectation value, the dynamical response function, can be used to compute scattering cross sections. In this work, we describe a near optimal quantum algorithm providing an approximation to the spectral density with energy resolution $\Delta$ and error $\epsilon$ using $O(\sqrt{\text{log}_2 (1/ε)[\text{log}_2 (1 / Δ) + \text{log}_2 (1/ε)]/ Δ)}$ operations. This is achieved without using expensive approximations to the time-evolution operator, but instead exploiting qubitization to implement an approximate Gaussian integral transform of the spectral density. Finally, we also describe appropriate error metrics to assess the quality of the spectral function approximations more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Validity of virial theorem in all-electron mixed basis density functional, Hartree–Fock, and GW calculations

In this paper, we calculate kinetic and potential energy contributions to the electronic ground-state total energy of several isolated atoms (He, Be, Ne, Mg, Ar, and Ca) by using the local density approximation (LDA) in density functional theory, the Hartree–Fock approximation (HFA), and the self-consistent GW approximation (GWA). To this end, we have implemented self-consistent HFA and GWA routines in our all-electron mixed basis code, TOMBO. We confirm that virial theorem is fairly well satisfied in all of these approximations, although the resulting eigenvalue of the highest occupied molecular orbital level, i.e., the negative of the ionization potential, is in excellent agreement only in the case of the GWA. We find that the wave function of the lowest unoccupied molecular orbital level of noble gas atoms is a resonating virtual bound state, and that of the GWA spreads wider than that of the LDA and thinner than that of the HFA.

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

Controlling dynamical thermal transport of biased bilayer graphene by impurity atoms

We address the dynamical thermal conductivity of biased bilayer graphene doped with acceptor impurity atoms for AA-stacking in the context of tight binding model Hamiltonian. The effect of scattering by dilute charged impurities is discussed in terms of the self-consistent Born approximation. Green’s function approach has been exploited to find the behavior of thermal conductivity of bilayer graphene within the linear response theory. We have found the frequency dependence of thermal conductivity for different values of concentration and scattering strength of dopant impurity. Also the dependence of thermal conductivity on the impurity concentration and bias voltage has been investigated in details.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2018

High-Resolution, Ground-Based Observations of the Lunar Sodium Exosphere During the Lunar Atmosphere and Dust Environment Explorer (LADEE) Mission

We present the first comprehensive set of lunar exospheric line width and line width derived effective temperatures as a function of lunar phase (66° waxing phase to 79° waning phase). Data were collected between November 2013 and May 2014 during six observing runs at the National Solar Observatory McMath-Pierce Solar Telescope by applying high-resolution Fabry-Perot spectroscopy (R ~ 180,000) to observe emission from exospheric sodium (5,889.9509 Å, D2 line). The 3-arc min field of view of the instrument, corresponding to ~336 km at the mean lunar distance (384,400 km), was positioned at several locations off the lunar limb; only equatorial observations taken out to 950 km are presented here. We find the sodium effective temperature distribution to be approximately a symmetric function of lunar phase with respect to full Moon. Within magnetotail passage we find temperatures in the range of 2500-9000 K. For phase angles greater than 40deg we find that temperatures flatten out to ~1700 K.

Subsolar↗
DOE OSTI2013

Computing the partition function, ensemble averages, and density of states for lattice spin systems by sampling the mean

An algorithm to approximately calculate the partition function (and subsequently ensemble averages) and density of states of lattice spin systems through non-Monte-Carlo random sampling is developed. This algorithm (called the sampling-the-mean algorithm) can be applied to models where the up or down spins at lattice nodes interact to change the spin states of other lattice nodes, especially non-Ising-like models with long-range interactions such as the biological model considered here. Because it is based on the Central Limit Theorem of probability, the sampling-the-mean algorithm also gives estimates of the error in the partition function, ensemble averages, and density of states. Easily implemented parallelization strategies and error minimizing sampling strategies are discussed. The sampling-the-mean method works especially well for relatively small systems, systems with a density of energy states that contains sharp spikes or oscillations, or systems with little a priori knowledge of the density of states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Theoretical and experimental study of the dynamic response of absorber-based, micro-scale, oscillatory probes for contact sensing applications

This paper presents two models for predicting the frequency response of micro-scale oscillatory probes. These probes are manufactured by attaching a thin fiber to the free end of one tine of a quartz tuning fork oscillator. In these studies, the attached fibers were either 75 μm diameter tungsten or 7 μm diameter carbon with lengths ranging from around 1 to 15 mm. The oscillators used in these studies were commercial 32.7 kHz quartz tuning forks. The first theoretical model considers lateral vibration of two beams serially connected and provides a characteristic equation from which the roots (eigenvalues) are extracted to determine the natural frequencies of the probe. A second, lumped model approximation is used to derive an approximate frequency response function for prediction of tine displacements as a function of a modal force excitation corresponding to the first mode of the tine in the absence of a fiber. These models are used to evaluate the effect of changes in both length and diameter of the attached fibers. Theoretical values of the natural frequencies of different modes show an asymptotic relationship with the length and a linear relationship with the diameter of the attached fiber. Similar results are observed from experiment, one with a tungsten probe having an initial fiber length of 14.11 mm incrementally etched down to 0.83 mm, and another tungsten probe of length 8.16 mm incrementally etched in diameter, in both cases using chronocoulometry to determine incremental volumetric material removal. The lumped model is used to provide a frequency response again reveals poles and zeros that are consistent with experimental measurements. Finite element analysis shows mode shapes similar to experimental microscope observations of the resonating carbon probes. This model provides a means of interpreting measured responses in terms of the relative motion of the tine and attached fibers. Of particular relevance is that, when a “zero” is observed in the response of the tine, one mode of the fiber is matched to the tine frequency and is acting as an absorber. This represents an optimal condition for contact sensing and for transferring energy to the fiber for fluid mixing, touch sensing, and surface modification applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1992

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗
DOE OSTI2014

Analytical solution of the Poisson-Nernst-Planck equations for an electrochemical system close to electroneutrality

Single charge densities and the potential are used to describe models of electrochemical systems. These quantities can be calculated by solving a system of time dependent nonlinear coupled partial differential equations, the Poisson-Nernst-Planck equations. Assuming small deviations from the electroneutral equilibrium, the linearized and decoupled equations are solved for a radial symmetric geometry, which represents the interface between a cell and a sensor device. The densities and the potential are expressed by Fourier-Bessels series. The system considered has a ratio between the Debye-length and its geometric dimension on the order of 10{sup −4} so the Fourier-Bessel series can be approximated by elementary functions. The time development of the system is characterized by two time constants, τ{sub c} and τ{sub g}. The constant τ{sub c} describes the approach to the stationary state of the total charge and the potential. τ{sub c} is several orders of magnitude smaller than the geometry-dependent constant τ{sub g}, which is on the order of 10 ms characterizing the transition to the stationary state of the single ion densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A primal–dual algorithm for risk minimization

In this paper, we develop an algorithm to efficiently solve risk-averse optimization problems posed in reflexive Banach space. Such problems often arise in many practical applications as, e.g., optimization problems constrained by partial differential equations with uncertain inputs. Unfortunately, for many popular risk models including the coherent risk measures, the resulting risk-averse objective function is nonsmooth. Here, this lack of differentiability complicates the numerical approximation of the objective function as well as the numerical solution of the optimization problem. To address these challenges, we propose a primal–dual algorithm for solving large-scale nonsmooth risk-averse optimization problems. This algorithm is motivated by the classical method of multipliers and by epigraphical regularization of risk measures. As a result, the algorithm solves a sequence of smooth optimization problems using derivative-based methods. We prove convergence of the algorithm even when the subproblems are solved inexactly and conclude with numerical examples demonstrating the efficiency of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Gaussian kernel density functions for compositional quantification in atom probe tomography

Highlights: • A methodology for compositional quantification from atom probe tomography is proposed. • A probability density function p(x) given discrete data points is constructed. • Each atom approximated by a Gaussian function with a particular smoothing parameter. • Density Kernel functions are used to find the optimum smoothing parameter. • The method is validated by using simulated data of semiconductor materials. - Abstract: Atom probe tomography (APT) has the ability to identify the nature and position of single atoms in a material with an almost 3D atomic resolution. However, the quantification of the material composition requires an appropriate treatment of the discrete APT data. When the amount of atoms selected to quantify the composition is relatively small, the spatial resolution is enhanced but the statistical error worsens. Conversely, the increase in the number of atoms considered reduces the spatial resolution, but improving the statistics. We have developed a methodology to reach an optimum equilibrium between the positional and the statistical error in an unsupervised form using Gaussian Kernel density functions. The validity of the method has been tested using APT simulated data of semiconductor materials. It has been proved that the chemical quantification in these materials requires the appropriate selection of the smoothing parameter, obtained without user intervention. In addition, the results have been compared to the usual techniques for composition measurement from APT data (voxelization and proximity histograms), showing better precision for high spatial resolution. This work supplies a data driven methodology based in Gaussian Kernel density functions for the accurate quantification of the composition from APT data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

Rain cell size statistics derived from radar observations at Wallops Island, Virginia

An investigation of two-dimensional rain cell size statistics has analyzed regression relations relating radar-determined rain rates to disdrometer data. This has yielded least-squares fits of radar reflectivity factors, and the application of a contouring program has generated 22,000 contours in which each isopleth belongs to predefined rain-rate intervals. An abundance of total and cell contours were observed belonging to all rain-rate categories. Both the computed number distributions and the conditional cumulative distributions as a function of contour diameter were found to be represented with good approximation by given exponential functions.

Goldhirsh, J.↗
DOE OSTI2024

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗
DOE OSTI2014

Description of plasmon-like band in silver clusters: The importance of the long-range Hartree-Fock exchange in time-dependent density-functional theory simulations

Absorption spectra of Ag{sub 20} and Ag{sub 55}{sup q} (q = +1, −3) nanoclusters are investigated in the framework of the time-dependent density functional theory in order to analyse the role of the d electrons in plasmon-like band of silver clusters. The description of the plasmon-like band from calculations using density functionals containing an amount of Hartree-Fock exchange at long range, namely, hybrid and range-separated hybrid (RSH) density functionals, is in good agreement with the classical interpretation of the plasmon-like structure as a collective excitation of valence s-electrons. In contrast, using local or semi-local exchange functionals (generalized gradient approximations (GGAs) or meta-GGAs) leads to a strong overestimation of the role of d electrons in the plasmon-like band. The semi-local asymptotically corrected model potentials also describe the plasmon as mainly associated to d electrons, though calculated spectra are in fairly good agreement with those calculated using the RSH scheme. Our analysis shows that a portion of non-local exchange modifies the description of the plasmon-like band.

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

Green’s function coupled cluster formulations utilizing extended inner excitations

Here, we analyze new approximations of the Green’s function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i( n, m ) method, where the excitation level of the inner auxiliary operators ( m ) used to describe the ionization potential and electron affinity effects in the N-1 and N+1 particle spaces is higher than the excitation level ( n ) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called “n + 1” rule in this category [or the GFCC-i( n, n + 1) method], which states that in order to maintain size-extensivity of the Green’s function matrix elements, the excitation level of inner auxiliary operators X p (ω) and Y q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N 2 and CO molecules in the inner and outer valence regimes. Lastly, in comparison with the Green’s function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

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

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1995

Hubble space telescope observations of young star clusters in NGC-4038/4039, 'the antennae' galaxies

New, high-resolution images of the disks of NGC 4038/4039 obtained with the Wide Field Camera of the Hubble Space Telescope (HST) are presented. NGC 4038/4039, nicknamed The Antennae, is a prototypical example of a pair of colliding galaxies believed to be at an early stage of a merger. Down to the limiting magnitude of V approximately 23 mag, the HST images reveal a population of over 700 blue pointlike objects within the disks. The mean absolute magnitude of these objects is M(sub V) = -11 mag, with the brightest objects reaching M(sub V) approximately -15. Their mean apparent color indices ar U - V = -0.7 mag and V - 1 = 0.8 mag on the Johnson UVI passband system, while their mean indices corrected for internal reddening are (u - v)(sub 0) = -1.0 mag and (V - I(sub 0) = 0.5. Their mean effective radius, determined from slightly resolved images, is 18 pc (for H(sub 0) = 50 km/s /Mpc). Based on their luminosities and resolution, most of these objects cannot be individual stars, but are likely young compact star clusters. The brighter ones are similar to the objects found in NGC 1275 and NGC 7252, which appear to be young globular clusters formed during recent galazy mergers. Based on their U - V and V - I colors, the brightest, bluest clusters of NGC 4038/4039 appear to be less than 10 Myr old. Most of these bright clusters are relatively tightly clustered themselves, with typically a dozen individual clusters belonging to a complex identified as a giant H II region from ground-based observations. The cluster luminosity function (LF) is approximately a power law, phi(L)dL proportional to L(exp -1.78+/-0.05)dL, with no hint of a turnover at fainter magnitudes. This power-law shape agrees with the LF of Magellanic Cloud clusters and Galactic open clusters, but differs from the LF of old globular cluster systems that is typically Gaussian with a Full Width at Half Maximum (FWHM) of approximately 3 mag. Besides the blue clusters, we also find about a dozen extremely red objects with V - I greater than 3.0. The highest number density of these red objects is found in the SE quadrant, where star formation appears to be most recent. We propose that these objects may be very young star clusters still embedded in their placental dust cocoons.

Whitmore, Bradley C.↗