Search NASASearch

SEARCH · Search NASA

Results for “Gradient estimation”

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 127 records · Page 7

Numerical experience with a class of algorithms for nonlinear optimization using inexact function and gradient information

For optimization problems associated with engineering design, parameter estimation, image reconstruction, and other optimization/simulation applications, low accuracy function and gradient values are frequently much less expensive to obtain than high accuracy values. Here, researchers investigate the computational performance of trust region methods for nonlinear optimization when high accuracy evaluations are unavailable or prohibitively expensive, and confirm earlier theoretical predictions when the algorithm is convergent even with relative gradient errors of 0.5 or more. The proper choice of the amount of accuracy to use in function and gradient evaluations can result in orders-of-magnitude savings in computational cost.

Carter, Richard G.

Constrained optimization using design of experiment surfaces

An algorithm for solving constrained optimization problems is presented. First, design of experiment techniques are used to survey the design space. After evaluating the objective and constraint functions, as specified by Taguchi orthogonal arrays, analytical models of these functions are generated using a least-squares regression analysis. Next, a nonlinear programming package is used to optimize the analytical model. Based on the optimization information, the design space is reduced so as to close in around the minimum, and the entire procedure is repeated until convergence. An important feature of the algorithm is that function gradients are not required; therefore, for problems in which gradients would have to be estimated using finite-differences the number of function evaluations required for the optimization is significantly reduced, when compared with traditional nonlinear programming techniques. In addition, there is no requirement that the gradients must be smooth and continuous.

Bolt, Marvin Vance

Thermoelastic response of metal matrix composites with large-diameter fibers subjected to thermal gradients

A new micromechanical theory is presented for the response of heterogeneous metal matrix composites subjected to thermal gradients. In contrast to existing micromechanical theories that utilize classical homogenization schemes in the course of calculating microscopic and macroscopic field quantities, in the present approach the actual microstructural details are explicitly coupled with the macrostructure of the composite. Examples are offered that illustrate limitations of the classical homogenization approach in predicting the response of thin-walled metal matrix composites with large-diameter fibers when subjected to thermal gradients. These examples include composites with a finite number of fibers in the thickness direction that may be uniformly or nonuniformly spaced, thus admitting so-called functionally gradient composites. The results illustrate that the classical approach of decoupling micromechanical and macromechanical analyses in the presence of a finite number of large-diameter fibers, finite dimensions of the composite, and temperature gradient may produce excessively conservative estimates for macroscopic field quantities, while both underestimating and overestimating the local fluctuations of the microscopic quantities in different regions of the composite. Also demonstrated is the usefulness of the present approach in generating favorable stress distributions in the presence of thermal gradients by appropriately tailoring the internal microstructure details of the composite.

Aboudi, Jacob

Fractionation of carbon isotopes by phytoplankton and estimates of ancient CO2 levels

Reports of the 13C content of marine particulate organic carbon are compiled and on the basis of GEOSECS data and temperatures, concentrations, and isotopic compositions of dissolved CO2 in the waters in which the related phytoplankton grew are estimated. In this way, the fractionation of carbon isotopes during photosynthetic fixation of CO2 is found to be significantly correlated with concentrations of dissolved CO2. Because ancient carbon isotopic fractionations have been determined from analyses of sedimentary porphyrins [Popp et al., 1989], the relationship between isotopic fractionation and concentrations of dissolved CO2 developed here can be employed to estimate concentrations of CO2 dissolved in ancient oceans and, in turn, partial pressures of CO2 in ancient atmospheres. The calculations take into account the temperature dependence of chemical and isotopic equilibria in the dissolved-inorganic-carbon system and of air-sea equilibria. Paleoenvironmental temperatures for each sample are estimated from reconstructions of paleogeography, latitudinal temperature gradients, and secular changes in low-latitude sea surface temperature. It is estimated that atmospheric partial pressures of CO2 were over 1000 micro atm 160 - 100 Ma ago, then declined to values near 300 micro atm during the next 100 Ma. Analysis of a high-resolution record of carbon isotopic fractionation at the Cenomanian-Turonian boundary suggests that the partial pressure of CO2 in the atmosphere was drawn down from values near 840 micro atm to values near 700 micro atm during the anoxic event.

NASA Discipline Exobiology

Flexural analysis of uplifted rift flanks on Venus

