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 1,297 records · Page 72

Approximate formulation of redistribution in the Ly-alpha, Ly-beta, H-alpha system

Simple approximate formulas are given for the coupled redistribution of Ly-alpha, Ly-beta, and H-alpha, by using well-defined approximations to an essentially exact formulation. These formulas incorporate all the essential physics including Raman scattering, lower state radiative decay, and correlated terms representing emission during a collision which must be retained in order that the emission coefficients are properly behaved in the line wings. Approximate expressions for the appropriate line broadening parameters are collected. Finally, practical expressions for the source functions are given. These are formulated through newly introduced nonimpact redistribution functions, which are shown to be reasonably approximated by existing (ordinary and generalized) redistribution functions.

Cooper, J.↗

A simple parameterization for quality factor as a function of linear energy transfer

The paper presents a simple analytic approximation of the radiation quality factor (Q) as a function of linear energy transfer for use in radiation protection calculations. The paper also presents estimated quality factors in water for protons over a broad range of incident energies. It is shown that the quality factors are less than unity for all proton energies greater than 13 MeV.

Townsend, Lawrence W.↗

Extending Parsimonious Bayesian Inference

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

97 MATHEMATICS AND COMPUTING↗

The response function of modulated grid Faraday cup plasma instruments

Modulated grid Faraday cup plasma analyzers are a very useful tool for making in situ measurements of space plasmas. One of their great attributes is that their simplicity permits their angular response function to be calculated theoretically. An expression is derived for this response function by computing the trajectories of the charged particles inside the cup. The Voyager Plasma Science (PLS) experiment is used as a specific example. Two approximations to the rigorous response function useful for data analysis are discussed. The theoretical formulas were tested by multi-sensor analysis of solar wind data. The tests indicate that the formulas represent the true cup response function for all angles of incidence with a maximum error of only a few percent.

Barnett, A.↗

Seasonal Variability of Middle Latitude Ozone in the Lowermost Stratosphere Derived from Probability Distribution Functions

We present a study of the distribution of ozone in the lowermost stratosphere with the goal of characterizing the observed variability. The air in the lowermost stratosphere is divided into two population groups based on Ertel's potential vorticity at 300 hPa. High (low) potential vorticity at 300 hPa indicates that the tropopause is low (high), and the identification of these two groups is made to account for the dynamic variability. Conditional probability distribution functions are used to define the statistics of the ozone distribution from both observations and a three-dimensional model simulation using winds from the Goddard Earth Observing System Data Assimilation System for transport. Ozone data sets include ozonesonde observations from northern midlatitude stations (1991-96) and midlatitude observations made by the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite (UARS) (1994- 1998). The conditional probability distribution functions are calculated at a series of potential temperature surfaces spanning the domain from the midlatitude tropopause to surfaces higher than the mean tropical tropopause (approximately 380K). The probability distribution functions are similar for the two data sources, despite differences in horizontal and vertical resolution and spatial and temporal sampling. Comparisons with the model demonstrate that the model maintains a mix of air in the lowermost stratosphere similar to the observations. The model also simulates a realistic annual cycle. Results show that during summer, much of the observed variability is explained by the height of the tropopause. During the winter and spring, when the tropopause fluctuations are larger, less of the variability is explained by tropopause height. This suggests that more mixing occurs during these seasons. During all seasons, there is a transition zone near the tropopause that contains air characteristic of both the troposphere and the stratosphere. The relevance of the results to the assessment of the environmental impact of aircraft effluence is also discussed.

Rood, Richard B.↗

The Unsteady Lift of a Finite Wing

Unsteady lift function for wings of finite aspect ratio have been calculated by approximate methods involving corrections of the aerodynamic inertia and of the angle of the infinite wing. The starting lift of the finite wing is found to be only slightly less than that of the infinite wing; whereas the final lift may be considerably less. The calculations indicate that the distribution of lift near the start is similar to the final distribution. Both the indicia and the oscillating lift functions are given. Approximate operational equivalents of the functions have been devised to facilitate the calculation of lift under various conditions of motion.

Jones, Robert T↗

Small-angle approximation to the transfer of narrow laser beams in anisotropic scattering media

