Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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 811 records · Page 45

Gravity waves

Gravity waves contributed to the establishment of the thermal structure, small scale (80 to 100 km) fluctuations in velocity (50 to 80 m/sec) and density (20 to 30%, 0 to peak). Dominant gravity wave spectrum in the middle atmosphere: x-scale, less than 100 km; z-scale, greater than 10 km; t-scale, less than 2 hr. Theorists are beginning to understand middle atmosphere motions. There are two classes: Planetary waves and equatorial motions, gravity waves and tidal motions. The former give rise to variability at large scales, which may alter apparent mean structure. Effects include density and velocity fluctuations, induced mean motions, and stratospheric warmings which lead to the breakup of the polar vortex and cooling of the mesosphere. On this scale are also equatorial quasi-biennial and semi-annual oscillations. Gravity wave and tidal motions produce large rms fluctuations in density and velocity. The magnitude of the density fluctuations compared to the mean density is of the order of the vertical wavelength, which grows with height. Relative density fluctuations are less than, or of the order of 30% below the mesopause. Such motions may cause significant and variable convection, and wind shear. There is a strong seasonal variation in gravity wave amplitude. Additional observations are needed to address and quantify mean and fluctuation statistics of both density and mean velocity, variability of the mean and fluctuations, and to identify dominant gravity wave scales and sources as well as causes of variability, both temporal and geographic.

Fritts, David↗

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

Satellite Image Mosaic Engine

A computer program automatically builds large, full-resolution mosaics of multispectral images of Earth landmasses from images acquired by Landsat 7, complete with matching of colors and blending between adjacent scenes. While the code has been used extensively for Landsat, it could also be used for other data sources. A single mosaic of as many as 8,000 scenes, represented by more than 5 terabytes of data and the largest set produced in this work, demonstrated what the code could do to provide global coverage. The program first statistically analyzes input images to determine areas of coverage and data-value distributions. It then transforms the input images from their original universal transverse Mercator coordinates to other geographical coordinates, with scaling. It applies a first-order polynomial brightness correction to each band in each scene. It uses a data-mask image for selecting data and blending of input scenes. Under control by a user, the program can be made to operate on small parts of the output image space, with check-point and restart capabilities. The program runs on SGI IRIX computers. It is capable of parallel processing using shared-memory code, large memories, and tens of central processing units. It can retrieve input data and store output data at locations remote from the processors on which it is executed.

Plesea, Lucian↗

Determination of the Limiting Magnitude

The limiting magnitude of an optical camera system is an important property to understand since it is used to find the completeness limit of observations. Limiting magnitude depends on the hardware and software of the system, current weather conditions, and the angular speed of the objects observed. If an object exhibits a substantial angular rate during the exposure, its light spreads out over more pixels than the stationary stars. This spreading causes the limiting magnitude to be brighter when compared to the stellar limiting magnitude. The effect, which begins to become important when the object moves a full width at half max during a single exposure or video frame. For targets with high angular speeds or camera systems with narrow field of view or long exposures, this correction can be significant, up to several magnitudes. The stars in an image are often used to measure the limiting magnitude since they are stationary, have known brightness, and are present in large numbers, making the determination of the limiting magnitude fairly simple. In order to transform stellar limiting magnitude to object limiting magnitude, a correction must be applied accounting for the angular velocity. This technique is adopted in meteor and other fast-moving object observations, as the lack of a statistically significant sample of targets makes it virtually impossible to determine the limiting magnitude before the weather conditions change. While the weather is the dominant factor in observing satellites, the limiting magnitude for meteors also changes throughout the night due to the motion of a meteor shower or sporadic source radiant across the sky. This paper presents methods for determining the limiting stellar magnitude and the conversion to the target limiting magnitude.

Kingery, Aaron↗

Landslide Hazard is Projected to Increase Across High Mountain Asia

