Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 1,171 records · Page 65

Radiation effects on atherosclerosis in atomic bomb survivors: a cross-sectional study using structural equation modeling

Past reports indicated that total-body irradiation at low to moderate doses could be responsible for cardiovascular disease risks, but the mechanism remains unclear. The purpose of this study was to investigate the association between radiation exposure and atherosclerosis, an underlying pathology of cardiovascular diseases, in the Japanese atomic bomb survivors. We performed a cross-sectional study measuring 14 clinical-physiological atherosclerosis indicators during clinical exams from 2010 to 2014 in 3274 participants of the Adult Health Study cohort. Multivariable analyses were performed by using a structural equation model with latent factors representing underlying atherosclerotic pathologies: (1) arterial stiffness, (2) calcification, and (3) plaque as measured with indicators chosen a priori on the basis of clinical-physiological knowledge. Radiation was linearly associated with calcification (standardized coefficient per Gy 0.15, 95 % confidence interval: CI [0.070, 0.23]) and plaque (0.11, 95 % CI [0.029, 0.20]), small associations that were comparable to about 2 years of aging per Gy of radiation exposure, but not with arterial stiffness (0.036, 95 % CI [-0.025, 0.095]). The model fitted better and had narrower confidence intervals than separate ordinary regression models explaining individual indicators independently. The associations were less evident when the dose range was restricted to a maximum of 2 or 1 Gy. By combining individual clinical-physiological indicators that are correlated because of common, underlying atherosclerotic pathologies, we found a small, but significant association of radiation with atherosclerosis.

59 BASIC BIOLOGICAL SCIENCES↗

A behavior and physiology-based decision support tool to predict thermal comfort and stress in non-pregnant, mid-gestation, and late-gestation sows

Background: Although thermal indices have been proposed for swine, none to our knowledge differentiate by reproductive stage or predict thermal comfort using behavioral and physiological data. The study objective was to develop a behavior and physiology-based decision support tool to predict thermal comfort and stress in multiparous (3.28 ± 0.81) non-pregnant (n = 11), mid-gestation (n = 13), and late-gestation (n = 12) sows. Results: Regression analyses were performed using PROC MIXED in SAS 9.4 to determine the optimal environmental indicator [dry bulb temperature (T DB ) and dew point] of heat stress (HS) in non-pregnant, mid-gestation, and late-gestation sows with respiration rate (RR) and body temperature (T B ) successively used as the dependent variable in a cubic function. A linear relationship was observed for skin temperature (T S ) indicating that TDB rather than the sow HS response impacted TS and so TS was excluded from further analyses. Reproductive stage was significant for all analyses (P < 0.05). Heat stress thresholds for each reproductive stage were calculated using the inflections points of RR for mild HS and T B for moderate and severe HS. Mild HS inflection points differed for non-pregnant, mid-gestation, and late gestation sows and occurred at 25.5, 25.1, and 24.0 °C, respectively. Moderate HS inflection points differed for non-pregnant, mid-gestation, and late gestation sows and occurred at 28.1, 27.8, and 25.5 °C, respectively. Severe HS inflection points were similar for non-pregnant and mid-gestation sows (32.9 °C) but differed for late-gestation sows (30.8 °C). These data were integrated with previously collected behavioral thermal preference data to estimate the T DB that non-pregnant, mid-gestation, and late-gestation sows found to be cool (T DB < T DB preference range), comfortable (T DB = T DB preference range), and warm (T DB preference range < T DB < mild HS). Conclusions: The results of this study provide valuable information about thermal comfort and thermal stress thresholds in sows at three reproductive stages. The development of a behavior and physiology-based decision support tool to predict thermal comfort and stress in non-pregnant, mid-gestation, and late-gestation sows is expected to provide swine producers with a more accurate means of managing sow environments.

59 BASIC BIOLOGICAL SCIENCES↗

Relationships between CH4 emission, biomass, and CO2 exchange in a subtropical grassland

