Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Extreme value analysis”

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 55 records · Page 3

3-D Surface Depression Profiling Using High Frequency Focused Air-Coupled Ultrasonic Pulses

Surface topography is an important variable in the performance of many industrial components and is normally measured with diamond-tip profilometry over a small area or using optical scattering methods for larger area measurement. This article shows quantitative surface topography profiles as obtained using only high-frequency focused air-coupled ultrasonic pulses. The profiles were obtained using a profiling system developed by NASA Glenn Research Center and Sonix, Inc (via a formal cooperative agreement). (The air transducers are available as off-the-shelf items from several companies.) The method is simple and reproducible because it relies mainly on knowledge and constancy of the sound velocity through the air. The air transducer is scanned across the surface and sends pulses to the sample surface where they are reflected back from the surface along the same path as the incident wave. Time-of-flight images of the sample surface are acquired and converted to depth/surface profile images using the simple relation (d = V*t/2) between distance (d), time-of-flight (t), and the velocity of sound in air (V). The system has the ability to resolve surface depression variations as small as 25 microns, is useable over a 1.4 mm vertical depth range, and can profile large areas only limited by the scan limits of the particular ultrasonic system. (Best-case depth resolution is 0.25 microns which may be achievable with improved isolation from vibration and air currents.) The method using an optimized configuration is reasonably rapid and has all quantitative analysis facilities on-line including 2-D and 3-D visualization capability, extreme value filtering (for faulty data), and leveling capability. Air-coupled surface profilometry is applicable to plate-like and curved samples. In this article, results are shown for several proof-of-concept samples, plastic samples burned in microgravity on the STS-54 space shuttle mission, and a partially-coated cylindrical ceramic composite sample. Impressive results were obtained for all samples when compared with diamond-tip profiles and measurements from micrometers. The method is completely nondestructive, noninvasive, non-contact and does not require light-reflective surfaces.

Roth, Don J.

Strain-energy-release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.

Strain energy release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtual crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.

Flow field analysis

The average mean wind speed integrated over a disk is shown to be extremely close to the mean value of wind speed which would be measured at the center of a disk for most geometries in which a WECS (Wind Energy Conversion System) would operate. Field test results are presented which compare instantaneous records of wind speed integrated over a disk with the wind speed measured at the center of the disk. The wind field that a rotating element would experience is presented which was synthesized from the outputs of an array of anemometers.

Cliff, W. C.

An analysis of particle track effects on solid mammalian tissues

Relative biological effectiveness (RBE) and quality factor (Q) at extreme values of linear energy transfer (LET) have been determined on the basis of experiments with single-cell systems and specific tissue responses. In typical single-cell systems, each heavy particle (Ar or Fe) passes through a single cell or no cell. In experiments on animal tissues, however, each heavy particle passes through several cells, and the LET can exceed 200 keV micrometers-1 in every cell. In most laboratory animal tissue systems, however, only a small portion of the hit cells are capable of expressing the end-point being measured, such as cell killing, mutation or carcinogenesis. The following question was therefore addressed: do RBEs and Q factors derived from single-cell experiments properly account for the damage at high LET when multiple cells are hit by HZE tracks? A review is offered in which measured radiation effects and known tissue properties are combined to estimate on the one hand, the number of cells at risk, p3n, per track, where n is the number of cells per track based on tissue and organ geometry, and p3 is the probability that a cell in the track is capable of expressing the experimental end-point. On the other hand, the tissue and single-cell responses are compared by determining the ratio RBE in tissue/RBE in corresponding single cells. Experimental data from the literature indicate that tissue RBEs at very high LET (Fe and Ar ions) are higher than corresponding single-cell RBEs, especially in tissues in which p3n is high.

Non-NASA Center

Evaluation of Global Climate Models for Use in Energy Analysis

The interplay between energy, climate, and weather is becoming more complex due to increasing contributions of renewable energy generation, energy storage, electrified end uses, and the increasing frequency of extreme weather events. Energy system analyses commonly rely on meteorological inputs to estimate renewable energy generation and energy demand; however, these inputs rarely represent the estimated impacts of future climate change. Climate models and publicly available climate change datasets can be used for this purpose, but the selection of inputs from the myriad of available models and datasets is a nuanced and subjective process. In this work, we assess datasets from various global climate models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). We present evaluations of their skills with respect to the historical climate and comparisons of their future projections of climate change for two climate change scenarios. We present the results for different climatic and energy system regions and include interactive figures in the accompanying software repository. Previous work has presented similar GCM evaluations, but none have presented variables and metrics specifically intended for comprehensive energy systems analysis including impacts on energy demand, thermal cooling, hydropower, water availability, solar energy generation, and wind energy generation. We focus on GCM output meteorological variables that directly affect these energy system components including the representation of extreme values that can drive grid resilience events. The objective of this work is not to recommend the best climate model and dataset for a given analysis, but instead to provide a reference to facilitate the selection of climate models and scenarios in subsequent work.

14 SOLAR ENERGY

