Search NASA⌕ Search

SEARCH · Search NASA

Results for “integrated gradients”

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

Efficient On-Orbit Singularity-Free Geopotential Estimation

For vehicles in orbital cruise, either in orbit around a body or in a transfer orbit, gravity is the primary external force; therefore, the complexity of the geopotential model heavily influences the accuracy of the navigation state. As this type of spacecraft requires a highly accurate geopotential model and performs the computation numerous times a second, this calculation needs to be efficient. This calculation needs to include the first and second derivatives of the geopotential because navigation systems can require both the gravitational acceleration and gradient for inertial integration and state filtering. The most efficient method for calculating the geopotential, the forward column method, is discussed in detail, as well as a method to avoid the singularities that exist when using this method. This is shown to decrease the computation time of the geopotential compared to other popular methods. In addition, methods for first and second order propagation are discussed, which decrease the rate at which the full geopotential model needs to be calculated while maintaining the accuracy required of navigation systems. These estimation methods are shown to decrease the necessary computation of the gravity model by multiple orders of magnitude. This method decreases the computation time of the Exploration Upper Stage navigation system geopotential model by potentially an order of magnitude compared with the previous model without affecting the navigation error. A properly implemented geopotential model can have the accuracy of an 8 by 8 model while approaching the computational requirement as using only the J234 coefficients.

Amert, Joel↗

A finite element: Boundary integral method for electromagnetic scattering

A method that combines the finite element and boundary integral techniques for the numerical solution of electromagnetic scattering problems is presented. The finite element method is well known for requiring a low order storage and for its capability to model inhomogeneous structures. Of particular emphasis in this work is the reduction of the storage requirement by terminating the finite element mesh on a boundary in a fashion which renders the boundary integrals in convolutional form. The fast Fourier transform is then used to evaluate these integrals in a conjugate gradient solver, without a need to generate the actual matrix. This method has a marked advantage over traditional integral equation approaches with respect to the storage requirement of highly inhomogeneous structures. Rectangular, circular, and ogival mesh termination boundaries are examined for two-dimensional scattering. In the case of axially symmetric structures, the boundary integral matrix storage is reduced by exploiting matrix symmetries and solving the resulting system via the conjugate gradient method. In each case several results are presented for various scatterers aimed at validating the method and providing an assessment of its capabilities. Important in methods incorporating boundary integral equations is the issue of internal resonance. A method is implemented for their removal, and is shown to be effective in the two-dimensional and three-dimensional applications.

Collins, J. D.↗

A displacement gradient BEM for accurate stress computation near boundaries in 2-D anisotropic problems

A displacement gradient method for 2D anisotropic elasticity problems is presented, which effectively minimizes the boundary layer effect through a two-step procedure. First, the boundary integral equations are solved for the unknown boundary displacements and tractions. Second, a direct integral equation for displacement gradients is developed in terms of boundary tractions. Three methods based on different evaluation procedures and locations for determining the displacement gradients are proposed. In the first method the displacement gradients are averaged at nodes common to adjacent elements. The second method stores the gradients element-wise. In the third method, the gradients are evaluated at the nodes of discontinuous elements. The three methods are applied to near-isotropic plates with circular and elliptic cutouts. It is concluded that all three methods can yield accurate stress distributions.

Sistla, R.↗

Two-Dimensional High-Lift Aerodynamic Optimization Using Neural Networks

The high-lift performance of a multi-element airfoil was optimized by using neural-net predictions that were trained using a computational data set. The numerical data was generated using a two-dimensional, incompressible, Navier-Stokes algorithm with the Spalart-Allmaras turbulence model. Because it is difficult to predict maximum lift for high-lift systems, an empirically-based maximum lift criteria was used in this study to determine both the maximum lift and the angle at which it occurs. The 'pressure difference rule,' which states that the maximum lift condition corresponds to a certain pressure difference between the peak suction pressure and the pressure at the trailing edge of the element, was applied and verified with experimental observations for this configuration. Multiple input, single output networks were trained using the NASA Ames variation of the Levenberg-Marquardt algorithm for each of the aerodynamic coefficients (lift, drag and moment). The artificial neural networks were integrated with a gradient-based optimizer. Using independent numerical simulations and experimental data for this high-lift configuration, it was shown that this design process successfully optimized flap deflection, gap, overlap, and angle of attack to maximize lift. Once the neural nets were trained and integrated with the optimizer, minimal additional computer resources were required to perform optimization runs with different initial conditions and parameters. Applying the neural networks within the high-lift rigging optimization process reduced the amount of computational time and resources by 44% compared with traditional gradient-based optimization procedures for multiple optimization runs.

Greenman, Roxana M.↗

High-Lift Optimization Design Using Neural Networks on a Multi-Element Airfoil

