Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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,531 records · Page 85

Model for rotating and nonuniform planetary exospheres.

The model neutral exosphere with a uniformly rotating exobase is generalized by allowing variations in exobase density and temperature which characterize the thermosphere just below the base. The corresponding velocity distribution function, satisfying the collisionless Boltzmann equation, is constructed and used to form a general expression for the velocity moments. Resulting density profiles of rotating exospheres with nonuniform densities and temperatures on the exobase are compared with corresponding nonrotating exospheres. Density enhancements due to rotation are found to be greatest above regions of exobase density or temperature minima. Equatorial density enhancements of terrestrial hydrogen, resulting from rotation, are estimated to be 15%-17% at altitudes of 10-20 earth radii. Corresponding increases in terrestrial helium are 30%-50% on the equator at altitudes of 0.7 to 1 earth radii even when there is a polar density bulge in the barosphere.

Hartle, R. E.↗

Quantum Nature of Dielectric Laser Accelerators

Dielectric laser accelerators (DLAs) hold great promise for producing economic and compact on-chip radiation sources. On-chip DLAs benefit from fabrication capabilities of the silicon industry and from breakthroughs in silicon-photonic nanostructures to enhance the interaction between particles and laser fields. Seemingly unrelated recent advances in the quantum interactions of electrons and light have raised interest in the underlying classical-quantum correspondence principle at the foundations of electron acceleration. Here, we present the observation of the underlying quantum nature of DLAs: observing quantized peaks in the electron-energy spectra. Our findings demonstrate quasi-phase-matching between an electron wave function and a light wave, which also demonstrates the role of the quantum wave function in the inverse Smith-Purcell effect. We harness the capabilities of an ultrafast transmission electron microscope (UTEM) to maintain a long electron-light interaction length extending over hundreds of periods of the laser pulse, mediated by a silicon-photonic nanograting DLA. The UTEM is shown as a new platform for characterization of future DLA concepts. The results raise fundamental questions regarding the role of quantum mechanics in DLA design, and more generally about the prospects of manipulating particles’ quantum wave functions in accelerator physics.

43 PARTICLE ACCELERATORS↗

A hill-sliding strategy for initialization of Gaussian clusters in the multidimensional space

A hill sliding technique was devised to extract Gaussian clusters from the multivariate probability density estimate of sample data for the first step of iterative unsupervised classification. Each cluster was assumed to posses a unimodal normal distribution. A clustering function proposed distinguished elements of a cluster under formation from the rest in the feature space. Initial clusters were extracted one by one according to the hill sliding tactics. A dimensionless cluster compactness parameter was proposed as a universal measure of cluster goodness and used satisfactorily in test runs with LANDSAT multispectral scanner data. The normalized divergence, defined by the cluster divergence divided by the entropy of the entire sample data, was utilized as a general separability measure between clusters. An overall clustering objective function was set forth in terms of cluster covariance matrices, from which the cluster compactness measure could be deduced. Minimal improvement of initial data partitioning was evaluated by this objective function in eliminating scattered sparse data points. The hill sliding clustering technique developed herein has the potential applicability to decomposition any multivariate mixture distribution into a number of unimodal distributions when an appropriate distribution function to the data set is employed.

Park, J. K.↗

Analysis of high-aspect-ratio jet-flap wings of arbitrary geometry

An analytical technique to compute the performance of an arbitrary jet-flapped wing is developed. The solution technique is based on the method of Maskell and Spence in which the well-known lifting-line approach is coupled with an auxiliary equation providing the extra function needed in jet-flap theory. The present method is generalized to handle straight, uncambered wings of arbitrary planform, twist, and blowing (including unsymmetrical cases). An analytical procedure is developed for continuous variations in the above geometric data with special functions to exactly treat discontinuities in any of the geometric and blowing data. A rational theory for the effect of finite wing thickness is introduced as well as simplified concepts of effective aspect ratio for rapid estimation of performance.

Lissaman, P. B. S.↗

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗

Photonics on the Mission to Mars

Human missions to Mars present some unique challenges for photonics devices. These devices will have exposure to many different space environments. During assembly they will be exposed to the Earth orbiting environment. Upon departure they will need to function through the Earth's Van Allen Radiation Belt. While the general interplanetary environment is less challenging than the radiation belt, they will operate in this environment for 18 months, subject to sudden saturation from solar flares. These components must continue to function properly through these saturation events presenting quite a challenge to photonic components, both optical and electronic. At Mars, the orbital environment is more benign than the Earth's. Components used as part of the landing vehicles must also deal with the pervasive dust environment for 3 - 6 months. These assembly and mission execution environments provide every form of space environmental challenges to photonic components. This paper will briefly discuss each environment and the expectations on the components for successful operation over the life of the mission.

