Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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 577 records · Page 32

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three dimension

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Anderson, B. H.↗

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three-dimension

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Anderson, B. H.↗

Space radiation interaction mechanisms in materials

Models of charged-particle impact under conditions typical of the space environment are reported, with a focus on impact excitation and nuclear reactions, especially for heavy ions. Impact excitation is studied by using a global model for electronic excitation based on formal relations through the classical dielectric function to derive an approximation related to the local plasma (electron density distribution) within the atoms and molecules; and corrections to the model resulting from the nonfluid nature of this plasma are discussed. Nuclear reactions are studied by reducing quantum-mechanical treatment of this general N-body problem to an equivalent two-body problem that is solvable, and by comparing the results with experimental data. The equations for heavy-charged-particle transport are derived and solution techniques demonstrated. Finally, these methods of analysis are applied to study the change in the electrical properties of a GaAs semiconductor for photovoltaic applications; proton damage to GaAs crystals is found to arise from stable replacement defects and to be nonannealable, in contrast to electron-induced damage.

Wilson, J. W.↗

Parameter estimation for the Euler-Bernoulli-beam

An approximation involving cubic spline functions for parameter estimation problems in the Euler-Bernoulli-beam equation (phrased as an optimization problem with respect to the parameters) is described and convergence is proved. The resulting algorithm was implemented and several of the test examples are documented. It is observed that the use of penalty terms in the cost functional can improve the rate of convergence.

Graif, E.↗

Computing approximate random Delta v magnitude probability densities

This paper describes the development and use of an algorithm to compute approximate statistics of the magnitude of a single random trajectory correction maneuver (TCM) Delta v vector. The TCM Delta v vector is modeled as a three component Cartesian vector each of whose components is a random variable having a normal (Gaussian) distribution with zero mean and possibly unequal standard deviations. The algorithm uses these standard deviations as input to produce approximations to (1) the mean and standard deviation of the magnitude of Delta v, (2) points of the probability density function of the magnitude of Delta v, and (3) points of the cumulative and inverse cumulative distribution functions of Delta v. The approximates are based on Monte Carlo techniques developed in a previous paper by the author and extended here. The algorithm described is expected to be useful in both pre-flight planning and in-flight analysis of maneuver propellant requirements for space missions.

Chadwick, C.↗

On the relative scattering of P- and S-waves

Using a single-scattering approximation, equations for the scattering attenuation coefficients of P-body and S-body waves are derived. The results are discussed in the light of the energy-renormalization approaches of Wu (1980, 1982) and Sato (1982) to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P-waves to S-waves and S-waves to P-waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average-wavenumber power spectra. The power spectra are approximated with piecewise constant functions, each segment of which contributes to the net frequency-dependent scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behavior associated with each spectral segment, including a transition peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observe such a peak in the apparent attenuation of seismic waves. Both the frequency and distance limits on the accuracy of the theoretical results are given.

Malin, P. E.↗

Surface photometry of comet P/Encke

A scheme to clean cometary digital images from offending background-star trails, and this technique has been applied to a pair of deep Kitt Peak 4-m plates of comet P/Encke, taken in October 1980. Simultaneous and subsequent digital spectra have been obtained at Lick Observatory. The non-polluted coma images show a strong asymmetric sunward-oriented fan/jet, and an extended and rounder (mostly gaseous) main coma, out to approximately 100,000 km radius. The stellar-trail point-spread function has a narrow width (sigma approximately 0.6 arcsec), so that spatial resolution better than approximately 300 km is achieved at the comet. The photometric gradient near the nucleus is very steep, strongly suggesting an icy-grain component which evaporates quickly (at radii equal to or less than 500 km) in the sunlight. Further from the nucleus, the profile becomes shallower, bluer, and more gas dominated. The effect of solar radiation pressure on C2, CN, and other molecules is probably responsible for the rounding of the outer, fainter isophotes. The source of the molecules is likely to be larger than the nucleus itself, and a substantial fraction may originate in the jet. The technique described here may also be applicable in surface photometry of galaxies, in cases where the heavy image pollution by foreground stars is present.

Djorgovski, S.↗

Theoretical models of single-element LTSA antenna elements