Knowledge of the thermal structure of a planet is vital to a thorough understanding of its general scheme of tectonics. Since no direct measurements of heat flow or thermal gradient are available for Venus, most estimates have been derived from theoretical considerations or by analog with the Earth. The flexural response of the lithosphere to applied loads is sensitive to regional thermal structure. Under the assumption that the yield strength as a function of depth can be specified, the temperature gradient can be inferred from the effective elastic plate thickness. Previous estimates of the effective elastic plate thickness of Venus range from 11-18 km for the foredeep north of Uorsar Rupes to 30-60 km for the annular troughs around several coronae. Thermal gradients inferred for these regions are 14-23 K km(exp -1) and 4-9 K km(exp -1) respectively. In this study, we apply the same techniques to investigate the uplifted flanks of an extensional rift. Hypotheses for the origin of uplifted rift flanks on Earth include lateral transport of heat from the center of the rift, vertical transport of heat by small-scale convection, differential thinning of the lithosphere, dynamical uplift, and isostatic response to mechanical uploading of the lithosphere. The 1st hypothesis is considered the dominant contributor to terrestrial rift flanks lacking evidence for volcanic activity, particularly for rift structures that are no longer active. In this study, we model the uplifted flanks of a venusian rift as the flexural response to a vertical end load.

Evans, Susan A.

Atmospheric Modeling Using Accelerometer Data During Mars Atmosphere and Volatile Evolution (MAVEN) Flight Operations

The Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft is a NASA orbiter designed to explore the Mars upper atmosphere, typically from 140 to 160 km altitude. In addition to the nominal science mission, MAVEN has performed several Deep Dip campaigns in which the orbit's closest point of approach, also called periapsis, was lowered to an altitude range of 115 to 135 km. MAVEN accelerometer data were used during mission operations to estimate atmospheric parameters such as density, scale height, along-track gradients, and wave structures. Density and scale height estimates were compared against those obtained from the Mars Global Reference Atmospheric Model and used to aid the MAVEN navigation team in planning maneuvers to raise and lower periapsis during Deep Dip operations. This paper describes the processes used to reconstruct atmosphere parameters from accelerometers data and presents the results of their comparison to model and navigation-derived values.

Tolson, Robert H.

Estimating tropical cyclone central pressure and outer winds from satellite microwave data

A new technique is presented for estimating central pressure and surface winds outside of the radius of maximum winds from brightness temperatures measured by the SCAMS on Nimbus 6. The problem of estimating maximum surface winds is not treated because the spatial resolution of the radiometer is not high enough for this purpose. The technique is developed on the basis of measurements over eight typhoons and five hurricanes during 1975. It is shown that the warm 55.45-GHz brightness temperature anomaly previously observed over Typhoon June can be found over other tropical cyclones. Maximum brightness temperature anomaly correlates with central pressure with a correlation coefficient of -0.859. Outer surface wind speeds are calculated by assuming gradient balance and using the regression between brightness temperature anomaly and central pressure to estimate pressure gradients.

Kidder, S. Q.

Potential vorticity structure of Titan's polar vortices from Cassini CIRS observations

The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.

Conor Andrew Nixon

Estimation of optical flow in airborne electro-optical sensors by stochastic approximation

The essence of motion or range estimation by passive electrooptical means is the ability to determine the correspondence of picture elements in pairs of image frames and to estimate their coordinates and their disparity (relative shifts) in the image plane of an electrooptical imaging sensor. The disparity can be in successive frames due to self-motion or in simultaneous frames of a stereo pair. A key issue is to provide these estimates on-line. This paper describes the theoretical background of such an interframe shift estimator. It is based on a stochastic gradient algorithm, specifically implementing a form of stochastic approximation, which can achieve rapid convergence of the shift estimate. Analytical and numerical simulation examples for random texture and isolated features validate the feasibility and the effectiveness of the estimator.

Merhav, S. J.

A variation of the Davis-Smith method for in-flight determination of spacecraft magnetic fields.

A variation of a procedure developed by Davis and Smith (1968) is presented for the in-flight determination of spacecraft magnetic fields. Both methods take statistical advantage of the observation that fluctuations in the interplanetary magnetic field over short periods of time are primarily changes in direction rather than in magnitude. During typical solar wind conditions between 0.8 and 1.0 AU, a statistical analysis of 2-3 days of continuous interplanetary field measurements yields an estimate of a constant spacecraft field with an uncertainty of plus or minus 0.25 gamma in the direction radial to the sun and plus or minus 15 gammas in the directions transverse to the radial. The method is also of use in estimating variable spacecraft fields with gradients of the order of 0.1 gamma/day and less and in other special circumstances.

Belcher, J. W.

Submillimeter observations of W3

