Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Gradient Calculation Methods on Arbitrary Polyhedral Unstructured Meshes for Cell-Centered CFD Solvers

A survey of gradient reconstruction methods for cell-centered data on unstructured meshes is conducted within the scope of accuracy assessment. Formal order of accuracy, as well as error magnitudes for each of the studied methods, are evaluated on a complex mesh of various cell types through consecutive local scaling of an analytical test function. The tests highlighted several gradient operator choices that can consistently achieve 1st order accuracy regardless of cell type and shape. The tests further offered error comparisons for given cell types, leading to the observation that the "ideal" gradient operator choice is not universal. Practical implications of the results are explored via CFD solutions of a 2D inviscid standing vortex, portraying the discretization error properties. A relatively naive, yet largely unexplored, approach of local curvilinear stencil transformation exhibited surprisingly favorable properties

Meshes↗
NASA NTRS1987

Accurate pressure gradient calculations in hydrostatic atmospheric models

A method for the accurate calculation of the horizontal pressure gradient acceleration in hydrostatic atmospheric models is presented which is especially useful in situations where the isothermal surfaces are not parallel to the vertical coordinate surfaces. The present method is shown to be exact if the potential temperature lapse rate is constant between the vertical pressure integration limits. The technique is applied to both the integration of the hydrostatic equation and the computation of the slope correction term in the horizontal pressure gradient. A fixed vertical grid and a dynamic grid defined by the significant levels in the vertical temperature distribution are employed.

Carroll, John J.↗
NASA NTRS1991

Airfoil optimization with efficient gradient calculations

The viscous airfoil design analysis code XFOIL was extended to allow optimization using conformal mapping coefficients as design variables. The optimization technique used was the Steepest Descent method applied to a Penalty Function. The gradients of the aerodynamic variables with respect to the design variables were cheaply calculated as by-products of XFOIL's integral boundary layer Newton solver. The speed of the optimization process further increased by updating the Newton system boundary layer variables after each optimization step using the available gradient information. Two examples are presented.

Sorensen, Thomas↗
NASA NTRS1997

Arterial Pressure Gradients during Upright Posture and 30 deg Head Down Tilt

Gravity alters local blood pressure within the body so that arterial pressures in the head and foot are lower and higher, respectively, than that at heart level. Furthermore, vascular responses to local alterations of arterial pressure are probably important to maintain orthostatic tolerance upon return to the Earth after space flight. However, it has been difficult to evaluate the body's arterial pressure gradient due to the lack of noninvasive technology. This study was therefore designed to investigate whether finger arterial pressure (FAP), measured noninvasively, follows a normal hydrostatic pressure gradient above and below heart level during upright posture and 30 deg head down tilt (HDT). Seven healthy subjects gave informed consent and were 19 to 52 years old with a height range of 158 to 181 cm. A Finapres device measured arterial pressure at different levels of the body by moving the hand from 36 cm below heart level (BH) to 72 cm above heart level (AH) in upright posture and from 36 cm BH to 48 cm AH during HDT in increments of 12 cm. Mean FAP creased by 85 mmHg transitioning from BH to AH in upright posture, and the pressure gradient calculated from hydrostatic pressure difference (rho(gh)) was 84 mmHg. In HDT, mean FAP decreased by 65 mmHg from BH to AH, and the calculated pressure gradient was also 65 mmHg. There was no significant difference between the measured FAP gradient and the calculated pressure gradient, although a significant (p = 0.023) offset was seen for absolute arterial pressure in upright posture. These results indicate that arterial pressure at various levels can be obtained from the blood pressure at heart level by calculating rho(gh) + an offset. The offset equals the difference between heart level and the site of measurement. In summary, we conclude that local blood pressure gradients can be measured by noninvasive studies of FAP.

Sanchez, E. R↗
NASA NTRS1983

Computational methods for aerodynamic design using numerical optimization

Five methods to increase the computational efficiency of aerodynamic design using numerical optimization, by reducing the computer time required to perform gradient calculations, are examined. The most promising method consists of drastically reducing the size of the computational domain on which aerodynamic calculations are made during gradient calculations. Since a gradient calculation requires the solution of the flow about an airfoil whose geometry was slightly perturbed from a base airfoil, the flow about the base airfoil is used to determine boundary conditions on the reduced computational domain. This method worked well in subcritical flow.

Peeters, M. F.↗
NASA NTRS2003

Optimization Based Efficiencies in First Order Reliability Analysis