The spectral domain computations for slotline guide wavelength and characteristic impedance provide the basic parameters that are need to calculate the radiation pattern of an Linearly Tapered Slot Antenna (LTSA) by using the half-plane Green's function and a stepped approximation to the LTSA structure. However, the spectral domain computations are somewhat time consuming and not suitable for computer aided design application. Therefore, the results of several computation were fitted by formulas that are readily evaluated on a calculator or computer. Antenna measurements were continued, and an attempt was made to predict the truncation effect by three different methods.

Source record↗

Effects of multiple scattering on time- and depth-resolved signals in airborne lidar systems

A semianalytic Monte Carlo radiative transfer model (SALMON) is employed to probe the effects of multiple-scattering events on the time- and depth-resolved lidar signals from homogeneous aqueous media. The effective total attenuation coefficients in the single-scattering approximation are determined as functions of dimensionless parameters characterizing the lidar system and the medium. Results show that single-scattering events dominate when these parameters are close to their lower bounds and that when their values exceed unity multiple-scattering events dominate.

Punjabi, A.↗

A continued program of planetary study at the University of Texas McDonald Observatory

The beginning of eclipses of the Pluto-Charon system was detected. The onset of coma formation of P/Halley at 5.4 au was detected and evidence of sublimation at 4.8 au when CN emission was detected. Extensive spatial maps of the gas in the comae of comets Halley and Giacobini-Zinner were obtained in fall 1985. Halley was time variable, and Giacobini-Zinner was depleted in C2 and C3 relative to CN. Comet Kopff was shown to have a pre-perihelion brightness maximum of its gas, consistent with mantle development if the comet is a high obliquity object. New Haser model scale lengths for CN, C3, and C2 were determined using results from the Faint Comet Survey. Spectra of 12 asteroids in unusual orbits showed no evidence of any comet-like emission features. In particular, 3200 Phaethon (1983 TB) has no gas or dust coma, in spite of the similarity of its orbit with the Geminid meteor stream. Data were analyzed on Saturn's H2 and CH4 bands for the recent southern summer using a Tomasko-Doose type of haze distribution. This haze model fits the data moderately well, giving a CH4 mixing ratio of (4.2 + or 0.4)x003. Simple functions were found to approximate the collision-induced rotation-translation thermal opacity of H2.

Smith, H. J.↗

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗

On the numerical performance of three-dimensional thick shell elements using a hybrid/mixed formulation

Three-dimensional thick shell elements with 8, 16, and 18 nodes are formulated by using the hybrid/mixed method. In bending applications, these elements are free from locking effect and give improved stress predictions. Finite element equations are derived from the Hellinger-Reissner variational principle in which both the displacement and stress fields are approximated by independent interpolation functions. For the assumption of stress parameters, three guidelines are followed: (1) suppression of kinematic deformation modes, (2) invariant element property, and (3) the constraint index exhibited by the element, when applied to constrained-media problems, must be greater than or equal to one. Numerical results are presented to show the element's behavior characteristics regarding sensitivity to locking, distortion effect (patch tests), mesh convergence and the accuracy of stress evaluation.

Graf, W.↗

Statistical-mechanical selection of the shapes of disk galaxies

A new method is proposed for selecting steady state shapes of disk galaxies as 'most probable states' of a large number of stars, given only the total energy and total angular momentum. A partial differential equation is derived for the mean gravitational potential; it is closely related to the 'sinh-Poisson' equation for the mean-field description of a line vortex system or electrostatic guiding-center plasma. A 'water bag' approximation to the distribution function for bound stars renders the equation analytically tractable, but accurate solution of it may require numerical integration, in view of its general nonlinearity.

Montgomery, David↗

Solutions to bi-Maxwellian transport equations for radial solar wind beyond 28 R(S)

This paper presents solar wind solutions for radial flow between 28 solar radii and 1 AU using the bi-Maxwellian-based 16-moment set of transport equations. In addition to the number density, drift velocity, and parallel and perpendicular temperatures, the 16-moment equations account for the transport of both longitudinal and transverse thermal energies as well as stress. Also, using the 16-moment approximation for the distribution function and assuming plasma parameter values characteristic of the solar wind, contour plots are generated for the proton velocity distribution function. It is shown how the shape of these plots depends on various macroscopic plasma parameters.

