Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical Testing Approach”

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 379 records · Page 21

Dynamic Modeling using Output-Error Parameter Estimation based on Frequency Responses Estimated with Multisine Inputs

A method is developed for estimating model parameters, such as nondimensional stability and control derivatives, by fitting transfer function or state-space models to empirical frequency response data using the output-error approach. The frequency response data were computed using Fourier transforms of measured input and output data. The control surfaces were excited with periodic multisine inputs which facilitated time-ecient estimation of multiple-input multiple-output frequency responses. The method was applied to lateral data from a nonlinear flight dynamics simulation of the F-16 aircraft, and to longitudinal data from multiple repeated flight test maneuvers of the NASA T-2 subscale aircraft. Results using simulation data showed the frequency response method compared well to other standard methods for parameter estimation. In addition to including all the available inputs, outputs, and harmonic frequencies in the estimation, relatively small subsets of the measured data could also be used to focus on identifying specific parts of the model. Results from flight test data showed that parameter estimates and uncertainties determined from repeated maneuvers were accurate and in statistical agreement with each other.

Grauer, Jared A.↗

Wind Tunnel Test of A Scale Model of A Venus Probe to Determine Aeroacoustics Environment

The Zephyr probe for NASA’s DAVINCI project is being designed to fly to Venus and collect measurements to characterize the Venusian atmosphere while it descends towards the surface. Despite a slow descent speed as the probe approaches the surface, the very dense atmosphere means the vibro-acoustic environment that the probe must withstand could be very harsh. A wind tunnel test campaign was conducted to characterize the surface-pressure fluctuations (acoustics) on the external surface of the probe during the descent. The goal was to provide forcing functions for vibro-acoustic analysis of all instruments inside the probe. A 25% scale-model of the Zephyr probe was manufactured, instrumented with 18 microphones, and tested in two wind tunnels located in the Fluid Mechanics Lab at NASA Ames Research Center at velocities ranging from 13 m/s to 48 m/s. Trip strips were added to the model to ensure a turbulent boundary layer on the model. Microphones were located appropriately to provide spectra of pressure fluctuations in regions with different flow characteristics as well as for calculating various two-point statistics. The acoustic spectra measured in the wind tunnel tests were scaled up to flight-scale using the estimated flight trajectory information. Besides fluctuating surface-pressure measurements, hot-film velocity measurements as well as smoke-laser flow visualization were conducted to better understand the flow around the probe. Lastly, several different geometry configurations were tested to see whether any reduction in acoustic levels could be achieved, but none of the configurations that were tested offered significant improvements. The flow visualization showed 3 primary flow regimes on the model: an attached turbulent boundary layer on the front of the model, a region of separated flow upstream of the flared lip, and a large, separated wake at the rear of the model. The acoustic spectra generally show three different shapes corresponding to these flow regions, with the highest levels seen at the rear of the model, and directly in front of the drag plate. The scaled overall levels of pressure fluctuations were found to be high: in the range of 130dB to 155dB.

planetary probe↗

Automated Classification of Thermal Infrared Spectra Using Self-organizing Maps

Existing and planned space missions to a variety of planetary and satellite surfaces produce an ever increasing volume of spectral data. Understanding the scientific informational content in this large data volume is a daunting task. Fortunately various statistical approaches are available to assess such data sets. Here we discuss an automated classification scheme based on Kohonen Self-organizing maps (SOM) we have developed. The SUM process produces an output layer were spectra having similar properties lie in close proximity to each other. One major effort is partitioning this output layer into appropriate regions. This is prefonned by defining dosed regions based upon the strength of the boundaries between adjacent cells in the SOM output layer. We use the Davies-Bouldin index as a measure of the inter-class similarities and intra-class dissimilarities that determines the optimum partition of the output layer, and hence number of SOM clusters. This allows us to identify the natural number of clusters formed from the spectral data. Mineral spectral libraries prepared at Arizona State University (ASU) and John Hopkins University (JHU) are used to test and evaluate the classification scheme. We label the library sample spectra in a hierarchical scheme with class, subclass, and mineral group names. We use a portion of the spectra to train the SOM, i.e. produce the output layer, while the remaining spectra are used to test the SOM. The test spectra are presented to the SOM output layer and assigned membership to the appropriate cluster. We then evaluate these assignments to assess the scientific meaning and accuracy of the derived SOM classes as they relate to the labels. We demonstrate that unsupervised classification by SOMs can be a useful component in autonomous systems designed to identify mineral species from reflectance and emissivity spectra in the therrnal IR.

Roush, Ted L.↗

