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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Analysis of Airport Ground Delay Program Decisions Using Data Mining Techniques

Air traffic service providers have to make decisions regarding changes to air traffic flow in the event of major weather disturbances and traffic congestions to maintain safety of the system. The behavior of the air traffic management system will be more predictable if consistent decisions are made under similar traffic and weather conditions. Consistency of deciding on control action depends on the weather and traffic conditions as well as accuracy in predicting these conditions. Weather parameters (defined in terms of forecast and actual weather and traffic conditions) on different days can be used to categorize days into days with little decision consistency, days with moderate decision consistency and days with high decision consistency. Four years of traffic, weather and ground delay program decisions data at major airports in the United States are used in the analysis. This paper examines performance of different data mining methods in the three regions of decision consistency. Not surprisingly, data mining methods have the best performance in the region of most decision consistency and have the poorest performance in the region of little decision consistency. In applications where data mining methods have differing performance in differing regions, it would be more useful to characterize the region specific performance instead of characterizing performance by a single parameter. Finally, the results show no significant variation in the performance of different data mining methods for this particular problem. The fact that different mining methods show no significant variation also provides further confidence in the results of data mining methods. Work in this abstract discusses initial results. This paper describes the results in terms of both forecast and actual environmental conditions and discusses how prediction errors impact decision consistency.

Kulkarni, Deepak↗
NASA NTRS1990

Variational formulation of high-performance finite elements - Parametrized variational principles

High-performance (HP) elements are simple finite elements constructed to deliver engineering accuracy with coarse arbitrary grids. This paper is part of a series on the variational basis of HP elements, with emphasis on those constructed with the free formulation (FF) and assumed natural strain (ANS) methods. The present paper studies parametrized variational principles that provide a foundation for the FF and ANS methods, as well as for a combustion of both methods.

Felippa, C. A.↗
DOE OSTI2021

Variational Information Planning for Sequential Decision Making

We consider the setting of sequential decision making where, at each stage, potential actions are evaluated based on expected reduction in posterior uncertainty, given by mutual information (MI). As MI typically lacks a closed form, we propose an approach which maintains variational approximations of, both, the posterior and MI utility. Our planning objective extends an established variational bound on MI to the setting of sequential planning. The result, variational information planning (VIP), is an efficient method for sequential decision making. We further establish convexity of the variational planning objective and, under conditional exponential family approximations, we show that the optimal MI bound arises from a relaxation of the well-known exponential family moment matching property. Here, we demonstrate VIP for sensor selection, experiment design, and active learning, where it meets or exceeds methods requiring more computation, or those specialized to the task.

Pacheco, Jason↗
NASA NTRS1984

Abnormal Circulation Changes in the Winter Stratosphere, Detected Through Variations of D Region Ionospheric Absorption

A method to detect stratospheric warmings using ionospheric absorption records obtained by an Absorption Meter (method A3) is introduced. The activity of the stratospheric circulation and the D region ionospheric absorption as well as other atmospheric parameters during the winter anomaly experience an abnormal variation. A simultaneity was found in the beginning of abnormal variation in the mentioned parameters, using the absorption records for detecting the initiation of the stratospheric warming. Results of this scientific experience of forecasting in the El Arenosillo Range, are presented.

Delamorena, B. A.↗
DOE OSTI2018

Microsegregation Measurement: Methods and Applications

Microsegregation is chemical variation at the microstructural level in a solidified ingot, which can cause undesirable effects on material properties. Its formation is driven by different levels of solute elemental solubility between solid and liquid phases, resulting in solute buildup or depletion in the liquid phase during solidification. Characterization of microsegregation can be essential for understanding its formation and removing it through heat treatment. Measurement is frequently carried out via scanning electron microscope (SEM), though metallographic analysis, mechanical testing and optical radiography have been used. SEM analysis can consist of point measurements on single dendrites, line scans across or along dendrites, elemental maps, and extended line-based or grid-based analyses extending over many dendrites, which can generate large amounts of data. Line-based or grid-based data can be analyzed using scatterplots of element concentrations to visualize segregation behavior and sorting/ranking of data by a single element or a combination of elements, to give compositional profiles as a function of extent of solidification. For both qualitative and quantitative analyses, effects of phase transformations and secondary phases need to be considered as well as sorting of measurement errors when generating compositional profiles. Applications of analysis include determination of partition ratios, examination of solute–element interactions, determination of effects of cooling rate and back diffusion, applications to freckling studies, validations of modeling, and development of heat-treatment profiles.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Hamilton's principle - Finite-element methods and flexible body dynamics

A variational formulation is given for the equations of motion for an unconstrained elastic body, and Hamilton's principle is used to derive the equations of motion and deformation of the body. Finite-element approximations are developed for these dynamical equations with respect to a body axis system satisfying the mean axis conditions. The free-body influence matrix for the body then is developed in terms of the finite-element model parameters.

Cavin, R. K., III↗
NASA NTRS1989

Retrieval of total precipitable water using radiometric measurements near 92 and 183 GHz