Methane flux was linearly correlated with plant biomass (r = 0.97, n = 6 and r = 0.95, n = 8) at two locations in a Florida Everglades Cladium marsh. One location, which had burned 4 months previously, exhibited a greater increase in methane flux as a function of biomass relative to sites at an unburned location. However, methane flux data from both sites fit a single regression (r = 0.94, n = 14) when plotted against net CO2 exchange suggesting that either methanogenesis in Everglades marl sediments is fueled by root exudation below ground, or that factors which enhance photosynthetic production and plant growth are also correlated with methane production and flux in this oligotrophic environment. The data presented are the first to show a direct relationship between spatial variability in plant biomass, net ecosystem production, and methane emission in a natural wetland.

Whiting, Gary J.↗

Calibration and Data Analysis Recommendations for Three-Component Moment Balances

Fundamental characteristics of design, calibration, and application of three-component moment balances are investigated in great detail. These balances are typically used to determine loads on control surfaces, canards, or other parts that are attached to a wind tunnel model. First, three different descriptions of the load state of a moment balance are reviewed. Then, load transformations between different load formats and the combined load diagram for two of the three load components are discussed. An error analysis showed that it is critical to maximize the product of the distance between the bending moment gages and their sensitivities in order to minimize the overall error in the normal force prediction. In addition, it is important to apply a sufficient number of calibration loadings near the first bending moment gage. Then, unwanted near-linear dependencies between the two bending moment gage outputs can be avoided. The error in the bending moment prediction is also investigated that results from the elastic deformation of the metric part of the balance under load. Finally, the application of the Non-Iterative Method to three-component moment balance calibration data is described in order to obtain regression models that can be used to predict loads from measured outputs during a wind tunnel test.

Strain-Gage Balance↗

Monte Carlo simulation of the effect of soil moisture variation on the microwave emission from soils

In this paper, results of a Monte Carlo simulation of the effect of noise on the relationship between the microwave emissivity of soil and its moisture content are presented. It is found that whenever the magnitude of the noise for the independent variable, in this case the soil moisture, is increased, both the slope of the regression and the correlation coefficient decrease. In particular, when the noise has magnitude equivalent to a coefficient of variation of 0.25, the slope and correlation coefficient are in good agreement with those obtained from the data of a 21-cm airborne microwave radiometer which was flown over a test site in hand county, South Dakota. The comparison was made using a linear relationship to determine the estimated emissivity from the ground measurements of soil moisture. The linear relationship was derived from a radiative transfer model calculation of the microwave emissivities using realistic soil-moisture profiles. The effect of surface roughness was included in the relationship, and the variability of the surface roughness was also simulated by a Monte Carlo technique.

Mo, T.↗

Prediction of elemental creep

Cyclic and steady-state creep tests were performed to provide data which were used to develop predictive equations. These equations, describing creep as a function of stress, temperature, and time, were developed through the use of a least squares regression analyses computer program for both the steady-state and cyclic data sets. Comparison of the data from the two types of tests, revealed that there was no significant difference between the cyclic and steady-state creep strains for the L-605 sheet under the experimental conditions investigated (for the same total time at load). Attempts to develop a single linear equation describing the combined steady-state and cyclic creep data resulted in standard errors of estimates higher than obtained for the individual data sets. A proposed approach to predict elemental creep in metals uses the cyclic creep equation and a computer program which applies strain and time hardening theories of creep accumulation.

Davis, J. W.↗

Dynamic Density: An Air Traffic Management Metric