Statistical Signal Analysis for Systems with Interferenced Inputs

A new approach is introduced, based on statistical signal analysis, which overcomes the error due to input signal interference. The model analyzed is given. The input signals u sub 1 (t) and u sub 2 (t) are assumed to be unknown. The measurable signals x sub 1 (t) and x sub 2 (t) are interferened according to the frequency response functions, H sub 12 (f) and H sub 21 (f). The goal of the analysis was to evaluate the power output due to each input, u sub 1 (t) and u sub 2 (t), for the case where both are applied to the same time. In addition, all frequency response functions are calculated. The interferenced system is described by a set of five equations with six unknown functions. An IBM XT Personal Computer, which was interfaced with the FFT, was used to solve the set of equations. The software was tested on an electrical two-input, one-output system. The results were excellent. The research presented includes the analysis of the acoustic radiation from a rectangular plate with two force inputs and the sound pressure as an output signal.

Bai, R. M.↗

On percolation as a cosmological test

Difficulties in the use of percolation as a complementary statistic for the galaxy clustering pattern are pointed out by studying simple toy models and dynamical N-body models that represent the competing clustering scenarios. The percolation properties are found not to be very sensitive to the presence of pancakes and strings once they are clumpy, and hence they do not distinguish properly between models that are very different. In the case of very smooth pancakes, the ability to percolate depends on sampling parameters, such as the mean number density and the volume, in a way which is unknown a priori because it depends on the same properties that the test ought to measure. This problem could, in principle, be eased by using volume-limited samples of high mean number density (an order of magnitude denser than the CfA redshift survey volume limited at 4000 km/s) and by comparing to models of identical number density and volume. An alternative approach, based on the sampling effects themselves, may provide a qualitative test for pancakes in samples of lower densities.

Dekel, A.↗

Using Two-Stream Theory to Capture Fluctuations of Satellite-Perceived TOA SW Radiances Reflected from Clouds Over Ocean

Shortwave (SW) fluxes estimated from broadband radiometry rely on empirically gathered and hemispherically resolved fields of outgoing top-of-atmosphere (TOA) radiances. This study aims to provide more accurate and precise fields of TOA SW radiances reflected from clouds over ocean by introducing a novel semiphysical model predicting radiances per narrow sun-observer geometry. This model was statistically trained using CERES-measured radiances paired with MODIS-retrieved cloud parameters as well as reanalysis-based geophysical parameters. By using radiative transfer approximations as a framework to ingest the above parameters, the new approach incorporates cloud-top effective radius and above-cloud water vapor in addition to traditionally used cloud optical depth, cloud fraction, cloud phase, and surface wind speed. A two-stream cloud albedo – serving to statistically incorporate cloud optical thickness and cloud-top effective radius – and Cox–Munk ocean reflectance were used to describe an albedo over each CERES footprint. Effective-radius-dependent asymmetry parameters were obtained empirically and separately for each viewing-illumination geometry. A simple equation of radiative transfer, with this albedo and attenuating above-cloud water vapor as inputs, was used in its log-linear form to allow for statistical optimization. We identified the two-stream functional form that minimized radiance residuals calculated against CERES observations and outperformed the state-of-the-art approach for most observer geometries outside the sun-glint and solar zenith angles between 20 and 70∘, reducing the median SD of radiance residuals per solar geometry by up to 13.2 % for liquid clouds, 1.9 % for ice clouds, and 35.8 % for footprints containing both cloud phases. Geometries affected by sun glint (constituting between 10 % and 1 % of the discretized upward hemisphere for solar zenith angles of 20 and 70∘, respectively), however, often showed weaker performance when handled with the new approach and had increased residuals by as much as 60 % compared to the state-of-the-art approach. Overall, uncertainties were reduced for liquid-phase and mixed-phase footprints by 5.76 % and 10.81 %, respectively, while uncertainties for ice-phase footprints increased by 0.34 %. Tested for a variety of scenes, we further demonstrated the plausibility of scene-wise predicted radiance fields. This new approach may prove useful when employed in angular distribution models and may result in improved flux estimates, in particular dealing with clouds characterized by small or large droplet/crystal sizes.

Two- stream theory↗

Jimsphere wind and turbulence exceedance statistic