A Technique Using Principal Component Analysis to Compare Seasonal Cycles of Earth Radiation from CERES and Model Computations

A method for quantitatively comparing the seasonal cycles of two global data sets is presented. The seasonal cycles of absorbed solar radiation (ASR) and outgoing longwave radiation (OLR) have been computed from an eight-year data set from the Clouds and Earth's Radiant Energy System (CERES) scanning radiometers and from a model data set produced by the NASA Goddard Space Flight Center Global Modeling and Assimilation Office. To compare the seasonal cycles from these two data sets, principal component (PC) analysis is used, where the PCs express the time variations and the corresponding empirical orthogonal functions (EOFs) describe the geographic variations. Ocean has a long thermal response time compared to land, so land and ocean are separated for the analysis. The root-mean square values for the seasonal cycles of ASR and OLR are extremely close for the two data sets. The first three PCs are quite close, showing that the time responses and magnitudes over the globe are very similar. The agreement between the two sets of PCs is quantified by computing the matrix of inner products of the two sets. For ASR over land, the first PCs of CERES and the model agree to better than 99.9%. The EOF maps are similar for most of the globe, but differ in a few places, and the agreement of the EOF maps is likewise quantified. Maps of differences between the annual cycles show regions of agreement and disagreement.

Smith, G. Louis

Monitoring physical and chemical parameters of Delaware Bay waters with an ERTS-1 data collection platform

The author has identified the following significant results. Results of the analysis of data collected during the summer of 1974 demonstrate that the ERTS Data Collection Platform (DCP) is quite responsive to changing water parameters and that this information can be successfully transmitted under all weather conditions. The monitoring of on-site probe outputs reveals a rapid response to changing water temperature, salinity, and turbidity conditions on incoming tides as the tidal salt wedge passes the probe location. The changes in water properties were corroborated by simultaneously sampling the water for subsequent laboratory analysis. Fluctuations observed in the values of salinity, conductivity, temperature and water depth over short time intervals were extremely small. Due to the nature of the probe, 10% to 20% fluctuations were observed in the turbidity values. The use of the average of the values observed during an overpass provided acceptable results. Good quality data was obtained from the satellite on each overpass regardless of weather conditions. Continued use of the DCP will help provide an indication of the accuracy of the probes and transmission system during long term use.

Klemas, V.

Guidelines for the Selection of Near-Earth Thermal Environment Parameters for Spacecraft Design

Thermal analysis and design of Earth orbiting systems requires specification of three environmental thermal parameters: the direct solar irradiance, Earth's local albedo, and outgoing longwave radiance (OLR). In the early 1990s data sets from the Earth Radiation Budget Experiment were analyzed on behalf of the Space Station Program to provide an accurate description of these parameters as a function of averaging time along the orbital path. This information, documented in SSP 30425 and, in more generic form in NASA/TM-4527, enabled the specification of the proper thermal parameters for systems of various thermal response time constants. However, working with the engineering community and SSP-30425 and TM-4527 products over a number of years revealed difficulties in interpretation and application of this material. For this reason it was decided to develop this guidelines document to help resolve these issues of practical application. In the process, the data were extensively reprocessed and a new computer code, the Simple Thermal Environment Model (STEM) was developed to simplify the process of selecting the parameters for input into extreme hot and cold thermal analyses and design specifications. In the process, greatly improved values for the cold case OLR values for high inclination orbits were derived. Thermal parameters for satellites in low, medium, and high inclination low-Earth orbit and with various system thermal time constraints are recommended for analysis of extreme hot and cold conditions. Practical information as to the interpretation and application of the information and an introduction to the STEM are included. Complete documentation for STEM is found in the user's manual, in preparation.

Anderson, B. J.

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS

Non-Linear Slosh Damping Model Development and Validation

Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can lead to significant savings by reducing the number and size of slosh baffles in liquid propellant tanks.

Yang, H. Q.

Gravitational influences on the behavior of confined diffusion flames

Closed form, mathematical representations for the shapes of laminar diffusion flames stabilized on two-dimensional slot burners are derived for arbitrary gravitational levels. The solution technique employs the Schvab-Zeldovich conserved species variable. Effects of gravity on the velocity field are incorporated through the momentum equation; these effects, in turn, influence the species field. The results reduce smoothly to the classical Burke-Schumann solution as gravity approaches zero. It is shown that both flame height and width increase with decreasing gravity levels. Several interesting cases are considered explicitly; these include effects of negative and extremely high values of gravity. In the latter case, it is shown that the flame collapses onto the burner surface. The analysis finds application in the modeling of steady and unsteady gas and solid propellant diffusion flames.

Hegde, U.

Impact of survey spatial variability on galaxy redshift distributions and the cosmological 3 × 2-point statistics for the Rubin Legacy Survey of Space and Time (LSST)