The broadening and the signal power detected of a laser beam traversing an anisotropic scattering medium were examined using the small-angle approximation to the radiative transfer equation in which photons suffering large-angle deflections are neglected. To obtain tractable answers, simple Gaussian and non-Gaussian functions for the scattering phase functions are assumed. Two other approximate approaches employed in the field to further simplify the small-angle approximation solutions are described, and the results obtained by one of them are compared with those obtained using small-angle approximation. An exact method for obtaining the contribution of each higher order scattering to the radiance field is examined but no results are presented.

Box, M. A.↗

Tomography of the rho meson in the QCD instanton vacuum: Transverse momentum dependent parton distribution functions

We analyze the rho meson unpolarized and polarized transverse momentum dependent parton distribution functions (TMDPDFs) in the instanton liquid model (ILM). The corresponding TMDs in ILM are approximated by a constituent quark beam function in the leading Fock state multiplied by a rapidity-dependent factor resulting from the staple-shaped Wilson lines, for fixed longitudinal momentum, transverse separation, and rapidity. At the resolution of the ILM, all of the rho meson TMDs are symmetric in parton 𝑥 for fixed transverse momentum, and Gaussian-like in the transverse momentum for fixed parton 𝑥. The latter is a direct consequence of the profiling of the quark zero modes in the ILM. The evolved TMDs at higher rapidity using the Collins-Soper kernel, and higher resolution using the renormalization group, show substantial skewness towards low parton 𝑥.

Gluons↗

Local thermodynamic mapping for effective liquid density-functional theory

The structural-mapping approximation introduced by Lutsko and Baus (1990) in the generalized effective-liquid approximation is extended to include a local thermodynamic mapping based on a spatially dependent effective density for approximating the solid phase in terms of the uniform liquid. This latter approximation, called the local generalized effective-liquid approximation (LGELA) yields excellent predictions for the free energy of hard-sphere solids and for the conditions of coexistence of a hard-sphere fcc solid with a liquid. Moreover, the predicted free energy remains single valued for calculations with more loosely packed crystalline structures, such as the diamond lattice. The spatial dependence of the weighted density makes the LGELA useful in the study of inhomogeneous solids.

Kyrlidis, Agathagelos↗

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗

Galerkin finite element method for time-fractional stochastic diffusion equations

In this paper, Galerkin finite element method for solving the time-fractional stochastic diffusion equations with multiplicative noise is proposed and investigated. The pathwise regularity properties of solutions to the semidiscrete Galerkin approximations are demonstrated and the convergence of optimal rates are derived. And also we construct the fully discrete scheme which is based on the approximations of the Mittag–Leffler function and analyze the error estimates of convergence in L{sub 2}-norm space. Finally, numerical results are conducted to confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Hybrid functional pseudopotentials

The consistency between the exchange-correlation functional used in pseudopotential construction and in the actual density functional theory calculation is essential for the accurate prediction of fundamental properties of materials. However, routine hybrid density functional calculations at present still rely on generalized gradient approximation pseudopotentials due to the lack of hybrid functional pseudopotentials. Here, we present a scheme for generating hybrid functional pseudopotentials, and we analyze the importance of pseudopotential density functional consistency for hybrid functionals. For the PBE0 hybrid functional, we benchmark our pseudopotentials for structural parameters and fundamental electronic gaps of the Gaussian-2 (G2) molecular dataset and some simple solids. Our results demonstrate that using our PBE0 pseudopotentials in PBE0 calculations improves agreement with respect to all-electron calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Aperture impedance of flared horns

The method of moments is often used when solving for the mutual coupling in arrays of aperture antennas. For elements that are waveguides or gradually flared horns the aperture fields can be approximated by a finite sum of waveguide modal functions. To solve for the flared horn case an approximation for the aperture impedances of the modes in the horn is needed. The WKB approach can be used to find these impedances, but this technique has an important limitation. It is known to fail in the vicinity of its turning points. The turning point is the cutoff point of the mode being considered. To overcome this limitation a different technique, the spherical mode approach, is discussed. This approach has no cutoff problems and works well for conical and pyramidal horns. Comparisons between the impedances and the resulting dominant mode reflection coefficient found using the two techniques are presented to illustrate this point.