This paper describes a method for the retrieval of total precipitable water (W) in dry atmospheres, which relies on the strong water vapor absorption line at 183 GHz (for W less than 0.6 g/sq cm) and on the absorption near 90 GHz (for W above 0.6 g/sq cm). The method is very sensitive to the variations of W less than 0.5 g/sq cm and is complementary to the established methods that use the weak 22 GHz water vapor absorption line to retrieve W in the 1-6 g/sq cm. The technique was demonstrated by the analysis of two Advanced Microwave Moisture Sounder observations of dry atmospheres following cold-air outbreaks on March 13, 1983, and February 23, 1986.

Wang, J. R.↗
DOE OSTI2025

Modeling Charge Transfer Excitation with Variation After Response Quantum Monte Carlo (Final Technical Report)

The electronically excited states of molecules play critical roles both in spectroscopy and in technological application, ranging from UV/Vis absorption experiments to the production of solar fuels. However, certain types of excitations, such as charge transfer excitations, remain difficult to treat with existing methods such as time dependent density functional theory or linear response coupled cluster. Under this Early Career Research Award, we developed a number of different theoretical methods to address this issue, which began with Variational Monte Carlo methods and branched into multiple types of excited-state-specific quantum chemistry.

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

Completed Beltrami-Michell formulation for analyzing mixed boundary value problems in elasticity

In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the Beltrami-Michell formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the more prevalent displacement of mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted application, could not become a true alternative to the Navier's displacement method, which can solve all three types of boundary value problems. The restrictions in the BMF have been alleviated by augmenting the classical formulation with a novel set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force formulation, has been termed the completed Beltrami-Michell formulation (CBMF). The CBMF can solve general elasticity problems with stress, displacement, and mixed boundary conditions in terms of stresses as the primary unknowns. The CBMF is derived from the stationary condition of the variational functional of the integrated force method. In the CBMF, stresses for kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. This paper presents the CBMF and its derivation from the variational functional of the integrated force method. Several examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value problems under thermomechanical loads. Selected example problems include a cylindrical shell wherein membrane and bending responses are coupled, and a composite circular plate.

Patnaik, Surya N.↗
NASA NTRS1980

Determination of variations of the solar radius from solar eclipse observations

This paper describes the method to determine the solar radius and its variations from observations made during total solar eclipses. In particular, the procedure to correct the spherical moon predictions for the effects of lunar mountains and valleys on the width and location of the path of totality is addressed in detail. The errors affecting this technique are addressed, a summary of the results of its application to three solar eclipses are presented, and the implications of the results on the constancy of the solar constant are described.

Sofia, S.↗
NASA NTRS1992

The variation of reconnection rate at the dayside magnetopause and cusp ion precipitation

A method is presented which allows estimation of the variation of the rate of magnetic reconnection at the dayside magnetopause. This is achieved using observations of the cusp particle precipitation made by low-altitude polar-orbiting spacecraft. In this paper we apply the technique to a previously published example of a cusp intersection by the DMSP F7 satellite. It is shown that the cusp signature in this case was produced by three separate bursts of reconnection which were of the order of 10 min apart, each lasting roughly 1 min. This is similar to the variation of reconnection rate which is required to explain typical flux transfer event signatures at the magnetopause.

Lockwood, M.↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2023

Dimensionality Reduction with Variational Encoders Based on Subsystem Purification

Efficient methods for encoding and compression are likely to pave the way toward the problem of efficient trainability on higher-dimensional Hilbert spaces, overcoming issues of barren plateaus. Here, we propose an alternative approach to variational autoencoders to reduce the dimensionality of states represented in higher dimensional Hilbert spaces. To this end, we build a variational algorithm-based autoencoder circuit that takes as input a dataset and optimizes the parameters of a Parameterized Quantum Circuit (PQC) ansatz to produce an output state that can be represented as a tensor product of two subsystems by minimizing $Tr(ρ^2)$. The output of this circuit is passed through a series of controlled swap gates and measurements to output a state with half the number of qubits while retaining the features of the starting state in the same spirit as any dimension-reduction technique used in classical algorithms. The output obtained is used for supervised learning to guarantee the working of the encoding procedure thus developed. We make use of the Bars and Stripes (BAS) dataset for an 8 × 8 grid to create efficient encoding states and report a classification accuracy of 95% on the same. Thus, the demonstrated example provides proof for the working of the method in reducing states represented in large Hilbert spaces while maintaining the features required for any further machine learning algorithm that follows.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Planetary atmosphere representation for Monte-Carlo simulation of aerospace vehicle trajectories

Methods of simulating large scale random atmospheric variations are examined for use in Monte Carlo simulations of hypervelocity vehicles. Random density, temperature, and pressure variations have a significant impact on the performance of such vehicles. Methods for generating random atmosphere profiles are described and compared in the context of a Monte Carlo simulation of the Aeroassist Flight Experiment (AFE). It is seen that these techniques result in statistically significant differences in key simulation point statistics.

Dutton, K.↗
DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

54 ENVIRONMENTAL SCIENCES↗