We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry data set degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-z estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i\lt 25.3$ for 10 yr LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd i-band depth for 1-yr (Y1), 3-yr (Y3), and 5-yr (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10\!-\!40~{{\rm per\ cent}}$ at extreme depth values. Using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times 2$ pt analysis, and assess the impact on cosmological parameters via a Fisher forecast. We find that galaxy clustering is most susceptible to variable depth, and non-uniformity needs to be mitigated below 3 per cent to recover unbiased cosmological constraints. There is little impact on galaxy–shear and shear–shear power spectra, given the expected LSST Y3 noise.

cosmology

RINA: An interactive system for the rapid interpretation of NASTRAN results

RINA is a general post-processor system which performs data reduction, post processing and post analysis of NASTRAN results. RINA can be run in either batch or time-sharing mode. The system, which is an expansion of the NASDAT program, performs the following tasks: (1) scan of extreme values for displacements, forces, stresses and margins of safety, (2) computations of envelopes for displacements, forces, stresses and margins of safety, (3) computations of allowables and margins of safety; and (4) generation of NASTRAN tables needed for reanalysis computations. Usage of all these options provides the analyst with an efficient tool for the study and interpretation of NASTRAN analysis results and their presentation for project documentation.

Raibstone, A. I.

Efficient Jitter Analysis for Spacecraft

Typically in space missions, the science instruments require a specific degree of pointing accuracy as well as dynamical quietness. This dynamical quietness, which is needed to allow the instruments to make measurements (remote sensing applications) or perform other functions, is usually characterized in terms of jitter and stability specifications. In order to insure that the spacecraft meets the requirements of its instruments, several jitter analyses are performed throughout the design phase of the spacecraft and beyond as the models of the spacecraft, its components, and disturbances mature. Each such analysis involves the simulation of the spacecraft and instrument dynamical response to all known disturbance scenarios, followed by the computation of jitter values for each instrument based on the specified jitter time windows. The direct approach for computing jitter values by sweeping maxima and minima throughout the time history may be costly in the computational sense as the size and number of the time histories involved could be quite large. Keeping in mind that typical spacecraft simulation time histories may easily involve hundreds of thousands or millions of points, it is imperative that a jitter analysis algorithm be developed which is more efficient than the direct approach. This paper presents a vectorized algorithm for efficient computation of spacecraft jitter values. The algorithm identifies the extreme points in the time history, which are the points that may dominate the jitter values depending on the location of the jitter window along the time history. The span of influence of each extremum is then computed by the algorithm and used in an efficient and vectorized fashion to obtain the jitter values. The algorithm deals with the multiple jitter windows 2 sequentially, first computing jitter values for the smallest time window, then looping over the remaining time windows until all jitter values are computed. A numerical example is carried out to demonstrate the efficiency and feasibility of the proposed jitter analysis technique.

Maghami, Peiman G.

The Physics of Turbulence in the Boundary Layer

The geometry of the velocity field in a numerically simulated incompressible turbulent boundary layer over a flat plate at Re theta=670 has been studied using the invariants of the velocity gradient tensor. These invariants are computed at every grid point in the flow and used to form the discriminant. Of primary interest are those regions in the flow where the discriminant is positive; regions where, according to the characteristic equation, the eigenvalues of the velocity gradient tensor are complex. An observer moving with a frame of reference which is attached to a fluid particle lying within such a region would see a local flow pattern of the type stable-focus-stretching or unstable-focus-compressing. When the flow is visualized this way, continuous, connected, large-scale structures are revealed that extend from the point just below the buffer layer out to the beginning of the wake region. These structures are aligned with the mean shear close to the wall and arch in the cross-stream direction away from the wall. In some cases the structures observed are very similar to to the hairpin eddy vision of boundary layer structure proposed by Theodorsen. That the structure of the flow is revealed more effectively by the discriminant rather than by the vorticity is important and adds support to recent observations of the discriminant in a channel flow simulation. Of particular importance is the fact that the procedure does not require the use of an arbitrary threshold in the discriminant. Further analysis using computer flow visualization shows a high degree of spatial correlation between regions of positive discriminant, extreme negative pressure fluctuations and large instantaneous values of Reynolds shear stress.

Kline, Stephen

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong

A preliminary assessment of the accuracy of selected meteorological parameters determined from Nimbus 6 satellite profile data

Published rms errors in rawinsonde data and discrepancies between satellite and rawinsonde profile data for temperature, dewpoint temperature, mixing ratio, and wind speed. Satellite rms errors were found to be 2 to 3 times as large as those for rawinsonde data. Gradients of the preceding parameters were computed for both rawinsonde and satellite data and compared with means and near extreme values computed from the AVE 2 and AVE 4 experiments. In all cases, it was found that satellite data can be used to determine with relatively good accuracy the near extreme gradients but not those whose value does not exceed the average. Synoptic charts were prepared to show that patterns of temperature could be determined with relatively good accuracy, while those of dew point were not as good as those for temperature. Winds represented by cloud motion vectors (satellite winds) were compared with rawinsonde winds, and it was found that large gaps exist in satellite values for a given pressure level and that errors in the satellite determined concluded that satellite profile data are very useful in synoptic analysis, particularly in data sparse regions as well as regions where near extreme gradients exist in the measured parameters.

Scoggins, J. R.