Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

The numerical evaluation of maximum-likelihood estimates of the parameters for a mixture of normal distributions from partially identified samples

Likelihood equations determined by the two types of samples which are necessary conditions for a maximum-likelihood estimate are considered. These equations, suggest certain successive-approximations iterative procedures for obtaining maximum-likelihood estimates. These are generalized steepest ascent (deflected gradient) procedures. It is shown that, with probability 1 as N sub 0 approaches infinity (regardless of the relative sizes of N sub 0 and N sub 1, i=1,...,m), these procedures converge locally to the strongly consistent maximum-likelihood estimates whenever the step size is between 0 and 2. Furthermore, the value of the step size which yields optimal local convergence rates is bounded from below by a number which always lies between 1 and 2.

Walker, H. F.↗
NASA NTRS1987

Inelastic strain analogy for piecewise linear computation of creep residues in built-up structures

An analogy between inelastic strains caused by temperature and those caused by creep is presented in terms of isotropic elasticity. It is shown how the theoretical aspects can be blended with existing finite-element computer programs to exact a piecewise linear solution. The creep effect is determined by using the thermal stress computational approach, if appropriate alterations are made to the thermal expansion of the individual elements. The overall transient solution is achieved by consecutive piecewise linear iterations. The total residue caused by creep is obtained by accumulating creep residues for each iteration and then resubmitting the total residues for each element as an equivalent input. A typical creep law is tested for incremental time convergence. The results indicate that the approach is practical, with a valid indication of the extent of creep after approximately 20 hr of incremental time. The general analogy between body forces and inelastic strain gradients is discussed with respect to how an inelastic problem can be worked as an elastic problem.

Jenkins, Jerald M.↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

The numerical evaluation of maximum-likelihood estimates of the parameters for a mixture of normal distributions from partially identified samples

Likelihood equations determined by the two types of samples which are necessary conditions for a maximum-likelihood estimate were considered. These equations suggest certain successive approximations iterative procedures for obtaining maximum likelihood estimates. The procedures, which are generalized steepest ascent (deflected gradient) procedures, contain those of Hosmer as a special case.

Walker, H. F.↗
NASA NTRS1974

First results of a general circulation model applied to the SST-NOx problem

Results of model runs, one using two-dimensional distribution of NO2 in the unperturbed stratosphere and another including an additional localized source of NO2 to approximate the effect of SSTs, are reported. The general circulation model and chemical reactions are described, and corrections in the previous latitudinal and seasonal gradients of total columnar ozone and the diffusion coefficient in the neighborhood of the tropopause are noted. Excellent agreement with previous observations was obtained. Global distributions of ozone and NO2 are described and represented in graphs. Results indicate ozone reduction of approximately 16% in the Northern Hemisphere and approximately 8% in the Southern Hemisphere, with the mid-latitude source of NO2 apparently having a blocking effect on the horizontal transport of ozone, resulting in larger reductions of ozone at high latitudes than at low ones.

Cunnold, D. M.↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS2012

From the Blazar Sequence to the Blazar Envelope: Revisiting the Relativistic Jet Dichotomy in Radio-Loud AGN

We revisit the concept of a blazar sequence that relates the synchrotron peak frequency (Vpeak) in blazars with synchrotron peak luminosity (Lpeak, in vLv) using a large sample of radio-loud AGN. We present observational evidence that the blazar sequence is formed from two populations in the synchrotron Vpeak - Lpeak plane, each forming an upper edge to an envelope of progressively misaligned blazars, and connecting to an adjacent group of radio galaxies having jets viewed at much larger angles to the line of sight. When binned by jet kinetic power (Lkin; as measured through a scaling relationship with extended radio power), we find that radio core dominance decreases with decreasing synchrotron Lpeak, revealing that sources in the envelope are generally more misaligned. We find population-based evidence of velocity gradients in jets at low kinetic powers (approximately 10(exp 42) - 10(exp 44.5) erg s(exp -1)), corresponding to FR I radio galaxies and most BL Lacs. These low jet power 'weak jet' sources, thought to exhibit radiatively inefficient accretion, are distinguished from the population of non-decelerating, low synchrotron-peaking (LSP) blazars and FR II radio galaxies ('strong' jets) which are thought to exhibit radiatively efficient accretion. The two-population interpretation explains the apparent contradiction of the existence of highly core-dominated, low-power blazars at both low and high synchrotron peak frequencies, and further implies that most intermediate synchrotron peak (ISP) sources are not intermediate in intrinsic jet power between LSP and high synchrotron-peaking (HSP) sources, but are more misaligned versions of HSP sources with similar jet powers.

Meyer, Eileen T.↗
NASA NTRS1980

