Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

Stabilized NADH as a Countermeasure for Jet Lag

Current remedies for jet lag (phototherapy, melatonin, stimulant, and sedative medications) are limited in efficacy and practicality. The efficacy of a stabilized, sublingual form of reduced nicotin amide adenine dinucleotide (NADH, ENADAlert, Menuco Corp.) as a countermeasure for jet lag was examined. Because NADH increases cellular production of ATP and facilitates dopamine synthesis, it may counteract the effects of jet lag on cognitive functioning and sleepiness. Thirty-five healthy, employed subjects participated in this double-blind, placebo-controlled study. Training and baseline testing were conducted on the West Coast before subjects flew overnight to the East Coast, where they would experience a 3-hour time difference. Upon arrival, individuals were randomly assigned to receive either 20 mg of sublingual stabilized ADH (n=18) or identical placebo tablets (n=17). All participants completed computer-administered tests (including CogScreen7) to assess changes in cognitive functioning, mood, and sleepiness in the morning and afternoon. Jet lag resulted in increased sleepiness for over half the participants and deterioration of cognitive functioning for approximately one third. The morning following the flight, subjects experienced lapses of attention in addition to disruptions in working memory, divided attention, and visual perceptual speed. Individuals who received NADH performed significantly better on 5 of 8 cognitive and psychomotor test measures (P less than or equal to 0.5) and showed a trend for better performance on the other three measures (P less than or equal to .l0). Subjects also reported less sleepiness compared with those who received placebo. No adverse effects were observed with NADH treatment. Stabilized NADH significantly reduced jet lag-induced disruptions of cognitive functioning, was easily administered, and was found to have no adverse side effects.

Kay, Gary G.↗

Feedback Resonating Control for a Wave Energy Converter

Through the use of advanced control techniques, wave energy converters (WECs) can achieve substantial increases in energy absorption. The motion of the WEC device is a significant contribution to the energy absorbed by the device. Reactive (complex conjugate) control maximizes the energy absorption due to the impedance matching. The issue with complex conjugate control is that, in general, the controller is noncausal, which requires prediction of the incoming waves. This article explores the potential of employing system identification techniques to build a causal transfer function that approximates the complex conjugate controller over a finite frequency band of interest. This approach is quite viable given the band-limited nature of ocean waves. The resulting controller is stable, and the average efficiency of the power captured by the causal controller in realistic ocean waves is 99%, when compared to the noncausal complex conjugate.

16 TIDAL AND WAVE POWER↗

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗

Spaceflight-Induced Cardiovascular Changes and Recovery During NASA's Functional Task Test

Microgravity-induced physiological changes could impair a crewmember s performance upon return to a gravity environment. The Functional Task Test (FTT) is designed to correlate these physiological changes to performance in mission-critical tasks. The Recovery from Fall/Stand Test (RFST) simulates one such task, measuring the ability to recover from a prone position and the cardiovascular response to orthostasis. The purpose of this study was to evaluate spaceflight-induced cardiovascular changes during the FTT. METHODS: Five astronauts participated in the FTT before 10-15 day missions, on landing day (R+0), and one (R+1), six (R+6) and thirty (R+30) days after landing. The RFST consisted of a 2-minute prone rest followed by a 3-minute stand during which heart rate (HR, Holter) and continuous blood pressure (BP, Finometer) were measured. Spectral heart rate variability (HRV) was calculated during the RFST to approximate autonomic function. Statistical analysis was performed with two-factor repeated measures ANOVA. RESULTS: During RFST, HR was higher on R+0 than preflight (p<0.004). This increase in HR persisted on R+1 and R+6 during the stand portion of RFST (p<0.026). BP was well-regulated on all test days. Parasympathetic activity was diminished on R+0 (p=0.035). Sympathovagal balance tended to be affected by spaceflight (main effect, p=0.072), appearing to be slightly elevated during postflight RFST except on R+30. Additionally, analysis of HR during the functional tasks yielded a higher HR on R+0 than preflight during 8 of 11 tasks analyzed, where all tasks had HR return to preflight values by R+30 (p<0.05). CONCLUSION: Spaceflight causes an increase in HR, decrease in parasympathetic activity, and increase in sympathovagal balance, which we confirmed during RFST. These spaceflight-induced changes seen in the RFST, along with the increased postflight HR in most functional tasks, can be used to assess functional performance after short-duration spaceflight.

Arzeno, N. M.↗

Correlations far from equilibrium in charged strongly coupled fluids subjected to a strong magnetic field

Within a holographic model, we calculate the time evolution of 2-point and 1-point correlation functions (of selected operators) within a charged strongly coupled system of many particles. That system is thermalizing from an anisotropic initial charged state far from equilibrium towards equilibrium while subjected to a constant external magnetic field. One main result is that thermalization times for 2-point functions are significantly (approximately three times) larger than those of 1-point functions. Magnetic field and charge amplify this difference, generally increasing thermalization times. However, there is also a competition of scales between charge density, magnetic field, and initial anisotropy, which leads to an array of qualitative changes on the 2- and 1-point functions. There appears to be a strong effect of the medium on 2-point functions at early times, but approximately none at later times. At strong magnetic fields, an apparently universal thermalization time emerges, at which all 2-point functions appear to thermalize regardless of any other scale in the system. Hence, this time scale is referred to as saturation time scale. As extremality is approached in the purely charged case, 2- and 1-point functions appear to equilibrate at infinitely late time. We also compute 2-point functions of charged operators. Our results can be taken to model thermalization in heavy ion collisions, or thermalization in selected condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multivariable quantum signal processing (M-QSP): prophecies of the two-headed oracle

