Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 973 records · Page 54

Landsat test of diffuse reflectance models for aquatic suspended solids measurement

Landsat radiance data were used to test mathematical models relating diffuse reflectance to aquatic suspended solids concentration. Digital CCT data for Landsat passes over the Bay of Fundy, Nova Scotia were analyzed on a General Electric Co. Image 100 multispectral analysis system. Three data sets were studied separately and together in all combinations with and without solar angle correction. Statistical analysis and chromaticity analysis show that a nonlinear relationship between Landsat radiance and suspended solids concentration is better at curve-fitting than a linear relationship. In particular, the quasi-single-scattering diffuse reflectance model developed by Gordon and coworkers is corroborated. The Gordon model applied to 33 points of MSS 5 data combined from three dates produced r = 0.98.

Munday, J. C., Jr.↗

Methodology for Determining Limit Torques for Threaded Fasteners

In aerospace design, where minimizing weight is always a priority, achieving the full capacity from fasteners is essential. To do so, the initial bolt preload must be maximized. The benefits of high preload are well documented and include improved fatigue resistance, a stiffer joint, and resistance to loosening. But many factors like elastic interactions and embedment tend to lower the initial preload placed on the bolt. These factors provide additional motivation to maximize the initial preload. But, to maximize bolt preload, you must determine what torque to apply. Determining this torque is greatly complicated by the large preload scatter generally seen with torque control. This paper presents a detailed methodology for generating limit torques for threaded fasteners. This methodology accounts for the large scatter in preload found with torque control, and therefore, addresses the statistical nature of the problem. It also addresses prevailing torque, a feature common in aerospace fasteners. Although prevailing torque provides a desired locking feature, it can also increase preload scatter. In addition, it can limit the amount of preload that can be generated due to the torsion it creates in the bolt. This paper discusses the complications of prevailing torque and how best to handle it. A wide range of torque-tension bolt testing was conducted in support of this research. The results from this research will benefit the design engineer as well as analyst involved in the design of bolted joints, leading to better, more optimized structural designs.

Hissam, Andy↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

A Survey of Logic Formalisms to Support Mishap Analysis

Mishap investigations provide important information about adverse events and near miss incidents. They are intended to help avoid any recurrence of previous failures. Over time, they can also yield statistical information about incident frequencies that helps to detect patterns of failure and can validate risk assessments. However, the increasing complexity of many safety critical systems is posing new challenges for mishap analysis. Similarly, the recognition that many failures have complex, systemic causes has helped to widen the scope of many mishap investigations. These two factors have combined to pose new challenges for the analysis of adverse events. A new generation of formal and semi-formal techniques have been proposed to help investigators address these problems. We introduce the term mishap logics to collectively describe these notations that might be applied to support the analysis of mishaps. The proponents of these notations have argued that they can be used to formally prove that certain events created the necessary and sufficient causes for a mishap to occur. These proofs can be used to reduce the bias that is often perceived to effect the interpretation of adverse events. Others have argued that one cannot use logic formalisms to prove causes in the same way that one might prove propositions or theorems. Such mechanisms cannot accurately capture the wealth of inductive, deductive and statistical forms of inference that investigators must use in their analysis of adverse events. This paper provides an overview of these mishap logics. It also identifies several additional classes of logic that might also be used to support mishap analysis.

Johnson, Chris↗

Far field effects of weak spatially periodic inhomogeneities.

Small variations in the properties of an acoustic medium produce cumulative effects that significantly alter wave forms propagating over long distances. To investigate some aspects of these effects, without statistical models, the one-dimensional wave equation with a weak spatially periodic variation in the sound speed is treated. An approximation is developed valid for such weak inhomogeneities over large distances based upon the approximate treatment of a linear differential equation with a periodic coefficient. It is applied to the propagation of an initial step wave and identifies the dispersive nature of the medium which accounts for the cumulative effects.

Friedman, M. B.↗

Development of a satellite SAR image spectra and altimeter wave height data assimilation system for ERS-1

The applicability of ERS-1 wind and wave data for wave models was studied using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR data. A series of global wave hindcasts is made for the surface stress and surface wind fields by assimilation of scatterometer data for the full 96-day SEASAT and also for two wind field analyses for shorter periods by assimilation with the higher resolution ECMWF T63 model and by subjective analysis methods. It is found that wave models respond very sensitively to inconsistencies in wind field analyses and therefore provide a valuable data validation tool. Comparisons between SEASAT SAR image spectra and theoretical SAR spectra derived from the hindcast wave spectra by Monte Carlo simulations yield good overall agreement for 32 cases representing a wide variety of wave conditions. It is concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra with confidence for wave studies if supported by realistic wave models and theoretical computations of the strongly nonlinear mapping of the wave spectrum into the SAR image spectrum. A closed nonlinear integral expression for this spectral mapping relation is derived which avoids the inherent statistical errors of Monte Carlo computations and may prove to be more efficient numerically.

Hasselmann, Klaus↗

Voids in the Ly-alpha forest