The high-lift performance of a multi-element airfoil was optimized by using neural-net predictions that were trained using a computational data set. The numerical data was generated using a two-dimensional, incompressible, Navier-Stokes algorithm with the Spalart-Allmaras turbulence model. Because it is difficult to predict maximum lift for high-lift systems, an empirically-based maximum lift criteria was used in this study to determine both the maximum lift and the angle at which it occurs. Multiple input, single output networks were trained using the NASA Ames variation of the Levenberg-Marquardt algorithm for each of the aerodynamic coefficients (lift, drag, and moment). The artificial neural networks were integrated with a gradient-based optimizer. Using independent numerical simulations and experimental data for this high-lift configuration, it was shown that this design process successfully optimized flap deflection, gap, overlap, and angle of attack to maximize lift. Once the neural networks were trained and integrated with the optimizer, minimal additional computer resources were required to perform optimization runs with different initial conditions and parameters. Applying the neural networks within the high-lift rigging optimization process reduced the amount of computational time and resources by 83% compared with traditional gradient-based optimization procedures for multiple optimization runs.

Greenman, Roxana M.↗

A similarity rule for compressibility and sidewall-boundary-layer effects in two-dimensional wind tunnels

The effect of the sidewall boundary layer on flow in two-dimensional wind tunnels is determined. The small-disturbance and isentropic approximations are made, and the sidewall-boundary-layer dynamics are modeled with the von Karman momentum-integral equation. The effects of the edge-velocity-gradient term in the sidewall momentum integral, which is usually dominant near the model, and the compressibility term are shown to be similar. It is shown that the effect of sidewall suction around the model is not similar to two-dimensional flow. Comparisons with experiment are made to verify the similarity rule.

Barnwell, R. W.↗

Electrophoretic mobility of cells in a vertical Ficoll gradiant

The upward migration of living cells and test particles under the influence of a constant electric field in a low conductivity Ficoll gradient occurs at nearly constant velocity. Viscosity and neutral polymer concentration affect migration rate. Decreasing viscosity speeds up the particle migration, decreasing neutral polymer (Ficoll) concentration, slows particle migration, since electrophoretic mobility increases approximately linearly with neutral polymer concentration. Neutral polymers interact with the cell surface to effectively raise its zeta potential. An analytic function was developed from the known dependence of these physical variables on migration distance; the analysis expresses migration velocity as an explicit function of position in the density gradient. It predicts an almost linear increase in velocity of about 12 to 16% over the working region of the gradient. It was numerically integrated and correctly predicts cell migration distance vs time curves without the use of any fitted parameters. The resulting migration curves follow the expected slowly varying exponential form that closely resembles a straight line. The ability to determine standard electrophoretic mobilities from such curves depends on knowledge of the effect of Ficoll on the zeta potential of the cell type that is separated.

Plank, L. D.↗

Convergent-tapered annular seals - Analysis and testing for rotordynamic coefficients

A combined analytical-computational method has been developed to calculate the pressure field and dynamic coefficients for tapered high-pressure annular seals. Completely developed turbulent flow is assumed in both the circumferential and axial directions, according to Hirs' (1973) bulk-flow turbulent-lubrication equations. In a numerical experiment with the method, a short bearing approximation is used to derive an analytical expression for the first-order (dynamic) pressure gradient. This expression is integrated numerically to define dynamic coefficients of the seal. Numerical results of the integration are compared with previous results for straight and tapered seals. It is found that direct stiffness and leakage coefficients increase in the present seal, while the remaining coefficients decrease. An optimal taper angle is shown to exist with respect to: (1) the direct stiffness and (2) the ratio of direct stiffness to leakage. On the basis of the optimal taper angle calculations, stiffness increases on the order of 40 to 50 percent are predicted. Experimental results for three taper angles are presented, and are found to be in good agreement with the theoretical predictions.

Childs, D. W.↗

'Downward control' of the mean meridional circulation and temperature distribution of the polar winter stratosphere

According to the 'downward control' principle, the extratropical mean vertical velocity on a given pressure level is approximately proportional to the meridional gradient of the vertically integrated zonal force per unit mass exerted by waves above that level. In this paper, a simple numerical model that includes parameterizations of both planetary and gravity wave breaking is used to explore the influence of gravity wave breaking in the mesosphere on the mean meridional circulation and temperature distribution at lower levels in the polar winter stratosphere. The results of these calculations suggest that gravity wave drag in the mesosphere can affect the state of the polar winter stratosphere down to altitudes below 30 km. The effect is most important when planetary wave driving is relatively weak: that is, during southern winter and in early northern winter. In southern winter, downwelling weakens by a factor of 2 near the stratospause and by 20% at 30 km when gravity wave drag is not included in the calculations. As a consequence, temperatures decrease considerably throughout the polar winter stratosphere (over 20 K above 40 km and as much as 8 K at 30 km, where the effect is enhanced by the long radiative relaxation timescale). The polar winter states obtained when gravity wave drag is omitted in this simple model resemble the results of simulations with some general circulation models and suggest that some of the shortcomings of the latter may be due to a deficit in mesospheric momentum deposition by small-scale gravity waves.

Garcia, Rolando R.↗

