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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Application of Linear Additive Conditions for Near-Infrared Diffuse Reflectance Absorption Spectroscopy

Determining the homogeneity of material mixing in real time during product processing is critical for quality control. According to the Kubelka–Munk (K-M) function of diffuse reflectance absorption spectrum, absorbance (A) is approximately linear with the content of the components when the sample scattering coefficient (S) is in a certain range. The S is determined by the particle size of powder samples. Therefore, this study determined particle size ranges that satisfy linear additivity in near-infrared diffuse reflectance spectroscopy (NIRDRS). Thus, the proposed NIRDRS analysis technique can be used to determine the homogeneity of material mixes or analyze the percentages of the components in the mixture. In this study, vitamin B3 and vitamin C were used for preparing mixed samples with varying percentages. The experimental results revealed that linear additivity is satisfied when the powder particle size is in the range of less than 280, 280–450, and 450–900 μm. When the confidence level is 0.01, the actual mixed spectra are not significantly different from the “simulated mixed spectra” constructed by linear addition, with their relative deviations less than 1.08%. The absolute errors of the actual and analytic percentages were within 2.98% for each component in the mixtures. The above conclusions also hold for sorghum, which has a complex material composition. Statistical models cannot analyze the percentages of components in the mixture. In contrast, linear addition and direct calibration approach avoids the use of a large number of samples for statistical modeling and analyze the percentages of mixed samples. Meanwhile, it can be used to discriminate and analyze the material mixing uniformity by building a mechanistic model.

Feng, Zhiyue↗
DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Competing Magnetic Interactions and Field-Induced Metamagnetic Transition in Highly Crystalline Phase-Tunable Iron Oxide Nanorods

The inherent existence of multi phases in iron oxide nanostructures highlights the significance of them being investigated deliberately to understand and possibly control the phases. Here, the effects of annealing at 250 °C with a variable duration on the bulk magnetic and structural properties of high aspect ratio biphase iron oxide nanorods with ferrimagnetic Fe3O4 and antiferromagnetic α-Fe2O3 are explored. Increasing annealing time under a free flow of oxygen enhanced the α-Fe2O3 volume fraction and improved the crystallinity of the Fe3O4 phase, identified in changes in the magnetization as a function of annealing time. A critical annealing time of approximately 3 h maximized the presence of both phases, as observed via an enhancement in the magnetization and an interfacial pinning effect. This is attributed to disordered spins separating the magnetically distinct phases which tend to align with the application of a magnetic field at high temperatures. The increased antiferromagnetic phase can be distinguished due to the field-induced metamagnetic transitions observed in structures annealed for more than 3 h and was especially prominent in the 9 h annealed sample. Our controlled study in determining the changes in volume fractions with annealing time will enable precise control over phase tunability in iron oxide nanorods, allowing custom-made phase volume fractions in different applications ranging from spintronics to biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2017

SUNSPOT AND STARSPOT LIFETIMES IN A TURBULENT EROSION MODEL

Quantitative models of sunspot and starspot decay predict the timescale of magnetic diffusion and may yield important constraints in stellar dynamo models. Motivated by recent measurements of starspot lifetimes, we investigate the disintegration of a magnetic flux tube by nonlinear diffusion. Previous theoretical studies are extended by considering two physically motivated functional forms for the nonlinear diffusion coefficient D : an inverse power-law dependence D ∝ B {sup −ν} and a step-function dependence of D on the magnetic field magnitude B . Analytical self-similar solutions are presented for the power-law case, including solutions exhibiting “super fast” diffusion. For the step-function case, the heat-balance integral method yields approximate solutions, valid for moderately suppressed diffusion in the spot. The accuracy of the resulting solutions is confirmed numerically, using a method which provides an accurate description of long-time evolution by imposing boundary conditions at infinite distance from the spot. The new models may allow insight into the differences and similarities between sunspots and starspots.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Rosseland and Flux Mean Opacities for Compton Scattering

Rosseland mean opacity plays an important role in theories of stellar evolution and X-ray burst models. In the high-temperature regime, when most of the gas is completely ionized, the opacity is dominated by Compton scattering. Our aim here is to critically evaluate previous works on this subject and to compute the exact Rosseland mean opacity for Compton scattering over a broad range of temperature and electron degeneracy parameter. We use relativistic kinetic equations for Compton scattering and compute the photon mean free path as a function of photon energy by solving the corresponding integral equation in the diffusion limit. As a byproduct we also demonstrate the way to compute photon redistribution functions in the case of degenerate electrons. We then compute the Rosseland mean opacity as a function of temperature and electron degeneracy and present useful approximate expressions. We compare our results to previous calculations and find a significant difference in the low-temperature regime and strong degeneracy. We then proceed to compute the flux mean opacity in both free-streaming and diffusion approximations, and show that the latter is nearly identical to the Rosseland mean opacity. We also provide a simple way to account for the true absorption in evaluating the Rosseland and flux mean opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Parameter distribution mapping in a gas turbine engine