High Mountain Asia has long been known as a hotspot for landslide risk, and studies have suggested that landslide hazard is likely to increase in this region over the coming decades. Extreme precipitation may become more frequent, with a nonlinear response relative to increasing global temperatures. However, these changes are geographically varied. This article maps probable changes to landslide hazard, as shown by a landslide hazard indicator (LHI) derived from downscaled precipitation and temperature. In order to capture the nonlinear response of slopes to extreme precipitation, a simple machine-learning model was trained on a database of landslides across High Mountain Asia to develop a regional LHI. This model was applied to statistically downscaled data from the 30 members of the Seamless System for Prediction and Earth System Research large ensembles to produce a range of possible outcomes under the Shared Socioeconomic Pathways 2-4.5 and 5-8.5. The LHI reveals that landslide hazard will increase in most parts of High Mountain Asia. Absolute increases will be highest in already hazardous areas such as the Central Himalaya, but relative change is greatest on the Tibetan Plateau. Even in regions where landslide hazard declines by year 2100, it will increase prior to the mid-century mark. However, the seasonal cycle of landslide occurrence will not change greatly across High Mountain Asia. Although substantial uncertainty remains in these projections, the overall direction of change seems reliable. These findings highlight the importance of continued analysis to inform disaster risk reduction strategies for stakeholders across High Mountain Asia.

Thomas A Stanley↗

Probing me Reverse Shock in an Oxygen-Rich Supernova Remnant

The aim of this project is to examine the O VI emission at three positions around the X- ray bright ring of the remnant in order to investigate the relation between the O VI emission, the X-ray O VII and O VIII emission, and the optical [OIII} emission, and how these vary around the rim of the remnant. All three pointings and the background pointing have now been observed; the archive notification for the most recent dataset was Oct 30, 2003. After reprocessing and screening, the net exposure time for the SE exposure is only 54 percent of the approved time (15 kilosec). for the SE exposure, the available statistics are not good enough for analysis. A request for reobservation to make up for the lost time in the SE pointing has been approved. Broad O VI 1032 and O VI 1038 emission is detected with velocity width of at least 800 km/s, and possibly exceeding 1000 km/s. The Flanagan et al. analysis of the Chandra grating data show bulk velocities in the X-ray gas of order +/- 1000 km/s. In the region of the FUSE E0lO2-SE pointing, the Chandra data indicate both blue-shifted and red-shifted emission. Analysis of the velocity structure of the O VI emission will provide additional constraints on the kinematics of the gas: is it emission from a tilted expanding barrel, or a more symmetric expansion? The O VI fluxes are also needed to assess whether the O VI is radiation from recombining O VII or instead fiom cooler gas ionizing toward O VII. The emission is faint, however, which complicates the analysis. Because the lines are so broad, absorption by intervening H2, CII, and foreground OVI must be considered. A number of stars in the SMC have been observed which provide information on foreground OVI absorption. The initial analyses have concentrated on the Li1F channel since the guidance is based on that channel. The Li2F channel is being examined exposure by exposure to see if any of the data can be used to improve the signal to noise in the Li1F data.

Gaetz, Terrance↗

Organometallic Polymeric Conductors

