Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Additional Trials & Tribulations in Synthesizing a Sulfide Standard

Well-characterized sulfide reference materials that can serve as matrix-matched calibrants for in-situ trace element analyses via laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) remain elusive. Here, we describe the creation an in-house sulfide standard at NASA JSC that will be used to measure siderophile and chalcophile trace elements in low pressure experimental products, specifically pentlandite ([FeNi]9S8) and pyrrhotite ([FeNi]1-xS), as well as natural sulfides in terrestrial and meteoritic mineral assemblages. Recent methods in creating a homogenous sulfide standard include pressed pellets, synthesized chips, or fused glass [e.g. 1-3]. Other studies, however, have had success in forming homogenous FeS standards via synthesization methods [e.g. 4-6]. We present a variation of the method described in [4], in which we create a pyrrhotite standard doped with a variety of trace elements (Zn, As, Se, Mo, Ru, Rh, Pd, Sn, Sb, Te, W, Os, Ir, Pt, Au) at ~10-40 ppm and Cu at ~200 ppm. Sulfide compositions were constructed using Fe and Ni metal and elemental sulfur powders. Trace elements were added to elemental sulfur from Atomic Absorption Spectroscopy (AAS) elemental standards as either nitrate or chloride solutions to prevent oxidation of the metal powders. The dried sulfur mixture was mixed with Fe and Ni powders and mechanically mixed before two aliquots were placed in separate SiO2 tubes. Each tube was held under vacuum for ~30 minutes, sealed under vacuum, and then heated at 800C for 48 hours. Like [4], the synthesis products were composed of pourous sulfide crystals. Major element analyses of both experimental aliquots, collected using an Electron Probe Microanalyser (EPMA), yield indistinguishable major element compositions (uncertainties in 2SE & 2RSE[%]), with an average of 57.90  0.09 (0.16 %), 4.95  0.03 (0.55 %), and 39.70  0.14 (0.34 %), for Fe, Ni, and S, respectively. Trace element data were measured using a Photon Machines 193nm laser ablation system coupled to a Thermo-Scientific Element-XR ICP-MS. Spot sizes were limited to 50 μm due to the porous nature of the sulfide target material. Trace element abundances, normalized to Fe as an internal standard, were also found to be homogenous between the two aliquots, with weighted mean 2RSE (%) values of <3.0 for all trace elements. Synthesized products were re-powdered and absolute concentrations measured via solution ICP-MS. Although the sulfide appears to be homogenous, sintering experiments will be performed to more closely match the standard density to natural sulfides and minimize differences in ablation behavior. Improved density also allows for higher sensitivity (i.e. more compact target material) and larger spot sizes or traverses, as void space is eliminated. Additional major and trace element analyses on the products of the sintering experiments will be undertaken.

Jacob B Setera↗
DOE OSTI2019

Managing variations among nodes in parallel system frameworks

Systems, apparatuses, and methods for managing variations among nodes in parallel system frameworks. Sensor and performance data associated with the nodes of a multi-node cluster may be monitored to detect variations among the nodes. A variability metric may be calculated for each node of the cluster based on the sensor and performance data associated with the node. The variability metrics may then be used by a mapper to efficiently map tasks of a parallel application to the nodes of the cluster. In one embodiment, the mapper may assign the critical tasks of the parallel application to the nodes with the lowest variability metrics. In another embodiment, the hardware of the nodes may be reconfigured so as to reduce the node-to-node variability.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

D-Wave Electron-H, -He+, and -Li2+ Elastic Scattering and Photoabsorption in P States of Two-Electron Systems

In previous papers [A. K. Bhatia, Phys. Rev. A 85, 052708 (2012); 86, 032709 (2012); 87, 042705 (2013)] electron-H, -He+, and -Li2+ P-wave scattering phase shifts were calculated using the variational polarized orbital theory. This method is now extended to the singlet and triplet D-wave scattering in the elastic region. The long-range correlations are included in the Schrodinger equation by using the method of polarized orbitals variationally. Phase shifts are compared to those obtained by other methods. The present calculation provides results which are rigorous lower bonds to the exact phase shifts. Using the presently calculated D-wave and previously calculated S-wave continuum functions, photoionization of singlet and triplet P states of He and Li+ are also calculated, along with the radiative recombination rate coefficients at various electron temperatures.

D-wave↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2018

