Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,063 records · Page 59

Analyzing the test process using structural coverage

A large, commercially developed FORTRAN program was modified to produce structural coverage metrics. The modified program was executed on a set of functionally generated acceptance tests and a large sample of operational usage cases. The resulting structural coverage metrics are combined with fault and error data to evaluate structural coverage. It was shown that in the software environment the functionally generated tests seem to be a good approximation of operational use. The relative proportions of the exercised statement subclasses change as the structural coverage of the program increases. A method was also proposed for evaluating if two sets of input data exercise a program in a similar manner. Evidence was provided that implies that in this environment, faults revealed in a procedure are independent of the number of times the procedure is executed and that it may be reasonable to use procedure coverage in software models that use statement coverage. Finally, the evidence suggests that it may be possible to use structural coverage to aid in the management of the acceptance test processed.

Ramsey, James↗

Finite-time arbitrary-motion unsteady cascade airfoil theory for helicopter rotors in hover

A complete and detailed derivation of finite-time arbitrary-motion cascade theory is presented for both Laplace and frequency domains. This theory includes the effect of returning wakes for both single- and multibladed rotors. The generalized cascade lift-deficiency function (CLDF) is shown to be consistent with the generalized Theodorsen lift-deficiency function when the wake spacing approaches infinity or when the reduced frequency tends to infinity. This function predicts correct zero-reduced-frequency limit. Accurate and efficient numerical procedures are presented for the evaluation of the CLDF. Numerical examples comparing the CLDF with Loewy's lift deficiency function in frequency domain are presented. Accurate Pade approximants of the CLDF are constructed using a Bode-plot approach which allows for complex poles.

Dinyavari, M. A. H.↗

Soft X-ray observations of pre-main sequence stars in the chamaeleon dark cloud

Einstein IPC observations of the nearby Chamaeleon I star forming cloud show 22 well-resolved soft X-ray sources in a 1x2 deg region. Twelve are associated with H-alpha emission line pre-main sequence (PMS) stars, and four with optically selected PMS stars. Several X-ray sources have two or more PMS stars in their error circles. Optical spectra were obtained at CTIO of possible stellar counterparts of the remaining X-ray sources. They reveal 5 probable new cloud members, K7-MO stars with weak or absent emission lines. These naked X-ray selected PMS stars are similar to those found in the Taurus-Auriga cloud. The spatial distributions and H-R diagrams of the X-ray and optically selected PMS stars in the cloud are very similar. Luminosity functions indicate the Chamaeleon stars are on average approximately 5 times more X-ray luminous than Pleiad dwarfs. A significant correlation between L sub x and optical magnitude suggests this trend may continue within the PMS phase of stellar evolution. The relation of increasing X-ray luminosity with decreasing stellar ages is thus extended to stellar ages as young as 1 million years.

Feigelson, Eric D.↗

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗

Calculation of Gallium-metal-Arsenic phase diagrams

Electrical contacts and metallization to GaAs solar cells must survive at high temperatures for several minutes under specific mission scenarios. The determination of which metallizations or alloy systems that are able to withstand extreme thermal excursions with minimum degradation to solar cell performance can be predicted by properly calculated temperature constitution phase diagrams. A method for calculating a ternary diagram and its three constituent binary phase diagrams is briefly outlined and ternary phase diagrams for three Ga-As-X alloy systems are presented. Free energy functions of the liquid and solid phase are approximated by the regular solution theory. Phase diagrams calculated using this method are presented for the Ga-As-Ge and Ga-As-Ag systems.

Scofield, J. D.↗

Undamped plasma waves

This paper describes small-amplitude nonlinear plasma wave solutions to the one-dimensional Vlasov-Maxwell equations. A sufficient condition for waves of a given phase velocity to exist arbitrarily close to a given spatially uniform Vlasov equilibrium is developed, and sufficient analytical information for the construction of approximate expressions for the electric potential and distribution functions is derived, with exact knowledge of the asymptotic behavior of the error terms. These results have a very surprising physical implication: the Landau damping of small-amplitude waves is not inevitable. Instead, there exist plasma waves that trap particles even at arbitrarily small amplitude and do not damp.

Holloway, James P.↗

Helicopter main-rotor speed effects: A comparison of predicted ranges of detection from the aural detection program ICHIN and the electronic detection program ARCAS

NASA LaRC personnel have conducted a strudy of the predicted acoustic detection ranges associated with reduced helicopter main rotor speeds. This was accomplished by providing identical input information to both the aural detection program ICHIN 6, (I Can Hear It Now, version 6) and the electronic acoustic detection program ARCAS (Assessment of Rotorcraft Detection by Acoustics Sensing). In this study, it was concluded that reducing the main rotor speed of the helicopter by 27 percent reduced both the predicted aural and electronic detection ranges by approximately 50 percent. Additionally, ARCAS was observed to function better with narrowband spectral input than with one-third octave band spectral inputs and the predicted electronic range of acoustic detection is greater than the predicted aural detection range.

Mueller, Arnold W.↗

Interpreting asteroid photometry and polarimetry using a model of shadowing and coherent backscattering