A 35 arcsec resolution, 400 micron continuum map of the core of the W3 molecular cloud complex is presented. The map shows prominent peaks centered approximately at the positions of the infrared sources IRS 4 and IRS 5 with stronger-emission at the western (IRS 4) peak. The visual extinction (inferred from the submillimeter measurements) through the center of both peaks is large (A sub v = 200). A total mass for the core of 900 solar masses is estimated, consistent with dynamical estimates if the observed velocity gradient is a result of rotation. The peak H2 density is approximately 2.8 x 10 to the 5th/cu cm. The observations exclude the possibility that the shell H II regions, W3A-C, could be in pressure equilibrium with the surrounding cloud. The observations establish that the core of W3 contains several massive condensations surrounding the youngest stars in the region.

Jaffe, D. T.

Validation of temperature retrievals obtained by the Limb Infrared Monitor of the Stratosphere (LIMS) experiment on Nimbus 7

An outline is given of LIMS temperature determinations (as a function of pressure) from measurements in two channels covering portions of the 15-micrometer band of carbon dioxide. The known sources of error from the radiometer and data reduction are used to estimate the systematic and random errors expected of the results. Observational determinations of the complete end-to-end precision are obtained by computing the standard deviation of six sequential temperature retrievals in regions where the atmosphere is horizontally uniform. This yields values of 0.2 to 0.6 K, in reasonable agreement with the estimates. A correction for horizontal gradients in the atmosphere leads to a large reduction in the differences between the stratospheric temperatures determined on the ascending and descending portions of the orbit. The temperatures agree in the mean with radiosondes and rocketsondes to within 1-2 K in most regions below 1 mbar. Several interesting, previously unseen features included cold regions in the mid-latitude mesosphere and wavelike vertical variations in the tropics.

Gille, J. C.

Navier-Stokes analysis of two- and three-dimensional flow field in solid rocket motors with segment joints

A multidimensional implicit Navier-Stokes analysis which uses numerical solution of ensemble-averaged Navier-Stokes equations in a nonorthogonal bodyfitted cylindrical-polar coordinate system has been applied to simulation of the internal flow field in solid-propellant rocket motor chambers with segment joints. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. Computations performed to simulate the axisymmetric flow field in the vicinity of the aft field joint in the Space Shuttle SRB using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot. Calculations were also performed for two cases involving asymmetric three-dimensional flow in the vicinity of the aft field joint in the SRB using 721,525 grid points to estimate circumferential velocities and pressure gradients at the joint.

Sabnis, J. S.

Analysis of Variation of Piston Temperature with Piston Dimensions and Undercrown Cooling

A theoretical analysis is presented that permits estimation of the changes in piston-temperature distribution induced by variations in the crown thickness, the ring-groove-pad thickness, and the undercrown surface heat-transfer coefficient. The analysis consists of the calculation of operating temperatures at various points in the piston body on the basis of the experimentally determined surface heat-transfer coefficients and boundary-region temperatures, as well as arbitrarily selected surface coefficients. Surface heat-transfer coefficients were estimated from the internal temperature gradients obtained by hardness surveys of aluminum pistons that had been operated under severe conditions in a liquid-cooled, single-cylinder, 5 1/2 by 6-inch test engine.

Sanders, J C

A High Order Finite Difference Scheme with Sharp Shock Resolution for the Euler Equations

We derive a high-order finite difference scheme for the Euler equations that satisfies a semi-discrete energy estimate, and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semi-discrete energy estimate is based on a symmetrization of the Euler equations that preserves the homogeneity of the flux vector, a canonical splitting of the flux derivative vector, and the use of difference operators that satisfy a discrete analogue to the integration by parts procedure used in the continuous energy estimate. Around discontinuities or sharp gradients, refined grids are created on which the discrete equations are solved after adding a newly constructed artificial viscosity. The positioning of the sub-grids and computation of the viscosity are aided by a detection algorithm which is based on a multi-scale wavelet analysis of the pressure grid function. The wavelet theory provides easy to implement mathematical criteria to detect discontinuities, sharp gradients and spurious oscillations quickly and efficiently.

Gerritsen, Margot

Assessment of Error in Synoptic-Scale Diagnostics Derived from Wind Profiler and Radiosonde Network Data

A topic of current practical interest is the accurate characterization of the synoptic-scale atmospheric state from wind profiler and radiosonde network observations. We have examined several related and commonly applied objective analysis techniques for performing this characterization and considered their associated level of uncertainty both from a theoretical and a practical standpoint. A case study is presented where two wind profiler triangles with nearly identical centroids and no common vertices produced strikingly different results during a 43-h period. We conclude that the uncertainty in objectively analyzed quantities can easily be as large as the expected synoptic-scale signal. In order to quantify the statistical precision of the algorithms, we conducted a realistic observing system simulation experiment using output from a mesoscale model. A simple parameterization for estimating the uncertainty in horizontal gradient quantities in terms of known errors in the objectively analyzed wind components and temperature is developed from these results.

Mace, Gerald G.