For aerospace applications, the use of polymers can result in tremendous weight savings over metals. Suitable polymeric materials for some applications like EMI shielding, spacecraft grounding, and charge dissipation must combine high electrical conductivity with long-term environmental stability, good processability, and good mechanical properties. Recently, other investigators have reported hybrid films made from an electrically conductive polymer combined with insulating polymers. In all of these instances, the films were prepared by infiltrating an insulating polymer with a precursor for a conductive polymer (either polypyrrole or polythiophene), and oxidatively polymerizing the precursor in situ. The resulting composite films have good electrical conductivity, while overcoming the brittleness inherent in most conductive polymers. Many aerospace applications require a combination of properties. Thus, hybrid films made from polyimides or other engineering resins are of primary interest, but only if conductivities on the same order as those obtained with a polystyrene base could be obtained. Hence, a series of experiments was performed to optimize the conductivity of polyimide-based composite films. The polyimide base chosen for this study was Kapton. 3-MethylThiophene (3MT) was used for the conductive phase. Three processing variables were identified for producing these composite films, namely time, temperature, and oxidant concentration for the in situ oxidation. Statistically designed experiments were used to examine the effects of these variables and synergistic/interactive effects among variables on the electrical conductivity and mechanical strength of the films. Multiple linear regression analysis of the tensile data revealed that temperature and time have the greatest effect on maximum stress. The response surface of maximum stress vs. temperature and time (for oxidant concentration at 1.2 M) is shown. Conductivity of the composite films was measured for over 150 days in air at ambient temperature. The conductivity of the films dropped only half an order of magnitude in that time. Films aged under vacuum at ambient temperature diminished slightly in conductivity in the first day, but did not change thereafter. An experimental design approach will be applied to maximize the efficiency of the laboratory effort. The material properties (initial and long term) will also be monitored and assessed. The experimental results will add to the existing database for electrically conductive polymer materials. Attachments: 1) Synthesis Crystal Structure, and Polymerization of 1,2:5,6:9,10-Tribenzo-3,7,11,13-tetradehydro(14) annulene. 2) Reinvestigation of the Photocyclization of 1,4-Phenylene Bis(phenylmaleic anhydride): Preparation and Structure of (5)Helicene 5,6:9,10-Dianhydride. 3) Preparation and Structure Charecterization of a Platinum Catecholate Complex Containing Two 3-Ethynyltheophone Groups. and 4) Rigid-Rod Polymers Based on Noncoplanar 4,4'-Biphenyldiamines: A Review of Polymer Properties vs Configuration of Diamines.

Youngs, Wiley J.↗

A statistical mechanical calculation of the thermodynamic properties of interstitial solid solutions involving second nearest neighbor interactions.

A previous calculation of the thermodynamic properties of interstitial solid solutions based on the technique of Kirkwood expansions has been extended to include the effects of second nearest neighbor solute atom mutual interactions. The error inherent in the first order (or quasi-chemical) counting of the degeneracy of the solution crystal is avoided. It is shown that, at high temperatures, even strong second nearest neighbor solute mutual interactions have a negligible effect on the entropy of the solution and a small, temperature-dependent effect on the solute partial enthalpy.

Alex, K.↗

Bayesian recursive image estimation.

The enhancement of images that are characterized only by statistical data where the picture contains additive noise is considered. The random process representing the output of the scanner is characterized by the output of a dynamic system with white noise input. The dynamic system describes a first-order vector-Markov process. The procedure of Kalman filtering is then utilized to recursively determine the minimum mean-square error estimate of the image. The result is then extended to obtain the smoothing of the data.

Nahi, N. E.↗

Linear and circular polarization of the cosmic background radiation

New data which consist of continued measurements of the linear polarization of the cosmic background radiation as well as the first measurement of the circular polarization are examined. Eleven declinations have been surveyed for linear polarization and one declination for circular polarization, all at 9 mm wavelength. No evidence is found for either a significant linear or circular component with statistical errors on the linear component of 20-60 microK for various models. For linear polarization, a 95 percent confidence level limit of 0.1 mK (0.00003) for an axisymmetric anisotropic model is achieved, while for spherical harmonics through third order, a corresponding limit of 0.2 mK is achieved. For a declination of 37 deg, a limit of 12 mK is placed on the time-varying component and 20 mK on the dc component of the circular polarization at the 95 percent confidence level. At 37 percent declination, the sensitivity per beam patch (7 deg) is 0.2 mK.

Lubin, P.↗

Coronal propagation of flare associated electrons and protons

A statistical study of characteristic times and intensities of 36 solar particle events observed between 1977 and 1979 by the Kiel Cosmic Ray Experiment on board HELIOS-1 and -2 has been carried out. For approx. 0.5 MeV electrons we order the times of maximum and the absolute intensities with respect to angular distance from the parent flare. Discussion of coronal parameters in terms of Reid's model leads to typical time constants for coronal diffusion and escape.