Exceedance statistics of winds and gusts observed over Cape Kennedy with Jimsphere balloon sensors are described. Gust profiles containing positive and negative departures, from smoothed profiles, in the wavelength ranges 100-2500, 100-1900, 100-860, and 100-460 meters were computed from 1578 profiles with four 41 weight digital high pass filters. Extreme values of the square root of gust speed are normally distributed. Monthly and annual exceedance probability distributions of normalized rms gust speeds in three altitude bands (2-7, 6-11, and 9-14 km) are log-normal. The rms gust speeds are largest in the 100-2500 wavelength band between 9 and 14 km in late winter and early spring. A study of monthly and annual exceedance probabilities and the number of occurrences per kilometer of level crossings with positive slope indicates significant variability with season, altitude, and filter configuration. A decile sampling scheme is tested and an optimum approach is suggested for drawing a relatively small random sample that represents the characteristic extreme wind speeds and shears of a large parent population of Jimsphere wind profiles.

Adelfang, S. I.↗

CANDELS: The Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, from z approx. 8 - 1.5. It will image > 250,000 distant galaxies using three separate cameras on the Hubble Space Tele8cope, from the mid-UV to near-IR, and will find and measure Type Ia supernovae beyond z > 1.5 to test their accuracy as standard candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10(exp 9) solar mass to z approx. 2, reaching the knee of the UV luminosity function of galaxies to z approx. 8. The survey covers approximately 800 square arc minutes and is divided into two parts. The CANDELS/Deep survey (5(sigma) point-source limit H =27.7mag) covers approx. 125 square arcminutes within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 50(sigma) point-source limit of H ? or approx. = 27.0 mag. Together with the Hubble Ultradeep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are non-proprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design.

Grogin, Norman A.↗

Associated with aerospace vehicles development of methodologies for the estimation of thermal properties

Thermal stress analyses are an important aspect in the development of aerospace vehicles at NASA-LaRC. These analyses require knowledge of the temperature distributions within the vehicle structures which consequently necessitates the need for accurate thermal property data. The overall goal of this ongoing research effort is to develop methodologies for the estimation of the thermal property data needed to describe the temperature responses of these complex structures. The research strategy undertaken utilizes a building block approach. The idea here is to first focus on the development of property estimation methodologies for relatively simple conditions, such as isotropic materials at constant temperatures, and then systematically modify the technique for the analysis of more and more complex systems, such as anisotropic multi-component systems. The estimation methodology utilized is a statistically based method which incorporates experimental data and a mathematical model of the system. Several aspects of this overall research effort were investigated during the time of the ASEE summer program. One important aspect involved the calibration of the estimation procedure for the estimation of the thermal properties through the thickness of a standard material. Transient experiments were conducted using a Pyrex standard at various temperatures, and then the thermal properties (thermal conductivity and volumetric heat capacity) were estimated at each temperature. Confidence regions for the estimated values were also determined. These results were then compared to documented values. Another set of experimental tests were conducted on carbon composite samples at different temperatures. Again, the thermal properties were estimated for each temperature, and the results were compared with values obtained using another technique. In both sets of experiments, a 10-15 percent off-set between the estimated values and the previously determined values was found. Another effort was related to the development of the experimental techniques. Initial experiments required a resistance heater placed between two samples. The design was modified such that the heater was placed on the surface of only one sample, as would be necessary in the analysis of built up structures. Experiments using the modified technique were conducted on the composite sample used previously at different temperatures. The results were within 5 percent of those found using two samples. Finally, an initial heat transfer analysis, including conduction, convection and radiation components, was completed on a titanium sandwich structural sample. Experiments utilizing this sample are currently being designed and will be used to first estimate the material's effective thermal conductivity and later to determine the properties associated with each individual heat transfer component.

Scott, Elaine P.↗

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a “Virtual Constellation” was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating “stable” calibration to within 5%—the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)–were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Student’s T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process.

Tuli, Fatima Tuz Zafrin↗

A Regression Method for Comparing Launch Zone Ranges Derived From Different Sources

In modern fighter aircraft air-to-air missile launch acceptable ranges (LARs) are provided to the pilot for his launch decision. These LARs are calculated in the avionic software, using very simple algebraic equations and algorithms, that attempt to model complex missile behavior. Frequently these algorithms are revised to reflect inaccuracies discovered, when comparisons are made with missile six degree-of-freedom (6-DOF) simulations. The process of modifying this software has become expensive. In the past, these simple LAR generators were modified based only on subjective judgmental approach, with little regard to the missile range design limits. With reduced budgets a more rational methodology should be taken to assess whether or not simple fighter LAR generators need to be modified. Any LAR generator can produce large amounts of range (minimum and maximum) data over the board engagement limits of modem fighter and air-to-air missiles. This infinite population lends itself to some statistical methodology to access simple LAR accuracy and the need for software modification. The methodology should produce statistical parameters for use in the decision making process of modifying a simple fighter LAR generator. This paper presents a regression method, based on a missile range design limit philosophy. Acceptance criteria from the regression method are given that provide simple statistical parameters to aid in this decision making process.

