Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear mixed models”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-resolution (30-m) urban land cover projections for Los Angeles California Urban Area: 2010 to 2100 under SSP3 and SSP5 [Updated simulations based on population-driven urban intensity transitions]

These data (v3) are updated from previous versions (1 and 2) in that they include consider the effects of population on transitions in urban land intensity. This leads to more reasonable differences in urban land projections under variant SSPs. For the present dataset, both SSP3 and SSP5 are provided. These data represent simulations of future land use and land cover (LULCC) for Los Angeles urban area (U.S. Census Bureau defined area) as raster tiff images at a 30-m pixel resolution and at decadal time steps from 2010 to 2100. LULCC classes in this product follow the National Land Cover Dataset (NLCD) classification. NLCD 21-24 correspond to open developed, low developed, medium developed, and high developed urban land classes, respectively. Only urban land cover classes (NLCD class 21, 22, 23, and 24) are dynamic over time; however, all NLCD classes are included in the final product. Therefore, NLCD classes that do not convert to an urban class will be similar to year 2000. The products were developed using a hybridized statistical and cellular automata approach. Linear mixed models (LMMs) were used to estimate future urban land budgets based on 1-km urban land fraction projections from Gao and Pesaresi (2021), whereas separate generalized linear mixed models (GLMMs) were used to estimate shifts in urban land intensities based on retrospective shifts in NLCD urban class intensities over a 20- year period. Based on urban land allocations from the statistical models, a cellular-automata and downscaling routine was used to simulate dynamic urban land expansion at a 30-m resolution based on suitability criteria. Scenarios of future urban landcover change projections include variant solutions for the Shared Socioeconomic Pathway 5 (SSP5) and SSP 3 based on different population assumptions, different land use intensification assumptions, variable land zoning constraints, and iterative adjustments to correct for over allocation of urban expansion across decadal time periods from 2010 to 2100. This results in 320 raster products.

Land↗
DOE OSTI2020

Assessment of Demographic, Genetic, and Imaging Variables Associated With Brain Resilience and Cognitive Resilience to Pathological Tau in Patients With Alzheimer Disease

Importance: Better understanding is needed of the degree to which individuals tolerate Alzheimer disease (AD)-like pathological tau with respect to brain structure (brain resilience) and cognition (cognitive resilience). Objective: To examine the demographic (age, sex, and educational level), genetic (APOE-ε4 status), and neuroimaging (white matter hyperintensities and cortical thickness) factors associated with interindividual differences in brain and cognitive resilience to tau positron emission tomography (PET) load and to changes in global cognition over time. Design, setting, an participants: In this cross-sectional, longitudinal study, tau PET was performed from June 1, 2014, to November 30, 2017, and global cognition monitored for a mean [SD] interval of 2.0 [1.8] years at 3 dementia centers in South Korea, Sweden, and the United States. The study included amyloid-β-positive participants with mild cognitive impairment or AD dementia. Data analysis was performed from October 26, 2018, to December 11, 2019. Exposures: Standard dementia screening, cognitive testing, brain magnetic resonance imaging, amyloid-β PET and cerebrospinal fluid analysis, and flortaucipir (tau) labeled with fluor-18 (18F) PET. Main outcomes and measures: Separate linear regression models were performed between whole cortex [18F]flortaucipir uptake and cortical thickness, and standardized residuals were used to obtain a measure of brain resilience. The same procedure was performed for whole cortex [18F]flortaucipir uptake vs Mini-Mental State Examination (MMSE) as a measure of cognitive resilience. Bivariate and multivariable linear regression models were conducted with age, sex, educational level, APOE-ε4 status, white matter hyperintensity volumes, and cortical thickness as independent variables and brain and cognitive resilience measures as dependent variables. Linear mixed models were performed to examine whether changes in MMSE scores over time differed as a function of a combined brain and cognitive resilience variable. Results: A total of 260 participants (145 [55.8%] female; mean [SD] age, 69.2 [9.5] years; mean [SD] MMSE score, 21.9 [5.5]) were included in the study. In multivariable models, women (standardized β = -0.15, P = .02) and young patients (standardized β = -0.20, P = .006) had greater brain resilience to pathological tau. Higher educational level (standardized β = 0.23, P < .001) and global cortical thickness (standardized β = 0.23, P < .001) were associated with greater cognitive resilience to pathological tau. Linear mixed models indicated a significant interaction of brain resilience × cognitive resilience × time on MMSE (β [SE] = -0.235 [0.111], P = .03), with steepest slopes for individuals with both low brain and cognitive resilience. Conclusions and relevance: Results of this study suggest that women and young patients with AD have relative preservation of brain structure when exposed to neocortical pathological tau. Interindividual differences in resilience to pathological tau may be important to disease progression because participants with both low brain and cognitive resilience had the most rapid cognitive decline over time.

