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 1,171 records · Page 65

Rolling-Element Fatigue Testing and Data Analysis - A Tutorial

In order to rank bearing materials, lubricants and other design variables using rolling-element bench type fatigue testing of bearing components and full-scale rolling-element bearing tests, the investigator needs to be cognizant of the variables that affect rolling-element fatigue life and be able to maintain and control them within an acceptable experimental tolerance. Once these variables are controlled, the number of tests and the test conditions must be specified to assure reasonable statistical certainty of the final results. There is a reasonable correlation between the results from elemental test rigs with those results obtained with full-scale bearings. Using the statistical methods of W. Weibull and L. Johnson, the minimum number of tests required can be determined. This paper brings together and discusses the technical aspects of rolling-element fatigue testing and data analysis as well as making recommendations to assure quality and reliable testing of rolling-element specimens and full-scale rolling-element bearings.

Vlcek, Brian L.↗

Review of the Aerodynamic Acceptance Test and Application to Anti-Icing Fluids Testing in the NRC Propulsion and Icing Wind Tunnel

In recent years, the FAA has worked with Transport Canada, National Research Council of Canada (NRC) and APS Aviation, Inc. to develop allowance times for aircraft operations in ice-pellet precipitation. These allowance times are critical to ensure safety and efficient operation of commercial and cargo flights. Wind-tunnel testing with uncontaminated anti-icing fluids and fluids contaminated with simulated ice-pellets had been carried out at the NRC Propulsion and Icing Wind Tunnel (PIWT) to better understand the flowoff characteristics and resulting aerodynamic effects. The percent lift loss on the thin, high-performance wing model tested in the PIWT was determined at 8 angle of attack and used as one of the evaluation criteria in determining the allowance times. Because it was unclear as to how performance degradations measured on this model were relevant to an actual airplane configuration, some means of interpreting the wing model lift loss was deemed necessary. In this report, the lift loss was related to the loss in maximum lift of a Boeing 737-200ADV airplane through the Aerodynamic Acceptance Test (AAT) performed for fluids qualification. This report provides a review of the research basis of the AAT in order to understand how this correlation was applied. A loss in maximum lift coefficient of 5.24 percent on the B737-200ADV airplane (which was adopted as the threshold in the AAT) corresponds to a lift loss of 7.3 percent on the PIWT model at 8 degrees angle of attack. There is significant scatter in the data used to develop the correlation related to varying effects of the various antiicing fluids that were tested and other factors. A statistical analysis indicated the upper limit of lift loss on the PIWT model was 9.2 percent. Therefore, for cases resulting in PIWT model lift loss from 7.3 to 9.2 percent, extra scrutiny of the visual observations is required in evaluating fluid performance with contamination. Additional research may result in future changes to this correlation.

Broeren, Andy P.↗

Detached-Eddy Simulation Based on the V2-F Model

Detached-eddy simulation (DES) based on the v(sup 2)-f Reynolds-averaged Navier-Stokes (RANS) model is developed and tested. The v(sup 2)-f model incorporates the anisotropy of near-wall turbulence which is absent in other RANS models commonly used in the DES community. The v(sup 2)-f RANS model is modified in order the proposed v(sup 2)-f-based DES formulation reduces to a transport equation for the subgrid-scale kinetic energy isotropic turbulence. First, three coefficients in the elliptic relaxation equation are modified, which is tested in channel flows with friction Reynolds number up to 2000. Then, the proposed v(sup 2)-f DES model formulation is derived. The constant, C(sub DES), required in the DES formulation was calibrated by simulating both decaying and statistically-steady isotropic turbulence. After C(sub DES) was calibrated, the v(sub 2)-f DES formulation is tested for flow around a circular cylinder at a Reynolds number of 3900, in which case turbulence develops after separation. Simulations indicate that this model represents the turbulent wake nearly as accurately as the dynamic Smagorinsky model. Spalart-Allmaras-based DES is also included in the cylinder flow simulation for comparison.

