SEARCH · Engineering Papers

Results for “maximum likelihood estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

569 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Predicting Maximum Temperatures of a Li-ion Battery on a Simulated Flight Profile using a Model-based Prognostics

One of the challenges in using Li-ion packs in aeronautics is their safety, and thermal runaway (TR) is a major concern. The current engineering solutions to prevent a Li-ion pack from a catastrophic TR require additional mass and volume to isolate cells. The excess mass could be reduced by improving detection and, thus, preventing a TR event. One of the possible early warning indicators of a TR is crossing a threshold temperature. We have developed an approach, based on the Unscented Kalman Filter (UKF), to predict the likelihood of reaching the threshold temperature for simulated flight profiles. The current battery prognostics algorithms for aerospace predict state-of-charge (SOC) and end-of-life (EOL) [1]. We extended this two-level algorithm to predict the maximum temperature during discharge. The amount of heat generated in a cell depends on factors such as cell chemistry, cell packaging, total cycles, operating temperature, and abuse history [2]. Our semi-empirical thermal model depends on three phenomenological parameters which account for those factors. In addition, a two-parameter reduced-order model is developed to predict the temperature rise for short bursts of “random-walk” (RW) discharge current sequence, which simulates a flight's current-loading profile. The performance of these models on different datasets and types of current loading will be presented. To predict the maximal temperatures for future cycles we must estimate the evolution of thermal parameters as the batteries age. It is found that the parameters of the 3-parametric thermal model cannot be estimated only from the RW data. To address the issue, we will present two alternative approaches: i) expanding the datasets to include discharge profiles beyond RWs; ii) model reduction to a two-parametric model. The two approaches will be illustrated by an application to the cycling data from a commercial LG 18650 cell. References: 1. M. Daigle, C.S. Kulkarni, End-of-discharge and End-of-life Prediction in Lithium-ion Batteries with Electrochemistry-based Aging Models, in: AIAA Infotech @ Aerospace, American Institute of Aeronautics and Astronautics, San Diego, California, USA, 2016. 2. M. Börner, et. al, Correlation of aging and thermal stability of commercial 18650-type lithium ion batteries, Journal of Power Sources. 342 (2017) 382–392.

Thermal runaway↗

NASA NTRS2017

Determination of the Risk of Radiation-Associated Circulatory and Cancer Disease Mortality in a NASA Early Astronaut Cohort

Of the many possible health challenges posed during extended exploratory missions to space, the effects of space radiation on cardiovascular disease and cancer are of particular concern. There are unique challenges to estimating those radiation risks; care and appropriate and rigorous methodology should be applied when considering small cohorts such as the NASA astronaut population. The objective of this work was to determine if there was sufficient evidence for excess risk of cardiovascular disease and cancer in early NASA astronaut cohorts. NASA astronauts in selection groups 1-7 were chosen; this relatively homogeneous cohort consists of 73 white males, who unlike today's astronauts, maintained similar smoking and drinking habits to the general US population, and have published radiation doses. The participants flew in space on missions Mercury through Shuttle and received space radiation doses between 0-74.1 milligrays. Cause of death information was obtained from the Lifetime Surveillance of Astronaut Health (LSAH) program at NASA Johnson Space Center. Mortality was compared with the US male population. Trends of mortality with dose were assessed using a logistic model, fitted by maximum likelihood. Only 32 (43.84 percent) of the 73 early astronauts have died. Standard mortality ratios (SMRs) for cancer (n=7, SMR=43.4, 95 percent CI 17.8, 84.9), all circulatory disease (n=7, SMR=33.2, 95 percent CI 13.7, 65.0), and ischemic heart disease (IHD) (n=5, SMR=40.1, 95 percent CI 13.2, 89.4) were significantly lower than for the US white male population. For cerebrovascular disease, the upper confidence interval for SMR included 100, indicating it was not significantly different from the US population (n=2, SMR = 77.0, 95 percent CI 9.4, 268.2). The power of the study is low and remains below 10 percent even when risks 10 times those reported in the literature are assumed. Due to small sample size, there is currently insufficient statistical power to evaluate space radiation exposure effects on mortality in NASA astronauts. In addition to a comprehensive longitudinal study of NASA astronauts, a research strategy of low dose epidemiology data integration with cell and animal studies should be utilized for space radiation risk assessment in the astronauts.

