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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

The method of models in variational scattering calculations.

It is shown that difficulties in atomic scattering calculations that stem from the use of inexact target wave functions can be overcome by employing the so-called method of models, provided that the projectile is distinguishable from the atomic electrons. The proposed method of models consists in replacing the target hamiltonian by a model hamiltonian, of which the appropriate target state is an eigenfunction. A connection with the positron scattering work of Peterkop and Rabik (1971) has been established.

Drachman, R. J.↗
NASA NTRS1976

A method to account for variation of average compressor inlet pressure during instantaneous distortion analyses

Instantaneous distortion analyses compare a time-varying value of an index (or 'surge margin used up') with a critical level (or 'available surge margin' of the compressor) to determine inlet-engine compatibility. Unless freestream conditions or propulsion system controls are changing, it is generally assumed that the available surge margin of the compressor is accurately determined from the steady-state operating point. Results are presented which show that variations of average compressor inlet pressure may occur without changes in freestream conditions or propulsion system controls. The volume dynamics of the compressor will cause these pressure variations to be attenuated and delayed by the time they reach the exit. This will cause the compressor pressure ratio (and available surge margin) to vary with time. A method is presented to calculate the available surge margin as a function of time and incorporate it into an instantaneous distortion analysis. Results show that inlet pressure variations which cause only a small change at the compressor exit can cause a significant variation in the available surge margin.

Burstadt, P. L.↗
NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗
NASA NTRS1998

Empirical correction for earth sensor horizon radiance variation

A major limitation on the use of infrared horizon sensors for attitude determination is the variability of the height of the infrared Earth horizon. This variation includes a climatological component and a stochastic component of approximately equal importance. The climatological component shows regular variation with season and latitude. Models based on historical measurements have been used to compensate for these systematic changes. The stochastic component is analogous to tropospheric weather. It can cause extreme, localized changes that for a period of days, overwhelm the climatological variation. An algorithm has been developed to compensate partially for the climatological variation of horizon height and at least to mitigate the stochastic variation. This method uses attitude and horizon sensor data from spacecraft to update a horizon height history as a function of latitude. For spacecraft that depend on horizon sensors for their attitudes (such as the Total Ozone Mapping Spectrometer-Earth Probe-TOMS-EP) a batch least squares attitude determination system is used. It is assumed that minimizing the average sensor residual throughout a full orbit of data results in attitudes that are nearly independent of local horizon height variations. The method depends on the additional assumption that the mean horizon height over all latitudes is approximately independent of season. Using these assumptions, the method yields the latitude dependent portion of local horizon height variations. This paper describes the algorithm used to generate an empirical horizon height. Ideally, an international horizon height database could be established that would rapidly merge data from various spacecraft to provide timely corrections that could be used by all.

Hashmall, Joseph A.↗
NASA NTRS1978

The development of methods for the prediction of primary creep behavior in metals

The applicability of a thermodynamic constitutive theory of deformation to the prediction of primary creep and creep strain relaxation behavior in metals is examined. Constitutive equations derived from the theory are subjected to a parametric analysis in order to determine the influence of several parameters on the curve forms generated by the equations. A computer program is developed which enables the solution of a generalized constitutive equation using experimental data as input. Several metals were tested to form a data base of primary creep and relaxation behavior. The extent to which these materials conformed to the constitutive equation showed wide variability, with the alloy Ti-6Al-4V exhibiting the most consistent results. Accordingly, most of the analysis is concentrated upon data from that alloy, although creep and relaxation data from all the materials tested are presented. Experimental methods are outlined as well as some variations in methods of analysis. Various theoretical and practical implications of the work are discussed.

Zerwekh, R. P.↗
NASA NTRS1978

Finite elements and finite differences for transonic flow calculations

The paper reviews the chief finite difference and finite element techniques used for numerical solution of nonlinear mixed elliptic-hyperbolic equations governing transonic flow. The forms of the governing equations for unsteady two-dimensional transonic flow considered are the Euler equation, the full potential equation in both conservative and nonconservative form, the transonic small-disturbance equation in both conservative and nonconservative form, and the hodograph equations for the small-disturbance case and the full-potential case. Finite difference methods considered include time-dependent methods, relaxation methods, semidirect methods, and hybrid methods. Finite element methods include finite element Lax-Wendroff schemes, implicit Galerkin method, mixed variational principles, dual iterative procedures, optimal control methods and least squares.

Hafez, M. M.↗
NASA NTRS1984

Instrument failure detection and isolation in a system with variable plant parameters