This paper develops a method for updating the gradient vector of the limit state function in reliability analysis using Broyden's rank one updating technique. In problems that use commercial code as a black box, the gradient calculations are usually done using a finite difference approach, which becomes very expensive for large system models. The proposed method replaces the finite difference gradient calculations in a standard first order reliability method (FORM) with Broyden's Quasi-Newton technique. The resulting algorithm of Broyden updates within a FORM framework (BFORM) is used to run several example problems, and the results compared to standard FORM results. It is found that BFORM typically requires fewer functional evaluations that FORM to converge to the same answer.

Peck, Jeffrey A.↗
NASA NTRSDate not supplied

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗
NASA NTRS1988

Mixing rates calculated from potential vorticity

Horizontal mixing coefficients K(yy) and K(yz) were computed using potential vorticity and gradients calculated from NMC data for geopotential heights and temperatures. It was found that the mid-latitude to polar mixing coefficients depend on the use of nongeostrophic winds, which yield K(yy) values which are mostly positive, whereas geostophic winds yield K(yy) values which were largely negative near the polar region. A comparison of the K(yy) values calculated from the potential vorticity fluxes and gradients with the K(yy) values calculated using the momentum balance between the flux of potential vorticity and the residual circulation has shown that the residual-circulation-balanced K(yy) values in the upper stratosphere were much larger than the values calculated by the present method. It is suggested that this difference is the result of breaking gravity waves in the upper stratosphere and lower mesosphere.

Newman, Paul A.↗
NASA NTRS1984

SSME thrust chamber simulation using Navier-Stokes equations

The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.

Przekwas, A. J.↗
NASA NTRS2016

Full Gradient Solution to Adaptive Hybrid Control

This paper focuses on the adaptation mechanisms in adaptive hybrid controllers. Most adaptive hybrid controllers update two filters individually according to the filtered-reference least mean squares (FxLMS) algorithm. Because this algorithm was derived for feedforward control, it does not take into account the presence of a feedback loop in the gradient calculation. This paper provides a derivation of the proper weight vector gradient for hybrid (or feedback) controllers that takes into account the presence of feedback. In this formulation, a single weight vector is updated rather than two individually. An internal model structure is assumed for the feedback part of the controller. The full gradient is equivalent to that used in the standard FxLMS algorithm with the addition of a recursive term that is a function of the modeling error. Some simulations are provided to highlight the advantages of using the full gradient in the weight vector update rather than the approximation.

Bean, Jacob↗
NASA NTRS2017

Full Gradient Solution to Adaptive Hybrid Control

This paper focuses on the adaptation mechanisms in adaptive hybrid controllers. Most adaptive hybrid controllers update two filters individually according to the filtered reference least mean squares (FxLMS) algorithm. Because this algorithm was derived for feedforward control, it does not take into account the presence of a feedback loop in the gradient calculation. This paper provides a derivation of the proper weight vector gradient for hybrid (or feedback) controllers that takes into account the presence of feedback. In this formulation, a single weight vector is updated rather than two individually. An internal model structure is assumed for the feedback part of the controller. The full gradient is equivalent to that used in the standard FxLMS algorithm with the addition of a recursive term that is a function of the modeling error. Some simulations are provided to highlight the advantages of using the full gradient in the weight vector update rather than the approximation.

Bean, Jacob↗
NASA NTRS1986

An analysis of mesoscale VAS retrievals using statistical structure functions

Mesoscale sounding data from a truth rawinsonde network and a VAS from March 6-7, 1982 are statistically studied. The parameters of temperature, geopotential height, and mixing ratio at 850, 700, 500, 300, and 200 mb levels, and thickness and precipitable water are examined. The VAS retrieval gradients derived from a physical algorithm and from a regression technique are analyzed; the analysis reveals that both techniques have similar error characteristics and capabilities for calculating gradients of mesoscale temperature and geopotential height. Ground-truth rawinsonde data and VAS retrievals are compared; it is observed that VAS data describes weaker gradients than rawinsonde data.

Fuelberg, H. E.↗
NASA NTRS2005

Traj_opt User's Guide

