Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simple linear regression”

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.

70 records · Page 4

Transient Modeling of Hybrid Rocket Low Frequency Instabilities

A comprehensive dynamic model of a hybrid rocket has been developed in order to understand and predict the transient behavior including instabilities. A linearized version of the transient model predicted the low-frequency chamber pressure oscillations that are commonly observed in hybrids. The source of the instabilities is based on a complex coupling of thermal transients in the solid fuel, wall heat transfer blocking due to fuel regression rate and the transients in the boundary layer that forms on the fuel surface. The oscillation frequencies predicted by the linearized theory are in very good agreement with 43 motor test results obtained from the hybrid propulsion literature. The motor test results used in the comparison cover a very wide spectrum of parameters including: 1) four separate research and development programs, 2) three different oxidizers (LOX, GOX, N2O), 3) a wide range of motor dimensions (i.e. from 5 inch diameter to 72 inch diameter) and operating conditions and 4) several fuel formulations. A simple universal scaling formula for the frequency of the primary oscillation mode is suggested.

Karabeyoglu, M. Arif↗

Non-Contact Thrust Stand Calibration Method for Repetitively-Pulsed Electric Thrusters

A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.

Wong, Andrea R.↗

Evaluation of the Williams-type model for barley yields in North Dakota and Minnesota

The Williams-type yield model is based on multiple regression analysis of historial time series data at CRD level pooled to regional level (groups of similar CRDs). Basic variables considered in the analysis include USDA yield, monthly mean temperature, monthly precipitation, soil texture and topographic information, and variables derived from these. Technologic trend is represented by piecewise linear and/or quadratic functions of year. Indicators of yield reliability obtained from a ten-year bootstrap test (1970-1979) demonstrate that biases are small and performance based on root mean square appears to be acceptable for the intended AgRISTARS large area applications. The model is objective, adequate, timely, simple, and not costly. It consideres scientific knowledge on a broad scale but not in detail, and does not provide a good current measure of modeled yield reliability.

Barnett, T. L.↗

Core Body Temperature Predictions Using Metabolic Energy Expenditure and Heart Rate During Simulated Extravehicular Activity

Long duration spaceflight missions will require crew to become more autonomous in conducting extravehicular activities (EVA) without direct communication with Mission Control for biomedical support. To enable such autonomy, we are developing a Crew State and Risk Model (CSRM) as a collection of key physiology domains that drive EVA crew capabilities and workloads. One model component of CSRM is human thermal regulation. In this paper, customized development of a baseline model to predict core body temperature is presented using physiology inputs of heart rate and metabolic rate. The model development dataset included a baseline study where participants (n=6, equal male and female) performed a 5-hour EVA in a two-part session while wearing a hybrid space suit simulator (HS3). The first session included an end-to-end EVA traversing 1500 meters to a geology site and traversing back to a habitat conducting geology, payload relocation, and maintenance operations every 500 meters. The second session included standalone tasks of a 2000-meter traverse followed by geology tasks. Thermal measures of core body temperature, local skin temperature, liquid cooling garment temperature, heart rate, and metabolic rate were collected through the test duration. Multiple regression was used to build a linear equation to predict core body temperature from inputs of heart rate and metabolic rate. Heat storage was calculated via predicted core body temperature plus LCG and skin temperatures. The model was tested against a dataset from a pressurized suited test (n=6, equal male and female) in the NASA Active Response Gravity Offload System (ARGOS) conducting similar end-to-end EVA and standalone tasks. Predicted error of the model against the raw test cases was 0.2±0.15 °C. The baseline prediction of core body temperature using heart rates and metabolic rates allows for simple real-time tracking from data collected in-flight to monitor crew consumables and thermal flight limits during EVA.

Bradley Hoffmann↗

Vegetation spatial variability and its effect on vegetation indices

