Search NASA⌕ Search

SEARCH · Search NASA

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A Model for Jet-Surface Interaction Noise Using Physically Realizable Upstream Turbulence Conditions

This paper is a continuation of previous work in which a generalized Rapid Distortion Theory (RDT) formulation was used to model low-frequency trailing-edge noise. The research was motivated by proposed next-generation aircraft configurations where the exhaust system is tightly integrated with the airframe. Data from recent experiments at NASA on the interaction between high-Reynolds-number subsonic jet flows and an external flat plate showed that the power spectral density (PSD) of the far-field pressure underwent considerable amplification at low frequencies. For example, at the 90deg observation angle, the low-frequency noise could be as much as 10 dB greater than the jet noise itself. In this paper, we present predictions of the noise generated by the interaction of a rectangular jet with the trailing edge of a semi-infinite flat plate. The calculations are based on a formula for the acoustic spectrum of this noise source derived from an exact formal solution of the linearized Euler equations involving (in this case) one arbitrary convected scalar quantity and a Rayleigh equation Green's function. A low-frequency asymptotic approximation for the Green's function based on a two-dimensional mean flow is used in the calculations along with a physically realizable upstream turbulence spectrum, which includes a finite decorrelation region. Numerical predictions of the sound field, based on three-dimensional RANS solutions to determine the mean flow, turbulent kinetic energy and turbulence length and time scales, for a range of subsonic acoustic Mach number jets and nozzle aspect ratios are compared with experimental data. Comparisons of the RANS results with flow data are also presented for selected cases. We find that a finite decorrelation region in the turbulence spectrum increases the low-frequency algebraic decay (the low frequency "roll-off") of the acoustic spectrum with angular frequency thereby producing much closer agreement with noise data for Strouhal numbers less than 0.1. Secondly, the large-aspect-ratio theory is able to predict the low-frequency amplification due to the jet-edge interaction reasonably well, even for moderate aspect ratio nozzles. We show also that the noise predictions for smaller aspect ratio jets can be fine-tuned using the appropriate RANS-based mean flow and turbulence properties.

Jet↗
NASA NTRSDate not supplied

Phased Array Characterization of Slat Noise Radiation from a High-Lift Common Research Model

A test campaign was conducted in the open-jet test section of the NASA Langley 14- by 22-Foot Subsonic Tunnel on a 10%-scale, semispan version of a High-Lift Common Research Model (CRM-HL) incorporating a leading-edge slat, trailing-edge flap and removable high-fidelity main landing gear. Far field acoustic measurements were obtained on baseline and acoustically treated model configurations using a traversing 97-microphone phased array that viewed the pressure side of the airframe. A modified version of the DAMAS deconvolution method incorporating corrections for shear layer acoustic wave decorrelation was employed to determine the locations and strengths of relevant noise sources along the span of the slat and in the vicinity of the flap edges and landing gear. The main goal of the test campaign was to use the array data to evaluate the effectiveness of various slat noise reduction concepts. It was found that slat-gap fillers produced significant broadband noise reduction for all Mach numbers and airfoil angles of attack that were investigated. Mach number scaling revealed an approximate velocity to the sixth power relationship for the slat noise at higher frequencies. Directivity measurements obtained from integration of DAMAS pressure-squared values over defined geometric zones around the model showed that baseline slat noise radiation approximated a dipole as a function of polar emission angle, becoming approximately omnidirectional with the introduction of the slat-gap filler. The test campaign clearly demonstrated the noise reduction benefits that can be obtained by applying appropriate treatments to leading edge slats on commercial transport-class aircraft.

William M Humphreys, Jr.↗
NASA NTRSDate not supplied

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗
DOE OSTI2022

Conditional probability density functional theory

Here we present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this paper we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Practical algorithms for multivariate rational approximation

We present two approaches for computing rational approximations to multivariate functions, motivated by their effectiveness as surrogate models for high-energy physics (HEP) applications. Our first approach builds on the Stieltjes process to efficiently and robustly compute the coefficients of the rational approximation. Our second approach is based on an optimization formulation that allows us to include structural constraints on the rational approximation (in particular, constraints demanding the absence of singularities), resulting in a semi-infinite optimization problem that we solve using an outer approximation approach. We present results for synthetic and real-life HEP data, and we compare the approximation quality of our approaches with that of traditional polynomial approximations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Probabilistic and deterministic aspects of linear estimation in geodesy

Recent advances in observational techniques related to geodetic work (VLBI, laser ranging) make it imperative that more consideration should be given to modeling problems. Uncertainties in the effect of atmospheric refraction, polar motion and precession-nutation parameters, cannot be dispensed with in the context of centimeter level geodesy. Even physical processes that have generally been previously altogether neglected (station motions) must now be taken into consideration. The problem of modeling functions of time or space, or at least their values at observation points (epochs) is explored. When the nature of the function to be modeled is unknown. The need to include a limited number of terms and to a priori decide upon a specific form may result in a representation which fails to sufficiently approximate the unknown function. An alternative approach of increasing application is the modeling of unknown functions as stochastic processes.

Dermanis, A.↗
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗
DOE OSTI2018

Survival of the most transferable at the top of Jacob’s ladder: Defining and testing the ω B97M(2) double hybrid density functional

A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resu|ting functional, ωB97M(2) is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 datq points clearly surpasses that of all of the tested density funstionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

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

Time correlation functions in classical dense fluids

