Search NASA⌕ Search

SEARCH · Search NASA

Results for “VARIATION 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 271 records · Page 15

Localized basis functions and other computational improvements in variational nonorthogonal basis function methods for quantum mechanical scattering problems involving chemical reactions

The Generalized Newton Variational Principle for 3D quantum mechanical reactive scattering is briefly reviewed. Then three techniques are described which improve the efficiency of the computations. First, the fact that the Hamiltonian is Hermitian is used to reduce the number of integrals computed, and then the properties of localized basis functions are exploited in order to eliminate redundant work in the integral evaluation. A new type of localized basis function with desirable properties is suggested. It is shown how partitioned matrices can be used with localized basis functions to reduce the amount of work required to handle the complex boundary conditions. The new techniques do not introduce any approximations into the calculations, so they may be used to obtain converged solutions of the Schroedinger equation.

Schwenke, David W.↗

Pulse-Echo Ultrasonic Imaging Method for Eliminating Sample Thickness Variation Effects

A pulse-echo, immersion method for ultrasonic evaluation of a material which accounts for and eliminates nonlevelness in the equipment set-up and sample thickness variation effects employs a single transducer and automatic scanning and digital imaging to obtain an image of a property of the material, such as pore fraction. The nonlevelness and thickness variation effects are accounted for by pre-scan adjustments of the time window to insure that the echoes received at each scan point are gated in the center of the window. This information is input into the scan file so that, during the automatic scanning for the material evaluation, each received echo is centered in its time window. A cross-correlation function calculates the velocity at each scan point, which is then proportionalized to a color or grey scale and displayed on a video screen.

Roth, Don J.↗

Pulse-echo ultrasonic imaging method for eliminating sample thickness variation effects

A pulse-echo, immersion method for ultrasonic evaluation of a material is discussed. It accounts for and eliminates nonlevelness in the equipment set-up and sample thickness variation effects employs a single transducer, automatic scanning and digital imaging to obtain an image of a property of the material, such as pore fraction. The nonlevelness and thickness variation effects are accounted for by pre-scan adjusments of the time window to insure that the echoes received at each scan point are gated in the center of the window. This information is input into the scan file so that, during the automatic scanning for the material evaluation, each received echo is centered in its time window. A cross-correlation function calculates the velocity at each scan point, which is then proportionalized to a color or grey scale and displayed on a video screen.

Roth, Don J.↗

A multivariate variational objective analysis-assimilation method. Part 2: Case study results with and without satellite data

The variational multivariate assimilation method described in a companion paper by Achtemeier and Ochs is applied to conventional and conventional plus satellite data. Ground-based and space-based meteorological data are weighted according to the respective measurement errors and blended into a data set that is a solution of numerical forms of the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation for a dry atmosphere. The analyses serve first, to evaluate the accuracy of the model, and second to contrast the analyses with and without satellite data. Evaluation criteria measure the extent to which: (1) the assimilated fields satisfy the dynamical constraints, (2) the assimilated fields depart from the observations, and (3) the assimilated fields are judged to be realistic through pattern analysis. The last criterion requires that the signs, magnitudes, and patterns of the hypersensitive vertical velocity and local tendencies of the horizontal velocity components be physically consistent with respect to the larger scale weather systems.

Achtemeier, Gary L.↗

Comparison of Artificial Compressibility Methods for Coupled Laminar Fluid Flow and Heat Transfer

Artificial compressibility (AC) methods for incompressible flows perturb the mass continuity equation by the time derivative of the pressure field to achieve direct coupling between pressure and velocity without the need for a projection method. Several variations of the original AC method have been proposed (e.g., the entropically damped artificial compressibility [EDAC] method) that can, in some scenarios, improve on the original model. In this paper, we methodically compare the results of the original AC model and the EDAC model for both steady-state and transient problems with coupling between the momentum and energy equations. We demonstrate that the EDAC model and the AC model yield similar numerical results, and the EDAC model is 10%-20% faster than the AC model for a given problem.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗

Development of New Modeling and Analysis Tools for Solar Sails

Existing finite-element-based structural analysis codes are ineffective in treating deployable gossamer space systems, including solar sails that are formed by long space-deployable booms and extremely large thin-film membrane apertures. Recognizing this, the NASA Space transportation Technology Program has initiated and sponsored a focused research effort to develop new and computationally efficient structural modeling and analysis tools for solar sails. The technical approach of this ongoing effort will be described. Two solution methods, the Distributed Transfer Function Method and the Parameter-Variation-Principle method, based on which the technical approach was formatted are also discussed.

boom↗

Variational Path Sampling of Rare Dynamical Events

This article reviews the concepts and methods of variational path sampling. These methods allow computational studies of rare events in systems driven arbitrarily far from equilibrium. Based upon a statistical mechanics of trajectory space and leveraging the theory of large deviations, they provide a perspective from which dynamical phenomena can be studied with the same types of ensemble reweighting ideas that have been used for static equilibrium properties. Applications to chemical, material, and biophysical systems are highlighted.

Singh, Aditya N↗

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

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

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↗

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

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

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

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

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

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