The general issue of using the Ly-alpha absorbers in quasar spectra as tracers of the large-scale matter distribution is discussed, emphasizing complications due to evolutionary effects and line blending. Then, a novel method of statistical analysis of the line data appropriate for studying large-scale structure is described. This method circumvents the problem of line blending by truncating the interval distribution at the low end. It also explicitly takes z-dependent effects into account. As an example, the method is applied to an 18-quasar spectrum sample of 417 lines.

Duncan, Robert C.↗

Large-eddy simulation of a plane wake

Previously the theoretical development leading to the dynamic localization model (DLM) for large-eddy simulation (LES) was presented. The method has been successfully applied to isotropic turbulence, channel flow, and the flow over a backward-facing step. Here we apply the model to the computation of the temporally developing place wake. The two main objectives of this project are: (1) Use the model to perform an LES of a time developing plane wake and compare the results with direction numerical simulation (DNS) data to see if important statistical measures can be readily predicted, and to provide a relative evaluation of the several versions of the model in terms of predictive capability and cost; and (2) If the tests in (1) show that the model generates reliable predictions, then use the LES to study various aspects of the physics of turbulent wakes and mixing layers.

Ghosal, Sandip↗

Extended Statistical Short-Range Guidance for Peak Wind Speed Analyses at the Shuttle Landing Facility: Phase II Results

This report describes the results from Phase II of the AMU's Short-Range Statistical Forecasting task for peak winds at the Shuttle Landing Facility (SLF). The peak wind speeds are an important forecast element for the Space Shuttle and Expendable Launch Vehicle programs. The 45th Weather Squadron and the Spaceflight Meteorology Group indicate that peak winds are challenging to forecast. The Applied Meteorology Unit was tasked to develop tools that aid in short-range forecasts of peak winds at tower sites of operational interest. A seven year record of wind tower data was used in the analysis. Hourly and directional climatologies by tower and month were developed to determine the seasonal behavior of the average and peak winds. Probability density functions (PDF) of peak wind speed were calculated to determine the distribution of peak speed with average speed. These provide forecasters with a means of determining the probability of meeting or exceeding a certain peak wind given an observed or forecast average speed. A PC-based Graphical User Interface (GUI) tool was created to display the data quickly.

Lambert, Winifred C.↗

Higher-Order Spectral Analysis of a Nonlinear Pitch and Plunge Apparatus

Simulated aeroelastic responses of a nonlinear pitch and plunge apparatus are analyzed using various statistical signal processing techniques including higher-order spectral methods. A MATLAB version of the Nonlinear Aeroelastic Testbed Apparatus (NATA) at the Texas A&M University is used to generate various aeroelastic response data including limit cycle oscillations (LCO). Traditional and higher-order spectral (HOS) methods are applied to the simulated aeroelastic responses. Higher-order spectral methods are used to identify critical signatures that indicate the transition from linear to nonlinear (LCO) aeroelastic behavior.

Silva, Walter A.↗

Can Remotely Sensed Snow Disappearance Explain Seasonal Water Supply?

Understanding the relationship between remotely sensed snow disappearance and seasonal water supply may become vital in coming years to supplement limited ground based, in situ measurements of snow in a changing climate. For the period 2001–2019, we investigated the relationship between satellite derived Day of Snow Disappearance (DSD)—the date at which snow has completely disappeared—and the seasonal water supply, i.e., the April—July total streamflow volume, for 15 snow dominated basins across the western U.S. A Monte Carlo framework was applied, using linear regression models to evaluate the predictive skill—defined here as a model’s ability to accurately predict seasonal flow volumes—of varied predictors, including DSD and in situ snow water equivalent (SWE), across a range of spring forecast dates. In all basins there is a statistically significant relationship between mean DSD and seasonal water supply (p ≤ 0.05), with mean DSD explaining roughly half of the variance. Satellite-based model skill improves later in the forecast season, surpassing the skill of in-situ-based (SWE) models in skill in 10 of the 15 basins by the latest forecast date. We found little to no correlation between model error and basin characteristics such as elevation and the ratio of snow water equivalent to total precipitation. Despite a relatively short data record, this exploratory analysis shows promise for improving seasonal water supply prediction, in particular for snow dominated basins lacking in situ observations.

snow remote sensing↗

System and Method for Finite Element Simulation of Helicopter Turbulence

The present invention provides a turbulence model that has been developed for blade-element helicopter simulation. This model uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades. for a total of twenty blade-element stations. The simulator system includes a software implementation of flight dynamics that adheres to the guidelines for turbulence set forth in military specifications. One of the features of the present simulator system is that it applies simulated turbulence to the rotor blades of the helicopter, rather than to its center of gravity. The simulator system accurately models the rotor penetration into a gust field. It includes time correlation between the front and rear of the main rotor, as well as between the side forces felt at the center of gravity and at the tail rotor. It also includes features for added realism, such as patchy turbulence and vertical gusts in to which the rotor disc penetrates. These features are realized by a unique real time implementation of the turbulence filters. The new simulator system uses two arrays one on either side of the main rotor to record the turbulence field and to produce time-correlation from the front to the rear of the rotor disc. The use of Gaussian Interpolation between the two arrays maintains the statistical properties of the turbulence across the rotor disc. The present simulator system and method may be used in future and existing real-time helicopter simulations with minimal increase in computational workload.