Elgart, S. R.↗

NASA NTRS2014

Mass and Volume Optimization of Space Flight Medical Kits

Resource allocation is a critical aspect of space mission planning. All resources, including medical resources, are subject to a number of mission constraints such a maximum mass and volume. However, unlike many resources, there is often limited understanding in how to optimize medical resources for a mission. The Integrated Medical Model (IMM) is a probabilistic model that estimates medical event occurrences and mission outcomes for different mission profiles. IMM simulates outcomes and describes the impact of medical events in terms of lost crew time, medical resource usage, and the potential for medically required evacuation. Previously published work describes an approach that uses the IMM to generate optimized medical kits that maximize benefit to the crew subject to mass and volume constraints. We improve upon the results obtained previously and extend our approach to minimize mass and volume while meeting some benefit threshold. METHODS We frame the medical kit optimization problem as a modified knapsack problem and implement an algorithm utilizing dynamic programming. Using this algorithm, optimized medical kits were generated for 3 mission scenarios with the goal of minimizing the medical kit mass and volume for a specified likelihood of evacuation or Crew Health Index (CHI) threshold. The algorithm was expanded to generate medical kits that maximize likelihood of evacuation or CHI subject to mass and volume constraints. RESULTS AND CONCLUSIONS In maximizing benefit to crew health subject to certain constraints, our algorithm generates medical kits that more closely resemble the unlimited-resource scenario than previous approaches which leverage medical risk information generated by the IMM. Our work here demonstrates that this algorithm provides an efficient and effective means to objectively allocate medical resources for spaceflight missions and provides an effective means of addressing tradeoffs in medical resource allocations and crew mission success parameters.

Monte Carlo Method↗

NASA NTRS2001

Evidence-Based Approach to the Analysis of Serious Decompression Sickness with Application to EVA Astronauts

It is important to understand the risk of serious hypobaric decompression sickness (DCS) in order to develop procedures and treatment responses to mitigate the risk. Since it is not ethical to conduct prospective tests about serious DCS with humans, the necessary information was gathered from 73 published reports. We hypothesize that a 4-hr 100% oxygen (O2) prebreathe results in a very low risk of serious DCS, and test this through analysis. We evaluated 258 tests containing information from 79,366 exposures in attitude chambers. Serious DCS was documented in 918 men during the tests. Serious DCS are signs and symptoms broadly classified as Type II DCS. A risk function analysis with maximum likelihood optimization was performed to identify significant explanatory variables, and to create a predictive model for the probability of serious DCS [P(serious DCS)]. Useful variables were Tissue Ratio, the planned time spent at altitude (T(sub alt)), and whether or not repetitive exercise was performed at altitude. Tissue Ratio is P1N2/P2, where P1N2 is calculated nitrogen (N2) pressure in a compartment with a 180-min half-time for N2 pressure just before ascent, and P2 is ambient pressure after ascent. A prebreathe and decompression profile Shuttle astronauts use for extravehicular activity (EVA) includes a 4-hr prebreathe with 100% O2, an ascent to P2 = 4.3 lb per sq. in. absolute, and a T(sub alt) = 6 hr. The P(serious DCS) is: 0.0014 (0.00096 - 0.00196, 95% confidence interval) with exercise and 0.00025 (0.00016 - 0.00035) without exercise. Given 100 Shuttle EVAs to date and no report of serious DCS, the true risk is less than 0.03 with 95% confidence (Binomial Theorem). It is problematic to estimate the risk of serious DCS since it appears infrequently, even if the estimate is based on thousands of altitude chamber exposures. The true risk to astronauts may lie between the extremes of the confidence intervals (0.00016 - 0.00196) since the contribution of other factors, particularly exercise, to the risk of serious DCS during EVA is unknown. A simple model that only accounts for four important variables in retrospective data is still helpful to increase our understanding about the risk of serious DCS.

Conkin, Johnny↗

NASA NTRS2012

Climate Change: A New Metric to Measure Changes in the Frequency of Extreme Temperatures using Record Data