Flight Testing↗

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗

Technology Benefit Estimator (T/BEST): User's Manual

The Technology Benefit Estimator (T/BEST) system is a formal method to assess advanced technologies and quantify the benefit contributions for prioritization. T/BEST may be used to provide guidelines to identify and prioritize high payoff research areas, help manage research and limited resources, show the link between advanced concepts and the bottom line, i.e., accrued benefit and value, and to communicate credibly the benefits of research. The T/BEST software computer program is specifically designed to estimating benefits, and benefit sensitivities, of introducing new technologies into existing propulsion systems. Key engine cycle, structural, fluid, mission and cost analysis modules are used to provide a framework for interfacing with advanced technologies. An open-ended, modular approach is used to allow for modification and addition of both key and advanced technology modules. T/BEST has a hierarchical framework that yields varying levels of benefit estimation accuracy that are dependent on the degree of input detail available. This hierarchical feature permits rapid estimation of technology benefits even when the technology is at the conceptual stage. As knowledge of the technology details increases the accuracy of the benefit analysis increases. Included in T/BEST's framework are correlations developed from a statistical data base that is relied upon if there is insufficient information given in a particular area, e.g., fuel capacity or aircraft landing weight. Statistical predictions are not required if these data are specified in the mission requirements. The engine cycle, structural fluid, cost, noise, and emissions analyses interact with the default or user material and component libraries to yield estimates of specific global benefits: range, speed, thrust, capacity, component life, noise, emissions, specific fuel consumption, component and engine weights, pre-certification test, mission performance engine cost, direct operating cost, life cycle cost, manufacturing cost, development cost, risk, and development time. Currently, T/BEST operates on stand-alone or networked workstations, and uses a UNIX shell or script to control the operation of interfaced FORTRAN based analyses. T/BEST's interface structure works equally well with non-FORTRAN or mixed software analysis. This interface structure is designed to maintain the integrity of the expert's analyses by interfacing with expert's existing input and output files. Parameter input and output data (e.g., number of blades, hub diameters, etc.) are passed via T/BEST's neutral file, while copious data (e.g., finite element models, profiles, etc.) are passed via file pointers that point to the expert's analyses output files. In order to make the communications between the T/BEST's neutral file and attached analyses codes simple, only two software commands, PUT and GET, are required. This simplicity permits easy access to all input and output variables contained within the neutral file. Both public domain and proprietary analyses codes may be attached with a minimal amount of effort, while maintaining full data and analysis integrity, and security. T/BESt's sotware framework, status, beginner-to-expert operation, interface architecture, analysis module addition, and key analysis modules are discussed. Representative examples of T/BEST benefit analyses are shown.

Generazio, Edward R.↗

Anomaly Detection for the Roman Space Telescope Wide Field Instrument’s Science Data Processing Pipeline

The Roman Space Telescope (RST) Wide Field Instrument (WFI) will be utilizing a preliminary Science Data Processing (SDP) pipeline during its Integration and Test, and to some extent during Operations, to track basic statistics and identify known features such as cosmic rays, snowballs as well as possible anomalies in raw detector data. In our detectors, these anomalies appear as jumps in the ramp of a readout and are classified as cosmic rays if they appear as a streak or snowballs if they’re more circular. The WFI employs an array of 18 H4RG-10 detectors that collect image samples. Each set of raw frames within a non-destructive exposure is packaged by the SDP pipeline into image cubes for each detector. Each cube is a time series of 4096 × 4096 accumulating pixel frames. The preliminary analysis pipeline is used to locate anomalies in these time-series accumulation frames and identify the type of anomaly, either natural phenomena or detector characteristic. To compare different methods, we’ve implemented both heuristic-based and data-driven methods to identify anomalies. For the heuristic-based approach, we identify snowballs and cosmic rays by the size and shape of outlier pixel clusters between consecutive frames. For data driven methods, we evaluated a Convolutional Neural Network (CNN) model, and more traditional methods like Principal Component Analysis (PCA). CNN is a supervised learning/classification method. Thus, we used a labeled dataset of anomalies to perform segmentation of the image and identify anomalies. We used previously identified cosmic rays and snowballs to measure the accuracy and efficiency of the mentioned approaches. In evaluating these methods, we aim to pick the best fit for the SDP pipeline’s anomaly detection in terms of both performance and runtime.

Paul Horton↗

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

Kepler Data Validation and Follow-up Programs