Nodal variational principle for excited states

It is proven that the exact excited-state wave function and energy may be obtained by minimizing the energy expectation value of trial wave functions that are constrained only to have the correct nodes of the state of interest. This excited-state nodal minimum principle has the advantage that it requires neither minimization with the constraint of wave-function orthogonality to all lower eigenstates nor the antisymmetry of the trial wave functions. It is also found that the minimization over the entire space can be partitioned into several interconnected minimizations within the individual nodal regions, and the exact excited-state energy may be obtained by a minimization in just one or several of these nodal regions. For the proofs of the theorem, it is observed that the many-electron eigenfunction (excited state as well as ground state), restricted to a nodal region, is equivalent to a ground-state wave function of one electron in a higher-dimensional space; and, alternatively, an explicit excited-state energy variational expression is utilized by generalizing the Jacobi method of multiplicative variation. In corollaries, error functions are constructed for cases for which the nodes are not necessarily exact. The exact nodes minimize the energy error functions with respect to nodal variations. In conclusion, simple numerical illustrations of the error functions are presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Marginal unbiased score expansion and application to CMB lensing

Here, we present the marginal unbiased score expansion (MUSE) method, an algorithm for generic high-dimensional hierarchical Bayesian inference. MUSE performs approximate marginalization over arbitrary non-Gaussian latent parameter spaces, yielding Gaussianized asymptotically unbiased and near-optimal constraints on global parameters of interest. It is computationally much cheaper than exact alternatives like Hamiltonian Monte Carlo (HMC), excelling on funnel problems which challenge HMC, and does not require any problem-specific user supervision like other approximate methods such as variational inference or many simulation-based inference methods. MUSE makes possible the first joint Bayesian estimation of the delensed Cosmic Microwave Background (CMB) power spectrum and gravitational lensing potential power spectrum, demonstrated here on a simulated data set as large as the upcoming South Pole Telescope 3G 1500 deg 2 survey, corresponding to a latent dimensionality of ~6 million and of order 100 global bandpower parameters. On a subset of the problem where an exact but more expensive HMC solution is feasible, we verify that MUSE yields nearly optimal results. We also demonstrate that existing spectrum-based forecasting tools which ignore pixel-masking underestimate predicted error bars by only ~10%. This method is a promising path forward for fast lensing and delensing analyses which will be necessary for future CMB experiments such as SPT-3G, Simons Observatory, or CMB-S4, and can complement or supersede existing HMC approaches. The success of MUSE on this challenging problem strengthens its case as a generic procedure for a broad class of high-dimensional inference problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

Seasonal variation of the vertical distribution of stratospheric ozone as observed with the Umkehr and BUV methods

Month-to-month variations evidenced by ozone profiles inferred from the classical Umkehr observations and from the back-scattered ultraviolet (BUV) satellite observations made from the Nimbus 4 satellite are examined. Upper stratospheric ozone profiles derived from BUV and Umkehr data display similar seasonal variations of about the same phase and magnitude for the 38 km to 50 km region. Between 28 km and 38 km, the seasonal variations are less marked, but the same rough picture emerges for both data sets. If both data sets indicate an increasing (or decreasing) trend over a period of years, it is not possible to conclude that a trend exists unless separate means exists for monitoring stratospheric dust. Because of the stratosphere well above the Junge layer seems less likely to be affected by volcanic debris, BUV data should be superior to Umkehr data for monitoring trends in the 38 km to 50 km range, provided that the calibration problems of flying such a monitoring instrument in space can be overcome.

Mateer, C. L.↗
DOE OSTI2023

Optimizing training trajectories in variational autoencoders via latent Bayesian optimization approach *

Unsupervised and semi-supervised ML methods such as variational autoencoders (VAE) have become widely adopted across multiple areas of physics, chemistry, and materials sciences due to their capability in disentangling representations and ability to find latent manifolds for classification and/or regression of complex experimental data. Like other ML problems, VAEs require hyperparameter tuning, e.g. balancing the Kullback–Leibler and reconstruction terms. However, the training process and resulting manifold topology and connectivity depend not only on hyperparameters, but also their evolution during training. Because of the inefficiency of exhaustive search in a high-dimensional hyperparameter space for the expensive-to-train models, here we have explored a latent Bayesian optimization (zBO) approach for the hyperparameter trajectory optimization for the unsupervised and semi-supervised ML and demonstrated for joint-VAE with rotational invariances. We have demonstrated an application of this method for finding joint discrete and continuous rotationally invariant representations for modified national institute of standards and technology database (MNIST) and experimental data of a plasmonic nanoparticles material system. The performance of the proposed approach has been discussed extensively, where it allows for any high dimensional hyperparameter trajectory optimization of other ML models.

