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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Space Shuttle Debris Impact Tool Assessment Using the Modern Design of Experiments

Complex computer codes are used to estimate thermal and structural reentry loads on the Shuttle Orbiter induced by ice and foam debris impact during ascent. Such debris can create cavities in the Shuttle Thermal Protection System. The sizes and shapes of these cavities are approximated to accommodate a code limitation that requires simple "shoebox" geometries to describe the cavities -- rectangular areas and planar walls that are at constant angles with respect to vertical. These approximations induce uncertainty in the code results. The Modern Design of Experiments (MDOE) has recently been applied to develop a series of resource-minimal computational experiments designed to generate low-order polynomial graduating functions to approximate the more complex underlying codes. These polynomial functions were then used to propagate cavity geometry errors to estimate the uncertainty they induce in the reentry load calculations performed by the underlying code. This paper describes a methodological study focused on evaluating the application of MDOE to future operational codes in a rapid and low-cost way to assess the effects of cavity geometry uncertainty.

DeLoach, Richard↗
DOE OSTI2024

Chapter 7: Learning Stable Local Volt/Var Controllers in Distribution Grids

This chapter describes a framework to synthesize provably stable local Volt/Var controllers for distributed energy resources (DERs) in power distribution grids (DGs). The goal is to control the reactive power injections of DERs to improve the system performance as quantified by a generic optimal reactive power flow (ORPF) problem. To achieve this, we jointly design for each DER the control function, which prescribes the reactive power update rule, and the equilibrium function, which approximates the ORPF solutions from local measurements of voltages and powers. We provide conditions on the equilibrium functions and the control parameters ensuring the stability of the closed-loop system. In particular, we discuss the trade-offs between each set of conditions accounting for practical considerations, like fully exploiting the DERs' generation capabilities and reducing the optimality gap. These conditions are then translated into learning constraints on the neural networks' parameters that are enforced in the training phase. We validate our framework with numerical simulations on the IEEE 37-bus network and through a comparison with an optimized version of standard piece wise linear control rules.

closed-loop asymptotic stability↗
NASA NTRS1990

Causal compensated perturbations in cosmology

A theoretical framework is developed to calculate linear perturbations in the gravitational and matter fields which arise causally in response to the presence of stiff matter sources in a FRW cosmology. It is shown that, in order to satisfy energy and momentum conservation, the gravitational fields of the source must be compensated by perturbations in the matter and gravitational fields, and the role of such compensation in containing the initial inhomogeneities in their subsequent evolution is discussed. A complete formal solution is derived in terms of Green functions for the perturbations produced by an arbitrary source in a flat universe containing cold dark matter. Approximate Green function solutions are derived for the late-time density perturbations and late-time gravitational waves in a universe containing a radiation fluid. A cosmological energy-momentum pseudotensor is defined to clarify the nature of energy and momentum conservation in the expanding universe.

Veeraraghavan, Shoba↗
DOE OSTI2019

First-principles study of the infrared spectrum in liquid water from a systematically improved description of H-bond network

An accurate ab initio theory of the H-bond structure of liquid water requires a high-level exchange correlation approximation from density functional theory. Based on the liquid structures modeled by ab initio molecular dynamics by using maximally localized Wannier functions as a basis, we study the infrared spectrum of water within the canonical ensemble. In particular, we employ both the Perdew-Burke-Ernzerhof (PBE) functional within the generalized gradient approximation (GGA) and the state-of-the-art meta-GGA level approximation provided by the strongly constrained and appropriately normed (SCAN) functional. In this work, we demonstrate that the SCAN functional improves not only the water structure but also the theoretical infrared spectrum of water. Our analyses show that the improvement in the stretching and bending bands can be mainly attributed to better descriptions of directional H bonding and the covalency at the inter- and intramolecular levels, respectively. On the other hand, better agreements in libration and hindered translation bands are due to the improved dynamics of the H-bond network enabled by a less structured liquid in the experimental direction. The spectrum predicted by SCAN shows much better agreement with experimental data than the conventionally widely adopted PBE functional at the GGA level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

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