Landsat MSS data were used to simulate low resolution satellite data, such as NOAA AVHRR, to quantify the fractional vegetation cover within a pixel and relate the fractional cover to the normalized difference vegetation index (NDVI) and the simple ratio (SR). The MSS data were converted to radiances from which the NDVI and SR values for the simulated pixels were determined. Each simulated pixel was divided into clusters using an unsupervised classification program. Spatial and spectral analysis provided a means of combining clusters representing similar surface characteristics into vegetated and non-vegetated areas. Analysis showed an average error of 12.7 per cent in determining these areas. NDVI values less than 0.3 represented fractional vegetated areas of 5 per cent or less, while a value of 0.7 or higher represented fractional vegetated areas greater than 80 per cent. Regression analysis showed a strong linear relation between fractional vegetation area and the NDVI and SR values; correlation values were 0.89 and 0.95 respectively. The range of NDVI values calculated from the MSS data agrees well with field studies.

Ormsby, J. P.↗

Development of the Ames Global Hyperspectral Synthetic Dataset

This study develops the surface BRDF (bidirectional reflectance distribution function) product of the Ames Global Hyperspectral Synthetic Dataset (AGHSD), based on the corresponding MODIS products, to support the NASA Surface Biology and Geology mission development. A main challenge in deriving a hyperspectral dataset from the multi-band satellite products is how to identify a succinct yet robust algorithm that allow us to infer BRDF at unobserved wavelengths based on the few observed bands. Using the theories of radiative transfer in vegetation canopies, we arrive at a simple equation that accurately approximates hyperspectral surface BRDF as the weighted sum of components from the soil and the vegetation. Each of the components is modeled by the product of the spectrally-dependent optical properties of a surface element (the spectra of the soil surface reflectance, the leaf single albedo, or the canopy scattering coefficient) and a spectrally-independent bidirectional scattering function. The optical properties of the soil and the vegetation can be obtained from existing spectral libraries or model simulations. The bidirectional scattering functions are represented by the Ross-Thick-Li-Sparse BRDF model, where the linear coefficients are estimated with regression analysis from the multi-band MODIS data. We validate the algorithm with simulations by Monte Carlo Ray Tracing model experiments, and the results are highly consistent with the theoretic derivation. We apply the algorithm to generate the AGHSD BRDF product at 1km and 8-day resolutions for the year of 2019. The results are biogeochemically and physically coherent and consistent, and thus serve the goal to support the science and application development of the SBG community.

Hyperspectral↗

Design and Development of a Model to Simulate 0-G Treadmill Running Using the European Space Agency's Subject Loading System

Develop a model that simulates a human running in 0 G using the European Space Agency s (ESA) Subject Loading System (SLS). The model provides ground reaction forces (GRF) based on speed and pull-down forces (PDF). DESIGN The theoretical basis for the Running Model was based on a simple spring-mass model. The dynamic properties of the spring-mass model express theoretical vertical GRF (GRFv) and shear GRF in the posterior-anterior direction (GRFsh) during running gait. ADAMs VIEW software was used to build the model, which has a pelvis, thigh segment, shank segment, and a spring foot (see Figure 1).the model s movement simulates the joint kinematics of a human running at Earth gravity with the aim of generating GRF data. DEVELOPMENT & VERIFICATION ESA provided parabolic flight data of subjects running while using the SLS, for further characterization of the model s GRF. Peak GRF data were fit to a linear regression line dependent on PDF and speed. Interpolation and extrapolation of the regression equation provided a theoretical data matrix, which is used to drive the model s motion equations. Verification of the model was conducted by running the model at 4 different speeds, with each speed accounting for 3 different PDF. The model s GRF data fell within a 1-standard-deviation boundary derived from the empirical ESA data. CONCLUSION The Running Model aids in conducting various simulations (potential scenarios include a fatigued runner or a powerful runner generating high loads at a fast cadence) to determine limitations for the T2 vibration isolation system (VIS) aboard the International Space Station. This model can predict how running with the ESA SLS affects the T2 VIS and may be used for other exercise analyses in the future.

Caldwell, E. C.↗

Comparison of CEAS and Williams-type models for spring wheat yields in North Dakota and Minnesota