The velocity autocorrelation functions and memory functions of dense classical fluids may be directly obtained from the static radial distribution function g(r) in an approximate way. Following the Mori projection operator formalism, the memory functions may be related to the fluctuating force correlation. At low densities, these functions may be evaluated by following the trajectories of particle pairs in the interatomic potential. At higher densities, the force correlation functions can be evaluated approximately from particle pair trajectories via the potential of the mean force. The results for argon fluids at various densities and temperatures agree satisfactorily with the molecular dynamics and the Enskog values. The decrease of the diffusion coefficient with density is partly due to the nature of g(r) which reflects the stronger clustering of atoms at higher densities.

Tsang, T.↗
DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗
DOE OSTI2020

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Smooth indentation of orthotropic beams

An approximate Green's function for surface displacements in an orthotropic beam is derived as the superposition of the half-plane solution for displacements and the beam theory deflections. The Green's function is used to formulate the integral equation for the problem of smooth contact between a rigid-cylinder and a simply supported orthotropic beam. The integral equation is solved using a least-squares approximation procedure. The contact stress distribution is presented for three materials with different degrees of orthotropy. Numerical results are given for contact force-contact length and contact force-indentation relations. The effects of curvature of the deflected beam, shear deformation and length-to-thickness ratio of the beam on the contact behavior are discussed.

Sankar, B. V.↗
DOE OSTI2018

Construction and Investigation of the Third-Order Approximation to the Solution of the Heat-Conduction Equation for Thin Shallow Shells by Using Legendre Polynomials in the Case of Stationary Heat Exchange

By using the N th order approximation of temperature function and its first derivative by Legendre polynomials, we construct the solution of problem of heat conduction for a thin-walled shallow isotropic shell and deduce the system of resolving equations for the N th approximation. For the first and third approximations, we solve this problem for the case of concentrated heat sources. In the third approximation, we also construct the plots of the dependences of temperature on the distance to the heat source and on the curvature of the shell under the conditions of symmetric or asymmetric stationary heat exchange.

42 ENGINEERING↗
DOE OSTI2022

Augmented Adam-Gibbs model for glass melt viscosity and configuration entropy as functions of temperature and composition

As the temperature of glass melt increases, its structure approaches the state of a simple liquid while the configuration entropy approaches a maximum value. We describe this gradual change using a power law function of inverse temperature. The Adam-Gibbs model for glass viscosity as a function of temperature and glass composition augmented in this way is greatly simplified when applied to common glass families occupying moderate composition regions, such as float glass or nuclear waste glasses, on which properties can be approximated as linear functions of composition. The parsimonious model thus obtained is preferable for use in optimizing glass formulation and mathematical modeling of glass melting and forming. For multicomponent glasses with N viscosity-affecting components, the augmented Adam-Gibbs model requires 2N + 3 adjustable parameters. The model efficacy is demonstrated by fitting the model to a viscosity-temperature-composition dataset for low-activity nuclear waste glasses.

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

Effect of 6 degrees head-down tilt on cardiopulmonary function: comparison with microgravity

BACKGROUND: Head-down tilt (HDT) of 6 degrees is a commonly used model of weightlessness, but there are few comparisons with actual microgravity. HYPOTHESIS: Our study was designed to prove that the changes in cardiopulmonary function seen in HDT would be similar to those seen in microgravity. METHODS: We compared measurements of cardiovascular and pulmonary function from three separate spaceflights of 14 to 17 d duration, with data collected during a 17-d period of HDT. RESULTS: HDT proved a good model of the cardiovascular response to microgravity, resulting in increases in cardiac output and stroke volume of a similar magnitude to those seen in microgravity, with a concomitant reduction in heart rate. By contrast, HDT was a poor model of the effects of microgravity on pulmonary ventilation and gas exchange. CONCLUSION: Pulmonary function in HDT approximated the changes seen in the 1-G supine posture, while in microgravity this was much closer to that seen in the 1-G upright position. The differences probably reflect the fact that changes in cardiovascular function result primarily from fluid shifts within the entire body, whereas changes in pulmonary ventilation are primarily a result of mechanical influences on the lung and chest and abdominal wall.

NASA Discipline Cardiopulmonary↗
NASA NTRS2004

Miniature Robotic Spacecraft for Inspecting Other Spacecraft

A report discusses the Miniature Autonomous Extravehicular Robotic Camera (Mini AERCam)-- a compact robotic spacecraft intended to be released from a larger spacecraft for exterior visual inspection of the larger spacecraft. The Mini AERCam is a successor to the AERCam Sprint -- a prior miniature robotic inspection spacecraft that was demonstrated in a space-shuttle flight experiment in 1997. The prototype of the Mini AERCam is a demonstration unit having approximately the form and function of a flight system. The Mini AERCam is approximately spherical with a diameter of about 7.5 in. (.19 cm) and a weight of about 10 lb (.4.5 kg), yet it has significant additional capabilities, relative to the 14-in. (36-cm), 35-lb (16-kg) AERCam Sprint. The Mini AERCam includes miniaturized avionics, instrumentation, communications, navigation, imaging, power, and propulsion subsystems, including two digital video cameras and a high-resolution still camera. The Mini AERCam is designed for either remote piloting or supervised autonomous operations, including station keeping and point-to-point maneuvering. The prototype has been tested on an air-bearing table and in a hardware-in-the-loop orbital simulation of the dynamics of maneuvering in proximity to the International Space Station.

Fredrickson, Steven↗