42 ENGINEERING↗
NASA NTRS1978

Drive train dynamic analysis

A method for parametric variations in drive train dynamic analysis is described. The method models the individual components of a drive system, forms the appropriate system interface coordinates and, calculates the system dynamic response at particular frequencies. Application of the method for prediction of the dynamic response characteristics of a helicopter transmission, and a comparison of results with test data are also included.

Giansante, N.↗
NASA NTRS1986

On the Maxwellian distribution, symmetric form, and entropy conservation for the Euler equations

The Euler equations of gas dynamics have some very interesting properties in that the flux vector is a homogeneous function of the unknowns and the equations can be cast in symmetric hyperbolic form and satisfy the entropy conservation. The Euler equations are the moments of the Boltzmann equation of the kinetic theory of gases when the velocity distribution function is a Maxwellian. The present paper shows the relationship between the symmetrizability and the Maxwellian velocity distribution. The entropy conservation is in terms of the H-function, which is a slight modification of the H-function first introduced by Boltzmann in his famous H-theorem. In view of the H-theorem, it is suggested that the development of total H-diminishing (THD) numerical methods may be more profitable than the usual total variation diminishing (TVD) methods for obtaining wiggle-free solutions.

Deshpande, S. M.↗
NASA NTRS2019

Elemental Variation and Petrogenesis of Pyroxenes in HED Meteorites by Multivariable Discriminant Analysis Method

HED (Howardite, Eucrite and Diogenite) are meteorites with mafic and ultramafic igneous composition. Previous studies suggested HED came from asteroid (4) Vesta and they were generated by magmatic melting followed by differentiation crystallization, metamorphic, and impact. Uniform oxygen isotopic composition of HED samples favors global magma ocean for (4) Vesta. However, the petrological diversity of HED may indicate more complex magma processes in the meteorites. In general, most of the geochemical studies are using traditional methods (e.g. element to element or ratio to ratio plots) to classify the different rock type and groups. With the traditional methods, only limited elements can be shown in a figure. As such, some meteorites may have been identified as HED by a traditional method but found different either in isotopic composition or other elemental characteristic. These anomalous HED meteorites may or may not come from asteroid (4) Vesta. In fact, magma is a unit system where any elemental change should affect to all other elements as a whole, it would be reasonable to consider all elements together to look for the systematical changes. Using multivariable discriminant analysis (MDA) method is one of such testing for their geochemical variation. The method may be able to help us to better understand the petrologic processes and the relationship among elements. This study is to test the MDA method by focus mainly on pyroxene composition of eucrite and diogenite.

Peng, Z. X.↗
NASA NTRS1985

Color measurement and discrimination

Theories of color measurement attempt to provide a quantative means for predicting whether two lights will be discriminable to an average observer. All color measurement theories can be characterized as follows: suppose lights a and b evoke responses from three color channels characterized as vectors, v(a) and v(b); the vector difference v(a) - v(b) corresponds to a set of channel responses that would be generated by some real light, call it *. According to theory a and b will be discriminable when * is detectable. A detailed development and test of the classic color measurement approach are reported. In the absence of a luminance component in the test stimuli, a and b, the theory holds well. In the presence of a luminance component, the theory is clearly false. When a luminance component is present discrimination judgements depend largely on whether the lights being discriminated fall in separate, categorical regions of color space. The results suggest that sensory estimation of surface color uses different methods, and the choice of method depends upon properties of the image. When there is significant luminance variation a categorical method is used, while in the absence of significant luminance variation judgments are continuous and consistant with the measurement approach.

Wandell, B. A.↗
DOE OSTI2021

Scalable Programming Workflows for Validation of Quantum Computers