McFarland, R. E.↗

Flare Statistics for Young Stars from a Convolutional Neural Network Analysis of TESS Data

All-sky photometric time-series missions have allowed for the monitoring of thousands of young (t(age) < 800 Myr) stars in order to understand the evolution of stellar activity. Here, we developed a convolutional neural network (CNN), stella, specifically trained to find flares in Transiting Exoplanet Survey Satellite (TESS) short-cadence data. We applied the network to 3200 young stars in order to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars t(age) > 50 Myr across all temperatures T(eff) ≥ 4000 K, while stars from 2300 ≤ T(eff) < 4000 K show no evolution across 800 Myr. Stars of T(eff) ≤ 4000 K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4%–7% more atmosphere over the first 1 Gyr. stella is an open-source Python toolkit hosted on GitHub and PyPI.

Adina D. Feinstein↗

An Automated Algorithm to Screen Massive Training Samples for a Global Impervious Surface Classification

An algorithm is developed to automatically screen the outliers from massive training samples for Global Land Survey - Imperviousness Mapping Project (GLS-IMP). GLS-IMP is to produce a global 30 m spatial resolution impervious cover data set for years 2000 and 2010 based on the Landsat Global Land Survey (GLS) data set. This unprecedented high resolution impervious cover data set is not only significant to the urbanization studies but also desired by the global carbon, hydrology, and energy balance researches. A supervised classification method, regression tree, is applied in this project. A set of accurate training samples is the key to the supervised classifications. Here we developed the global scale training samples from 1 m or so resolution fine resolution satellite data (Quickbird and Worldview2), and then aggregate the fine resolution impervious cover map to 30 m resolution. In order to improve the classification accuracy, the training samples should be screened before used to train the regression tree. It is impossible to manually screen 30 m resolution training samples collected globally. For example, in Europe only, there are 174 training sites. The size of the sites ranges from 4.5 km by 4.5 km to 8.1 km by 3.6 km. The amount training samples are over six millions. Therefore, we develop this automated statistic based algorithm to screen the training samples in two levels: site and scene level. At the site level, all the training samples are divided to 10 groups according to the percentage of the impervious surface within a sample pixel. The samples following in each 10% forms one group. For each group, both univariate and multivariate outliers are detected and removed. Then the screen process escalates to the scene level. A similar screen process but with a looser threshold is applied on the scene level considering the possible variance due to the site difference. We do not perform the screen process across the scenes because the scenes might vary due to the phenology, solar-view geometry, and atmospheric condition etc. factors but not actual landcover difference. Finally, we will compare the classification results from screened and unscreened training samples to assess the improvement achieved by cleaning up the training samples. Keywords:

Tan, Bin↗

Delineation of soil temperature regimes from HCMM data

Evaluation of LANDSAT and Heat Capacity Mapping Mission (HCMM) data as input into National Cooperative Soil Survey is discussed. Signature classification techniques were applied to 13 May 76 LANDSAT data. LANDSAT data was overlaid with HCMM data, revealing registration problems caused by a shortage of control points in LANDSAT data, and the WARP program developed to improve registration accuracy. Initial images for control point selection were produced using digital terrain elevation data. Statistical procedures for evaluating data classification and to describe spatial distribution of surface temperature and its correlation with soil surface conditions were investigated.

Day, R. L.↗

A second-order Budkyo-type parameterization of landsurface hydrology

A simple, second order parameterization of the water fluxes at a land surface for use as the appropriate boundary condition in general circulation models of the global atmosphere was developed. The derived parameterization incorporates the high nonlinearities in the relationship between the near surface soil moisture and the evaporation, runoff and percolation fluxes. Based on the one dimensional statistical dynamic derivation of the annual water balance, it makes the transition to short term prediction of the moisture fluxes, through a Taylor expansion around the average annual soil moisture. A comparison of the suggested parameterization is made with other existing techniques and available measurements. A thermodynamic coupling is applied in order to obtain estimations of the surface ground temperature.

Andreou, S. A.↗

Relationships between a potential field-source surface model of the coronal magnetic field and properties of the solar wind at 1 AU

A comparison is made between the properties at the source surface of a potential field model of the solar corona and solar wind properties observed with ISEE 3 at 1 AU and extrapolated back to the sun. This is done for three consecutive Carrington rotations in the summer of 1979. The already known properties of a velocity minimum at the interplanetary current sheet and a tendency for the average velocity to increase with longitudinal distance from the current sheet are recovered. A better correlation between source surface properties and velocity is demonstrated by using the field strength on the source surface rather than longitudinal distance from the current sheet. Applying a simple algorithm to distinguish between transient and corotating interplanetary variations substantially increases the correlation. However, given a reasonable estimate of the number of degrees of freedom for the sample, only the correlation between source surface field strength and flow speed after the transients have been eliminated is statistically significant at the 1 percent level.

Suess, S. T.↗

Origin of the Greenland Ice Sheet

A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.

North, G. R.↗