Jee, Sol Keun↗

An application of extremal field guidance in statistical low thrust guidance analysis

The subject of this paper is the application of extremal field guidance to the covariance analysis of low thrust retargeting guidance. The general procedure for low thrust navigation analysis is briefly reviewed. The basic equations of trajectory optimization and extremal field guidance are presented with emphasis on the low thrust problem. The adopted formulation includes thruster on-off times as control variables and accounts for the effects of model parameter variations. A method of solution of the guidance equations is outlined; and the application of the solution in a covariance analysis, including a procedure for the treatment of model parameter variations as either biases or first order Gauss-Markov stochastic processes, is discussed.

Jacobson, R. A.↗

A demonstration of a transportable radio interferometric surveying system with 3-cm accuracy on a 307-m base line

A precision geodetic measurement system (Aries, for Astronomical Radio Interferometric Earth Surveying) based on the technique of very long base line interferometry has been designed and implemented through the use of a 9-m transportable antenna and the NASA 64-m antenna of the Deep Space Communications Complex at Goldstone, California. A series of experiments designed to demonstrate the inherent accuracy of a transportable interferometer was performed on a 307-m base line during the period from December 1973 to June 1974. This short base line was chosen in order to obtain a comparison with a conventional survey with a few-centimeter accuracy and to minimize Aries errors due to transmission media effects, source locations, and earth orientation parameters. The base-line vector derived from a weighted average of the measurements, representing approximately 24 h of data, possessed a formal uncertainty of about 3 cm in all components. This average interferometry base-line vector was in good agreement with the conventional survey vector within the statistical range allowed by the combined uncertainties (3-4 cm) of the two techniques.

Ong, K. M.↗

Turbulent transport modelling of separating and reattaching shear flows

The improvement of capabilities for computer simulation of turbulent recirculating flows was investigated. Attention has been limited to two dimensional flows and principally to statistically stationary motion. Improvement of turbulence modeling explored the treatment of the near wall sublayer and of the exterior fully turbulent region, working within the framework of turbulence closures requiring the solution of transport equations for the turbulence energy and its dissipation rate. The work on the numerical procedure, based on the Gosman-Pun program TEACH, addressed the problems of incorporating the turbulence model as well as the extension to time dependent flows, the incorporation of a third order approximation of convective transport, and the treatment of non-orthogonal boundaries.

Launder, B. E.↗

The seasonal response of the Held-Suarez climate model to prescribed ocean temperature anomalies. II - Dynamical analysis

The heating associated with equatorial, subtropical, and midlatitude ocean temperature anamolies in the Held-Suarez climate model is analyzed. The local and downstream response to the anomalies is analyzed, first by examining the seasonal variation in heating associated with each ocean temperature anomaly, and then by combining knowledge of the heating with linear dynamical theory in order to develop a more comprehensive explanation of the seasonal variation in local and downstream atmospheric response to each anomaly. The extent to which the linear theory of propagating waves can assist the interpretation of the remote cross-latitudinal response of the model to the ocean temperature anomalies is considered. Alternative hypotheses that attempt to avoid the contradictions inherent in a strict application of linear theory are investigated, and the impact of sampling errors on the assessment of statistical significance is also examined.

Phillips, T. J.↗

Hemispheric differences in the morphology of the high latitude ionosphere measured at about 500 km

Ambient electron density measurements made by the radio frequency capacitance probe on Ariel 4 have been analyzed in order to study the structures present at about 500 km altitude at latitudes poleward of + or - 50 deg Lambda during the 1971 Northern Hemisphere winter (December solstice) and the 1972 Southern Hemisphere winter (June solstice); data obtained on 94 days centered on each solstice have been used. The general morphology of, and the extreme densities within, the midlatitude electron density trough, the polar cap depletions, and the electron density enhancements associated with the cusp-auroral zone were determined statistically. Analyses of data acquired during quiet magnetic conditions show that the Northern and Southern Hemisphere winter ionospheres were significantly different; in particular, the Southern Hemisphere densities were lower than those in the Northern Hemisphere. The maximum electron densities observed in the Northern Hemisphere occurred in a magnetic-local-time range symmetrical about the 02-14 MLT meridian, whereas in the Southern Hemisphere the maxima were symmetrical with respect to the midnight-noon magnetic meridian.