Actuator and aerodynamic modeling for high-angle-of-attack aeroservoelasticity

Accurate prediction of airframe/actuation coupling is required by the imposing demands of modern flight control systems. In particular, for agility enhancement at high angle of attack and low dynamic pressure, structural integration characteristics such as hinge moments, effective actuator stiffness, and airframe/control surface damping can have a significant effect on stability predictions. Actuator responses are customarily represented with low-order transfer functions matched to actuator test data, and control surface stiffness is often modeled as a linear spring. The inclusion of the physical properties of actuation and its installation on the airframe is therefore addressed using detailed actuator models which consider the physical, electrical, and mechanical elements of actuation. The aeroservoelastic analysis procedure is described in which the actuators are modeled as detailed high-order transfer functions and as approximate low-order transfer functions. The impacts of unsteady aerodynamic modeling on aeroservoelastic stability are also investigated by varying the order of approximation, or number of aerodynamic lag states, in the analysis. Test data from a thrust-vectoring configuration of an F/A-l8 aircraft are compared to predictions to determine the effects on accuracy as a function of modeling complexity.

Brenner, Martin J.↗
DOE OSTI2018

Exploring non-adiabatic approximations to the exchange–correlation functional of TDDFT

A decomposition of the exact exchange–correlation potential of time-dependent density functional theory into an interaction component and a kinetic component offers a new starting point for non-adiabatic approximations. The components are expressed in terms of the exchange–correlation hole and the difference between the one-body density matrix of the interacting and Kohn–Sham systems, which must be approximated in terms of quantities accessible from the Kohn–Sham evolution. We explore several preliminary approximations, evaluate their fulfillment of known exact conditions, and test their performance on simple model systems for which available exact solutions indicate the significance of going beyond the adiabatic approximation.

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

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Actuator and aerodynamic modeling for high-angle-of-attack aeroservoelasticity

Accurate prediction of airframe/actuation coupling is required by the imposing demands of modern flight control systems. In particular, for agility enhancement at high angle of attack and low dynamic pressure, structural integration characteristics such as hinge moments, effective actuator stiffness, and airframe/control surface damping can have a significant effect on stability predictions. Actuator responses are customarily represented with low-order transfer functions matched to actuator test data, and control surface stiffness is often modeled as a linear spring. The inclusion of the physical properties of actuation and its installation on the airframe is therefore addressed in this paper using detailed actuator models which consider the physical, electrical, and mechanical elements of actuation. The aeroservoelastic analysis procedure is described in which the actuators are modeled as detailed high-order transfer functions and as approximate low-order transfer functions. The impacts of unsteady aerodynamic modeling on aeroservoelastic stability are also investigated in this paper by varying the order of approximation, or number of aerodynamic lag states, in the analysis. Test data from a thrust-vectoring configuration of an F/A-18 aircraft are compared to predictions to determine the effects on accuracy as a function of modeling complexity.

Brenner, Martin J.↗
DOE OSTI2019

Scaling laws and stability of nano-sized defect clusters in niobium via atomistic simulations and statistical analysis

The predictions of scaling laws for the structure and properties of defect clusters are generally limited to small defect clusters in their ground-state configurations. We investigated the size and geometrical configuration dependence of nano-sized defect clusters in niobium (Nb) using molecular dynamics. We studied the structure and stability of large clusters of size up to fifty defects for vacancies and one hundred defects for interstitials, as well as the role of helium and metastable configurations on the stability of these clusters. We compared three different interatomic potentials in order to determine the relative stability of these clusters as a function of their size and geometrical configurations. Additionally, we conducted a statistical analysis to predict the formation and binding energies of interstitial clusters as a function of both their size and configuration. We find that the size dependence of vacancy and interstitial clusters can be approximated by functional forms that account for bulk and surface effects as well as some considerations of elastic interactions. We also find that helium and metastable configurations can make vacancy and interstitial clusters thermally stable depending on the configuration. Our parameterized functional forms for the formation and binding energies are valid for a very broad range of defect sizes and configurations making it possible to be used directly in a coarse-grained modeling strategy such as Monte Carlo, cluster dynamics or dislocation dynamics which look at defect accumulation and evolution in microstructures.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗
DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Soft mean spherical approximation for dusty plasma liquids: Level of accuracy and analytic expressions

