Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational methods”

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 865 records · Page 48

Estimation of surface moisture availability from remote temperature measurements

Analysis of moisture availability, determined from GOES IR surface temperature measurements and a boundary layer model, and antecedent precipitation are presented for data taken over Texas and Argentina. Correlations between antecedent precipitation and moisture availability are high where the spatial variation in moisture availability is much larger than an inherent uncertainty in the method, which is judged to be about + or - 0.2 for the moisture availability parameter. In regions where there are large gradients of precipitation, the IR method can provide useful information on variations in soil moisture availability and, by implication, in precipitation.

Flores, Alberto L.↗

On instruction sets and their formats

Central to instruction set design is the issue of the instruction format. Some common format encoding techniques are discussed, and a method of representing high-level language parse trees by means of formats that cover successively larger portions of a tree is introduced. Variations are introduced on the method that represent directed acyclic graphs as well as simple parse trees, and that encode constants in a special fashion. For a particular representation, the number of times each format is executed to run a sample program to completion is measured.

Flynn, M. J.↗

Support Vector Machines for Classification of Direct Energy Deposition Standoff Distance for Improved Process Control

A critical factor in the implementation of direct energy deposition is the ability to maintain the standoff distance between the nozzle and the build surface, as this influences powder capture efficiency and overall part quality. Due to process-related variations, layer height may vary, causing unintended variation in standoff distance and poor build quality. While prior work has utilized contact probing to qualify standoff distance during processing, in situ methods for qualification of standoff distance are of major interest. The present work seeks to understand efficacy of image-based methods for classifying standoff distance variation in real-time using support vector machines (SVMs). It was hypothesized that the size of the melt pool and the amount of spatter will have significant correlations with deviations in the standoff distance; thus, SVMs were used on a dataset that is comprised of morphological features of melt pool size and image entropy. The SVM model was used to classify melt pool images into categories according to standoff distance variation from nominal. K-folds cross validation was used to find the optimal hyperparameters for the SVM model. To understand the impact of the selected features on the classification performance and inference speed, multiple models were trained with differing numbers of included features. Results for classification score, inference time, and image preprocessing/feature extraction from these data are reported. The present results show that the SVM model was able to predict the standoff distance classification with an accuracy of 97 percent and a speed of 0.122 s per image, making it a viable solution for real-time control of standoff distance.

Klesmith, Zoe↗

Hybrid analytical technique for the nonlinear analysis of curved beams

The application of a two-step hybrid technique to the geometrically nonlinear analysis of curved beams is used to demonstrate the potential of hybrid analytical techniques in nonlinear structural mechanics. The hybrid technique is based on successive use of the perturbation method and a classical direct variational procedure. The functions associated with the various-order terms in the perturbation expansion of the fundamental unknowns, and their sensitivity derivatives with respect to material and geometric parameters of the beam, are first obtained by using the perturbation method. These functions are selected as coordinate functions (or modes) and the classical direct variational technique is then used to compute their amplitudes. The potential of the proposed hybrid technique for nonlinear analysis of structures is discussed. The effectiveness of the hybrid technique is demonstrated by means of numerical examples. The symbolic computation system Mathematica is used in the present study. The tasks performed on Mathematica include: (1) generation of algebraic expressions for the perturbation functions of the different response quantities and their sensitivity derivatives: and (2) determination of the radius of convergence of the perturbation series.

Noor, A. K.↗

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↗

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↗

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↗

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.↗

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↗

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.↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

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↗