The definition of a metric of air traffic controller workload based on air traffic characteristics is essential to the development of both air traffic management automation and air traffic procedures. Dynamic density is a proposed concept for a metric that includes both traffic density (a count of aircraft in a volume of airspace) and traffic complexity (a measure of the complexity of the air traffic in a volume of airspace). It was hypothesized that a metric that includes terms that capture air traffic complexity will be a better measure of air traffic controller workload than current measures based only on traffic density. A weighted linear dynamic density function was developed and validated operationally. The proposed dynamic density function includes a traffic density term and eight traffic complexity terms. A unit-weighted dynamic density function was able to account for an average of 22% of the variance in observed controller activity not accounted for by traffic density alone. A comparative analysis of unit weights, subjective weights, and regression weights for the terms in the dynamic density equation was conducted. The best predictor of controller activity was the dynamic density equation with regression-weighted complexity terms.

Laudeman, I. V.↗

Reinforcement Learning via Gaussian Processes with Neural Network Dual Kernels

While deep neural networks (DNNs) and Gaussian Processes (GPs) are both popularly utilized to solve problems in reinforcement learning, both approaches feature undesirable drawbacks for challenging problems. DNNs learn complex non-linear embeddings, but do not naturally quantify uncertainty and are often data-inefficient to train. GPs infer posterior distributions over functions, but popular kernels exhibit limited expressivity on complex and high-dimensional data. Fortunately, recently discovered conjugate and neural tangent kernel functions encode the behavior of overparameterized neural networks in the kernel domain. We demonstrate that these kernels can be efficiently applied to regression and reinforcement learning problems by analyzing a baseline case study.We apply GPs with neural network dual kernels to solve reinforcement learning tasks for the first time. We demonstrate, using the well understood mountain-car problem, that GPs empowered with dual kernels perform at least as well as those using the conventional radial basis function kernel. Finally, we conjecture that by inheriting the probabilistic rigor of GPs and the powerful embedding properties of DNNs, GPs using NN dual kernels will empower future reinforcement learning models on difficult domains.

97 MATHEMATICS AND COMPUTING↗

Cardiovascular responses of snakes to hypergravity

Snakes have provided useful vertebrate models for understanding circulatory adaptation to gravity, attributable to their elongate body shape and evolutionary diversificaton in terms of ecology and behavior. Recently we have studied cardiovascular responses of snakes to hypergravic acceleration forces produced acutely in the head-to-tail direction (+Gz) on a short-arm centrifuge. Snakes were held in a nearly straight position within a horizontal plastic tube and subjected to a linear force gradient during acceleration. Carotid blood flow provided an integrated measure of cardiovascular performance. Thus, cardiovascular tolerance of snakes to stepwise increments of Gz was measured as the caudal Gz force at which carotid blood flow ceased. Tolerance to increasing Gz varies according to adaptive evolutionary history inferred from the ecology and behavior of species. With respect to data for six species we investigated, multiple regression analysis demonstrates that Gz tolerance correlates with gravitational habitat, independently of body length. Relative to aquatic and non-climbing species, carotid blood flow is better maintained in arboreal or scansorial species, which tolerate hypergravic forces of +2 to +3.5 Gz. Additionally, semi-arboreal rat snakes (Elaphe obsoleta) exhibit plasticity of responses to long-term, intermittent +1.5 Gz stress. Compared to non-acclimated controls, acclimated snakes show greater increases of heart rate during head-up tilt or acceleration, greater sensitivity of arterial pressure to circulating catecholamines, higher blood levels of prostaglandin ratios favorable to maintenance of arterial blood pressure, and medial hypertrophy in major arteries and veins. As in other vertebrates, Gz tolerance of snakes is enhanced by acclimation, high arterial pressure, comparatively large blood volume, and body movements. Vascular studies of snakes suggest the importance to acclimation of local responses involving vascular tissue, in addition to centrally mediated responses to fluid shifts.

Review, Tutorial↗

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗

Helicopter stability and control modeling improvements and verification on two helicopters

