Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient approximations”

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

Implementation of Vigneron's streamwise pressure gradient approximation in the PNS equations

For single sweep parabolized Navier-Stokes solvers, the streamwise pressure gradient must be modified in the subsonic region to eliminate numerical instabilities. The accuracy of this modification on the solution of the parabolized Navier-Stokes equations with Vigneron's technique is shown to depend on how the numerical approximation of the pressure gradient is formed. A simple test case of supersonic laminar flow over a flat plate is computed with two different numerical methods for solving the PNS equations. Significant errors in the temperature profile and skin friction coefficient are demonstrated using a fully conservative differencing treatment of Vigneron's splitting for the pressure gradient typically used in parabolized Navier-Stokes solvers. The physical reason for this error is discussed. An alternate formulation is demonstrated which minimizes these errors.

Morrison, Joseph H.↗

Galileo/NIMS near-infrared thermal imagery of the surface of Venus

Numerous highland and lowland features on the surface of Venus are observed in multispectral imagery acquired at approximately 50 km spatial resolution by the Near-Infrared Mapping Spectrometer (NIMS) on board the Galileo spacecraft in Feb. 1990. Specifically, such features are observed at 1.18 microns, a wavelength particularly sensitive to thermal emission from the hot, lower atmosphere (less than 10 km) and surface, and show up particularly well when this image is 'de-clouded' using a simultaneously-acquired 2.3-microns image of the upper, cloudy atmosphere. Due to the steep atmospheric temperature gradient (approximately 8 degrees per kilometer), hot lowland areas appear relatively bright, while cooler, highland areas appear dark (due to the steep atmospheric temperature gradient - approximately 8 degrees per kilometer - surface temperatures span approximately 100 K over the 13 kilometer range of surface altitudes observed in this image). Prominent highland features include Maxwell Montes (approximately 12 km altitude), Alpha Regio (2.5 km), Eistla Regio (approximately 2.0 km), Bell Regio (2-3 km), and the western edge of Aphrodite Terra (2-2.5 km). Low-lying regions include Sedna Planitia (-1.0 km), Tinatin Planitia (-0.5 km), and the Bereghinya Planitia (0 km). From correlations with radar altimetry maps, such imagery may place useful constraints on surface emissivity and temperature variations, as well as on the nature of continuum opacity of CO2 in the 1-micron region.

Carlson, R. W.↗

Numerical optimization in Hilbert space using inexact function and gradient evaluations

Trust region algorithms provide a robust iterative technique for solving non-convex unstrained optimization problems, but in many instances it is prohibitively expensive to compute high accuracy function and gradient values for the method. Of particular interest are inverse and parameter estimation problems, since function and gradient evaluations involve numerically solving large systems of differential equations. A global convergence theory is presented for trust region algorithms in which neither function nor gradient values are known exactly. The theory is formulated in a Hilbert space setting so that it can be applied to variational problems as well as the finite dimensional problems normally seen in trust region literature. The conditions concerning allowable error are remarkably relaxed: relative errors in the gradient error condition is automatically satisfied if the error is orthogonal to the gradient approximation. A technique for estimating gradient error and improving the approximation is also presented.

Carter, Richard G.↗

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗

Warm-Core Intensification Through Horizontal Eddy Heat Transports into the Eye