60 APPLIED LIFE SCIENCES↗
NASA NTRSDate not supplied

Linear Mixed-Effects Models for Human-in-the-Loop Tracking Experiment Data

Linear mixed-effects models provide several benefits over more traditional statistical inference tests that are particularly useful for most human-in-the-loop tracking experiment data. However, surprisingly, mixed models are virtually not used for the analysis of tracking experiment data. This paper uses of linear mixed-effects models to analyze combined tracking data from two previous human-in-the-loop roll tracking experiments that compared control behavior metrics collected in both a research aircraft and a motion-base simulator. In the experiments, pilots' behavior under 10 different motion configurations with varying motion filter gains and break frequencies was evaluated and compared to that in the real aircraft. The linear mixed-effects model analysis on the combined dataset confirmed the main statistical outcomes of the individual experiments. The main benefits of mixed models for this type of data were demonstrated by successfully combining data from two experiments that used different experimental conditions and of which one had an additional apparatus and the other a missing participant. Finally, the mixed-model analysis was able to explicitly test for scientifically relevant statistical differences in the dependent measures between the aircraft and simulator, as well as both experiments.

manual control↗
NASA NTRSDate not supplied

Developing Rationale for Experimental Designs Employed During Forthcoming NASA Quesst Mission Community Noise Campaigns

Since the publication of R.A. Fisher’s “The Arrangement of Field Experiments” in 1926, the enumerated principles that contributed to the improvement of agricultural field experiments in the early 20th century (randomization, replication, blocking, and appropriate variation of factors) remain hallmarks of all well-designed clinical trials and scientific endeavors in the 21st century. Beginning in 2026, NASA will conduct community noise campaigns with the first-of-kind X-59 experimental aircraft to 1) demonstrate the possibility of quiet supersonic flight over land and 2) to collect live data about annoyance (a categorical response) in relation to estimated noise levels (an experimental factor) produced by the new noise phenomenon, a low-noise “sonic thump”. The resulting predicted relationship between perceptual response and estimated noise level, e.g., mixed logistic regression or related models in the class of generalized linear (mixed) models, is believed to be a useful policy tool that domestic and international regulators can use when deciding whether existing prohibitions of supersonic flight over land can be replaced with a noise-based limit. The X-59 is engineered with the expectation that annoyance in the tested ranges should be a rare outcome, and, consequently, NASA convened an interdisciplinary team to examine and refine the rationale underlying the experimental design for future community campaigns. In this presentation, we review relevant literature on experiments with rare binary outcomes, relate the deliberations of the NASA team to tenets of good experimental design, and highlight several challenges and operational realities of this ambitious campaign.

Design of Experiments; blocking; randomization; co↗
DOE OSTI2017

Effects of hydrokinetic turbine sound on the behavior of four species of fish within an experimental mesocosm

The development of hydrokinetic energy technologies (e.g., tidal turbines) has raised concern over the potential impacts of underwater sound produced by hydrokinetic turbines on fish species likely to encounter these turbines. To assess the potential for behavioral impacts, we exposed four species of fish to varying intensities of recorded hydrokinetic turbine sound in a semi-natural environment. Although we tested freshwater species (redhorse suckers [Moxostoma spp], freshwater drum [Aplondinotus grunniens], largemouth bass [Micropterus salmoides], and rainbow trout [Oncorhynchus mykiss]), these species are also representative of the hearing physiology and sensitivity of estuarine species that would be affected at tidal energy sites. Here, we evaluated changes in fish position relative to different intensities of turbine sound as well as trends in location over time with linear mixed-effects and generalized additive mixed models. We also evaluated changes in the proportion of near-source detections relative to sound intensity and exposure time with generalized linear mixed models and generalized additive models. Models indicated that redhorse suckers may respond to sustained turbine sound by increasing distance from the sound source. Freshwater drum models suggested a mixed response to turbine sound, and largemouth bass and rainbow trout models did not indicate any likely responses to turbine sound. Lastly, findings highlight the importance for future research to utilize accurate localization systems, different species, validated sound transmission distances, and to consider different types of behavioral responses to different turbine designs and to the cumulative sound of arrays of multiple turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2015