A linearized model of helicopter flight dynamics is developed which includes the flapping, lead-lag, and dynamic inflow degrees of freedom (DOF). The model is a combination of analytical terms and numerically determined stability derivatives, and is used to investigate the importance of the rotor DOF to stability and control modeling. The results show that the rotor DOF can have a significant impact on some of the natural modes in a linear model. The flap and dynamic inflow DOF show the greatest influence. Flapping exhibits strong coupling to the body, dynamic inflow, and to lead-lag to a lesser extent. Dynamic inflow tends to damp the high-frequency flapping modes, and reduces the damping on coupled body-flap motion. Dynamic inflow also couples to the flapping motion to produce complex roots. With body-flap and lag regressing modes as exceptions, the results show essentially similar behavior for most modes of articulated and hingeless rotor helicopters.

Schrage, D. P.↗

Convolutional Neural Network for Transition Modeling Based on Linear Stability Theory

Transition prediction is an important aspect of aerodynamic design because of its impact on skin friction and potential coupling with flow separation characteristics. Traditionally, the modeling of transition has relied on correlation-based empirical formulas based on integral quantities such as the shape factor of the boundary layer. However, in many applications of computational fluid dynamics, the shape factor is not straightforwardly available or not well-defined. We propose using the complete velocity profile along with other quantities (e.g., frequency, Reynolds number) to predict the perturbation amplification factor. While this can be achieved with regression models based on a classical fully connected neural network, such a model can be computationally more demanding. We propose a novel convolutional neural network inspired by the underlying physics as described by the stability equations. Specifically, convolutional layers are first used to extract integral quantities from the velocity profiles, and then fully connected layers are used to map the extracted integral quantities, along with frequency and Reynolds number, to the output (amplification ratio). Numerical tests on classical boundary layers clearly demonstrate the merits of the proposed method. More importantly, we demonstrate that, for Tollmien-Schlichting instabilities in two-dimensional, low-speed boundary layers, the proposed network encodes information in the boundary layer profiles into an integral quantity that is strongly correlated to a well-known, physically defined parameter – the shape factor.

Laminar-turbulent transition↗

GMLS-Nets: A Framework for Learning from Unstructured Data

Data fields sampled on irregularly spaced points arise in many applications in the sciences and engineering. For regular grids, Convolutional Neural Networks (CNNs) have been successfully used to gaining benefits from weight sharing and invariances. We generalize CNNs by introducing methods for data on unstructured point clouds based on Generalized Moving Least Squares (GMLS). GMLS is a non-parametric technique for estimating linear bounded functionals from scattered data, and has recently been used in the literature for solving partial differential equations. By parameterizing the GMLS estimator, we obtain learning methods for operators with unstructured stencils. In GMLS-Nets the necessary calculations are local, readily parallelizable, and the estimator is supported by a rigorous approximation theory. We show how the framework may be used for unstructured physical data sets to perform functional regression to identify associated differential operators and to regress quantities of interest. The results suggest the architectures to be an attractive foundation for data-driven model development in scientific machine learning applications

97 MATHEMATICS AND COMPUTING↗

Differential credibility assessment for statistical downscaling

Climate science is increasingly using (i) ensembles of climate projections from multiple models derived using different assumptions and/or scenarios and (ii) process-oriented diagnostics of model fidelity. Efforts to assign differential credibility to projections and/or models are also rapidly advancing. A framework to quantify and depict the credibility of statistically downscaled model output is presented and demonstrated. Here, the approach employs transfer functions in the form of robust and resilient generalized linear models applied to downscale daily minimum and maximum temperature anomalies at 10 locations using predictors drawn from ERA-Interim reanalysis and two global climate models (GCM; GFDL-ESM2M and MPI-ESM-LR). The downscaled time series are used to derive several impact relevant CLIMDEX temperature indices that are assigned credibility based on (1) the reproduction of relevant large-scale predictors by the GCMs (i.e. fraction of regression beta-weights derived from predictors that are well-reproduced) and (2) the degree of variance in the observations reproduced in the downscaled series following application of a new variance inflation technique. Credibility of the downscaled predictands varies across locations, between the two GCM and is generally higher for minimum temperature than maximum temperature. The differential credibility assessment framework demonstrated here is easy to use and flexible. It can be applied as is to inform decision makers regarding projection confidence, and/or extended to include other components of the transfer functions, and/or used to weight members of a statistically downscaled ensemble.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of constrained linear inversions to the selection of the Lagrange multiplier

