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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

A simplified model for the viscous crossflow in a slotted test section

A simplified physical model was constructed which simulates the viscous crossflow in a fluid layer near the slots at a fixed streamwise location in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical solutions of the two dimensional, incompressible Navier-Stokes equations in stream function and vorticity, which govern the model flow, were obtained. Fairly general slot geometry was incorporated by means of the Thompson-Thames-Mastin transformation. An approximate factorization scheme with cyclic acceleration parameters was employed to solve a finite difference analog of the stream function equation. The vorticity equation was numerically solved with a modified version of the classical alternating direction implicit scheme. Although no quantitative assessment of solution accuracy can be made, numerical results for variations in incremental wall pressure around the slat are at least qualitatively similar to some experimental results.

Cooke, C. H.↗
NASA NTRS1981

A simplified model for the viscous crossflow in a slotted test section

A simplified physical model is constructed which simulates the viscous crossflow in a fluid layer near the slots at a fixed streamwise location in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical solutions of the two-dimensional, incompressible Navier-Stokes equations in stream function and vorticity, which govern the model flow, are obtained. Fairly general slot geometry is incorporated by means of the Thompson-Thames-Mastin transformation. An approximate factorization scheme with cyclic acceleration parameters is employed to solve a finite difference analog of the stream function equation. The vorticity equation is numerically solved with a modified version of the classical alternating direction implicit (ADI) scheme. Although no quantitative assessment of solution accuracy can be made, numerical results for variations in incremental wall pressure around the slot are at least qualitatively similar to some experimental results of Berndt and Sorenson.

Cooke, C. H.↗
NASA NTRS1988

Estimating phase integrals - A generalization of Russell's law

An attempt is made to demonstrate, by means of Hapke's (1981) photometric function, that a simple and reliable method exists for the estimation of phase integrals from limited higher phase angle measurements. This method is a generalization of the approximation first proposed by Russell (1916) and more recently treated by Veverka (1971). It is shown that this generalization of Russell's law can employ observations anywhere in the range of phase angles from 40 to 90 deg; optimum estimates of q will be obtained if the data near 70 deg are obtainable.

Verbiscer, Anne J.↗
DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

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

Bloch and Bethe Ansätze for the Harper model: A butterfly with a boundary

Based on a recent generalization of Bloch's theorem, herein we present a Bloch Ansatz for the Harper model with an arbitrary rational magnetic flux in various geometries, and solve the associated Ansatz equations analytically. In the case of a cylinder and a particular boundary condition, the energy spectrum of edge states has no dependence on the length of the cylinder, which allows us to construct a quasi-one-dimensional edge theory that is exact and describes two edges simultaneously. We prove that energies of bulk states, generating the so-called Hofstadter's butterfly, depend on a single geometry-dependent spectral parameter and have exactly the same functional form for the cylinder and the torus with general twisted boundary conditions, and argue that the (edge) bulk spectrum of a semi-infinite cylinder in an irrational magnetic field is (the complement of) a Cantor set. Finally, realizing that the bulk projection of the Harper Hamiltonian is a linear form over a deformed Weyl algebra, we introduce a Bethe Ansatz valid for both cylinder and torus geometries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

General contraction of Gaussian basis sets. Part 2: Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) 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 outmost 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 (ANO) sets.

Almloef, Jan↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2023

Nonequilibrium design strategies for functional colloidal assemblies

Here, we use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.

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

Pairwise connected tensor network representation of path integrals

It has been recently shown how the tensorial nature of real-time path integrals (PIs) involving the Feynman-Vernon influence functional can be utilized with matrix product states, taking advantage of the finite length of the bath-induced memory. Tensor networks (TNs) promise to provide a unified language to express the structure of a PI. A generalized TN specifically incorporating the pairwise interaction structure of the influence functional and its invariance with respect to the average forward-backward position or the sojourn value in the form of the blip representation is derived and implemented. This pairwise connected TNPI (PC-TNPI) is illustrated through applications to typical spin-boson problems and explorations of the differences caused by the exact form of the spectral density. The storage and performance scalings are reported, showing the compactness of the representation and the efficiency of the contraction process. Finally, taking advantage of the compressed representation, the viability of using PC-TNPI for simulating multistate problems is demonstrated. The PC-TNPI structure can be shown to yield other TN algorithms currently in use. Consequently, it should be possible to use it as a starting point for deriving other optimized procedures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrolysis of Acetamide on Low-Index CeO 2 Surfaces: Ceria as a Deamidation and General De-esterification Catalyst