Demars, H. G.↗

Solar particle abundances at energies of greater than 1 MeV per nucleon and the role of interplanetary shocks

The abundances of elements in large solar energetic-particle events in the energy range of 2-12 MeV per nucleon are examined. It is confirmed that the abundances relative to mean values vary approximately monotonically as a function of mass, except for He-4; some events show a gradual depletion of heavy ions, whereas a small number displays a gradual increase. A further organization of abundance data is shown, which depends on the longitude of the source region. Enhancements in Fe/C and other heavy elements relative to C occur when source regions are near west 60 deg; the enhancements are attributed to the sampling of a flare-heated material. Depletions of these elements are found to be greatest for source regions near central meridian; they are matched by a steepening of the spectrum and can be understood in terms of diffusive shock acceleration.

Cane, H. V.↗

Compensator design for stability enhancement with collocated controllers

A continuum model rather than a finite element model is used. The optimal compensator design is formulated as a stochastic regulator problem and is shown to be solvable by the general infinite-dimensional theory developed by the author despite the lack of exponential stabilizability. Infinite-dimensional steady-state Riccati equations characterizing the feedback control gain and the Kalman filter gain operators can be solved explicitly. The associated performance indexes including the mean-square control-effort are calculated in closed form. As a first approximation, the compensator transfer function can be realized as a bank of bandpass filters in parallel centered at the undamped mode frequencies. Numerical calculations for the gain and bandwidths for a typical configuration are presented. The performance of the compensator is evaluated when in fact in the true model there is no actuator noise. The theoretical problem involved is to show that the infinite-dimensional stochastic process is asymptotically stationary. It is possible to calculate the steady-state covariance in closed form and thereby calculate performance indexes of interest explicitly, facilitating the choice of optimal design parameters.

Balakrishnan, A. V.↗

Heat flow vs. atmospheric greenhouse on early Mars

Researchers derived a quantitative relationship between the effectiveness of an atmospheric greenhouse and internal heat flow in producing the morphological differences between earlier and later Martian terrains. The derivation is based on relationships previously derived by other researchers. The reasoning may be stated as follows: the CO2 mean residence time in the Martian atmosphere is almost certainly much shorter than the total time span over which early climate differences are thought to have been sustained. Therefore, recycling of previously degassed CO2 quickly becomes more important than the ongoing supply of juvenile CO2. If so, then the atmospheric CO2 pressure, and thereby the surface temperature, may be approximated mathematically as a function of the total degassed CO2 in the atmosphere plus buried material and the ratio of the atmospheric and regolith mean residence times. The latter ratio can also be expressed as a function of heat flow. Hence, it follows that the surface temperature may be expressed as a function of heat flow and the total amount of available CO2. However, the depth to the water table can simultaneously be expressed as a function of heat flow and the surface temperature (the boundary condition). Therefore, for any given values of total available CO2 and regolith conductivity, there exist coupled independent equations which relate heat flow, surface temperature, and the depth to the water table. This means we can now derive simultaneous values of surface temperature and the depth of the water table for any value of the heat flow. The derived relationship is used to evaluate the relative importance of the atmospheric greenhouse effect and the internal regolith thermal gradient in producing morphological changes for any value of the heat flow, and to assess the absolute importance of each of the values of the heat flow which are thought to be reasonable on independent geophysical grounds.

Fanale, F. P.↗

Scattering behavior of Lunar Lake playa determined from parabola bidirectional reflectance data

Bidirectional reflectance data obtained with the Portable Apparatus for Rapid Acquisition of Bidirectional Observations of Land and Atmosphere instrument at Lunar Lake, Nevada, were analyzed to determine the scattering properties of playas. The data are approximated by a Lambertian function, except at high phase angles in the solar principal plane, where Fresnel reflectance appears to dominate. The data also depart slightly (15 percent) from a Lambertian function at high emission angles in all azimuthal planes. No published photometric model accounts for the observed data. It is hypothesized that the observed scattering behavior is the superposition of volume and surface (Fresnel) scattering mechanisms.

Shepard, Michael K.↗