Schellert, G.↗

Error modeling for differential GPS

Differential Global Positioning System (DGPS) positioning is used to accurately locate a GPS receiver based upon the well-known position of a reference site. In utilizing this technique, several error sources contribute to position inaccuracy. This thesis investigates the error in DGPS operation and attempts to develop a statistical model for the behavior of this error. The model for DGPS error is developed using GPS data collected by Draper Laboratory. The Marquardt method for nonlinear curve-fitting is used to find the parameters of a first order Markov process that models the average errors from the collected data. The results show that a first order Markov process can be used to model the DGPS error as a function of baseline distance and time delay. The model's time correlation constant is 3847.1 seconds (1.07 hours) for the mean square error. The distance correlation constant is 122.8 kilometers. The total process variance for the DGPS model is 3.73 sq meters.

Blerman, Gregory S.↗

Attitude Estimation Signal Processing: A First Report on Possible Algorithms and Their Utility

In this brief effort, time has been of the essence. The data had to be acquired from APL/Lincoln Labs, stored, and sorted out to obtain the pertinent streams. This has been a significant part of this effort and hardware and software problems have been addressed with the appropriate solutions to accomplish this part of the task. Passed this, some basic and important algorithms are utilized to improve the performance of the attitude estimation systems. These algorithms are an essential part of the signal processing for the attitude estimation problem as they are utilized to reduce the amount of the additive/multiplicative noise that in general may or may not change its structure and probability density function, pdf, in time. These algorithms are not currently utilized in the processing of the data, at least, we are not aware of their use in this attitude estimation problem. Some of these algorithms, like the variable thresholding, are new conjectures, but one would expect that someone somewhere must have utilized this kind of scheme before. The variable thresholding idea is a straightforward scheme to use in case of a slowly varying pdf, or statistical moments of the unwanted random process. The algorithms here are kept simple but yet effective for processing the data and removing the unwanted noise. For the most part, these algorithms can be arranged so that their consecutive and orderly execution would complement the preceding algorithm and improve the overall performance of the signal processing chain.

Riasati, Vahid R.↗

Improved Modeling in a Matlab-Based Navigation System

An innovative approach to autonomous navigation is available for low earth orbit satellites. The system is developed in Matlab and utilizes an Extended Kalman Filter (EKF) to estimate the attitude and trajectory based on spacecraft magnetometer and gyro data. Preliminary tests of the system with real spacecraft data from the Rossi X-Ray Timing Explorer Satellite (RXTE) indicate the existence of unmodeled errors in the magnetometer data. Incorporating into the EKF a statistical model that describes the colored component of the effective measurement of the magnetic field vector could improve the accuracy of the trajectory and attitude estimates and also improve the convergence time. This model is identified as a first order Markov process. With the addition of the model, the EKF attempts to identify the non-white components of the noise allowing for more accurate estimation of the original state vector, i.e. the orbital elements and the attitude. Working in Matlab allows for easy incorporation of new models into the EKF and the resulting navigation system is generic and can easily be applied to future missions resulting in an alternative in onboard or ground-based navigation.

Deutschmann, Julie↗

Alternative Test Methods for Electronic Parts

It is common practice within NASA to test electronic parts at the manufacturing lot level to demonstrate, statistically, that parts from the lot tested will not fail in service using generic application conditions. The test methods and the generic application conditions used have been developed over the years through cooperation between NASA, DoD, and industry in order to establish a common set of standard practices. These common practices, found in MIL-STD-883, MIL-STD-750, military part specifications, EEE-INST-002, and other guidelines are preferred because they are considered to be effective and repeatable and their results are usually straightforward to interpret. These practices can sometimes be unavailable to some NASA projects due to special application conditions that must be addressed, such as schedule constraints, cost constraints, logistical constraints, or advances in the technology that make the historical standards an inappropriate choice for establishing part performance and reliability. Alternate methods have begun to emerge and to be used by NASA programs to test parts individually or as part of a system, especially when standard lot tests cannot be applied. Four alternate screening methods will be discussed in this paper: Highly accelerated life test (HALT), forward voltage drop tests for evaluating wire-bond integrity, burn-in options during or after highly accelerated stress test (HAST), and board-level qualification.