Posterior propriety for hierarchical models with log-likelihoods that have norm bounds

Statisticians often use improper priors to express ignorance or to provide good frequency properties, requiring that posterior propriety be verified. Our paper addresses generalized linear mixed models, GLMMs, when Level I parameters have Normal distributions, with many commonly-used hyperpriors. It provides easy-to-verify sufficient posterior propriety conditions based on dimensions, matrix ranks, and exponentiated norm bounds, ENBs, for the Level I likelihood. Since many familiar likelihoods have ENBs, which is often verifiable via log-concavity and MLE finiteness, our novel use of ENBs permits unification of posterior propriety results and posterior MGF/moment results for many useful Level I distributions, including those commonly used with multilevel generalized linear models, e.g., GLMMs and hierarchical generalized linear models, HGLMs. Furthermore, those who need to verify existence of posterior distributions or of posterior MGFs/moments for a multilevel generalized linear model given a proper or improper multivariate F prior as in Section 1 should find the required results in Sections 1 and 2 and Theorem 3 (GLMMs), Theorem 4 (HGLMs), or Theorem 5 (posterior MGFs/moments).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Interpretable Machine Learning for Characterizing Electric Vehicle Charging Behavior: Insights from Real-World Data

As electric vehicle (EV) adoption rises globally, concerns about the impact on aging electrical grids grow, particularly regarding the charging behavior of EV drivers. This study analyzes real-world driving and charging data from Ford battery electric vehicles (BEVs) collected between 2018 and 2019 to develop interpretable models that characterize charging behavior and quantify influencing factors. Prior research has relied on assumptions regarding driver behavior, often overlooking actual charging patterns. By employing generalized linear mixed models (GLMMs), this work offers insights into how various elements, such as next trip distance and state of charge (SOC), influence charging decisions. The dataset comprises over three million park-trip pairs from 1,997 vehicles, revealing that features related to driving behavior significantly dictate charging behavior, while infrastructure and regional factors have lesser impacts. The findings suggest that existing simulation models may oversimplify EV charging behavior assumptions. This work utilizes real-world EV driving and charging data to train interpretable models that describe charging behavior and quantify the factors most associated with how drivers use charging infrastructure. This research underscores the need for interpretable, data-driven methodologies to inform future EV infrastructure planning and grid management.

29 - ENERGY PLANNING, POLICY AND ECONOMY↗
DOE OSTI2021

Breeding latitude is associated with the timing of nesting and migration around the annual calendar among Purple Martin (Progne subis) populations

Timing is critical for species that migrate seasonally, as timing mismatches can lead to reduced reproductive success and population declines. We investigated the effect of breeding latitude on the migration timing of a long-distance migratory songbird (Purple Martin, Progne subis subis) by using light-level geolocators deployed across breeding sites between 26° and 53° latitude (n = 329). We used nesting data (1st egg dates; n=28,165) to explore latitudinal patterns in the timing of breeding. Using generalized linear mixed models (GLMMs) we found that breeding latitude was highly associated with timing around the calendar, explaining population-level timing across both spring and fall migration. Populations breeding at more northern latitudes migrated later in fall (mean days/degree latitude =2), in spring (mean days/degree latitude =3), and nested later (days/degree latitude = 3) than populations breeding at more southern latitudes. Generalized additive mixed models (GAMMs) revealed significant non-linearity between breeding latitude and spring and fall migration timing, and the timing of nesting. In spring, the most northern breeders (> 40 degrees N) were more synchronous in their timing, departing their nonbreeding sites and arriving earlier at their breeding sites relative to their breeding latitude. These most northern breeders also had similar nest timing and fall migration departure timing. Breeding latitude was also associated with seasonal partitioning of migratory and sedentary periods. The most northern breeders spent 10% more time on fall migration, 2% more at wintering grounds, 1% more on spring migration, and 13% less time at their breeding sites. Our results support the notion that selection for timing at breeding latitudes can have year-round effects on the timing of seasonal events in long-distance migratory songbirds. Future studies should examine whether these timing patterns associated with latitude limit phenotypic responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Numerical Simulations of Two-Phase Reacting Flow in a Single-Element Lean Direct Injection (LDI) Combustor Using NCC

