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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Variance of phase fluctuations of waves propagating through a random medium

As an electromagnetic wave propagates through a random scattering medium, such as a forest, its energy is attenuated and random phase fluctuations are induced. The magnitude of the random phase fluctuations induced is important in estimating how well a Synthetic Aperture Radar (SAR) can image objects within the scattering medium. The two-layer random medium model, consisting of a scattering layer between free space and ground, is used to calculate the variance of the phase fluctuations induced between a transmitter located above the random medium and a receiver located below the random medium. The scattering properties of the random medium are characterized by a correlation function of the random permittivity fluctuations. The effective permittivity of the random medium is first calculated using the strong fluctuation theory, which accounts for large permittivity fluctuations of the scatterers. The distorted Born approximation is used to calculate the first-order scattered field. A perturbation series for the phase of the received field in the Rytov approximation is then introduced and the variance of the phase fluctuations is also calculated assuming that the transmitter and receiver are in the paraxial limit of the random medium, which allows an analytic solution to be obtained. Results are compared using the paraxial approximation, scalar Green's function formulation, and dyadic Green's function formulation. The effects studied are the dependence of the variance of the phase fluctuations on receiver location in lossy and lossless regions, medium thickness, correlation length and fractional volume of scatterers, depolarization of the incident wave, ground layer permittivity, angle of incidence, and polarization.

Chu, Nelson C.↗
NASA NTRS1980

Microstructure of the interplanetary magnetic field near 4 and 5 AU

Interplanetary magnetic field fluctuations near 4 and 5 AU are investigated on the basis of Pioneer 10 and 11 vector magnetic field measurements in order to identify extant wave modes and propagation directions as functions of frequency in the range 0.00004 to 0.009 Hz. Analysis of the approximately 500 spectral matrices as functions of frequency obtained from each spacecraft indicates the distributions of maximum variance directions to be peaked along the normal to the plane of the minimum variance direction and the mean magnetic field, as would be expected for a planar, small-amplitude Alfven wave with a minimum variance direction close to the wave vector. It is thus inferred that single small-amplitude Alfven waves are occasionally present in the temporal and frequency bands analyzed, however waves with wave vectors not well approximated by the minimum variance direction, Alfven wave ensembles, finite-amplitude waves and nonplanar fluctuations may also be present. Observations are also noted to suggest the predominance of the Alfven mode over the fast and slow modes, however do not support theoretical models of outward propagating plane Alfven waves of solar origin.

Parker, G. D.↗
DOE OSTI2023

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

First-principles studies of the atomic, elastic and vibrational properties of BaNbO{sub 3}

We have employed the ab initio pseudopotential method, within a generalized gradient approximation of the density functional theory, a linear response approach, to determine atomic, elastic and lattice dynamical properties of the BaNbO{sub 3}. The obtained structural and elastic parameter are in good agreement with previous experimental and theoretical results. The analysis of B/G ratio indicates that BaNbO{sub 3} is a ductile material. The calculated lattice dynamical properties compared with phonon spectrum of the similar structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

First-Principles Integrated Adsorption Modeling for Selective Capture of Uranium from Seawater by Polyamidoxime Sorbent Materials

Nuclear power is a relatively carbon-free energy source that has the capacity to be utilized today in an effort to stem the tides of global warming. The growing demand for nuclear energy, however, could put significant strain on our uranium ore resources, and the mining activities utilized to extract that ore can leave behind long-term environmental damage. A potential solution to enhance the supply of uranium fuel is to recover uranium from seawater using amidoximated adsorbent fibers. This technology has been studied for decades but is currently plagued by the material’s relatively poor selectivity of uranium over its main competitor vanadium. In this work, we investigate the binding schemes between uranium, vanadium, and the amidoxime functional groups on the adsorbent surface. Using quantum chemical methods, binding strengths are approximated for a set of complexation reactions between uranium and vanadium with amidoxime functionalities. Those approximations are then coupled with a comprehensive aqueous adsorption model developed in this work to simulate the adsorption of uranium and vanadium under laboratory conditions. Experimental adsorption studies with uranium and vanadium over a wide pH range are performed, and the data collected are compared against simulation results to validate the model. It was found that coupling ab initio calculations with process level adsorption modeling provides accurate predictions of the adsorption capacity and selectivity of the sorbent materials. Furthermore, this work demonstrates that this multiscale modeling paradigm could be utilized to aid in the selection of superior ligands or ligand compositions for the selective capture of metal ions. Furthermore, this first-principles integrated modeling approach opens the door to the in silico design of next-generation adsorbents with potentially superior efficiency and selectivity for uranium over vanadium in seawater.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Push it to the limit: comparing periodic and local approaches to density functional theory for intermolecular interactions

