SEARCH · Engineering Papers
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.
SIRTF guaranteed time observers have specified field surveys to be carried out with IRAC with integration times per position ranging from 90 to over 10,000 seconds. This paper describes the plans for these surveys and their application to the study of galaxy formation and evolution.
This study applies an indicators framework to investigate climate drivers of tundra vegetation trends and variability over the 1982–2019 period. Previously known indicators relevant for tundra productivity (summer warmth index (SWI), coastal spring sea-ice (SI) area, coastal summer open-water (OW)) and three additional indicators (continentality, summer precipitation, and the Arctic Dipole (AD): second mode of sea level pressure variability) are analyzed with maximum annual Normalized Difference Vegetation Index (MaxNDVI) and the sum of summer bi-weekly (time-integrated) NDVI (TI-NDVI) from the Advanced Very High Resolution Radiometer time-series. Climatological mean, trends, and correlations between variables are presented. Changes in SI continue to drive variations in the other indicators. As spring SI has decreased, summer OW, summer warmth, MaxNDVI, and TI-NDVI have increased. However, the initial very strong upward trends in previous studies for MaxNDVI and TI-NDVI are weakening and becoming spatially and temporally more variable as the ice retreats from the coastal areas. TI-NDVI has declined over the last decade particularly over High Arctic regions and southwest Alaska. The continentality index (CI) (maximum minus minimum monthly temperatures) is decreasing across the tundra, more so over North America than Eurasia. The relationship has weakened between SI and SWI and TI-NDVI, as the maritime influence of OW has increased along with total precipitation. The winter AD is correlated in Eurasia with spring SI, summer OW, MaxNDVI, TI-NDVI, the CI and total summer precipitation. This winter connection to tundra emphasizes the role of SI in driving the summer indicators. The winter (DJF) AD drives SI variations which in turn shape summer OW, the atmospheric SWI and NDVI anomalies. The winter and spring indicators represent potential predictors of tundra vegetation productivity a season or two in advance of the growing season.
It has long been recognized that frequency and phasing of structural modes in the presence of airflow play a fundamental role in the occurrence of flutter. Animation of simulation results for the long, slender Semi-Span Super-Sonic Transport (S4T) wind-tunnel model demonstrates that, for the case of mass-ballasted nacelles, the flutter mode can be described as a traveling wave propagating downstream. Such a characterization provides certain insights, such as (1) describing the means by which energy is transferred from the airflow to the structure, (2) identifying airspeed as an upper limit for speed of wave propagation, (3) providing an interpretation for a companion mode that coalesces in frequency with the flutter mode but becomes very well damped, (4) providing an explanation for bursts of response to uniform turbulence, and (5) providing an explanation for loss of low frequency (lead) phase margin with increases in dynamic pressure (at constant Mach number) for feedback systems that use sensors located upstream from active control surfaces. Results from simulation animation, simplified modeling, and wind-tunnel testing are presented for comparison. The simulation animation was generated using double time-integration in Simulink of vertical accelerometer signals distributed over wing and fuselage, along with time histories for actuated control surfaces. Crossing points for a zero-elevation reference plane were tracked along a network of lines connecting the accelerometer locations. Accelerometer signals were used in preference to modal displacement state variables in anticipation that the technique could be used to animate motion of the actual wind-tunnel model using data acquired during testing. Double integration of wind-tunnel accelerometer signals introduced severe drift even with removal of both position and rate biases such that the technique does not currently work. Using wind-tunnel data to drive a Kalman filter based upon fitting coefficients to analytical mode shapes might provide a better means to animate the wind tunnel data.
The entire procedure of incorporating the requirements and goals of a space flight project into integrated, time ordered sequences of spacecraft commands, is called the uplink process. The Uplink Process Control Task (UPCT) was created to examine the uplink process and determine ways to improve it. The Problem Statement Language/Problem Statement Analyzer (PSL/PSA) designed to assist the designer/analyst/engineer in the preparation of specifications of an information system is used as a supporting tool to aid in the analysis. Attention is given to a definition of the uplink process, the definition of PSL/PSA, the construction of a PSA database, the value of analysis to the study of the uplink process, and the PSL/PSA lessons learned.
Modifications to a three dimensional, implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations are described. The modifications involve the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. The paper presents a detailed description of the time marching aeroelastic procedure and presents comparisons with experimental data to provide an assessment of the capability. Flutter results are shown for an isolated 45 degree swept-back wing and a supersonic transport configuration with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement, giving confidence in the accuracy of the aeroelastic capability that was developed.
Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.
Temperature - fuel-flow and temperature-area feedback control systems were investigated as means of controlling tailpipe gas temperature of a turbojet engine during transient operation in the high-speed region. Proportional-plus-integral control was used in both systems, but in the temperature-area control system it was necessary to add nonlinear components to the basic proportional-plus-integral control to provide satisfactory transient response to a desired step increase in temperature. Time integral of temperature-error functions were used as criteria for determining optimum transient response. A description of engine dynamics was obtained from frequency-response data.
Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.
We present results from the Medium Deep Survey (MDS), a Key Project using the Hubble Space Telescope (HST). Wide Field Camera (WFC) images of random fields have been taken in 'parallel mode' with an effective resolution of 0.2 sec full width at half maximum (FWHM) in the V(F555W) and I(F785LP) filters. The exposures presented here were targeted on a field away from 3C 273, and resulted in approximately 5 hr integration time in each filter. Detailed morphological structure is seen in galaxy images with total integrated magnitudes down to V approximately = 22.5 and I approximately = 21.5. Parameters are estimated that best fit the observed galaxy images, and 143 objects are identified (including 23 stars) in the field to a fainter limiting magnitude of I approximately = 23.5. We outline the extragalactic goals of the HST Medium Deep Survey, summarize our basic data reduction procedures, and present number (magnitude) counts, a color-magnitude diagram for the field, surface brightness profiles for the brighter galaxies, and best-fit half-light radii for the fainter galaxies as a function of apparent magnitude. A median galaxy half-light radius of 0.4 sec is measured, and the distribution of galaxy sizes versus magnitude is presented. We observe an apparent deficit of galaxies with half-light radii between approximately 0.6 sec and 1.5 sec, with respect to standard no-evolution or mild evolution cosmological models. An apparent excess of compact objects (half-light radii approximately 0.1 sec) is also observed with respect to those models. Finally, we find a small excess in the number of faint galaxy pairs and groups with respect to a random low-redshift field sample.
Because of the flow environment associated with the SSME powerhead pressure (3000 psia), temperature (1800 R), and mechanical complexity and the high vibration test stand environment, detailed flow measurements are difficult to make. The feasibility of using laser anemometry techniques to map velocity flow fields in an SSME powerhead is studied. In the study three engine powerhead component flow environments: (1) the high pressure fuel turbopump preburner, (2) the fuel turbopump turbine rotor and stator region, and (3) the 180 deg turnaround duct - are being considered. Flow parameters measured by the anemometry techniques are time averaged values of the velocity magnitude and flow direction, turbulence intensity, velocity component correlation, integral time scale, and turbulence spectrum.
Impact responses of composite laminates with and without initial stresses are investigated using the finite element method. A nine-node isoparametric quadrilateral element based on the Mindlin plate theory and the von Karman large deflection assumptions is developed. An experimentally established contact law which accounts for the permanent indentation is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjunction with successive iterations within each time step. Numerical results, including the contact force histories, deflections and strains in the plate, are presented.
The annual cycle of sea surface temperature (SST) in the southwestern Atlantic Ocean was estimated using four years (July 1984-July 1988) of NOAA Advanced Very High Resolution Radiometer observations. High resolution satellite observations at 1-km space and daily time resolution were grided at 100-km space and 5-day time intervals to develop an analysis dataset for determination of low frequency SST variability. The integral time scale, a measure of serial correlation, was found to vary from 40 to 60 days in the domain of interest. The existence of superannual trends in the SST data was investigated, but conclusive results could not be obtained. The annual cycle (and, in particular, the annual harmonic) explains a large proportion of the SST variability. The estimated amplitude of the cycle ranges between 5 deg and 13 deg C throughout the study area, with minima in August-September and maxima in February. The resultant climatology is compared with an arbitrary 5-day satellite SST field, and with the COADS/ICE SST climatology. It was found that the higher resolution satellite-based SST climatology resolves boundary current structure and has significantly better structural agreement with the observed field.
A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.
A typical heterodyne remote sensing system contains three major elements: the antenna, the radiometer, and the spectrometer. The radiometer consists of the local oscillator, the mixer, and the intermediate frequency amplifiers. This subsystem performs the function of down converting the high frequency incident thermal emission signal to a lower intermediate frequency. The spectrometer measures the power spectrum of the down-converted signal simultaneously in many contiguous frequency channels. Typical spectrum analysis requirements involve measurement of signal bandwidths of 100 to 1000 MHz with a channel resolution of 0.5 to 10 MHz. Three general approaches are used for spectrometers: (1) filter banks, (2) Acousto-Optic Spectrometers (AOS's), and (3) digital autocorrelators. In contrast to the two frequency domain techniques, an autocorrelator works in the time domain. The autocorrelation function (ACF) of the incoming signal is computed and averaged over the integration time. The averaged ACF is then Fourier transformed to obtain the signal power spectrum. Significant progress was made in the development of sub mm antennas and radiometers. It is now time to begin research in the development of low power spaceborne spectrometers and to reduce their size and weight. The near-term research goal will be to develop a prototype digital autocorrelation spectrometer, using VLSI gate array technology, which will have a small size, low power requirements, and can be used in spacecraft mm and sub mm radiometer systems. The long-range objective of this technology development is to make extremely low power, less than 10 mW/channel, small and stable wideband spectrometers which can be used in future mm and sub mm wavelength space missions such as the Large Deployable Reflector.
We show that a simple modification of Formulation 1 of Farassat results in a new analytic expression that is highly suitable for broadband noise prediction when extensive turbulence simulation is available. This result satisfies all the stringent requirements, such as permitting the use of the exact geometry and kinematics of the moving body, that we have set as our goal in the derivation of useful acoustic formulas for the prediction of rotating blade and airframe noise. We also derive a simple analytic expression for the autocorrelation of the acoustic pressure that is valid in the near and far fields. Our analysis is based on the time integral of the acoustic pressure that can easily be obtained at any resolution for any observer time interval and digitally analyzed for broadband noise prediction. We have named this result as Formulation 2B of Farassat. One significant consequence of Formulation 2B is the derivation of the acoustic velocity potential for the thickness and loading terms of the Ffowcs Williams-Hawkings (FW-H) equation. This will greatly enhance the usefulness of the Fast Scattering Code (FSC) by providing a high fidelity boundary condition input for scattering predictions.
An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.
A sample of 30 galaxies selected for their intense IRAS flux at 60 and 100 micron using the Arecibo telescope at 21 cm to measure the continuum and HI line luminosities were observed. The centimeter wave continuum correlates very well with the far-infrared flux, with a correlation coefficient as high as that found for other samples, and the same ratio between FIR and radio luminosities. Weaker correlations are seen between the FIR and optical luminosity and between the FIR and radio continuum. There is very little correlation between the FIR and the HI mass deduced from the integral of the 21 cm line. The strength of the radio continuum correlation suggests that there is little contribution to either the radio and FIR from physical processes not affecting both. If they each reflect time integrals of the star formation rate then the time constants must be similar, or the star formation rate must change slowly in these galaxies.
Current radiometric earth and atmospheric sensing systems in the centimeter wavelength range generally employ a directive antenna connected through a single terminal pair to a Dicke receiver. It is shown that this approach does not lend itself to systems with greatly increased spatial resolution. Signal to noise considerations relating to antenna efficiency force the introduction of active elements at the subarray level; thus, if Dicke switching is to be used, it must be distributed throughout the system. Some possible approaches are suggested. The introduction of active elements at the subarray level is found to ease the design constraints on time delay elements, necessary for bandwidth, and on multiple beam generation, required in order to achieve sufficient integration time with high resolution.