A series of numerical simulations of Jet-A spray reacting flow in a single-element lean direct injection (LDI) combustor have been conducted by using the National Combustion Code (NCC). The simulations have been carried out using the time filtered Navier-Stokes (TFNS) approach ranging from the steady Reynolds-averaged Navier-Stokes (RANS), unsteady RANS (URANS), to the dynamic flow structure simulation (DFS). The sub-grid model employed for turbulent mixing and combustion includes the well-mixed model, the linear eddy mixing (LEM) model, and the filtered mass density function (FDF/PDF) model. The starting condition of the injected liquid spray is specified via empirical droplet size correlation, and a five-species single-step global reduced mechanism is employed for fuel chemistry. All the calculations use the same grid whose resolution is of the RANS type. Comparisons of results from various models are presented.

Liu, Nan-Suey↗
DOE OSTI2021

Environmental predictors of pulmonary nontuberculous mycobacteria (NTM) sputum positivity among persons with cystic fibrosis in the state of Florida

Nontuberculous mycobacteria (NTM) are opportunistic human pathogens that are commonly found in soil and water, and exposure to these organisms may cause pulmonary nontuberculous mycobacterial disease. Persons with cystic fibrosis (CF) are at high risk for developing pulmonary NTM infections, and studies have shown that prolonged exposure to certain environments can increase the risk of pulmonary NTM. It is therefore important to determine the risk associated with different geographic areas. Using annualized registry data obtained from the Cystic Fibrosis Foundation Patient Registry for 2010 through 2017, we conducted a geospatial analysis of NTM infections among persons with CF in Florida. A Bernoulli model in SaTScan was used to identify clustering of ZIP codes with higher than expected numbers of NTM culture positive individuals. Generalized linear mixed models with a binomial distribution were used to test the association of environmental variables and NTM culture positivity. We identified a significant cluster of M . abscessus and predictors of NTM sputum positivity, including annual precipitation and soil mineral levels.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

MINE: maximally informative next experiment—toward a new GWAS experimental design and methodology

Abstract The computational methodology of Genome Wide Association Studies (GWAS) currently has several limitations: (i) the number of observations (rows) on a quantitative trait tends to be smaller than the number of single nucleotide polymorphisms (SNPs) (columns) in the design matrix; (ii) each SNP is usually modeled separately, failing to acknowledge interaction between each other (ie epistasis); (iii) there is implicit linkage disequilibrium (LD) between neighboring SNPs due to their linkage. To overcome these issues, we developed a tool that uses ensemble methods to fit mixed linear models to GWAS data, and these ensemble methods include the development of a new experimental design approach in GWAS, which uses the resultant models and data to select the next informative experiment over time. This new adaptive and staged approach for GWAS experimental design was developed and tested in a 3 yr adaptive model-guided discovery experiment against a fixed classical design. In Sorghum bicolor a total of 79, 86, and 78 accessions were tested in years 1, 2, and 3, respectively out of 343 accessions available in the Bioenergy Association Panel (BAP) each identified for 232,303 SNPs, 1 every 2–3 kb in the genomes. We demonstrated the feasibility of MINE enacted with 8 people in the field per year over 3 yr vs in 1 large classical design enacted with 20 people in 1 yr. The MINE results for chromosomal regions identified controlling dry weight were confirmed against results from previous sorghum GWAS experiments and 1 large classical design for the BAP panel.

Genetics & Heredity↗
NASA NTRSDate not supplied

Mitigating the Impacts of Measurement Error in the Quesst Mission Community Noise Study

Beginning in 2025, the NASA Quesst mission will conduct a series of community response tests involving flyovers of the X-59 aircraft at select localities across the United States. Several waves of a longitudinal survey will be administered over approximately one month of testing in order to capture perceptual responses to low-amplitude sonic booms, or “sonic thumps”. Simultaneously, noise exposure levels will be estimated by fusing model-based predictions with measurements taken from a sparse network of monitors in the region. As one of the aims of the study is to produce a dose-response curve, a regression model relating perceptual response to noise exposure levels, it is important to acknowledge the potential attenuation bias that results from measurement error in the estimated noise exposure levels. In this presentation we review and compare several methods for dealing with measurement error in generalized linear mixed models. The methods are demonstrated on simulated data and real data collected during past NASA risk reduction studies.

measurement error↗
DOE OSTI2018

Ecological analysis of associations between groundwater quality and hypertension and cardiovascular disease in rural Saskatchewan, Canada using Bayesian hierarchical models and administrative health data