The shadow-hiding models for the opposition effect and negative polarization of atmosphereless solar system bodies do not explain some experimental findings, such as the enhancing opposition effect and negative polarization with decreasing particle size down to wavelength scales. The enhancement for laboratory photometric and polarimetric data on artificial glass samples with different particle size is shown. These results are in agreement with the so-called coherent backscattering or interference mechanism proposed for the interpretation of the opposition effect and negative polarization. Two different approaches for describing the opposition effect and negative polarization produced by the shadow-interference mechanism were developed. One is based on exact electromagnetic solutions for simple scattering systems that include dipole-dipole and dipole-surface coupling; The other is based on a point-scatterer approximation characterized by model photometric and polarimetric phase functions, and the mutual shadowing effect is derived using virtual volumes associated with the point-scatterers. Both approaches yield qualitatively similar results, although neither is entirely satisfactory. We regard them as prototypes for a future unified model of shadowing and coherent backscattering. The sharp opposition effect of 44 Nysa and the asteroid albedo-polarization rule are here explained using the point-scatterer approach.

Shkuratov, Y. G.↗

Investigation of a geodesy coexperiment to the Gravity Probe B relativity gyroscope program

Geodesy is the science of measuring the gravitational field of and positions on the Earth. Estimation of the gravitational field via gravitation gradiometry, the measurement of variations in the direction and magnitude of gravitation with respect to position, is this dissertation's focus. Gravity Probe B (GP-B) is a Stanford satellite experiment in gravitational physics. GP-B will measure the precession the rotating Earth causes on the space time around it by observing the precessions of four gyroscopes in a circular, polar, drag-free orbit at 650 km altitude. The gyroscopes are nearly perfect niobium-coated spheres of quartz, operating at 1.8 K to permit observations with extremely low thermal noise. The permissible gyroscope drift rate is miniscule, so the torques on the gyros must be tiny. A drag-free control system, by canceling accelerations caused by nongravitational forces, minimizes the support forces and hence torques. The GP-B system offers two main possibilities for geodesy. One is as a drag-free satellite to be used in trajectory-based estimates of the Earth's gravity field. We described calculations involving that approach in our previous reports, including comparison of laser only, GPS only, and combined tracking and a preliminary estimate of the possibility of estimating relativistic effects on the orbit. The second possibility is gradiometry. This technique has received a more cursory examination in previous reports, so we concentrate on it here. We explore the feasibility of using the residual suspension forces centering the GP-B gyros as gradiometer signals for geodesy. The objective of this work is a statistical prediction of the formal uncertainty in an estimate of the Earth's gravitation field using data from GP-B. We perform an instrument analysis and apply two mathematical techniques to predict uncertainty. One is an analytical approach using a flat-Earth approximation to predict geopotential information quality as a function of spatial wavelength. The second estimates the covariance matrix arising in a least-squares estimate of a spherical harmonic representation of the geopotential using GP-B gradiometer data. The results show that the GP-B data set can be used to create a consistent estimate of the geopotential up to spherical harmonic degree and order 60. The formal uncertainty of all coefficients between degrees 5 and 50 is reduced by factors of up to 30 over current satellite-only estimates and up to 7 over estimates which include surface data. The primary conclusion resulting from this study is that the gravitation gradiometer geodesy coexperiment to GP-B is both feasible and attractive.

Everitt, C. W. F.↗

Formulas for the Supersonic Loading, Lift, and Drag of Flat Swept-Back Wings with Leading Edges Behind the Mach Line

The method of superposition of linearized conical flows has been applied to the calculation of the aerodynamic properties, in supersonic flight, of thin flat, swept-back wings at an angle of attack. The wings are assumed to have rectilinear plan forms, with tips parallel to the stream, and to taper in the conventional sense. The investigation covers the moderately supersonic speed range where the Mach lines from the leading-edge apex lie ahead of the wing. The trailing edge may lie ahead of or behind the Mach lines from its apex. The case in which the Mach cone from one tip intersects the other tip is not treated. Formulas are obtained for the load distribution, the total lift, and the drag due to lift. For the cases in which the trailing edge is outside the Mach cone from its apex the formulas are complete. For wings with both leading and trailing edges behind their respective Mach lines, a degree of approximation is necessary. Charts of some of the functions derived are included to facilitate computing, and several examples are worked out in outline.

Cohen, Doris↗

Dynamic response at altitude of a turbojet engine with variable area exhaust nozzle

The dynamic characteristics of turbojet engine with variable exhaust nozzle area were investigated over a range of altitudes and flight Mach numbers. These characteristics generalize to standard static level condition. Data resulting from approximate step disturbances in either manipulated variable suggested functional relationships from which functions can be derived. The constants of the transfer functions are listed for five combinations of engine speed and exhaust nozzle area. The minimum data needed to completely define the transfer functions were: experimentally determined dynamic characteristics, such as engine time constant and initial rise ratio, resulting from either manipulated variable; static characteristics determined from steady-state performance curves.

Delio, Gene J↗

A superlinear interior points algorithm for engineering design optimization