The influence of the choice of the Lagrange multiplier on constrained linear inversions is explored, with reference made to applications in inferring the columnar aerosol size distributions from spectral aerosol optical depth measurements. A range of the Lagrange multiplier is examined to find all positive solutions for the solution vector, which represents modifying factors to the assumed form of the size distribution. An iterative method is devised to constrain the calculations to consideration of only positive quantities and a requirement that the regression fit to data be consistent with measurement errors. The determination of the variances and covariances is formulated and applied to existing data sets for optical depth. Variances in the solution are found to be large for particle radii when the information content of the data is small.

King, M. D.↗

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly

Abstract Background Looking older for one’s chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. Objectives To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. Methods We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5–87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). Results We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one’s age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62–0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77–0.95), less age-related hearing loss (model 2; B = −0.76, 95% CI −1.35 to −0.17) and fewer cataracts (OR 0.84, 95% CI 0.73–0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04–0.10). Conclusions PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

Mekić, Selma↗

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.↗

Greenhouse Gas Concentration Data Recovery Algorithm for a Low Cost, Laser Heterodyne Radiometer

The goal of a coordinated effort between groups at GWU and NASA GSFC is the development of a low-cost, global, surface instrument network that continuously monitors three key carbon cycle gases in the atmospheric column: carbon dioxide (CO2), methane (CH4), carbon monoxide (CO), as well as oxygen (O2) for atmospheric pressure profiles. The network will implement a low-cost, miniaturized, laser heterodyne radiometer (mini-LHR) that has recently been developed at NASA Goddard Space Flight Center. This mini-LHR is designed to operate in tandem with the passive aerosol sensor currently used in AERONET (a well established network of more than 450 ground aerosol monitoring instruments worldwide), and could be rapidly deployed into this established global network. Laser heterodyne radiometry is a well-established technique for detecting weak signals that was adapted from radio receiver technology. Here, a weak light signal, that has undergone absorption by atmospheric components, is mixed with light from a distributed feedback (DFB) telecommunications laser on a single-mode optical fiber. The RF component of the signal is detected on a fast photoreceiver. Scanning the laser through an absorption feature in the infrared, results in a scanned heterodyne signal io the RF. Deconvolution of this signal through the retrieval algorithm allows for the extraction of altitude contributions to the column signal. The retrieval algorithm is based on a spectral simulation program, SpecSyn, developed at GWU for high-resolution infrared spectroscopies. Variations io pressure, temperature, composition, and refractive index through the atmosphere; that are all functions of latitude, longitude, time of day, altitude, etc.; are modeled using algorithms developed in the MODTRAN program developed in part by the US Air Force Research Laboratory. In these calculations the atmosphere is modeled as a series of spherically symmetric shells with boundaries specified at defined altitudes. Temperature, pressure, and species mixing ratios are defined at these boundaries. Between the boundaries, temperature is assumed to vary linearly with altitude while pressure (and thus gas density) vary exponentially. The observed spectrum at the LHR instrument will be the integration of the contributions along this light path. For any absorption measurement the signal at a particular spectral frequency is a linear combination of spectral line contributions from several species. For each species that might absorb in a spectral region, we have pre-calculated its contribution as a function of temperature and pressure. The integrated path absorption spectrum can then by calculated using the initial sun angle (from location, date, and time) and assumptions about pressure and temperature profiles from an atmospheric model. The modeled spectrum is iterated to match the experimental observation using standard multilinear regression techniques. In addition to the layer concentrations, the numerical technique also provides uncertainty estimates for these quantities as well as dependencies on assumptions inherent in the atmospheric models.

Miller, J. Houston↗