The CEAS and Williams-type yield models are both based on multiple regression analysis of historical time series data at CRD level. The CEAS model develops a separate relation for each CRD; the Williams-type model pools CRD data to regional level (groups of similar CRDs). Basic variables considered in the analyses are USDA yield, monthly mean temperature, monthly precipitation, and variables derived from these. The Williams-type model also used soil texture and topographic information. Technological trend is represented in both by piecewise linear functions of year. Indicators of yield reliability obtained from a ten-year bootstrap test of each model (1970-1979) demonstrate that the models are very similar in performance in all respects. Both models are about equally objective, adequate, timely, simple, and inexpensive. Both consider scientific knowledge on a broad scale but not in detail. Neither provides a good current measure of modeled yield reliability. The CEAS model is considered very slightly preferable for AgRISTARS applications.

Barnett, T. L.↗

Understanding software faults and their role in software reliability modeling

This study is a direct result of an on-going project to model the reliability of a large real-time control avionics system. In previous modeling efforts with this system, hardware reliability models were applied in modeling the reliability behavior of this system. In an attempt to enhance the performance of the adapted reliability models, certain software attributes were introduced in these models to control for differences between programs and also sequential executions of the same program. As the basic nature of the software attributes that affect software reliability become better understood in the modeling process, this information begins to have important implications on the software development process. A significant problem arises when raw attribute measures are to be used in statistical models as predictors, for example, of measures of software quality. This is because many of the metrics are highly correlated. Consider the two attributes: lines of code, LOC, and number of program statements, Stmts. In this case, it is quite obvious that a program with a high value of LOC probably will also have a relatively high value of Stmts. In the case of low level languages, such as assembly language programs, there might be a one-to-one relationship between the statement count and the lines of code. When there is a complete absence of linear relationship among the metrics, they are said to be orthogonal or uncorrelated. Usually the lack of orthogonality is not serious enough to affect a statistical analysis. However, for the purposes of some statistical analysis such as multiple regression, the software metrics are so strongly interrelated that the regression results may be ambiguous and possibly even misleading. Typically, it is difficult to estimate the unique effects of individual software metrics in the regression equation. The estimated values of the coefficients are very sensitive to slight changes in the data and to the addition or deletion of variables in the regression equation. Since most of the existing metrics have common elements and are linear combinations of these common elements, it seems reasonable to investigate the structure of the underlying common factors or components that make up the raw metrics. The technique we have chosen to use to explore this structure is a procedure called principal components analysis. Principal components analysis is a decomposition technique that may be used to detect and analyze collinearity in software metrics. When confronted with a large number of metrics measuring a single construct, it may be desirable to represent the set by some smaller number of variables that convey all, or most, of the information in the original set. Principal components are linear transformations of a set of random variables that summarize the information contained in the variables. The transformations are chosen so that the first component accounts for the maximal amount of variation of the measures of any possible linear transform; the second component accounts for the maximal amount of residual variation; and so on. The principal components are constructed so that they represent transformed scores on dimensions that are orthogonal. Through the use of principal components analysis, it is possible to have a set of highly related software attributes mapped into a small number of uncorrelated attribute domains. This definitively solves the problem of multi-collinearity in subsequent regression analysis. There are many software metrics in the literature, but principal component analysis reveals that there are few distinct sources of variation, i.e. dimensions, in this set of metrics. It would appear perfectly reasonable to characterize the measurable attributes of a program with a simple function of a small number of orthogonal metrics each of which represents a distinct software attribute domain.

Munson, John C.↗

A Note on the Relationship between Temperature and Water Vapor in Quasi-Equilibrium and Climate States