Using DFT calculations and acetamide as the main example, we show that ceria is a potential catalyst for the hydrolysis of amide and similar bonds. The overall reaction is endergonic in the gas phase, yielding acetic acid and ammonia, but is slightly exergonic in the aqueous phase, which facilitates ionization of the products (CH 3 COO – and NH 4 + ). Neighboring Ce and O sites on the CeO 2 (111), (110), and (100) facets are conducive to the formation of an activated metastable tetrahedral intermediate (TI) complex, followed by C–N bond scission. With van der Waals and solvation effects taken into account, the overall reaction energetics is found to be most favorable on the (111) facet as desorption of acetic acid is much more uphill energetically on (110) and (100). We further suggest that the Ce–O–Ce sites on ceria surfaces can activate X(=Y)–Z type bonds in amides, amidines, and carboxylate and phosphate esters, among many others that we term “generalized esters”. A Brønsted-Evans–Polanyi relationship is identified correlating the stability of the transition and final states of the X–Z generalized ester bond scission. A simple descriptor (ΣΔχ) based on the electronegativity of the atoms that constitute the bond (X, Y, Z) versus those of the catalytic site (O, Ce, Ce) captures the trend in the stability of the transition state of generalized ester bond scission and suggests a direction for modifying ceria for targeting specific organic substrates.

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

Theoretical Insights for Improving the Schottky-Barrier Height at the Ga2O3/Pt Interface

In this work, we study the Schottky-barrier height (SBH) at the junction between ..beta..-Ga2O3 and platinum, a system of great importance for the next generation of high-power and high-temperature electronic devices. Specifically, we obtain interfacial atomic structures at different orientations using our structure-matching algorithm and compute their SBH using electronic structure calculations based on hybrid density-functional theory.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,P↗
DOE OSTI2021

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

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

Bootstrapping the String Kawai-Lewellen-Tye Kernel

We show that a generalized version of the 4-point string theory KLT double-copy map is the most general solution to the minimal-rank double-copy bootstrap in effective field theory. This follows from significant restrictions of the 4-point map resulting from the 6-point bootstrap analysis. The generalized 4-point double-copy map is defined by a function with only two parameters times a simple function that is symmetric in s,t,u. The two parameters can be interpreted as independent choices of α', one for each of the two sets of amplitudes double-copied with the map. Specifically, each of those two sets of amplitudes must obey either the string monodromy relations or the field theory KK & BCJ relations; there are no other options. We propose a closed form of the new double-copy map that interpolates between the original KLT string double-copy and the open & closed string period integrals. The construction clarifies the "single-valued projection" property of the Riemann zeta-function values for the 4-point string theory double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Landau damping in space plasmas

The Landau damping of electrostatic Langmuir waves and ion-acoustic waves in a hot, isotropic, nonmagnetized, generalized Lorentzian plasma is analyzed using the modified plasma dispersion function. Numerical solutions for the real and imaginary parts of the wave frequency omega sub 0 - (i)(gamma) have been obtained as a function of the normalized wave number (k)(lambda sub D), where lambda sub D is the electron Debye length. For both particle distributions the electrostatic modes are found to be strongly damped at short wavelengths. At long wavelengths, this damping becomes less severe, but the attenuation of Langmuir waves is much stronger for a generalized Lorentzian plasma than for a Maxwellian plasma. It is concluded that Landau damping of ion-acoustic waves is only slightly affected by the presence of a high energy tail, but is strongly dependent on the ion temperature.

Thorne, Richard M.↗