Hybrid quantum-classical workflows have become standard methods for executing variational algorithms and other quantum simulation techniques, which are key applications for noisy intermediate scale quantum (NISQ) computers. Validating these simulations is an important task which helps gauge the progress of quantum computer development, and classical simulation can serve as a tool to this end. Both exact and more scalable approximate methods with quantifiable error bounds can be used in validation tasks where the applicable metrics include the distance from a calculable ground truth, the quality of an error model fit to data, etc. Here we present a library extension that includes methods for validation of quantum simulations based on scalable hybrid workflows executable on high performance computers. We provide examples that use approximate methods based on tensor networks and stabilizer simulators to bound the error of quantum simulations on NISQ hardware.

Nguyen, Thien↗
DOE OSTI2019

The Dominant Influencer of Voltage Fluctuation (DIVF) for Power Distribution System

Variability in power generation at multiple consumer owned renewable energy sources can cause random voltage fluctuations, which may result in voltage violations at multiple nodes in a power distribution system. Traditional methods of voltage control are ineffective to address fast voltage variations. New dynamic methods of voltage control using end user devices can significantly benefit from the knowledge of the dominant influencer of voltage fluctuation (DIVF). This paper introduces information theoretic metrics to identify DIVF for every node in the power distribution system. A node can be DIVF due to its location in the distribution system and its capacity to generate and change its power. Use of traditional power flow algorithms to identify the DIVF is computationally complex which limits their use in real-time applications. In this work, information theoretic indicators of the DIVF are computed analytically in a computationally efficient manner, and their effectiveness to identify the DIVF is tested using the IEEE 69 node test system. From simulations, it is illustrated that differential entropy, Kullback-Leibler distance and Frechet distance are effective indicators of the DIVF and can successfully identify the most influential nodes providing the foundation for intelligent voltage control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Comparison of One- and Two-Variable Linear Regression Models and Classic Energy Intensity for Energy Performance Tracking of Two Manufacturing Sectors

We report manufacturing facilities consumed about 32% of total domestic energy in the United States in 2016. To evaluate the energy savings achieved through the implementation of energy conservation projects and to establish powerful arguments for future projects, accurate energy performance tracking methods are necessary. The classic energy intensity method (i.e., the ratio of annual total energy over annual total production) is the means most widely used to measure savings because it can be understood and calculated easily. This method considers the variation in production rates to some extent; however, it fundamentally assumes facilities' base energy consumption (energy consumption at zero production) to be zero, which rarely holds true. Furthermore, this method does not consider variations in other relevant parameters, such as weather conditions. Therefore, the regression models approach is commonly recommended to track energy performance improvements. However, because it requires more data and statistical expertise, the regression models approach has been adopted by only a few facilities, and many others suspect it is not worth the effort for their specific cases. For this reason, the improved accuracy this approach offers needs to be demonstrated. By analyzing 477 monthly energy (electricity and natural gas) data sets, this study quantitatively compared the accuracy of classic energy intensity, one-variable (production only) linear regression models, and two-variable (i.e., production and weather) linear regression models for two manufacturing sectors (primary metal and transportation equipment manufacturing). Results showed that significant improvements in accuracy were achieved with one-variable regression models when compared with the classic energy intensity method and with two-variable regression models when compared with one-variable regression models. These improvements were demonstrated using the p values of intercept, cooling degree days, and heating degree days. Based on these results, to achieve a good balance between accuracy improvements and resources requirements, one-variable (production only) linear regression models for electricity consumption and two-variable (production and heating degree days) linear regression models for natural gas consumption are recommended for both sectors in the case of limited resources. Lastly, facilities can also use the results to decide which approaches fit better in their cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2007

Projection Method for Flows with Large Density Variations

Numerical models of solidification including a mushy zone are notoriously inefficient; most of them are based on formulations that require the coupled solution of the velocity components in the momentum equation greatly restricting the range of applicability of the models. There are only two models known to the authors that have used a projection or fractional step formulation, but none of these were used to model problems of any significant size. A third model was only applied to a partial mushy zone with no all-fluid region. Our initial attempts at modeling directional solidification in the presence of a developing mushy zone using a projection formulation encountered very serious difficulties once solidification starts. These difficulties were traced to the inability of the method to deal with large local density differences in the vicinity of the fluid-mush interface. As a result, a modified formulation of the projection method has been developed, that maintains the coupling between the body force and the pressure gradient and is presented in this work. The new formulation is shown to be robust and efficient, and can be applied to problems involving very large meshes. This is illustrated in this work through its application to simulations involving Pb-Sb and Pb-Sn alloys.

Heinrich, Juan C.↗
DOE OSTI2024

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