The approach that the Kepler Mission uses to remove false positive events and to validate the discoveries consists of two parts; data validation (DPI) and follow up observations program (FOP). DV consists of several methods of examining the data from the spacecraft observations. First, to rule out statistical fluctuations in the data, accept only signals that show 3 or more transits and that have a total signal-to-noise ratio that exceeds 7 sigma. Second. to identify small stellar companions to the target star, we check for secondary eclipses and determine if the transit characteristics are appropriate for a planetary companion. Third, check for background binaries that are in the target aperture. Here we measure the movement of the image centroid before, during, and after the transit. If the target is producing the signal, a dimming wi11 move the image centroid in a known direction and magnitude. If the signal comes from a nearby star, the amplitude and direction of the motion wi11 be different, This test is expected to rule out the hundreds of binary signals expected from background stars. The precision of the measurement depends on the stellar fluxes and positions but can be better than 0.01 pixel; i.e., 0.04". Those candidates that pass these tests are examined using ground-based telescopes and radial velocity spectrometers. First medium precision RV is used to rule out any remaining stellar companions. Then high spatial resolution imaging is used to check for nearby stars that are in the aperture- (The Kepler apertures depend on magnitude but are of order 36 sq are sec in area.) If no stars are present that quid generate the observed signal, then the candidate goes to a large telescope such as Keck, HET, or Wi1lilam Herschel for high precision observations to get the planet mass or an upper limit to it, if there are some stars in the aperture, then the photometric observations are employed to look for the transit by cane of the confounding stars. If none are seen, then the candidate is sent to a large telescope such as Keck

William, Borucki J.↗

Soil moisture estimation using GOES-VISSR infrared data - A case study with a simple statistical method

Five days of clear sky observations of Kansas and Nebraska are used to examine the statistical relationship between soil moisture and infrared surface temperature observations taken from a geosynchronous satellite. Linear regression is used to relate soil moisture to surface temperature and other variables that represent wind speed, vegetation cover, and low-level temperature advection. Results show good agreement between estimated and observed soil moisture features on each of the 5 days. The average coefficient of determination for five pseudoindependent tests in which the test day is held out of the regression is 0.71. It is shown that a depletion coefficient of 0.92, when used to compute antecedent precipitation index (API), produces the best correlation between API and soil moisture as inferred from GOES thermal infrared data. By averaging daily predicted values over the 5-day rain-free case study period, 92 percent of the variance of the morning surface temperature change is explained by a simple multiple linear regression with all independent variables, or, alternatively, 85 percent of the observed variance in API is explained. It is concluded that this approach can distinguish at least four classes of soil wetness, but the necessity for measurement of surface advection may limit its usefulness in remote areas.

Wetzel, Peter J.↗

Implementing Software Safety in the NASA Environment

Until recently, NASA did not consider allowing computers total control of flight systems. Human operators, via hardware, have constituted the ultimate safety control. In an attempt to reduce costs, NASA has come to rely more and more heavily on computers and software to control space missions. (For example. software is now planned to control most of the operational functions of the International Space Station.) Thus the need for systematic software safety programs has become crucial for mission success. Concurrent engineering principles dictate that safety should be designed into software up front, not tested into the software after the fact. 'Cost of Quality' studies have statistics and metrics to prove the value of building quality and safety into the development cycle. Unfortunately, most software engineers are not familiar with designing for safety, and most safety engineers are not software experts. Software written to specifications which have not been safety analyzed is a major source of computer related accidents. Safer software is achieved step by step throughout the system and software life cycle. It is a process that includes requirements definition, hazard analyses, formal software inspections, safety analyses, testing, and maintenance. The greatest emphasis is placed on clearly and completely defining system and software requirements, including safety and reliability requirements. Unfortunately, development and review of requirements are the weakest link in the process. While some of the more academic methods, e.g. mathematical models, may help bring about safer software, this paper proposes the use of currently approved software methodologies, and sound software and assurance practices to show how, to a large degree, safety can be designed into software from the start. NASA's approach today is to first conduct a preliminary system hazard analysis (PHA) during the concept and planning phase of a project. This determines the overall hazard potential of the system to be built. Shortly thereafter, as the system requirements are being defined, the second iteration of hazard analyses takes place, the systems hazard analysis (SHA). During the systems requirements phase, decisions are made as to what functions of the system will be the responsibility of software. This is the most critical time to affect the safety of the software. From this point, software safety analyses as well as software engineering practices are the main focus for assuring safe software. While many of the steps proposed in this paper seem like just sound engineering practices, they are the best technical and most cost effective means to assure safe software within a safe system.

Wetherholt, Martha S.↗