Plante, Jeannette↗

Entry, Descent, and Landing technological barriers and crewed MARSvehicle performance analysis

Mars has been explored historically only by robotic crafts, but a crewed mission encompasses several new engineering challenges – high ballistic coefficient entry, hypersonic decelerators, guided entry for reaching intended destinations within acceptable margins for error in the landing ellipse, and payload mass are all critical factors for evaluation. A comprehensive EDL parametric analysis has been conducted in support of a high mass landing architecture by evaluating three types of vehicles −70° Sphere Cone, Ellipsled and SpaceX hybrid architecture called Red Dragon as potential candidate options for crewed entry vehicles. Aerocapture at the Martian orbit of about 400 km and subsequent Entry-from-orbit scenarios were investigated at velocities of 6.75 km/s and 4 km/s respectively. A study on aerocapture corridor over a range of entry velocities (6–9 km/s) suggests that a hypersonic L/D of 0.3 is sufficient for a Martian aerocapture. Parametric studies conducted by varying aeroshell diameters from 10 m to 15 m for several entry masses up to 150 mt are summarized and results reveal that vehicles with entry masses in the range of about 40–80 mt are capable of delivering cargo with a mass on the order of 5–20 mt. For vehicles with an entry mass of 20 mt to 80 mt, probabilistic Monte Carlo analysis of 5000 cases for each vehicle were run to determine the final landing ellipse and to quantify the statistical uncertainties associated with the trajectory and attitude conditions during atmospheric entry. Strategies and current technological challenges for a human rated Entry, Descent, and Landing to the Martian surface are presented in this study.

Prabhakar Subrahmanyam↗

Statistical Properties of SGR 1900+14 Bursts

We study the statistics of soft gamma repeater (SGR) bursts using a database of 187 events detected with BATSE and 837 events detected with the Rossi X-Ray Timing Explorer Proportional Counter Array: all events are from SGR 1900+14 during its 1998-1999 active phase. We find that the fluence or energy distribution of bursts is consistent with a power law of index 1.66, over 4 orders of magnitude. This scale-free distribution resembles the Gutenberg-Richter law for earthquakes and gives evidence for self-organized criticality in SGRS. The distribution of time intervals between successive bursts from SGR 1900+14 is consistent with a lognormal distribution. There is no correlation between burst intensity and the waiting times till the next burst, but there is some evidence for a correlation between burst intensity and the time elapsed since the previous burst. We also find a correlation between the duration and the energy of the bursts, but with significant scatter. In all these statistical properties, SGR bursts resemble earthquakes and solar flares more closely than they resemble any known accretion-powered or nuclear-powered phenomena. Thus, our analysis lends support to the hypothesis that the energy source for SGR bursts is internal to the neutron star and plausibly magnetic.

Gogus, Ersin↗

Diffusion flame in homologous turbulent shear flows.

The simplifed statistical theory developed previously is employed to analyze the equilibrium and near-equilibrium combustion of initially unmixed reactants. It is found that the flame zone in the limit of large Damkohler numbers is very thick and is of the order of the local integral scale of turbulence. This is in contrast to the existing phenomenological theories which predict the infinitesimally thin flame sheet, in the same limit, as it is with the laminar diffusion flame. Qualitative agreements with the available experimental results are shown. It is found that singularities exist at the edges of the flame which are removed as Damkohler number is reduced. Also, it is found that the heat transfer may take place against the local mean temperature gradient in certain regions within the flame.

Chung, P. M.↗