Highlights: • Hardness and magnesium were associated with decreased risk of cardiovascular disease. • Arsenic concentrations were not associated with the risk of cardiovascular disease. • Summarized water health standards were not associated with cardiovascular disease. Associations between groundwater quality and the prevalence of hypertension and ischemic heart disease were investigated in rural areas of the Canadian province of Saskatchewan. The partially ecological study was analyzed using Bayesian hierarchical models to account for spatial variability in risk. Exposure measures and health outcomes were estimated based on previously collected water quality surveillance data from public water supplies and private wells and administrative health data. Water quality exposures for each study region were estimated by applying geostatistical techniques to arsenic concentrations and principal component (PC) scores. The PC scores summarized groups of parameters measuring either health standards or aesthetic objectives described by the province. Generalized linear mixed models with a log link assessed associations between water quality and observed count of health outcomes relative to the expected value. The Bayesian models contained uncorrelated and spatially correlated random effects for each geographic region. Effect estimates were controlled for sex and age by stratification of case and expected case counts, for smoking by inclusion of sex- and age-specific prevalence of chronic obstructive pulmonary disease as a surrogate covariate, and for education and income by use of census data. There was no evidence for associations between groundwater arsenic concentrations in public or private water supplies and increased risk of hypertension or cardiovascular disease. An increase in the second aesthetic objectives PC score from public supplies was associated with a protective effect against ischemic heart disease. This PC value summarized hardness and magnesium Similarly, an increase in the second aesthetic objectives PC in private supplies was associated with decreased prevalence of hypertension. The results of this study are consistent with others demonstrating a relationship between elevated hardness and magnesium concentrations in drinking water and reduced risk of hypertension and cardiovascular disease. Further investigation is warranted with individual exposure history, particularly with respect to the potential beneficial effect of hard water on the prevalence of hypertension and cardiovascular disease.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Organic geochemistry of sediments from the continental margin off southern New England, U.S.A.--Part I. Amino acids, carbohydrates and lignin

Total organic carbon (TOC), lignin, amino acids, sugars and amino sugars were measured in recent sediments for the continental margin off southern New England. The various organic carbon fractions decreased in concentration with increasing distance from shore. The fraction of the TOC that was accounted for by these major components also decreased with increasing distance from shore. The concentration of lignin indicated that only about 3-5% of the organic carbon in the nearshore sediment was of terrestrial origin. The various fractions were highly correlated, which was consistent with a simple linear mixing model of shelf organic matter with material form the slope and rise and indicated a significant transport of sediment from the continental shelf to the continental slope and rise.

Non-NASA Center↗
DOE OSTI2019

Early risk stratification using Rubidium-82 positron emission tomography in STEMI patients

Background: Assessment of infarct size after myocardial infarction is predictive of subsequent morphological changes and clinical outcome. This study aimed to assess subacute post-intervention Rubidium-82 ({sup 82}Rb)-PET imaging in predicting left ventricle ejection fraction, regional wall motion, and final infarct size by CMR at 3-months after STEMI. Methods: STEMI patients undergoing percutaneous coronary intervention were included prospectively. Rest-only {sup 82}Rb-PET perfusion imaging was performed at median 36 hours [IQR: 22 to 50] after the treatment. The extent of hypoperfusion and absolute blood flow (mL·min·g) were estimated on a global and a 17-segment model with dedicated software. At 3-months follow-up patients completed the CMR functional and late gadolinium enhancement imaging. Results: 42 patients were included, but only 35 had follow-up CMR and constituted the study population. Absolute blood flow was significantly lower in the infarct-related territory compared to remote myocardium, P < .005. Extent of hypoperfusion correlated with final infarct size, r = 0.58, P < .001, while blood flow correlated with ejection fraction, r = 0.41, P < .05. In linear mixed models, higher subacute absolute blood flow (β = 4.6, confidence interval [3.5; 5.2], P < .001, R{sup 2} = 0.67) was associated with greater wall motion. Segmental extent of subacute hypoperfusion (β = 0.43 [0.38; 0.49], P < .001, R{sup 2} = 0.58) was associated with the degree of late gadolinium enhancement at 3-months. Conclusions: Subacute rest-only {sup 82}Rb-PET is feasible following STEMI and seems predictive of myocardial function and infarct size at 3-months.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Modeling Pancreatic Tumor Motion Using 4-Dimensional Computed Tomography and Surrogate Markers