This paper considers the problem of detecting instrument failures in a system in which there are plant parameter variations. The method described is based on the use of Kalman filters to provide a detection comparison in a manner similar to Clark's dedicated observer method. To improve the sensitivity to sensor variation and reduce the sensitivity to plant parameter variation the Kalman filters are designed to be sub-optimal. A procedure is described to design the sub-optimal filters so as to improve a measure of the ratio of sensor to plant sensitivity.

Alexandro, F. J., Jr.↗
NASA NTRS1989

On the accuracy of shape sensitivity derivatives

Six methods for calculating shape sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the finite-element mesh is refined, displacement sensitivity derivatives converge more slowly than the displacement themselves. Five of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer, rather than reduced accuracy.

Haftka, Raphael T.↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS2008

Signal Analysis Algorithms for Optimized Fitting of Nonresonant Laser Induced Thermal Acoustics Damped Sinusoids

This study seeks a numerical algorithm which optimizes frequency precision for the damped sinusoids generated by the nonresonant LITA technique. It compares computed frequencies, frequency errors, and fit errors obtained using five primary signal analysis methods. Using variations on different algorithms within each primary method, results from 73 fits are presented. Best results are obtained using an AutoRegressive method. Compared to previous results using Prony s method, single shot waveform frequencies are reduced approx.0.4% and frequency errors are reduced by a factor of approx.20 at 303K to approx. 0.1%. We explore the advantages of high waveform sample rates and potential for measurements in low density gases.

Balla, R. Jeffrey↗
NASA NTRS2022

NASA Engineering and Safety Center Technical Bulletin No. 22-04: Uncertainty Quantification of Reduced Order Structural Dynamic Models

Uncertainty quantification (UQ) provides statistical bounds on prediction accuracy based on finite element model (FEM) uncertainty. An alternate method for UQ, called the Hybrid Parametric Variation (HPV) combines a parametric variation of the Hurty/Craig-Bampton (HCB) fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method. This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available.

Uncertainty Quantification↗
NASA NTRS1993

Surface-induced brightness temperature variations and their effects on detecting thin cirrus clouds using IR emission channels in the 8-12 micrometer region

A method for detecting cirrus clouds in terms of brightness temperature differences between narrow bands at 8, 11, and 12 mu m has been proposed by Ackerman et al. (1990). In this method, the variation of emissivity with wavelength for different surface targets was not taken into consideration. Based on state-of-the-art laboratory measurements of reflectance spectra of terrestrial materials by Salisbury and D'Aria (1992), we have found that the brightness temperature differences between the 8 and 11 mu m bands for soils, rocks and minerals, and dry vegetation can vary between approximately -8 K and +8 K due solely to surface emissivity variations. We conclude that although the method of Ackerman et al. is useful for detecting cirrus clouds over areas covered by green vegetation, water, and ice, it is less effective for detecting cirrus clouds over areas covered by bare soils, rocks and minerals, and dry vegetation. In addition, we recommend that in future the variation of surface emissivity with wavelength should be taken into account in algorithms for retrieving surface temperatures and low-level atmospheric temperature and water vapor profiles.

Gao, Bo-Cai↗
NASA NTRS2000

A Simple Approach for Calibrating Imaging Systems with a Solid-State Sensor

Solid-state array imaging has become an indispensable tool for various modem scientific instrumentation and industrial applications. Often in experiments, use of an imaging system with linear response is very desirable. The linearity of commercial products, however, has yet to be accurately determined in order to assess their nonlinearity effects. Linearity is affected both by sensors and electronics. It is also very desirables to employ low-cost products if we can provide the performance comparable of expensive ones through calibration. Here, we present a very simple approach for approximately checking linearity of imaging systems and also a simple calibration method for precisely measuring intensity ratio or relative intensity variation. These methods can be readily implemented with existing equipment. The calibration method works well with an intermediate-grade analog sensor interfaced by a frame grabber. In our initial attempt for calibration, the nonlinearity of the camera can be reduced from about five percent to much less than one percent in standard deviation.

Leslie, Fred W.↗
NASA NTRS2010

Terahertz Mapping of Microstructure and Thickness Variations

A noncontact method has been devised for mapping or imaging spatial variations in the thickness and microstructure of a layer of a dielectric material. The method involves (1) placement of the dielectric material on a metal substrate, (2) through-the-thickness pulse-echo measurements by use of electromagnetic waves in the terahertz frequency range with a raster scan in a plane parallel to the substrate surface that do not require coupling of any kind, and (3) appropriate processing of the digitized measurement data.

Roth, Donald J.↗