Tulunay, Y. K.↗

The NASA MSFC Earth Global Reference Atmospheric Model-2007 Version

Reference or standard atmospheric models have long been used for design and mission planning of various aerospace systems. The NASA/Marshall Space Flight Center (MSFC) Global Reference Atmospheric Model (GRAM) was developed in response to the need for a design reference atmosphere that provides complete global geographical variability, and complete altitude coverage (surface to orbital altitudes) as well as complete seasonal and monthly variability of the thermodynamic variables and wind components. A unique feature of GRAM is that, addition to providing the geographical, height, and monthly variation of the mean atmospheric state, it includes the ability to simulate spatial and temporal perturbations in these atmospheric parameters (e.g. fluctuations due to turbulence and other atmospheric perturbation phenomena). A summary comparing GRAM features to characteristics and features of other reference or standard atmospheric models, can be found Guide to Reference and Standard Atmosphere Models. The original GRAM has undergone a series of improvements over the years with recent additions and changes. The software program is called Earth-GRAM2007 to distinguish it from similar programs for other bodies (e.g. Mars, Venus, Neptune, and Titan). However, in order to make this Technical Memorandum (TM) more readable, the software will be referred to simply as GRAM07 or GRAM unless additional clarity is needed. Section 1 provides an overview of the basic features of GRAM07 including the newly added features. Section 2 provides a more detailed description of GRAM07 and how the model output generated. Section 3 presents sample results. Appendices A and B describe the Global Upper Air Climatic Atlas (GUACA) data and the Global Gridded Air Statistics (GGUAS) database. Appendix C provides instructions for compiling and running GRAM07. Appendix D gives a description of the required NAMELIST format input. Appendix E gives sample output. Appendix F provides a list of available parameters to enable the user to generate special output. Appendix G gives an example and guidance on incorporating GRAM07 as a subroutine in other programs such as trajectory codes or orbital propagation routines.

Leslie, F.W.↗

Time-Frequency Analysis of Rocket Nozzle Wall Pressures During Start-up Transients

Surveys of the fluctuating wall pressure were conducted on a sub-scale, thrust- optimized parabolic nozzle in order to develop a physical intuition for its Fourier-azimuthal mode behavior during fixed and transient start-up conditions. These unsteady signatures are driven by shock wave turbulent boundary layer interactions which depend on the nozzle pressure ratio and nozzle geometry. The focus however, is on the degree of similarity between the spectral footprints of these modes obtained from transient start-ups as opposed to a sequence of fixed nozzle pressure ratio conditions. For the latter, statistically converged spectra are computed using conventional Fourier analyses techniques, whereas the former are investigated by way of time-frequency analysis. The findings suggest that at low nozzle pressure ratios -- where the flow resides in a Free Shock Separation state -- strong spectral similarities occur between fixed and transient conditions. Conversely, at higher nozzle pressure ratios -- where the flow resides in Restricted Shock Separation -- stark differences are observed between the fixed and transient conditions and depends greatly on the ramping rate of the transient period. And so, it appears that an understanding of the dynamics during transient start-up conditions cannot be furnished by a way of fixed flow analysis.

Baars, Woutijn J.↗

Stochastic Reduced Order Models with Python (SROMPy)