A simple approximate result for the maximum growth rate of baroclinic instabilities

The Charney problem for baroclinic instability involves the quasi-geostrophic instability of a zonal flow on a beta plane where the zonal flow is characterized by a constant vertical shear. The atmosphere is non-Boussinesq and continuous. The solution of this problem involves confluent hypergeometric functions, and the mathematical difficulty of the problem has precluded extracting simple results of generality. It is shown that there exists very simple, powerful approximate result for the growth rate of the most rapidly growing instability, viz., that this growth rate is linearly proportional to the surface meridional temperature gradient. The coefficient of proportionality is also easily determined. The result extends to substantially more general profiles than those in the Charney problem.

Lindzen, R. S.↗
NASA NTRS1973

Long-term motion in a restricted problem of rotational motion

A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the long-term effects of gravity-gradient torque and orbital evolution on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is constrained to precess and spin with constant rates. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite to second order in n/omega sub 0, where n is the orbital mean motion of the center of mass and omega sub 0 is a reference value of the magnitude of the satellite's rotational angular velocity. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.

Cochran, J. E.↗
DOE OSTI2024

Identifying Controlling Variables for Mercury Vapors in Alpha-4 at Y-12: Two Year Data Collection Update

Multiple sensor packages were deployed at Alpha-4 by SRNL, in collaboration with United Cleanup Oak Ridge LLC (UCOR), to monitor mercury vapor concentrations and meteorological parameters. These sensors collected data, inside and outside of the legacy-use facility, for approximately two years. Though data gaps still exist, particularly in colder months, several controlling variables were identified that govern mercury vapor concentrations within Alpha-4. Temperature, barometric pressure gradients, humidity, and wind speed have been identified as controlling variables. A strong positive correlation was seen between mercury vapor concentrations and temperature which generally followed diurnal fluctuations. Temperatures below approximately 10 degrees Celsius did not show any spikes above the PEL, indicating more work can be performed at any time during the winter months – more data should be collected to confirm consistency in this finding. Additionally, spikes in mercury vapor concentrations above that of the permissible exposure limit (PEL; 100 µg/m 3 ) occurred primarily in late afternoon or evening/overnight hours (between 3 PM and 6 AM), which suggests D&D operations might be best scheduled during morning or daytime hours prior to the late afternoon. However, a limited number of spikes did occur outside of the identified window, although this may be attributed to disturbances in air flow and mercury vapor release from work activities performed inside of the Alpha-4 building. The analysis conducted allows for a strong predictive capability for estimating mercury vapor concentrations based upon accurate meteorological parameters. Still, additional data collection, particularly in the winter months, could help to strengthen the predictive power and validate the identified data trends. Further, increased temporal resolution could also help to better characterize the incipient stages of the increases and decreases in the mercury vapor concentration. Continued monitoring support by SRNL at Y-12 is underway at Alpha-4 to further close remaining data gaps and support deactivation and decommissioning work. Within a collaborative effort with UCOR, the SRNL team is collecting mercury vapor data to study the efficacy of a novel mercury suppressant, FerroBlack® which was recently deployed at Alpha-4. In addition, modifications to the current monitoring setup to increase measurement resolution is also being investigated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗
NASA NTRS1995

The anomalous nuclear component in the three-dimensional heliosphere

The anomalous nuclear component is neither of solar nor galactic cosmic ray origin. Its presence in the heliosphere is an independent probe for both interplanetary electrodynamical investigations--especially solar modulation--and probably the most direct means for determining the elemental and isotopic composition of those neutral atoms in the local interstellar medium that have high first ionization potentials (e.g., He, N, O, Ne, Ar, etc.). This report is a brief account of the evolution of our knowledge of this component for readers not specializing in the subject. Included are the initial discoveries of the component, its composition, spectra, heliospheric radial and latitudinal intensity gradients, modulation over the approximately 22 year solar magnetic cycle, trapping in the magnetosphere and its use for estimating the location of a heliospheric termination shock. Recent measurements from the ULYSSES spacecraft have provided conclusive evidence that incoming neutral atoms, after photoionization by solar uv, are picked up by the solar wind, thus lending support for the concept that after their acceleration--probably at a termination shock--they return to the inner heliosphere as pseudo-cosmic rays. ULYSSES spacecraft investigations extending to approximately 56 deg south latitude reveal, for both the anomalous nuclear component and the galactic cosmic rays, that there is a surprisingly small latitudinal intensity gradient. Thus, for the current phase of the solar cycle, modulation is much more spherically symmetric in the inner solar system than had generally been believed. A further surprise is the continual presence of approximately 26 day recurrent modulation at hight latitudes, without corresponding magnetic field compressions. These results are changing our ideas and requiring modification of our models for solar modulation in three dimensions.