The Epoch of Disk Settling: Z Approximately Equal to 1 to Now

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies over 0.2 < z < 1.2. DEEP2 provides a large sample of high resolution galaxy spectra and dual-band Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties, in order to examine trends. At a fixed stellar mass, galaxies systematically decrease in disturbed motions and increase in rotation velocity and potential well depth with time. The most massive galaxies are the most well-ordered at all times, with higher rotation velocities and less disturbed motions compared to less massive galaxies. We quantify disturbed motions with an integrated gas velocity dispersion (sigma(sub g)), which is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Due to finite slit width and seeing, sigma(sub g) integrates over unresolved velocity gradients which can correspond to non-ordered gas kinematics such as small-scale velocity gradients, gas motions due to star-formation, or super-imposed clumps along the line-of-sight. We compile surveys of galaxy kinematics over 1.2 < z < 3.8 and do not find any trends with redshift, likely because these studies are biased toward the most highly star-forming systems. In summary, over the last approx 8 billion years since z = 1.2, blue galaxies evolve from disturbed to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies always being the most evolved at any time.

Kassin, Susan A.↗

Analytic Formulation and Numerical Implementation of an Acoustic Pressure Gradient Prediction

The scattering of rotor noise is an area that has received little attention over the years, yet the limited work that has been done has shown that both the directivity and intensity of the acoustic field may be significantly modified by the presence of scattering bodies. One of the inputs needed to compute the scattered acoustic field is the acoustic pressure gradient on a scattering surface. Two new analytical formulations of the acoustic pressure gradient have been developed and implemented in the PSU-WOPWOP rotor noise prediction code. These formulations are presented in this paper. The first formulation is derived by taking the gradient of Farassat's retarded-time Formulation 1A. Although this formulation is relatively simple, it requires numerical time differentiation of the acoustic integrals. In the second formulation, the time differentiation is taken inside the integrals analytically. The acoustic pressure gradient predicted by these new formulations is validated through comparison with the acoustic pressure gradient determined by a purely numerical approach for two model rotors. The agreement between analytic formulations and numerical method is excellent for both stationary and moving observers case.

Lee, Seongkyu↗

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗

The Stokes force on a droplet in an unbounded fluid medium due to capillary effects

The Stokes force on a fluid droplet is obtained when the droplet is placed in an unbounded fluid medium and motion ensues due to an arbitrary interfacial-tension gradient on the droplet surface. The force, derived here for a spherical droplet, is proportional to the integral of the interfacial-tension gradient over the droplet surface. When the interfacial-tension gradients are caused by temperature or concentration variations, the result for the force may be further specialized when convective transport effects are negligible. In this case, it is possible to express the force in terms of the gradient of the undisturbed temperature (or concentration) field evaluated at the location of the droplet center in a form analogous to Faxen's force law.

Subramanian, R. Shankar↗

An integrated approach to the synthesis of geometrically non-linear structures

An integrated approach to the minimum weight design of geometrically nonlinear three-dimensional truss structures with geometric imperfections, subject to inequality constraints on static displacements, stresses, local buckling and cross sectional areas, is investigated. The integrated structural synthesis problem involves design and response quantities as independent variables and equilibrium equations, describing the finite element model, as equality constraints. The nonlinear structural analysis and the optimization are thus merged together into a single process. A computer program developed to compute the constraint values and analytical gradients is coupled with a generalized reduced gradient algorithm to solve the integrated problem. Numerical results for a geometrically nonlinear shallow dome example problem are presented for various types of imperfections. Furthermore, it is found that the algorithm is capable of detecting and guarding against system as well as element elastic instability using equilibrium information only, that is, without imposing system and local buckling inequality constraints.

Smaoui, H.↗

Cosmic-ray gradient measurements and modulation beyond the inner solar wind termination shock

Data provided by the IMP, Voyager, and Pioneer spacecraft for the cosmic-ray particles with E greater than 60 MeV show the existence of integral radial cosmic-ray-density gradients which were nearly constant from 1977 to 1982. In this paper, these measurements are explained using a model in which significant modulation occurs in the turbulent shocked plasma flow between the inner solar wind termination and the outer contact discontinuity separating the interstellar medium flow from the heliospheric plasma.

Quenby, J. J.↗

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Simple and efficient numerical procedures for evaluating the gradient of Newton-type potentials are presented. Convergences of both normal and tangential components of the gradient are examined. The convergence of the vector potential is also examined, and it is shown that the scheme for handling near-hypersingular integrals also is effective for the nearly singular potential terms.

Fink, Patricia W.↗

Chemical potential effects on neutrino diffusion in supernovae

The validity of imposing a zero chemical potential for neutrinos in hydrodynamic calculations of collapsing supernovae is investigated in the diffusion approximation of neutrino transport. A coupled system of equations is solved for neutrino and energy diffusion fluxes as well as lepton diffusion in a collapsing supernovae ambient medium, and the results indicate a substantial growth in the neutrino chemical potential for densities greater than 10 to the 12th power gm/cu cm. The rate of energy transport is shown to be significantly affected by increases in Fermi integrals and chemical-potential gradients accompanied by decreases in temperature, and the extent of neutrino particle/antiparticle reactions is found also to affect energy diffusion rates. It is concluded that the photon-like behavior usually assumed for neutrinos may be incorrect and that an extension of the Sn transport approximation to include lepton characteristics is necessary for a definitive answer to the question of neutrino transport in supernovae.

Mazurek, T. J.↗