Stochastic Reduced Order Models with Python (SROMPy) is a software package developed to enable user-friendly utilization of the stochastic reduced order model (SROM) approach for uncertainty quantification. A SROM is a low dimensional, discrete approximation to a random quantity that enables efficient and non-intrusive stochastic computations. With SROMPy, a user can easily generate a SROM to approximate a random variable or vector described by several different types of probability distributions using the Python programming language. Once a SROM is constructed, the software can be used to propagate uncertainty through a user-defined computational model to estimate statistics of a given quantity of interest. This report is meant to introduce the SROMPy module and brie y demonstrate its capabilities. A simple example of a spring-mass system with a random input is included to illustrate the practicality of the SROM approach to uncertainty quantification and relative ease of applying it with SROMPy. The example includes a comparison with a solution obtained using classical Monte Carlo simulation, demonstrating the similarities and advantages of using the SROM approach.

Warner, James E.↗

The central engine of quasars and AGNs: A relativistic proton radiative shock

Active galactic nuclei (AGNs) and quasars (QSOs) appear to emit roughly equal energy per decade from radio to gamma-ray energies (e.g. Ramaty and Ligenfelter 1982). This argues strongly for a nonthermal radiation mechanism (see Rees 1984). In addition, statistical studies have indicated that the spectra of these objects in the IR-UV and 2 to 50 keV X-ray band, can be fitted very well with power laws of specific indices. These spectral indices do not seem to depend on the luminosity or morphology of the objects (Rothschild et al. 1983; Malkan 1984), and any theory should account for them in a basic and model independent way. If shocks accelerate relativistic protons via the first-order Fermi mechanism (e.g. Axfor 1981), the radiating electrons can be produced as secondaries throughout the source by proton-proton (p-p) collisions and pion decay, thus eliminating Compton losses (Protheroe and Kazanas 1983). As shown by Kazanas (1984), if relativistic electrons are injected at high energies, e+-e- pair production results in a steady state electron distribution that is very similar to that observed in AGNs, independent of the details of injection and the dynamics of the source. The conditions required by this mechanism are met in the shock model of Eichler (1984) and Ellison and Eichler (1984) which allows the self-consistent calculation of the shock acceleration efficiency.

Kazanas, D.↗

Application of solar max ACRIM data to analyze solar-driven climatic variability on Earth

Terrestrial climatic effects associated with solar variability have been proposed for at least a century, but could not be assessed quantitatively owing to observational uncertainities in solar flux variations. Measurements from 1980 to 1984 by the Active Cavity Radiometer Irradiance Monitor (ACRIM), capable of resolving fluctuations above the sensible atmosphere less than 0.1% of the solar constant, permit direct albeit preliminary assessments of solar forcing effects on global temperatures during this period. The global temperature response to ACRIM-measured fluctuations was computed from 1980 to 1985 using the NYU transient climate model including thermal inertia effects of the world ocean; and compared the results with observations of recent temperature trends. Monthly mean ACRIM-driven global surface temperature fluctuations computed with the climate model are an order of magnitude smaller, of order 0.01 C. In constrast, global mean surface temperature observations indicate an approx. 0.1 C increase during this period. Solar variability is therefore likely to have been a minor factor in global climate change during this period compared with variations in atmospheric albedo, greenhouse gases and internal self-inducedoscillations. It was not possible to extend the applicability of the measured flux variations to longer periods since a possible correlation of luminosity with solar annual activity is not supported by statistical analysis. The continuous monitoring of solar flux by satellite-based instruments over timescales of 20 years or more comparable to timescales for thermal relaxation of the oceans and of the solar cycle itself is needed to resolve the question of long-term solar variation effects on climate.

Hoffert, M. I.↗

An improved approach for flight readiness assessment

An improved methodology for quantitatively evaluating failure risk for a spaceflight system in order to assess flight readiness is presented. This methodology is of particular value when information relevant to failure prediction, including test experience and knowledge of parameters used in engineering analyses of failure phenomena, is limited. In this approach, engineering analysis models that characterize specific failure modes based on the physics and mechanics of the failure phenomena are used in a prescribed probabilistic structure to generate a failure probability distribution that is modified by test and flight experience in a Bayesian statistical procedure. The probabilistic structure and statistical methodology are generally applicable to any failure mode for which quantitative engineering analysis can be employed to characterize the failure phenomenon and are particularly well suited for use under the constraints on information availability that are typical of such spaceflight systems as the Space Shuttle and planetary spacecraft.