Silvestro, J. W.↗

Simplifying activations with linear approximations in neural networks

A key step in Neural Networks is activation. Among the different types of activation functions, sigmoid, tanh, and others involve the usage of exponents for calculation. From a hardware perspective, exponential implementation implies the usage of Taylor series or repeated methods involving many addition, multiplication, and division steps, and as a result are power-hungry and consume many clock cycles. We implement a piecewise linear approximation of the sigmoid function as a replacement for standard sigmoid activation libraries. This approach provides a practical alternative by leveraging piecewise segmentation, which simplifies hardware implementation and improves computational efficiency. In this paper, we detail piecewise functions that can be implemented using linear approximations and their implications for overall model accuracy and performance gain. Our results show that for the DenseNet, ResNet, and GoogLeNet architectures, the piecewise linear approximation of the sigmoid function provides faster execution times compared to the standard TensorFlow sigmoid implementation while maintaining comparable accuracy. Specifically, for MNIST with DenseNet, accuracy reaches 99.91% (Piecewise) vs. 99.97% (Base) with up to 1.31x speedup in execution time. For CIFAR-10 with DenseNet, accuracy improves to 98.97% (Piecewise) vs. 99.40% (Base) while achieving 1.24x faster execution. Similarly, for CIFAR-100 with DenseNet, the accuracy is 97.93% (Piecewise) vs. 98.39% (Base), with a 1.18x execution time reduction. These results confirm the proposed method’s capability to efficiently process large-scale datasets and computationally demanding tasks, offering a practical means to accelerate deep learning models, including LSTMs, without compromising accuracy.

Activation function↗

Electronic Structure and Linhard Function of Iron-Based Superconductor La{sub 0.5−x}Na{sub 0.5 + x}Fe{sub 2}As{sub 2}

We perform theoretical studies of the iron-based superconductor La{sub 0.5−x}Na{sub 0.5 + x}Fe{sub 2}As{sub 2}. Using density functional calculation with virtual crystal approximation, band structures, density of states, and Fermi surfaces are obtained for different La/Na composition by varying x value. Furthermore, the nesting property is also investigated using the Linhard function. Our main findings are as follows: For undoped (x = 0.0) case, the magnetic and electronic structure share great similarities with extensively studied BaFe{sub 2}As{sub 2}. Upon hole doping (x = 0.3), the Q = (π,π) nesting property are preserved. However, upon electron doping (x = -0.3), the Fermi surfaces undergo a dramastic change into 3D and complex structure. We find this change is attributed to the comprehensive influence of band filling and chemical pressure of La{sup 3+}. We propose the asymmetric change of Fermi surfaces upon doping may in some sense explain the electron-hole asymmetry in phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural tailoring of advanced turboprops

A computer program has been developed for the performance of numerical optimizations of highly swept propfan blades by minimizing an objective function that is defined either as direct operating cost or the aeroelastic difference between a blade and its scaled model. Three component analysis categories are employed: an optimization algorithm, approximate analysis procedures for objective function and constraint evaluation, and refined analysis procedures for optimum design validation. The analyses conducted by the program encompass aerodynamic efficiency evaluation, finite element stress and vibration analysis, acoustics, flutter, and forced response life prediction.

Brown, K. W.↗

Flexure-torsion behavior of prismatic beams. I - Section properties via power series

The behavior of a tip-loaded cantilever beam with an arbitrary cross section is studied using Saint-Venant's semi-inverse method along with a power series solution for the out-of-plane flexure and torsion warping functions. The power series coefficients are determined by solving a set of variationally derived linear algebraic equations. For complex cross sections, the calculated coefficients represented a 'best-fit approximation' to the exact warping function. The resulting warping functions are used to determine the cross-sectional properties (torsion constant, shear correction factors, shear deformation coefficients, and shear center location). A new linear relation is developed for locating the shear center, where the twist rate is zero about the line of shear centers. Moreover, the kinematic relations for a new fully compatible one-dimensional beam theory are developed. Numerical results are presented first to verify the approach and second to provide section data on NACA four-series airfoils not currently found in the literature.

Kosmatka, J. B.↗