Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 271 records · Page 15

Code for paper "Reconstructing an Implosion with Space-Time Smoothing of Velocimetry Traces"

New diagnostics to obtain high-temporal resolution observation of surface velocity provide a rich collection of data from high-explosive (HE) experiments. Using this data to estimate the temporal evolution of an imploding, HE-driven surface is of great interest for validating simulation of the physical system and for design of future experiments. This is the code for a paper in which we propose a method for using a Gaussian process to estimate the surface position over time and provide uncertainty estimates on that reconstruction accounting for measurement error, uncertainty to sparse spatial sampling velocity data, and robustness to outlying observations. In the paper, the code is applied on data from a hemispherical HE implosion experiment. This code is for release with the paper for peer-review, replicating the results and making the figures from the paper.

Grosskopf, Michael↗

Detailed Uncertainty Analysis of the Ares I A106 Liftoff/Transition Database

The Ares I A106 Liftoff/Transition Force and Moment Aerodynamics Database describes the aerodynamics of the Ares I Crew Launch Vehicle (CLV) from the moment of liftoff through the transition from high to low total angles of attack at low subsonic Mach numbers. The database includes uncertainty estimates that were developed using a detailed uncertainty quantification procedure. The Ares I Aerodynamics Panel developed both the database and the uncertainties from wind tunnel test data acquired in the NASA Langley Research Center s 14- by 22-Foot Subsonic Wind Tunnel Test 591 using a 1.75 percent scale model of the Ares I and the tower assembly. The uncertainty modeling contains three primary uncertainty sources: experimental uncertainty, database modeling uncertainty, and database query interpolation uncertainty. The final database and uncertainty model represent a significant improvement in the quality of the aerodynamic predictions for this regime of flight over the estimates previously used by the Ares Project. The maximum possible aerodynamic force pushing the vehicle towards the launch tower assembly in a dispersed case using this database saw a 40 percent reduction from the worst-case scenario in previously released data for Ares I.

Hanke, Jeremy L.↗

Updates to the Regional Seismic Travel Time (RSTT) Model: 2. Path-dependent Travel-time Uncertainty

Abstract The regional seismic travel time (RSTT) model and software were developed to improve travel-time prediction accuracy by accounting for three-dimensional crust and upper mantle structure. Travel-time uncertainty estimates are used in the process of associating seismic phases to events and to accurately calculate location uncertainty bounds (i.e. event location error ellipses). We improve on the current distance-dependent uncertainty parameterization for RSTT using a random effects model to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The random effects model separates the error between observed slowness and model predicted slowness into bias and random components. The path-specific travel-time uncertainty is calculated by integrating these mean squared errors along a seismic-phase ray path. We demonstrate that event location error ellipses computed for a 90% coverage ellipse metric (used by the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (IDC)), and using the path-specific travel-time uncertainty approach, are more representative (median 82.5% ellipse percentage) of true location error than error ellipses computed using distance-dependent travel-time uncertainties (median 70.1%). We also demonstrate measurable improvement in location uncertainties using the RSTT method compared to the current station correction approach used at the IDC (median 74.3% coverage ellipse).

58 GEOSCIENCES↗

1,2-Dichlorohexafluoro-Cyclobutane (1,2-c-C4F6Cl2, R-316c) a Potent Ozone Depleting Substance and Greenhouse Gas: Atmospheric Loss Processes, Lifetimes, and Ozone Depletion and Global Warming Potentials for the (E) and (Z) stereoisomers