Consensus on global warming is the result of multiple and varying lines of evidence, and one key ramification is the increase in frequency of extreme climate events including record high temperatures. Here we develop a metric- called "record equivalent draws" (RED)-based on record high (low) temperature observations, and show that changes in RED approximate changes in the likelihood of extreme high (low) temperatures. Since we also show that this metric is independent of the specifics of the underlying temperature distributions, RED estimates can be aggregated across different climates to provide a genuinely global assessment of climate change. Using data on monthly average temperatures across the global landmass we find that the frequency of extreme high temperatures increased 10-fold between the first three decades of the last century (1900-1929) and the most recent decade (1999-2008). A more disaggregated analysis shows that the increase in frequency of extreme high temperatures is greater in the tropics than in higher latitudes, a pattern that is not indicated by changes in mean temperature. Our RED estimates also suggest concurrent increases in the frequency of both extreme high and extreme low temperatures during 2002-2008, a period when we observe a plateauing of global mean temperature. Using daily extreme temperature observations, we find that the frequency of extreme high temperatures is greater in the daily minimum temperature time-series compared to the daily maximum temperature time-series. There is no such observable difference in the frequency of extreme low temperatures between the daily minimum and daily maximum.

Munasinghe, L.↗

NASA NTRS2008

Application of an Ensemble Smoother to Precipitation Assimilation

Assimilation of precipitation in a global modeling system poses a special challenge in that the observation operators for precipitation processes are highly nonlinear. In the variational approach, substantial development work and model simplifications are required to include precipitation-related physical processes in the tangent linear model and its adjoint. An ensemble based data assimilation algorithm "Maximum Likelihood Ensemble Smoother (MLES)" has been developed to explore the ensemble representation of the precipitation observation operator with nonlinear convection and large-scale moist physics. An ensemble assimilation system based on the NASA GEOS-5 GCM has been constructed to assimilate satellite precipitation data within the MLES framework. The configuration of the smoother takes the time dimension into account for the relationship between state variables and observable rainfall. The full nonlinear forward model ensembles are used to represent components involving the observation operator and its transpose. Several assimilation experiments using satellite precipitation observations have been carried out to investigate the effectiveness of the ensemble representation of the nonlinear observation operator and the data impact of assimilating rain retrievals from the TMI and SSM/I sensors. Preliminary results show that this ensemble assimilation approach is capable of extracting information from nonlinear observations to improve the analysis and forecast if ensemble size is adequate, and a suitable localization scheme is applied. In addition to a dynamically consistent precipitation analysis, the assimilation system produces a statistical estimate of the analysis uncertainty.

Zhang, Sara↗

NASA NTRS2023

Using Regionalized Air Quality Model Performance and Bayesian Maximum Entropy Data Fusion to Map Global Surface Ozone Concentration and Associated Uncertainty

Estimates of ground-level ozone concentrations have been improved through data fusion of observations and atmospheric chemistry models. Our previous global ozone estimates for the Global Burden of Disease study corrected for bias uniformly across continents and then corrected near monitoring stations using the Bayesian Maximum Entropy (BME) framework for data fusion. Here, we use the Regionalized Air Quality Model Performance (RAMP) framework to correct model bias over a much larger spatial range than BME can, accounting for the spatial inhomogeneity of bias and nonlinearity as a function of modeled ozone. RAMP bias correction is applied to a composite of 9 global chemistry-climate models, based on the nearest set of monitors. These estimates are then fused with observations using BME, which matches observations at measurement stations, with the influence of observations declining with distance in space and time. We create global ozone maps for each year from 1990 to 2017 at fine spatial resolution. RAMP is shown to create unrealistic discontinuities due to the spatial clustering of ozone monitors, which we overcome by applying a weighting for RAMP based on the number of monitors nearby. Incorporating RAMP before BME has little effect on model performance near stations, but strongly increases R 2 by 0.15 at locations farther from stations, shown through a checkerboard cross-validation. Corrections to estimates differ based on location in space and time, confirming heterogeneity. We quantify the likelihood of exceeding selected ozone levels, finding that parts of the Middle East, India, and China are most likely to exceed 55 parts per billion (ppb) in 2017. About 96% of the global population was exposed to ozone levels above the World Health Organization guideline of 60 µg m −3 (30 ppb) in 2017. Our annual fine-resolution ozone estimates may be useful for several applications including epidemiology and assessments of impacts on health, agriculture, and ecosystems.