The soft mean spherical approximation is employed for the study of the thermodynamics of dusty plasma liquids, the latter treated as Yukawa one-component plasmas. Within this integral theory method, the only input necessary for the calculation of the reduced excess energy stems from the solution of a single non-linear algebraic equation. Consequently, thermodynamic quantities can be routinely computed without the need to determine the pair correlation function or the structure factor. The level of accuracy of the approach is quantified after an extensive comparison with numerical simulation results. The approach is solved over a million times with input spanning the whole parameter space and reliable analytic expressions are obtained for the basic thermodynamic quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Fragment-based {sup 13}C nuclear magnetic resonance chemical shift predictions in molecular crystals: An alternative to planewave methods

We assess the quality of fragment-based ab initio isotropic {sup 13}C chemical shift predictions for a collection of 25 molecular crystals with eight different density functionals. We explore the relative performance of cluster, two-body fragment, combined cluster/fragment, and the planewave gauge-including projector augmented wave (GIPAW) models relative to experiment. When electrostatic embedding is employed to capture many-body polarization effects, the simple and computationally inexpensive two-body fragment model predicts both isotropic {sup 13}C chemical shifts and the chemical shielding tensors as well as both cluster models and the GIPAW approach. Unlike the GIPAW approach, hybrid density functionals can be used readily in a fragment model, and all four hybrid functionals tested here (PBE0, B3LYP, B3PW91, and B97-2) predict chemical shifts in noticeably better agreement with experiment than the four generalized gradient approximation (GGA) functionals considered (PBE, OPBE, BLYP, and BP86). A set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided based on these benchmark calculations. Statistical cross-validation procedures are used to demonstrate the robustness of these fits.

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

Dynamical systems and nonlinear transient rheology of the far-from-equilibrium Bjorken flow

In relativistic kinetic theory, the one-particle distribution function is approximated by an asymptotic perturbative power series in the Knudsen number which is divergent. For the Bjorken flow, we expand the distribution function in terms of its moments and study their nonlinear evolution equations. The resulting coupled dynamical system can be solved for each moment consistently using a multiparameter transseries which makes the constitutive relations inherit the same structure. A new nonperturbative dynamical renormalization scheme is born out of this formalism that goes beyond the linear response theory.We show that there is a Lyapunov function, also known as dynamical potential, which is, in general, a function of the moments and time satisfying Lyapunov stability conditions along renormalization group flows connected to the asymptotic hydrodynamic fixed point. As a result, the transport coefficients get dynamically renormalized at every order in the time-dependent perturbative expansion by receiving nonperturbative corrections present in the transseries. The connection between the integration constants and the UV data is discussed using the language of dynamical systems. Furthermore, we show that the first dissipative correction in the Knudsen number to the distribution function is not only determined by the known effective shear viscous term but also a new high-energy nonhydrodynamic mode. It is demonstrated that the survival of this new mode is intrinsically related to the nonlinear mode-to-mode coupling with the shear viscous term. Finally, we comment on some possible phenomenological applications of the proposed nonhydrodynamic transport theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1976

Atmospheric refraction errors in laser ranging systems

The effects of horizontal refractivity gradients on the accuracy of laser ranging systems were investigated by ray tracing through three dimensional refractivity profiles. The profiles were generated by performing a multiple regression on measurements from seven or eight radiosondes, using a refractivity model which provided for both linear and quadratic variations in the horizontal direction. The range correction due to horizontal gradients was found to be an approximately sinusoidal function of azimuth having a minimum near 0 deg azimuth and a maximum near 180 deg azimuth. The peak to peak variation was approximately 5 centimeters at 10 deg elevation and decreased to less than 1 millimeter at 80 deg elevation.

Gardner, C. S.↗
DOE OSTI2015

Image reconstruction in higher dimensions: myocardial perfusion imaging of tracer dynamics with cardiac motion due to deformation and respiration

Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. We find these results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Finally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.

60 APPLIED LIFE SCIENCES↗