Simpson, J. A.↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
NASA NTRS1985

The response of turbulent boundary layers to sudden pertubations

The processes present in a turbulent boundary layer (BL) which experiences a sudden change in the surface roughness are discussed in terms of numerical models in comparison with empirical data. Consideration is limited to flow perturbations which preserve the BL approximation. Generalized continuity, momentum and turbulent energy transport equations are reviewed. The perturbed BL has two regions: one with large velocity gradients and an outer zone with smaller gradients. Expressions are defined for the velocity and temperature profiles and the roughness length after a step change in surface roughness, which causes the disturbance. A pertrurbation could also arise from a heat flux input of a pressure gradient. The effects of the shape of the perturbing step and the presence of several disturbing edges are explored. Finally, it is shown that closure forms and constants for self-preserving flows are worthy of further development because of their demonstrable abilities to describe simple perturbed flows.

Smits, A. J.↗
NASA NTRS2008

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↗
DOE OSTI2024

Symmetry-Based Structured Matrices for Efficient Approximately Equivariant Networks

There has been much recent interest in designing symmetry-aware neural networks (NNs) exhibiting relaxed equivariance. Such NNs aim to interpolate between being exactly equivariant and being fully flexible, affording consistent performance benefits. In a separate line of work, certain structured parameter matrices -- those with displacement structure, characterized by low displacement rank (LDR) -- have been used to design small-footprint NNs. Displacement structure enables fast function and gradient evaluation, but permits accurate approximations via compression primarily to classical convolutional neural networks (CNNs). In this work, we propose a general framework -- based on a novel construction of symmetry-based structured matrices -- to build approximately equivariant NNs with significantly reduced parameter counts. Our framework integrates the two aforementioned lines of work via the use of so-called Group Matrices (GMs), a forgotten precursor to the modern notion of regular representations of finite groups. GMs allow the design of structured matrices -- resembling LDR matrices -- which generalize the linear operations of a classical CNN from cyclic groups to general finite groups and their homogeneous spaces. We show that GMs can be employed to extend all the elementary operations of CNNs to general discrete groups. Further, the theory of structured matrices based on GMs provides a generalization of LDR theory focussed on matrices with cyclic structure, providing a tool for implementing approximate equivariance for discrete groups. We test GM-based architectures on a variety of tasks in the presence of relaxed symmetry. We report that our framework consistently performs competitively compared to approximately equivariant NNs, and other structured matrix-based compression frameworks, sometimes with a one or two orders of magnitude lower parameter count.

Samudre, Ashwin↗
NASA NTRS1988

Analytical approximation of a distorted reflector surface defined by a discrete set of points

Reflector antennas on Earth orbiting spacecrafts generally cannot be described analytically. The reflector surface is subjected to a large temperature fluctuation and gradients, and is thus warped from its true geometrical shape. Aside from distortion by thermal stresses, reflector surfaces are often purposely shaped to minimize phase aberrations and scanning losses. To analyze distorted reflector antennas defined by discrete surface points, a numerical technique must be applied to compute an interpolatory surface passing through a grid of discrete points. In this paper, the distorted reflector surface points are approximated by two analytical components: an undistorted surface component and a surface error component. The undistorted surface component is a best fit paraboloid polynomial for the given set of points and the surface error component is a Fourier series expansion of the deviation of the actual surface points, from the best fit paraboloid. By applying the numerical technique to approximate the surface normals of the distorted reflector surface, the induced surface current can be obtained using physical optics technique. These surface currents are integrated to find the far field radiation pattern.

Acosta, Roberto J.↗
NASA NTRS1979

Convective instability when the temperature gradient and rotation vector are oblique to gravity. I - Fluids without diffusion

A linear stability analysis of fluid layers under uniform rotation (generally oblique to gravity) which possess both vertical and horizontal temperature gradients is made by considering ideal fluids without diffusion within a Boussinesq approximation. This simplified configuration is used to assess the preferred convective modes as a function of latitude on a planet like Jupiter. The tilted rotation vector introduces a preference for roll-like disturbances with north-south orientations, while the horizontal temperature gradient produces a thermal wind shear which favors convective rolls oriented parallel to the flow in an east-west direction. It is found that the horizontal temperature gradient needed to produce a preference for the axisymmetric or east-west rolls increases with an increasing rotation rate and a decreasing latitude. The parameter values for Jupiter are estimated with the use of a simple radiative convective model, indicating a preference for axisymmetric rolls at nearly all latitudes if the convection zone depth is greater than about 200 km below the one atmosphere pressure level and convective roll characteristics which contribute to an equatorial acceleration.

Hathaway, D. H.↗