With the aim of systematically comparing two popular approaches to density functional theory–all-electron calculations with local basis sets, and periodic calculations employing plane wave basis sets and norm-conserving pseudopotentials–we have computed complete-basis binding energies across the S22 set of intermolecular interactions, a dataset consisting of noncovalent interactions of small- and medium-sized molecules containing first- and second-row atoms, using the Troullier-Martins norm-conserving pseudopotentials with SPW92, a local spin-density approximation; and PBE, a generalised gradient approximation. Here, we have found that it is challenging to reach the basis set limit with these periodic calculations; for the methods and systems examined, a minimum vacuum distance of 30 Å between a system and its nearest images is necessary–unless some form of dipole correction is employed–as is a kinetic energy cutoff of at least 80 Ry. The trends in convergence with respect to vacuum size and kinetic energy cutoff are largely independent of the level of density functional approximation employed. A sense of the impact of each hyperparameter on basis set error provides a foundation for ensuring quality calculations in future studies and allows us to quantify the basis set errors incurred in existing studies on similar systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Dynamics of water confined in lyotropic liquid crystals: Molecular dynamics simulations of the dynamic structure factor

The properties of water under confinement are of practical and fundamental interest. Here in this work we study the properties of water in the self-assembled lyotropic phases of gemini surfactants with a focus on testing the standard analysis of quasi-elastic neutron scattering (QENS) experiments. In QENS experiments the dynamic structure factor is measured and fit to models to extract the translational diffusion constant, D T , and rotational relaxation time, τ R . We test this procedure by using simulation results for the dynamic structure factor, extracting the dynamic parameters from the fit as is typically done in experiments, and comparing the values to those directly measured in the simulations. We find that the decoupling approximation, where the intermediate scattering function is assumed to be a product of translational and rotational contributions, is quite accurate. The jump-diffusion and isotropic rotation models, however, are not accurate when the degree of confinement is high. In particular, the exponential approximations for the intermediate scattering function fail for highly confined water and the values of D T and τ R can differ from the measured value by as much as a factor of two. Other models have more fit parameters, however, and with the range of energies and wave-vectors accessible to QENS, the typical analysis appears to be the best choice. In the most confined lamellar phase, the dynamics are sufficiently slow that QENS does not access a large enough time scale and neutron spin echo measurements would be a valuable technique in addition to QENS.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Experimentally-Derived Bidirectional Reflectance Distribution Function Data in Support of the Orbital Debris Program Office