Recent work shows that quantum signal processing (QSP) and its multi-qubit lifted version, quantum singular value transformation (QSVT), unify and improve the presentation of most quantum algorithms. QSP/QSVT characterize the ability, by alternating ansätze, to obliviously transform the singular values of subsystems of unitary matrices by polynomial functions; these algorithms are numerically stable and analytically well-understood. That said, QSP/QSVT require consistent access to a single oracle, saying nothing about computing joint properties of two or more oracles; these can be far cheaper to determine given an ability to pit oracles against one another coherently. This work introduces a corresponding theory of QSP over multiple variables: M-QSP. Surprisingly, despite the non-existence of the fundamental theorem of algebra for multivariable polynomials, there exist necessary and sufficient conditions under which a desired stable multivariable polynomial transformation is possible. Moreover, the classical subroutines used by QSP protocols survive in the multivariable setting for non-obvious reasons, and remain numerically stable and efficient. Up to a well-defined conjecture, we give proof that the family of achievable multivariable transforms is as loosely constrained as could be expected. The unique ability of M-QSP to obliviously approximate joint functions of multiple variables coherently leads to novel speedups incommensurate with those of other quantum algorithms, and provides a bridge from quantum algorithms to algebraic geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reconstructing high-dimensional Hilbert-valued functions via compressed sensing