Ozone↗

NASA NTRS1996

Early Energetic Particle Irradiation of the HED Parent Body Regolith

Previous studies have shown that many individual grains within the dark phase of the Kapoeta howardite were irradiated with energetic particles while residing on the surface of the early HED regolith. Particle tracks in these grains vary in density by more than an order of magnitude and undoubtedly were formed by energetic heavy (Fe) ions associated with early solar flares. Early Irradiation of HED Regolith: Concentrations of excess Ne alone are not sufficient to decide between competing galactic and solar irradiation models. However, from recent studies of depth samples of oriented lunar rocks, we have shown that the cosmogenic 21-Ne/22-Ne ratio produced in feldspar differs substantially between Galactic Cosmic Radiation (GCR) and solar protons, and that this difference is exactly that predicted from cross-section data. Using Ne literature data and new isotopic data we obtained on acid-etched, separated feldspar from both the light and dark phases of Kapoeta, we derive 21-Ne/22-Ne = 0.80 for the recent GCR irradiation and 21-Ne/22-Ne = 0.68 for the early regolith irradiation. This derived ratio indicates that the early Ne production in the regolith occurred by both galactic and solar protons. If we adopt a likely one-component regolith model in which all grains were exposed to galactic protons but individual grains had variable exposure to solar protons, we estimate that this early GCR irradiation lasted for about 3-6 m.y. More complex two-component regolith models involving separate solar and galactic irradiation would permit this GCR age to be longer. Higher-energy solar protons would permit the GCR to be longer. Higher-energy solar protons would permit the GCR age to be shorter. Further, cosmogenic 126(Xe) in Kapoeta dark is no more than a factor of about 2 higher than that observed in Kapoeta light. Because 126(Xe) can only be formed by galactic protons and not solar protons, these data support a short GCR irradiation for the HED regolith. This would also be the maximum time peRiod for the solar irradiation. Various asteroidal regolith models, based on Monte Carlo modeling of impact rates as a function of size and on irradiation features of meteorites, predict surface exposure times of about 0.1 to 10 m.y., and depend on such factors as gravity, rock mechanical properties, and micrometeoroid flux. Because the depth at which solar Fe tracks are produced (is much less than 1 micrometer) is much less than the depth at which Solar Cosmic Rays (SCR) Ne is produced (about 1 cm), for a reasonably well-stirred HED regolith the "surface exposure time" for SCR 21-Ne production should be significantly longer than that for solar tracks and some other surface irradiation features. Enhanced Solar Proton Irradiation: For bulk samples of Kapoeta dark feldspar and a one-component regolith model, the derived ratio of 21-Ne/22-Ne = 0.68 implies that the early production ratio of SCR 21-Ne to GCR 21-Ne was about 0.5-1.5. This ratio is independent of any assumptions about the fraction of dark grains that are irradiated or of the variability in the degree of solar irradiation among grains. The 21-Ne SCR/GCR ratio indirectly derived from bulk Kapoeta pyroxene is somewhat larger, as is the ratio derived for simple two-component regolith models. Individual feldspar grains that were extensively solar irradiated would require even larger 21-Ne SCR/GCR production ratios. In contrast, the theoretical SCR/GCR production ratio for lunar feldspar with 0 g/CM2 shield ing is is less than or equal to 2, and the lowest ratio observed in near-surface samples of lunar anorthosites is less than or equal to 1. Considering the greater solar distance of Vesta (compared to the Moon), the likelihood that SCR 21-Ne was acquired under some shielding where production rates are lower, and the likelihood that the exposure time to galactic protons exceeded the exposure time to solar protons because of their very different penetration depths, the 21-Ne SCR/GCR production ratio on the HED parent body was probably < 0.1. The relatively large difference between the derived 21-Ne SCR/GCR ratio in Kapoeta dark feldspar and the estimated production ratio strongly indicates that the early solar irradiation involved a flux -20-50x the recent solar flux. This enhanced proton flux was probably associated with an overall greater solar activity in the first approximately 10(exp 7) to 10(exp 8) years of solar history.

Bogard, D. D.↗

NASA NTRS1995