Trajectory optimization program Traj_opt was developed at Ames Research Center to help assess the potential benefits of ultrahigh temperature ceramic materials applied to reusable space vehicles with sharp noses and wing leading edges. Traj_opt loosely couples the Ames three-degrees-of-freedom trajectory package Traj (see NASA-TM-2004-212847) with the SNOPT optimization package (Stanford University Technical Report SOL 98-1). Traj_opt version January 22, 2003 is covered by this user guide. The program has been applied extensively to entry and ascent abort trajectory calculations for sharp and blunt crew transfer vehicles. The main optimization variables are control points for the angle of attack and bank angle time histories. No propulsion options are provided, but numerous objective functions may be specified and the nonlinear constraints implemented include a distributed surface heating constraint capability. Aero-capture calculations are also treated with an option to minimize orbital eccentricity at apoapsis. Traj_opt runs efficiently on a single processor, using forward or central differences for the gradient calculations. Results may be displayed conveniently with Gnuplot scripts. Control files recommended for five standard reentry and ascent abort trajectories are included along with detailed descriptions of the inputs and outputs.

Saunders, David A.↗
NASA NTRS1994

Multidisciplinary design optimization using genetic algorithms

Multidisciplinary design optimization (MDO) is an important step in the conceptual design and evaluation of launch vehicles since it can have a significant impact on performance and life cycle cost. The objective is to search the system design space to determine values of design variables that optimize the performance characteristic subject to system constraints. Gradient-based optimization routines have been used extensively for aerospace design optimization. However, one limitation of gradient based optimizers is their need for gradient information. Therefore, design problems which include discrete variables can not be studied. Such problems are common in launch vehicle design. For example, the number of engines and material choices must be integer values or assume only a few discrete values. In this study, genetic algorithms are investigated as an approach to MDO problems involving discrete variables and discontinuous domains. Optimization by genetic algorithms (GA) uses a search procedure which is fundamentally different from those gradient based methods. Genetic algorithms seek to find good solutions in an efficient and timely manner rather than finding the best solution. GA are designed to mimic evolutionary selection. A population of candidate designs is evaluated at each iteration, and each individual's probability of reproduction (existence in the next generation) depends on its fitness value (related to the value of the objective function). Progress toward the optimum is achieved by the crossover and mutation operations. GA is attractive since it uses only objective function values in the search process, so gradient calculations are avoided. Hence, GA are able to deal with discrete variables. Studies report success in the use of GA for aircraft design optimization studies, trajectory analysis, space structure design and control systems design. In these studies reliable convergence was achieved, but the number of function evaluations was large compared with efficient gradient methods. Applicaiton of GA is underway for a cost optimization study for a launch-vehicle fuel-tank and structural design of a wing. The strengths and limitations of GA for launch vehicle design optimization is studied.

Unal, Resit↗
NASA NTRS1971

Study to determine and improve design for lithium-doped solar cells

Solar cell experiments show that a single lithium density parameter, the lithium density gradient, calculated from nondestructive capacitance measurements, provides the basis for accurate predictions of lithium cell behavior in a 1-MeV electron environment for fluences ranging between 3 X 10 to the 13th power e/sq cm and 3 X 10 to the 15th power/e sq cm. The oxygen-rich (quartz crucible) lithium cell with phosphorous starting dopant and lithium gradient between approximately 5 X 10 to the 18th power and 1.5 x 10 to the 19th power/cm to the 4th power was found superior in performance to the commercial 10 ohm-cm n/p control cells. Post-recovery stability of oxygen-rich cells was satisfactory. An average post-recovery current drop of approximately 1 mA was observed for 70 crucible cells after 1 year-equivalent storage time at 80 C. In contrast the oxygen-poor (float zone and Lopex) lithium cells displayed spotty initial performance and stability problems at room temperature.

Brucker, G.↗
NASA NTRS1990

Shape design sensitivities using fully automatic 3-D mesh generation

Previous work in three dimensional shape optimization involved specifying design variables by associating parameters directly with mesh points. More recent work has shown the use of fully-automatic mesh generation based upon a parameterized geometric representation. Design variables have been associated with a mathematical model of the part rather than the discretized representation. The mesh generation procedure uses a nonuniform grid intersection technique to place nodal points directly on the surface geometry. Although there exists an associativity between the mesh and the geometrical/topological entities, there is no mathematical functional relationship. This poses a problem during certain steps in the optimization process in which geometry modification is required. For the large geometrical changes which occur at the beginning of each optimization step, a completely new mesh is created. However, for gradient calculations many small changes must be made and it would be too costly to regenerate the mesh for each design variable perturbation. For that reason, a local remeshing procedure has been implemented which operates only on the specific edges and faces associated with the design variable being perturbed. Two realistic design problems are presented which show the efficiency of this process and test the accuracy of the gradient computations.

Botkin, M. E.↗