Watson, Michael D.↗

Propagation of priors for more accurate and efficient spectroscopic functional fits and their application to ferroelectric hysteresis

Multi-dimensional spectral-imaging is a mainstay of the scanning probe and electron microscopies, micro-Raman, and various forms of chemical imaging. In many cases, individual spectra can be fit to a specific functional form, with the model parameter maps, providing direct insight into material properties. Since spectra are often acquired across a spatial grid of points, spatially adjacent spectra are likely to be similar to one another; yet, this fact is almost never used when considering parameter estimation for functional fits. On datasets tried here, we show that by utilizing proximal information, whether it be in the spatial or spectral domains, it is possible to improve the reliability and increase the speed of such functional fits by ~2-3x, as compared to random priors. We explore and compare three distinct new methods: (1) spatially averaging neighborhood spectra, and propagating priors based on functional fits to the averaged case, (2) hierarchical clustering-based methods where spectra are grouped hierarchically based on response, with the priors propagated progressively down the hierarchy, and (3) regular clustering without hierarchical methods with priors propagated from fits to cluster means. Our results highlight that utilizing spatial and spectral neighborhood information is often critical for accurate parameter estimation in noisy environments, which we show for ferroelectric hysteresis loops acquired on a prototypical PbTiO3 thin film with piezoresponse spectroscopy. This method is general and applicable to any spatially measured spectra where functional forms are available. Examples include exploring the superconducting gap with tunneling spectroscopy, using the Dynes formula, or current-voltage curve fits in conductive atomic force microscopy mapping. Here we explore the problem for ferroelectric hysteresis, which, given its large parameter space, constitutes a more difficult task than, for example, fitting current-voltage curves with a Schottky emission formula.

42 ENGINEERING↗

Transition Metal Activation And Functionalization Of Carbon-Hydrogen And Carbon-Carbon Bonds

This project has as its overall goal improvement in the intelligent use of our energy resources, specifically petroleum derived products, and is aimed at the development of new routes for the manipulation of C-H and C-C bonds. During this project period, our research is focused on the following general goals: (1) C-H bond activation reactions in functionalized substrates, including alkylnitriles, ketones, esters, ethers, fluoromethanes, alcohols, amines, and alkynes, and fundamental studies of their effects on metal-carbon bond strengths using trispyrazolylborate-rhodium (TpRh) complexes; (2) Hydrogenation reactions using first-row transition metals; and (3) C–CN bond activation in aryl-nitriles and allyl-nitriles.

02 PETROLEUM↗

General contraction of Gaussian basis sets. II - Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.

Almlof, Jan↗

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗

Atomic Fingerprinting of Heteroatoms Using Noncontact Atomic Force Microscopy

We report immense strides have been made in increasing the resolution of scanning probe microscopy. Noncontact atomic force microscopy (n-AFM) now offers one the ability to characterize and visualize single molecules with subatomic resolution. Specifically, nc-AFM with a carbon monoxide (CO) functionalized tip has the ability to discriminate functional groups (–C≡C–, –CH 2 , –C$=$O, …), although discriminating atomic species often remains as an ongoing challenge. Here, real-space pseudopotentials constructed within density functional theory are employed to accurately simulate nc-AFM images of molecules containing heteroatoms (S, I, and N) within dibenzothiophene (DBT), 2-iodotriphenylene (ITP), acridine (ACR) and ferrous phthalocyanine (FePc). It is found that S and I atoms can be easily identified from C based on their unique features. For N atoms, a use of tip functionalization is proposed to effectively discriminate them from C atoms.

36 MATERIALS SCIENCE↗

Some results on contractive mappings as related to pattern recognition

Several of the techniques used in pattern recognition are reformulated as the problem of determining fixed points of a function. If x sub 0 is a fixed point of f and if f is contractive at x sub 0, then, for any y belonging to a sufficiently small neighborhood of x sub 0 the orbit of y will converge to x sub 0. Several general results regarding contractive mappings are developed with emphasis on functions.

Brown, D. R.↗

Analytical characterization of AOTV perigee aerothermodynamic regime