A profile map of parameter values in a region is determined based on average values along linear paths through the region. In one example, a temperature map of a region of a gas turbine engine is created based on time-of-flight measurements from acoustic transceiver pairs arranged circumferentially around the region. A speed of sound for each transceiver pair is determined based on the time-of-flight measurements. An average temperature along each path is estimated from the time of flight, and the profile map is computed from the average temperatures. The profile map may be computed using a polynomial approximation technique, a grid optimization technique or a basis function technique.

42 ENGINEERING↗
DOE OSTI2025

Wiener Filter Deconvolution for Analog Signals

This project explores the use of Wiener deconvolution to recover an original signal that has been distorted by a known transfer function and by noise. A simulated Gaussian pulse was used as the test signal, and a transfer function was applied in the frequency domain to model the system distortion. Controlled noise was then introduced to approximate real-world signal degradation. A Wiener filter was implemented to reverse the effects of the transfer function while minimizing the influence of noise. The recovered signal was compared with the original pulse to evaluate the effectiveness of the filter. The results demonstrate that Wiener deconvolution offers a stable and effective approach to signal recovery, balancing complete transfer function inversion and noise suppression.

Campbell, Savanna [Fermilab]↗
DOE OSTI2018

Monte Carlo calculation and verification of the geometrical factors for the NPDGamma experiment

The NPDGamma experiment measures the parity-violating asymmetry in gamma-ray emission in the capture of polarized neutrons on liquid parahydrogen. The sensitivity to the asymmetry for each detector in the array is used as a parameter in the extraction of the physics asymmetry from the measured data. The detector array is approximately cylindrically symmetric around the target and a step-wise sinusoidal function has been used for the sensitivity in the previous iteration of the NPDGamma experiment, but deviations from cylindrical symmetry necessitate the use of a Monte Carlo model to determine corrections to the geometrical factors. For the calculations, source code modifications to MCNPX were done in order to calculate the sensitivity of each cesium iodide detector to the physics asymmetry. We describe the MCNPX model and results from calculations and how the results are validated through measurement of the parity violating asymmetry of gamma-rays from neutron capture on chlorine.

Grammer, K. B.↗
DOE OSTI2017

Approximate analytic solutions to coupled nonlinear Dirac equations

Here, we consider the coupled nonlinear Dirac equations (NLDEs) in 1+11+1 dimensions with scalar–scalar self-interactions g 1 2 /2($\bar{ψ}$ψ) 2 + g 2 2 /2($\bar{Φ}$Φ) 2 + g 2 3 ($\bar{ψ}$ψ)($\bar{Φ}$Φ) as well as vector–vector interactions g 1 2 /2($\bar{ψ}$γμψ)($\bar{ψ}$γμψ) + g 2 2 /2($\bar{Φ}$γμΦ)($\bar{Φ}$γμΦ) + g 2 3 ($\bar{ψ}$γμψ)($\bar{Φ}$γμΦ). Writing the two components of the assumed rest frame solution of the coupled NLDE equations in the form ψ=e –iω1t R 1 cosθ,R 1 sinθΦ=e –iω2t R 2 cosη,R 2 sinη, and assuming that θ(x),η(x) have the same functional form they had when g3 = 0, which is an approximation consistent with the conservation laws, we then find approximate analytic solutions for Ri(x) which are valid for small values of g 3 2 /g 2 2 and g 3 2 /g 1 2 . In the nonrelativistic limit we show that both of these coupled models go over to the same coupled nonlinear Schrödinger equation for which we obtain two exact pulse solutions vanishing at x → ±∞.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

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

A comparison of QTP functionals against coupled-cluster methods for EAs of small organic molecules

EA-EOM-CCSD electron affinities and LUMO energies of various Kohn–Sham density functional theory (DFT) methods are calculated for an a priori IP benchmark set of 64 small, closed-shell molecules. The purpose of these calculations was to investigate whether the QTP KS-DFT functionals can emulate EA-EOM-CC with only a mean-field approximation. We show that the accuracy of DFT—relative to CCSD—improves significantly when elements of correlated orbital theory are introduced into the parameterization to define the QTP family of functionals. In particular, QTP(02), which has only a single range separation parameter, provides results accurate to a MAD of <0.15 eV for the whole set of 64 molecules compared to EA-EOM-CCSD, far exceeding the results from the non-QTP family of density functionals.

Chemistry↗
DOE OSTI2016

Optimization of an exchange-correlation density functional for water

We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the “correct” parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems.

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

Measuring neutrino mass imprinted on the anisotropic galaxy clustering