Purpose: To assess intrafractional positional variations of pancreatic tumors using 4-dimensional computed tomography (4D-CT), their impact on gross tumor volume (GTV) coverage, the reliability of biliary stent, fiducial seeds, and the real-time position management (RPM) external marker as tumor surrogates for setup of respiratory gated treatment, and to build a correlative model of tumor motion. Methods and Materials: We analyzed the respiration-correlated 4D-CT images acquired during simulation of 36 patients with either a biliary stent (n=16) or implanted fiducials (n=20) who were treated with RPM respiratory gated intensity modulated radiation therapy for locally advanced pancreatic cancer. Respiratory displacement relative to end-exhalation was measured for the GTV, the biliary stent, or fiducial seeds, and the RPM marker. The results were compared between the full respiratory cycle and the gating interval. Linear mixed model was used to assess the correlation of GTV motion with the potential surrogate markers. Results: The average ± SD GTV excursions were 0.3 ± 0.2 cm in the left-right direction, 0.6 ± 0.3 cm in the anterior-posterior direction, and 1.3 ± 0.7 cm in the superior-inferior direction. Gating around end-exhalation reduced GTV motion by 46% to 60%. D95% was at least the prescribed 56 Gy in 76% of patients. GTV displacement was associated with the RPM marker, the biliary stent, and the fiducial seeds. The correlation was better with fiducial seeds and with biliary stent. Conclusions: Respiratory gating reduced the margin necessary for radiation therapy for pancreatic tumors. GTV motion was well correlated with biliary stent or fiducial seed displacements, validating their use as surrogates for daily assessment of GTV position during treatment. A patient-specific internal target volume based on 4D-CT is recommended both for gated and not-gated treatment; otherwise, our model can be used to predict the degree of GTV motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRSDate not supplied

Teal Cycle Ergometer Vibration Isolation and Stabilization System Human-in-the-Loop Ground Evaluation for Long-Duration Spaceflight

Purpose: To counter the deleterious effects of weightlessness on the cardiopulmonary system, astronauts living on the International Space Station exercise on a variety of countermeasure equipment including the Cycle Ergometer Vibration Isolation and Stabilization System (CEVIS). Operational since 2001, the onboard CEVIS will be replaced by a new model, known as Teal CEVIS (TC). As a part of ground evaluation, TC hardware underwent human-in-the-loop (HITL) testing to verify the TC hardware produces workloads that elicit physiologic responses comparable to a laboratory cycle ergometer (LAB). Methods: Seven subjects (5 M/2 F) performed submaximal cycle ergometer testing with indirect calorimetry measures on TC and LAB on separate test days. Testing consisted of graded 30 watt increases in workload until subjects reached 85% of age-predicted max heart rate (HR). Exercise outcomes included rate of oxygen uptake (VO2; liters/min), rate of energy expenditure (REE; kcal/min), and HR (beats/min). Linear mixed models (LMM) were fitted to compare VO2, REE, and HR responses between devices across power outputs with fixed effects for power (P) and device (D) and with random effects for subject. LMM effect coefficients (β), std errors (SE), pseudo-partial R2 of effects (pR2), and model likelihood ratio statistics (χ2, Pr(> chisq); α<.05) are provided. Results: LMM main effects for P and D were observed for VO2 (βp=.0107, SE = .0002, pR2=0.97; βd = -0.075, SE = 0.019, pR2=0.14; χ2(1)= 13.57, p<.001), such that VO2 was higher across stages on TC. Main effects of P and D were observed for REE (βp= .0586, SE = 0.001, pR2= 0.97; βd = -0.33208, SE = 0.099, pR2= 0.056; χ2(1)=10.481, p < .01), such that REE was higher across power outputs on TC. A P x D interaction was observed for HR, along with main effects for P and D (βp x d =-0.048, SE = 0.02, pR2=0.048; βp= 0.42, SE = 0.016, pR2=0.87; χ2(1) = 5.398, p = 0.02), such that higher workloads elicited a greater difference in HR between devices. Conclusions: HITL results show, TC elicits greater physiologic responses across power outputs compared to LAB. However, pR2 for device effects show small differences between devices. Therefore, TC can be expected to provide appropriate physiological stimulus across workloads and be considered a reliable tool to mitigate the effects of weightlessness.

Cycle Ergometer Vibration and Isolation System↗
DOE OSTI2021

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