An ideal and simple formulation is successfully derived that well represents a quasi-linear relationship found between the domain-averaged water vapor, q (mm), and temperature, T (K), fields obtained from a series of quasi-equilibrium (long-term) simulations for the Tropics using the two-dimensional Goddard Cumulus Ensemble (GCE) model. Earlier model work showed that the forced maintenance of two different wind profiles in the Tropics leads to two different equilibrium states. Investigating this finding required investigation of the slope of the moisture-temperature relations, which turns out to be linear in the Tropics. The extra-tropical climate equilibriums become more complex, but insight on modeling sensitivity can be obtained by linear stepwise regression of the integrated temperature and humidity. A globally curvilinear moisture-temperature distribution, similar to the famous Clausius-Clapeyron curve (i.e., saturated water vapor pressure versus temperature), is then found in this study. Such a genuine finding clarifies that the dynamics are crucial to the climate (shown in the earlier work) but the thermodynamics adjust. The range of validity of this result is further examined herein. The GCE-modeled tropical domain-averaged q and T fields form a linearly-regressed "q-T" slope that genuinely resides within an ideal range of slopes obtained from the aforementioned formulation. A quantity (denoted as dC2/dC1) representing the derivative between the static energy densities due to temperature (C2) and water vapor (C1) for various quasi-equilibrium states can also be obtained. A dC2/dC1 value near unity obtained for the GCE-modeled tropical simulations implies that the static energy densities due to moisture and temperature only differ by a pure constant for various equilibrium states. An overall q-T relation also including extra-tropical regions is, however, found to have a curvilinear relationship. Accordingly, warm/moist regions favor change in water vapor faster than temperature, while cold/dry regions favor an increase in temperature quicker than water vapor.

Shie, C.-L.↗

The effect of clouds on the earth's radiation budget

The radiative fluxes from the Earth Radiation Budget Experiment (ERBE) and the cloud properties from the International Satellite Cloud Climatology Project (ISCCP) over Indonesia for the months of June and July of 1985 and 1986 were analyzed to determine the cloud sensitivity coefficients. The method involved a linear least squares regression between co-incident flux and cloud coverage measurements. The calculated slope is identified as the cloud sensitivity. It was found that the correlations between the total cloud fraction and radiation parameters were modest. However, correlations between cloud fraction and IR flux were improved by separating clouds by height. Likewise, correlations between the visible flux and cloud fractions were improved by distinguishing clouds based on optical depth. Calculating correlations between the net fluxes and either height or optical depth segregated cloud fractions were somewhat improved. When clouds were classified in terms of their height and optical depth, correlations among all the radiation components were improved. Mean cloud sensitivities based on the regression of radiative fluxes against height and optical depth separated cloud types are presented. Results are compared to a one-dimensional radiation model with a simple cloud parameterization scheme.

Ziskin, Daniel↗

Comparison of Likelihood Methods for Generalized Linear Mixed Models with Application to Quiet Supersonic Flights 2018 Data

Repeated measurement will be a feature of the survey data collected during the Quesst missionX-59 community response tests (CRT). Since each participant will report his or her categorical level of annoyance in response to multiple events, the responses from any single individual may be correlated with one another. Several models within the class of generalized linear mixed models (GLMM) are pertinent to the analysis of correlated categorical outcomes; the random intercept logistic regression model is one example. Both Bayesian and frequentist methods for fitting these models are available, with frequentist methods relying on some form of approximation (of either an integral or the integrand) that appears in the marginal likelihood function. Given several anticipated similarities of the X-59 CRT data to data collected during a past risk reduction, Quiet Supersonic Flights 2018 (QSF18), this short note is intended to create awareness. It documents an instance in which a reported population average dose-response relationship derived from QSF18 single event data was distorted by the integral approximation applied in likelihood-based methods. We review some of the available literature on the topic, compare the outputs of several different computational approaches implemented in available statistical software, and present simple corrective actions that may be useful during the Quesst mission.

dose-response model↗

Assessment of the Uniqueness of Wind Tunnel Strain-Gage Balance Load Predictions