The mechanism for the formation and intensification of the hurricane warm core is not well understood. The generally accepted explanation is that the warm core forms as a result of gentle subsidence of air within the eye that warms as a result of adiabatic compression. Malkus suggested that this subsidence is part of a deep circulation in which air begins descent at high levels in the eye, acquires cyclonic angular momentum as it descends to lower levels, and then diverges at low levels, where it is entrained back into the eyewall. Inward mixing from the eyewall is hypothesized to force the subsidence and maintain the moisture and momentum budgets of the subsiding air. Willoughby suggested that air within the eye has remained so since it was first enclosed during the formation of the eyewall and that it subsides at most only a few kilometers rather than through the depth of the troposphere. He relates the subsidence to the low-level divergence and entrainment into the eyewall noted by Malkus, but suggests that shrinkage of the eye's volume is more than adequate to account for the air lost to the eyewall or converted to cloudy air by turbulent mixing across the eye boundary. Smith offered an alternative view of the subsidence forcing, suggesting that vertical motion in a mature hurricane eye is generated largely by imbalances between the downward vertical pressure gradient force and the upward buoyancy force. The vertical pressure gradient force is associated with the decay and/or radial spread of the tangential wind field with height at those levels were the winds are in approximate gradient wind balance. The rate of subsidence is just that required to warm the air sufficiently such that the buoyancy remains in close hydrostatic balance with an increasing vertical pressure gradient force. In this study, a very high-resolution simulation of Hurricane Bob using a cloud-resolving grid scale of 1.3 km is used to examine the heat budget within the storm with particular emphasis on the mechanisms for warming of the eye.

Braun, Scott A.↗

Galactic cosmic ray gradients, field-aligned and latitudinal, among Voyagers 1/2 and IMP-8

The present investigation represents a summary of a comprehensive analysis of the same subject conducted by Roelof et al. (1981). It is pointed out that the tandem earth-Jupiter trajectories of the Voyager 1/2 spacecraft, combined with baseline measurements from the earth-orbiting IMP 7/8 spacecraft, provide the first opportunity for unambiguously separating latitude from radial or field-aligned effects in galactic cosmic ray gradients. Attention is given to the method of data analysis, and the separation of field-aligned and latitudinal gradients. It is found that latitudinal gradients approximately equal to or greater than 1 percent per deg in the cosmic ray intensity were a common feature of the interplanetary medium between 1 and 5 AU in 1977-78. Except in the most disturbed periods, cosmic ray intensities are well-ordered in field-aligned structures.

Roelof, E. C.↗

High-resolution submillimeter-wave radiometry of supersonic flow

The recent development of a high-resolution submillimeter-wave heterodyne radiometer has made possible the first measurements of H2O molecule rotational line excitation temperatures and detailed profiles in supersonic flow. Absorption signals were measured across the flow for the 2/11/ from 2//02/ (752 GHz) para-H2O rotational transition against a hot background. These signals decrease downstream owing to the volume expansion of the gas away from the sonic nozle exit in the high-vacuum chamber. Radiative transfer calculations based on the large-velocity-gradient approximation and multilevel statistical equilibrium agree with these results and with the measured spectral line shapes. The data reveal nearly isentropic gas expansion and cooling. These studies have shown that submillimeter-wave heterodyne radiometry can be useful for remote sensing of supersonic flow with low mass flux, provided the signal transmission is through a dry or thin atmosphere.

Dionne, G. F.↗

Density-functional calculations of the surface tension of liquid Al and Na

Calculations of the surface tensions of liquid Al and Na are described using the full ionic density functional formalism of Wood and Stroud (1983). Surface tensions are in good agreement with experiment in both cases, with results substantially better for Al than those found previously in the gradient approximation. Preliminary minimization with respect to surface profile leads to an oscillatory profile superimposed on a nearly steplike ionic density disribution; the oscillations have a wavellength of about a hardsphere diameter.

Stroud, D.↗

C3H2 observations as a diagnostic probe for molecular clouds