We present and analyze a novel sparse polynomial technique for approximating high-dimensional Hilbert-valued functions, with application to parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our theoretical framework treats the function approximation problem as a joint sparse recovery problem, where the set of jointly sparse vectors is possibly infinite. To achieve the simultaneous reconstruction of Hilbert-valued functions in both parametric domain and Hilbert space, we propose a novel mixed-norm based `1 regularization method that exploits both energy and sparsity. Our approach requires extensions of concepts such as the restricted isometry and null space properties, allowing us to prove recovery guarantees for sparse Hilbert-valued function reconstruction. We complement the enclosed theory with an algorithm for Hilbert-valued recovery, based on standard forward-backward algorithm, meanwhile establishing its strong convergence in the considered infinite-dimensional setting. Finally, we demonstrate the minimal sample complexity requirements of our approach, relative to other popular methods, with numerical experiments approximating the solutions of high-dimensional parameterized elliptic PDEs.

Dexter, Nick↗

Electronic structure and thermoelectric properties of CoAsSb with post-DFT approaches

Here, we study the electronic structure and thermoelectric properties of recently synthesized CoAsSb. The calculated bandgap becomes more accurate for increasingly-complex electronic structure methods: generalized gradient approximation, hybrid functionals, selfconsistent linearized quasiparticle GW method (LQSGW), and LQSGW with simplified vertex corrections. By equating the bandgaps of each method from a rigid shift of the bands, we evaluate the contributions made to thermoelectric properties beyond the bandgap. In doing so, we evaluate the efficacy of a commonpractice: a rigid shift applied to less-costly electronic structure methods in order to gain some accuracy of the more-costly methods. We find that while the shift made the Seebeck coefficients much closer to one another than from the original electronic structures, there remain differences in the Goldsmid-Sharp (thermoelectric) bandgap between the methods and from the intended electronic bandgap. Additionally, some lasting differences in temperature dependence remain between the methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydrodynamic instability and coalescence of binary neutron stars

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configuration, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of confirgurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, and spins. The influence of the stiffness of the equation of state is also explored by varying the effective polytopic index. Typically, we find that the radial infall velocity just prior to contact is about 10% of the tangential orbital velocity. Once the stability limit is reached, the final evolution only takes another orbit. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibiltiy of mass transfer form one neutron star to the other. We show that stable mass transfer is unlikely except when the mass of one of the components is very small (M is less than or approximately 0.4 solar mass) and the viscosity is high enough to maintain corotation. Otherwise, either the two stars come into contact or the dynamical instability sets in before a Roche limit can be reached.

Lai, Dong↗

Skin and Proximity Effects in Electrodes and Furnace Shells

A review of two-dimensional (2D) analytical models of skin and proximity effects in large industrial furnaces with three electrodes arranged in an equilateral triangle is given. The models cover three different cases: one electrode only, three electrodes where two are approximated by line currents, and induced shell currents where all electrodes are approximated by line currents. The first two models show how the skin and proximity effects depend on electrode material properties and size, and the distance between the electrodes. The third model shows how the strength of the induced shell currents will depend on electrode position and furnace size. These models are compared to numerical studies including distributed electrodes and shell currents. The analytical models are accurate when induced shell currents can be disregarded. However, strong shell currents may have a significant impact on the current distribution within the electrodes. This electrode-shell proximity effect competes with the electrode-electrode proximity effect. Finally, the 2D models have been compared with three-dimensional (3D) case studies of large industrial furnaces. In 3D, the shell currents are significantly smaller than what are predicted by the 2D models, but they are sufficiently strong to cause a significant correction of the electrode current density.

36 MATERIALS SCIENCE↗

Overview of the High Performance Antiproton (HiPAT) Experiment

The annihilation of matter with antimatter represents the highest energy density of any known reaction, producing 10(exp 8) MJ/g, approximately 10 orders of magnitude more energy per unit mass than chemical based combustion. To take the first step towards using this energy for propulsion applications the NASA MSFC Propulsion Research Center (PRC) has initiated a research activity examining the storage of low energy antiprotons. Storage was identified as a key enabling technology since it builds the experience base necessary to understand the handling of antiprotons for virtually all utilization and high-density storage concepts. To address this need, a device referred to as the High Performance Antiproton Trap (HiPAT) is under development at the NASA MSFC PRC. The HiPAT is an electromagnetic system (Penning-Malmberg design) consisting of a 4 Tesla superconductor, a high voltage confinement electrode system (operation up to 20 KV), and an ultra high vacuum test section (operating in the 10(exp -12) torr range). The system was designed to be portable with an ultimate goal of maintaining 10(exp 12) charged particles with a half-life of 18 days. Currently, this system is being experimentally evaluated using normal matter ions which are cheap to produce and relatively easy to handle. These normal ions provide a good indication of overall trap behavior, with the exception of assessing annihilation losses. The ions are produced external to HiPAT using two hydrogen ion sources, with adjustable beam energy and current. Ion are transported in a beam line and controlled through the use of electrostatic optics. These optics serve to both focus and gate the incoming ions, providing microsecond-timed pulses that are dynamically captured by cycling the HiPAT electric containment field like a 'trap door'. The layout of this system more closely simulates the operations expected at an actual antiproton production facility where 'packets' of antiprotons with pulse widths measured in 100's of nanoseconds could be provided. Initial dynamic capture experiments have been performed with both trap and ton source system functioning at approximately 10% of maximum levels, minimizing the potential for extraneous effects. Dynamic trapping techniques demonstrated the successful capture of millions of hydrogen ions with good agreement with the predicted loading based on the timing sequence, trap electric field, and ion beam current. These techniques will be expanded to examine multiple shot capture or 'stacking' to increase the total number of ions stored within HiPAT.

Martin, James J.↗

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗

Probabilistic partition of unity networks for high–dimensional regression problems

We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a fixed-degree polynomial. We present a training strategy that leverages the expectation maximization (EM) algorithm. During the training, we alternate between (i) applying gradient descent to update the DNN coefficients; and (ii) using closed-form formulae derived from the EM algorithm to update the mixture of experts model parameters. Under the probabilistic formulation, step (ii) admits the form of embarrassingly paralleliazable weighted least-squares solves. The PPOU-Nets consistently outperform the baseline fully-connected neural networks of comparable sizes in numerical experiments of various data dimensions. Here, we also explore the proposed model in applications of quantum computing, where the PPOU-Nets act as surrogate models for cost landscapes associated with variational quantum circuits.

97 MATHEMATICS AND COMPUTING↗

Improved Architectures and Training Algorithms for Deep Operator Networks

Operator learning techniques have recently emerged as a powerful tool for learning maps between infinite-dimensional Banach spaces. Trained under appropriate constraints, they can also be effective in learning the solution operator of partial differential equations (PDEs) in an entirely self-supervised manner. In this work we analyze the training dynamics of deep operator networks (DeepONets) through the lens of Neural Tangent Kernel theory, and reveal a bias that favors the approximation of functions with larger magnitudes. To correct this bias we propose to adaptively re-weight the importance of each training example, and demonstrate how this procedure can effectively balance the magnitude of back-propagated gradients during training via gradient descent. We also propose a novel network architecture that is more resilient to vanishing gradient pathologies. Taken together, our developments provide new insights into the training of DeepONets and consistently improve their predictive accuracy by a factor of 10-50x, demonstrated in the challenging setting of learning PDE solution operators in the absence of paired input-output observations.

97 MATHEMATICS AND COMPUTING↗

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gradient-based optimization for regression in the functional tensor-train format

Predictive analysis of complex computational models, such as uncertainty quantification (UQ), must often rely on using an existing database of simulation runs. In this paper we consider the task of performing low-multilinear-rank regression on such a database. Specifically we develop and analyze an efficient gradient computation that enables gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT). We compare our algorithms with 22 other nonparametric and parametric regression methods on 10 real-world data sets and show that for many physical systems, exploiting low-rank structure facilitates efficient construction of surrogate models. Here, we use a number of synthetic functions to build insight into behavior of our algorithms, including the rank adaptation and group-sparsity regularization procedures that we developed to reduce overfitting. Finally we conclude the paper by building a surrogate of a physical model of a propulsion plant on a naval vessel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