A new test was developed to assess the uniqueness of wind tunnel strain-gage balance load predictions that are obtained from regression models of calibration data. The test helps balance users to gain confidence in load predictions of non-traditional balance designs. It also makes it possible to better evaluate load predictions of traditional balances that are not used as originally intended. The test works for both the Iterative and Non-Iterative Methods that are used in the aerospace testing community for the prediction of balance loads. It is based on the hypothesis that the total number of independently applied balance load components must always match the total number of independently measured bridge outputs or bridge output combinations. This hypothesis is supported by a control volume analysis of the inputs and outputs of a strain-gage balance. It is concluded from the control volume analysis that the loads and bridge outputs of a balance calibration data set must separately be tested for linear independence because it cannot always be guaranteed that a linearly independent load component set will result in linearly independent bridge output measurements. Simple linear math models for the loads and bridge outputs in combination with the variance inflation factor are used to test for linear independence. A highly unique and reversible mapping between the applied load component set and the measured bridge output set is guaranteed to exist if the maximum variance inflation factor of both sets is less than the literature recommended threshold of five. Data from the calibration of a six{component force balance is used to illustrate the application of the new test to real-world data.

strain-gage balance↗

Human sympathetic and vagal baroreflex responses to sequential nitroprusside and phenylephrine

We evaluated a method of baroreflex testing involving sequential intravenous bolus injections of nitroprusside followed by phenylephrine and phenylephrine followed by nitroprusside in 18 healthy men and women, and we drew inferences regarding human sympathetic and vagal baroreflex mechanisms. We recorded the electrocardiogram, photoplethysmographic finger arterial pressure, and peroneal nerve muscle sympathetic activity. We then contrasted least squares linear regression slopes derived from the depressor (nitroprusside) and pressor (phenylephrine) phases with 1) slopes derived from spontaneous fluctuations of systolic arterial pressures and R-R intervals, and 2) baroreflex gain derived from cross-spectral analyses of systolic pressures and R-R intervals. We calculated sympathetic baroreflex gain from integrated muscle sympathetic nerve activity and diastolic pressures. We found that vagal baroreflex slopes are less when arterial pressures are falling than when they are rising and that this hysteresis exists over pressure ranges both below and above baseline levels. Although pharmacological and spontaneous vagal baroreflex responses correlate closely, pharmacological baroreflex slopes tend to be lower than those derived from spontaneous fluctuations. Sympathetic baroreflex slopes are similar when arterial pressure is falling and rising; however, small pressure elevations above baseline silence sympathetic motoneurons. Vagal, but not sympathetic baroreflex gains vary inversely with subjects' ages and their baseline arterial pressures. There is no correlation between sympathetic and vagal baroreflex gains. We recommend repeated sequential nitroprusside followed by phenylephrine doses as a simple, efficientmeans to provoke and characterize human vagal and sympathetic baroreflex responses.

NASA Discipline Cardiopulmonary↗

Right and left ventricular volumes in vitro by a new nongeometric method

We present an evaluation of a new nongeometric technique for calculating right and left ventricular volumes. This method calculates ventricular chamber volumes from multiple cross-sectional echocardiographic views taken from a single point as the echo beam is tilted progressively through the ventricle. Right and left ventricular volumes are calculated from both the approximate short axis and approximate apical position on 20 in vitro human hearts and compared with the actual chamber volumes. The results for both ventricles from both positions are excellent. Correlation coefficients are > 0.95 for all positions; the standard errors are in the range of 5 to 7 mL and the slopes and intercepts for the regression lines are not significantly different from 1 and 0, respectively (except for the left ventricular short-axis intercept). For all positions, approximately 6 to 8 views are needed for peak accuracy (7.5 degrees to 10 degrees separation). This approach offers several advantages. No geometric assumptions about ventricular shape are made. All images are acquired from a single point (or window), and the digitized points can be used to make a three-dimensional reconstruction of the ventricle. Also, during the calculations a volume distribution curve for the ventricle is produced. The shape of this curve can be characteristic for certain situations (ie, right ventricle, short axis) and can be used to make new simple equations for calculating volume. We conclude that this is an accurate nongeometric method for determining both right and left ventricular volumes in vitro.

NASA Discipline Cardiopulmonary↗

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