Recently the three-membered ring molecule, cyclopropenylidene, C3H2, has been identified in the laboratory and detected in molecular clouds by Thaddeus, Vrtilek and Gottlieb (1985). This molecule is wide-spread throughout the Galaxy and has been detected in 25 separate sources including cold dust clouds, circumstellar envelopes, HII regions, and the spiral arms observed against the Cas supernova remnant. In order to evaluate the potential of C3H2 as a diagnostic probe for molecular clouds, and to attempt to identify the most useful transitions, statistical equilibrium calculations were carried out for the lowest 24 levels of the ortho species and the lowest 10 levels of the para species. Many of the sources observed by Matthews and Irvine (1985) show evidence of being optically thick in the 1(10)-1(01) line. Consequently, the effects of radiative trapping should be incorporated into the equilibrium calculations. This was done using the Large Velocity Gradient approximation for a spherical cloud of uniform density. Some results of the calculations for T(K)=10K are given. Figures are presented which show contours of the logarithm of the ratio of peak line brightness temperatures for ortho-para pairs of lines at similar frequencies. It appears that the widespread nature of C3H2, the relatively large strength of its spectral lines, and their sensitivity to density and molecular abundance combine to make this a useful molecule for probing physical conditions in molecular clouds. The 1(10)-1(01) and 2(20)-2(11) K-band lines may be especially useful in this regard because of the ease with which they are observed and their unusual density-dependent emission/absorption properties.

Avery, L. W.↗

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Canuto, V. M.↗

New Third-Order Moments for the PBL

Turbulent convection is inherently non-local and a primary condition for a successful treatment of the PBL is a reliable model of non-locality. In the dynamic equations governing the convective flux, turbulent kinetic energy, etc, non-locality enters through the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA , severely underestimates the TOMs (by up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use. Here, we present a new analytic expression for the TOMs which is considerably simpler than the 1994 expression and which at the same time yields a much better fit to the LES data.

Canuto, V. M.↗

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗

Three Isotopes of Oxygen in Lunar Samples - The Same as Earth or Different?

One of the most important measurements that can be made for any extraterrestrial sample is determination of the relative abundance of O-16, O-17 and O-18. To make the comparison, investigators report their results as delta (δ) values compared to terrestrial delta O-17 vs delta O-18 for a suite of specimens because for any given reservoir, i.e. a planetary parent body such as the Earth or Mars, the numbers should define a line of gradient approximately one half. Different bodies should be identified from their delta O-17 which has traditionally been defined by the equation delta O-17 = delta O-17 - (a constant, ca. 0.52) × delta O-18.

Pillinger, C. T.↗

A rotor optimization using regression analysis

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Giansante, N.↗

Differential Equation Approximation Using Gradient-Boosted Quantile Regression

The operation of cyber-physical-human (CPH) systems is subject to various epistemic and aleatory uncertainties. Overall trustworthiness of CPH systems relies on the trustworthiness of its components and their interactions. It is important that computational models comprising the cyber component of CPH provide predictions accompanied by a measure of confidence in model outcomes. Uncertainty quantification (UQ) and propagation are especially important in safety critical CPH systems. Gradient-boosted trees is a modeling approach capable both of learning the dynamics of a system and performing UQ. In this paper, we devise a method for using gradient boosting to learn the dynamics of a second order differential equation and estimate uncertainty at the same time. We do this by creating a custom loss function that trains the model to approximate the second derivative of a noisy time series, and to penalize based on a parameter that corresponds to the desired quantile. The resulting gradient boosting model can simulate stochastic trajectories of the system given a single starting point, that is, it can estimate both the expected trajectory and its uncertainty. We show that the uncertainty estimation is well calibrated and that the model can learn the dynamics even in the presence of noise. We demonstrate the approach on a simple cartpole system.

Autonomous systems↗

The Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

An approximate method for the calculation of the compressible laminar boundary layer with heat transfer and arbitrary pressure gradient, based on Thwaites' correlation concept, is presented. With the definition of dimensionless shear and heat-transfer parameters and an assumed correlation of these parameters in terms of a momentum parameter, a complete system of relations for calculating skin friction and heat transfer results. Knowledge of velocity or temperature profiles is not necessary in using this calculation method. When the method is applied to a convergent-divergent, axially symmetric rocket nozzle, it shows that high rates of heat transfer are obtained at the initial stagnation point and at the throat of the nozzle. Also indicated are negative displacement thicknesses in the convergent portion of the nozzle; these occur because of the high density within the lower portions of the cooled boundary layer. (author)

Cohen, Clarence B↗