The atmospheric processing of (E)- and (Z)-1,2-dichlorohexafluorocyclobutane (1,2-c-C4F6Cl2, R-316c) was examined in this work as the ozone depleting (ODP) and global warming (GWP) potentials of this proposed replacement compound are presently unknown. The predominant atmospheric loss processes and infrared absorption spectra of the R-316c isomers were measured to provide a basis to evaluate their atmospheric lifetimes and, thus, ODPs and GWPs. UV absorption spectra were measured between 184.95 to 230 nm at temperatures between 214 and 296 K and a parametrization for use in atmospheric modeling is presented. The Cl atom quantum yield in the 193 nm photolysis of R- 316c was measured to be 1.90 +/- 0.27. Hexafluorocyclobutene (c-C4F6) was determined to be a photolysis co-product with molar yields of 0.7 and 1.0 (+/-10%) for (E)- and (Z)-R-316c, respectively. The 296 K total rate coefficient for the O(1D) + R-316c reaction, i.e., O(1D) loss, was measured to be (1.56 +/- 0.11) × 10(exp −10)cu cm/ molecule/s and the reactive rate coefficient, i.e., R-316c loss, was measured to be (1.36 +/- 0.20) × 10(exp −10)cu cm/molecule/s corresponding to a approx. 88% reactive yield. Rate coefficient upper-limits for the OH and O3 reaction with R-316c were determined to be <2.3 × 10(exp −17) and <2.0 × 10(exp −22)cu cm/molecule/s, respectively, at 296 K. The quoted uncertainty limits are 2(sigma) and include estimated systematic errors. Local and global annually averaged lifetimes for the (E)- and (Z)-R-316c isomers were calculated using a 2-D atmospheric model to be 74.6 +/- 3 and 114.1 +/-10 years, respectively, where the estimated uncertainties are due solely to the uncertainty in the UV absorption spectra. Stratospheric photolysis is the predominant atmospheric loss process for both isomers with the O(1D) reaction making a minor, approx. 2% for the (E) isomer and 7% for the (Z) isomer, contribution to the total atmospheric loss. Ozone depletion potentials for (E)- and (Z)-R-316c were calculated using the 2-D model to be 0.46 and 0.54, respectively. Infrared absorption spectra for (E)- and (Z)-R-316c were measured at 296 K and used to estimate their radiative efficiencies (REs) and GWPs; 100-year time-horizon GWPs of 4160 and 5400 were obtained for (E)- and (Z)-R-316c, respectively. Both isomers of R-316c are shown in this work to be long-lived ozone depleting substances and potent greenhouse gases.

UV cross section↗

Construction of Polarimetric Radar-Based Reference Rain Maps for the Iowa Flood Studies Campaign

The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.

Radar↗

Evaluation of cloud height, optical thickness, and phase retrievals from the CHROMA algorithm applied to Sentinel-3 OLCI data

We previously developed the Cloud Height Retrieval from O 2 Molecular Absorption (CHROMA) algorithm for the Ocean Color Instrument (OCI) on the new NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission. Here, we apply CHROMA to observations from the Ocean Land Colour Instrument (OLCI) to guide expectations for PACE, as it will take some time to obtain large-scale validation data for OCI. We use cloud top height (CTH), phase, and (for liquid clouds) cloud optical thickness (COT) data from the ground-based Atmospheric Radiation Measurement (ARM) network to evaluate the OLCI retrievals. We found that OLCI and Moderate Resolution Imaging Spectroradiometer (MODIS) CTH compare similarly well to the ARM reference. OLCI has a tendency to underestimate CTH as CTH increases, and algorithm assumptions about cloud geometric thickness may contribute to this. ARM COT from multifilter shadowband radiometers (MFRSR) and Sun photometers are well-correlated with one another, albeit with a roughly 30 % offset on average; OLCI and MODIS COT agree more closely with the MFRSR data. OLCI retrieval uncertainty estimates show skill at telling low-uncertainty cases from high-uncertainty ones, although CTH uncertainties are underestimated. Additionally, we compare the OLCI data to satellite retrievals based on thermal infrared measurements from MODIS and Sea and Land Surface Temperature Radiometer (SLSTR) data. Differences are broadly consistent with physical expectations based on the A-band vs. thermal techniques, although one key challenge in such aggregated comparisons is different cloud masking sensitivities and algorithm failure rates meaning additional sampling differences are introduced. We conclude by discussing the transition to and possible enhancements for PACE OCI.