In preliminary design studies of AOTV (Aeroassisted Orbital Transfer Vehicle) configurations, it is desirable to have a simple analytical method for characterizing the aerothermodynamic regime at skip trajectory perigee as a function of vehicle L/D and m/CDA (ballistic coefficient). The present study derives an approximate perigee solution from the generalized equations of motion. For a prescribed entry velocity vector, the solution determines a Zeta function (proportional to density divided by m/CDA) at perigee as a parametric function of the exit velocity vector and (constant) L/D. The perigee Zeta function then determines perigee density (or altitude) as a parametric function of m/CDA. The solution allows the following classic aerothermodynamic parameters to be determined at perigee as parametric functions of m/CDA, L/D, and the exit velocity vector: Reynolds Number (viscous effects), Mach Number (compressibility effects), Knudsen Number (rarefaction effects), sphere and disk bow shock standoff distance, Damkohler Number (relaxation effects), viscous correlation parameter (viscid-inviscid interactions), and Stanton Number (convective heat transfer). Results of the analysis are given for low L/D and mid L/D AOTV configurations on return from geosynchronous and L5 orbits. It is concluded the method successfully provides preliminary estimates of the aerothermodynamic parameters through the use of simple algebraic equations and plots.

Desautel, D.↗

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of the thermodynamic properties of a mixture of gases as a function of temperature and pressure

The evaluation of the thermodynamic properties of a gas mixture can be performed using a generalized correlation which makes use of the second virial coefficient. This coefficient is based on statistical mechanics and is a function of temperature and composition, but not of pressure. The method provides results accurate to within 3 percent for gases which are nonpolar or only slightly polar. When applied to highly polar gases, errors of 5 to 10 percent may result. For gases which associate, even larger errors are possible. The sequences of calculations can be routinely programmed for a digital computer. The thermodynamic properties of a mixture of neon, argon and ethane were calculated by such a program. The result will be used for the design of the gas replenishment system for the Energetic Gamma Ray Experiment Telescope.

Colon, G.↗

Descriptor-Based Analysis of Atomic Layer Deposition Mechanisms on Spinel LiMn 2 O 4 Lithium-Ion Battery Cathodes

Protective coatings have been shown to effectively suppress Mn ion dissolution from the spinel LiMn 2 O 4 lithium-ion battery cathode by stabilizing the surface against undesired side reactions with the electrolyte. In spite of extensive study, however, there remains a lack of atomic-scale understanding of how such coatings are deposited, and no molecular-level descriptor to predict trends in deposition mechanisms has been identified. We have recently shown that Al 2 O 3 coatings grown by atomic layer deposition (ALD) with alternating trimethylaluminum (TMA) and water exposures exhibit submonolayer growth because of precursor decomposition on the lithium manganate spinel (LMO) surface during early ALD pulses. In the present work, we elucidate the underlying mechanisms of this Al 2 O 3 ALD process using density functional theory (DFT) calculations and X-ray photoelectron spectroscopy (XPS) experiments, and we introduce a generalized descriptor-based framework to understand the resulting trends across a spectrum of surface structures and functionalities. We demonstrate that all decomposition products, including CH 3 -aluminum adducts and dissociated CH 3 groups, are Lewis bases and are coordinated to oxygen atoms on the LMO surface, leading to charge transfer to Lewis acidic Mn 3d states. Inert-transfer XPS supports these theoretical predictions, showing an increase in near-surface Mn 3+ content following TMA exposure and shifts in C is spectra consistent with C-O bond formation. We extend the DFT studies to various low- and high-index LMO surface facets, as a proxy for tuning the Lewis acid-base interactions between surface-bound CH 3 * and near-surface Mn ions. The thermochemistry for TMA reactions on these chemically distinct LMO surfaces demonstrates that ALD is structure-sensitive and that there is higher reactivity for TMA decomposition and Al 2 O 3 nucleation near LMO steps and defects. Motivated by the Lewis basic character of the decomposition products, we introduce the oxygen vacancy formation energy as a descriptor for decomposition energetics, and we demonstrate that all energetics are correlated to this quantity through the number of electrons that are transferred along the reaction coordinate. Finally, based on these findings, we hypothesize that improved electrochemical cycling with only 1-2 ALD cycles may be due to selective passivation of defect sites on the LMO surface that are more susceptible to Mn dissolution, and we suggest that similar descriptor-based analyses could be useful for the study of other ALD coatings on oxide substrates.

25 ENERGY STORAGE↗