Moore, N. R.↗

Three-dimensional position and shape of the bow shock and their variation with Alfvenic, sonic and magnetosonic Mach numbers and interplanetary magnetic field orientation

A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Peredo, M.↗

Polysilicon Gate Enhancement of the Random Dopant Induced Threshold Voltage Fluctuations in Sub-100 nm MOSFET's with Ultrathin Gate Oxide

In this paper, we investigate various aspects of the polysilicon gate influence on the random dopant induced threshold voltage fluctuations in sub-100 nm MOSFET's with ultrathin gate oxides. The study is done by using an efficient statistical three-dimensional (3-D) "atomistic" simulation technique described else-where. MOSFET's with uniform channel doping and with low doped epitaxial channels have been investigated. The simulations reveal that even in devices with a single crystal gate the gate depletion and the random dopants in it are responsible for a substantial fraction of the threshold voltage fluctuations when the gate oxide is scaled-in the range of 1-2 nm. Simulation experiments have been used in order to separate the enhancement in the threshold voltage fluctuations due to an effective increase in the oxide thickness associated with the gate depletion from the direct influence of the random dopants in the gate depletion layer. The results of the experiments show that the both factors contribute to the enhancement of the threshold voltage fluctuations, but the effective increase in the oxide-thickness has a dominant effect in the investigated range of devices. Simulations illustrating the effect or the polysilicon grain boundaries on the threshold voltage variation are also presented.

Asenov, Asen↗

Mathematical Inversion of Lightning Data: Techniques and Applications

A survey of some interesting mathematical inversion studies dealing with radio, optical, and electrostatic measurements of lightning are presented. A discussion of why NASA is interested in lightning, what specific physical properties of lightning are retrieved, and what mathematical techniques are used to perform the retrievals are discussed. In particular, a relatively new multi-station VHF time-of-arrival (TOA) antenna network is now on-line in Northern Alabama and will be discussed. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The LMA supports on-going ground-validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. The LMA also provides detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and offers interesting comparisons with other meteorological/geophysical datasets. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. A new channel mapping retrieval algorithm is introduced for this purpose. To characterize the spatial distribution of retrieval errors, the algorithm has been applied to analyze literally tens of millions of computer-simulated lightning VHF point sources that have been placed at various ranges, azimuths, and altitudes relative to the LMA network. Statistical results are conveniently summarized in high-resolution, color-coded, error maps.

Koshak, William↗

Statistical Analysis of CFD Solutions from 2nd Drag Prediction Workshop

In June 2001, the first AIAA Drag Prediction Workshop was held to evaluate results obtained from extensive N-Version testing of a series of RANS CFD codes. The geometry used for the computations was the DLR-F4 wing-body combination which resembles a medium-range subsonic transport. The cases reported include the design cruise point, drag polars at eight Mach numbers, and drag rise at three values of lift. Although comparisons of the code-to-code medians with available experimental data were similar to those obtained in previous studies, the code-to-code scatter was more than an order-of-magnitude larger than expected and far larger than desired for design and for experimental validation. The second Drag Prediction Workshop was held in June 2003 with emphasis on the determination of installed pylon-nacelle drag increments and on grid refinement studies. The geometry used was the DLR-F6 wing-body-pylon-nacelle combination for which the design cruise point and the cases run were similar to the first workshop except for additional runs on coarse and fine grids to complement the runs on medium grids. The code-to-code scatter was significantly reduced for the wing-body configuration compared to the first workshop, although still much larger than desired. However, the grid refinement studies showed no sign$cant improvement in code-to-code scatter with increasing grid refinement.

Hemsch, M. J.↗