We present a quasi-Newton interior points algorithm for nonlinear constrained optimization. It is based on a general approach consisting of the iterative solution in the primal and dual spaces of the equalities in Karush-Kuhn-Tucker optimality conditions. This is done in such a way to have primal and dual feasibility at each iteration, which ensures satisfaction of those optimality conditions at the limit points. This approach is very strong and efficient, since at each iteration it only requires the solution of two linear systems with the same matrix, instead of quadratic programming subproblems. It is also particularly appropriate for engineering design optimization inasmuch at each iteration a feasible design is obtained. The present algorithm uses a quasi-Newton approximation of the second derivative of the Lagrangian function in order to have superlinear asymptotic convergence. We discuss theoretical aspects of the algorithm and its computer implementation.

Herskovits, J.↗

Stochastic solution to quantum dynamics

The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. The technique, which is developed and validated in two- and three- dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space.

John, Sarah↗

M dwarfs, microlensing, and the mass budget of the Galaxy

We show that faint red stars do not contribute significantly to the mass budget of the Galaxy or to microlensing statistics. Our results are obtained by analyzing two long exposures of a high-latitude field taken with the Wide Field Camera (WFC) on the newly repaired Hubble Space Telescope (HST). Stars are easily distinguished from galaxies essentially to the limiting magnitudes of the images. We find five stars with 2.0 less than V - I less than 3.0 and I less than 25.3 and no stars with V - I greater than 3.0. Therefore, main-sequence stars with M(sub I) greater than 10 that are above the hydrogen-burning limit in the dark halo or the spheroid contribute less than 6% of the unseen matter. Faint red disk stars, M-dwarfs, contribute at most 15% to the mass of the disk. We parameterize the faint end of the cumulative distribution of stars, Phi, as a function of luminosity L(sub V), d Phi/d ln L(sub V) proportional to L(sub V exp -gamma). For spheroid stars, gamma less than 0.32 over the range 6 less than M(sub V) less than 17, with 98% confidence. The disk luminosity function falls, gamma less than 0, for 15 approximately less than M(sub V) approximately less than 19. Faint red stars in the disk or thick disk, and stars with M(sub V) less than 16 in the spheroid contribute tau less than 10(exp -8) to the optical depth to microlensing toward the Large Magellanic Cloud.

Bahcall, John N.↗

AROTEL Temperature Retrievals Within PSC's Using Raman Scattering

NASA Goddard Space Flight Center's Airborne Raman Ozone Temperature Lidar (AROTEL) made temperature retrievals within Polar Stratospheric Clouds (PSCs) on several flights during the SAGE III Ozone Loss and Validation Experiment (SOLVE) campaign. The location of the PSCs was confirmed using simultaneously acquired data from the NASA Langley Aerosol lidar. Retrievals were made on flight dates 991207, 991210 and 000127 from just above the aircraft to 25 kilometers geometric altitude. Raman temperature retrievals are, to first order, insensitive to Mie interference because the Raman signals are red shifted by 2331 cm(exp -1) from the initial laser wavelength. Backscattering from clouds and aerosols is consequently not observed by the detector; however, extinction does impact the measurement and limits retrievals to optically thin clouds. Comparisons between retrievals employing Rayleigh and Raman scattering show the Raman temperatures to be significantly warmer than those employing Rayleigh scattering. Uncertainties are a function of altitude: at 25 km they were approximately 3 K. Temperature profiles could not be retrieved for optically thick clouds.

Burris, John↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Analytical Derivation of the Vegetation Optical Depth from the Microwave Polarization Difference Index

A numerical solution for the canopy optical depth in an existing microwave-based land surface parameter retrieval model is presented. The optical depth is derived from the microwave polarization difference index and the dielectric constant of the soil. The original procedure used an approximation in the form of a logarithmic decay function to define this relationship, and was derived through a series of lengthy polynomials. These polynomials had to be recalculated when the scattering albedo or antenna incidence angle changes. The new procedure is computationally more efficient and accurate.

Meesters, Antoon G. C. A.↗

Atomic Data and Spectral Line Intensities for CA XVII

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ca XVII. The configurations used are 2s(sup 2), 2s2p, 2p(sup 2), 2l3l', 214l' and 2s5l', with l = s,p and l' = s,p, d giving rise to 92 fine-structure levels in intermediate coupling. Collision strengths are calculated at seven incident energies (15, 30, 75, 112.5, 150, 187.5 and 225 Ry) for the transitions within the three lowest configurations corresponding to the 10 lowest energy levels, and five incident energies (75, 112.5, 150, 187.5 and 225 Ry) for transitions between the lowest five levels and the n = 3,4,5 configurations. Calculations have been carried out using the distorted wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates of the present work, and R-Matrix results for the 2s2, 2s2p, 2p2 configurations available in the literature, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8)-10(exp 14)/cu cm at an electron temperature of log Te(K)=6.7, corresponding to the maximum abundance of Ca XVII. Spectral line intensities are calculated, and their diagnostic relevance L; discussed. This dataset will be made available in the next version of the CHIANTI database.

Bhatia, A.K.↗