Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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

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

Materials Science↗
DOE OSTI2024

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

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

HVEM structures and mutants reveal distinct functions of binding to LIGHT and BTLA/CD160

HVEM is a TNF (tumor necrosis factor) receptor contributing to a broad range of immune functions involving diverse cell types. It interacts with a TNF ligand, LIGHT, and immunoglobulin (Ig) superfamily members BTLA and CD160. Assessing the functional impact of HVEM binding to specific ligands in different settings has been complicated by the multiple interactions of HVEM and HVEM binding partners. To dissect the molecular basis for multiple functions, we determined crystal structures that reveal the distinct HVEM surfaces that engage LIGHT or BTLA/CD160, including the human HVEM–LIGHT–CD160 ternary complex, with HVEM interacting simultaneously with both binding partners. Based on these structures, we generated mouse HVEM mutants that selectively recognized either the TNF or Ig ligands in vitro. Knockin mice expressing these muteins maintain expression of all the proteins in the HVEM network, yet they demonstrate selective functions for LIGHT in the clearance of bacteria in the intestine and for the Ig ligands in the amelioration of liver inflammation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

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

Solar p-mode oscillations as a tracer of radial differential rotation

Photoelectric observations of solar p-modes obtained with improved wavenumber and frequency resolution are presented. The observations are compared with model calculations of the p-modes, and the degree of spatial and temporal coherence of the observed wave pattern is investigated. It is found that the p-mode oscillations pervade the visible surface of the sun with a high degree of coherence in space and time, so that the whole complex pattern of standing waves with its nodes and antinodes can be regarded as a fixed pattern corotating with the solar surface layers. The p-modes are introduced as a tracer of solar rotational flow velocities. The equatorial differential rotation is estimated as a function of effective depth on the basis of the theoretical contribution functions for the p-modes recently derived by Ulrich et al. (1978). The results strongly indicate that the angular speed of rotation is not uniform even in the relatively shallow layer extending about 20,000 km below the photosphere.

Deubner, F.-L.↗
DOE OSTI2020

Crossing symmetry, transcendentality and the Regge behaviour of 1d CFTs

We develop the technology for Polyakov-Mellin (PM) bootstrap in one- dimensional conformal field theories (CFT 1 ). By adding appropriate contact terms, we bootstrap various effective field theories in AdS 2 and analytically compute the CFT data to one loop. The computation can be extended to higher orders in perturbation theory, if we ignore mixing, for any external dimension. We develop PM bootstrap for O(N) theories and derive the necessary contact terms for such theories (which also involves a new higher gradient contact term absent for N = 1). We perform cross-checks which include considering the diagonal limit of the 2d Ising model in terms of the 1d PM blocks. As an independent check of the validity of the results obtained with PM bootstrap, we propose a suitable basis of transcendental functions, which allows to fix the four-point correlators of identical scalar primaries completely, up to a finite number of ambiguities related to the number of contact terms in the PM basis. We perform this analysis both at tree level (with and without exchanges) and at one loop. We also derive expressions for the corresponding CFT data in terms of harmonic sums. Finally, we consider the Regge limit of one-dimensional correlators and derive a precise connection between the latter and the large-twist limit of CFT data. Exploiting this result, we study the crossing equation in the three OPE limits and derive some universal constraints for the large-twist limit of CFT data in Regge-bounded theories with a finite number of exchanges.

1/N expansion↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

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

HALE UAS Concept of Operations. Version 3.0

This document is a system level Concept of Operations (CONOPS) from the perspective of future High Altitude Long Endurance (HALE) Unmanned Aircraft Systems (UAS) service providers and National Airspace System (NAS) users. It describes current systems (existing UAS), describes HALE UAS functions and operations to be performed (via sample missions), and offers insight into the user s environment (i.e., the UAS as a system of systems). It is intended to be a source document for NAS UAS operational requirements, and provides a construct for government agencies to use in guiding their regulatory decisions, architecture requirements, and investment strategies. Although it does not describe the technical capabilities of a specific HALE UAS system (which do, and will vary widely), it is intended to aid in requirements capture and to be used as input to the functional requirements and analysis process. The document provides a basis for development of functional requirements and operational guidelines to achieve unrestricted access into the NAS. This document is an FY06 update to the FY05 Access 5 Project-approved Concept of Operations document previously published in the Public Domain on the Access 5 open website. This version is recommended to be approved for public release also. The updates are a reorganization of materials from the previous version with the addition of an updated set of operational requirements, inclusion of sample mission scenarios, and identification of roles and responsibilities of interfaces within flight phases.

Source record↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Giant gravitons, Harish-Chandra integrals, and BPS states in symplectic and orthogonal $\mathcal{N}$ = 4 SYM

We find generating functions for half BPS correlators in $\mathcal{N}$ = 4 SYM theories with gauge groups Sp(2N), SO(2N + 1), and SO(2N) by computing the norms of a class of BPS coherent states. These coherent states are built from operators involving Harish-Chandra integrals. Such operators have an interpretation as localized giant gravitons in the bulk of anti-de-Sitter space. This extends the analysis of [1] to Sp(2N), SO(2N + 1), and SO(2N) gauge theories. We show that we may use ordinary Schur functions as a basis for the sector of states with no cross-caps in these theories. This is consistent with the construction of these theories as orientifold projections of an SU(2N) theory. We make note of some relations between the symmetric functions that appear in the expansion of these coherent states and symplectic Schur functions. We also comment on some connections to Schubert calculus and Gromov-Witten invariants, which suggest that the Harish-Chandra integral may be extended to such problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