Sayer, Andrew M. [Univ. of Maryland Baltimore Coun↗

Evidence of Temporal Variation of Titan Atmospheric Density in 2005-2013

One major science objective of the Cassini mission is an investigation of Titan's atmosphere constituent abundances. Titan's atmospheric density is of interest not only to planetary scientists but also to mission design and mission control engineers. Knowledge of the dependency of Titan's atmospheric density with altitude is important because any unexpectedly high atmospheric density has the potential to tumble the spacecraft during a flyby. During low-altitude Titan flyby, thrusters are fired to counter the torque imparted on the spacecraft due to the Titan atmosphere. The denser the Titan's atmosphere is, the higher are the duty cycles of the thruster firings. Therefore thruster firing telemetry data could be used to estimate the atmospheric torque imparted on the spacecraft. Since the atmospheric torque imparted on the spacecraft is related to the Titan's atmospheric density, atmospheric densities are estimated accordingly. In 2005-2013, forty-three low-altitude Titan flybys were executed. The closest approach altitudes of these Titan flybys ranged from 878 to 1,074.8 km. Our density results are also compared with those reported by other investigation teams: Voyager-1 (in November 1980) and the Huygens Atmospheric Structure Instrument, HASI (in January 2005). From our results, we observe a temporal variation of the Titan atmospheric density in 2005-2013. The observed temporal variation is significant and it isn't due to the estimation uncertainty (5.8%, 1 sigma) of the density estimation methodology. Factors that contributed to this temporal variation have been conjectured but are largely unknown. The observed temporal variation will require synergetic analysis with measurements made by other Cassini science instruments and future years of laboratory and modeling efforts to solve. The estimated atmospheric density results are given in this paper help scientists to better understand and model the density structure of the Titan atmosphere.

Atmosphere↗

On the representation and estimation of spatial uncertainty

This paper describes a general method for estimating the nominal relationship and expected error (covariance) between coordinate frames representing the relative locations of objects. The frames may be known only indirectly through a series of spatial relationships, each with its associated error, arising from diverse causes, including positioning errors, measurement errors, or tolerances in part dimensions. This estimation method can be used to answer such questions as whether a camera attached to a robot is likely to have a particular reference object in its field of view. The calculated estimates agree well with those from an independent Monte Carlo simulation. The method makes it possible to decide in advance whether an uncertain relationship is known accurately enough for some task and, if not, how much of an improvement in locational knowledge a proposed sensor will provide. The method presented can be generalized to six degrees of freedom and provides a practical means of estimating the relationships (position and orientation) among objects, as well as estimating the uncertainty associated with the relationships.

Smith, Randall C.↗

Validation of Regional CO2 Concentrations in the ECMWF Real-Time Analysis and Carbon-Tracker Reanalysis with Airborne Observations from ACT-A Field Campaign

Through verifying against hundreds of hours of airborne in-situ measurements from the NASA-sponsored Atmospheric Carbon and Transport – America (ACT-A) field campaign, this study systematically examines the regional uncertainties and biases of the carbon dioxide (CO2) concentrations from two of the state-of-the-art global analysis products, namely the real-time analysis from the European Center (EC) for Medium Range Forecasting and NOAA’s near real-time Carbon Tracker (CT) reanalysis. It is found that both the EC and CT-NRT analyses agree reasonably well with the independent ACT-A flight-level CO2 measurements in the free troposphere but the uncertainties are considerably larger in the boundary layer during both the summer months of 2016 and the winter months of 2017. There are also strong variabilities in accuracy and bias between seasons, and across three different subregions in the United States (Mid-Atlantic, Midwest and South). Overall, the analysis uncertainties of the EC and CT-NRT analyses in terms of root-mean square deviations against airborne data are comparable to each other, both of which are between 1-2 ppm in the free troposphere but can be as large as 10 ppm near the surface, which are grossly consistent with the difference between the two analyses. The current study not only provides systematic uncertainty estimates for both analysis products over North America but also demonstrated that these two independent estimates can be used to approximate the overall regional CO2 analysis uncertainties. Both statistics are important in future studies in quantifying the uncertainties of regional carbon concentration and flux estimates, as well as in assessing the impact of regional transport through more refined regional modeling and analysis systems.

Chen, Hans W.↗

Allowable Trajectory Variations for Space Shuttle Orbiter Entry-Aeroheating CFD

Reynolds-number criteria are developed for acceptable variations in Space Shuttle Orbiter entry trajectories for use in computational aeroheating analyses. The criteria determine if an existing computational fluid dynamics solution for a particular trajectory can be extrapolated to a different trajectory. The criteria development begins by estimating uncertainties for seventeen types of computational aeroheating data, such as boundary layer thickness, at exact trajectory conditions. For each type of datum, the allowable uncertainty contribution due to trajectory variation is set to be half of the value of the estimated exact-trajectory uncertainty. Then, for the twelve highest-priority datum types, Reynolds-number relations between trajectory variation and output uncertainty are determined. From these relations the criteria are established for the maximum allowable trajectory variations. The most restrictive criterion allows a 25% variation in Reynolds number at constant Mach number between trajectories.

Wood, William A.↗

Optimal Estimation Retrievals and Their Uncertainties: What Every Atmospheric Scientist Should Know

Remote sensing instruments are heavily used to provide observations for both the operational and research communities. These sensors do not provide direct observations of the desired atmospheric variables, but instead, retrieval algorithms are necessary to convert the indirect observations into the variable of interest. It is critical to be aware of the underlying assumptions made by many retrieval algorithms, including that the retrieval problem is often ill posed and that there are various sources of uncertainty that need to be treated properly. In short, the retrieval challenge is to invert a set of noisy observations to obtain estimates of atmospheric quantities. The problem is often complicated by imperfect forward models, by imperfect prior knowledge, and by the existence of nonunique solutions. Optimal estimation (OE) is a widely used physical retrieval method that combines measurements, prior information, and the corresponding uncertainties based on Bayes’s theorem to find an optimal solution for the atmospheric state. Furthermore, OE also allows the relative contributions of the different sources of error to the uncertainty in the final retrieved atmospheric state to be understood. Here, we provide a novel Python library to illustrate the use of OE for inverse problems in the atmospheric sciences. We introduce two example problems: how to retrieve drop size distribution parameters from radar observations and how to retrieve the temperature profile from ground-based microwave sensors. Using these examples, we discuss common pitfalls, how the various error sources impact the retrieval, and how the quality of the retrieval results can be quantified.

54 ENVIRONMENTAL SCIENCES↗

Decision Science for Machine Learning (DeSciML)

The increasing use of machine learning (ML) models to support high-consequence decision making drives a need to increase the rigor of ML-based decision making. Critical problems ranging from climate change to nonproliferation monitoring rely on machine learning for aspects of their analyses. Likewise, future technologies, such as incorporation of data-driven methods into the stockpile surveillance and predictive failure analysis for weapons components, will all rely on decision-making that incorporates the output of machine learning models. In this project, our main focus was the development of decision scientific methods that combine uncertainty estimates for machine learning predictions, with a domain-specific model of error costs. Other focus areas include uncertainty measurement in ML predictions, designing decision rules using multiobjecive optimization, the value of uncertainty reduction, and decision-tailored uncertainty quantification for probability estimates. By laying foundations for rigorous decision making based on the predictions of machine learning models, these approaches are directly relevant to every national security mission that applies, or will apply, machine learning to data, most of which entail some decision context.

97 MATHEMATICS AND COMPUTING↗

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗

Methodology for deriving the telescope focus function and its uncertainty for a heterodyne pulsed Doppler lidar

Doppler lidars provide two measured parameters, radial velocity and signal-to-noise ratio, from which winds and turbulent properties are routinely derived. Attenuated backscatter, which gives quantitative information on aerosols, clouds, and precipitation in the atmosphere, can be used in conjunction with the winds and turbulent properties to create a sophisticated classification of the state of the atmospheric boundary layer. Calculating attenuated backscatter from the signal-to-noise ratio requires accurate knowledge of the telescope focus function, which is usually unavailable. Inaccurate assumptions of the telescope focus function can significantly deform attenuated backscatter profiles, even if the instrument is focused at infinity. Here, we present a methodology for deriving the telescope focus function using a co-located ceilometer for pulsed heterodyne Doppler lidars. The method was tested with Halo Photonics StreamLine and StreamLine XR Doppler lidars but should also be applicable to other pulsed heterodyne Doppler lidar systems. The method derives two parameters of the telescope focus function, the effective beam diameter and the effective focal length of the telescope. Additionally, the method provides uncertainty estimates for the retrieved attenuated backscatter profile arising from uncertainties in deriving the telescope function, together with standard measurement uncertainties from the signal-to-noise ratio. The method is best suited for locations where the absolute difference in aerosol extinction at the ceilometer and Doppler lidar wavelengths is small.

54 ENVIRONMENTAL SCIENCES↗

Deep learning forecasts the spatiotemporal evolution of fluid-induced microearthquakes

Microearthquakes generated by subsurface fluid injection record the evolving stress state and permeability of reservoirs. Forecasting their spatiotemporal evolution is therefore critical for applications such as enhanced geothermal systems, carbon dioxide sequestration and other geoengineering applications. Here we propose a transformer neural network model that ingests hydraulic stimulation history and prior microearthquake observations to forecast four key quantities: cumulative microearthquake count, cumulative logarithmic seismic moment, and the 50th- and 95th-percentile extents of the microearthquake cloud. Applied to the EGS Collab Experiment 1 dataset, the model achieves R2 > 0.98 for the 1-s forecast horizon and R2 > 0.88 for the 15-s forecast horizon across all targets, and supplies uncertainty estimates through a learned standard deviation term. These accurate, uncertainty-quantified forecasts enable real-time inference of fracture propagation and permeability evolution, demonstrating the strong potential of deep-learning approaches to improve seismic-risk assessment and guide mitigation strategies in future fluid-injection operations.

Chung, Jaehong↗

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

A Side by Side Comparison of Filter-Based PM(sub 2.5) Measurements at a Suburban Site: A Closure Study

Reliable determination of the effects of air quality on public health and the environment requires accurate measurement of PM(sub 2.5) mass and the individual chemical components of fine aerosols. This study seeks to evaluate PM(sub 2.5) measurements that are part of a newly established national network by comparing them with a more conventional sampling system. Experiments were carried out during 2002 at a suburban site in Maryland, United States, where two samplers from the U.S. Environmental Protection Agency (USEPA) Speciation Trends Network: Met One Speciation Air Sampling System STNS and Thermo Scientific Reference Ambient Air Sampler STNR, two Desert Research Institute Sequential Filter Samplers DRIF, and a continuous TEOM monitor (Thermo Scientific Tapered Element Oscillating Microbalance) were sampling air in parallel. These monitors differ not only in sampling configuration but also in protocol-specific sample analysis procedures. Measurements of PM(sub 2.5) mass and major contributing species were well correlated among the different methods with r-values > 0.8. Despite the good correlations, daily concentrations of PM(sub 2.5) mass and major contributing species were significantly different at the 95% confidence level from 5 to 100% of the time. Larger values of PM(sub 2.5) mass and individual species were generally reported from STNR and STNS. The January STNR average PM(sub 2.5) mass (8.8 (micro)g/per cubic meter) was 1.5 (micro)g/per cubic meter larger than the DRIF average mass. The July STNS average PM(sub 2.5) mass (27.8 (micro)g/per cubic meter) was 3.8 (micro)g/per cubic meter larger than the DRIF average mass. These differences can only be partially accounted for by known random errors. Variations in flow control, face velocity, and sampling artifacts likely influence the measurement of PM(sub 2.5) speciation and mass closure. Simple statistical tests indicate that the current uncertainty estimates used in the STN network may underestimate the actual uncertainty.

Haines, Jennifer C.↗

A Satellite Mortality Study to Support Space Systems Lifetime Prediction

Estimating the operational lifetime of satellites and spacecraft is a complex process. Operational lifetime can differ from mission design lifetime for a variety of reasons. Unexpected mortality can occur due to human errors in design and fabrication, to human errors in launch and operations, to random anomalies of hardware and software or even satellite function degradation or technology change, leading to unrealized economic or mission return. This study focuses on data collection of public information using, for the first time, a large, publically available dataset, and preliminary analysis of satellite lifetimes, both operational lifetime and design lifetime. The objective of this study is the illustration of the relationship of design life to actual lifetime for some representative classes of satellites and spacecraft. First, a Weibull and Exponential lifetime analysis comparison is performed on the ratio of mission operating lifetime to design life, accounting for terminated and ongoing missions. Next a Kaplan-Meier survivor function, standard practice for clinical trials analysis, is estimated from operating lifetime. Bootstrap resampling is used to provide uncertainty estimates of selected survival probabilities. This study highlights the need for more detailed databases and engineering reliability models of satellite lifetime that include satellite systems and subsystems, operations procedures and environmental characteristics to support the design of complex, multi-generation, long-lived space systems in Earth orbit.

Fox, George↗