The anisotropic galaxy clustering of large scale structure observed by the Baryon Oscillation Spectroscopic Survey Data Release 11 is analyzed to probe the sum of neutrino masses in the small m {sub ν} ∼< 1 eV limit in which the early broadband shape determined before the last scattering surface is immune from the variation of m {sub ν}. The signature of m {sub ν} is imprinted on the altered shape of the power spectrum at later epoch, which provides an opportunity to access the non-trivial m {sub ν} through the measured anisotropic correlation function in redshift space (hereafter RSD instead of Redshift Space Distortion). The non-linear RSD corrections with massive neutrinos in the quasi linear regime are approximately estimated using one-loop order terms. We suggest an approach to probe m {sub ν} simultaneously with all other distance measures and coherent growth functions, exploiting this deformation of the early broadband shape of the spectrum at later epoch. If the origin of cosmic acceleration is unknown, m {sub ν} is poorly determined after marginalizing over all other observables. However, we find that the measured distances and coherent growth functions are minimally affected by the presence of mild neutrino mass. Although the standard model of cosmic acceleration is assumed to be the cosmological constant, the constraint on m {sub ν} is little improved. Interestingly, the measured Cosmic Microwave Background (hereafter CMB) distance to the last scattering surface sharply slices the degeneracy between the matter content and m {sub ν}, and the m {sub ν} is observed to be m {sub ν} = 0.19{sup +0.28}{sub −0.17} eV which is different from massless neutrino at 68% confidence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Properties of a potential energy matrix in oscillator basis

Matrix elements of potential energy are examined in detail. We consider a model problem — a particle in a central potential. The most popular forms of central potential are taken up, namely, square-well potential, Gaussian, Yukawa and exponential potentials. We study eigenvalues and eigenfunctions of the potential energy matrix constructed with oscillator functions. It is demonstrated that eigenvalues coincide with the potential energy in coordinate space at some specific discrete points. We establish approximate values for these points. It is also shown that the eigenfunctions of the potential energy matrix are the expansion coefficients of the spherical Bessel functions in a harmonic oscillator basis. We also demonstrate a close relation between the separable approximation and L{sup 2} basis (J-matrix) method for the quantum theory of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A 3D approximate maximum likelihood localization solver

A robust three-dimensional solver was needed to accurately and efficiently estimate the time sequence of locations of fish tagged with acoustic transmitters and vocalizing marine mammals to describe in sufficient detail the information needed to assess the function of dam-passage design alternatives and support Marine Renewable Energy. An approximate maximum likelihood solver was developed using measurements of time difference of arrival from all hydrophones in receiving arrays on which a transmission was detected. Field experiments demonstrated that the developed solver performed significantly better in tracking efficiency and accuracy than other solvers described in the literature.

Source record↗
DOE OSTI2025

The Langdon effect in laser plasmas: Absorption and conduction

A plasma heated by inverse bremsstrahlung absorption of laser light develops a non-Maxwellian electron distribution function, called the Langdon effect [A. B. Langdon, Phys. Rev. Lett. 44, 575 (1980)]. These non-Maxwellian distributions are sufficiently long-lived to impact the absorption processes itself as well as the transport of heat by electrons. The theory of the Langdon effect in a homogeneous plasma is reviewed to clarify some aspects of Langdon's derivation as well as to confirm that the widely used super-Gaussian approximation works fairly well to describe the shape of the distribution function and reduction of the absorption rate. The Langdon effect on thermal conduction in an inhomogeneous plasma is developed by considering perturbations in a homogeneous absorbing plasma, which develops a heat flux due to both temperature and density gradients. A practical theory of the heat flux is developed by fitting the results of Vlasov–Fokker–Planck simulations, which avoids several approximations that compromised the usefulness of past theoretical predictions, most critically, the effect of electron–electron collisions on the fluxes. The present fits parameterize the coefficients of the temperature gradient (thermal conductivity) and the density gradient for a plasma of any ionization state and for any laser intensity where the theory of the Langdon effect remains locally valid. It is expected that this generalized theory of heat flow in an absorbing plasma will improve the predictive capability of radiation-hydrodynamics simulations of laser-produced plasmas, especially those formed in inertial confinement fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations

We study Chebyshev filter diagonalization as a tool for the computation of many interior eigenvalues of very large sparse symmetric matrices. In this technique the subspace projection onto the target space of wanted eigenvectors is approximated with filter polynomials obtained from Chebyshev expansions of window functions. After the discussion of the conceptual foundations of Chebyshev filter diagonalization we analyze the impact of the choice of the damping kernel, search space size, and filter polynomial degree on the computational accuracy and effort, before we describe the necessary steps towards a parallel high-performance implementation. Because Chebyshev filter diagonalization avoids the need for matrix inversion it can deal with matrices and problem sizes that are presently not accessible with rational function methods based on direct or iterative linear solvers. To demonstrate the potential of Chebyshev filter diagonalization for large-scale problems of this kind we include as an example the computation of the 10{sup 2} innermost eigenpairs of a topological insulator matrix with dimension 10{sup 9} derived from quantum physics applications.

97 MATHEMATICS AND COMPUTING↗