Do Bond Functions Help for the Calculation of Accurate Bond Energies?

The bond energies of 8 chemically bound diatomics are computed using several basis sets with and without bond functions (BF). The bond energies obtained using the aug-pVnZ+BF basis sets (with a correction for basis set superposition error, BSSE) tend to be slightly smaller that the results obtained using the aug-pV(n+I)Z basis sets, but slightly larger than the BSSE corrected aug-pV(n+I)Z results. The aug-cc-pVDZ+BF and aug-cc-pVTZ+BF basis sets yield reasonable estimates of bond energies, but, in most cases, these results cannot be considered highly accurate. Extrapolation of the results obtained with basis sets including bond functions appears to be inferior to the results obtained by extrapolation using atom-centered basis sets. Therefore bond functions do not appear to offer a path for obtaining highly accurate results for chemically bound systems at a lower computational cost than atom centered basis sets.

Bauschlicher, Charles W., Jr.↗
NASA NTRS2013

RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers

RUSHMAPS is a new onboard data reduction scheme that gives real-time access to key science parameters (e.g. moments) of a class of heliophysics science and/or solar system exploration investigation that includes plasma particle spectrometers (PPS), but requires moments reporting (density, bulk-velocity, temperature, pressure, etc.) of higher-level quality, and tolerates a lowpass (variable quality) spectral representation of the corresponding particle velocity distributions, such that telemetry use is minimized. The proposed methodology trades access to the full-resolution velocity distribution data, saving on telemetry, for real-time access to both the moments and an adjustable-quality (increasing quality increases volume) spectral representation of distribution functions. Traditional onboard data storage and downlink bandwidth constraints severely limit PPS system functionality and drive cost, which, as a consequence, drives a limited data collection and lower angular energy and time resolution. This prototypical system exploit, using high-performance processing technology at GSFC (Goddard Space Flight Center), uses a SpaceCube and/or Maestro-type platform for processing. These processing platforms are currently being used on the International Space Station as a technology demonstration, and work is currently ongoing in a new onboard computation system for the Earth Science missions, but they have never been implemented in heliospheric science or solar system exploration missions. Preliminary analysis confirms that the targeted processor platforms possess the processing resources required for realtime application of these algorithms to the spectrometer data. SpaceCube platforms demonstrate that the target architecture possesses the sort of compact, low-mass/power, radiation-tolerant characteristics needed for flight. These high-performing hybrid systems embed unprecedented amounts of onboard processing power in the CPU (central processing unit), FPGAs (field programmable gate arrays), and DSP (digital signal processing) elements. The fundamental computational algorithm de constructs 3D velocity distributions in terms of spherical harmonic spectral coefficients (which are analogous to a Fourier sine-cosine decomposition), but uses instead spherical harmonics Legendre polynomial orthogonal functions as a basis for the expansion, portraying each 2D angular distribution at every energy or, geometrically, spherical speed-shell swept by the particle spectrometer. Optionally, these spherical harmonic spectral coefficients may be telemetered to the ground. These will provide a smoothed description of the velocity distribution function whose quality will depend on the number of coefficients determined. Successfully implemented on the GSFC-developed processor, the capability to integrate the proposed methodology with both heritage and anticipated future plasma particle spectrometer designs is demonstrated (with sufficiently detailed design analysis to advance TRL) to show specific science relevancy with future HSD (Heliophysics Science Division) solar-interplanetary, planetary missions, sounding rockets and/or CubeSat missions.

Figueroa-Vinas, Adolfo↗
DOE OSTI2025

Optimizing stochastic algorithms for hadron correlation function computations in lattice QCD using a localized distillation basis

Distillation is a quark-smearing method for the construction of a broad class of hadron operators useful in lattice QCD computations and defined via a projection operator into a vector space of smooth gauge-covariant fields. A new orthonormal basis for this space is constructed which builds in locality. This basis is useful for the construction of stochastic methods to estimate the correlation functions computed in Monte Carlo calculations relevant for hadronic physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

PhenoProfiling: Mapping phenotypic outcomes to molecular determinants of biochemical activity

The aggregate genomes of the trillions of microorganisms within soil, animal hosts, and aquatic systems encode for an extensive functional capacity for myriad biochemical activities. C, N, P, and S metabolism, synthesis of signaling molecules and vitamins, mineralization, and other activities are essential to microbe, community, and plant physiology, and more broadly to climate, water, animal, flora, and human health. The current understanding of the molecular basis for the function of microbial communities stems primarily from comparative metagenomic and metatranscriptomic studies. These same tools are employed to ascertain the impacts to the community resulting from perturbations, such as climate change, emerging pollutants, fires, and seawater infiltration for environmental communities, and dietary changes, xenobiotic exposure, and various disease states for the human gut microbiome. Such studies can identify the potential for a specific function, but they cannot determine that a particular cell is functionally active, nor can they determine the molecular architecture required for function. In short, genes and transcripts alone fail to reveal the complex subcellular arrangement of proteins and molecules that elicit a given phenotype in a microbial cell.

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

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