The NASA Orbital Debris Program Office (ODPO) has used various optical assets to acquire photometric data of Earth-orbiting objects to define the orbital debris environment. To better characterize and model optical data acquired from ground-based telescopes, the Optical Measurements Center (OMC) at NASA Johnson Space Center emulates illumination conditions seen in space by using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. One of the OMC goals is to improve the size calculation used for optical data by developing an optical-based Size Estimation Model. The current size estimation requires applying a Lambertian phase function, a set albedo value, and range to the observed magnitude. The first step to improving the sampled brightness of laboratory targets is to remove aspect-angle dependencies. Then, the volume of possible object viewing angles is sampled at 21 combinations of azimuth and zenith angles for each solar phase angle. Finally, the acquired images are input into an image processing program that generates approximations for the object’s Bidirectional Reflectance Distribution Function (BRDF) and phase function. The BRDF is a radiometric concept that identifies an object’s material composition by matching a BRDF approximated with photometric data collected by ground-based telescopes with a BRDF generated experimentally from a known object in the laboratory. This paper discusses the validation of experimental BRDF and phase function approximations produced in the OMC and how the findings will be incorporated into ODPO models. A Lambertian sphere is imaged and the subsequent experimental functions are scrutinized to confirm that they correspond to an object that has an isotropic luminance. With the image processing algorithm validated, test objects with varying optical properties are then imaged to confirm that the produced photometric functions are both unique and repeatable. Once the validation is complete, the OMC will be used to evaluate a subset of fragments from a hypervelocity impact test of a mock-up satellite and assess the appropriate phase function and size estimates using BRDF measurements for a large volume of targets composed of various shapes, sizes, and materials.

Orbital Debris Program Office (ODPO)↗
DOE OSTI2018

Ab Initio Investigation of Structural Stability and Electronic and Magnetic Properties of the Half-Heusler Alloys: MTiSb (M = Fe, Co, and Ni)

The structural, electronic, and magnetic properties of the half-Heusler MTiSb (M = Fe, Co, Ni) compounds have been examined using the first-principles calculations with the full-potential linear augmented plane wave (FP-LAPW) method within the density functional theory. The electronic structures and magnetic properties of MTiSb (M = Fe, Co, Ni) are studied. It is found that the calculated lattice parameter is in good agreement with the theoretical values. Using the general gradient approximation (GGA), we observe that for the three compounds, the α phase is more stable than β and γ phases. CoTiSb is a semiconductor, NiTiSb is a metal, while the FeTiSb compound is a half-metallic, and the total magnetic moments are 0.00085, − 0.26074, and 0.99353 μ B, respectively, which is in agreement with the Slater–Pauling rule, Furthermore, the origin for the appearance of the half-metallic band gap in the FeTiSb compound was also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗
NASA NTRS1924

Static stability of seaplane floats and hulls

Values of lateral and longitudinal metacentric heights for various seaplanes were calculated by means of approximate formulae derived here. The data are given in tabular form. Upon plotting these metacentric heights against the corresponding gross weights, it appears that the metacentric height is approximately a straight line function of the gross weight. For the lateral metacentric height GM = 13 + .002 W and for longitudinal metacentric height GM = 15 + .002 W, GM is in feet and the gross weight (W) is in pounds. Although only approximate, it is thought that the values indicated here are a reliable guide to current practice. It is recommended that the longitudinal and lateral metacentric heights be made equal and of the value given by GM = 15 = .002 W. The proper length or spacing required to satisfy the indicated value may then be obtained from substitution in the approximate formulae for metacentric height.

Diehl, W S↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Exchange-correlation energies of atoms from efficient density functionals: influence of the electron density

The primary goal of Kohn–Sham density functional theory is to evaluate the exchange-correlation contribution to electronic properties. However, the accuracy of a density functional can be affected by the electron density. Here we apply the nonempirical Tao–Mo (TM) semilocal functional to study the influence of the electron density on the exchange and correlation energies of atoms and ions, and compare the results with the commonly used nonempirical semilocal functionals local spin-density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE), Tao–Perdew–Staroverov–Scuseria (TPSS), and hybrid functional PBE0. We find that the spin-restricted Hartree–Fock density yields the exchange and correlation energies in good agreement with the Optimized Effective Potential method, particularly for spherical atoms and ions. However, the errors of these semilocal and hybrid functionals become larger for self-consistent densities. We further find that the quality of the electron density have greater effect on the exchange-correlation energies of kinetic energy density-dependent meta-GGA functionals TPSS and TM than on those of the LSDA and GGA, and therefore, should have greater influence on the performance of meta-GGA functionals. Lastly, we show that the influence of the density quality on PBE0 is slightly reduced, compared to that of PBE, due to the exact mixing.

74 ATOMIC AND MOLECULAR PHYSICS↗