On the uncertainty in single molecule fluorescent lifetime and energy emission measurements

Time-correlated single photon counting has recently been combined with mode-locked picosecond pulsed excitation to measure the fluorescent lifetimes and energy emissions of single molecules in a flow stream. Maximum likelihood (ML) and least square methods agree and are optimal when the number of detected photons is large however, in single molecule fluorescence experiments the number of detected photons can be less than 20, 67% of those can be noise and the detection time is restricted to 10 nanoseconds. Under the assumption that the photon signal and background noise are two independent inhomogeneous poisson processes, we derive the exact joint arrival time probably density of the photons collected in a single counting experiment performed in the presence of background noise. The model obviates the need to bin experimental data for analysis, and makes it possible to analyze formally the effect of background noise on the photon detection experiment using both ML or Bayesian methods. For both methods we derive the joint and marginal probability densities of the fluorescent lifetime and fluorescent emission. the ML and Bayesian methods are compared in an analysis of simulated single molecule fluorescence experiments of Rhodamine 110 using different combinations of expected background nose and expected fluorescence emission. While both the ML or Bayesian procedures perform well for analyzing fluorescence emissions, the Bayesian methods provide more realistic measures of uncertainty in the fluorescent lifetimes. The Bayesian methods would be especially useful for measuring uncertainty in fluorescent lifetime estimates in current single molecule flow stream experiments where the expected fluorescence emission is low. Both the ML and Bayesian algorithms can be automated for applications in molecular biology.

Brown, Emery N.↗

NASA NTRS1996

On the Uncertainty in Single Molecule Fluorescent Lifetime and Energy Emission Measurements

Time-correlated single photon counting has recently been combined with mode-locked picosecond pulsed excitation to measure the fluorescent lifetimes and energy emissions of single molecules in a flow stream. Maximum likelihood (ML) and least squares methods agree and are optimal when the number of detected photons is large, however, in single molecule fluorescence experiments the number of detected photons can be less than 20, 67 percent of those can be noise, and the detection time is restricted to 10 nanoseconds. Under the assumption that the photon signal and background noise are two independent inhomogeneous Poisson processes, we derive the exact joint arrival time probability density of the photons collected in a single counting experiment performed in the presence of background noise. The model obviates the need to bin experimental data for analysis, and makes it possible to analyze formally the effect of background noise on the photon detection experiment using both ML or Bayesian methods. For both methods we derive the joint and marginal probability densities of the fluorescent lifetime and fluorescent emission. The ML and Bayesian methods are compared in an analysis of simulated single molecule fluorescence experiments of Rhodamine 110 using different combinations of expected background noise and expected fluorescence emission. While both the ML or Bayesian procedures perform well for analyzing fluorescence emissions, the Bayesian methods provide more realistic measures of uncertainty in the fluorescent lifetimes. The Bayesian methods would be especially useful for measuring uncertainty in fluorescent lifetime estimates in current single molecule flow stream experiments where the expected fluorescence emission is low. Both the ML and Bayesian algorithms can be automated for applications in molecular biology.

Brown, Emery N.↗

NASA NTRS2003

Bayesian estimation of dose thresholds

An example is described of Bayesian estimation of radiation absorbed dose thresholds (subsequently simply referred to as dose thresholds) using a specific parametric model applied to a data set on mice exposed to 60Co gamma rays and fission neutrons. A Weibull based relative risk model with a dose threshold parameter was used to analyse, as an example, lung cancer mortality and determine the posterior density for the threshold dose after single exposures to 60Co gamma rays or fission neutrons from the JANUS reactor at Argonne National Laboratory. The data consisted of survival, censoring times and cause of death information for male B6CF1 unexposed and exposed mice. The 60Co gamma whole-body doses for the two exposed groups were 0.86 and 1.37 Gy. The neutron whole-body doses were 0.19 and 0.38 Gy. Marginal posterior densities for the dose thresholds for neutron and gamma radiation were calculated with numerical integration and found to have quite different shapes. The density of the threshold for 60Co is unimodal with a mode at about 0.50 Gy. The threshold density for fission neutrons declines monotonically from a maximum value at zero with increasing doses. The posterior densities for all other parameters were similar for the two radiation types.

NASA Discipline Radiation Health↗